资源预览内容
第1页 / 共15页
第2页 / 共15页
第3页 / 共15页
第4页 / 共15页
第5页 / 共15页
第6页 / 共15页
第7页 / 共15页
第8页 / 共15页
第9页 / 共15页
第10页 / 共15页
亲,该文档总共15页,到这儿已超出免费预览范围,如果喜欢就下载吧!
资源描述
注意:冒号后面都是解释,如果步骤。没说更改目录就一直是前一步的目录。每个命令都有一个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主要通过调用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 dateleader_file输入参数轨道类型(ODR)、数据头文件( SARLEADERyyyymmdd)和影像获取时间 yyyymmdd输出主要结果作用SAR原始数据预处理,主要完成1)解析数据产品文件,将其转为ROIPAC的内部格式;2)填补原始中缺失的行; 3)调整采样窗口的起始时间;4)提取成像参数Roi_prep.命令使用roi_prep.pl date OrbitTypepl输入参数影像获取时间()、轨道类型( ODR)yyyymmdd输出主要结果.slc.rsc、 .roi.in、.roi.out作用生成 slc 参数文件Roi.pl命令使用roi.pl date no 1 no输入参数影像获取时间( yyyymmdd)输出主要结果.slc作用生成 slc 图像roipac2do命令使用Roipac2doris date sris输入参数输出主要结果作用影像获取时间( yyyymmdd)master.res或 slave.res生成 Doris 处理所需的影像信息1.3 选择主影像兴趣区域重新成像I yao观测主影像成像图像,确定研究区域,按照研究区域边界扩展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_s 所需参数文master_crop.inetup件辅影像列表: Make_slcs.list生成结果作用确定裁剪区域,建立主影像,确定辅影像列表1.5 辅影像成像在 SLC目录下,运行 make_slcs_ers 命令,对所有辅影像按照兴趣范围成像。运行过程为逐个进行辅影像目录,执行 step_slc_ers 命令成像。 step_slc_ers 用法同 1.2./ 注意其他辅影像成像的时候是按照粗裁的范围裁剪的,他不进行精裁,由于主影像又进行了精裁,一般辅影像比较大。对应 1 级产品不用成像,直接仿照说明书,读取即可1 建立连接并读取主影像link_slcs_路径(初始数据的)cd master_date(即是主影像日期文件夹)/进入主影像日期文件夹step readwhole XXX (where XXX isERS,Envisat ,RSAT, orTSX)第一次对主影像读取,全部读取,不裁剪2 裁剪主影像cp $MY SCR/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_INTERV AL1/ 星历时间间隔( s)M_ORB_EXTRATIME 6/Time in seconds before first andlast 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 before first andlast line to output ephemerides.输出结果提出的精密轨道信息保存在slave.res中。作用提取辅影像精密轨道信息2.2 主辅影像配准主辅影像粗配准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_WINSIZE1024 256/ 匹配窗口的大小,处理的时候被改变为奇数。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文件,记录主影像上待配准点点的坐标输出结果作用COREGPM主要参数FC_WINSIZE64 64/匹配窗口的大小FC_ACC8 8/在方位向和距离向的搜索偏移量FC_INITOFFcoarsecorr/使用像素级配准后的偏移量FC_OSFACTOR32/过采样的内插倍数配准结果保存在INSAR_date 文件夹下coreg 文件夹中的coreg.out和文件里计算主辅影像每个配准点的偏移量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_
网站客服QQ:2055934822
金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号