StaMPS软件操作流程

StaMPS软件操作流程注意：冒号后面都是解释，如果步骤。没说更改目录就一直是前一步的目录。每个命令都有一个log对应1成像处理（R0I_PAC）SAR影像原始数据准备建立SLC文件夹，将SAR原数据拷贝到SLC目录里，为每景影像按照获取时间分别建立yyyymmdd文件夹，每景SAR影像原始数据对应的数据和头文件分别命名为IMAGERYyymmdd和SARLEADERyyyymmdd。主影像成像（对应的是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命令介绍

slc_ersraw.p1用orbit_typeleader_filedate输入参数轨道类曲(ODR)、数据头文件(SARLEADERyyyymmdd)和影像获取时间yyyymmdd输出主要结果<date>.raw 、<date>.raw.rsc 、<date>_parse_lines.out、shift,out、shift,out.rsc作用SAR原始数据预处理，主要完成1)解析数据产品文件，将其转为R0I_PAC的内部格式；2)填补原始中缺失的行；3)调整采样窗口的起始时间；4)提取成像参数Roi_prep.p1命令使用roi_prep.pl dateOrbitType输入参数影像获取时间(yyyymmdd)>轨道类型(ODR)输出主<date>.sic.rsc、

要结果<date>.roi.in 、<date>.roi.out作用生成sic参数文件Roi.p1命令使用roi.pldateno1no输入参数影像获取时间(yyyymmdd)输出主要结果<date>.sic作用生成sic图像roipac2doris命令使用Roipac2dorisdate[s]输入参数影像获取时间(yyyymmdd)输出主要结果master,res 或slave,res作用生成Doris处理所需的影像信息1.3选择主影像兴趣区域重新成像Iyao观测主影像成像图像，确定研究区域，按照研

究区域边界扩展1000像元的范围来编辑

co重新运行step_slc_ers,按照兴趣

区域重新成像。step_slc_ers命令用法同

L2.〃粗裁。直接在SLC文件夹下生成有这个

文件，在里面修改，别忘记去掉#号。1.4建立主影像精裁，观察主影像成像图像，编辑master_crop.in文件，再次确定裁剪区域。运行step_master_setup命令建立主影像。//从ROI_PAC_SCR中复制master_crop.in到主影像成像的文件夹中。Step_master_setup所需参数文件生成结果master_crop.in还要注意精裁的裁剪范围是在粗裁后的影像中得范围不是在初始成像后的主影像的行列。Step_master_setup所需参数文件生成结果master_crop.in辅影像列表：Make_slcs.list在与SLC同目录下建立文件夹INSAR_master_date,将masterdate.slc.rsc、master_crop.slcdem.dorisingeocode.dorisntiming.dorisinmaster.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）提取精密轨道信息（冒号后面的是此步骤的解释说明，不是让你这么操作）在insarmaster文件夹下step_master_orbit_ODR：主影像精密轨道信息提取。〃没有精轨不运行这步，直接下步make_orbits：为每幅辅影像建立相应的文件夹，在每个文件夹内执行step_orbit命令，读取相应辅影像精密轨道信息。step_master_orbit_ODR运行Doris中的M_PORBITS步骤：M_PORBI TS主要参数M_ORBDIR $SAR_ODR_DIR/ERS1〃精密轨道文件路径//M_ORB_INTERVAL1//星历时间间隔（s）MORBEXTRATIME//Timeinsecondsbeforefirstandlastlinetooutputephemerides.输出结果提出的精密轨道信息保存在

master.res中。作用提取主影像精密轨道信息作用step_orbit执行orbit_ers.dorisin命令卡,运彳亍Doris中S_PORBITS步骤：S_PORBITS主要参数输出结果作用S_PORBITS主要参数输出结果作用$SAR_ODR_DIR/ERS1〃精密轨道文件路径S_ORB_INTERVAL1星历时间间隔（S）SORBEXTRATIME6////Timeinsecondsbeforefirstandlastline////Timetooutputephemerides.提出的精密轨道信息保存在slave.res中。提取辅影像精密轨道信息 主辅影像配准主辅影像粗配准make_coarse：进入每幅辅影像相应的文件夹，在每个文件夹内执行step_coarse命令，(coarse.dorisin)进行主辅影像一级轨道配准和二级像素级配

