DorisInSAR处理基本步骤和基本参数设定Word文档下载推荐.docx
- 文档编号:5670953
- 上传时间:2023-05-05
- 格式:DOCX
- 页数:21
- 大小:21.82KB
DorisInSAR处理基本步骤和基本参数设定Word文档下载推荐.docx
《DorisInSAR处理基本步骤和基本参数设定Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《DorisInSAR处理基本步骤和基本参数设定Word文档下载推荐.docx(21页珍藏版)》请在冰点文库上搜索。
第五步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处理比较关键的参数进行了一个集中设置。
因此若要进行整体处理时,必须也先要运行第一步和第二步,然后再运行第八步,即可。
1,命令介绍.
生成模板:
编辑模板:
编辑第一步模板:
run–e1
编辑第二步模板:
run–e2
编辑第三步模板:
run–e3
编辑第四步模板:
run–e4
编辑第五步模板:
run–e5
编辑第六步模板:
run–e6
编辑第七步模板:
run–e7
编辑第八步模板:
run–e9
(每行前加C表示屏蔽该行。
打开模板文件后,按“i“表示可以对进行插入操作,如果删去时,直接用delete或后退键。
按ESC键时则退出当前的编辑状态。
若要退出文件时,则先按ESC,然后:
q!
,或:
wq,其中:
表示退出不保存修改,:
wq为保存修改。
)
运行处理:
运行第一步:
run–s1
运行第二步:
run–s2
运行第三步:
run–s3
运行第四步:
run–s4
运行第五步:
run–s5
运行第六步:
run–s6
运行第七步:
run–s7
运行第八步:
run–s9
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格式的参数文件{
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_NWIN21
像素级配准划分窗口的数目,大于5。
当CC_IN_POS被选择时,就不需要设置。
CC_WINSIZE
256256
匹配窗口的大小,处理的时候被改变为奇数。
CC_ACC
6432
在空间域需要用户设置,在方位向和距离向的搜索偏移量。
在频率域则自动设置为CC_WINSIZE的一半。
CC_INITOFF00|“orbit”
初始的偏移量。
“orbot”直接读入使用轨道计算的偏移量。
00
用户之际给点的初始偏移量
m_filtazi(主影像方位向滤波)s_filtazi(辅影像方位向滤波)
AF-BLOCKSIZE1024
FFT变换中每次距离向中的行数。
一般是越大越好。
AF_OVERLAP
HANMMING
0.75
光谱滤波的权值。
滤波器的alpha系数
AF_OUT_MASTERmaster.afiltered
滤波后的主影像文件名
AF_OUT_SLAVE
slave.afileered
滤波后的辅影像文件名
AF_OUT_FORMAT
cr4
输出文件数据的格式,cr4:
complexreal4
ci2:
complexshorts
filtrange(距离向滤波)
fine(三级子像素精配准)
FC_METHOD
magfft|magspace|oversample
精配准计算方法。
Magfft:
幅度的fft变换后计算相关系数法
Magspace:
幅度直接计算系数法
Overspace:
最好的方法,过采样处理。
FC_NWIN
601
子像素级配准时划分窗口的数目,比较大才好。
FC-IN_POS
filename
FC_WINSIZE
6464
匹配窗口的大小。
FC_ACC88
FC_INITOFF00|COARSECORR
同CC_INITOFF,COARSECORR表明使用像素级配准后的偏移量。
FC_OSFACTOR16
过采样的内插倍数。
2.4,第四步:
coregpm估算配准多项式系数。
resample辅影像重采样。
coregpm(估算配准多项式系数)
CPM_THRESHOLD
0.4
经配准后生成的大量匹配点,当相关系数大于CPM_THRESHOLD,才参与平差运算。
阈值的大小与FILE中窗口的大小有关。
小窗口情况下阈值趋于1.0,阈值越大越好。
对于64×
64的窗口,阈值选择0.2即可。
CPM_DEGREE
2
多项式的阶数。
2阶即可。
CPM_WEIGHTBAMLER|NONE|LINEAR|QUADRATIC
最小二乘中权矩阵方法选择。
BAMLER:
权矩阵为相关系数的2倍。
LINEAR:
权矩阵为相关系数的1倍.
CPM_MAXITER
20
最小二乘运算的最大迭代次数。
resample(辅影像重采样)
RS_METHODRECT|TRI|CC4P|CC6P|TS6P|TS8P|。
。
重采样的内插核.默认RC6P
RS_OUT_FILE
s_resampled.raw
重采样后生成的文件名。
RS_OUT_FORMATCR4|CI2
输出文件数据类型。
CR4:
comlex_real4;
CI2:
complex_short
RS_DBOW格式(a,a,a,a)
重采样的区域大小。
C掉表示整副采样。
(这里的区域坐标是针对主影像而言的)。
RS_SHIFTAZION|OFF
ON,表示采样时,采样核随着多普勒中心频率变化而变化,
OFF,则表示不随便变化而采用固定采样核。
一般情况下,如果多普勒中心频率变化很小,可以OFF掉,
默认为ON。
filtrange(距离向滤波)
2.5,第五步:
干涉图生成。
interfero干涉图生成。
coherence生成相干系数图
括号里也是必须的,但必须是四选二,或选一和二,或选三和四,不能全选。
{
comprefpha用轨道来计算平地相位
Subtrrefpha减去平地相位
comprefdem利用提供的粗DEM来反算求得平地相位
subtrrefdem减去由粗DEM反算得到的平地相位
}
interfero(干涉图生成)
INT_OUT_CINT(或INT_OUT_INT)
生成的干涉图文件(INT_OUT_CINT为复数float型,INT_OUT_INT为float型,)。
必需,但二选一。
INT_METHOD
生成的干涉图的方法,有(old,oversample)两种,
必需,一般选oversample。
INT_MULTILOOK
多视处理系数。
必需。
如果没有多视要求,则默认(11),如果在这一步参考相位没有去除,则多视参数不应太大,一般默认为:
(102)或者(51)、(82).
comprefpha(用轨道来计算平地相位)
FE_METHOD
采用计算平地相位的方法,只有精轨(porbits),如果C掉此行,则为粗轨。
可选。
FE_GEGREE
位置与平地相位之间的多项式系数,一般二次
必需,默认5。
FE_NPOINTS
最小二乘法来计算平地相位所用来拟合的点数
必需,默认501。
Subtrrefpha(减去平地相位)
SRP_METHOD
采用减去平地相位的方法,有(polynomial(多项式)和exact)两种
多项式方式表示减去采用前一步采用拟合的多项式得到的平地相位,exact表示采用三个核心方程逐点计算平地相位并减去。
默认为polynomial。
SRP_OUT_CINT
去平地效应后的干涉图(复数float型)。
SRP_MULTILOOK
多视系数
必需,默认(22),注意,这里要与前一步的多视系数相乘等于真实多视系数。
coherence(成相干系数图)
COH_OUT_CCOH(或COH_OUT_COH)
生成的相干图文件(COH_OUT_CCOH为复数float型,COH_OUT_COH为float型,),二者只能选一。
必需,但二选一,给定其路径。
COH_MULTILOOK
多视系数。
必需,默认(102)
COH_WINSIZE
用于计算相干系数的窗口。
必需,(1515)
COH_METHOD
用于计算相干系数的方法
可选LD;
comprefdem(利用提供的粗DEM来反算求得平地相位)
CRD_METHOD
求取平地相位方法,有(NEAREST和TRILINEAR)两种。
由于是粗DEM,与影像的分辨率不一,需要过采样到同一分辩率下,过采样方式两种,一种是最近邻元方法(NEAREST),一种是内插方式,
必需,默认为TRILINEAR
CRD_IN_FORMAT
输入粗DEM的格式,有(I2(shortint),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.003(弧度)
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得到的平地相位与原始干涉合成相位之间的间隔差异,由前面得
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- DorisInSAR 处理 基本 步骤 基本参数 设定