INSAR处理ROIPAC应用指南

1、ROI_PAC应用指南1 软件下载及安装1.1 软件下载及在线帮助Download the software from: Open Channel FoundationTo get help from: ROI_PAC Wiki1.2 ROI_PAC及其一些必要软件的安装ROI_PAC安装参考：，可根据提示按步骤进行安装；或阅读安装目录下文档AAREADME_BUILD_ROIPAC。下面的安装步骤介绍单独安装各个组件的方法(与网站上的方法略有不同)。1.2.1 FFTWFFTW是由MIT的Matteo Frigo和Steven G Johnson开发的,可用来进行快速傅立叶变换,目前最新版本

2、v3.2.2, 相关网址 解压FFTW安装包： tar -xzf fftw-3.2.2.tar.gz进入 FFTW 代码所在目录，编译并安装FFTW软件到指定目录(设为usr/local/fftw)， 依次执行： ./configure -prefix=/usr/local/fftw -enable-float make make install1.2.2 ROI_PAC安装解压ROI_PAC安装包，并将解压后文件夹剪切到/usr/local下： tar -xzf ROI_PAC_3_0_1.tgz进入ROI_PAC代码所在目录(ROI_PAC_3_0_1/ROI_PAC)，首先设置FFTW的

3、安装环境，然后利用软件中自带的安装脚本进行安装： export FFTW_LIB_DIR=/usr/local/fftw/lib export FFTW_INC_DIR=/usr/local/fftw/include ./contrib./multibuild.sh此时，ROI_PAC的可执行代码在ROI_PAC_3_0_1/ROI_PAC/multibuild*/installs/defaults/bin下，在ROI_PAC_3_0_1/ROI_PAC下创建INT_BIN文件夹并将可执行代码拷贝到该文件夹： mkdir ROI_PAC_3_0_1/ROI_PAC/INT_BINexport

4、FFTW_INC_DIR=/usr/local/fftw/includecd ROI_PAC_3_0_1/ROI_PAC/multibuild*/installs/defaults/bincp * ROI_PAC_3_0_1/ROI_PAC/INT_BIN/将ROI_PAC_3_0_1/ROI_PAC/DOC/SAR_CONFIG拷贝至主目录(/usr/local/ ROI_PAC_3_0_1/ROI_PAC)，并编辑以设置环境变量，SAR_CONFIG文件中设置的环境变量如下：ROIPAC="/usr/local/ ROI_PAC_3_0_1/ROI_PAC"export

5、 INT_BIN="$ROIPAC/INT_BIN"export INT_SCR="$ROIPAC/INT_SCR"export MY_SCR="$ROIPAC/MY_SCR"export MY_BIN="$ROIPAC/MY_BIN"export FFTW_LIB="/usr/local/lib"export FFTW_INC="/usr/local/include"export PATH="$PATH:$INT_BIN:$INT_SCR"export

6、SAR_PRC_DIR="$ROIPAC/PRC"export SAR_ODR_DIR="$ROIPAC/ODR"启动软件前运行命令source设置环境变量：source /usr/local/ROI_PAC_3_0_1/ROI_PAC/SAR_CONFIG1.2.3 SNAPHU安装解压SNAPHU安装包，进入目录，输入make，编译成二进制可执行文件，并在/ROI_PAC/MY_BIN目录下创建snaphu程序的快捷方式：cdsnaphu-v1.4.2/src makeln /usr/local/ROI_PAC_3_0_1/ROI_PAC/MY_BI

7、N1.2.4 MDX安装解压MDX安装包，修改mdx_176_75目录下的Makemdx_g77文件，对于fedora来说，将第9行和第10行修改为：XINC=-I/usr/include/X11 I/usr/include/XmXLIB=-lXm L/usr/lib lXt lX11Xm库需事先安装grace软件；在终端输入g77 v，如果没有该命令，还需安装gcc3.4.6。运行以下命令完成安装：make f Makemdx_g772 ROI_PAC二轨DInSAR处理步骤2.1 前期准备：SAR Data: ERS-1/2, Envisat ASAR, ALOS PALSAROrbit:

8、 PRC(DPAF)/ODR(Delft)/HDR(header) for ERS DOR(ESA DORIS)/ODR(Delft) for Envisat raw data; HDR for int HDR for ALOSDEM: SRTM DEM 90m or 30m Resolution创建工作目录(如/TESTDIR)，在该目录下创建主副影像目录(以成像日期命名，如/TESTDIR/070225)，若需要单独的轨道数据，将所需的轨道文件按一定的格式存储于所对应的文件夹(SAR_CONFIG中设置)。2.2 文件类型： ROI_PAC File Suffix &#

9、160; ROI_PAC File Structure   File Format  Data Type*.slccpxBSQcomplex*.intcpxBSQcomplex*.corrmgBILfloating point*.unwrmgBILfloating point*.hgtrmgBILfloating point*.ampcpxBSQcomplexcpx = complex values, alternating real and imaginary pixels rmg = alternating lines of amplitude an

10、d phase BSQ = band sequentialBIL = band interleave2.3 数据处理 2.3.1 数据整理，生成RAW数据For ERS1/2:make_raw.plUsage: make_raw.pl orbit_type leader_file dateorbit type can be: HDR (header file), PRC (DPAF) or ODR (Delft)将数据放入以成像日期命名的文件夹。每景影像(CEOS格式)包括三个文件，例如：volume directory: ROI_PAC命令示例：(当前目录为TESTDIR/930110)ma

11、ke_raw.pl PRC S930110For Envisat:make_raw_envi.plUsage:make_raw_envi.pl asa_file_root orbit_type outnameOrbit type can be: DOR ( ESA DORIS) or ODR (Delft) (use “Orbit Type=HDR” in processing)将数据放入以成像日期命名的文件夹。每景影像包括一个文件：ASA_IM_0CNPDE20040623_031040_000000752028_00018_12096_0013.N1 (ESA name)orENV1_4_

