DORIS操作流程详解.docx
- 文档编号:3931086
- 上传时间:2022-11-26
- 格式:DOCX
- 页数:16
- 大小:22.22KB
DORIS操作流程详解.docx
《DORIS操作流程详解.docx》由会员分享,可在线阅读,更多相关《DORIS操作流程详解.docx(16页珍藏版)》请在冰豆网上搜索。
DORIS操作流程详解
1,命令介绍.
生成模板:
run–g
编辑模板:
编辑第一步模板:
run–e1
编辑第二步模板:
run–e2
编辑第三步模板:
run–e3
编辑第四步模板:
run–e4
编辑第五步模板:
run–e5
编辑第六步模板:
run–e6
编辑第七步模板:
run–e7
编辑第八步模板:
run–e9
(每行前加C表示屏蔽该行。
打开模板文件后,按“i“表示可以对进行插入操作,如果删去时,直接用delete或后退键。
按ESC键时则退出当前的编辑状态。
若要退出文件时,则先按ESC,然后:
q!
,或:
wq,其中:
q!
表示退出不保存修改,:
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格式的参数文件{
名称:
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_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
子像素级配准时划分窗口的数目,比较大才好。
当CC_IN_POS被选择时,就不需要设置。
FC-IN_POS filename
ASCII文件,记录主影像上待匹配的点的坐标。
(不采用)
FC_WINSIZE 6464
匹配窗口的大小。
FC_ACC88
在空间域需要用户设置,在方位向和距离向的搜索偏移量。
FC_INITOFF00|COARSECORR
同CC_INITOFF,COARSECORR表明使用像素级配准后的偏移量。
FC_OSFACTOR16
过采样的内插倍数。
2.4,第四步:
影像重采样。
(1)需要处理步骤
必需:
coregpm估算配准多项式系数。
必需:
resample辅影像重采样。
可选:
filtrange距离向滤波
(2)小步骤参数设置说明
############################################################
分步骤:
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,第五步:
干涉图生成。
(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
用处:
多视处理系数。
说明:
必需。
如果没有多视要求,则默认(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得到的平地相位与原始干涉合成相位之间的间隔差异,由前面得到,正值表示合成相位向右和向上的一个偏移。
说明:
可选
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(相位高程转换)
##########
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- DORIS 操作 流程 详解