StaMPS软件操作作业流程.docx

StaMPS软件操作作业流程.docx

《StaMPS软件操作作业流程.docx》由会员分享，可在线阅读，更多相关《StaMPS软件操作作业流程.docx（26页珍藏版）》请在冰豆网上搜索。

StaMPS软件操作作业流程

StaMPS软件操作作业流程

注意：

冒号后面都是解释，如果步骤。

没说更改目录就一直是前一步的目录。

每个命令都有一个log对应

1成像处理（ROI_PAC）

SAR影像原始数据准备

建立SLC文件夹，将SAR原数据拷贝到SLC目录里，为每景影像按照获取时间分别建立yyyymmdd文件夹，每景SAR影像原始数据对应的数据和头文件分别命名为IMAGERYyymmdd和SARLEADERyyyymmdd。

主影像成像（对应的是0级）

按照主影像选取原则选取主影像，进入到主影像文件夹master_date，对于ERS-1/2数据，运行step_slc_ers命令成像。

,，及其StaMPS本身的roipac2doris命令来完成。

step_slc_ers命令介绍

Step_slc_ers

Make_raw.pl

命令使用

make_raw.plorbit_typeleader_filedate

输入参数

轨道类型（ODR）、数据头文件（SARLEADERyyyymmdd）和影像获取时间yyyymmdd

输出主要结果

.raw、.raw.rsc、_parse_lines.out、shift.out、shift.out.rsc

作用

SAR原始数据预处理，主要完成1）解析数据产品文件，将其转为ROI＿PAC的内部格式；2）填补原始中缺失的行；3）调整采样窗口的起始时间；4）提取成像参数

Roi_prep.pl

命令使用

roi_prep.pldateOrbitType

输入参数

影像获取时间（yyyymmdd）、轨道类型（ODR）

输出主要结果

.slc.rsc、.roi.in、.roi.out

作用

生成slc参数文件

Roi.pl

命令使用

roi.pldateno1no

输入参数

影像获取时间（yyyymmdd）

输出主要结果

.slc

作用

生成slc图像

roipac2doris

命令使用

Roipac2dorisdate[s]

输入参数

影像获取时间（yyyymmdd）

输出主要结果

master.res或slave.res

作用

生成Doris处理所需的影像信息

选择主影像兴趣区域重新成像Iyao

观测主影像成像图像，确定研究区域，。

重新运行step_slc_ers，按照兴趣区域重新成像。

.//粗裁。

直接在SLC文件夹下生成有这个文件，在里面修改，别忘记去掉#号。

建立主影像

精裁，观察主影像成像图像，，再次确定裁剪区域。

运行step_master_setup命令建立主影像。

//。

还要注意精裁的裁剪范围是在粗裁后的影像中得范围不是在初始成像后的主影像的行列。

Step_master_setup

所需参数文件

master_crop.in

生成结果

辅影像列表：

Make_slcs.list

在与SLC同目录下建立文件夹INSAR_master_date,将master_date.slc.rsc、master_crop.slc、dem.dorisin、geocode.dorisn、timing.dorisin、master.res等文件拷入此文件夹

作用

确定裁剪区域，建立主影像，确定辅影像列表

辅影像成像

在SLC目录下，运行make_slcs_ers命令，对所有辅影像按照兴趣范围成像。

运行过程为逐个进行辅影像目录，执行step_slc_ers命令成像。

.

//注意其他辅影像成像的时候是按照粗裁的范围裁剪的，他不进行精裁，由于主影像又进行了精裁，一般辅影像比较大。

对应1级产品不用成像，直接仿照说明书，读取即可

1建立连接并读取主影像

link_slcs_路径（初始数据的）

cdmaster_date（即是主影像日期文件夹）//进入主影像日期文件夹

stepreadwholeXXX（whereXXXis‘ERS’,‘Envisat’,‘RSAT’,or‘TSX’）

第一次对主影像读取，全部读取，不裁剪

2裁剪主影像

cp$MYSCR/或者从安装文件里考

并编辑

3裁剪后主影像再次读取

step_master_read

4读取各个辅影像

返回到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_INTERVAL1//星历时间间隔（s）

M_ORB_EXTRATIME6//Timeinsecondsbeforefirstandlastlinetooutputephemerides.

