说明:按顺序运行时,该参数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(平滑地形区域测量),NOSTATCOSTS(非统计学模式测量)) 说明:必需,默认TOPO(DEM生成时用)。 名称:UW_SNAPHU_INIT
用处:解缠所用的网络流生成方式。有(MST(最小生成树),MCF(最小费用流))。 说明:必需,默认MST。 }
############################################################ 分步骤:dinsar (差分干涉处理)
以下括号内是在做差分处理的形变对操作时才用到{ 名称:DI_OUT_FILE 用处:生成的差分干涉图。 说明:必需。
名称:DI_IN_TOPOSLAVE
用处:地形对的辅影像信息文件(地形对的slave.res文件)。 说明:必需。 名称:DI_IN_TOPOINT
用处: 地形对的生成干涉图步骤的信息文件(地形对的products.res)。 说明:必需。
名称:DI_IN_TOPOMASTER
用处:地形对的主影像信息文件(地形对的master.res文件)。 说明:只针对四轨法有效,对三轨无效。 }
2.7,第七步:DEM生成。 (1)需要处理步骤
必需:slant2h 相位高程转换。 必需:geocode地理编码。 (2)小步骤参数设置说明
############################################################ 分步骤:slant2h(相位高程转换)
########################################################### 名称:S2H_METHOD
用处:高程转换方法,有(ambiguity,schwabish(最小二乘拟合),rodriguez)三种选择。 说明:必需。 名称:S2H_OUT_HEI
用处:输出的高程值文件(hei.raw)。 说明:必需。 名称:S2H_OUT_PHI
用处:输出的经度值文件(phi.raw)。 说明:只针对ambiguity方法。
名称:S2H_OUT_lAM
用处:输出的纬度值文件(lam.raw)。 说明:只针对ambiguity方法。 名称:S2H_NPOINT
用处:点数。计算N个点在不同高程处的相位与高程关系,进而拟合二者之间的关系多项式,。
说明:只针对schwabisch方法,默认200。 名称:S2H_DEGREE1D
用处:拟合多项式的次数(相位对高程的系数)(1次时为2,二次为5)。 说明:只针对schwabisch方法.默认为2。 名称:S2H_NHEIGHTS 用处:给定的控制点数目。
说明:只针对schwabisch方法,一般为S2H_DEGREE1D的值+1。 名称:S2H_DEGREE2D
用处:拟合多项式的次数(高程对位置的系数)(1次时为2,二次为5)。 说明:只针对schwabisch方法.默认为5。
############################################################ 分步骤:geocode(地理编码)
########################################################### 名称:GEO_OUT_PHI
用处:输出的地理编码后的经度值文件(lam.raw)。 说明:必需。 名称:GEO_OUT_LAM
用处:输出的地理编码后的纬度值文件(phi.raw)。
说明:必需。 2.8,整体处理。 (1)需要处理步骤 Run –s1 Run –s2 Run –s9
第八步处理参数在第三到第七步中有介绍。
在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_TOPOSLAVE 和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_TOPOSLAVE 和DI_IN_TOPINT的参数值. (7)运行run -s6