Has anyone tried MintPy for SBAS using SNAP?

s1tbx
#21

Here is how you prep the data. You need to do some coding in Python

Prepare data like this directory structure.

  1. Install hyp3-lib from the Github repo using the setup.py file
  2. Clip raster data to the same extent:
import glob
from osgeo import gdal
import os
from hyp3lib.cutGeotiffs import cutFiles


# Path that contains the directories with the unzipped Hyp3 data
images_folder = r'E:\London\Hyp3'

glob_items = []
target_files = ['*corr.tif', '*unw_phase.tif', '*inc_map.tif', '*dem.tif']
for item in target_files:
    for file in glob.glob(os.path.join(images_folder, '*', item)):
        glob_items.append(file)

print('Clipping images')
cutFiles(glob_items)
  1. Prepare config file using glob paths
mintpy.load.processor        = hyp3
##---------interferogram datasets:
mintpy.load.unwFile          = E:/London/Hyp3/*/*unw_phase_clip.tif
mintpy.load.corFile          = E:/London/Hyp3/*/*corr_clip.tif
##---------geometry datasets:
mintpy.load.demFile          = E:/London/Hyp3/*/*dem_clip.tif
mintpy.load.incAngleFile     = E:/London/Hyp3/*/*inc_map_clip.tif
  1. Run MintPy smallbaselineApp using config file as input.
4 Likes
#22

Is this for SBAS?
How can I select my study area?

#23

@suribabu

Yes this is for SBAS.

If you want to subset your Hyp3 data in your config file you should add:

mintpy.subset.lalo = S:N,W:E

S,N,W,E being the extent of the bounding box. This is based on the projection of the Hyp3 data which is UTM.

2 Likes
#25

A sidenote:

Please get acquainted with the open source tool PyRate for SBAS, SNAP9 will support it in the same manner as StaMPS today:

PyRate documentation — PyRate 0.6.1 documentation (geoscienceaustralia.github.io)

GitHub - GeoscienceAustralia/PyRate: A Python tool for estimating velocity and time-series from Interferometric Synthetic Aperture Radar (InSAR) data.

3 Likes
#26

Hey, just curious. why 1master and 3slaves. Is this a general network for sbas, or any procedure set the requriement

#27

It’s not required. You could draw your connected network whatever you want. Redundancy should be better.

#28

is there any rule for pairing the images or any experience we can use. Maybe we just remove the pairs with low quality inferogram, but I try to save some time…

#29

There are two typical ways to do it:

  1. For each image, calculate interferograms with n closest images, 3 for instance.
  2. Set a threshold for temporal and spatial baselines, something like 130m and 36 days.

You can use the SBAS tool at ASF Vertex in order to draw the network by yourself, or downloading a network based on thresholds.

#30

hey there,
I’m also preparing datafiles for sbas using snap.
If I understand correctly, all data should be co-registrated to a chosen master data firstly and then generate inteferogram pairs. This is not possible in snap for now because multi-master is not supported.
Your method is to clip all interferograms into same area and dimensions. This mainly rely on the geocoding accuracy. I guess it may have some problems. Apply co-regsitration at early phase maybe more reasonable

#31

You have to process them all individually one by one. This is a very labor intensive process if you do it manually in SNAP. I really recommend you create a script for it.

I just do standard InSAR processing that you can find on RUS Training. I automated that process to fit the SBAS requirement where each date needs 3 slaves.

2 Likes
#32

If it possible, please explain how to process SBAS in SNAP and StaMPS?

#33

A modified version of StaMPS has been provided here: Please test the SBAS_snap2stamps
You can try to follow along these instructions.

1 Like
#34

Hi dear ABraun
I want to calculate the closure_phase_bias.py using the MintPy software that I per-processed on ISCE,
but I cant calculate with (MintPy/closure_phase_bias.py at main · insarlab/MintPy · GitHub)
I would be grateful if you could guide me

#35

I haven’t tried MintPy myself yet, sorry.

#36

Hi dear panjibrotoisworo
I want to calculate the closure_phase_bias.py using the MintPy software that I per-processed on ISCE,
but I cant calculate with (MintPy/closure_phase_bias.py at main · insarlab/MintPy · GitHub)
I would be grateful if you could guide me

#37

thanks dear ABraun

#38

Hello, thanks for your sharing. but I wonder how to Subset all four MintPy products by region, to make all the input data have the same row and column.