doris使用InSAR处理过程及参数说明Word格式文档下载.docx

1、运行第三步：run s3运行第四步：run s4运行第五步：run s5运行第六步：run s6运行第七步：run s7运行第八步：run s92，参数解说及设置。2.1，第一步：读取主影像信息。（1）需要处理步骤。必需：m_readfiles（读取主影像头文件信息）.m_crop（读取主影像数据体信息）.不需：m_ovs（主影像过采样步骤）：作干涉处理时，一般不运行。可选：m_prorbits（读取主影像精轨文件）:若要使用提供的精密轨道数据时，就执行该行，否则直接C掉。（2）小步骤参数说明及设置。#INSAR通用参数设置，每步都有的。#名称：LOGFILE用处： 用以记载DORIS处理过程

2、中所有信息。文件名：doris.log.M_RESFILE 用以记载主影像文件所有信息。master.resS_RESFILE 用以记载辅影像文件所有信息。slave.resI_RESFILE 用以记载主辅影像配准及以后步骤处理的所有信息。master_slave.res分步骤：M_readfiles（读取主影像头文件信息）M_IN_METHOD使用的原始数据的传感器类型。有（ERS（CEOS格式），ASAR（ENVISAT格式），JERS（CEOS格式），RSAT（RSAT格式），ATLANTIS（RSAT格式）。说明：必需针对CEOS和RSAT格式的参数文件M_IN_VOL SLC的VOL

3、UME文件M_IN_LEA SLC的LEADER文件M_IN_DAT SLC的DAT文件M_IN_NULL SLC的null文件，默认为dummy针对ENVISAT格式的参数文件 SLC的DAT文件（ASAR只有一个文件）M_crops（读取主影像数据体信息）M_CROP_IN输入数据体的路径。必需 .dat文件。M_CROP_OUT输出的raw数据体路径。M_DBOW要处理的数据区的大小，直接按影像大小来计数。可选 ，C掉时表示默认为整块影像大小。M_DBOW_GEO要处理的数据区的大小，按按经纬度来计数。一般不用这个，它与M_DBOW不能同选。M_ovr（主影像过采样，针对PSINSAR处

4、理的）M_porbits（读取精密轨道信息）M_ORBDIR精密轨道文件路径。针对精轨处理才给定路径。2.2，第二步：读辅影像信息。设置基本与第一步类同。这里不作详细介绍。2.3，第三步：影像配准。（1）需要处理步骤coarseorb 一级轨道级配准。coarsecorr 二级像素级粗配准m_filtazi 主影像方位向滤波 s_filtazi 辅影像方位向滤波filtrange 距离向滤波fine 三级子像素精配准（2）小步骤参数设置说明coarseorb（一级轨道级配准）coarsecorr（二级像素级粗配准）CC_METHOD magfft |magspace相关系数计算的方法，magf

5、ft 是先使用fft变换到频率域，然后计算相关系数。推荐使用magspace 是直接在空间域计算相关系数。CC_IN_POS filenameASCII 文件，记录主影像上待匹配的点的坐标。（不采用）CC_NWIN 21像素级配准划分窗口的数目，大于。当CC_IN_POS被选择时，就不需要设置。CC_WINSIZE 256 256匹配窗口的大小，处理的时候被改变为奇数。CC_ACC 64 32在空间域需要用户设置，在方位向和距离向的搜索偏移量。在频率域则自动设置为 CC_WINSIZE的一半。CC_INITOFF 0 0| “orbit”初始的偏移量。“orbot” 直接读入使用轨道计算的偏移

6、量。0 0 用户之际给点的初始偏移量m_filtazi（主影像方位向滤波）s_filtazi（辅影像方位向滤波）AF-BLOCKSIZE 1024FFT变换中每次距离向中的行数。一般是越大越好。AF_OVERLAPHANMMING 0.75光谱滤波的权值。滤波器的alpha 系数AF_OUT_MASTER master.afiltered滤波后的主影像文件名AF_OUT_SLAVE slave.afileered滤波后的辅影像文件名AF_OUT_FORMAT cr4输出文件数据的格式， cr4 ：complex real4 ci2 ：complex shortsfiltrange（距离向滤波）

7、fine （三级子像素精配准）FC_METHOD magfft | magspace | oversample精配准计算方法。Magfft： 幅度的fft变换后计算相关系数法Magspace ：幅度直接计算系数法Overspace ：最好的方法，过采样处理。FC_NWIN 601子像素级配准时划分窗口的数目，比较大才好。FC-IN_POS file nameFC_WINSIZE 64 64匹配窗口的大小。FC_ACC 8 8FC_INITOFF 0 0 |COARSECORR同CC_INITOFF ，COARSECORR表明使用像素级配准后的偏移量。FC_OSFACTOR 16过采样的内插倍数

8、。2.4，第四步：影像重采样。coregpm 估算配准多项式系数。resample 辅影像重采样。coregpm（估算配准多项式系数）CPM_THRESHOLD 0.4经配准后生成的大量匹配点，当相关系数大于CPM_THRESHOLD，才参与平差运算。阈值的大小与FILE 中窗口的大小有关。小窗口情况下阈值趋于1.0 ，阈值越大越好。对于64的窗口，阈值选择0.2即可。CPM_DEGREE 2多项式的阶数。2阶即可。CPM_WEIGHT BAMLER | NONE| LINEAR | QUADRATIC最小二乘中权矩阵方法选择。BAMLER：权矩阵为相关系数的2倍。LINEAR：权矩阵为相关系

9、数的1倍.CPM_MAXITER 20最小二乘运算的最大迭代次数。resample （辅影像重采样）RS_METHOD RECT | TRI |CC4P| CC6P| TS6P | TS8P | 。重采样的内插核.默认RC6PRS_OUT_FILE s_resampled.raw重采样后生成的文件名。RS_OUT_FORMAT CR4 | CI2输出文件数据类型。CR4： comlex_real 4；CI2：complex_shortRS_DBOW 格式（a，a，a，a）重采样的区域大小。C掉表示整副采样。（这里的区域坐标是针对主影像而言的）。RS_SHIFTAZI ONOFFON，表示采样时

10、，采样核随着多普勒中心频率变化而变化，OFF，则表示不随便变化而采用固定采样核。一般情况下，如果多普勒中心频率变化很小，可以OFF掉，默认为ON。filtrange （距离向滤波）以下转自：2.5，第五步：干涉图生成。interfero 干涉图生成。coherence生成相干系数图括号里也是必须的，但必须是四选二，或选一和二，或选三和四，不能全选。comprefpha 用轨道来计算平地相位Subtrrefpha 减去平地相位comprefdem 利用提供的粗DEM来反算求得平地相位subtrrefdem 减去由粗DEM反算得到的平地相位interfero （干涉图生成）INT_OUT_CINT

11、（或INT_OUT_INT）生成的干涉图文件（INT_OUT_CINT为复数float型， INT_OUT_INT为float型，）。必需，但二选一。INT_METHOD生成的干涉图的方法，有（old，oversample）两种，必需，一般选oversample。INT_MULTILOOK多视处理系数。必需。如果没有多视要求，则默认（1 1），如果在这一步参考相位没有去除，则多视参数不应太大，一般默认为：（10 2）或者（5 1）、（8 2）.comprefpha （用轨道来计算平地相位）FE_METHOD采用计算平地相位的方法，只有精轨（porbits），如果C掉此行，则为粗轨。可选。FE_

12、GEGREE位置与平地相位之间的多项式系数，一般二次必需，默认5。FE_NPOINTS最小二乘法来计算平地相位所用来拟合的点数必需，默认501。Subtrrefpha （减去平地相位）SRP_METHOD采用减去平地相位的方法，有（polynomial（多项式）和exact）两种多项式方式表示减去采用前一步采用拟合的多项式得到的平地相位，exact表示采用三个核心方程逐点计算平地相位并减去。默认为polynomial。SRP_OUT_CINT去平地效应后的干涉图（复数float型）。cint.minrefpha.raw。SRP_MULTILOOK多视系数必需，默认（2 2），注意，这里要与前一

13、步的多视系数相乘等于真实多视系数。coherence（成相干系数图）COH_OUT_CCOH（或COH_OUT_COH）生成的相干图文件（COH_OUT_CCOH为复数float型，COH_OUT_COH为float型，），二者只能选一。必需，但二选一，给定其路径。COH_MULTILOOK多视系数。必需，默认（10 2）COH_WINSIZE用于计算相干系数的窗口。必需，（15 15）COH_METHOD用于计算相干系数的方法可选 LD；comprefdem （利用提供的粗DEM来反算求得平地相位）CRD_METHOD求取平地相位方法，有（NEAREST和TRILINEAR）两种。由于是粗D

14、EM，与影像的分辨率不一，需要过采样到同一分辩率下，过采样方式两种，一种是最近邻元方法（NEAREST），一种是内插方式，必需，默认为TRILINEARCRD_IN_FORMAT输入粗DEM的格式，有（I2（short int），R4（float），R8（double），I2_BEGENDIAN）,Gtopo30 和SRTM都是I2型。可选，C掉时系统默认为I2.CRD_IN_DEM输入DEM的路径名称。是已经拼好的DEM。CRD_IN_SIZE输入DEM的大小，分经纬向。按（纬度，经度）输入CRD_IN_DELTA输入DEM的间隔，一般经纬向间隔是一样的。必需，SRTM的间隔为 0.0083

15、3333333（弧度）CRD_IN_UL输入DEM的左上角的经纬度。CRD_DENSE过采样输入DEM的倍数。对NEAREST为2适合，对TRILINEAR为01适合.越小越耗内存。CRD_OUT_FILE由DEM反算得到的平地相位文件。CRD_OUT_DEM_LP原DEM得到的当前雷达坐标系下的新DEM（已采样）。视用途来定，一般情况，做PS时才用到。subtrrefdem （减去由粗DEM反算得到的平地相位）SRD_OUT_CINT去平地效应后的复干涉图文件。SRD_OFFSET基于粗DEM得到的平地相位与原始干涉合成相位之间的间隔差异，由前面得到，正值表示合成相位向右和向上的一个偏移。可

16、选2.6，第六步：相位解缠。filtphase 干涉图滤波。unwrap 相位解缠。dinsar 差分干涉处理（只对DINSAR有用）filtphase（干涉图滤波）PF_IN_FILE待滤波的干涉图文件。按顺序运行时，该参数C掉，如果只运行该步时必须，且文件名后，要加上影像的宽度。PF_OUT_FILE滤波后的干涉图文件。PF_METHOD滤波方法,（有goldstein（Goldstein滤波），spatialconv（空间回旋滤波）两种）。必需，默认spatialconv。unwrap（相位解缠）UW_METHOD解缠方法，（有SNAPHU（网络流解缠方法）和RAMON（支切解缠）。必需

17、，默认SNAPHU。UW_SEEDS RAMON（支切解缠）时用到的参数，用以给定初始解缠的种子数个数。UW_OUT_FILE解缠输出文件。UW_OUT_FORMAT解缠输出文件格式，有（HGT（振幅，相位，振幅，相位。）和REAL4（float）两种格式）。以下括号内为SNAPHU解缠参数UW_SNAPHU_LOG解缠日志文件。UW_SNAPHU_coh解缠所需干涉图。UW_SNAPHU_MODE解缠模式。有（TOPO（用于地形测量），DEFO（差分干涉测量），SMOOTH（平滑地形区域测量），NOSTATCOSTS（非统计学模式测量）必需，默认TOPO（DEM生成时用）。UW_SNAPHU_INIT解缠所用的网