准。Step_coarse命令执行Doris软件的COARSEORB和COARSECORR两步骤。COARSEORB主要参数不需要参数，休…取士锚皇’傍的的f和依靠读取王辅影像的的master.res和slave.res来兀成输出结果轨道配准结果保存在coreg.out文件中作用基于主辅影像轨道数据配准，精度达30个像兀COARSECORR主要参数CC_METHOD magfft〃相关系数计算方法，magfft是先使用fft变换到频率域，然后计算相关系数。magspace是直接在空间域计算相关系数CC_NWIN 21 〃像素级配准划分窗口的数目，大于5.CC_WINSIZE 1024256 〃匹配窗口的大小，处理的时候被改变为奇数。CC_INITOFF orbit 〃利用轨道计算的偏移量作为初始偏移量， “ 一输出结果计 算 出Coarse_correlation_translation_lines与coarse_correlation_translation_pixels，保存在coreg.out文件中。作用利用相关系数法实现像素级配准。主辅影像精配准make_coreg：精配准，精配准过程中所有基线距小于设定的bperp_max值的直接与主影像配准，基线距比较大的和离它最近的三个辅影像配准。然后利用Hooper提出的配准方法，利用

中间过渡的方法实现所有主辅影像的配准。make_coreg命令运行Doris软件中的FINE和COREGPM步骤，即三级子像素精配准和估计配准多项式系数。FINE主要参数FC_METHOD oversample 〃精配准计算方法，过采样处理FC_IN_POS fc_pos.in //ASCII文件，记录主影像上待配准点点的坐标FC_WINSIZE 6464 //匹配窗口的大小FC_ACC 88 〃在方位向和距离向的搜索偏移量FC_INITOFF coarsecorr 〃使用像素级配准后的偏移量FC_OSFACTOR 32 //过采样的内插倍数输出结果配准结果保存在INSAR_date文件夹下coreg文件夹中的coreg.out和CPM_DATA.n1.n2文件里作用计算主辅影像每个配准点的偏移量COREGPM主要参数CPM_THRESHOLD0.30 〃经配准后生成大量的匹配点，当相关系数大于CPM_THRESHOLD时，才参与平差运算。阈值的大小与File中窗口的大小有关，小窗口情况下阈值趋于1.0，阈值越大越好。CPM_DEGREE 2 〃多项式的阶数，2阶即可。CPM_WEIGHT bamler //最小二乘中权矩阵方法选择，bamler：权矩阵为相关系数的2倍，linear：权矩阵为相关系数的1倍。CPM_MAXITER(、80叩，। (/俚小二乘运算的最大迭代举二输出结果计算出方位向和距离向配准多项式系数，输出的结果保存在coreg文件夹的coreg.out文件里。作用计算配准多项式系数主影像强度图模拟与地形相位计算(平坦地区地形可不执行)cdINSAR_master_date//在INSAR_master_date文件夹修改timing.dorisin文件参数step_master_timing：模拟主影像强度图，计算DEM与主影像的TimingError。make_dems：进入到每个辅影像目录，执行step_dei命令，计算每个干涉对的地形相位。step_master_timing运行Doris软件的M_SIMAMP和M_TIMING步骤：MP主要参数SAM_IN_FORMATI2/r型的主要参数在后乂格式，主要包括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输出结果mastersam.raw模拟的主影像幅度图主影像裁剪的DEMdem_和根据samraw或者demcrop.raw（行列号在作用 M_TIMING主要参数step_master_timing.logNcolsDEM和NrowsDEM模拟主影像幅度图MTEMETHODmagfft//Magfft是先使用fft变换到频率域，然后计算相关系数。Magspace是直接在空间域计算相关系数。Magfft计算快于magspace。MTE_NWIN 30//划分窗口的数目MTE_INITOFF 00〃初始偏移量MTE_WINSIZE40962048//匹配窗口大小输出结果模拟幅度图与主影像方位向和距离向的TimingError：Master_azimuth_timing_error和Master_range_timing_error模拟图与主影像的粗配准偏移量Coarse_correlation_translation_lines 与Coarse_correlation_translation_pixels；作用计算模拟幅度图与主影像的TimingError以及粗配准偏移量。/DEM偏移量计算。steP_dem运行加示软件中的comprefdem步骤：COMPREFDEM主要参数CRD_METHOD trilinear 〃求取地形相位方法，有neares4Dtrilinear两种。由于是粗DEM，与影像的分辨率不一致，需要过采样到同一分辨率下。CRD_INCLUDE_FEOFF //phasew.r.t.ellipsoidCRD_OUT_FILErefdem_1l.raw〃由DEM反算得到的地形相位文件

