StaMPS软件操作流程

1、注意：冒号后面都是解释，如果步骤。没说更改目录就一直是前一步的目录。每个命令都有一个log对应1成像处理(ROI_PAC)SAR影像原始数据准备建立SLC文件夹，将SAR原数据拷贝到SLC目录里，为每景影像按照获取时间分别建立yyyymmdd文件夹，每景SAR影像原始数据对应的数据和头文件分别命名为IMAGERYyymmdd和SARLEADERyyyymmdd。主影像成像(对应的是0级)按照主影像选取原则选取主影像，进入到主影像文件夹master_date,对于ERS-1/2数据，运行step_slc_ers命令成像。step_slc_ers主要通过调用ROI_PAC软件中的,命令，及其Sta

2、MPS本身的roipac2doris命令来完成。stepslcers命令介绍Step_slc_ers命令使用orbittypeleaderfiledate输入参数轨道类型(ODR)、数据头文件(SARLEADERyyyymmdd)和影像获取时间yyyymmdd输出主要结果作用SAR原始数据预处理，主要完成1)解析数据产品文件，将其转为ROI_PAC的内部格式；2)填补原始中缺失的行；3)调整采样窗口的起始时间；4)提取成像参数命令使用dateOrbitType输入参数影像获取时间(yyyymmdd)、轨道类型(ODR)输出主要结果<date>.、<date>.、<

3、date>.作用生成slc参数文件命令使用dateno1no输入参数影像获取时间(yyyymmdd)输出主要结果<date>.slc作用生成slc图像roipac2doris命令使用Roipac2dorisdates输入参数影像获取时间(yyyymmdd)输出主要结果或作用生成Doris处理所需的影像信息选择主影像兴趣区域重新成像Iyao观测主影像成像图像，确定研究区域，按照研究区域边界扩展1000像元的范围来编辑。重输出结果提出的精密轨道信息保存在中。作用提取主影像精密轨道信息新运行stepslcers,按照兴趣区域重新成像。stepslcers命令用法同.steporbi

4、t执行命令卡，运行Doris中SPORBITS步骤:S_PORBITS主要参数SORBDIR$SARODRDIR/ERS1输出结果提出的精密轨道信息保存在中。作用提取辅影像精密轨道信息主辅影像配准主辅影像粗配准make_coarse:进入每幅辅影像相应的文件夹，在每个文件夹内执行step_coarse命令，（）进行主辅影像一级轨道配准和二级像素级配准。Stepcoarse命令执行Doris软件的COARSEORB和COARSECORR两步骤。COARSEORB主要参数不需要参数，依靠读取主辅影像的的和来完成输出结果轨道配准结果保存在文件中作用基于主辅影像轨道数据配准，精度达30个像兀COARS

5、ECORR主要参数CCMETHODmagfftCCWINSIZE1024256stamps(2,2)：调用，迭代计算每个干涉图中每个候选点的噪声相位，计算每个候选点的值。psestgammaquick主要参数max_topo_err5PixelswithuncorrelatedDEMerrorgreaterthanthiswillnotbepicked。filter_grid_size50Candidatepixelsareresampledtoagridwiththisspacingbeforefilteringtodeterminethespatiallycorrelatedphase.f

6、ilter_weighting'P-square'Candidatepixelsareweightedduringresamplingaccordingtothisscheme.clap_win32Wavelengthslongerthanthisarepassed反算出频率，因为是低通频域滤波。clap_alpha1aTogetherwiththe卩term,determinestherelativecontributionofthelow-passandadaptivephaseelementstotheCLAPfilter.clap_beta卩YDetermineswhe

7、nconvergenceisreachedanditerationceases.收敛阈值Lowcohthresh31Ydensityrand20点太多的话这个值可以小些如2percentrand20输出结果作用基于计算出的每个候选点的丫值，同时考虑幅度离差选取PS像兀stamps(4,4):调用，去除那些只在一些干涉图中相位稳定的PS点及其受临近PS点影响而表示为PS点特征的像元。同时根据噪声标准差去除噪声点。psweed主要参数weedalpha8ThetimeseriesphaseforeachpairissmoothedusingaGaussianwindowwithstandardde

8、viation1/weedalpha.Theoriginalphaseminusthesmoothedphaseisassumedtobenoise.weedstandarddevIftheminimumstandarddeviationisgreaterthanthethreshold,thepixelisdropped.输出结果作用去除那些只在一些干涉图中相位稳定的PS点及其受临近PS点影响而表示为PS点特征的像元stamps(5,5):调用和(计算每一个干涉图噪声标准差，作为权重)，对选取的PS像元的缠绕相位去除spatiall-uncorrelatedlookangel(DEM)err

