StaMPS软件操作流程

注意：冒号后面都是解释，如果步骤。没说更改目录就一直是前一步的目录。每个命令都有一个log对应1成像处理（ROI_PAC）1.1SAR影像原始数据准备建立SLC文件夹，将SAR原数据拷贝到SLC目录里，为每景影像按照获取时间分别建立yyyymmdd文件夹，每景SAR影像原始数据对应的数据和头文件分别命名为IMAGERYyymmdd和SARLEADERyyyymmdd。1.2主影像成像(对应的是0级)按照主影像选取原则选取主影像，进入到主影像文件夹master_date，对于ERS-1/2数据，运行step_slc_ers命令成像。step_slc_ers主要通过调用ROI_PAC软件中的make_raw.pl,roi_prep.pl,roi.pl命令，及其StaMPS本身的roipac2doris命令来完成。step_slc_ers命令介绍Step_slc_ersMake_raw.pl命令使用make_raw.plorbit_typeleader_filedate输入参数轨道类型（ODR）、数据头文件（SARLEADERyyyymmdd）和影像获取时间yyyymmdd输出主要结果date.raw、date.raw.rsc、date_parse_lines.out、shift.out、shift.out.rsc作用SAR原始数据预处理，主要完成1）解析数据产品文件，将其转为ROI＿PAC的内部格式；2）填补原始中缺失的行；3）调整采样窗口的起始时间；4）提取成像参数Roi_prep.pl命令使用roi_prep.pldateOrbitType输入参数影像获取时间（yyyymmdd）、轨道类型（ODR）输出主要结果date.slc.rsc、date.roi.in、date.roi.out作用生成slc参数文件Roi.pl命令使用roi.pldateno1no输入参数影像获取时间（yyyymmdd）输出主要结果date.slc作用生成slc图像roipac2doris命令使用Roipac2dorisdate[s]输入参数影像获取时间（yyyymmdd）输出主要结果master.res或slave.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.in到主影像成像的文件夹中。还要注意精裁的裁剪范围是在粗裁后的影像中得范围不是在初始成像后的主影像的行列。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等文件拷入此文件夹作用确定裁剪区域，建立主影像，确定辅影像列表1.5辅影像成像在SLC目录下，运行make_slcs_ers命令，对所有辅影像按照兴趣范围成像。运行过程为逐个进行辅影像目录，执行step_slc_ers命令成像。step_slc_ers用法同1.2.//注意其他辅影像成像的时候是按照粗裁的范围裁剪的，他不进行精裁，由于主影像又进行了精裁，一般辅影像比较大。对应1级产品不用成像，直接仿照说明书，读取即可1建立连接并读取主影像link_slcs_路径（初始数据的）cdmaster_date（即是主影像日期文件夹）//进入主影像日期文件夹stepreadwholeXXX(whereXXXis‘ERS’,‘Envisat’,‘RSAT’,or‘TSX’)第一次对主影像读取，全部读取，不裁剪2裁剪主影像cp$MYSCR/master_crop.in或者从安装文件里考并编辑3裁剪后主影像再次读取step_master_read4读取各个辅影像返回到SLC目录make_read2差分干涉处理（DORIS）2.1提取精密轨道信息（冒号后面的是此步骤的解释说明，不是让你这么操作）在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中。作用提取主影像精密轨道信息step_orbit执行orbit_ers.dorisin命令卡,运行Doris中S_PORBITS步骤：S_PORBITS主要参数S_ORBDIR$SAR_ODR_DIR/ERS1//精密轨道文件路径S_ORB_INTERVAL1//星历时间间隔（s）S_ORB_EXTRATIME6//Timeinsecondsbeforefirstandlastlinetooutputephemerides.输出结果提出的精密轨道信息保存在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文件中作用基于主辅影像轨道数据配准，精度达30个像元COARSECORR主要参数CC_METHODmagfft//相关系数计算方法，magfft是先使用fft变换到频率域，然后计算相关系数。magspace是直接在空间域计算相关系数CC_NWIN21//像素级配准划分窗口的数目，大于5.CC_WINSIZE1024256//匹配窗口的大小，处理的时候被改变为奇数。CC_INITOFForbit//利用轨道计算的偏移量作为初始偏移量输出结果计算出Coarse_correlation_translation_lines与coarse_correlation_translation_pixels，保存在coreg.out文件中。作用利用相关系数法实现像素级配准。2.2.2主辅影像精配准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文件里作用计算主辅影像每个配准点的偏移量COREGPM主要参数CPM_THRESHOLD0.30//经配准后生成大量的匹配点，当相关系数大于CPM_THRESHOLD时，才参与平差运算。阈值的大小与File中窗口的大小有关，小窗口情况下阈值趋于1.0，阈值越大越好。CPM_DEGREE2//多项式的阶数，2阶即可。CPM_WEIGHTbamler//最小二乘中权矩阵方法选择，bamler：权矩阵为相关系数的2倍，linear：权矩阵为相关系数的1倍。CPM_MAXITER8000//最小二乘运算的最大迭代次数输出结果计算出方位向和距离向配准多项式系数，输出的结果保存在coreg文件夹的coreg.out文件里。作用计算配准多项式系数2.3主影像强度图模拟与地形相位计算(平坦地区地形可不执行)cdINSAR_master_date//在INSAR_master_date文件夹，修改timing.dorisin文件参数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_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//求取地形