CRD_OUT_DEM_LPdem_radar.raw〃原DEM得到的当前雷达坐标系统下的新罂后采样）、./年— “ 如H4输出结果地形相位文件refdem_1l.raw和雷达坐标系下的DEM：dem_radar.raw（在各辅影像文件夹中），dem_radar_i.raw（在主影像文件夹中）这几个文件和主影像行列号一样，float格式。作用计算地形相位2.4辅影像重采样make_resample：对辅影像按照精配准算出的多项式来重采样。进入到每个辅影像目录，执行step_resample。（resample.dorisin,自动力口入RS_DBOW10001099915003019）step_resample运行Dgs软件中的RESAMPLE步骤：RESAMPLE主要参数RS_METHOD rc12P 〃重采样的内插核RS_OUT_FILE slave_res.slc 〃重采样后生成的文件名RS_OUT_FQRMAT__cr4，“一，，一〃输出文件数据类型输出结果重米样后的辅影像slave_res.slc作用辅影像重采样2.5差分干涉处理make_ifgs：主辅影像进行干涉处理，去除平地相位和地形相位。进入到每个辅影像目录，运

行step_ifg命令。（integram.dorisin）（make_ifgs_nodem对于平坦地区不用去地形）stepifg运行Doris软件中的■—d SUBTRREFDEM、COMPREFPHA和SUBTRREFPHA步骤，即 生成干涉图，利用轨道来计算平地相位，减去平地相位，减去由粗DEM反算得到的地形相位：inter主要参数INT_OUT_INTMULTCINTcint.raw 〃生成的干涉图文件，为复数float型ILOOK11 //多视处理系数输出结果干涉冬：cint.raw作用干涉处理，生成干涉图COMPREFPHA主要参数输出结果FE_METHOFE_DEGREFE_NPOIN平地InterDporbits 〃采用计算平地相位的方法，精轨E 3 〃位置与平地相位之间的多项式系数TS相2位/终小端法执计算学过所程来拟存的于ferogram.out中数作用根据精轨信息计算平地相位SUBTRREFPHA(SUBTRREFDEM).,，3巳要参数SRP_METHOD polynomial 〃采用减去平地相位的方法SRP_OUT_CINTcint.minrefpha.raw 〃去除平地效应后的干涉图输出结果去除平地相位的干涉图：cint.minrefpha.raw（去地形后没有了）作用去除平地相位影响SUBTRREFDEM主要参数SRD_OUT_CINTcint.minrefdem.raw//去除地形相位后的复干涉图文件SRD_OFFSET 00 〃基于粗DEM得到的地形相位与原始干涉合成相位之间的间隔差异。输出结果去除地形相位后的干涉图文件cint.minrefdem.raw作用去除地形相位2.6地理编码step_geo：从雷达坐标系转换到地图坐标系//只

在INSAR_masterdate文件夹下任意一个辅文件（就是日期那个文件夹）夹下运行一次就行。SteP_ge0运行Doris软件中的GEO.步骤:GEOCODE主要参数GEO_OUT_LAMGEO_OUT_PHIlon.raw 〃经度输出文件lat.raw 〃纬度输出文件输出结果像元经就是累纬度信息lon.raw和lat.raw（就生成个结果在INSAR文件夹，积的形变。）八一-一一一，」…作用地出码，参考坐标系转换3PS处理3.1准备注意小基线不是在此步初选点，是在确定小基线对后再初选点，别弄错了〃命令是mt_prep0.43250200〃命令是mt_prep0.43250200还是在INSARmasterdate目录下mt_prep：按照设定的幅度离差阈值、距离向方位向分块和重叠范围等参数进行处理，选定初始PS候选点。mt_prep在运行过程中还调用了 壬mt_extract_info和mtextractcandsmt_prep主要参数da_thresh 〃幅度离差阈值，一般取值0.4-0.42，小于此幅度离差的选为PS候选点，默认为0.4rg_patches 〃距离向上分块，默认为1az_patches 〃方位向上分块，默认为1rg_overlap 〃距离向上块间重叠像元数，默认50az_overlap 〃方位向上块间重叠像元数，默认200输出结果建立PATCH（分几个块，产生几个文件）文件夹，（直接在INSARmasterdate目录下）主影像行数和列数

