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生成相干系数图
括号里也是必须的，但必须是四选二，或选一和二，或选三和四，不能全选。