StaMPS软件操作的流程纲要纲要.docx
《StaMPS软件操作的流程纲要纲要.docx》由会员分享,可在线阅读,更多相关《StaMPS软件操作的流程纲要纲要.docx(40页珍藏版)》请在冰豆网上搜索。
StaMPS软件操作的流程纲要纲要
注意:
冒号后边都是解说,假如步骤。
没说改正目录就向来是前一步的目录。
每个命令都有一个log对应
1成像办理(ROI_PAC)
SAR影像原始数据准备
成立SLC文件夹,将yyyymmdd文件夹,每景
SAR原数据拷贝到SLC目录里,为每景影像依据获取时间分别成立
SAR影像原始数据对应的数据和头文件分别命名为IMAGERYyymmdd
SARLEADERyyyymmdd。
主影像成像(对应的是0级)
依据主影像选用原则选用主影像,进入到主影像文件夹
master_date,关于ERS-1/2数据,
运转step_slc_ers
命令成像。
step_slc_ers
主要经过调用
ROI_PAC软件中的make_raw.pl,
roi_prep.pl,
命令,及其StaMPS自己的roipac2doris
命令来达成。
step_slc_ers命令介绍
命令使用
orbit_type
leader_file
date
rs
输入参数
轨道种类(ODR)、数据头文件
(SARLEADERyyyymmdd)和影像获取时间
yyyymmdd
输出主要结果
.raw
、
、
、
、
作用
SAR原始数据预办理,主要达成
1)分析数
据产品文件,将其转为
ROI_PAC的内部格
式;2)填充原始中缺失的行;
3)调整采样
窗口的开端时间;
4)提取成像参数
命令使用
dateOrbitType
输入参数
影像获取时间(yyyymmdd)、轨道种类(ODR)
输出主要结果
、
、
作用
生成slc参数文件
命令使用
dateno1no
输入参数
影像获取时间(yyyymmdd)
输出主要结果
.slc
作用
生成slc图像
roipac2dori
命令使用
Roipac2dorisdate[s]
s
输入参数
影像获取时间(yyyymmdd)
输出主要结果
或
作用
生成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
$SAR_ODR_DIR/ERS1//精细轨道文件路径
S_ORB_INTERVAL
1
//星历时间间隔(s)
S_ORB_EXTRATIME6
//Timeinsecondsbeforefirstandlastline
tooutputephemerides.
输出结果
提出的精细轨道信息保留在
中。
作用
提取辅影像精细轨道信息
主辅影像配准
主辅影像粗配准
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:
权矩阵为有关系数的
1
倍。
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
mt_prep
主要参数
da_thresh
//
幅度离差阈值,一般取值
,小于此幅
度离差的选为
PS候选点,默以为
rg_patches
//
距离向上分块,默以为
1
az_patches
//
方向向上分块,默以为
1
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
image
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'];
%1
float
value
per
PS
candidate
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),
改正某个参数值。
参数结果保留在文件里。
//在INSAR
masterdate
下,运转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