width.txt、len.txt，基线分布bperp.1.in,辅影像日期day.1.in,主影像日期master_day.1.in，干涉对日期ifgday.1.in，文件头信息heading.1.in，主影像雷达波长lambda.1.in，视角文件look_angele.1.in，幅度定标信息calamp.out，DEM参数demparm.in分块情况patch.list，平均幅度图mean_amp.flt，分块范围patch.in和patch_noover.in，PS候选点pscands.1.ij和pscands.Lda；每个候选点的经纬度信息pscands.Lll，每个候选点的高程信息pscands.1.hgt，每个候选点的相位信息pscands.1.ph等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']; %YYYYMMDDYYYYMMDD,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']; % satelliteheadinglambdaname=['lambda.1.in']; % wavelengthcalname=['calamp.out']; % amplitidecalibrationswidthname=['width.txt']; % widthofinterferogramslenname=['len.txt'>];_ _^、%length.ofinterferograms作用提取所需信息，幅度校正，建立分块，选取初始PS候选点3.2stamp处理〃在INSARmasterdate下先启动matla)getparm(parmname),无参数时显示所有参数，默认参数从ps_parms_default.m提取默认参数从ps_parms_default.m提取setparm('parmname’,value)，修改某个参数值。参数结果保存在parms.mat文件里。〃在INSARmasterdate下，运行stamps〃在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_gamma_parm)。stamps(1,1):调用ps_load_initial.m，将PS处理所需各种文件加载进matlab工作空间ps_load_inital所需文件pscands.1.ph，pscands.1.ij，bperp.1.in， day.1.in，ifgday.1.in，master_day.1.in，pscands.1.ll，pscands.1.da，

pscands.l.hgt ,look_angle.1.in ,heading.l.in,lambda.I.in,calamp.out,len.txt,width.txt输出结果干涉图中每个候选点的基线分布bpl.mat,每个候选点的幅度离差dal.mat,每个候选点的高度值hg1.mat，每个候选点的视角la1.mat,每个候选点的干涉相位ph1.mat，每个候选点的点位、行列号、经纬度等信息psLmat作用 将PS处理所需要的各种文件加载到Matlab工作空间，生成相应的.mat文件stamps(2,2)：调用ps_est_gamma_quick.m，迭代计算每个干涉图中每个候选点的噪声相位，计算每个候选点的y值。ps_est_gamma_quick主要参数max_topo_err5主要参数max_topo_err5/MaximumuncorrelatedDEMerror(inm).PixelswithuncorrelatedDEMerrorgreaterthanthiswillnotbepicked。filter_grid_size50//Pixelsize

filter_grid_size50//Pixelsizeofgrid(inm).Candidatepixelsareresampledtoagridwiththisspacingbeforefilteringtodeterminethespatiallycorrelatedphase.filter_weighting 'P-square'//Weightingscheme(PSprobabilitysquared),theotherpossibilitybeing‘SNR’.Candidatepixelsareweightedduringresamplingaccordingtothisscheme.clap_win 32 //CLAP(CombinedLow-passandAdaptivePhase)filterwindowsize。窗口大小单位不是米，是32*32的矩阵clap_low_pass_wavelength800//CLAPfilterlow-passcontributioncut-offspatialwavelength(inm).Wavelengthslongerthanthisarepassed.反算出频率，因为是低通频域滤波。clap_alpha1 //CLAPaterm.Togetherwiththepterm,determinestherelativecontributionofthelow-passandadaptivephaseelementstotheCLAPfilter.clap_beta0.3 //CLAPptermgamma_change_convergence0.005//ThresholdforchangeingammachangeinmeanvalueofY.Determineswhenconvergenceisreachedanditerationceases.收敛阈值Low_coh_thresh 31//就是0.31，较低y的阈值，小于次值时没有ps点，应用于在P-square方法算weighing时。

输出结果pm1.mat作用迭代计算每个干涉图中每个候选点的噪声相位，计算每个候选点的y值stamps(3,3)：调用ps_select.m，基于计算出的每个候选点的y值，同时考虑幅度离差选取PS像元。ps_select主要参数selectmethod\DENSITY, //Otheroption\PERCENT’.densityrand 20 //Maximumacceptablespatialdensity(perkm2)ofselectedpixelswithrandomphase.点太多的话这个值可以小些如2percentrand 20 //Maximumacceptablepercentageofselectedpixelswithrandomphase.输出结果select1.mat作用基于计算出的每个候选点的y值，同时考虑幅度离差选取PS像元ps_weed主要参数stamps(4,4):调用ps_weed.m,去除那些只在一些干涉图中相位稳定的PS点及其受临近PS点影响而表示为PS点特征的像元。同时根据噪声标准差ps_weed主要参数weed_alpha8 //Smoothingparameterforestimatingphasenoisedistributionforeachpairofneighbouringpixels.ThetimeseriesphaseforeachpairissmoothedusingaGaussianwindowwithstandarddeviation1/weedalpha.Theoriginalphaseminusthesmoothedphaseisassumedtobenoise.weed_standard_dev1.0//Thresholdstandarddeviation.Iftheminimumstandarddeviationisgreaterthanthethreshold,thepixelis

dropped.输出结果psweed.1.node,psweed2.node,weed1.mat,pm2.mat,ph2.mat,ps2.mat,hgt2.mat,la2.mat,bp2.mat作用去除那些只在些干涉图中相位稳定的PS点及其受临近PS点影响而表示为PS点特征的像元stamps(5,5):调用ps_correct_phase,m和ps_calc_ifg_std.m(计算每一个干涉图噪声标准差，作为权重)，对选取的PS像元的缠绕相位去除spatia^-uncorrelatedlookangel(DEM)error.还要调用error.还要调用ps_merge_patchesps_correct_phaseps_calc_if主要参数

输出结果主要参数

输出结果merge_resample_size0 //Coarserposting（inm）toresampleto.Ifsetto0,noresamplingisapplied.rC2Tm0t（去除SULA和masternoise后相位），ifgstd2.mat（噪声标准差）作用 对选取的PS像元的缠绕相位纠spatiall-uncorrelatedlookangel(DEM)error.stamps(6,6):调用ps_unwrap.m，对选取的PS像元进行相位解缠。ps_unwrap.m主要参数Unwrap_method‘3D' //Unwrappingmethod.Unwrap_ifg_index‘all'//Indextointerferogramstobeunwrapped.Unwrap_prefilter_flag‘y'//Prefilterphasebeforeunwrappingtoreducenoise.Unwrap_patch_phase‘n’ //UsethepatchphasefromStep3asprefilteredphase.Ifsetto‘n’(recommended),PSphaseisfilteredusingaGoldsteinadaptivephasefilter.Unwrap_grid_size100 //Resamplinggridspacing.Unwrap_gold_n_win32 //WindowsizeforGoldsteinfilter.Unwrap_time_win180 //Smoothingwindow(indays)forestimatingphasenoisedistributionforeachpairofneighbouringpixels.Unwrapgoldalpha0.8//Valueof°forGoldsteinfilter. _输出结果phuw2.mat, unwrap.1.node,unwrap.2.node作用对选取的PS像元的缠绕相位进行相位解缠stamps(7,7):调用ps_calc_scla.m、ps_smooth_scla.m(减去SCLA后重新解缠时用)，估计spatially-correlatedlookangleerror，同时估计出主影像大气影响相位和轨道误差相位(AOE)

stamps（8,8）:ps_scn_filt.m,滤除其他空间相关项如辅影像AOE。ps_calc_scla.mps_smooth_scla.mps_scn_filt.m主要参数Recalc_index‘all'//IndextointerferogramstobeusedintheSCLAestimation.Scla_deramp‘y' //Ifsetto‘y',aphaserampisestimatedforeachMterferogrOm.0ther0Pti0nisIn；-. 工-「 、目，小[一输出结果scla2.mat（SCLA和C主影像大气和轨道误差也在里面）,scla_smooth2.mat（滤波后的…）,scn2.mat（滤波后的各项误差相位）作用估计spatially-correlatedlookangleerror，同时估计出主影像大气影响相位和轨道误差相位（AOE）最后输出显示ps_plot（'v'，4）4小基线处理如果没有运行前面的PS处理，则需要先在INSARmasterdate目录下读取基线信息加载到matlab工作空间，命令如下

mt_extract_infomatlab>>ps_load_info确定小基线对sb_find主要参数RHO_MIN//最小相关系数（默认0.50）可以减小这个参数获得更多小基线对设置成0.1或0.2DDIFF_MAX〃完全失相关的时间间隔（默认1500天）BDIFF_MAX〃临界基线值（默认1070）输出结果small_baselines.list作用以相天系数为条件寻找小基线干涉对〃终端在INSARmasterdate目录下运行在matlab里sb_find:寻找小基线干涉对//sb_find（rho_min,ddiff_max,bdiff_max）运行后显示基线图。也可以使用命令调出基线图plotsbbaselinesmake_small_baselines：方位向、距离向滤波，按照选择的小基线对，生成差分干涉图。运行Doris软件中的

m_filtazi，s_filtazi，FILTRANGE与cpxfiddle命令。INTERFERO步骤，及其Doris的cpxsumcpxfiddle命令。并在INSARmasterdate下^^成SMALLBASELINES文件夹将原来配准后的文件slave_res.slc考入到相应的干涉对组成的文件夹内分别进行主辅影像方位向滤波与距离向滤波，生成干涉图。利用PS处理中每个干涉对中的干涉图和去除地形相位的差分干涉图，结合现有的干涉图，用Doris命令cpxsum处理，生成新的差分干涉图。目前，对此过程的执行过程还不是很清楚。m_filtazi，s_filtazi主要参数AF_BLOCKSIZE4096 //fftlengtheachcolumnFFT变换中每次距离向中的行数，一般是越大越好。AF_OVERLAP 64 //Halfoftheoverlapbetweenconsecutivebufferesinazimuthdirection.Partiallythesamedataisusedtoestimatethespectrum,whichmighthavecertainadvantages.cAF_HAMMING0.75 〃光谱滤波的权值。滤波器的alpha系数。AF_OUT_MASTER master_azifilt.slc 〃滤波后的主影像文件名AF_OUT_SLAVE slave_azifilt.slc 〃滤波后的辅影像文件名AF_OUT_FORMAT cr4 〃输出文件数据的格式输出结果方位向滤波后的主影像master_azifilt.slc方位向滤波后的辅影像slave_azifilt.slc作用主辅影像方位向滤波FILTRANGE主要参数RF_METHOD adaptive 〃距离向滤波的方法RF_FFTLENGTH256 //LengthofwindowforadaptivemethodRF_NLMEAN 15 //oddTakemeanoverRF_NLMEANlinestoreducenoiseforpeakestimation.Hastobeodd.

RF_THRESHOLD 5 //ThresholdonSNRofpeakestimationtoperformrangefilteringRF_HAMMING 0.75 //alphaRF_OVERSAMPLE 4 //Oversamplemasterandslavewiththisfactorbeforecomputingthecomplexinterferogramforpeakestimation.RF_WEIGHTCORROFF//Inpeakestimation,weightvaluestobiashigherfrequencies.RF_OUT_MASTERmaster_filtrg.slc〃距离向滤波后的主影像文件RF_OUT_SLAVEslave_filtrg.slc〃距离向滤波后的辅影像文件RF_OUT_FORMATcr4〃输出文件数据的格式输出结果距离向滤波后的主影像文件master_filtrg.slc距离向滤波后的辅影像文件slave_filtrg.slc作用主辅影像距离向滤波INTERFERO主要参数INT_OUT_CINT cint.raw //干涉图INT_MULTILOOK 11 //多视输出结果（干涉图cint.raw）去除平地和地形的差分干涉图cint.minrefdem作用主辅影像生成干涉图mt_prepmt_prep：按照设定的幅度差分离差阈值、距离向方位向分块和重叠范围等参数进行处理，选定初始高相干候选点。//cdSMALLBASELINES//cdSMALLBASELINES进入这个文件夹mt_prep0.63250200mt_prep在运行过程中还调用了mtextractinfo、sb_parms_initial.m和mt_extract_candsmt_prep主要参数da_thresh 〃幅度差分离差阈值，一般取值0.6，小于此幅度差分离差的选为高相干候选点，默认为0.6rg_patches 〃距离向上分块，默认为1az_patches 〃方位向上分块，默认为1rg_overlap 〃距离向上块间重叠像元数，默认50az_overlap 〃方位向上块间重叠像元数，默认200

输出结果建立PATCH文件夹，主影像行数和列数width.txt、len.txt，基线分布bperp.Lin,辅影像日期day.Lin,主影像日期master_day.1.in，干涉对日期ifgday.Lin，文件头信息heading.Lin，主影像雷达波长lambda.1.in，视角文件look_angele.1.in，幅度定标信息calamp.out，平均幅度图mean_amp.flt，DEM参数demparm.in，分块情况patch.list，分块范围patch.in和patch_noover.in，PS候选点pacands.Lij和pscands.1.da；每个候选点的经纬度信息pscands.Lll，每个候选点的高程信息pscands.1.hgt，每个候选点的相位信息、pscands..1,ph等_ r_-r、 .• 一r.、 ，、作用提取所需信息，幅度校正，块，选取初始高相干候选点建立分S处理运行过程同PS中stamps的使用，可参考上面。先启动matlab先启动matlab还是在SMALLBASELINES下getparm(parmname),无参数时显示所有参数，默认参数从ps_parms_default.m默认参数从ps_parms_default.m提取setparm('parmname’,value)，修改某个参数值。参数结果保存在parms.mat文件里。Stamps(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_gamma_parm)。stamps(1,1)：调用sb_load_initial.m,将小基线处理所需各种文件加载进matlab工作空间sb_load_inital所需文件 pscands.1.ph，，pscands.1.ij，bperp.1.in，day.1.in，ifgday.1.in，master_day.1.in，pscands.1.ll，pscands.1.da，pscands.1.hgt ，look_angle.1.in ，heading.1.in，lambda.1.in，calamp.out，len.txt,width.txt输出结果干涉图中每个候选点的基线分布bpLmat，每个候选点的幅度离差daLmat，每个候选点的高度值hg1.mat，每个候选点的视角la1.mat,每个候选点的相位ph1.mat，每个候选点的点位、行列号、经纬度等信息psLmat作用 将小基线处理所需要的各种文件加载到Matlab工作空间，生成相应的.mat文件

stamps(2,2)：调用ps_est_gamma_quick.m,迭代计算每个干涉图中每个候选点的噪声相位，计算每个候选点的y值。ps_est_gamma_quick u l//N育 ■ 主要参数max_topo_err 5//MaximumuncorrelatedDEMerror(inm).PixelswithuncorrelatedDEMerrorgreaterthanthiswillnotbepicked。filter_grid_size 50//Pixelsizeofgrid(inm).Candidatepixelsareresampledtoagridwiththisspacingbeforefilteringtodeterminethespatiallycorrelatedphase.filter_weighting 'P-square'//Weightingscheme(PSprobabilitysquared),theotherpossibilitybeing‘SNR’.Candidatepixelsareweightedduringresamplingaccordingtothisscheme.

clap_win 32 //CLAP(CombinedLow-passandAdaptivePhase)filterwindowsize。clap_low_pass_wavelength800//CLAPfilterlow-passcontributioncut-offspatialwavelength(inm).Wavelengthslongerthanthisarepassedclap_alpha1//CLAPaterm.Togetherwiththepterm,determinestherelativecontributionofthelow-passandadaptivephaseelementstotheCLAPfilter.clap_beta0.3 //CLAPptermgamma_change_convergence0.005 //Thresholdforchangeinchangeinmeanvalueof丫.Determineswhenconvergenceisreachedanditerationceases.Small_bseline_flag'y’

Low_coh_thresh 15输出结果pm1.mat作用迭代计算每个干涉图中每个候选点的噪声相位，计算每个候选点的y值stamps(3,3)：调用ps_select.m，基于计算出的每个候选点的y值，同时考虑幅度离差选取高相干像元。ps_select主要参数selectmethod\DENSITY, //Otheroption\PERCENT’.densityrand 2 //Maximumacceptablespatialdensity(perkm2)ofselectedpixelswithrandomphase.percentrand 1 //Maximumacceptablepercentageofselectedpixelswithrandomphase.输出结果select1.mat作用基于计算出的每个候选点的y值，同时考虑幅度离差选取高相干像元stamps(4,4):调用ps_weed.m,去除那些只在一些干涉图中相位稳定的高相干点及其受临近高相干点影响而表示为高相干点特征的像元。ps_weed主要参数weed_alpha8 //Smoothingparameterforestimatingphasenoisedistributionforeachpairofneighbouringpixels.ThetimeseriesphaseforeachpairissmoothedusingaGaussianwindowwithstandarddeviation1/weedalpha.Theoriginalphaseminusthesmoothedphaseisassumedtobenoise.weed_standard_dev10//Thresholdstandarddeviation.Iftheminimumstandarddeviationisgreaterthanthethreshold,thepixelisdropped.设置成1比较收紧。weedl.mat,pm2.mat,ph2.mat,ps2.mat,hgt2.mat,la2.mat,bp2.mat作用去除那些只在一些干涉图中相位稳定的高相干点及其受临近高相干点影响而表示为高相干点特征的像元作用stamps(5,5):调用ps_c