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處理所需的影像信息選擇主影像興趣區(qū)域重新成像Iyao觀測主影像成像圖像，確定研究區(qū)域，按照研究區(qū)域邊界擴展1000像元的范圍來編輯。重輸出結果提出的精密軌道信息保存在中。作用提取主影像精密軌道信息新運行stepslcers,按照興趣區(qū)域重新成像。stepslcers命令用法同.steporbi

4、t執(zhí)行命令卡，運行Doris中SPORBITS步驟:S_PORBITS主要參數SORBDIR$SARODRDIR/ERS1輸出結果提出的精密軌道信息保存在中。作用提取輔影像精密軌道信息主輔影像配準主輔影像粗配準make_coarse:進入每幅輔影像相應的文件夾，在每個文件夾內執(zhí)行step_coarse命令，（）進行主輔影像一級軌道配準和二級像素級配準。Stepcoarse命令執(zhí)行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):調用，去除那些只在一些干涉圖中相位穩(wěn)定的PS點及其受臨近PS點影響而表示為PS點特征的像元。同時根據噪聲標準差去除噪聲點。psweed主要參數weedalpha8ThetimeseriesphaseforeachpairissmoothedusingaGaussianwindowwithstandardde

8、viation1/weedalpha.Theoriginalphaseminusthesmoothedphaseisassumedtobenoise.weedstandarddevIftheminimumstandarddeviationisgreaterthanthethreshold,thepixelisdropped.輸出結果作用去除那些只在一些干涉圖中相位穩(wěn)定的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):調用，去除那些只在一些干涉圖中相位穩(wěn)定的高相干點及其受臨近高相干點影響而表示為高相干點特征的像元。psweed主要參數weed_alpha8ThetimeseriesphaseforeachpairissmoothedusingaGaussianwindowwithstandarddeviation1/weedalpha.Theoriginalphaseminusthesmoothedphaseisassumedtobe

16、noise.weed_standard_dev10Iftheminimumstandarddeviationisgreaterthanthethreshold,thepixelisdropped.設置成1比較收緊。輸出結果作用去除那些只在一些干涉圖中相位穩(wěn)定的高相干點及其受臨近高相干點影響而表示為咼相干點特征的像元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，同時估計出主影像大氣影