bug introduced by subset.py PR #1369
Description
Hi @yunjunz ,
I am afraid the PR for subset.py to support radar coordinates introduced a bug. I get:
smallbaselineApp.py ./unittestGalapagosSenDT128.template --dir mintpy
multilook method : nearest
Traceback (most recent call last):
File "/Users/famelung/code/rsmas_insar/tools/miniforge3/bin/smallbaselineApp.py", line 8, in <module>
sys.exit(main())
File "/Users/famelung/code/rsmas_insar/tools/MintPy/src/mintpy/cli/smallbaselineApp.py", line 209, in main
run_smallbaselineApp(inps)
File "/Users/famelung/code/rsmas_insar/tools/MintPy/src/mintpy/smallbaselineApp.py", line 1155, in run_smallbaselineApp
app.run(steps=inps.runSteps)
File "/Users/famelung/code/rsmas_insar/tools/MintPy/src/mintpy/smallbaselineApp.py", line 908, in run
self.run_load_data(sname)
File "/Users/famelung/code/rsmas_insar/tools/MintPy/src/mintpy/smallbaselineApp.py", line 176, in run_load_data
mintpy.cli.load_data.main(iargs)
File "/Users/famelung/code/rsmas_insar/tools/MintPy/src/mintpy/cli/load_data.py", line 121, in main
load_data(inps)
File "/Users/famelung/code/rsmas_insar/tools/MintPy/src/mintpy/load_data.py", line 825, in load_data
iDict = read_subset_box(iDict)
File "/Users/famelung/code/rsmas_insar/tools/MintPy/src/mintpy/load_data.py", line 197, in read_subset_box
pix_box = coord.bbox_geo2radar(geo_box)
File "/Users/famelung/code/rsmas_insar/tools/MintPy/src/mintpy/objects/coord.py", line 502, in bbox_geo2radar
y, x, y_res, x_res = self.geo2radar(lat, lon, print_msg=print_msg)
File "/Users/famelung/code/rsmas_insar/tools/MintPy/src/mintpy/objects/coord.py", line 364, in geo2radar
az[i], rg[i] = self._get_lookup_row_col(
File "/Users/famelung/code/rsmas_insar/tools/MintPy/src/mintpy/objects/coord.py", line 229, in _get_lookup_row_col
raise RuntimeError(f'No corresponding coordinate found for y/x: {y}/{x}')
RuntimeError: No corresponding coordinate found for y/x: -0.81/-91.13
The previous version (git checkout cc2cc9a8) works fine.
You can reproduce the error using:
wget http://149.165.154.65:/data/HDF5EOS/unittestGalapagosSenDT128/unittestGalapagosSenDT128.tar
tar xvf unittestGalapagosSenDT128.tar
cd unittestGalapagosSenDT128
smallbaselineApp.py ./unittestGalapagosSenDT128.template --dir mintpy
smallbaselineApp.py --version
The version I commited for the PR works fine (https://github.com/falkamelung/MintPy).
It might be related to buf. At
MintPy/src/mintpy/objects/coord.py
Lines 364 to 370 in 3f134ca
the az[i], rg[i] = self._get_lookup_row_col(...) function works for i=2 and the error occurs for i=3. This is how far I came.