12、442_2871_2889_18031.baq (WInSAR name)数据命名规则：For ALOS:make_raw_alos.pl (modified by Li, pre-process program in MY_BIN)Usage: make_raw_alos.pl alos_file_root outname fbd2fbs_conversionOrbit type: HDR将数据放入以成像日期命名的文件夹。每景影像包括三个文件(FBS)或四个文件(FBD)：data: IMG-HH-ALPSRP111820640-H1.0_A(FBS&FBD)IMG-HV-ALPSR

13、P111820640-H1.0_A(FBD)leader: LED-ALPSRP111820640-H1.0_Avolume directory (not used): VOL-ALPSRP111820640-H1.0_A数据命名规则：ROI_PAC命令示例: (当前目录为TESTDIR/070225)make_raw_alos.pl IMG-HH-ALPmake_raw_alos.pl IMG-HH-ALP 070225 FBD2FBS2.3.2由RAW数据到SLC，由SLC生成干涉图这个过程由process_2pass.pl完成：Usage: process_2pass infile Do

14、ItFrom EndItAt (Start-End point detail)首先在TESTDIR文件夹下创建并编辑输入信息文件c (parameter detail)ROI_PAC命令示例(当前目录为TESTDIR)：process_2crawslcsprocess_2cslcs full_res运行过程中的问题参见：2.3.3 ROI_PAC处理流程图InSAR Processing (Parameter Details)Create a file named "c" that cont

15、ains the input parameters for the main processing script.  The file must define the following parameters: Required ParametersDescriptionSarDir1Directory containing raw SAR data for first sceneSarDir2Directory containing raw SAR data for second sceneIntDirDirectory containing interferometric out

16、put filesIn addition, the following optional parameters can be specified: Optional ParametersDescriptionOptional ValuesDefault Value SimDirDirectory containing topographic simulation output files SIM DEM Full path to Digital Elevation Model (DEM) used to remove topography NU

17、LL GeoDir Directory containing tectonic simulation output files GEO FilterStrength Strength of the filter used in the processing. 0-10.75 UnwrappedThreshold Determines the extent of unwrapping.   Connected regions with a magnitude greater than the th

18、reshold as found in the file *.msk are included in the unwrapping process.0-10.1 OrbitType Type of orbit file, either precise orbits (PRC) or Delft orbits (ODR)PRC, ODRODR BaselineType Given to discriminate baseline estimation technique OrbitType Rlooks_sim Specifi

19、es the number of looks to take when averaging the raw SAR data. 4 Rlooks_int  Rlooks_sim Rlooks_unwRlooks_sim Rlooks_sml16 pixel_ratio pixel size ratio of azimuth/range 5 Alooks_sml Defined as (pixel_ratio X Rlooks_sml)  80 before

20、_z_ext   1500 after_z_ext   1500 near_rng_ext   700 far_rng_ext   0 usergivendop1 User defined doppler value for scene 1 0 usergivendop2 User defined doppler value for scene 2  0 unw_se

21、edx x coordinate where the phase unwrapping begins. The default places the seed in the middle of the image.  -9999 unw_seedy y coordinate where the phase unwrapping begins.  The default places the seed in the middle of the image. -9999 x_start Specifies t

22、he range offset in pixels between the two SAR scenes.  The default value signals the processing to use the estimated range offset calculated from the a priori orbit information. 0.01 y_start Specifies the azimuth offset in pixels between the two SAR scenes.  The default valu

23、e signals the processing to use the estimated azimuth offset calculated from the a priori orbit information. 0.01 flattening Used to specify the method of flattening. Set this variable to "orbit" if you do not have a DEM or other model simulation.  Set to "topo&quo

24、t; to re-estimate the baseline information using a DEM and/or other model information.orbit, topoorbit Threshold_mag Defines an amplitude threshold when producing a mask during the unwrapping.  This helps to reduce side effects on the edges of the radar scene where amplitudes are smal

25、l. (0=unwrap everywhere)  0-0.000055.0e-5 Threshold_ph_grd   0 sigma_thresh   1.0 smooth_width  Specifies the size of the filter window in range and azimuth used when creating the phase gradient mask 5 slope_width  

26、 5 concurrent_roi Allows SAR processor to work on both SAR scenes simultaneously.yes, nono mapping  dem_baseddem_based cleanup Remove intermediate files when finishedyes, nono CO_MODEL Specifies a Coseismic model used in flattening. Make sure to set

27、"flattening=topo".  Read the instructions for INTER_MODEL for details about naming input files. NULL INTER_MODEL Specifies an Interseismic model used in flattening. Make sure to set "flattening=topo" for this option.  Provide three data files for each com

28、ponent of displacement along with three resource files.  ROI_PAC will expect the files to be named XE, XE.rsc, XN, XN.rsc, XU, and XU.rsc for the east, north, and vertical components, respectively, where the "X" can be a user defined prefix specified by the INTER_MODEL parameter. 

29、;NULL Filt_method Specify type of filter usedpsfilt, adapt_filtpsfilt unw_methodSpecify unwrapping method.  Bridge.pl is only compatible with the "old" unwrapping method. old, icuoldExample:   SarDir1=920610 SarDir2=980509 IntDir=920610-980509 DEM=/a/computer/data/sar/DEM/SF/SF.dem flattening=topo REFERENCE: /das/ROIPAC/proc.htmlInSAR Processing (Start-End Point Details)START/END pointCommentsraw*Always start here after running make_raw.plroi_prepCreate input files for S