ionospheric corrections in MintPy

[EJFielding]

Is the new ionospheric correction analysis in MintPy designed only for the ISCE2 topsStack (stackSentinel) ionospheric corrections, or does it support the alosStack, and stripmapStack ionospheric corrections?

[yunjunz]

We have tested it on alosStack and topsStack. But it should work for all isce2 stack processors. If that's not the case, please don't hesitate to open an issue and let us know.

[Qori88]

Hi Yunjun and Eric,
I implement correct ionosphere for alos1 from stackstripmap -W ionosphere. For mintpy.ionosphericDelay.method I use split_spectrum. But I found this error :

Estimate ionospheric delay time-series via ifgram_inversion.py ...
ifgram_inversion.py /media/imrsl/HDD2/TInSAR/Kalimantan/Pre/12_ion/mintpy3/inputs/ionStack.h5 -t /media/imrsl/HDD2/TInSAR/Kalimantan/Pre/12_ion/mintpy3/smallbaselineApp.cfg -w no --update
read input option from template file: /media/imrsl/HDD2/TInSAR/Kalimantan/Pre/12_ion/mintpy3/smallbaselineApp.cfg
phase unwrapping error correction "bridging+phase_closure" is turned ON
use dataset "unwrapPhase_bridging_phaseClosure" by default
Traceback (most recent call last):
File "/home/imrsl/anaconda3/envs/mintpy/bin/smallbaselineApp.py", line 10, in
sys.exit(main())
^^^^^^
File "/home/imrsl/anaconda3/envs/mintpy/lib/python3.11/site-packages/mintpy/cli/smallbaselineApp.py", line 209, in main
run_smallbaselineApp(inps)
File "/home/imrsl/anaconda3/envs/mintpy/lib/python3.11/site-packages/mintpy/smallbaselineApp.py", line 1155, in run_smallbaselineApp
app.run(steps=inps.runSteps)
File "/home/imrsl/anaconda3/envs/mintpy/lib/python3.11/site-packages/mintpy/smallbaselineApp.py", line 932, in run
self.run_ionospheric_delay_correction(sname)
File "/home/imrsl/anaconda3/envs/mintpy/lib/python3.11/site-packages/mintpy/smallbaselineApp.py", line 611, in run_ionospheric_delay_correction
mintpy.cli.iono_split_spectrum.main(iargs)
File "/home/imrsl/anaconda3/envs/mintpy/lib/python3.11/site-packages/mintpy/cli/iono_split_spectrum.py", line 145, in main
run_iono_split_spectrum(inps)
File "/home/imrsl/anaconda3/envs/mintpy/lib/python3.11/site-packages/mintpy/iono_split_spectrum.py", line 42, in run_iono_split_spectrum
ifgram_inversion.main(cmd.split()[1:])
File "/home/imrsl/anaconda3/envs/mintpy/lib/python3.11/site-packages/mintpy/cli/ifgram_inversion.py", line 264, in main
inps = cmd_line_parse(iargs)
^^^^^^^^^^^^^^^^^^^^^
File "/home/imrsl/anaconda3/envs/mintpy/lib/python3.11/site-packages/mintpy/cli/ifgram_inversion.py", line 178, in cmd_line_parse
raise ValueError(msg)
ValueError: input dataset name "unwrapPhase_bridging_phaseClosure" not found in file: /media/imrsl/HDD2/TInSAR/Kalimantan/Pre/12_ion/mintpy3/inputs/ionStack.h5

any guidance, I appreciate

[Qori88]

Hi @yunjunz,
I add your code into estimate iono time-series part :
# 3. estimate iono time-series
print('\n'+'-'*80)
print('Estimate ionospheric delay time-series via ifgram_inversion.py ...')
cmd = f'ifgram_inversion.py {inps.iono_stack_file} -t {inps.template_file} -w no --update'
cmd += ' --dset unwrapPhase --mask-dset no '
print(cmd)
ifgram_inversion.main(cmd.split()[1:])

And the error is changed to :

[Qori88]

BTW, I am not sure if your input iono files are correct or not (have not tried iono correction within stripmapStack for many years), please ensure.

I also tried to change the ionosphere stack with interferogram stack that also have unw interferogram, spatial coh and conncomp but the error is the same

[yunjunz]

Could you try changing the following again?

from

MintPy/src/mintpy/iono_split_spectrum.py

Line 40 in d2e551a

cmd = f'ifgram_inversion.py {inps.iono_stack_file} -t {inps.template_file} -w no --update'

to

    cmd = f'ifgram_inversion.py {inps.iono_stack_file} --dset unwrapPhase --weight-func no --update'

[Qori88]

Yunjun, it works with your last change suggestion. Thank you very much.
I continue to the last step but the final result is still messy. I will check my .cfg file again.

[yunjunz]

The fix is now in #1239

[Qori88]

here is my .txt file
alos_pre12_ion.txt