9、or.还要调用psmergepatchesps_correct_phasepscalcifgstd主要参数mergeresamplesize0Ifsetto0,noresamplingisapplied.输出结果(去除SULA和masternoise后相位),(噪声标准差)作用对选取的PS像元的缠绕相位纠正spatiall-uncorrelatedlookangel(DEM)error.stamps(6,6):调用，对选取的PS像元进行相位解缠。主要参数Unwrapmethod3D'Unwrapifgindexall'Unwrapprefilterflagy'Unwra

10、ppatchphasen'Ifsetton'(recommended),PSphaseisfilteredusingaGoldsteinadaptivephasefilter.Unwrapgridsize100Unwrapgoldnwin32Unwraptimewin180Unwrapgoldalpha输出结果作用对选取的PS像元的缠绕相位进行相位解缠stamps（7,7）:调用、（减去SCLA后重新解缠时用），估计spatially-correlatedlookangleerror，同时估计出主影像大气影响相位和轨道误差相位（AOE）stamps（8,8）:，滤除其他空间相关

11、项如辅影像AOE。ps_smooth_scla.m主要参数Recalcindexall'Scladerampy'输出结果（SCLA和C主影像大气和轨道误差也在里面），（滤波后的），（滤波后的各项误差相位）作用估计spatially-correlatedlookangleerror，同时估计出主影像大气影响相位和轨道误差相位（AOE）最后输出显示ps_plot（v'，4）4小基线处理如果没有运行前面的PS处理，则需要先在INSARmasterdate目录下读取基线信息加载到matlab工作空间，命令如下mt_extract_infomatlab>>ps_loa

12、d_info确定小基线对在matlab里sbfind:寻找小基线干涉对Partiallythesamedataisusedtoestimatethespectrum,whichmighthavecertainadvantages.cAFHAMMINGHastobeodd.RFTHRESHOLD5RFWEIGHTCORROFFRFOUTMASTERat文件stamps(2,2)：调用，迭代计算每个干涉图中每个候选点的噪声相位，计算每个候选点的丫值。psestgammaquick主要参数max_topo_err5PixelswithuncorrelatedDEMerrorgreaterthanth

13、iswillnotbepicked。filter_grid_size50Candidatepixelsareresampledtoagridwiththisspacingbeforefilteringtodeterminethespatiallycorrelatedphase.filter_weightingP-square'Candidatepixelsareweightedduringresamplingaccordingtothisscheme.clap_win32Wavelengthslongerthanthisarepassedclap_alpha1aTogetherwith

14、the卩term,determinestherelativecontributionofthelow-passandadaptivephaseelementstotheCLAPfilter.clap_beta卩YDetermineswhenconvergenceisreachedanditerationceases.Small_bseline_flagy'Lowcohthresh15输出结果作用迭代计算每个干涉图中每个候选点的噪声相位，计算每个候选点的丫值stamps(3,3):调用，基于计算出的每个候选点的丫值，同时考虑幅度离差选取高相干像元。psselect主要参数selectme

15、thodDENSITY'densityrand2percentrand1输出结果作用基于计算出的每个候选点的丫值，同时考虑幅度离差选取高相干像元stamps(4,4):调用，去除那些只在一些干涉图中相位稳定的高相干点及其受临近高相干点影响而表示为高相干点特征的像元。psweed主要参数weed_alpha8ThetimeseriesphaseforeachpairissmoothedusingaGaussianwindowwithstandarddeviation1/weedalpha.Theoriginalphaseminusthesmoothedphaseisassumedtobe

16、noise.weed_standard_dev10Iftheminimumstandarddeviationisgreaterthanthethreshold,thepixelisdropped.设置成1比较收紧。输出结果作用去除那些只在一些干涉图中相位稳定的高相干点及其受临近高相干点影响而表示为咼相干点特征的像元stamps(5,5):调用和，对选取的高相干像元的缠绕相位去除spatiall-uncorrelatedlookangel(DEM)error.ps_correct_phasepscalcifgstd主要参数mergeresamplesize100Ifsetto0,noresamp

17、lingisapplied.输出结果作用选取的高相干像元的缠绕相位去除spatiall-uncorrelatedlookangel(DEM)error.stamps(6,6):调用，对选取的高相干像元的缠绕相位进行相位解缠，retrievethephaserespecttotheoriginalmasterbyleast-squaresinversion.ps_unwrapsbinvertuw主要参数Unwrapmethod3DQUICK'Unwrapifgindexall'Unwrapprefilterflagy'Unwrappatchphasen'Ifset

18、ton'(recommended),PSphaseisfilteredusingaGoldsteinadaptivephasefilter.Unwrapgridsize200Unwrapgoldnwin32Unwraptimewin180aUnwrapgoldalpha输出结果作用对选取的高相干像元的缠绕相位进行相位解缠retrievethephaserespecttotheoriginalmasterbyleast-squaresinversion.stamps(7,7)和stamps(8,8):调用、和，估计小基线干涉对和PS处理干涉图中的spatially-correlatedlookangleerror，同时估计出主影像大气影