InSAR处理doris使用详细步骤
- 格式：doc
- 大小：214.50 KB
- 文档页数：28
doris使用--DORIS干涉处理过程简介。
2009-06-11 21:441，Doris处理可以分步处理，也可以整体处理。
Doris处理要先生成模板文件，doris每步处理都对应一个模板文件，这些模板文件包含了doris每步处理所需要的参数，因此用户可以直接打开模板文件，用以编辑更改然后保存，也可以在linux终端中采用Vi命令，或nedit（相当于window下的写字板）来进行编辑保存。
2，查看版本信息，终端键入：doris 即可查看。
3，生成模板。
命令：run –g说明：会在当前工作目录下生成模板文件，它包含三个文件夹目录：Inputfiles，Outdata ，Outinfo。
Outdata用于存放每步生成的结果数据文件。
Outinfo用于存放每步处理的log文件，记录中间处理过程。
Inputfiles中存放的就是模板文件。
它包含的文件：对这些模板进行编辑时，文件中每行前面加c和空格，表示程序处理时不读取该行参数，否则都会执行到的，编辑参数选项时直接更改后面的的数值即可。
4，doris处理步骤总体介绍。
Doris的INSAR处理共分八步：每步又可以分为许多小步，可以通过编辑来选择执行或不执行某些小步骤（小步骤在模板文件中的标识符为process）第一步run –s1：读入主影像的头文件及数据体信息以及轨道信息（可以使用粗轨），其对应的模板文件input.m_initial.，这步包含：读取主影像头文件信息，读取主影像数据体信息，主影像过采样，读取主影像精轨文件。
（如下：加c 表示屏蔽该行，以下都同）第二步run –s2：读入辅影像的头文件及数据体信息以及轨道信息（可以使用粗轨），其对应的模板文件：input.s_initial。
这步包含：读取辅影像头文件信息，读取辅影像数据体信息，辅影像过采样，读取辅影像精轨文件。
第三步run –s3：影像配准。
对应的模板文件：input.coregistration.包含步骤有：一级轨道级配准，二级像素级粗配准，三级子像素级精配准，主影像方位向滤波，辅影像方位向滤波，以及影像距离向滤波。
第四步run –s4：影像重采样。
对应的模板文件：input.resample。
包含步骤有：估计配准多项式系数，影像重采样，以及距离向滤波。
第五步run –s5：生成干涉图。
对应的模板文件：input.products.。
包含步骤有：干涉图生成，除平地效应，以及生成相干系数图。
其中去平地效应有两种方式：一种是基于轨道去平地效应（如果在第一二步提供精轨，那么就是基于精轨去平地效应，否则是基于粗轨去平地效应）。
其中comprefpha步骤为用轨道来计算平地相位。
Subtrrefpha步骤为减去平地相位。
第二种是基于粗DEM去平地效应。
其中comprefdem 为利用提供的粗DEM来反算求得平地相位，subtrrefdem为减去由粗DEM反算得到的平地相位。
第六步run –s6：干涉图滤波及解缠。
对应的模板文件：input.filter_unwrap.包含的步骤有：干涉图滤波，相位解缠，以及差分干涉处理（针对Dinsar）注意的是doris没有提供强有力的解缠程序，它只提供一种简单的枝切法解缠，这对于复杂地形下的干涉图解缠会存在很大误差，因此它一般调用比较有名的snaphu解缠算法来完成解缠过程。
Snaphu解缠也可以单独拿出来进行处理，第七步run –s7：DEM生成。
对应的模板文件：input.s2h_geocode. 。
它包含的步骤有：相位高程转换，地理编码第八步run –s9（整体的步骤包括第三至第七）：doris分步处理前七步即完成，这里的第八步实际上是整体处理，它直接包含第三至第七步的全过程，对参数设置进行了简化，它对应的模板设置文件：input.quicklook.。
该模板文件是将前三到七步模板文件中认为对insar处理比较关键的参数进行了一个集中设置。
因此若要进行整体处理时，必须也先要运行第一步和第二步，然后再运行第八步，即可。
doris使用--InSAR处理过程及参数说明（1）2009-06-11 21:511，命令介绍.生成模板：run –g编辑模板：编辑第一步模板：run –e1编辑第二步模板：run –e2编辑第三步模板：run –e3编辑第四步模板：run –e4编辑第五步模板：run –e5编辑第六步模板：run –e6编辑第七步模板：run –e7编辑第八步模板：run –e9（每行前加 C 表示屏蔽该行。
打开模板文件后，按“i“表示可以对进行插入操作，如果删去时，直接用delete或后退键。
按ESC键时则退出当前的编辑状态。
若要退出文件时，则先按ESC，然后：q！，或：wq，其中：q！表示退出不保存修改，：wq为保存修改。
）运行处理：运行第一步：run –s1运行第二步：run –s2运行第三步：run –s3运行第四步：run –s4运行第五步：run –s5运行第六步：run –s6运行第七步：run –s7运行第八步：run –s92，参数解说及设置。
2.1，第一步：读取主影像信息。
1)在文件input.m_initial中处理（1）需要处理步骤。
必需：m_readfiles(读取主影像头文件信息).必需：m_crop(读取主影像数据体信息).不需：m_ovs（主影像过采样步骤）：作干涉处理时，一般不运行。
可选：m_prorbits（读取主影像精轨文件）:若要使用提供的精密轨道数据时，就执行该行，否则直接C掉。
（2）小步骤参数说明及设置。
############################################################INSAR通用参数设置，每步都有的。
###########################################################名称：LOGFILE用处：用以记载DORIS处理过程中所有信息。
文件名：doris.log.名称：M_RESFILE用处：用以记载主影像文件所有信息。
文件名：master.res名称：S_RESFILE用处：用以记载辅影像文件所有信息。
doris使用--InSAR处理过程及参数说明（2）2009-06-11 21:522.5，第五步：干涉图生成。
5)在文件input.products中处理（1）需要处理步骤必需：interfero 干涉图生成。
必需：coherence生成相干系数图括号里也是必须的，但必须是四选二，或选一和二，或选三和四，不能全选。
{必需：comprefpha 用轨道来计算平地相位必需：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_MULTILOOK用处：多视处理系数。
说明：必需。
如果没有多视要求，则默认(1 1)，如果在这一步参考相位没有去除，则多视参数不应太大，一般默认为：(10 2)或者(5 1)、(8 2).############################################################分步骤：comprefpha （用轨道来计算平地相位）###########################################################名称：FE_METHOD用处：采用计算平地相位的方法，只有精轨(porbits)，如果C掉此行，则为粗轨。
说明：可选。
名称：FE_DEGREE用处：位置与平地相位之间的多项式系数，一般二次说明：必需，默认5。
名称：FE_NPOINTS用处：最小二乘法来计算平地相位所用来拟合的点数说明：必需，默认501。
############################################################分步骤：Subtrrefpha （减去平地相位）###########################################################名称：SRP_METHOD用处：采用减去平地相位的方法，有（polynomial（多项式）和exact）两种多项式方式表示减去采用前一步采用拟合的多项式得到的平地相位，exact表示采用三个核心方程逐点计算平地相位并减去。
说明：默认为polynomial。
名称：SRP_OUT_CINT用处：去平地效应后的干涉图（复数float型）。
说明：cint.minrefpha.raw。
名称：SRP_MULTILOOK用处：多视系数说明：必需，默认（2 2），注意，这里要与前一步的多视系数相乘等于真实多视系数。
############################################################分步骤：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）两种。
