Skip to content

view.py fails when using large multilook factor #950

@EJFielding

Description

@EJFielding

Description of the problem

The view.py script is not correctly adjusting the location of the reference point when large multilooking factor is used for the display and then fails when original location is outside the multilooked image.

Full script that generated the error

view.py inputs/ifgramStack.h5 unwrapPhase --save

Full error message

run view.py in MintPy version 1.5.1.post6, date 2023-01-22
input file is ifgramStack file: /u/sar-r0/fielding/Calif/SF_Bay/UAVSAR/Stacks/Haywrd_15302_03/s4/MintPy/inputs/ifgramStack.h5 in float32 format
file size in y/x: (7345, 4950)
input dataset: "['unwrapPhase']"
num of datasets in file ifgramStack.h5: 330
num of datasets to exclude: 0
num of datasets to display: 198
data   coverage in y/x: (0, 0, 4950, 7345)
subset coverage in y/x: (0, 0, 4950, 7345)
data   coverage in lat/lon: None
subset coverage in lat/lon: None
------------------------------------------------------------------------
colormap: jet
figure title: unwrapPhase
figure size : [15.00, 8.00]
dataset number: 198
row     number: 6
column  number: 17
figure  number: 2
total number of pixels: 3.7E+09
* multilook 12 by 12 with nearest interpolation for display to save memory
consider reference pixel in y/x: [200, 1800]
mark interferograms with 'dropIfgram=False' in red colored title
----------------------------------------
Figure 1 - unwrapPhase_1.png
reading data as a list of 2D matrices ...
Traceback (most recent call last):
  File "/u/vento-r0/ericf/packages/Anaconda/anaconda3/envs/insar2022/bin/view.py", line 8, in <module>
    sys.exit(main())
  File "/u/pez-z2/fielding/tools/MintPy/src/mintpy/cli/view.py", line 172, in main
    obj.plot()
  File "/u/pez-z2/fielding/tools/MintPy/src/mintpy/view.py", line 1611, in plot
    plot_figure(j, self, metadata=self.atr)
  File "/u/pez-z2/fielding/tools/MintPy/src/mintpy/view.py", line 1286, in plot_figure
    data = read_data4figure(i_start, i_end, inps, metadata)
  File "/u/pez-z2/fielding/tools/MintPy/src/mintpy/view.py", line 1074, in read_data4figure
    d[d!=0] -= d[ref_y, ref_x]
IndexError: index 1800 is out of bounds for axis 1 with size 412

System information

  • Operating system: Linux-CentOS
  • Python environment: conda
  • MintPy version: MintPy version 1.5.1.post6, date 2023-01-22
  • Your custom / default template file (if the bug is related to a specific dataset):

Metadata

Metadata

Assignees

No one assigned

    Labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions