InSAR处理doris使用详细步骤Word文档下载推荐.docx

1、对应的模板文件：input.coregistration.包含步骤有：一级轨道级配准，二级像素级粗配准，三级子像素级精配准 ，主影像方位向滤波，辅影像方位向滤波，以及影像距离向滤波。 第四步run Cs4：影像重采样。input.resample。包含步骤有：估计配准多项式系数，影像重采样，以及距离向滤波。 第五步run Cs5：生成干涉图。input.products.。干涉图生成，除平地效应，以及生成相干系数图。 其中去平地效应有两种方式： 一种是基于轨道去平地效应（如果在第一二步提供精轨，那么就是基于精轨去平地效应，否则是基于粗轨去平地效应）。其中comprefpha步骤为用轨道来计算平

2、地相位。Subtrrefpha步骤为减去平地相位。 第二种是基于粗DEM去平地效应。其中comprefdem 为利用提供的粗DEM来反算求得平地相位，subtrrefdem为减去由粗DEM反算得到的平地相位。 第六步run Cs6：干涉图滤波及解缠。input.filter_unwrap.包含的步骤有：干涉图滤波，相位解缠，以及差分干涉处理（针对Dinsar） 注意的是doris没有提供强有力的解缠程序，它只提供一种简单的枝切法解缠，这对于复杂地形下的干涉图解缠会存在很大误差，因此它一般调用比较有名的snaphu解缠算法来完成解缠过程。Snaphu解缠也可以单独拿出来进行处理， 第七步run

3、Cs7：DEM生成。input.s2h_geocode. 。它包含的步骤有：相位高程转换，地理编码 第八步run Cs9（整体的步骤包括第三至第七）：doris分步处理前七步即完成，这里的第八步实际上是整体处理，它直接包含第三至第七步的全过程，对参数设置进行了简化，它对应的模板设置文件：input.quicklook.。该模板文件是将前三到七步模板文件中认为对insar处理比较关键的参数进行了一个集中设置。 因此若要进行整体处理时，必须也先要运行第一步和第二步，然后再运行第八步，即可。 doris使用-InSAR处理过程及参数说明（1） 2009-06-11 21:51 1，命令介绍. 生成模

4、板：run Cg 编辑模板： 编辑第一步模板：run Ce1 编辑第二步模板：run Ce2 编辑第三步模板：run Ce3 编辑第四步模板：run Ce4 编辑第五步模板：run Ce5 编辑第六步模板：run Ce6 编辑第七步模板：run Ce7 编辑第八步模板：run Ce9 （每行前加 C 表示屏蔽该行。 打开模板文件后，按“i“表示可以对进行插入操作，如果删去时，直接用delete或后退键。按ESC键时则退出当前的编辑状态。若要退出文件时，则先按ESC，然后：q！，或：wq，其中：表示退出不保存修改，：wq为保存修改。） 运行处理： 运行第一步：run Cs1 运行第二步：run

5、Cs2 运行第三步：run Cs3 运行第四步：run Cs4 运行第五步：run Cs5 运行第六步：run Cs6 运行第七步：run Cs7 运行第八步：run Cs9 2，参数解说及设置。 2.1，第一步：读取主影像信息。 1）在文件input.m_initial中处理 （1）需要处理步骤。 必需：m_readfiles（读取主影像头文件信息）. 必需：m_crop（读取主影像数据体信息）. 不需：m_ovs（主影像过采样步骤）：作干涉处理时，一般不运行。 可选：m_prorbits（读取主影像精轨文件）:若要使用提供的精密轨道数据时，就执行该行，否则直接C掉。 （2）小步骤参数说明及

6、设置。 # INSAR通用参数设置，每步都有的。 # 名称：* 用处： 用以记载DORIS处理过程中所有信息。 文件名：doris.log. 名称：M_* 用以记载主影像文件所有信息。master.res 名称：S_* 用以记载辅影像文件所有信息。 doris使用-InSAR处理过程及参数说明（2） 2009-06-11 21:52 2.5，第五步：干涉图生成。 5）在文件input.products中处理 （1）需要处理步骤 必需：interfero 干涉图生成。coherence生成相干系数图 括号里也是必须的，但必须是四选二，或选一和二，或选三和四，不能全选。 comprefpha 用轨

7、道来计算平地相位 必需：Subtrrefpha 减去平地相位comprefdem 利用提供的粗DEM来反算求得平地相位 必需：subtrrefdem 减去由粗DEM反算得到的平地相位 （2）小步骤参数设置说明 # 分步骤：interfero （干涉图生成）INT_OUT_CINT（或INT_OUT_INT）生成的干涉图文件（INT_OUT_CINT为复数float型， INT_OUT_INT为float型，）。必需，但二选一。 名称：INT_METHOD生成的干涉图的方法，有（old，oversample）两种， 说明：必需，一般选oversample。INT_*OK 用处：多视处理系数。必需

