StaMPS软件操作作业流程.doc
《StaMPS软件操作作业流程.doc》由会员分享,可在线阅读,更多相关《StaMPS软件操作作业流程.doc(17页珍藏版)》请在咨信网上搜索。
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关键经过调用ROI_PAC软件中mak
2、e_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输出关键结果.raw、.raw.rsc、_parse_lines.out、shift.out、shift.out.rsc作用SAR原始数据预处理,关键完成1)解析数据产品文件,将其转为ROIPAC内部格式;2)
3、填补原始中缺失行;3)调整采样窗口起始时间;4)提取成像参数Roi_prep.pl命令使用roi_prep.pl date OrbitType输入参数影像获取时间(yyyymmdd)、轨道类型(ODR)输出关键结果.slc.rsc、.roi.in、.roi.out作用生成slc参数文件Roi.pl命令使用roi.pl date no 1 no输入参数影像获取时间(yyyymmdd)输出关键结果.slc作用生成slc图像roipac2doris命令使用Roipac2doris date s输入参数影像获取时间(yyyymmdd)输出关键结果master.res或slave.res作用生成Dori
4、s处理所需影像信息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到主影像成像文件夹中。还要注意精裁裁剪范围是在粗裁后影像
5、中得范围不是在初始成像后主影像行列。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命令,对全部辅影像根据爱好范围成像。运行过程为逐一进行辅影像目录,实施st
6、ep_slc_ers命令成像。step_slc_ers使用方法同1.2./注意其它辅影像成像时候是根据粗裁范围裁剪,她不进行精裁,因为主影像又进行了精裁,通常辅影像比较大。对应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裁
7、剪后主影像再次读取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_D
8、IR/ERS1 /精密轨道文件路径M_ORB_INTERVAL 1 /星历时间间隔(s)M_ORB_EXTRATIME 6 /Time in seconds before first 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_O
9、RB_EXTRATIME 6 /Time in seconds before first and last line to output ephemerides.输出结果提出精密轨道信息保留在slave.res中。作用提取辅影像精密轨道信息2.2 主辅影像配准2.2.1 主辅影像粗配准make_coarse:进入每幅辅影像对应文件夹,在每个文件夹内实施step_coarse命令,(coarse.dorisin)进行主辅影像一级轨道配准和二级像素级配准。Step_coarse命令实施Doris软件COARSEORB和COARSECORR两步骤。COARSEORB关键参数不需要参数,依靠读取主辅影
10、像master.res和slave.res来完成输出结果轨道配准结果保留在coreg.out文件中作用基于主辅影像轨道数据配准,精度达30个像元COARSECORR关键参数CC_METHOD magfft /相关系数计算方法,magfft是先使用fft变换到频率域,然后计算相关系数。magspace是直接在空间域计算相关系数CC_NWIN 21 /像素级配准划分窗口数目,大于5. CC_WINSIZE 1024 256 /匹配窗口大小,处理时候被改变为奇数。CC_INITOFF orbit /利用轨道计算偏移量作为初始偏移量输出结果计算出Coarse_correlation_translati
11、on_lines和coarse_correlation_translation_pixels,保留在coreg.out文件中。作用利用相关系数法实现像素级配准。2.2.2 主辅影像精配准make_coreg:精配准,精配准过程中全部基线距小于设定bperp_max值直接和主影像配准,基线距比较大和离它最近三个辅影像配准。然后利用Hooper提出配准方法,利用中间过渡方法实现全部主辅影像配准。make_coreg命令运行Doris软件中FINE和COREGPM步骤,即三级子像素精配准和估量配准多项式系数。FINE关键参数FC_METHOD oversample /精配准计算方法,过采样处理FC_
12、IN_POS fc_pos.in / ASCII文件,统计主影像上待配准点点坐标FC_WINSIZE 64 64 / 匹配窗口大小FC_ACC 8 8 /在方位向和距离向搜索偏移量FC_INITOFF coarsecorr /使用像素级配准后偏移量FC_OSFACTOR 32 / 过采样内插倍数输出结果配准结果保留在INSAR_date文件夹下coreg文件夹中coreg.out和CPM_DATA.n1.n2文件里作用计算主辅影像每个配准点偏移量COREGPM关键参数CPM_THRESHOLD 0.30 /经配准后生成大量匹配点,当相关系数大于CPM_THRESHOLD时,才参与平差运算。阈值
13、大小和File中窗口大小相关,小窗口情况下阈值趋于1.0,阈值越大越好。CPM_DEGREE 2 /多项式阶数,2阶即可。CPM_WEIGHT bamler / 最小二乘中权矩阵方法选择,bamler:权矩阵为相关系数2倍,linear:权矩阵为相关系数1倍。CPM_MAXITER 8000 /最小二乘运算最大迭代次数输出结果计算出方位向和距离向配准多项式系数,输出结果保留在coreg文件夹coreg.out文件里。作用计算配准多项式系数2.3 主影像强度图模拟和地形相位计算(平坦地域地形可不实施 )cd INSAR_master_date/在INSAR_master_date文件夹,修改ti
14、ming.dorisin文件参数step_master_timing:模拟主影像强度图,计算DEM和主影像Timing Error。make_dems:进入到每个辅影像目录,实施step_dem命令,计算每个干涉正确地形相位。step_master_timing运行Doris软件M_SIMAMP和M_TIMING步骤:M_SIMAMP关键参数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 360
15、1 / DEM行数和列数SAM_IN_DELTA 0. / DEM采样间隔 SAM_IN_UL 39.99583333 115.9958333 / 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.logNcols
16、DEM和NrowsDEM)作用模拟主影像幅度图M_TIMING关键参数MTE_METHOD magfft / Magfft是先使用fft变换到频率域,然后计算相关系数。Magspace是直接在空间域计算相关系数。Magfft计算快于magspace。MTE_NWIN 30 / 划分窗口数目MTE_INITOFF 0 0 /初始偏移量MTE_WINSIZE 4096 2048 / 匹配窗口大小输出结果模拟幅度图和主影像方位向和距离向Timing Error:Master_azimuth_timing_error和Master_range_timing_error模拟图和主影像粗配准偏移量Coar
17、se_correlation_translation_lines和Coarse_correlation_translation_pixels;作用计算模拟幅度图和主影像Timing Error和粗配准偏移量。/DEM偏移量计算。step_dem运行Doris软件中COMPREFDEM步骤:COMPREFDEM关键参数CRD_METHOD trilinear /求取地形相位方法,有nearest和trilinear两种。因为是粗DEM,和影像分辨率不一致,需要过采样到同一分辨率下。CRD_INCLUDE_FE OFF / phase w.r.t. ellipsoidCRD_OUT_FILE re
18、fdem_1l.raw /由DEM反算得到地形相位文件CRD_OUT_DEM_LP dem_radar.raw /原DEM得到目前雷达坐标系统下新DEM(已采样)输出结果地形相位文件refdem_1l.raw和雷达坐标系下DEM :dem_radar.raw(在各辅影像文件夹中), dem_radar_i.raw(在主影像文件夹中)这多个文件和主影像行列号一样,float格式。作用计算地形相位2.4 辅影像重采样make_resample:对辅影像根据精配准算出多项式来重采样。进入到每个辅影像目录,实施step_resample。(resample.dorisin,自动加入RS_DBOW 10
19、00 10999 1500 3019)step_resample运行Doris软件中RESAMPLE步骤:RESAMPLE关键参数RS_METHOD rc12p /重采样内插核RS_OUT_FILE slave_res.slc /重采样后生成文件名RS_OUT_FORMAT cr4 /输出文件数据类型输出结果重采样后辅影像slave_res.slc作用辅影像重采样2.5差分干涉处理make_ifgs:主辅影像进行干涉处理,去除平地相位和地形相位。进入到每个辅影像目录,运行step_ifg命令。(integram.dorisin)(make_ifgs_nodem对于平坦地域不用去地形)step_
20、ifg运行Doris软件中SUBTRREFDEM、COMPREFPHA和SUBTRREFPHA步骤,即 生成干涉图,利用轨道来计算平地相位,减去平地相位,减去由粗DEM反算得到地形相位:inter关键参数INT_OUT_CINT cint.raw /生成干涉图文件,为复数float型INT_MULTILOOK 1 1 / 多视处理系数输出结果干涉图:cint.raw作用干涉处理,生成干涉图COMPREFPHA关键参数FE_METHOD porbits /采取计算平地相位方法,精轨FE_DEGREE 3 /位置和平地相位之间多项式系数FE_NPOINTS 201 /最小二乘法来计算平地相位所用来
21、拟合点数输出结果平地相位终端实施过程存于Interferogram.out中作用依据精轨信息计算平地相位SUBTRREFPHA(SUBTRREFDEM)关键参数SRP_METHOD polynomial /采取减去平地相位方法SRP_OUT_CINT cint.minrefpha.raw /去除平地效应后干涉图输出结果去除平地相位干涉图:cint.minrefpha.raw(去地形后没有了)作用去除平地相位影响SUBTRREFDEM关键参数SRD_OUT_CINT cint.minrefdem.raw/去除地形相位后复干涉图文件SRD_OFFSET 0 0 /基于粗DEM得到地形相位和原始干涉
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- StaMPS 软件 操作 作业 流程
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【天****】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【天****】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。