《StaMPS软件操作流程》由会员分享,可在线阅读,更多相关《StaMPS软件操作流程(16页珍藏版)》请在金锄头文库上搜索。
1、注意:冒号后面都是解释,如果步骤。没说更改目录就一直是前一步的目录。每个命令都有一个log对应1成像处理(R0I_PAC)1.1 SAR 影像原始数据准备建立 SLC 文件夹,将 SAR 原数据拷贝到 SLC 目录里,为每景影像按照获取时间分别建立 yyyymmdd 文件夹,每景 SAR 影像原始数据对应的数据和头文件分别命名为 IMAGERYyymmdd 和 SARLEADERyyyymmdd。1.2 主影像成像(对应的是 0级)按照主影像选取原则选取主影像,进入到主影像文件夹mas ter_da te,对于ERST/2数据, 运行step_slc_ers命令成像。step_slc_ers主
2、要通过调用ROI_PAC软件中的make_raw.pl, roi_prep.pl, roi.pl 命令,及其 StaMPS 本身的 roipac2doris 命令来完成。step slc ers 命令介绍Step_slc_ersMake_raw.pl命令使用make raw.pl orbit type leader 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 、da te.ro i.out作用生成slc参数文件Roi.pl命令使用roi.pl date no 1 no输入参数影像获取时间(yyyymmdd)输出主要结果da te.slc作用生成slc图像roipac2doris命令使用Roipac2doris date s输入参数影像获取时间(yyyymmdd)输出主要结果master.res 或 slave.res作用
4、生成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.in 到主影像 成像的文
5、件夹中。还要注意精裁的裁剪范围是在粗裁后的影像中得范围不是在初始成像后的主影像的行列。Step_master_set up所需参数文件mas ter crop.in生成结果辅影像列表:Make slcs.list作用确定裁剪区域,建立主影像,确定辅影像列表1.5 辅影像成像在SLC目录下,运行make_slcs_ers命令,对所有辅影像按照兴趣范围成像。运行过程为逐 个进行辅影像目录,执行 step_slc_ers 命令成像。 step_slc_ers 用法同 1.2./注意其他辅影像成像的时候是按照粗裁的范围裁剪的,他不进行精裁,由于主影像又进 行了精裁,一般辅影像比较大。对应 1 级产品不
6、用成像,直接仿照说明书,读取即可1 建立连接并读取主影像link_slcs_路径(初始数据的)cd master_date (即是主影像日期文件夹)/进入主影像日期文件夹s t ep r e a d 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_EXTRATIME 6/Time in
8、seconds before first and last lineto 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 before first and last lineto output e
9、phemerides.输出结果提出的精密轨道信息保存在slave.res中。作用提取辅影像精密轨道信息2.2 主辅影像配准2.2.1 主辅影像粗配准make_coarse:进入每幅辅影像相应的文件夹,在每个文件夹内执行 step_coarse命令,(coarse.dorisin)进行主辅影像一级轨道配准和二级像素级配准。Step coarse 命令执行 Doris 软件的 COARSEORB 和 COARSECORR 两步骤。COARSEORB主要参数不需要参数,依靠读取主辅影像的的master.res和slave.res来兀成输出结果轨道配准结果保存在coreg.out文件中作用基于主辅影像
10、轨道数据配准,精度达30个像兀COARSECORR主要参数CC_METHODmagfft 相关系数计算方法,magfft是先使用fft变换到频率域,然后计算相关系数。magspace是直接在空间域计算相关 系数CC_NWIN21像素级配准划分窗口的数目,大于5.CC_WINSIZE1024 256匹配窗口的大小,处理的时候被改变为奇数。CC INITOFForbit利用轨道计算的偏移量作为初始偏移量输出结果计算出 Coarse_correlation_translationines 与coarse correlation translation pixels, 保存在 coreg.out 文件
11、中。作用利用相关系数法实现像素级配准。2.2.2 主辅影像精配准make_coreg:精配准,精配准过程中所有基线距小于设定的bperp_max值的直接与主影像配 准,基线距比较大的和离它最近的三个辅影像配准。然后利用Hooper提出的配准方法,利 用中间过渡的方法实现所有主辅影像的配准。make_coreg命令运行Doris软件中的FINE和COREGPM步骤,即三级子像素精配准和估计配 准多项式系数。FINE主要参数FC_METH0Doversample/精配准计算方法,过采样处理FC_IN_POSfjpos.in/ ASCII文件,记录主影像上待配准点点的坐标FC_WINSIZE64 6
12、4/ 匹配窗口的大小FC_ACC8 8/在方位向和距离向的搜索偏移量FC_INITOFFcoarsecorr/使用像素级配准后的偏移量FC OSFACTOR32/过采样的内插倍数输出结果配准结果保存在INSAR_da te文件夹下coreg文件夹中的coreg.out和文 件里作用计算主辅影像每个配准点的偏移量COREGPM主要参数CPM_THRESHOLD 0.30/经配准后生成大量的匹配点,当相关系数大于CPM_THRESHOLD时,才参与平差运算。阈值的大小与File 中窗口的大小有关,小窗口情况下阈值趋于1.0,阈值越大越好。 CPM_DEGREE2/多项式的阶数,2阶即可。CPM_W
13、EIGHTbamler/最小一乘中权矩阵方法选择,bamler:权矩阵为相关系数的2倍,linear:权矩阵为相关系数的1 倍。CPM MAXITER8000/最小二乘运算的最大迭代次数输出结果计算出方位向和距离向配准多项式系数,输出的结果保存在coreg文 件夹的coreg.out文件里。作用计算配准多项式系数2.3 主影像强度图模拟与地形相位计算(平坦地区地形可不执行 )cd INSAR_master_date在 INSAR_master_date 文件夹,修改 timing.dorisin 文件参数 step_master_timing :模拟主影像强度图,计算DEM与主影像的Timin
14、g Error。 make_dems:进入到每个辅影像目录,执行step_dem命令,计算每个干涉对的地形相位。step master timing 运行 Doris 软件的 M SIMAMP 和 M TIMING 步骤:M_SIMAMP主要参数SAM_IN_FORMAT I2/r型的 /DEM格式,主要包括I2 (整),R4 (float)和只8 (double)SAM_IN_DEM/sungt/suzhou_test/DEM/suzhou.dem /DEM路径SAM_IN_SIZE3601 3601/ DEM 行数和列数SAM_IN_DELTA 0.000833333/ DEM 采样间隔S
15、AM IN UL3 115.9958333 / DEM左上角纬度和精度,可以在hdr文件里看到,用写字板打开可看 SAM_IN_NODATA -32768/-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主要参数MTE METHODmagfft/ Magfft是先使用fft