Unwrapping failed and zero amplitude after runing seuential.py

[SiriusXCH0103]

I process 16 Sentinel-1A images following the "workflows.md", but it fails every time when unwrapping. No files generated in "./PS_DS/unwrap" dir. The "run_unwrap_ps_ds.sh" command outputs the following message when unwrapping each wrapped file:

This is the Open Source version of ISCE.
Some of the workflows depend on a separate licensed package.
To obtain the licensed package, please make a request for ISCE
through the website: https://download.jpl.nasa.gov/ops/request/index.cfm.
Alternatively, if you are a member, or can become a member of WinSAR
you may be able to obtain access to a version of the licensed sofware at
https://winsar.unavco.org/software/isce
length, width:  3764   15603

snaphu v1.4.2
Reading wrapped phase from file PS_DS/20170708_20180727.int
No weight file specified.  Assuming uniform weights
Reading correlation data from file PS_DS/tcorr_ds_ps.bin
NaN or infinity found in correlation data
Abort

I checked output files of each step using "mdx.py". I found the output files "./Sequential/miniStacks/_/EVD/*.slc" are all zeros, and the maps are all black. As a reault, the "integratePS.py" outputs are also black.

What is the reason of this problem?

[otirmizi]

Hello,
Snaphu does not do well with NaN values. The first thing to do would be to double check and see if all your data is actually blank or if you are just zoomed in too far. You can add -z -8 to the end of the mdx command or use the zoom button from the toolbar to check.

If the files are actually all blank, you will need to recreate them using previous steps. If there is data there, you can edit the integratePS.py script and add ifgram_ds_ps = np.nan_to_num(ifgram_ds_ps) at line 129 and re-run it. This will replace all NaN values with zeroes and you should be able to unwrap after that.

[hy935]

您好，Snaphu 在 NaN 值方面表现不佳。首先要做的是仔细检查，看看你的所有数据是否真的是空白的，或者你是否只是放大得太远了。您可以添加到 mdx 命令的末尾，也可以使用工具栏中的缩放按钮进行检查。-z -8

如果文件实际上都是空白的，则需要使用前面的步骤重新创建它们。如果那里有数据，您可以编辑 integratePS.py 脚本并在第 129 行添加并重新运行它。这会将所有 NaN 值替换为零，之后您应该能够解包。ifgram_ds_ps = np.nan_to_num(ifgram_ds_ps)

Hello, I'm encountering the same issue, but after using your method, the error still persists.

sh run_unwrap_ps_ds.sh
This is the Open Source version of ISCE.
Some of the workflows depend on a separate licensed package.
To obtain the licensed package, please make a request for ISCE
through the website: https://download.jpl.nasa.gov/ops/request/index.cfm.
Alternatively, if you are a member, or can become a member of WinSAR
you may be able to obtain access to a version of the licensed sofware at
https://winsar.unavco.org/software/isce
/home/InSAR01/anaconda3/envs/insar/lib/python3.9/site-packages/osgeo/gdal.py:312: FutureWarning: Neither gdal.UseExceptions() nor gdal.DontUseExceptions() has been explicitly called. In GDAL 4.0, exceptions will be enabled by default.
warnings.warn(
length, width: 1784 4200

snaphu v1.4.2
Reading wrapped phase from file PS_DS/20230109_20230121.int
No weight file specified. Assuming uniform weights
Reading correlation data from file PS_DS/tcorr_ds_ps.bin
NaN or infinity found in correlation data
Abort