How to prepare InSAR coregistered stack from ISCE2 --> dolphin --> Mintpy. geometryGeo.h5 creation gives error.

[Karthi J]

Hi, everyone,

                     I am getting error while trying to process the dolphin's output in mintpy. Can you please give me insights and solution.

Can we process the all workflows with the Mintpy?

I used ISCE 2.6.3 & Mintpy 1.6.2 & 1.6.3 to make CSLC stack.

I used 

"stackSentinel.py -s .../slc -o .../orbits -d .../demLat_N19_N22_Lon_E084_E087.dem.wgs84 -b '19.760 19.917 85.725 85.935' -a .../orbits -c 5 -C NESD -W slc -m 20160101 --num_proc 32 --num_proc4topo 32 -e 0.30" 

I then ran dolphin with the 

"cslc_file_list:
  - .../merged/SLC/20160101/20160101.slc.full.vrt
  - .../merged/SLC/20160113/20160113.slc.full.vrt
  - .../merged/SLC/20160125/20160125.slc.full.vrt"

The thing is, I ran dolphin with full resolution. the strides dolphin run dolpin_config.yaml strides 1 1

I imported the files using  

"########## 1. Load Data (--load to exit after this step)
## load_data.py -H to check more details and example inputs.
mintpy.load.processor        = isce
##---------for ISCE only:
mintpy.load.metaFile         = .../coregistered/reference/IW*.xml
mintpy.load.baselineDir      = .../coregistered/baselines
##---------interferogram datasets:
mintpy.load.unwFile          = .../dolphin/unwrapped/*.unw.tif
mintpy.load.corFile          = .../dolphin/interferograms/*.int.cor.tif
mintpy.load.connCompFile     = .../dolphin/unwrapped/*.unw.conncomp.tif
##---------geometry datasets:
mintpy.load.demFile          = .../coregistered/merged/geom_reference/hgt.rdr
mintpy.load.lookupYFile      = .../coregistered/merged/geom_reference/lat.rdr
mintpy.load.lookupXFile      = .../coregistered/merged/geom_reference/lon.rdr
mintpy.load.incAngleFile     = .../coregistered/merged/geom_reference/los.rdr
mintpy.load.azAngleFile      = .../coregistered/merged/geom_reference/los.rdr
mintpy.load.shadowMaskFile   = .../coregistered/merged/geom_reference/shadowMask.rdr" 

I get error like "Traceback (most recent call last):
File "/home/incois/mohanty/.conda/envs/isce2/bin/smallbaselineApp.py", line 10, in
sys.exit(main())
File "/home/incois/mohanty/.conda/envs/isce2/lib/python3.9/site-packages/mintpy/cli/smallbaselineApp.py", line 209, in main
run_smallbaselineApp(inps)
File "/home/incois/mohanty/.conda/envs/isce2/lib/python3.9/site-packages/mintpy/smallbaselineApp.py", line 1155, in run_smallbaselineApp
app.run(steps=inps.runSteps)
File "/home/incois/mohanty/.conda/envs/isce2/lib/python3.9/site-packages/mintpy/smallbaselineApp.py", line 908, in run
self.run_load_data(sname)
File "/home/incois/mohanty/.conda/envs/isce2/lib/python3.9/site-packages/mintpy/smallbaselineApp.py", line 182, in run_load_data
stack_file, geom_file, _, ion_file = ut.check_loaded_dataset(self.workDir, print_msg=True)[:4]
File "/home/incois/mohanty/.conda/envs/isce2/lib/python3.9/site-packages/mintpy/utils/utils.py", line 91, in check_loaded_dataset
raise FileNotFoundError(errno.ENOENT, os.strerror(errno.ENOENT), geom_file + msg_aria)
FileNotFoundError: [Errno 2] No such file or directory: '/home/incois/mohanty/karthi/SBAS/dolphin/mintpy_1/inputs/geometryGeo.h5'"" 

Here's the detailed error log. 

Even if I directly geocode the geometryRadar.h5 into geometryGeo.h5 using geocode.py the size mismatch occurs. like this 

"IndexError: boolean index did not match indexed array along dimension 0; dimension is 1389 but corresponding boolean dimension is 266"