1、输入参数轨道类型(ODR)、数据头文件(SARLEADERyyyymmdd)和影像获取时间yyyymmdd输出主要结果.raw、.raw.rsc、_parse_lines.out、shift.out、shift.out.rsc作用SAR原始数据预处理,主要完成1)解析数据产品文件,将其转为ROIPAC的内部格式;2)填补原始中缺失的行;3)调整采样窗口的起始时间;4)提取成像参数Roi_prep.plroi_prep.pl date OrbitType影像获取时间(yyyymmdd)、轨道类型(ODR).slc.rsc、.roi.in、.roi.out生成slc参数文件Roi.plroi.pl
2、 date no 1 no影像获取时间(yyyymmdd).slc生成slc图像roipac2dorisRoipac2doris date smaster.res或slave.res生成Doris处理所需的影像信息 选择主影像兴趣区域重新成像I yao观测主影像成像图像,确定研究区域,。重新运行step_slc_ers,按照兴趣区域重新成像。./粗裁。直接在SLC文件夹下生成有这个文件,在里面修改,别忘记去掉#号。 建立主影像精裁,观察主影像成像图像,再次确定裁剪区域。运行step_master_setup命令建立主影像。/。还要注意精裁的裁剪范围是在粗裁后的影像中得范围不是在初始成像后的主影
3、像的行列。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命令成像。./注
4、意其他辅影像成像的时候是按照粗裁的范围裁剪的,他不进行精裁,由于主影像又进行了精裁,一般辅影像比较大。对应1级产品不用成像,直接仿照说明书,读取即可1建立连接并读取主影像link_slcs_路径 (初始数据的)cd master_date(即是主影像日期文件夹)/进入主影像日期文件夹step read whole XXX (where XXX is ERS, Envisat, RSAT, or TSX)第一次对主影像读取,全部读取,不裁剪2裁剪主影像cp $MY SCR/ 或者从安装文件里考并编辑3裁剪后主影像再次读取step_master_read4读取各个辅影像返回到SLC目录make_r
5、ead2 差分干涉处理(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_
6、ORB_EXTRATIME 6 /Time in seconds before first and last line to output ephemerides.输出结果提出的精密轨道信息保存在master.res中。提取主影像精密轨道信息,运行Doris中S_PORBITS步骤:S_PORBITSS_ORBDIR $SAR_ODR_DIR/ERS1 /精密轨道文件路径S_ORB_INTERVAL 1 /星历时间间隔(s)S_ORB_EXTRATIME 6 /Time in seconds before first and last line to output ephemerides.提出
7、的精密轨道信息保存在slave.res中。提取辅影像精密轨道信息 主辅影像配准 主辅影像粗配准make_coarse:进入每幅辅影像相应的文件夹,在每个文件夹内执行step_coarse命令,()进行主辅影像一级轨道配准和二级像素级配准。Step_coarse命令执行Doris软件的COARSEORB和COARSECORR两步骤。COARSEORB不需要参数,依靠读取主辅影像的的master.res和slave.res来完成轨道配准结果保存在coreg.out文件中基于主辅影像轨道数据配准,精度达30个像元COARSECORRCC_METHOD magfft /相关系数计算方法,magfft是
8、先使用fft变换到频率域,然后计算相关系数。magspace是直接在空间域计算相关系数CC_NWIN 21 /像素级配准划分窗口的数目,大于5. CC_WINSIZE 1024 256 /匹配窗口的大小,处理的时候被改变为奇数。CC_INITOFF orbit /利用轨道计算的偏移量作为初始偏移量计算出Coarse_correlation_translation_lines与coarse_correlation_translation_pixels,保存在coreg.out文件中。利用相关系数法实现像素级配准。 主辅影像精配准make_coreg:精配准,精配准过程中所有基线距小于设定的bpe
9、rp_max值的直接与主影像配准,基线距比较大的和离它最近的三个辅影像配准。然后利用Hooper提出的配准方法,利用中间过渡的方法实现所有主辅影像的配准。make_coreg命令运行Doris软件中的FINE和COREGPM步骤,即三级子像素精配准和估计配准多项式系数。FINEFC_METHOD oversample /精配准计算方法,过采样处理FC_IN_POS fc_pos.in / ASCII文件,记录主影像上待配准点点的坐标FC_WINSIZE 64 64 / 匹配窗口的大小FC_ACC 8 8 /在方位向和距离向的搜索偏移量FC_INITOFF coarsecorr /使用像素级配准
10、后的偏移量FC_OSFACTOR 32 / 过采样的内插倍数配准结果保存在INSAR_date文件夹下coreg文件夹中的coreg.out和CPM_DATA.n1.n2文件里SAM_IN_FORMAT I2/r型的 /DEM格式,主要包括I2(整),R4(float)和R8(double)SAM_IN_DEM /sungt/suzhou_test/DEM/suzhou.dem /DEM路径SAM_IN_SIZE 3601 3601 / DEM行数和列数SAM_IN_DELTA 0.000833333 / DEM采样间隔 SAM_IN_UL 39.99583333 115.9958333 /
11、DEM左上角纬度和精度,可以在hdr文件里看到,用写字板打开可看SAM_IN_NODATA -32768/-9999 /无值区的取值SAM_OUT_FILE master_sam.raw /模拟主影像幅度图SAM_OUT_DEM dem_sam.raw /裁剪的DEM模拟的主影像幅度图master_sam.raw和根据主影像裁剪的DEM dem_sam_raw 或者demcrop.raw(行列号在step_master_timing.log的NcolsDEM和NrowsDEM)模拟主影像幅度图M_TIMINGMTE_METHOD magfft / Magfft是先使用fft变换到频率域,然后计
12、算相关系数。Magspace是直接在空间域计算相关系数。Magfft计算快于magspace。MTE_NWIN 30 / 划分窗口的数目MTE_INITOFF 0 0 /初始偏移量MTE_WINSIZE 4096 2048 / 匹配窗口大小模拟幅度图与主影像方位向和距离向的Timing Error:Master_azimuth_timing_error和Master_range_timing_error模拟图与主影像的粗配准偏移量Coarse_correlation_translation_lines与Coarse_correlation_translation_pixels;计算模拟幅度图与
13、主影像的Timing Error以及粗配准偏移量。/DEM偏移量计算。step_dem运行Doris软件中的COMPREFDEM步骤:COMPREFDEMCRD_METHOD trilinear /求取地形相位方法,有nearest和trilinear两种。由于是粗DEM,与影像的分辨率不一致,需要过采样到同一分辨率下。CRD_INCLUDE_FE OFF / phase w.r.t. ellipsoidCRD_OUT_FILE refdem_1l.raw /由DEM反算得到的地形相位文件CRD_OUT_DEM_LP dem_radar.raw /原DEM得到的当前雷达坐标系统下的新DEM(已采样)地形相位文件refdem_1l.raw和雷达坐标系下的DEM :dem_radar.raw(在各辅影像文件夹中), dem_radar_i.ra
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1