Skip to content

Spatial subset not working for ARIA data #469

@rzinke

Description

@rzinke

Description of the problem
I have a stack of ARIA interferograms I am tyring to process. I am trying to subset the data set according to an area of interest, so I specify the lat/lon bounds using mintpy.subset.lalo. However, when I process the data set, it processes the entire spatial extent, leading to frequent memory overloads. Does MintPy support spatial subsetting for ARIA data?

*Loading data

Running load_data.py -t smallbaselineapp.cfg Parameters.txt --enforce produces the following output:

...
mintpy.load.lookupYFile    = auto  #[path of latitude /row   /y coordinate file], not required for geocoded data
mintpy.load.lookupXFile    = auto  #[path of longitude/column/x coordinate file], not required for geocoded data
mintpy.load.incAngleFile   = auto  #[path of incidence angle file], optional but recommend
mintpy.load.azAngleFile    = auto  #[path of azimuth   angle file], optional
mintpy.load.shadowMaskFile = auto  #[path of shadow mask file], optional but recommend
mintpy.load.waterMaskFile  = auto  #[path of water  mask file], optional but recommend
mintpy.load.bperpFile      = auto  #[path pattern of 2D perpendicular baseline file], optional
##---------subset (optional):
## if both yx and lalo are specified, use lalo option unless a) no lookup file AND b) dataset is in radar coord
mintpy.subset.yx           = auto    #[y0:y1,x0:x1 / no], auto for no
mintpy.subset.lalo         = 36.0:44.1,86.5:91.5    #[lat0:lat1,lon0:lon1 / no], auto for no


########## 2. modify_network
## reference: Yunjun et al. (2019, section 4.2 and 5.3.1); Chaussard et al. (2015, GRL)
## 1) Coherence-based network modification = (threshold + MST) by default
## It calculates a average coherence value for each interferogram using spatial coherence and input mask (with AOI)
## Then it finds a minimum spanning tree (MST) network with inverse of average coherence as weight (keepMinSpanTree)
(MintPy) [rzinke@sar A1_ML3_ERA5_Topo_PxMask]$ load_data.py -t smallbaselineApp.cfg Parameters41A.txt --enforce
processor : isce
SAR platform/sensor : unknown from project name "None"
--------------------------------------------------
prepare metadata files for isce products
/u/sarh0/rzinke/tools/MintPy/mintpy/load_data.py:588: UserWarning: No input metadata file found: auto
  warnings.warn('No input metadata file found: {}'.format(inpsDict['mintpy.load.metaFile']))
WARNING: mintpy.subset.lalo is not supported if 1) no lookup file AND    2) radar/unkonwn coded dataset
        ignore it and continue.
--------------------------------------------------
searching interferometric pairs info
input data files:
WARNING: No reqired ['unwrapPhase', 'azimuthOffset'] data files found!
--------------------------------------------------
searching geometry files info
input data files:
WARNING: No reqired height data files found!
--------------------------------------------------
updateMode : True
compression: None

Does the message

WARNING: mintpy.subset.lalo is not supported if 1) no lookup file AND    2) radar/unkonwn coded dataset
        ignore it and continue.

indicate that a lookup file is needed for ARIA data to be spatially subset?

Metadata

Metadata

Assignees

Labels

No labels
No labels

Type

No type

Projects

No projects

Milestone

No milestone

Relationships

None yet

Development

No branches or pull requests

Issue actions