StaMPS软件操作的流程纲要纲要Word文件下载.docx
- 文档编号:21994695
- 上传时间:2023-02-02
- 格式:DOCX
- 页数:40
- 大小:34.39KB
StaMPS软件操作的流程纲要纲要Word文件下载.docx
《StaMPS软件操作的流程纲要纲要Word文件下载.docx》由会员分享,可在线阅读,更多相关《StaMPS软件操作的流程纲要纲要Word文件下载.docx(40页珍藏版)》请在冰豆网上搜索。
dateOrbitType
影像获取时间(yyyymmdd)、轨道种类(ODR)
生成slc参数文件
dateno1no
影像获取时间(yyyymmdd)
.slc
生成slc图像
roipac2dori
Roipac2dorisdate[s]
s
或
生成Doris
办理所需的影像信息
辅影像成像
选择主影像兴趣地区从头成像Iyao
观察主影像成像图像,确立研究地区,依据研究地区界限扩展1000像元的范围来编写
roi.proc。
从头运转step_slc_ers,依据兴趣地区从头成像。
step_slc_ers命令用法同
1.2.//粗裁。
直接在SLC文件夹下生成有这个文件,在里面改正,别忘掉去掉#号。
成立主影像
精裁,察看主影像成像图像,编写文件,再次确立裁剪地区。
运转
step_master_setup命令成立主影像。
//从ROI_PAC_SCR中复制到主影像
成像的文件夹中。
还要注意精裁的裁剪范围是在粗裁后的影像中得范围不是在初始成像后的主影像的队列。
Step_master_set所需参数文件
up
生成结果
辅影像列表:
在与SLC
同目录下成立文件夹
INSAR_master_date,将
等文件拷入此文件夹
确立裁剪地区,成立主影像,确立辅影像列表
在SLC目录下,运转make_slcs_ers命令,对所有辅影像依据兴趣范围成像。
运转过程为逐一进行辅影像目录,履行step_slc_ers命令成像。
step_slc_ers用法同1.2.
注意其余辅影像成像的时候是依据粗裁的范围裁剪的,他不进行精裁,因为主影像又进行了精裁,一般辅影像比较大。
对应1级产品不用成像,直接模仿说明书,读取即可
成立连结并读取主影像
link_slcs_
路径
(初始数据的)
cdmaster_date
(即是主影像日期文件夹)
//进入主影像日期文件夹
stepread
whole
XXX(whereXXXis‘ERS’,‘Envisat’,‘RSAT’,or‘TSX’)
4第一次对主影像读取,所有读取,不裁剪
5裁剪主影像
6cp$MY或许从安装文件里考
7并编写
8裁剪后主影像再次读取
9step_master_read
10读取各个辅影像
返回到SLC目录
make_read
2差分干预办理(DORIS)
提取精细轨道信息(冒号后边的是此步骤的解说说明,不是让你
这么操作)
insarmaster文件夹下
step_master_orbit_ODR:
主影像精细轨道信息提取。
//没有精轨不运转这步,直接下步
make_orbits:
为每幅辅影像成立相应的文件夹,在每个文件夹内履行
step_orbit命令,读取
相应辅影像精细轨道信息。
step_master_orbit_ODR运转Doris中的M_PORBITS步骤:
M_PORBITS
主要参数
M_ORBDIR
$SAR_ODR_DIR/ERS1//精细轨道文件路径
M_ORB_INTERVAL
1
//星历时间间隔(s)
M_ORB_EXTRATIME6
//Timeinsecondsbeforefirstandlastline
tooutputephemerides.
输出结果
提出的精细轨道信息保留在
中。
提取主影像精细轨道信息
step_orbit履行命令卡,运转Doris中S_PORBITS步骤:
S_PORBITS
S_ORBDIR
S_ORB_INTERVAL
S_ORB_EXTRATIME6
提取辅影像精细轨道信息
主辅影像配准
主辅影像粗配准
make_coarse:
进入每幅辅影像相应的文件夹,在每个文件夹内履行step_coarse命令,
)进行主辅影像一级轨道配准和二级像素级配准。
Step_coarse命令履行Doris软件的COARSEORB和COARSECORR两步骤。
COARSEORB
不需要参数,依赖读取主辅影像的的和来达成轨道配准结果保留在文件中
鉴于主辅影像轨道数据配准,精度达30个像元
COARSECORR
CC_METHOD
magfft
//有关系数计算方法,magfft是先使用fft
变换到频次域,而后计算有关系数。
magspace是直接在空间域计算有关
系数
CC_NWIN
21
//像素级配准区分窗口的数量,大于5.
CC_WINSIZE
1024256
//般配窗口的大小,办理的时候被改变
为奇数。
CC_INITOFF
orbit
//利用轨道计算的偏移量作为初始偏移量
计算出Coarse_correlation_translation_lines与
coarse_correlation_translation_pixels,保留在文件中。
作用利用有关系数法实现像素级配准。
主辅影像精配准
make_coreg:
精配准,精配准过程中所有基线距小于设定的
bperp_max值的直接与主影像配
准,基线距比较大的和离它近来的三个辅影像配准。
而后利用用中间过渡的方法实现所有主辅影像的配准。
Hooper提出的配准方法,利
make_coreg命令运转Doris软件中的FINE和COREGPM步骤,即三级子像素精配准和预计配
准多项式系数。
FINE
FC_METHOD
oversample
//
精配准计算方法,过采样办理
FC_IN_POS
//ASCII
文件,记录主影像上待配准
点点的坐标
FC_WINSIZE
6464
般配窗口的大小
FC_ACC88
在方向向和距离向的搜寻偏移量
FC_INITOFF
coarsecorr
使用像素级配准后的偏移量
FC_OSFACTOR
32
过采样的内插倍数
配准结果保留在
INSAR_date
文件夹下
coreg文件夹中的
和
文件里
计算主辅影像每个配准点的偏移量
COREGPM
CPM_THRESHOLD
经配准后生成大批的般配点,当
有关系数大于CPM_THRESHOLD时,才参加平差运算。
阈值的大小与
File
中窗口的大小有关,小窗口状况下阈值趋于
,阈值越大越好。
CPM_DEGREE
2
多项式的阶数,2阶即可。
CPM_WEIGHT
bamler
最小二乘中权矩阵方法选择,
bamler:
权矩阵为有关系数的
2倍,linear:
倍。
CPM_MAXITER
8000
最小二乘运算的最大迭代次数
计算出方向向和距离向配准多项式系数,输出的结果保留在
coreg文
件夹的文件里。
计算配准多项式系数
主影像强度图模拟与地形相位计算(平展地域地形可不履行)
cdINSAR_master_date//在INSAR_master_date文件夹,改正文件参数step_master_timing:
模拟主影像强度图,计算DEM与主影像的TimingError。
make_dems:
进入到每个辅影像目录,履行step_dem命令,计算每个干预对的地形相位。
step_master_timing运转Doris软件的M_SIMAMP和M_TIMING步骤:
M_SIMAMP
主要参数SAM_IN_FORMATI2/r型的//DEM格式,主要包含I2(整),R4
float)和R8(double)
SAM_IN_DEM/sungt/suzhou_test/DEM/suzhou.dem//DEM路径
SAM_IN_SIZE36013601//DEM行数和列数
SAM_IN_DELTA//DEM采样间隔
SAM_IN_UL115.9958333//DEM左上角纬度和精
度,能够在hdr文件里看到,用写字
板翻开可看
M_TIMING
SAM_IN_NODATA-32768/-9999//无值区的取值
SAM_OUT_FILE//模拟主影像幅度图
SAM_OUT_DEM//裁剪的DEM
模拟的主影像幅度图和依据主影像裁剪的DEM
dem_sam_raw或许(队列号在的NcolsDEM和NrowsDEM)
模拟主影像幅度图
MTE_METHODmagfft//Magfft是先使用fft变换到频
率域,而后计算有关系数。
Magspace是直接在空间域计算有关系数。
Magfft计算快于magspace。
MTE_NWIN30//区分窗口的数量
MTE_INITOFF00//初始偏移量
MTE_WINSIZE40962048//般配窗口大小
输出结果模拟幅度图与主影像方向向和距离向的TimingError:
Master_azimuth_timing_error和Master_range_timing_error
模拟图与主影像的粗配准偏移量Coarse_correlation_translation_lines
Coarse_correlation_translation_pixels;
计算模拟幅度图与主影像的
Timing
Error以及粗配准偏移量。
/DEM
偏移量计算。
step_dem运转Doris软件中的COMPREFDEM步骤:
COMPREFDEM
CRD_METHOD
trilinear
求取地形相位方法,有nearest
trilinear两种。
因为是粗DEM,与影像的分辨率不一致,需要过采样到
同一分辨率下。
CRD_INCLUDE_FEOFF
//phasew.r.t.ellipsoid
CRD_OUT_FILE
由DEM反算获取的地形相位文件
CRD_OUT_DEM_LP//
原DEM获取的目前雷达坐标系统下
的新DEM(已采样)
输出结果地形相位文件和雷达坐标系下的DEM:
RESAMPLE
RS_METHODrc12p//
重采样的内插核
RS_OUT_FILE
重采样后生成的文件名
RS_OUT_FORMATcr4
输出文件数据种类
重采样后的辅影像
辅影像重采样
差分干预办理
make_ifgs:
主辅影像进行干预办理,去除平川相位和地形相位。
进入到每个辅影像目录,运转step_ifg命令。
()(make_ifgs_nodem关于平展地域不用去地形)
step_ifg运转Doris软件中的SUBTRREFDEM、COMPREFPHA和SUBTRREFPHA步骤,即生
成干预图,利用轨道来计算平川相位,减去平川相位,减去由粗DEM反算获取的地形相位:
inter
INT_OUT_CINT//
生成的干预图文件,为复数float
型
INT_MULTILOOK11
多视办理系数
干预图:
干预办理,生成干预图
COMPREFPHA
FE_METHOD
porbits
采纳计算平川相位的方法,精轨
FE_DEGREE
3//
地点与平川相位之间的多项式系数
FE_NPOINTS
201//
最小二乘法来计算平川相位所用来拟合的点数
平川相位终端履行过程存于
中
依据精轨信息计算平川相位
SUBTRREFPHA
(SUBTRREFDEM)
SRP_METHOD
polynomial
采纳减去平川相位的方法
SRP_OUT_CINT
去除平川效应后的干预图
去除平川相位的干预图:
(去地形后没有了)
去除平川相位影响
SUBTRREFDEM
SRD_OUT_CINT
去除地形相位后的复干预图文件
SRD_OFFSET
00
鉴于粗DEM获取的地形相位与原始干预合
成相位之间的间隔差别。
去除地形相位后的干预图文件
去除地形相位
地理编码
step_geo:
从雷达坐标系变换到地图坐标系//只在INSAR_masterdate文件夹下随意一个辅文件(就是日期那个文件夹)夹下运转一次就行。
Step_geo运转Doris软件中的GEOCODE步骤:
GEOCODE
GEO_OUT_LAM
经度输出文件
GEO_OUT_PHI
纬度输出文件
像元经纬度信息
和(就生成一个结果在
INSAR文件夹,
就是积累的形变。
)
地理编码,参照坐标系变换
3PS办理
准备
注意小基线不是在此步初选点,是在确立小基线对后再初选点,别弄错了
//命令是mt_prep3250200仍是在INSARmasterdate目录下
mt_prep:
依据设定的幅度离差阈值、距离向方向向分块和重叠范围等参数进行办理,选定初始PS候选点。
mt_prep
在运转过程中还调用了
mt_extract_info
和mt_extract_cands
da_thresh
幅度离差阈值,一般取值
,小于此幅
度离差的选为
PS候选点,默以为
rg_patches
距离向上分块,默以为
az_patches
方向向上分块,默以为
rg_overlap
距离向上块间重叠像元数,默认
50
az_overlap
方向向上块间重叠像元数,默认
200
成立PATCH(分几个块,产生几个文件)文件夹,
(直接在INSAR
master
date目录下)主影像行数和列数
,基线散布bperp.1.in,
辅影像日期day.1.in,
主影
像日期
,干预对日期
,文件头信息
,主影像雷达波长
,视角文件
,幅度定标信息
,DEM参数
分块状况
,
均匀幅度图
,分块范围
,PS
候选点
;
每个候选点的经纬度信息
,每个候选点的高程信息
,每个候选点的
相位信息
等
phname=['
pscands.1.ph'
];
%foreachPScandidate,a
floatcomplexvalueforeachifg
ijname=['
pscands.1.ij'
%
ID#Azimuth#Range#1line
perPScandidate
bperpname=['
bperp.1.in'
%inmeters1line
perslave
image
dayname=['
day.1.in'
YYYYMMDD,1lineperslave
ifgdayname=['
ifgday.1.in'
%YYYYMMDDYYYYMMDD,1line
perifg
masterdayname=['
master_day.1.in'
%YYYYMMDD
llname=['
pscands.1.ll'
%2floatvalues(lonand
lat)perPScandidate
daname=['
pscands.1.da'
%1
float
value
per
PS
candidate
hgtname=['
pscands.1.hgt'
laname=['
look_angle.1.in'
%gridoflookanglevalues
headingname=['
heading.1.in'
%satelliteheading
lambdaname=['
lambda.1.in'
%wavelength
calname=['
calamp.out'
%amplitidecalibrations
widthname=['
width.txt'
%widthofinterferograms
lenname=['
len.txt'
%lengthofinterferograms
提取所需信息,幅度校订,成立分块,选用初始
PS候选点
stamps办理
//在INSAR
masterdate
下先启动matlab
getparm(parmname),无参数时显示所有参数,默认参数从
提取;
setparm(‘parmname’,value),
改正某个参数值。
参数结果保留在文件里。
下,运转stamps
stamps(start_step,end_step,patches_flag
,est_gamma_parm)
Patches_flag
默以为‘y’,假如设置为‘n’把所有数据作为一块办理。
Est_gamma_parm是一个可选参数
passed
tops_est_gamma_quick。
选
y则利用
ps_est_gamma_quick(est_gamma_parm),n利用ps_est_gamma(est
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- StaMPS 软件 操作 流程 纲要