1、StaMPS软件操作的流程纲要纲要注意:冒号后边都是解说,假如步骤。没说改正目录就向来是前一步的目录 。每个命令都有一个 log对应1成像办理(ROI_PAC)SAR影像原始数据准备成立SLC文件夹,将yyyymmdd文件夹,每景SAR原数据拷贝到 SLC目录里,为每景影像依据获取时间分别成立SAR影像原始数据对应的数据和头文件分别命名为 IMAGERYyymmddSARLEADERyyyymmdd。主影像成像(对应的是0级)依据主影像选用原则选用主影像,进入到主影像文件夹master_date,关于ERS-1/2数据,运转step_slc_ers命令成像。step_slc_ers主要经过调用
2、ROI_PAC软件中的make_raw.pl,roi_prep.pl,命令,及其StaMPS自己的roipac2doris命令来达成。step_slc_ers命令介绍命令使用orbit_typeleader_filedaters输入参数轨道种类(ODR)、数据头文件(SARLEADERyyyymmdd)和影像获取时间yyyymmdd输出主要结果.raw、作用SAR原始数据预办理,主要达成1)分析数据产品文件,将其转为ROIPAC的内部格式;2)填充原始中缺失的行;3)调整采样窗口的开端时间;4)提取成像参数命令使用dateOrbitType输入参数影像获取时间(yyyymmdd)、轨道种类(O
3、DR)输出主要结果、作用生成slc参数文件命令使用dateno1no输入参数影像获取时间(yyyymmdd)输出主要结果.slc作用生成slc图像roipac2dori命令使用Roipac2dorisdatess输入参数影像获取时间(yyyymmdd)输出主要结果或作用生成Doris办理所需的影像信息辅影像成像选择主影像兴趣地区从头成像 Iyao观察主影像成像图像,确立研究地区,依据研究地区界限扩展 1000像元的范围来编写roi.proc 。从头运转 step_slc_ers ,依据兴趣地区从头成像。 step_slc_ers 命令用法同1.2./ 粗裁。直接在 SLC文件夹下生成有这个文件
4、,在里面改正,别忘掉去掉 #号。成立主影像精裁,察看主影像成像图像,编写 文件,再次确立裁剪地区。运转step_master_setup 命令成立主影像。 /从ROI_PAC_SCR中复制 到主影像成像的文件夹中。还要注意精裁的裁剪范围是在粗裁后的影像中得范围不是在初始成像后的主影像的队列。Step_master_set 所需参数文件 up生成结果辅影像列表:在与SLC同目录下成立文件夹INSAR_master_date,将、等文件拷入此文件夹作用确立裁剪地区,成立主影像,确立辅影像列表在SLC目录下,运转make_slcs_ers命令,对所有辅影像依据兴趣范围成像。运转过程为逐一进行辅影像目
5、录,履行step_slc_ers命令成像。step_slc_ers用法同1.2.注意其余辅影像成像的时候是依据粗裁的范围裁剪的,他不进行精裁,因为主影像又进行了精裁,一般辅影像比较大。对应1级产品不用成像,直接模仿说明书,读取即可成立连结并读取主影像link_slcs_路径(初始数据的)cdmaster_date(即是主影像日期文件夹)/进入主影像日期文件夹stepreadwholeXXX(whereXXXisERS,Envisat,RSAT,orTSX)4第一次对主影像读取,所有读取,不裁剪5裁剪主影像6cp$MY 或许从安装文件里考7并编写8裁剪后主影像再次读取9step_master_r
6、ead10读取各个辅影像返回到SLC目录make_read2差分干预办理( 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_I
7、NTERVAL1/星历时间间隔(s)M_ORB_EXTRATIME6/Timeinsecondsbeforefirstandlastlinetooutputephemerides.输出结果提出的精细轨道信息保留在中。作用提取主影像精细轨道信息step_orbit履行命令卡,运转Doris中S_PORBITS步骤:S_PORBITS主要参数S_ORBDIR$SAR_ODR_DIR/ERS1/精细轨道文件路径S_ORB_INTERVAL1/星历时间间隔(s)S_ORB_EXTRATIME6/Timeinsecondsbeforefirstandlastlinetooutputephemerides
8、.输出结果提出的精细轨道信息保留在中。作用提取辅影像精细轨道信息主辅影像配准 主辅影像粗配准make_coarse:进入每幅辅影像相应的文件夹,在每个文件夹内履行 step_coarse命令,)进行主辅影像一级轨道配准和二级像素级配准。Step_coarse命令履行 Doris软件的COARSEORB和COARSECORR两步骤。COARSEORB主要参数输出结果作用不需要参数,依赖读取主辅影像的的和来达成轨道配准结果保留在文件中鉴于主辅影像轨道数据配准,精度达 30个像元COARSECORR主要参数CC_METHODmagfft/有关系数计算方法,magfft是先使用fft变换到频次域,而后
9、计算有关系数。magspace是直接在空间域计算有关系数CC_NWIN21/像素级配准区分窗口的数量,大于5.CC_WINSIZE1024256/般配窗口的大小,办理的时候被改变为奇数。CC_INITOFForbit/利用轨道计算的偏移量作为初始偏移量输出结果计算出Coarse_correlation_translation_lines与coarse_correlation_translation_pixels,保留在 文件中。作用 利用有关系数法实现像素级配准。 主辅影像精配准make_coreg:精配准,精配准过程中所有基线距小于设定的bperp_max值的直接与主影像配准,基线距比较大的
10、和离它近来的三个辅影像配准。而后利用用中间过渡的方法实现所有主辅影像的配准。Hooper提出的配准方法,利make_coreg命令运转Doris软件中的FINE和COREGPM步骤,即三级子像素精配准和预计配准多项式系数。FINE主要参数FC_METHODoversample/精配准计算方法,过采样办理FC_IN_POS/ASCII文件,记录主影像上待配准点点的坐标FC_WINSIZE6464/般配窗口的大小FC_ACC88/在方向向和距离向的搜寻偏移量FC_INITOFFcoarsecorr/使用像素级配准后的偏移量FC_OSFACTOR32/过采样的内插倍数输出结果配准结果保留在INSAR
11、_date文件夹下coreg文件夹中的和文件里作用计算主辅影像每个配准点的偏移量COREGPM主要参数CPM_THRESHOLD/经配准后生成大批的般配点,当有关系数大于CPM_THRESHOLD时,才参加平差运算。阈值的大小与File中窗口的大小有关,小窗口状况下阈值趋于,阈值越大越好。CPM_DEGREE2/多项式的阶数,2阶即可。CPM_WEIGHTbamler/最小二乘中权矩阵方法选择,bamler:权矩阵为有关系数的2倍,linear:权矩阵为有关系数的1倍。CPM_MAXITER8000/最小二乘运算的最大迭代次数输出结果计算出方向向和距离向配准多项式系数,输出的结果保留在core
12、g文件夹的文件里。作用计算配准多项式系数主影像强度图模拟与地形相位计算 (平展地域地形可不履行 )cdINSAR_master_date/在INSAR_master_date文件夹,改正文件参数step_master_timing:模拟主影像强度图,计算DEM与主影像的TimingError。make_dems:进入到每个辅影像目录,履行 step_dem命令,计算每个干预对的地形相位。step_master_timing运转Doris软件的M_SIMAMP和M_TIMING 步骤:M_SIMAMP主要参数 SAM_IN_FORMAT I2/r型的 /DEM格式,主要包含 I2(整),R4fl
13、oat)和R8(double)SAM_IN_DEM /sungt/suzhou_test/DEM/suzhou.dem /DEM路径SAM_IN_SIZE 36013601 /DEM行数和列数SAM_IN_DELTA /DEM采样间隔SAM_IN_UL 115.9958333 /DEM左上角纬度和精度,能够在hdr文件里看到,用写字板翻开可看输出结果作用M_TIMING主要参数SAM_IN_NODATA -32768/-9999 /无值区的取值SAM_OUT_FILE /模拟主影像幅度图SAM_OUT_DEM /裁剪的DEM模拟的主影像幅度图 和依据主影像裁剪的 DEMdem_sam_raw或
14、许(队列号在的NcolsDEM和NrowsDEM)模拟主影像幅度图MTE_METHOD magfft /Magfft是先使用fft变换到频率域,而后计算有关系数。Magspace是直接在空间域计算有关系数。Magfft计算快于magspace。MTE_NWIN 30 /区分窗口的数量MTE_INITOFF 00 /初始偏移量MTE_WINSIZE 40962048 /般配窗口大小输出结果 模拟幅度图与主影像方向向和距离向的 TimingError:Master_azimuth_timing_error 和Master_range_timing_error模拟图与主影像的粗配准偏移量 Coars
15、e_correlation_translation_linesCoarse_correlation_translation_pixels;作用计算模拟幅度图与主影像的TimingError以及粗配准偏移量。/DEM偏移量计算。step_dem运转Doris软件中的COMPREFDEM步骤:COMPREFDEM主要参数CRD_METHODtrilinear/求取地形相位方法,有nearest和trilinear两种。因为是粗DEM,与影像的分辨率不一致,需要过采样到同一分辨率下。CRD_INCLUDE_FEOFF/phasew.r.t.ellipsoidCRD_OUT_FILE/由DEM反算获取
16、的地形相位文件CRD_OUT_DEM_LP/原DEM获取的目前雷达坐标系统下的新DEM(已采样)输出结果 地形相位文件和雷达坐标系下的 DEM :RESAMPLE主要参数RS_METHODrc12p/重采样的内插核RS_OUT_FILE/重采样后生成的文件名RS_OUT_FORMATcr4/输出文件数据种类输出结果重采样后的辅影像作用辅影像重采样差分干预办理make_ifgs:主辅影像进行干预办理,去除平川相位和地形相位。进入到每个辅影像目录,运转step_ifg命令。()(make_ifgs_nodem关于平展地域不用去地形)step_ifg运转Doris软件中的SUBTRREFDEM、CO
17、MPREFPHA和SUBTRREFPHA步骤,即 生成干预图,利用轨道来计算平川相位,减去平川相位,减去由粗 DEM反算获取的地形相位:inter主要参数INT_OUT_CINT/生成的干预图文件,为复数float型INT_MULTILOOK11/多视办理系数输出结果干预图:作用干预办理,生成干预图COMPREFPHA主要参数FE_METHODporbits/采纳计算平川相位的方法,精轨FE_DEGREE3/地点与平川相位之间的多项式系数FE_NPOINTS201/最小二乘法来计算平川相位所用来拟合的点数输出结果平川相位终端履行过程存于中作用依据精轨信息计算平川相位SUBTRREFPHA(SU
18、BTRREFDEM)主要参数SRP_METHODpolynomial/采纳减去平川相位的方法SRP_OUT_CINT/去除平川效应后的干预图输出结果去除平川相位的干预图:(去地形后没有了)作用去除平川相位影响SUBTRREFDEM主要参数SRD_OUT_CINT去除地形相位后的复干预图文件SRD_OFFSET00/鉴于粗DEM获取的地形相位与原始干预合成相位之间的间隔差别。输出结果去除地形相位后的干预图文件作用去除地形相位地理编码step_geo:从雷达坐标系变换到地图坐标系/只在INSAR_masterdate文件夹下随意一个辅文件(就是日期那个文件夹)夹下运转一次就行。Step_geo运转
19、Doris软件中的GEOCODE步骤:GEOCODE主要参数GEO_OUT_LAM/经度输出文件GEO_OUT_PHI/纬度输出文件输出结果像元经纬度信息和(就生成一个结果在INSAR文件夹,就是积累的形变。)作用地理编码,参照坐标系变换3PS办理准备注意小基线不是在此步初选点,是在确立小基线对后再初选点,别弄错了/命令是mt_prep 3250200 仍是在INSAR masterdate 目录下mt_prep:依据设定的幅度离差阈值、距离向方向向分块和重叠范围等参数进行办理,选定初始PS候选点。mt_prep在运转过程中还调用了mt_extract_info和mt_extract_cand
20、smt_prep主要参数da_thresh/幅度离差阈值,一般取值,小于此幅度离差的选为PS候选点,默以为rg_patches/距离向上分块,默以为1az_patches/方向向上分块,默以为1rg_overlap/距离向上块间重叠像元数,默认50az_overlap/方向向上块间重叠像元数,默认200输出结果成立PATCH(分几个块,产生几个文件)文件夹,(直接在INSARmasterdate目录下)主影像行数和列数、,基线散布bperp.1.in,辅影像日期day.1.in,主影像日期,干预对日期,文件头信息,主影像雷达波长,视角文件,幅度定标信息,DEM参数分块状况,均匀幅度图,分块范围
21、和,PS候选点和;每个候选点的经纬度信息,每个候选点的高程信息,每个候选点的相位信息等phname=pscands.1.ph;%foreachPScandidate,afloatcomplexvalueforeachifgijname=pscands.1.ij;%ID#Azimuth#Range#1lineperPScandidatebperpname=bperp.1.in;%inmeters1lineperslaveimagedayname=day.1.in;%YYYYMMDD,1lineperslaveimageifgdayname=ifgday.1.in;%YYYYMMDDYYYYMMD
22、D,1lineperifgmasterdayname=master_day.1.in;%YYYYMMDDllname=pscands.1.ll;%2floatvalues(lonandlat)perPScandidatedaname=pscands.1.da;%1floatvalueperPScandidatehgtname=pscands.1.hgt;%1floatvalueperPScandidatelaname=look_angle.1.in;%gridoflookanglevaluesheadingname=heading.1.in;%satelliteheadinglambdanam
23、e=lambda.1.in;%wavelengthcalname=calamp.out;%amplitidecalibrationswidthname=width.txt;%widthofinterferogramslenname=len.txt;%lengthofinterferograms作用提取所需信息,幅度校订,成立分块,选用初始PS候选点stamps办理/在INSARmasterdate下先启动matlabgetparm(parmname),无参数时显示所有参数,默认参数从提取;setparm(parmname,value),改正某个参数值。参数结果保留在文件里。/在INSARmasterdate下,运转stampsstamps(start_step,end_step,patches_flag,est_gamma_parm)Patches_flag默以为y,假如设置为n把所有数据作为一块办理。Est_gamma_parm是一个可选参数passedtops_est_gamma_quick。选y则利用ps_est_gamma_quick(est_gamma_parm),n利用ps_est_gamma(est
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1