Processing Hyp3 data

[Benjy Marks]

Hi,

I am attempting to build a work flow for my students to process insar data via hyp3. I have set up the following docker container to try to streamline the process: https://github.com/benjym/insar

I am using the following configuration file:

mintpy.load.processor        = hyp3

##---------interferogram datasets:

mintpy.load.unwFile          = /home/work/data/interferograms/*/*unw_phase_clip.tif

mintpy.load.corFile          = /home/work/data/interferograms/*/*corr_clip.tif

##---------geometry datasets:

mintpy.load.demFile          = /home/work/data/DEM/dem_clip.tif

After downloading some processed hyp3 insar files and running smallbaselineApp.py, everything works fine up to the tropospheric correction stage where I get the following error:

INFO: You are using the latest ECMWF platform for downloading datasets: https://cds.climate.copernicus.eu/api/v2

Downloading 1 of 3: /home/work/ERA5/ERA5_N6124550_N6297610_E292790_E390330_20150303_19.grb 

{'product_type': 'reanalysis', 'format': 'grib', 'variable': ['geopotential', 'temperature', 'specific_humidity'], 'pressure_level': ['1', '2', '3', '5', '7', '10', '20', '30', '50', '70', '100', '125', '150', '175', '200', '225', '250', '300', '350', '400', '450', '500', '550', '600', '650', '700', '750', '775', '800', '825', '850', '875', '900', '925', '950', '975', '1000'], 'year': '2015', 'month': '03', 'day': '03', 'time': '19:00', 'area': '6297610.00/292790.00/6124550.00/390330.00'}

I believe this is related to the area being passed in the wrong units, which stems from the coordinate system being wrong in ifgramStack.h5 . When I run info.py inputs/ifgramStack.h5 I get the following (where I believe the coordinates are in fact RADAR but are listed as GEO?). Can you please offer some advice as to what I should try next?

******************** Basic File Info ************************

file name: /home/work/inputs/ifgramStack.h5

file type: ifgramStack

coordinates : GEO

******************** HDF5 File Structure ********************

Attributes in / level:

  ALOOKS                   4

  ANTENNA_SIDE             -1

  CENTER_LINE_UTC          68837.5

  DATA_TYPE                float32

  DATE12                   150303-190728

  EARTH_RADIUS             6371000.0

  FILE_LENGTH              2164

  FILE_PATH                /home/work/data/interferograms/S1AA_20150303T190658_20190728T190737_VVP1608_INT80_G_ueF_EA73/S1AA_20150303T190658_20190728T190737_VVP1608_INT80_G_ueF_EA73_unw_phase_clip.tif

  FILE_TYPE                ifgramStack

  HEADING                  193.1769

  HEIGHT                   693000.0

  LAT_REF1                 6124480.0

  LAT_REF2                 6124480.0

  LAT_REF3                 6297600.0

  LAT_REF4                 6297600.0

  LENGTH                   2164

  LON_REF1                 292800.0

  LON_REF2                 390400.0

  LON_REF3                 292800.0

  LON_REF4                 390400.0

  ORBIT_DIRECTION          ASCENDING

  PLATFORM                 Sen

  PROCESSOR                hyp3

  PROJECT_NAME             default

  P_BASELINE_BOTTOM_HDR    49.6686

  P_BASELINE_TOP_HDR       49.6686

  REF_LAT                  6249720.0

  REF_LON                  324680.0

  REF_X                    398

  REF_Y                    598

  RLOOKS                   20

  UTM_ZONE                 56S

  WAVELENGTH               0.055465764662349676

  WIDTH                    1220

  X_FIRST                  292800.0

  X_STEP                   80.0

  X_UNIT                   meters

  Y_FIRST                  6297600.0

  Y_STEP                   -80.0

  Y_UNIT                   meters

  mintpy.load.corFile      /home/work/data/interferograms/*/*corr_clip.tif

  mintpy.load.demFile      /home/work/data/DEM/dem_clip.tif

  mintpy.load.processor    hyp3

  mintpy.load.unwFile      /home/work/data/interferograms/*/*unw_phase_clip.tif

  number_bands             1

  scheme                   BSQ

HDF5 dataset "/bperp               ": shape (3,)                , dtype <float32>

HDF5 dataset "/coherence           ": shape (3, 2164, 1220)     , dtype <float32>

  MODIFICATION_TIME    1622598769.8476129

HDF5 dataset "/date                ": shape (3, 2)              , dtype <|S8>

HDF5 dataset "/dropIfgram          ": shape (3,)                , dtype <bool>

HDF5 dataset "/unwrapPhase         ": shape (3, 2164, 1220)     , dtype <float32>

  MODIFICATION_TIME    1622598769.315683

thanks,

benjy

[Zhang Yunjun]

Hi benjy,

The tropospheric delay support of hyp3 in mintpy is just added (you would need to load an incidence angle file for this purpose), you could follow the progress on the github issue here: https://github.com/insarlab/MintPy/issues/540.

Your GitHub description and setup looks very nice! I would encourage you to put a link to the mintpy github issue above, it might help the ASF folks to be aware of it and foster communications and problem-solving.

Regards,

Yunjun