Insar操作过程

Insar操作过程（一）

仰天长啸 发表于2010年03月29日 18:28 阅读(9) 评论(0) 分类： InSAR 举报

1，命令介绍.

生成模板：run –g 编辑模板：

编辑第一步模板：run –e1 编辑第二步模板：run –e2 编辑第三步模板：run –e3 编辑第四步模板：run –e4 编辑第五步模板：run –e5 编辑第六步模板：run –e6 编辑第七步模板：run –e7 编辑第八步模板：run –e8

（每行前加 C 表示屏蔽该行。 打开模板文件后，按“i“表示可以对进行插入操作，如果删去时，直接用delete或后退键。按ESC键时则退出当前的编辑状态。若要退出文件时，则先按ESC，然后：q！，或：wq，其中：q！表示退出不保存修改，：wq为保存修改。） 运行处理：

运行第一步：run –s1 运行第二步：run –s2 运行第三步：run –s3 运行第四步：run –s4 运行第五步：run –s5 运行第六步：run –s6 运行第七步：run –s7 运行第八步：run –s8 2，参数解说及设置。

2.1，第一步：读取主影像信息。 （1）需要处理步骤。

必需：m_readfiles(读取主影像头文件信息). 必需：m_crop(读取主影像数据体信息).

不需：m_ovs（主影像过采样步骤）：作干涉处理时，一般不运行。

可选：m_prorbits（读取主影像精轨文件）:若要使用提供的精密轨道数据时，就执行该行，否则直接C掉。 （2）小步骤参数说明及设置。

############################################################ INSAR通用参数设置，每步都有的。

###########################################################

名称：LOGFILE

用处： 用以记载DORIS处理过程中所有信息。 文件名：doris.log. 名称：M_RESFILE

用处： 用以记载主影像文件所有信息。 文件名：master.res 名称：S_RESFILE

用处： 用以记载辅影像文件所有信息。 文件名：slave.res 名称：I_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的VOLUME文件 名称：M_IN_LEA SLC的LEADER文件 名称：M_IN_DAT SLC的DAT文件

名称：M_IN_NULL SLC的null文件，默认为dummy ｝

针对ENVISAT格式的参数文件｛

名称：M_IN_DAT SLC的DAT文件(ASAR只有一个文件) ｝

############################################################ 分步骤：M_crops（读取主影像数据体信息）

########################################################### 名称：M_CROP_IN

用处：输入数据体的路径。 说明：必需 .dat文件。 名称：M_CROP_OUT

用处：输出的raw数据体路径。 说明：必需 名称：M_DBOW

用处：要处理的数据区的大小，直接按影像大小来计数。 说明：可选 ，C掉时表示默认为整块影像大小。 名称：M_DBOW_GEO

用处：要处理的数据区的大小，按按经纬度来计数。 说明：一般不用这个，它与M_DBOW不能同选。

############################################################

分步骤：M_ovr（主影像过采样，针对PSINSAR处理的）

########################################################### ############################################################ 分步骤：M_porbits（读取精密轨道信息）

########################################################### 名称：M_ORBDIR

用处：精密轨道文件路径。

说明：针对精轨处理才给定路径。 2.2，第二步：读辅影像信息。

设置基本与第一步类同。这里不作详细介绍。 2.3，第三步：影像配准。 （1）需要处理步骤

必需：coarseorb 一级轨道级配准。 必需：coarsecorr 二级像素级粗配准 必需：m_filtazi 主影像方位向滤波 必需： s_filtazi 辅影像方位向滤波 不需：filtrange 距离向滤波

必需：fine 三级子像素精配准 (2)小步骤参数设置说明

############################################################ 分步骤：coarseorb（一级轨道级配准）

########################################################### ############################################################ 分步骤：coarsecorr（二级像素级粗配准）

########################################################### CC_METHOD magfft |magspace 相关系数计算的方法，

magfft 是先使用fft变换到频率域，然后计算相关系数。——推荐使用 magspace 是直接在空间域计算相关系数。 CC_IN_POS filename

ASCII 文件，记录主影像上待匹配的点的坐标。（不采用） CC_NWIN 21

像素级配准划分窗口的数目，大于５。当CC_IN_POS被选择时，就不需要设置。 CC_WINSIZE 256 256

匹配窗口的大小，处理的时候被改变为奇数。 CC_ACC 64 32

在空间域需要用户设置，在方位向和距离向的搜索偏移量。 在频率域则自动设置为 CC_WINSIZE的一半。 CC_INITOFF 0 0| “orbit” 初始的偏移量。

