How to use stackSentinel to handle only a small area of a single burst？

isce-framework / isce2

InSAR Scientific Computing Environment version 2

Other

508 stars 250 forks source link

How to use stackSentinel to handle only a small area of a single burst？ #807

Closed ZGHHGZ closed 7 months ago

ZGHHGZ commented 9 months ago

The smallest unit of stackSentinel for data alignment is the burst, but my study area is often very small, maybe 1/100 of a single burst, but the program needs to align the entire burst, which usually takes a long time and generates huge files, is there a modification to align only a certain area in the burst? I looked at the program code in topsStack and found that it reads the width and length of the entire burst. It is not clear to me from which file this reading is taken and whether I can modify the data that the program needs to read to achieve my purpose. I tried changing the length and width values within the .vrt file or .xml, but it still aligns the entire burst!

comma-never-coma commented 9 months ago

A feasible method is to crop combined interferograms and geom files, thus u can obtain smaller filtered and unwrapped interferograms.

Nirdesh31 commented 8 months ago

Hi @comma-never-coma, Can you elaborate what do you mean by 'crop combined interferograms and geom files' here? And what files are you referring here?

ZGHHGZ commented 8 months ago

Hi @comma-never-coma, Can you elaborate what do you mean by 'crop combined interferograms and geom files' here? And what files are you referring here?

https://github.com/isce-framework/isce2/blob/main/contrib/timeseries/prepStackToStaMPS/bin/run_SLCcropStack.csh

This might help you. It can be used to clip the registered SLC.

comma-never-coma commented 8 months ago

This refers to the files located under the 'merged/interferograms' and 'merged/geom_reference' directories. In fact, you can also see this section of cropped data in the sample data Sentinel-1 on Fernandina with ISCE from mintpy. BTW, I haven't tried to crop SLC with 'run_SLCcropStack.csh'. I just use cropped merged interferograms and geom_reference files for filtering and unwrapping later. Generating interferograms with cropped SLC looks better.

ZGHHGZ commented 7 months ago

Hi @comma-never-coma, Can you elaborate what do you mean by 'crop combined interferograms and geom files' here? And what files are you referring here?

This refers to the files located under the 'merged/interferograms' and 'merged/geom_reference' directories. In fact, you can also see this section of cropped data in the sample data Sentinel-1 on Fernandina with ISCE from mintpy. BTW, I haven't tried to crop SLC with 'run_SLCcropStack.csh'. I just use cropped merged interferograms and geom_reference files for filtering and unwrapping later. Generating interferograms with cropped SLC looks better.

Based on isce2, I have implemented the process of only downloading a single burst for processing, first cropping and then registering. For hundreds of images, if your research area is very small, it can take 1-2 hours and tens of GB of space to complete all the processes from data downloading to SBAS-InSAR processing. https://github.com/ZGHHGZ/ML-TS-InSAR

HN-Jiang commented 1 month ago

Hi @comma-never-coma A quick question about cropped merged interferograms and geom_reference files. I use stackSentinel.py -W interferogram and get 01-11 step in run_files. So if I want to use run_SLCcropStack.csh to crop interferograms to save times for unwrapping. Can I directly copy this script into my run_files folder and run it?

ZGHHGZ commented 1 month ago

Hi @comma-never-coma A quick question about cropped merged interferograms and geom_reference files. I use stackSentinel.py -W interferogram and get 01-11 step in run_files. So if I want to use run_SLCcropStack.csh to crop interferograms to save times for unwrapping. Can I directly copy this script into my run_files folder and run it?

‘isce2/contrib/timeseries/prepStackToStaMPS/bin The SLC cropping function in ‘/run_SLCcropStack.csh’ is specifically designed to prepare data for the StaMPS software and is not adapted to the stackSentinel interferometry process. isce's stackSentinel script interferometry utilises unmerged SLC data. If you want to make the cropped merged SLC available for interferometry, you may need to modify some of the stackSentinel code. 006de0c7cc4d68685eab37374591abb5

Of course, you can also try the following method for fast processing of time series InSAR. https://github.com/ZGHHGZ/ESIS

https://github.com/ZGHHGZ/ML-TS-InSAR