输出结果

提出的精密轨道信息保存在master.res中。

作用

提取主影像精密轨道信息

运行Doris中S_PORBITS步骤：

S_PORBITS

主要参数

S_ORBDIR$SAR_ODR_DIR/ERS1//精密轨道文件路径

S_ORB_INTERVAL1//星历时间间隔（s）

S_ORB_EXTRATIME6//Timeinsecondsbeforefirstandlastlinetooutputephemerides.

输出结果

提出的精密轨道信息保存在slave.res中。

作用

提取辅影像精密轨道信息

主辅影像配准

主辅影像粗配准

make_coarse：

进入每幅辅影像相应的文件夹，在每个文件夹内执行step_coarse命令，（）进行主辅影像一级轨道配准和二级像素级配准。

Step_coarse命令执行Doris软件的COARSEORB和COARSECORR两步骤。

COARSEORB

主要参数

不需要参数，依靠读取主辅影像的的master.res和slave.res来完成

输出结果

轨道配准结果保存在coreg.out文件中

作用

基于主辅影像轨道数据配准，精度达30个像元

COARSECORR

主要参数

CC_METHODmagfft//相关系数计算方法，magfft是先使用fft变换到频率域，然后计算相关系数。

magspace是直接在空间域计算相关系数

CC_NWIN21//像素级配准划分窗口的数目，大于5.CC_WINSIZE1024256//匹配窗口的大小，处理的时候被改变为奇数。

CC_INITOFForbit//利用轨道计算的偏移量作为初始偏移量

输出结果

计算出Coarse_correlation_translation_lines与

coarse_correlation_translation_pixels，保存在coreg.out文件中。

作用

利用相关系数法实现像素级配准。

主辅影像精配准

make_coreg：

精配准，精配准过程中所有基线距小于设定的bperp_max值的直接与主影像配准，基线距比较大的和离它最近的三个辅影像配准。

然后利用Hooper提出的配准方法，利用中间过渡的方法实现所有主辅影像的配准。

make_coreg命令运行Doris软件中的FINE和COREGPM步骤，即三级子像素精配准和估计配准多项式系数。

FINE

主要参数

FC_METHODoversample//精配准计算方法，过采样处理

FC_IN_POSfc_pos.in//ASCII文件，记录主影像上待配准点点的坐标

FC_WINSIZE6464//匹配窗口的大小

FC_ACC88//在方位向和距离向的搜索偏移量

FC_INITOFFcoarsecorr//使用像素级配准后的偏移量

FC_OSFACTOR32//过采样的内插倍数

输出结果

配准结果保存在INSAR_date文件夹下coreg文件夹中的coreg.out和CPM_DATA.n1.n2文件里

主要参数

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_DELTA0.000833333//DEM采样间隔

SAM_IN_UL39.99583333115.9958333//DEM左上角纬度和精度，可以在hdr文件里看到，用写字板打开可看

SAM_IN_NODATA-32768/-9999//无值区的取值

SAM_OUT_FILEmaster_sam.raw//模拟主影像幅度图

SAM_OUT_DEMdem_sam.raw//裁剪的DEM

输出结果

模拟的主影像幅度图master_sam.raw和根据主影像裁剪的DEMdem_sam_raw或者demcrop.raw（行列号在step_master_timing.log的NcolsDEM和NrowsDEM）

作用

模拟主影像幅度图

M_TIMING

主要参数

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；

作用

计算模拟幅度图与主影像的TimingError以及粗配准偏移量。

/DEM偏移量计算。

step_dem运行Doris软件中的COMPREFDEM步骤：

COMPREFDEM

主要参数

CRD_METHODtrilinear//求取地形相位方法，有nearest和trilinear两种。

由于是粗DEM，与影像的分辨率不一致，需要过采样到同一分辨率下。

CRD_INCLUDE_FEOFF//phasew.r.t.ellipsoid

CRD_OUT_FILErefdem_1l.raw//由DEM反算得到的地形相位文件

CRD_OUT_DEM_LPdem_radar.raw//原DEM得到的当前雷达坐标系统下的新DEM（已采样）

输出结果

地形相位文件refdem_1l.raw和雷达坐标系下的DEM：

dem_radar.raw

（在各辅影像文件夹中），dem_radar_i.ra