“orbot” 直接读入使用轨道计算的偏移量。 0 0 用户之际给点的初始偏移量

############################################################

分步骤：m_filtazi（主影像方位向滤波）s_filtazi（辅影像方位向滤波）

########################################################### AF-BLOCKSIZE 1024

FFT变换中每次距离向中的行数。一般是越大越好。 AF_OVERLAP

HANMMING 0.75

光谱滤波的权值。滤波器的alpha 系数 AF_OUT_MASTER master.afiltered 滤波后的主影像文件名

AF_OUT_SLAVE slave.afileered 滤波后的辅影像文件名 AF_OUT_FORMAT cr4

输出文件数据的格式， cr4 ：complex real4 ci2 ：complex shorts ############################################################ 分步骤：filtrange（距离向滤波）

########################################################### ############################################################ 分步骤：fine （三级子像素精配准）

########################################################### FC_METHOD magfft | magspace | oversample 精配准计算方法。

Magfft： 幅度的fft变换后计算相关系数法 Magspace ：幅度直接计算系数法

Overspace ：最好的方法，过采样处理。 FC_NWIN 601

子像素级配准时划分窗口的数目，比较大才好。当CC_IN_POS被选择时，就不需要设置。

FC-IN_POS file name

ASCII 文件，记录主影像上待匹配的点的坐标。（不采用） FC_WINSIZE 64 64 匹配窗口的大小。 FC_ACC 8 8

在空间域需要用户设置，在方位向和距离向的搜索偏移量。 FC_INITOFF 0 0 |COARSECORR

同CC_INITOFF ，COARSECORR表明使用像素级配准后的偏移量。 FC_OSFACTOR 16 过采样的内插倍数。

2.4，第四步：影像重采样。 （1）需要处理步骤

必需：coregpm 估算配准多项式系数。 必需：resample 辅影像重采样。 可选：filtrange 距离向滤波 (2)小步骤参数设置说明

############################################################ 分步骤：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：权矩阵为相关系数的1倍. CPM_MAXITER 20

最小二乘运算的最大迭代次数。

############################################################ 分步骤：resample （辅影像重采样）

########################################################### RS_METHOD RECT | TRI |CC4P| CC6P| TS6P | TS8P | 。。。。。 重采样的内插核.默认RC6P

RS_OUT_FILE s_resampled.raw 重采样后生成的文件名。 RS_OUT_FORMAT CR4 | CI2

输出文件数据类型。CR4： comlex_real 4；CI2：complex_short RS_DBOW 格式（a，a，a，a）

重采样的区域大小。C掉表示整副采样。 （这里的区域坐标是针对主影像而言的）。 RS_SHIFTAZI ON｜OFF

ON，表示采样时，采样核随着多普勒中心频率变化而变化， OFF，则表示不随便变化而采用固定采样核。

一般情况下，如果多普勒中心频率变化很小，可以OFF掉， 默认为ON。

############################################################ 分步骤：filtrange （距离向滤波）

###########################################################

2.5，第五步：干涉图生成。 （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_GEGREE

用处：位置与平地相位之间的多项式系数，一般二次 说明：必需，默认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）两种。

由于是粗DEM，与影像的分辨率不一，需要过采样到同一分辩率下，过采样方式两种，一种是最近邻元方法（NEAREST），一种是内插方式， 说明：必需，默认为TRILINEAR

名称：CRD_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.00833333333（弧度） 名称：CRD_IN_UL

用处：输入DEM的左上角的经纬度。按（纬度，经度）输入 说明：必需。 名称：CRD_DENSE

用处：过采样输入DEM的倍数。对NEAREST为2适合，对TRILINEAR为0—1适合.越小越耗内存。 说明：必需

名称：CRD_OUT_FILE

用处：由DEM反算得到的平地相位文件。 说明：必需

名称：CRD_OUT_DEM_LP

用处：原DEM得到的当前雷达坐标系下的新DEM（已采样）。 说明：视用途来定，一般情况，做PS时才用到。

############################################################ 分步骤：subtrrefdem （减去由粗DEM反算得到的平地相位） ########################################################### 名称：SRD_OUT_CINT

用处：去平地效应后的复干涉图文件。 说明：必需 名称：SRD_OFFSET

用处：基于粗DEM得到的平地相位与原始干涉合成相位之间的间隔差异，由前面得到，正值表示合成相位向右和向上的一个偏移。 说明：可选

2.6，第六步：相位解缠。 （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（平滑地形区域测量），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

（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

本文来源：https://www.bwwdw.com/article/7ju6.html