suesun0890 / roi_pac-sentinel1 Goto Github PK
View Code? Open in Web Editor NEWThis project forked from raphaelgrandin/roi_pac-sentinel1
InSAR processing of Sentinel-1 using ROI_PAC
License: Other
This project forked from raphaelgrandin/roi_pac-sentinel1
InSAR processing of Sentinel-1 using ROI_PAC
License: Other
#################################################### ### Sentinel-1 pre-processor for ROI_PAC ### Raphael Grandin (IPGP) -- April 2016 ### [email protected] ### http://www.ipgp.fr/~grandin #################################################### ### ### README FILE ### ### Content : ### - installation guidelines ### - how it works ### - bibliography ### - license ### ### #################################################### ---------------------------------------------------------------- ---------------------------------------------------------------- INSTALLATION GUIDELINES 1/ Install the following software / packages : - ROI_PAC Version 3.0.1 - python 2.6 => will be called thru command "python" - a few python libraries should be installed, including: mpmath, gdal - optionally : snaphu (http://web.stanford.edu/group/radar/softwareandlinks/sw/snaphu/) - optionally : generic mapping tools (http://gmt.soest.hawaii.edu/) 2 / Modify a few ROI_PAC routines and re-compile ROI_PAC. The modified routines are located in directory "ROI_PAC_changes" The modifications are there to enable ROI_PAC to deal with wide Sentinel-1 images and to handle correctly ESA orbit state vectors. 3 / Copy the content of "MY_SCR" directory to the user's own MY_SCR dir. This directory contains a few perl scripts that are occasionally called during the process. ---------------------------------------------------------------- ---------------------------------------------------------------- WARNING ROI_PAC is a copyrighted software that requires a license. Licenses are available at no charge for non-commercial use from the Open Channel Foundation. Read the license terms. ---------------------------------------------------------------- ---------------------------------------------------------------- HOW IT WORKS The principle of Sentinel-1 InSAR processing is described in Grandin (2015). A prerequisite is that the user downloads Sentinel-1 data in SLC format. Then, the user fills a parameter file that tells the package what data should be used, and where to find the data. The parameter file should look like this : ***************************************** BEGIN ******************************************* PATHDIR /home/rgrandin/S1-ROI_PAC_V2.4 WORKINGDIR /media/Quadra_4/Nepal_SD_test/Frame3/20150409_20150503 DIR_ARCHIVE /media/Quadra_3/Nepal2015/Sentinel/data DEM /media/Quadra_3/Nepal2015/DEM/srtmV4_nepal_corrpatch_shift_destrip_crop_Nepal2015_Frame3_resamp4x.dem DEM_low /media/Quadra_3/Nepal2015/DEM/srtmV4_nepal_corrpatch_shift_destrip_crop_Nepal2015_Frame3_resamp4x_4rlks.dem LABEL_ante 20150409 LABEL_post 20150503 DIR_IMG_ante S1A_IW_SLC__1SSV_20150409T122127_20150409T122154_005407_006E01_6F00.SAFE DIR_IMG_ante S1A_IW_SLC__1SSV_20150409T122152_20150409T122224_005407_006E01_0098.SAFE DIR_IMG_post S1A_IW_SLC__1SSV_20150503T122118_20150503T122145_005757_007649_5EE7.SAFE DIR_IMG_post S1A_IW_SLC__1SSV_20150503T122143_20150503T122211_005757_007649_F885.SAFE DIR_IMG_post S1A_IW_SLC__1SSV_20150503T122208_20150503T122235_005757_007649_A437.SAFE LOOKS_RANGE 12 LOOKS_AZIMUTH 4 ***************************************** END ******************************************* Explanation : ************ MANDATORY ************* PATHDIR = path to the directory where the Sentinel-1 ROI_PAC package is located WORKINGDIR = path to the working directory, i.e. where the processing will take place DIR_ARCHIVE = path to the directory containing Sentinel-1 SAFE directories. DEM = path and name of a dense DEM (resolution should be ~ 20 m or better – oversampled SRTM is OK) DEM_low = path and name of a coarser DEM (it is OK to multilook the above DEM by a factor 4) LABEL_ante = name of master date, in YYYYMMDD format LABEL_post = name of slave date, in YYYYMMDD format DIR_IMG_ante = name of Sentinel-1 data directory for the master image DIR_IMG_ante = another directory to append and create a long strip DIR_IMG_ante = yet another one... they should be provided in chronological order. DIR_IMG_post = same, but for the slave image ************ OPTIONAL ************* LOOKS_RANGE = (integer) number of looks in range (12 is OK) LOOKS_AZIMUTH = (integer) number of looks in azimuth (4 is OK) SKIP_BEG_ante = (integer) skip this number of bursts at the top of first master image SKIP_END_ante = (integer) skip this number of bursts at the bottom of last master image SKIP_BEG_post = (integer) skip this number of bursts at the top of first slave image SKIP_END_post = (integer) skip this number of bursts at the bottom of last slave image SPECTRAL_DIV = (yes or no) run (yes) or skip (no) spectral diversity. Default: yes. Finally, the user copies the file named "S1-ROI_PAC_00_batch.sh" (located in directory "scripts") into his working directory. This script is run from command line, and takes the parameter file as the single (required) argument. For example : ./S1-ROI_PAC_00_batch.sh topsar_param.in This script runs the successive processing steps. Each step corresponds to a script, numbered from 01 to 09 in directory "scripts". All the steps are run successively for the three sub-swaths. By default, the process is run using polarization VV, but this can be changed by editing "S1-ROI_PAC_00_batch.sh". The user is advised to run the steps manually one after the other in the first place in order to become familiar with the procesing. The scripts take the parameter file as first argument, followed by sub-swath name and polarization. For instance : ./S1-ROI_PAC_01_stitch.csh topsar_param.in iw2 vv ./S1-ROI_PAC_02_baseline.csh topsar_param.in iw2 vv ./S1-ROI_PAC_03_correl.csh topsar_param.in iw2 vv ...... etc ...... The successive steps are the following : STEP 1 : S1-ROI_PAC_01_stitch.csh 01 / From a list of provided SAFE directories, order them in chronological order 02 / Decode the metadata from .xml files in SAFE directories. If necessary, the zipped SAFE directories are unzipped. 03 / Compare length of master / slave images and only keep bursts in the intersection of the two acquisitions 04 / Download orbits (restituted and precise) from ESA QC web site (S1-ROI_PAC_01b_download_orbits.sh) 05 / Correct azimuth phase drift burst-by-burst to bring azimuth spectrum to baseband, for the master and slave SLCs 06 / Stitch the bursts to produce continuous images 07 / Save two alternative SLCs for the master image containing forward-looking and backward-looking phases in burst overlaps. STEP 2 : S1-ROI_PAC_02_baseline.csh 08 / Start ROI_PAC processing to create file tree structure 09 / Compute baselines STEP 3 : S1-ROI_PAC_03_correl.csh 10 / Compute dense pixel offsets between master and slave images STEP 4 : S1-ROI_PAC_04_sync.csh 11 / Re-run azimuth phase drift correction for slave image using time offset deduded from pixel offsets. 12 / Also save two alternative SLCs for the slave image containing forward-looking and backward-looking phases in burst overlaps. STEP 5 : S1-ROI_PAC_05_interf.csh 13 / Run ROI_PAC from SLCs to interferogram corrected from simulation STEP 5b (optional) : S1-ROI_PAC_05_interf.csh 13b / Re-run interferogram calculation only STEP 6i : S1-ROI_PAC_06_sd.csh 14 / Compute double-difference interferogram for spectral diversity STEP 7 : S1-ROI_PAC_07_sd_finish.csh 15 / Re-run azimuth phase drift correction for slave image using time offset deduded from pixel offsets and refinement from spectral diversity 16 / Re-calculate interferogram STEP 7b (optional) : S1-ROI_PAC_07b_sd_refine.csh 16b / Update across-track burst overlap interferogram to account for correction deduced from spectral diversity. STEP 8 : S1-ROI_PAC_08_unwrap.csh 17 / Filter and unwrap interferogram STEP 9 : S1-ROI_PAC_09_done.csh 18 / Geocode unwrapped interferogram using coarse DEM STEP 10 (unstable) : S1-ROI_PAC_10_geocode_beta.csh 19 / Stich the geocoded sub-swaths together to produce seamless interogram. ---------------------------------------------------------------- ---------------------------------------------------------------- REFERENCES Grandin, R., (2015). "Interferometric processing of SLC Sentinel-1 TOPS data." In Proceedings of the European Space Agency, SP-371. FRINGE 2015, 23-27 March 2015, Frascati (Rome), Italy. Grandin, R., E. Klein, M. Métois, and C. Vigny (2016). "3D displacement field of the 2015 Mw8.3 Illapel earthquake (Chile)" "from across- and along-track Sentinel-1 TOPS interferometry" Geophys. Res. Lett., 43, doi:10.1002/2016GL067954. PDFs can be found in directory "bibliography" ---------------------------------------------------------------- ---------------------------------------------------------------- LICENSE * Copyright (C) 2016 R.GRANDIN * [email protected] * "Sentinel-1 pre-processor for ROI_PAC". * "Sentinel-1 pre-processor for ROI_PAC" is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version. * "Sentinel-1 pre-processor for ROI_PAC" is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. * You should have received a copy of the GNU General Public License along with "Sentinel-1 pre-processor for ROI_PAC". If not, see <http://www.gnu.org/licenses/>.
A declarative, efficient, and flexible JavaScript library for building user interfaces.
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
An Open Source Machine Learning Framework for Everyone
The Web framework for perfectionists with deadlines.
A PHP framework for web artisans
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
Some thing interesting about web. New door for the world.
A server is a program made to process requests and deliver data to clients.
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
Some thing interesting about visualization, use data art
Some thing interesting about game, make everyone happy.
We are working to build community through open source technology. NB: members must have two-factor auth.
Open source projects and samples from Microsoft.
Google ❤️ Open Source for everyone.
Alibaba Open Source for everyone
Data-Driven Documents codes.
China tencent open source team.