InSAR处理doris使用详细步骤.docx
《InSAR处理doris使用详细步骤.docx》由会员分享，可在线阅读，更多相关《InSAR处理doris使用详细步骤.docx（16页珍藏版）》请在冰豆网上搜索。
InSAR处理doris使用详细步骤
InSAR处理doris使用详细步骤
doris使用--DORIS干涉处理过程简介。
2009-06-1121:
44
1，Doris处理可以分步处理，也可以整体处理。
Doris处理要先生成模板文件，doris每步处理
都对应一个模板文件，这些模板文件包含了doris每步处理所需要的参数，因此用户可以直接打开模板文件，用以编辑更改然后保存，也可以在linux终端中采用Vi命令，或nedit（相当于window下的写字板）来进行编辑保存。
2，查看版本信息，终端键入：
doris即可查看。
3，生成模板。
命令：
runCg
说明：
会在当前工作目录下生成模板文件，它包含三个文件夹目录：
Inputfiles，Outdata，Outinfo。
Outdata用于存放每步生成的结果数据文件。
Outinfo用于存放每步处理的log文件，记录中间处理过程。
Inputfiles中存放的就是模板文件。
它包含的文件：
对这些模板进行编辑时，文件中每行前面加c和空格，表示程序处理时不读取该行参数，否则都会执行到的，编辑参数选项时直接更改后面的的数值即可。
4，doris处理步骤总体介绍。
Doris的INSAR处理共分八步：
每步又可以分为许多小步，可以通过编辑来选择执行或不执行某些小步骤（小步骤在模板文件中的标识符为process）第一步runCs1：
读入主影像的头文件及数据体信息以及轨道信息（可以使用粗轨），其对应的模板文件input.m_initial.，这步包含：
读取主影像头文件信息，读取主影像数据体信息，主影像过采样，读取主影像精轨文件。
（如下：
加c表示屏蔽该行，以下都同）
第二步runCs2：
读入辅影像的头文件及数据体信息以及轨道信息（可以使用粗轨），其对应的模板文件：
input.s_initial。
这步包含：
读取辅影像头文件信息，读取辅影像数据体信息，辅影像过采样，读取辅影像精轨文件。
第三步runCs3：
影像配准。
对应的模板文件：
input.coregistration.包含步骤有：
一级轨道级配准，二级像素级粗配准，三级子像素级精配准，主影像方位向滤波，辅影像方位向滤波，以及影像距离向滤波。
第四步runCs4：
影像重采样。
对应的模板文件：
input.resample。
包含步骤有：
估计配准多项式系数，影像重采样，以及距离向滤波。
第五步runCs5：
生成干涉图。
对应的模板文件：
input.products.。
包含步骤有：
干涉图生成，除平地效应，以及生成相干系数图。
其中去平地效应有两种方式：
一种是基于轨道去平地效应（如果在第一二步提供精轨，那么就是基于精轨去平地效应，否则是基于粗轨去平地效应）。
其中comprefpha步骤为用轨道来计算平地相位。
Subtrrefpha步骤为减去平地相位。
第二种是基于粗DEM去平地效应。
其中comprefdem为利用提供的粗DEM来反算求得平地相位，subtrrefdem为减去由粗DEM反算得到的平地相位。
第六步runCs6：
干涉图滤波及解缠。
对应的模板文件：
input.filter_unwrap.包含的步骤有：
干涉图滤波，相位解缠，以及差分干涉处理（针对Dinsar）
注意的是doris没有提供强有力的解缠程序，它只提供一种简单的枝切法解缠，这对于复杂地形下的干涉图解缠会存在很大误差，因此它一般调用比较有名的snaphu解缠算法来完成解缠过程。
Snaphu解缠也可以单独拿出来进行处理，
第七步runCs7：
DEM生成。
对应的模板文件：
input.s2h_geocode.。
它包含的步骤有：
相位高程转换，地理编码第八步runCs9（整体的步骤包括第三至第七）：
doris分步处理前七步即完成，这里的第八步实际上是整体处理，它直接包含第三至第七步的全过程，对参数设置进行了简化，它对应的模板设置文件：
input.quicklook.。
该模板文件是将前三到七步模板文件中认为对insar处理比较关键的参数进行了一个集中设置。
因此若要进行整体处理时，必须也先要运行第一步和第二步，然后再运行第八步，即可。
doris使用--InSAR处理过程及参数说明
（1）2009-06-1121:
51
1，命令介绍.
生成模板：
runCg编辑模板：
编辑第一步模板：
runCe1编辑第二步模板：
runCe2编辑第三步模板：
runCe3编辑第四步模板：
runCe4编辑第五步模板：
runCe5编辑第六步模板：
runCe6编辑第七步模板：
runCe7编辑第八步模板：
runCe9
（每行前加C表示屏蔽该行。
打开模板文件后，按“i“表示可以对进行插入操作，如果删去时，直接用delete或后退键。
按ESC键时则退出当前的编辑状态。
若要退出文件时，则先按ESC，然后：
q！
，或：
wq，其中：
q！
表示退出不保存修改，：
wq为保存修改。
）运行处理：
运行第一步：
runCs1运行第二步：
runCs2运行第三步：
runCs3运行第四步：
runCs4运行第五步：
runCs5运行第六步：
runCs6运行第七步：
runCs7运行第八步：
runCs92，参数解说及设置。
2.1，第一步：
读取主影像信息。
1）在文件input.m_initial中处理
（1）需要处理步骤。
必需：
m_readfiles（读取主影像头文件信息）.必需：
m_crop（读取主影像数据体信息）.
不需：
m_ovs（主影像过采样步骤）：
作干涉处理时，一般不运行。
可选：
m_prorbits（读取主影像精轨文件）:
若要使用提供的精密轨道数据时，就执行该行，否则直接C掉。
（2）小步骤参数说明及设置。
############################################################INSAR通用参数设置，每步都有的。
###########################################################名称：
*****
用处：
用以记载DORIS处理过程中所有信息。
文件名：
doris.log.名称：
M_*****
用处：
用以记载主影像文件所有信息。
文件名：
master.res名称：
S_*****
用处：
用以记载辅影像文件所有信息。
doris使用--InSAR处理过程及参数说明
（2）2009-06-1121:
52
2.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_*****OK用处：
多视处理系数。
说明：
必需。
如果没有多视要求，则默认（11），如果在这一步参考相位没有去除，则多视参数不应太大，一般默认为：
（102）或者（51）、（82）.############################################################分步骤：
comprefpha（用轨道来计算平地相位）
###########################################################名称：
FE_METHOD
用处：
采用计算平地相位的方法，只有精轨（porbits），如果C掉此行，则为粗轨。
说明：
可选。
名称：
FE_DEGREE
用处：
位置与平地相位之间的多项式系数，一般二次说明：
必需，默认5。
名称：
FE_*****
用处：
最小二乘法来计算平地相位所用来拟合的点数说明：
必需，默认501。
############################################################分步骤：
Subtrrefpha（减去平地相位）
###########################################################名称：
SRP_METHOD
用处：
采用减去平地相位的方法，有（polynomial（多项式）和exact）两种多项式方式表示减去采用前一步采用拟合的多项式得到的平地相位，exact表示采用三个核心方程逐点计算平地相位并减去。
说明：
默认为polynomial。
名称：
SRP_OUT_CINT
用处：
去平地效应后的干涉图（复数float型）。
说明：
cint.minrefpha.raw。
名称：
SRP_*****OK用处：
多视系数
说明：
必需，默认（22），注意，这里要与前一步的多视系数相乘等于真实多视系数。
############################################################分步骤：
coherence（成相干系数图）
###########################################################名称：
COH_OUT_CCOH（或COH_OUT_COH）
用处：
生成的相干图文件（COH_OUT_CCOH为复数float型，COH_OUT_COH为float型，），二者只能选一。
说明：
必需，但二选一，给定其路径。
名称：
COH_*****OK用处：
多视系数。
说明：
必需，默认（102）名称：
COH_*****
用处：
用于计算相干系数的窗口。
说明：
必需，（1515）名称：
COH_METHOD
用处：
用于计算相干系数的方法说明：
可选LD；
############################################################分步骤：
comprefdem（利用提供的粗DEM来反算求得平地相位）###########################################################名称：
CRD_METHOD
用处：
求取平地相位方法，有（*****和*****AR）两种。
由于是粗DEM，与影像的分辨率不一，需要过采样到同一分辩率下，过采样方式两种，一种是最近邻元方法（*****），一种是内插方式，说明：
必需，默认为*****AR名称：
CRD_IN_FORMAT
用处：
输入粗DEM的格式，有（I2（shortint），R4（float），R8（double），I2_*****AN）,Gtopo30和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为0―1适合.越小越耗内存。
说明：
必需
名称：
CRD_OUT_FILE
用处：
由DEM反算得到的平地相位文件。
说明：
必需
名称：
CRD_OUT_DEM_LP
用处：
原DEM得到的当前雷达坐标系下的新DEM（已采样）。
说明：
视用途来定，一般情况，做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
用处：
待滤波的干涉图文件。
说明：
按顺序运行时，该参数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
用处：
解缠输出文件格式，有（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。
｝
############################################################
分步骤：
*****NTIALINSAR（差分干涉处理）以下括号内是在做差分处理的形变对操作时才用到｛名称：
DI_OUT_FILE用处：
生成的差分干涉图。
说明：
必需。
名称：
DI_IN_*****VE
用处：
地形对的辅影像信息文件（地形对的slave.res文件）。
说明：
必需。
名称：
DI_IN_*****
用处：
地形对的生成干涉图步骤的信息文件（地形对的products.res）。
说明：
必需。
名称：
DI_IN_*****TER
用处：
地形对的主影像信息文件（地形对的master.res文件）。
说明：
只针对四轨法有效，对三轨无效。
｝
2.7，第七步：
DEM生成。
7）在文件input.s2h_geocode中处理
（1）需要处理步骤
必需：
slant2h相位高程转换。
必需：
geocode地理编码。
（2）小步骤参数设置说明
############################################################分步骤：
slant（相位高程转换）
###########################################################
名称：
S2H_METHOD
用处：
高程转换方法，有（ambiguity，schwabish（最小二乘拟合），rodriguez）三种选择。
说明：
必需。
名称：
S2H_OUT_HEI
用处：
输出的高程值文件（hei.raw）。
说明：
必需。
名称：
S2H_OUT_PHI
用处：
输出的经度值文件（phi.raw）。
说明：
只针对ambiguity方法。
名称：
S2H_OUT_LAM
用处：
输出的纬度值文件（lam.raw）。
说明：
只针对ambiguity方法。
名称：
S2H_*****
用处：
点数。
计算N个点在不同高程处的相位与高程关系，进而拟合二者之间的关系多项式，。
说明：
只针对schwabisch方法，默认200。
名称：
S2H_*****D
用处：
拟合多项式的次数（相位对高程的系数）（1次时为2，二次为5）。
说明：
只针对schwabisch方法.默认为2。
名称：
S2H_*****S用处：
给定的控制点数目。
说明：
只针对schwabisch方法，一般为S2H_*****D的值＋1。
名称：
S2H_*****D
用处：
拟合多项式的次数（高程对位置的系数）（1次时为2，二次为5）。
说明：
只针对schwabisch方法.默认为5。
############################################################分步骤：
geocoding（地理编码）
###########################################################名称：
GEO_OUT_PHI
用处：
输出的地理编码后的经度值文件（phi.raw）。
说明：
必需。
名称：
GEO_OUT_LAM
用处：
输出的地理编码后的纬度值文件（lambda.raw）。
说明：
必需。
2.8，整体处理。
（1）需要处理步骤RunCs1RunCs2RunCs9
第八步处理参数在第三到第七步中有介绍。
doris使用--DInSAR三轨差分与四轨差分法实施过程2009-06-1121:
53
D－InSAR处理过程。
在DORIS中，D－INSAR有三轨和四轨两种方式。
不论在三轨还是四轨都要保证是在同一坐标下，并且都配准同一坐标系，即所有图像可以重合的。
认真看了下官网的设置...一行一行找..终于发现...其中resample步有窗口大小的设置,通过这一步来保证两者的重叠.
因为在doris处理中，resample步设用的窗口坐标是以主影像为基准的，而且三轨法差分两次用的都是同一主影像，因此经过配准后，在topo对和defo对的resample步设置RS_DBOW为同样值，是可以保证两者重叠的。
（RB_DBOW设置时必须保证区域都在重叠区域内，不超限）
至此,我以为三轨差分的步骤:
（1）.先做defo对,至解缠一步..记下resample时的slave采样时的窗口文件大小,可从master_slave.res中得到..（这个大小可以是默认值（即屏掉该行），即全副配准，也可以是在resample这一步给定一个区域，要注意的是，给定的区域不要超限，建议手工设置比较好，这样更有可能保证不超限）
（2）.再新建工作目录，做topo对,至resample时,打开其模板文件,起用RS_DBOW参数,将其窗口
大小设为第一步记下的defo的master_slave.res中的RS_DBOW值,然后执行,一直到解缠步（这里解缠只执行滤波步）.
（3）.在defo工作目录中,run-e6,启用D－Insar步,屏掉滤波和解缠两步.将topo对的解缠结果拷至defo的outdata下,然后给出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中的
重采样文件，放到（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,启用D－Insar步,屏掉滤波和解缠两步.将topo对的解缠结果拷至defo的outdata下,然后给出topo的slave.res和master-slave.res的路径,即DI_IN_*****VE和DI_IN_TOPINT的参数值.（7）运行run-s6
用doris进行INSAR的操作过程
****操作时屏幕上显示的处理过程****
默认分类2009-06-0909:
58:
03阅读289评论2字号：
大中小
Xiaoyan@hp-47c03a5dcc20/cygdrive/e/insar/doristest_Data/Your_work
$cde:
Xiaoyan@hp-47c03a5dcc20/cygdrive/e
$cdinsar/
Xiaoyan@hp-47c03a5dcc20/cygdrive/e/insar
$cdtest
Xiaoyan@hp-47
升级会员 