8、。如果没有多视要求，则默认（1 1），如果在这一步参考相位没有去除，则多视参数不应太大，一般默认为：（10 2）或者（5 1）、（8 2）. # 分步骤：comprefpha （用轨道来计算平地相位）FE_METHOD采用计算平地相位的方法，只有精轨（porbits），如果C掉此行，则为粗轨。可选。FE_DEGREE位置与平地相位之间的多项式系数，一般二次 说明：必需，默认5。FE_*最小二乘法来计算平地相位所用来拟合的点数 说明：必需，默认501。Subtrrefpha （减去平地相位）SRP_METHOD采用减去平地相位的方法，有（polynomial（多项式）和exact）两种 多项式方

9、式表示减去采用前一步采用拟合的多项式得到的平地相位，exact表示采用三个核心方程逐点计算平地相位并减去。默认为polynomial。SRP_OUT_CINT去平地效应后的干涉图（复数float型）。cint.minrefpha.raw。SRP_*OK 用处：多视系数必需，默认（2 2），注意，这里要与前一步的多视系数相乘等于真实多视系数。coherence（成相干系数图）COH_OUT_CCOH（或COH_OUT_COH）生成的相干图文件（COH_OUT_CCOH为复数float型，COH_OUT_COH为float型，），二者只能选一。必需，但二选一，给定其路径。COH_*OK 用处：多视

10、系数。必需，默认（10 2） 名称：COH_*用于计算相干系数的窗口。必需，（15 15） 名称：COH_METHOD用于计算相干系数的方法 说明：可选 LD；comprefdem （利用提供的粗DEM来反算求得平地相位） # 名称：CRD_METHOD求取平地相位方法，有（*和*AR）两种。 由于是粗DEM，与影像的分辨率不一，需要过采样到同一分辩率下，过采样方式两种，一种是最近邻元方法（*），一种是内插方式， 说明：必需，默认为*AR 名称：CRD_IN_FORMAT输入粗DEM的格式，有（I2（short int），R4（float），R8（double），I2_*AN）,Gtopo30

11、 和SRTM都是I2型。可选 ，C掉时系统默认为I2. 名称：CRD_IN_DEM输入DEM的路径名称。是已经拼好的DEM。必需 名称：CRD_IN_SIZE输入DEM的大小，分经纬向。按（纬度，经度）输入 说明：必需CRD_IN_DELTA输入DEM的间隔，一般经纬向间隔是一样的。必需，SRTM的间隔为 0.*-*333（弧度） 名称：CRD_IN_UL输入DEM的左上角的经纬度。CRD_DENSE过采样输入DEM的倍数。对*为2适合，对*AR为01适合.越小越耗内存。CRD_OUT_FILE由DEM反算得到的平地相位文件。CRD_OUT_DEM_LP原DEM得到的当前雷达坐标系下的新DEM

12、（已采样）。视用途来定，一般情况，做PS时才用到。subtrrefdem （减去由粗DEM反算得到的平地相位） # 名称：SRD_OUT_CINT去平地效应后的复干涉图文件。SRD_OFFSET基于粗DEM得到的平地相位与原始干涉合成相位之间的间隔差异，由前面得到，正值表示合成相位向右和向上的一个偏移。可选 2.6，第六步：相位解缠。 6）在文件input.filter_unwrap中处理 （1）需要处理步骤filtphase 干涉图滤波。unwrap 相位解缠。dinsar 差分干涉处理（只对DINSAR有用） （2）小步骤参数设置说明filtphase（干涉图滤波）PF_IN_FILE 待

13、滤波的干涉图文件。按顺序运行时，该参数C掉，如果只运行该步时必须，且文件名后，要加上影像的宽度。PF_OUT_FILE滤波后的干涉图文件。PF_METHOD滤波方法,（有goldstein（Goldstein滤波），spatialconv（空间回旋滤波）两种）。必需，默认spatialconv。unwrap（相位解缠）UW_METHOD解缠方法，（有SNAPHU（网络流解缠方法）和RAMON（支切解缠）。必需，默认SNAPHU。UW_SEEDS RAMON（支切解缠）时用到的参数，用以给定初始解缠的种子数个数。UW_OUT_FILE 用处：解缠输出文件。UW_OUT_FORMAT 解缠输出文件

14、格式，有（HGT（振幅，相位，振幅，相位。）和REAL4（float）两种格式）。 以下括号内为SNAPHU解缠参数 名称：UW_SNAPHU_LOG 用处：解缠日志文件。UW_SNAPHU_coh 用处：解缠所需干涉图。UW_SNAPHU_MODE解缠模式。有（TOPO（用于地形测量），DEFO（差分干涉测量），SMOOTH（平滑地形区域测量），*OSTS（非统计学模式测量） 说明：必需，默认TOPO（DEM生成时用）。UW_SNAPHU_INIT解缠所用的网络流生成方式。有（MST（最小生成树），MCF（最小费用流）。必需，默认MST。 # 分步骤：*NTIAL INSAR （差分干涉处理

