UPDATE : New optionnal output (ortho geometry) for DiapOTB chain

python_src/diapOTB.py

@@ -143,6 +143,7 @@ if __name__ == "__main__":
     ml_geoGrid_range = ml_range
     ml_geoGrid_azimut = ml_azimut
     gain_interfero = float(ConfigSectionMap(Config, "DIn_SAR")['interferogram_gain'])
+    activateOrthoInterferogram = str2bool(ConfigSectionMap(Config, "DIn_SAR")['interferogram_ortho'])
 
     if (geoGrid_threshold < 0) or (geoGrid_threshold > 1) :
         print("Wrong Threshold for fine deformation grid")
@@ -380,6 +381,7 @@ if __name__ == "__main__":
     ######## SARRobustInterferogram Application (interf step) #######
     print("\n SARRobustInterferogram Application \n")
     interferogram = "interferogram.tif"
+    interferogram_ortho = "interferogram_ortho.tif"
     appInterferogram = otb.Registry.CreateApplication("SARRobustInterferogram")
     appInterferogram.SetParameterString("insarmaster", master_Image)
     appInterferogram.SetParameterString("insarslave", os.path.join(output_dir, slave_Image_CoRe))
@@ -395,4 +397,8 @@ if __name__ == "__main__":
     appInterferogram.SetParameterFloat("gain", gain_interfero)
     appInterferogram.SetParameterString("ram", "4000")
     appInterferogram.SetParameterString("out", os.path.join(output_dir, interferogram))
+    if activateOrthoInterferogram :
+        appInterferogram.SetParameterEmpty("ortho", True)
+        appInterferogram.SetParameterString("indem", dem)
+        appInterferogram.SetParameterString("outopt", os.path.join(output_dir, interferogram_ortho))
     appInterferogram.ExecuteAndWriteOutput()

python_src/ex_config.ini

@@ -21,4 +21,5 @@ GridStep_range = 150
 GridStep_azimut = 150
 Grid_Threshold = 0.3
 Grid_Gap = 0.7
-Interferogram_gain = 0.
\ No newline at end of file
+Interferogram_gain = 0.
+Interferogram_ortho = yes
\ No newline at end of file