AutoRIFT + Geogrid Issue with Sentinel-1 TOPS Mode Processing (posted on the AutoRIFT github channel too, looking for some advice)

[Sci123-lab]

Dear InSAR Community,

I have a question regarding the use of AutoRIFT with Geogrid.

I have processed Sentinel-1 TOPS-mode data (subswaths 2 and 3) with denseoffsets set to true in the input file to generate full SLCs which are inputs for AutoRIFT.

(Note: I previously generated dense offsets using AMPCOR, and I observed azimuth and range offsets. Now, I am trying to reproduce these offsets using AutoRIFT to compare the performance of both methods. I assume that the error I am encountering with AutoRIFT is due to my lack of knowledge of the tool.)

I would like to share the error message I am receiving and ask if anyone has encountered and resolved this issue. Any insights or suggestions would be greatly appreciated!

AutoRIFT
When I run autoRIFT directly,here is what I get... autoRIFT starts and ends after a while but no files are being stored.

testautoRIFT_ISCE.py -m reference.slc.full -s secondary.slc.full

2025-03-13 22:23:12,352 - root - DEBUG - Creating readonly ISCE mmap with
file = reference.slc.full
bands = 1
width = 46747
length = 14049
scheme = BIP
dtype = CFLOAT

2025-03-13 22:23:12,356 - root - DEBUG - Creating readonly ISCE mmap with
file = secondary.slc.full
bands = 1
width = 46747
length = 14049
scheme = BIP
dtype = CFLOAT

Preprocessing filter width 5
Using preprocessing methods ['hps', 'hps']
Pre-process Start!!!
Using Wallis Filter Width: 5
Pre-process Done!!!
4.201286554336548
Uniform Data Type Done!!!
13.169963836669922

AutoRIFT Start!!!
/home/hassan/miniconda3/envs/isce2/lib/python3.8/site-packages/isce/components/contrib/geo_autoRIFT/autoRIFT/autoRIFT.py:1468: RuntimeWarning: divide by zero encountered in true_divide
DxMadmin = np.ones(Dx.shape) / OverSampleRatio / SearchLimitX * 2;
/home/hassan/miniconda3/envs/isce2/lib/python3.8/site-packages/isce/components/contrib/geo_autoRIFT/autoRIFT/autoRIFT.py:1469: RuntimeWarning: divide by zero encountered in true_divide
DyMadmin = np.ones(Dy.shape) / OverSampleRatio / SearchLimitY * 2;
/home/hassan/miniconda3/envs/isce2/lib/python3.8/site-packages/numpy/lib/nanfunctions.py:1119: RuntimeWarning: All-NaN slice encountered
r, k = function_base._ureduce(a, func=_nanmedian, axis=axis, out=out,
/home/hassan/miniconda3/envs/isce2/lib/python3.8/site-packages/isce/components/contrib/geo_autoRIFT/autoRIFT/autoRIFT.py:1399: RuntimeWarning: Mean of empty slice
C = np.nanmean(B,axis=0)
AutoRIFT Done!!!
79.58786702156067
/home/hassan/miniconda3/envs/isce2/lib/python3.8/site-packages/isce/applications/testautoRIFT_ISCE.py:396: DeprecationWarning: np.bool is a deprecated alias for the builtin bool. To silence this warning, use bool by itself. Doing this will not modify any behavior and is safe. If you specifically wanted the numpy scalar type, use np.bool_ here.
Deprecated in NumPy 1.20; for more details and guidance: https://numpy.org/devdocs/release/1.20.0-notes.html#deprecations
noDataMask = noDataMask.astype(np.bool)

What dos this output log mean

GeoGrid

I tried Geogrid command line code and specified reference secondary coregistered slcs that I transformed to acceptable format of Geogrid, and also a dem in xy cartesian coordinates (m) to do geogrid first but it fails with the error.
(I tried before transformation to, still the same error)

testGeogrid_ISCE.py -m /home/hassan/ISCE2-Tops/insar/merged/reference_grd.tif -s /home/hassan/ISCE2-Tops/insar/merged/secondary_grd.tif -d /home/hassan/ISCE2-Tops/insar/merged/dem_nztm.tif

The error i get is as follows:

Traceback (most recent call last):
File "/home/hassan/miniconda3/envs/isce2/lib/python3.8/site-packages/isce/applications/testGeogrid_ISCE.py", line 392, in
main()
File "/home/hassan/miniconda3/envs/isce2/lib/python3.8/site-packages/isce/applications/testGeogrid_ISCE.py", line 386, in main
metadata_m = loadMetadata(inps.indir_m)
File "/home/hassan/miniconda3/envs/isce2/lib/python3.8/site-packages/isce/applications/testGeogrid_ISCE.py", line 131, in loadMetadata
info.sensingStart = min([x.sensingStart for x in frames])
ValueError: min() arg is an empty sequence

Looking ahead to advices and people working currently on autoRIFT or have done in past. Thanks a bunch.