《StaMPS软件操作流程》由会员分享，可在线阅读，更多相关《StaMPS软件操作流程（25页珍藏版）》请在金锄头文库上搜索。
1、注意：冒号后面都是解释，如果步骤。没说更改目录就一直是前一步的目录。每个命令都有一个 log 对应1成像处理（ROI_PAC ）1.1 SAR 影像原始数据准备建立SLC文件夹，将SAR原数据拷贝到SLC目录里，为每景影像按照获取时 间分别建立 yyyymmdd 文件夹，每景 SAR 影像原始数据对应的数据和头文件 分别命名为 IMAGERYyymmdd 和 SARLEADERyyyymmdd。1.2 主影像成像（对应的是 0级）按照主影像选取原则选取主影像，进入到主影像文件夹master_date，对于 ERS-1/2数据，运行step_slc_ers命令成像。step_slc_ers主要通
2、过调用 ROI_PAC 软件中的 make_raw.pl, roi_prep.pl, roi.pl 命令，及其 StaMPS 本身的 roipac2doris 命令来完成。step_slc_ers 命令介绍Step_sl c_ersMake_raw .pl命令使用make_raw.plorbit_typeleader_file date输入参数轨道类型（ODR ）、数据头文件 （SARLEADERyyyymmdd）和影像获 取时间yyyymmdd输出主要结 果作用SAR原始数据预处理，主要完成 1）解析数据产品文件，将其转为 ROI_PAC的内部格式；2）填补 原始中缺失的行；3）调整米样窗
3、口的起始时间；4）提取成像参数Roi_prep .pl命令使用roi_prep.pl date OrbitType输入参数影像获取时间（yyyymmdd）、轨 道类型（ODR）输出主要结 果da te.slc.rsc、da te.roi.in、date.roi. out作用生成slc参数文件Roi.pl命令使用roi.pl date no 1 no输入参数影像获取时间（yyyymmdd）输出主要结 果da te.slc作用生成slc图像roipac2d oris命令使用Roipac2doris date s输入参数影像获取时间（yyyymmdd）输出主要结 果master.res 或 slav
4、e.res作用生成Doris处理所需的影像信息1.3 选择主影像兴趣区域重新成像 Iyao 观测主影像成像图像，确定研究区域，按照研究区域边界扩展 1000 像元的范围来编辑roi.proc。重新运行step_slc_ers，按照兴趣区域重新成像。 step_slc_ers 命令用法同 1.2./粗裁。直接在 SLC 文件夹下生成有这个文 件，在里面修改，别忘记去掉#号。1.4 建立主影像精裁，观察主影像成像图像，编辑 master_crop.in 文件，再次确定裁剪区 域。运行 step_master_setup 命令建立主影像。 /从 ROI_PAC_SCR 中复制 master_crop
5、.in 到主影像成像的文件夹中。还要注意精裁的裁剪范围是在粗裁后的影像中得范围不是在初始成像后的主影像的行列。Step_master_setup所需参数 文件mas ter_crop.in生成结果辅影像列表：Make slcs.list作用确定裁剪区域，建立主影像，确定辅 影像列表1.5 辅影像成像在 SLC 目录下，运行 make_slcs_ers 命令，对所有辅影像按照兴趣范围成 像。运行过程为逐个进行辅影像目录，执行 step_slc_ers 命令成像。 step_slc_ers 用法同 1.2./注意其他辅影像成像的时候是按照粗裁的范围裁剪的，他不进行精裁， 由于主影像又进行了精裁，一
6、般辅影像比较大。对应 1 级产品不用成像，直接仿照说明书，读取即可1 建立连接并读取主影像link_slcs_路径（初始数据的）cd master_date （即是主影像日期文件夹）/进入主影像日期文件夹step read whole XXX （where XXX is ERS, Envisat, RSAT, or TSX） 第一次对主影像读取，全部读取，不裁剪2 裁剪主影像cp $MY SCR/master_crop.in 或者从安装文件里考并编辑3 裁剪后主影像再次读取step_master_read4 读取各个辅影像返回到SLC目录make_read2 差分干涉处理（ DORIS）2.1
7、 提取精密轨道信息（冒号后面的是此步骤的解释说明，不是让你这么操 作）在 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_EXTRATI
8、ME 6/Time in seconds beforefirst and last line to output ephemerides.输出结果提出的精密轨道信息保存在master.res中。作用提取主影像精密轨道信息step_orbit 执行 orbit_ers.dorisin 命令卡,运行 Doris 中 S_PORBITS 步骤:S_PORBITS主要参数S_ORBDIR$SAR_ODR_DIR/ERS1 精密轨道文件路径S_ORB_INTERVAL 1星历时间间隔（s）S_ORB_EXTRATIME 6/Time in seconds beforefirst and last lin
9、e to output ephemerides.输出结果提出的精密轨道信息保存在slave.res中。作用提取辅影像精密轨道信息2.2 主辅影像配准2.2.1 主辅影像粗配准make_coarse :进入每幅辅影像相应的文件夹，在每个文件夹内执行 step_coarse命令，(coarse.dorisin)进行主辅影像一级轨道配准和二级像素级配准。Step_coarse 命令执行 Doris 软件的 COARSEORB 和 COARSECORR 两步骤。COARSEORB主要参数不需要参 数，依靠读取主辅影像的的master.res和slave.res来兀成输出结果轨道配准结果保存在coreg
10、.out文件中作用基于主辅影像轨道数据配准，精度达30个像兀COARSECORR主要参数CC_METHODmagfft 相关系数计算方法，magfft是先使用fft变换到频率域，然后计算相关系数。magspace是直接在空间域计算相关 系数CC_NWIN21像素级配准划分窗口的数目，大于5.CC_WINSIZE1024 256匹配窗口的大小，处理的时候被改变为奇数。CC INITOFForbit利用轨道计算的偏移量作为初始偏移量输出结果计算出 Coarse_correlation_translation_lines 与 coarse_correlation_translation_pixels
11、 ， 保存在 coreg.out 文件中。作用利用相关系数法实现像素级配准。2.2.2 主辅影像精配准make_coreg :精配准，精配准过程中所有基线距小于设定的bperp_max值的 直接与主影像配准，基线距比较大的和离它最近的三个辅影像配准。然后 利用Hooper提出的配准方法，利用中间过渡的方法实现所有主辅影像的配 准。make_coreg命令运行Doris软件中的FI NE和COREGPM步骤，即三级子像素精配准和估计配准多项式系数。FINE主要参数FC_METHODoversample/精配准计算方法，过采样处理FC_IN_POSfc_pos.in/ ASCII文件，记录主影像上
12、待配准点点的坐标FC_WINSIZE64 64/ 匹配窗口的大小FC_ACC8 8/在方位向和距离向的搜索偏移量FC_INITOFFcoarsecorr/使用像素级配准后的偏移量FC OSFACTOR32/过采样的内插倍数输出结果配准结果保存在INSAR_date文件夹下coreg文件夹中的coreg.out和文 件里作用计算主辅影像每个配准点的偏移量COREGPM主要参数CPM_THRESHOLD 0.30/经配准后生成大量的匹配点，当相关系数大于CPM_THRESHOLD时，才参与平差运算。阈值的大小与File 中窗口的大小有关，小窗口情况下阈值趋于1.0，阈值越大越好。 CPM_DEGR
13、EE2/多项式的阶数，2阶即可。CPM_WEIGHTbamler/最小一乘中权矩阵方法选择，bamler:权矩阵为相关系数的2倍，linear:权矩阵为相关系数的1 倍。CPM MAXITER8000/最小二乘运算的最大迭代次数输出结果计算出方位向和距离向配准多项式系数，输出的结 果保存在coreg文件夹的coreg.out文件里。作用计算配准多项式系数2.3 主影像强度图模拟与地形相位计算（平坦地区地形可不执行 ）cd INSAR_master_date在 INSAR_master_date 文件夹，修改 timing.dorisin 文件参数step_master_timing:模拟主影像
14、强度图，计算DEM与主影像的TimingError。make_dems:进入到每个辅影像目录，执行step_dem命令，计算每个干涉 对的地形相位。step_master_timing 运行 Doris 软件的 M_SIMAMP 和 M_TIMING 步骤:M_SIMAMP主要参数SAM_IN_FORMAT I2/r 型的 /DEM 格式，主要 包括 12 （整），R4（float）和只8 （double） SAM_IN_DEM/sungt/suzhou_test/DEM/suzhou.dem /DEM 路径SAM_IN_SIZE3601 3601/ DEM 行数和列数SAM_IN_DELTA
15、 0.000833333/ DEM 采样间隔SAM IN UL3 115.9958333 / DEM左上角纬度和精度，可以在hdr文件 里看到，用写字板打开可 看SAM_IN_NODATA -3276弘9999无值区的取值SAM_OUT_FILEmaster_sam.raw 模拟主影像幅度图SAM_OUT_DEMdem_sam.raw裁剪的DEM输出结果模拟的主影像幅度图master_sam.raw和根据主影像裁剪的 DEM dem_sam_raw 或者 demcrop.raw （行列号在 step master timing.log 的 NcolsDEM 和 NrowsDEM）作用模拟主影像幅度图M_TIMING