15、） 以下括号内是在做差分处理的形变对操作时才用到 名称：DI_OUT_FILE 用处：生成的差分干涉图。DI_IN_*VE地形对的辅影像信息文件（地形对的slave.res文件）。DI_IN_* 地形对的生成干涉图步骤的信息文件（地形对的products.res）。DI_IN_*TER地形对的主影像信息文件（地形对的master.res文件）。只针对四轨法有效，对三轨无效。 2.7，第七步： 7）在文件input.s2h_geocode中处理 （1）需要处理步骤slant2h 相位高程转换。geocode地理编码。slant（相位高程转换） #S2H_METHOD高程转换方法，有（ambigu

16、ity，schwabish（最小二乘拟合），rodriguez）三种选择。S2H_OUT_HEI输出的高程值文件（hei.raw）。S2H_OUT_PHI输出的经度值文件（phi.raw）。只针对ambiguity方法。S2H_OUT_LAM输出的纬度值文件（lam.raw）。S2H_*点数。计算N个点在不同高程处的相位与高程关系，进而拟合二者之间的关系多项式，。只针对schwabisch方法，默认200。S2H_*D拟合多项式的次数（相位对高程的系数）（1次时为2，二次为5）。只针对schwabisch方法.默认为2。S2H_*S 用处：给定的控制点数目。只针对schwabisch方法，一般

17、为S2H_*D的值1。拟合多项式的次数（高程对位置的系数）（1次时为2，二次为5）。只针对schwabisch方法.默认为5。geocoding（地理编码）GEO_OUT_PHI输出的地理编码后的经度值文件（phi.raw）。GEO_OUT_LAM输出的地理编码后的纬度值文件（lambda.raw）。 2.8，整体处理。 （1）需要处理步骤 Run Cs1 Run Cs2 Run Cs9 第八步处理参数在第三到第七步中有介绍。 doris使用-DInSAR三轨差分与四轨差分法实施过程 2009-06-11 21:53 DInSAR处理过程。 在DORIS中，DINSAR有三轨和四轨两种方式。

18、不论在三轨还是四轨都要保证是在同一坐标下，并且都配准同一坐标系，即所有图像可以重合的。 认真看了下官网的设置.一行一行找.终于发现.其中resample步有窗口大小的设置,通过这一步来保证两者的重叠. 因为在doris处理中，resample步设用的窗口坐标是以主影像为基准的，而且三轨法差分两次用的都是同一主影像，因此经过配准后，在topo对和defo对的resample步设置RS_DBOW为同样值，是可以保证两者重叠的。（RB_DBOW设置时必须保证区域都在重叠区域内，不超限） 至此,我以为三轨差分的步骤: （1）.先做defo对,至解缠一步.记下resample时的slave采样时的窗口文

19、件大小,可从master_slave.res中得到.（这个大小可以是默认值（即屏掉该行），即全副配准，也可以是在resample这一步给定一个区域，要注意的是，给定的区域不要超限，建议手工设置比较好，这样更有可能保证不超限） （2）.再新建工作目录，做topo对,至resample时,打开其模板文件,起用RS_DBOW参数,将其窗口 大小设为第一步记下的defo的master_slave.res中的RS_DBOW值,然后执行,一直到解缠步（这里解缠只执行滤波步）. （3）.在defo工作目录中,run -e6,启用DInsar步,屏掉滤波和解缠两步. 将topo对的解缠结果拷至defo的out

20、data下,然后给出topo的slave.res和master-slave.res的路径,即DI_IN_*VE 和DI_IN_TOPINT的参数值. （4）运行run -s6 四轨差分步骤： 如1，2，3，4四影像，1，2为topo对，3，4为defo对， （1） 选定一个主影像，（以1为公共主影像吧）。新建三个工作目录，作三次 配准，（1，2），（1，3），（1，4）;至重采样步，注意在resmaple步设置所有的RS_DBOW为相同值，且保证不超限。 （2） （1,2）继续做到解缠为止。 （3） 将（1，3）中相关3的文件，如，slave.res,以及生成在其outdata中的 重采样文件

21、，放到（1.4）中对应的位置，记住要用重采样文件替换掉（1，4）中1的raw文件，用slave.res替换掉（1，4）中的slave.res,相应的slave.res,并更改其中的Data_output_file:将其改为新导入的那个重采样文件的路径。 （4） 这时，（1，2）就为topo目录，新的（1，4）就为defo目录 （5） 对defo对继续执行，一直到解缠步（这里解缠只执行滤波步） （6）在defo工作目录中,run -e6,启用DInsar步,屏掉滤波和解缠两步. 将topo对的解缠结果拷至defo的outdata下,然后给出topo的slave.res和master-slave.

22、res的路径,即DI_IN_*VE 和DI_IN_TOPINT的参数值. （7）运行run -s6 用doris进行INSAR的操作过程 *操作时屏幕上显示的处理过程* 默认分类 2009-06-09 09:58:03 阅读289 评论2 字号：大中小 Xiaoyanhp-47c03a5dcc20 /cygdrive/e/insar/doristest_Data/Your_work $ cd e: Xiaoyanhp-47c03a5dcc20 /cygdrive/e $ cd insar/ Xiaoyanhp-47c03a5dcc20 /cygdrive/e/insar $ cd test Xiaoyanhp-47