INSAR 处理 ROI_PAC应用指南 (2).doc

ROI_PAC应用指南1 软件下载及安装1.1 软件下载及在线帮助Download the software from: Open Channel FoundationTo get help from: ROI_PAC Wiki1.2 ROI_PAC及其一些必要软件的安装ROI_PAC安装参考：/installation，可根据提示按步骤进行安装；或阅读安装目录下文档AAREADME_BUILD_ROIPAC。下面的安装步骤介绍单独安装各个组件的方法(与网站上的方法略有不同)。1.2.1 FFTWFFTW是由MIT的Matteo Frigo和Steven G Johnson开发的,可用来进行快速傅立叶变换,目前最新版本v3.2.2, 相关网址 /解压FFTW安装包： tar -xzf fftw-3.2.2.tar.gz进入 FFTW 代码所在目录，编译并安装FFTW软件到指定目录(设为usr/local/fftw)， 依次执行： cd fftw-3.2.2 ./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的安装环境，然后利用软件中自带的安装脚本进行安装： 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 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_PACexport INT_BIN=$ROIPAC/INT_BINexport INT_SCR=$ROIPAC/INT_SCRexport MY_SCR=$ROIPAC/MY_SCRexport MY_BIN=$ROIPAC/MY_BINexport FFTW_LIB=/usr/local/libexport FFTW_INC=/usr/local/includeexport PATH=$PATH:$INT_BIN:$INT_SCRexport SAR_PRC_DIR=$ROIPAC/PRCexport SAR_ODR_DIR=$ROIPAC/ODR启动软件前运行命令source设置环境变量：source /usr/local/ROI_PAC_3_0_1/ROI_PAC/SAR_CONFIG1.2.3 SNAPHU安装解压SNAPHU安装包，进入snaphu-v1.4.2/src目录，输入make，编译成二进制可执行文件，并在/ROI_PAC/MY_BIN目录下创建snaphu程序的快捷方式：cdsnaphu-v1.4.2/src makeln snaphu-v1.4.2/src /usr/local/ROI_PAC_3_0_1/ROI_PAC/MY_BIN1.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: 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 ROI_PAC File Structure File FormatData Type*.slccpxBSQcomplex*.intcpxBSQcomplex*.corrmgBILfloating point*.unwrmgBILfloating point*.hgtrmgBILfloating point*.ampcpxBSQcomplexcpx = complex values, alternating real and imaginary pixels rmg = alternating lines of amplitude and 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格式)包括三个文件，例如：data: IMAGERY1993011018252739T1Of3leader: SARLEADER1993011018252739T1Of1volume directory: VDF1993011018252739T1ROI_PAC命令示例：(当前目录为TESTDIR/930110)make_raw.pl PRC SARLEADER1995052318253409T1Of1 930110For 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_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-ALPSRP111820640-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 DoItFrom EndItAt (Start-End point detail)首先在TESTDIR文件夹下创建并编辑输入信息文件c (parameter detail)ROI_PAC命令示例(当前目录为TESTDIR)：process_2crawslcsprocess_2cslcs full_res运行过程中的问题参见：/eas/PeoplePlaces/Faculty/matt/2a.html2.3.3 ROI_PAC处理流程图InSAR Processing (Parameter Details)Create a file named c that contains 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 output filesIn addition, the following optional parameters can be specified: Optional ParametersDescriptionOptional ValuesDefault ValueSimDirDirectory containing topographic simulation output filesSIMDEMFull path to Digital Elevation Model (DEM) used to remove topographyNULLGeoDirDirectory containing tectonic simulation output filesGEOFilterStrengthStrength of the filter used in the processing.0-10.75UnwrappedThresholdDetermines the extent of unwrapping. Connected regions with a magnitude greater than the threshold as found in the file *.msk are included in the unwrapping process.0-10.1OrbitTypeType of orbit file, either precise orbits (PRC) or Delft orbits (ODR)PRC, ODRODRBaselineTypeGiven to discriminate baseline estimation techniqueOrbitTypeRlooks_simSpecifies the number of looks to take when averaging the raw SAR data.4Rlooks_intRlooks_simRlooks_unwRlooks_simRlooks_sml16pixel_ratiopixel size ratio of azimuth/range5Alooks_smlDefined as (pixel_ratio X Rlooks_sml)80before_z_ext1500after_z_ext1500near_rng_ext700far_rng_ext0usergivendop1User defined doppler value for scene 10usergivendop2User defined doppler value for scene 20unw_seedxx coordinate where the phase unwrapping begins. The default places the seed in the middle of the image.-9999unw_seedyy coordinate where the phase unwrapping begins. The default places the seed in the middle of the image.-9999x_startSpecifies the 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.01y_startSpecifies the azimuth offset in pixels between the two SAR scenes. The default value signals the processing to use the estimated azimuth offset calculated from the a priori orbit information.0.01flatteningUsed 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 to re-estimate the baseline information using a DEM and/or other model information.orbit, topoorbitThreshold_magDefines 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 small. (0=unwrap everywhere)0-0.000055.0e-5Threshold_ph_grd0sigma_thresh1.0smooth_widthSpecifies the size of the filter window in range and azimuth used when creating the phase gradient mask5slope_width5concurrent_roiAllows SAR processor to work on both SAR scenes simultaneously.yes, nonomappingdem_baseddem_basedcleanupRemove intermediate files when finishedyes, nonoCO_MODELSpecifies a Coseismic model used in flattening. Make sure to set flattening=topo. Read the instructions for INTER_MODEL for details about naming input files.NULLINTER_MODELSpecifies an Interseismic model used in flattening. Make sure to set flattening=topo for this option. Provide three data files for each component 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.NULLFilt_methodSpecify type of filter usedpsfilt, adapt_filtpsfiltunw_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 SAR processor orbbaseCalculate baselines slcsCreate single-l