StaMPS軟件操作流程.doc

精品文檔注意：冒號后面都是解釋，如果步驟。沒說更改目錄就一直是前一步的目錄。每個命令都有一個log對應(yīng)1成像處理（ROI_PAC）1.1 SAR影像原始數(shù)據(jù)準(zhǔn)備建立SLC文件夾，將SAR原數(shù)據(jù)拷貝到SLC目錄里，為每景影像按照獲取時間分別建立yyyymmdd文件夾，每景SAR影像原始數(shù)據(jù)對應(yīng)的數(shù)據(jù)和頭文件分別命名為IMAGERYyymmdd和SARLEADERyyyymmdd。1.2 主影像成像(對應(yīng)的是0級)按照主影像選取原則選取主影像，進入到主影像文件夾master_date，對于ERS-1/2數(shù)據(jù)，運行step_slc_ers命令成像。step_slc_ers主要通過調(diào)用ROI_PAC軟件中的make_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輸入?yún)?shù)軌道類型（ODR）、數(shù)據(jù)頭文件（SARLEADERyyyymmdd）和影像獲取時間yyyymmdd輸出主要結(jié)果.raw、.raw.rsc、_parse_lines.out、shift.out、shift.out.rsc作用SAR原始數(shù)據(jù)預(yù)處理，主要完成1）解析數(shù)據(jù)產(chǎn)品文件，將其轉(zhuǎn)為ROIPAC的內(nèi)部格式；2）填補原始中缺失的行；3）調(diào)整采樣窗口的起始時間；4）提取成像參數(shù)Roi_prep.pl命令使用roi_prep.pl date OrbitType輸入?yún)?shù)影像獲取時間（yyyymmdd）、軌道類型（ODR）輸出主要結(jié)果.slc.rsc、.roi.in、.roi.out作用生成slc參數(shù)文件Roi.pl命令使用roi.pl date no 1 no輸入?yún)?shù)影像獲取時間（yyyymmdd）輸出主要結(jié)果.slc作用生成slc圖像roipac2doris命令使用Roipac2doris date s輸入?yún)?shù)影像獲取時間（yyyymmdd）輸出主要結(jié)果master.res或slave.res作用生成Doris處理所需的影像信息1.3 選擇主影像興趣區(qū)域重新成像I yao觀測主影像成像圖像，確定研究區(qū)域，按照研究區(qū)域邊界擴展1000像元的范圍來編輯c。重新運行step_slc_ers，按照興趣區(qū)域重新成像。step_slc_ers命令用法同1.2./粗裁。直接在SLC文件夾下生成有這個文件，在里面修改，別忘記去掉#號。1.4 建立主影像精裁，觀察主影像成像圖像，編輯master_crop.in文件，再次確定裁剪區(qū)域。運行step_master_setup命令建立主影像。/從ROI_PAC_SCR中復(fù)制master_crop.in到主影像成像的文件夾中。還要注意精裁的裁剪范圍是在粗裁后的影像中得范圍不是在初始成像后的主影像的行列。Step_master_setup所需參數(shù)文件master_crop.in生成結(jié)果輔影像列表：Make_slcs.list在與SLC同目錄下建立文件夾INSAR_master_date,將master_date.slc.rsc、master_crop.slc、dem.dorisin、geocode.dorisn、timing.dorisin、master.res等文件拷入此文件夾作用確定裁剪區(qū)域，建立主影像，確定輔影像列表1.5 輔影像成像在SLC目錄下，運行make_slcs_ers命令，對所有輔影像按照興趣范圍成像。運行過程為逐個進行輔影像目錄，執(zhí)行step_slc_ers命令成像。step_slc_ers用法同1.2./注意其他輔影像成像的時候是按照粗裁的范圍裁剪的，他不進行精裁，由于主影像又進行了精裁，一般輔影像比較大。對應(yīng)1級產(chǎn)品不用成像，直接仿照說明書，讀取即可1建立連接并讀取主影像link_slcs_路徑 （初始數(shù)據(jù)的）cd master_date（即是主影像日期文件夾）/進入主影像日期文件夾step read whole XXX (where XXX is ERS, Envisat, RSAT, or TSX)第一次對主影像讀取，全部讀取，不裁剪2裁剪主影像cp $MY SCR/master_crop.in 或者從安裝文件里考并編輯3裁剪后主影像再次讀取step_master_read4讀取各個輔影像返回到SLC目錄make_read2 差分干涉處理（DORIS）2.1 提取精密軌道信息（冒號后面的是此步驟的解釋說明，不是讓你這么操作）在insarmaster文件夾下step_master_orbit_ODR：主影像精密軌道信息提取。/沒有精軌不運行這步，直接下步make_orbits：為每幅輔影像建立相應(yīng)的文件夾，在每個文件夾內(nèi)執(zhí)行step_orbit命令，讀取相應(yīng)輔影像精密軌道信息。step_master_orbit_ODR運行Doris中的M_PORBITS步驟：M_PORBITS主要參數(shù)M_ORBDIR $SAR_ODR_DIR/ERS1 /精密軌道文件路徑M_ORB_INTERVAL 1 /星歷時間間隔（s）M_ORB_EXTRATIME 6 /Time in seconds before first and last line to output ephemerides.輸出結(jié)果提出的精密軌道信息保存在master.res中。作用提取主影像精密軌道信息step_orbit執(zhí)行orbit_ers.dorisin命令卡,運行Doris中S_PORBITS步驟：S_PORBITS主要參數(shù)S_ORBDIR $SAR_ODR_DIR/ERS1 /精密軌道文件路徑S_ORB_INTERVAL 1 /星歷時間間隔（s）S_ORB_EXTRATIME 6 /Time in seconds before first and last line to output ephemerides.輸出結(jié)果提出的精密軌道信息保存在slave.res中。作用提取輔影像精密軌道信息2.2 主輔影像配準(zhǔn)2.2.1 主輔影像粗配準(zhǔn)make_coarse：進入每幅輔影像相應(yīng)的文件夾，在每個文件夾內(nèi)執(zhí)行step_coarse命令，（coarse.dorisin）進行主輔影像一級軌道配準(zhǔn)和二級像素級配準(zhǔn)。Step_coarse命令執(zhí)行Doris軟件的COARSEORB和COARSECORR兩步驟。COARSEORB主要參數(shù)不需要參數(shù)，依靠讀取主輔影像的的master.res和slave.res來完成輸出結(jié)果軌道配準(zhǔn)結(jié)果保存在coreg.out文件中作用基于主輔影像軌道數(shù)據(jù)配準(zhǔn)，精度達30個像元COARSECORR主要參數(shù)CC_METHOD magfft /相關(guān)系數(shù)計算方法，magfft是先使用fft變換到頻率域，然后計算相關(guān)系數(shù)。magspace是直接在空間域計算相關(guān)系數(shù)CC_NWIN 21 /像素級配準(zhǔn)劃分窗口的數(shù)目，大于5. CC_WINSIZE 1024 256 /匹配窗口的大小，處理的時候被改變?yōu)槠鏀?shù)。CC_INITOFF orbit /利用軌道計算的偏移量作為初始偏移量輸出結(jié)果計算出Coarse_correlation_translation_lines與coarse_correlation_translation_pixels，保存在coreg.out文件中。作用利用相關(guān)系數(shù)法實現(xiàn)像素級配準(zhǔn)。2.2.2 主輔影像精配準(zhǔn)make_coreg：精配準(zhǔn)，精配準(zhǔn)過程中所有基線距小于設(shè)定的bperp_max值的直接與主影像配準(zhǔn)，基線距比較大的和離它最近的三個輔影像配準(zhǔn)。然后利用Hooper提出的配準(zhǔn)方法，利用中間過渡的方法實現(xiàn)所有主輔影像的配準(zhǔn)。make_coreg命令運行Doris軟件中的FINE和COREGPM步驟，即三級子像素精配準(zhǔn)和估計配準(zhǔn)多項式系數(shù)。FINE主要參數(shù)FC_METHOD oversample /精配準(zhǔn)計算方法，過采樣處理FC_IN_POS fc_pos.in / ASCII文件，記錄主影像上待配準(zhǔn)點點的坐標(biāo)FC_WINSIZE 64 64 / 匹配窗口的大小FC_ACC 8 8 /在方位向和距離向的搜索偏移量FC_INITOFF coarsecorr /使用像素級配準(zhǔn)后的偏移量FC_OSFACTOR 32 / 過采樣的內(nèi)插倍數(shù)輸出結(jié)果配準(zhǔn)結(jié)果保存在INSAR_date文件夾下coreg文件夾中的coreg.out和CPM_DATA.n1.n2文件里作用計算主輔影像每個配準(zhǔn)點的偏移量COREGPM主要參數(shù)CPM_THRESHOLD 0.30 /經(jīng)配準(zhǔn)后生成大量的匹配點，當(dāng)相關(guān)系數(shù)大于CPM_THRESHOLD時，才參與平差運算。閾值的大小與File中窗口的大小有關(guān)，小窗口情況下閾值趨于1.0，閾值越大越好。CPM_DEGREE 2 /多項式的階數(shù)，2階即可。CPM_WEIGHT bamler / 最小二乘中權(quán)矩陣方法選擇，bamler：權(quán)矩陣為相關(guān)系數(shù)的2倍，linear：權(quán)矩陣為相關(guān)系數(shù)的1倍。CPM_MAXITER 8000 /最小二乘運算的最大迭代次數(shù)輸出結(jié)果計算出方位向和距離向配準(zhǔn)多項式系數(shù)，輸出的結(jié)果保存在coreg文件夾的coreg.out文件里。作用計算配準(zhǔn)多項式系數(shù)2.3 主影像強度圖模擬與地形相位計算(平坦地區(qū)地形可不執(zhí)行 )cd INSAR_master_date/在INSAR_master_date文件夾，修改timing.dorisin文件參數(shù)step_master_timing：模擬主影像強度圖，計算DEM與主影像的Timing Error。make_dems：進入到每個輔影像目錄，執(zhí)行step_dem命令，計算每個干涉對的地形相位。step_master_timing運行Doris軟件的M_SIMAMP和M_TIMING步驟：M_SIMAMP主要參數(shù)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 3601 / DEM行數(shù)和列數(shù)SAM_IN_DELTA 0.000833333 / DEM采樣間隔 SAM_IN_UL 39.99583333 115.9958333 / DEM左上角緯度和精度，可以在hdr文件里看到，用寫字板打開可看SAM_IN_NODATA -32768/-9999 /無值區(qū)的取值SAM_OUT_FILE master_sam.raw /模擬主影像幅度圖SAM_OUT_DEM dem_sam.raw /裁剪的DEM輸出結(jié)果模擬的主影像幅度圖master_sam.raw和根據(jù)主影像裁剪的DEM dem_sam_raw 或者demcrop.raw（行列號在step_master_timing.log的NcolsDEM和NrowsDEM）作用模擬主影像幅度圖M_TIMING主要參數(shù)MTE_METHOD magfft / Magfft是先使用fft變換到頻率域，然后計算相關(guān)系數(shù)。Magspace是直接在空間域計算相關(guān)系數(shù)。Magfft計算快于magspace。MTE_NWIN 30 / 劃分窗口的數(shù)目MTE_INITOFF 0 0 /初始偏移量MTE_WINSIZE 4096 2048 / 匹配窗口大小輸出結(jié)果模擬幅度圖與主影像方位向和距離向的Timing Error：Master_azimuth_timing_error和Master_range_timing_error模擬圖與主影像的粗配準(zhǔn)偏移量Coarse_correlation_translation_lines與Coarse_correlation_translation_pixels；作用計算模擬幅度圖與主影像的Timing Error以及粗配準(zhǔn)偏移量。/DEM偏移量計算。step_dem運行Doris軟件中的COMPREFDEM步驟：COMPREFDEM主要參數(shù)CRD_METHOD trilinear /求取地形相位方法，有nearest和trilinear兩種。由于是粗DEM，與影像的分辨率不一致，需要過采樣到同一分辨率下。CRD_INCLUDE_FE OFF / phase w.r.t. ellipsoidCRD_OUT_FILE refdem_1l.raw /由DEM反算得到的地形相位文件CRD_OUT_DEM_LP dem_radar.raw /原DEM得到的當(dāng)前雷達坐標(biāo)系統(tǒng)下的新DEM（已采樣）輸出結(jié)果地形相位文件refdem_1l.raw和雷達坐標(biāo)系下的DEM ：dem_radar.raw（在各輔影像文件夾中）， dem_radar_i.raw（在主影像文件夾中）這幾個文件和主影像行列號一樣，float格式。作用計算地形相位2.4 輔影像重采樣make_resample：對輔影像按照精配準(zhǔn)算出的多項式來重采樣。進入到每個輔影像目錄，執(zhí)行step_resample。(resample.dorisin,自動加入RS_DBOW 1000 10999 1500 3019)step_resample運行Doris軟件中的RESAMPLE步驟：RESAMPLE主要參數(shù)RS_METHOD rc12p /重采樣的內(nèi)插核RS_OUT_FILE slave_res.slc /重采樣后生成的文件名RS_OUT_FORMAT cr4 /輸出文件數(shù)據(jù)類型輸出結(jié)果重采樣后的輔影像slave_res.slc作用輔影像重采樣2.5差分干涉處理make_ifgs：主輔影像進行干涉處理，去除平地相位和地形相位。進入到每個輔影像目錄，運行step_ifg命令。（integram.dorisin）（make_ifgs_nodem對于平坦地區(qū)不用去地形）step_ifg運行Doris軟件中的SUBTRREFDEM、COMPREFPHA和SUBTRREFPHA步驟，即 生成干涉圖，利用軌道來計算平地相位，減去平地相位，減去由粗DEM反算得到的地形相位：inter主要參數(shù)INT_OUT_CINT cint.raw /生成的干涉圖文件，為復(fù)數(shù)float型INT_MULTILOOK 1 1 / 多視處理系數(shù)輸出結(jié)果干涉圖：cint.raw作用干涉處理，生成干涉圖COMPREFPHA主要參數(shù)FE_METHOD porbits /采用計算平地相位的方法，精軌FE_DEGREE 3 /位置與平地相位之間的多項式系數(shù)FE_NPOINTS 201 /最小二乘法來計算平地相位所用來擬合的點數(shù)輸出結(jié)果平地相位終端執(zhí)行過程存于Interferogram.out中作用根據(jù)精軌信息計算平地相位SUBTRREFPHA(SUBTRREFDEM)主要參數(shù)SRP_METHOD polynomial /采用減去平地相位的方法SRP_OUT_CINT cint.minrefpha.raw /去除平地效應(yīng)后的干涉圖輸出結(jié)果去除平地相位的干涉圖：cint.minrefpha.raw（去地形后沒有了）作用去除平地相位影響SUBTRREFDEM主要參數(shù)SRD_OUT_CINT cint.minrefdem.raw/去除地形相位后的復(fù)干涉圖文件SRD_OFFSET 0 0 /基于粗DEM得到的地形相位與原始干涉合成相位之間的間隔差異。輸出結(jié)果去除地形相位后的干涉圖文件cint.minrefdem.raw作用去除地形相位2.6 地理編碼step_geo：從雷達坐標(biāo)系轉(zhuǎn)換到地圖坐標(biāo)系/只在INSAR_masterdate文件夾下任意一個輔文件（就是日期那個文件夾）夾下運行一次就行。Step_geo運行Doris軟件中的GEOCODE步驟：GEOCODE主要參數(shù)GEO_OUT_LAM lon.raw /經(jīng)度輸出文件GEO_OUT_PHI lat.raw /緯度輸出文件輸出結(jié)果像元經(jīng)緯度信息lon.raw和lat.raw（就生成一個結(jié)果在INSAR文件夾，就是累積的形變。）作用地理編碼，參考坐標(biāo)系轉(zhuǎn)換3 PS處理3.1準(zhǔn)備注意小基線不是在此步初選點，是在確定小基線對后再初選點，別弄錯了/命令是mt_prep 0.4 3 2 50 200還是在INSAR master date目錄下mt_prep：按照設(shè)定的幅度離差閾值、距離向方位向分塊和重疊范圍等參數(shù)進行處理，選定初始PS候選點。mt_prep在運行過程中還調(diào)用了mt_extract_info和mt_extract_candsmt_prep主要參數(shù)da_thresh /幅度離差閾值，一般取值0.4-0.42，小于此幅度離差的選為PS候選點，默認(rèn)為0.4rg_patches /距離向上分塊，默認(rèn)為1az_patches /方位向上分塊，默認(rèn)為1rg_overlap /距離向上塊間重疊像元數(shù)，默認(rèn)50az_overlap /方位向上塊間重疊像元數(shù)，默認(rèn)200輸出結(jié)果建立PATCH（分幾個塊，產(chǎn)生幾個文件）文件夾，（直接在INSAR master date目錄下）主影像行數(shù)和列數(shù)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，幅度定標(biāo)信息calamp.out，DEM參數(shù)demparm.in分塊情況patch.list，平均幅度圖mean_amp.flt，分塊范圍patch.in和patch_noover.in，PS候選點pscands.1.ij和pscands.1.da；每個候選點的經(jīng)緯度信息pscands.1.ll，每個候選點的高程信息pscands.1.hgt，每個候選點的相位信息pscands.1.ph等phname=pscands.1.ph; % for each PS candidate, a float complex value for each ifgijname=pscands.1.ij; % ID# Azimuth# Range# 1 line per PS candidatebperpname=bperp.1.in; % in meters 1 line per slave imagedayname=day.1.in; % YYYYMMDD, 1 line per slave imageifgdayname=ifgday.1.in; % YYYYMMDD YYYYMMDD, 1 line per ifgmasterdayname=master_day.1.in; % YYYYMMDDllname=pscands.1.ll; % 2 float values (lon and lat) per PS candidatedaname=pscands.1.da; % 1 float value per PS candidatehgtname=pscands.1.hgt; % 1 float value per PS candidatelaname=look_angle.1.in; % grid of look angle valuesheadingname=heading.1.in; % satellite headinglambdaname=lambda.1.in; % wavelengthcalname=calamp.out; % amplitide calibrationswidthname=width.txt; % width of interferogramslenname=len.txt; % length of interferograms作用提取所需信息，幅度校正，建立分塊，選取初始PS候選點3.2 stamps處理/在INSAR master date下先啟動matlabgetparm(parmname),無參數(shù)時顯示所有參數(shù)，默認(rèn)參數(shù)從ps_parms_default.m提??；setparm(parmname,value),修改某個參數(shù)值。參數(shù)結(jié)果保存在parms.mat文件里。/在INSAR master date下，運行stampsstamps（start_step，end_step，patches_flag，est_gamma_parm）Patches_flag 默認(rèn)為y，如果設(shè)置為n把所有數(shù)據(jù)作為一塊處理。Est_gamma_parm是一個可選參數(shù)passed to ps_est_gamma_quick。選y則利用ps_est_gamma_quick（est_gamma_parm），n利用ps_est_gamma（est_gamma_parm）。stamps(1,1)：調(diào)用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.1.hgt，look_angle.1.in，heading.1.in，lambda.1.in，calamp.out，len.txt, width.txt輸出結(jié)果干涉圖中每個候選點的基線分布bp1.mat，每個候選點的幅度離差da1.mat，每個候選點的高度值hg1.mat，每個候選點的視角la1.mat,每個候選點的干涉相位ph1.mat，每個候選點的點位、行列號、經(jīng)緯度等信息ps1.mat作用將PS處理所需要的各種文件加載到Matlab工作空間，生成相應(yīng)的.mat文件stamps(2,2)：調(diào)用ps_est_gamma_quick.m，迭代計算每個干涉圖中每個候選點的噪聲相位，計算每個候選點的值。ps_est_gamma_quick主要參數(shù)max_topo_err 5 / Maximum uncorrelated DEM error (in m). Pixels with uncorrelated DEM error greater than this will not be picked 。filter_grid_size 50 / Pixel size of grid (in m). Candidate pixels are resampled to a grid with this spacing before filtering to determine the spatially correlated phase.filter_weighting P-square/ Weighting scheme (PS probability squared), the other possibility being SNR. Candidate pixels are weighted during resampling according to this scheme.clap_win 32 / CLAP (Combined Low-pass and Adaptive Phase) filter window size。 窗口大小單位不是米，是32*32的矩陣clap_low_pass_wavelength 800 / CLAP filter low-pass contribution cut-off spatial wavelength (in m). Wavelengths longer than this are passed.反算出頻率，因為是低通頻域濾波。clap_alpha 1 /CLAP term. Together with the term, determines the relative contribution of the low-pass and adaptive phase elements to the CLAP filter.clap_beta 0.3 / CLAPtermgamma_change_convergence 0.005 /Threshold for change in gammachange in mean value of . Determines when convergence is reached and iteration ceases. 收斂閾值Low_coh_thresh 31/就是0.31,較低的閾值，小于次值時沒有ps點,應(yīng)用于在P-square方法算weighing時。輸出結(jié)果pm1.mat作用迭代計算每個干涉圖中每個候選點的噪聲相位，計算每個候選點的值stamps(3,3)：調(diào)用ps_select.m，基于計算出的每個候選點的值，同時考慮幅度離差選取PS像元。ps_select主要參數(shù)select method DENSITY /Other option PERCENT.density rand 20 / Maximum acceptable spatial density (per km2) of selected pixels with random phase. 點太多的話這個值可以小些如2percent rand 20 /Maximum acceptable percentage of selected pixels with random phase. 輸出結(jié)果select1.mat作用基于計算出的每個候選點的值，同時考慮幅度離差選取PS像元stamps(4,4):調(diào)用ps_weed.m，去除那些只在一些干涉圖中相位穩(wěn)定的PS點及其受臨近PS點影響而表示為PS點特征的像元。同時根據(jù)噪聲標(biāo)準(zhǔn)差去除噪聲點。ps_weed主要參數(shù)weed_alpha 8 /Smoothing parameter for estimating phase noise distribution for each pair of neighbouring pixels. The time series phase for each pair is smoothed using a Gaussian window with standard deviation 1/weed alpha. The original phase minus the smoothed phase is assumed to be noise.weed_standard_dev 1.0 /Threshold standard deviation. If the minimum standard deviation is greater than the threshold, the pixel is dropped. 輸出結(jié)果psweed.1.node, psweed2.node, weed1.mat, pm2.mat, ph2.mat, ps2.mat, hgt2.mat, la2.mat, bp2.mat作用去除那些只在一些干涉圖中相位穩(wěn)定的PS點及其受臨近PS點影響而表示為PS點特征的像元stamps(5,5):調(diào)用ps_correct_phase.m和ps_calc_ifg_std.m（計算每一個干涉圖噪聲標(biāo)準(zhǔn)差，作為權(quán)重），對選取的PS像元的纏繞相位去除spatiall-uncorrelated look angel (DEM) error. 還要調(diào)用ps_merge_patchesps_correct_phaseps_calc_ifg_std主要參數(shù)merge_resample_size 0 / Coarser posting (in m) to resample to. If set to 0, no resamplingis applied.輸出結(jié)果rc2.mat(去除SULA和master noise后相位), ifgstd2.mat(噪聲標(biāo)準(zhǔn)差)作用對選取的PS像元的纏繞相位糾正spatiall-uncorrelated look angel (DEM) error.stamps(6,6):調(diào)用ps_unwrap.m，對選取的PS像元進行相位解纏。ps_unwrap.m主要參數(shù)Unwrap_method 3D /Unwrapping method.Unwrap_ifg_index all /Index to interferograms to be unwrapped.Unwrap_prefilter_flag y /Prefilter phase before unwrapping to reduce noise. Unwrap_patch_phase n /Use the patch phase from Step 3 as prefiltered phase. If setto n (recommended), PS phase is filtered using a Goldsteinadaptive phase filter.Unwrap_grid_size 100 /Resampling grid spacing. Unwrap_gold_n_win 32 /Window size for Goldstein filter.Unwrap_time_win 180 / Smoothing window (in days) for estimating phase noise distribution for each pair of neighbouring pixels. Unwrap_gold_alpha 0.8 /Value of for Goldstein filter.輸出結(jié)果phuw2.mat, unwrap.1.node, unwrap.2.node作用對選取的PS像元的纏繞相位進行相位解纏stamps(7,7): 調(diào)用ps_calc_scla.m、ps_smooth_scla.m(減去SCLA后重新解纏時用)，估計spatially-correlated look angle error，同時估計出主影像大氣影響相位和軌道誤差相位（AOE）stamps(8,8): ps_scn_filt.m，濾除其他空間相關(guān)項如輔影像AOE。ps_calc_scla.mps_smooth_scla .mps_scn_filt.m主要參數(shù)Recalc_index all /Index to interferograms to be used in the SCLA estimation.Scla_deramp y /If set to y, a phase ramp is estimated for each interferogram.Other option is n.輸出結(jié)果scla2.mat（SCLA和C主影像大氣和軌道誤差也在里面）, scla_smooth2.mat（濾波后的）, scn2.mat（濾波后的各項誤差相位）作用估計spatially-correlated look angle error，同時估計出主影像大氣影響相位和軌道誤差相位（AOE）最后輸出顯示ps_plot(v，4)4 小基線處理如果沒有運行前面的PS處理，則需要先在INSAR master date目錄下讀取基線信息加載到matlab工作空間，命令如下mt_extract_infomatlabps_load_info4.1 確定小基線對在matlab里sb_find:尋找小基線干涉對/sb_find(rho_min,ddiff_max,bdiff_max)運行后顯示基線圖。也可以使用命令plot_sb_baselines調(diào)出基線圖sb_find主要參數(shù)RHO_MIN /最小相關(guān)系數(shù)(默認(rèn)0.50)可以減小這個參數(shù)獲得更多小基線對 設(shè)置成0.1或0.2DDIFF_MAX /完全失相關(guān)的時間間隔(默認(rèn)1500天)BDIFF_MAX /臨界基線值(默認(rèn)1070)輸出結(jié)果small_baselines.list作用以相關(guān)系數(shù)為條件尋找小基線干涉對/終端在INSAR master date目錄下運行make_small_baselines：方位向、距離向濾波，按照選擇的小基線對，生成差分干涉圖。運行Doris軟件中的m_filtazi，s_filtazi，F(xiàn)ILTRANGE與INTERFERO步驟，及其Doris的cpxsum與cpxfiddle命令。并在INSAR master date下生成SMALL BASELINES文件夾將原來配準(zhǔn)后的文件slave_res.slc考入到相應(yīng)的干涉對組成的文件夾內(nèi)分別進行主輔影像方位向濾波與距離向濾波，生成干涉圖。利用PS處理中每個干涉對中的干涉圖和去除地形相位的差分干涉圖，結(jié)合現(xiàn)有的干涉圖，用Doris命令cpxsum處理，生成新的差分干涉圖。目前，對此過程的執(zhí)行過程還不是很清楚。m_filtazi，s_filtazi主要參數(shù)AF_BLOCKSIZE 4096 / fftlength each column FFT變換中每次距離向中的行數(shù)，一般是越大越好。AF_OVERLAP 64 / Half of the overlap between consecutive bufferes in azimuth direction. Partially the same data is used to estimate the spectrum ,which might have certain advantages.c AF_HAMMING 0.75 /光譜濾波的權(quán)值。濾波器的alpha系數(shù)。 AF_OUT_MASTER master_azifilt.slc /濾波后的主影像文件名AF_OUT_SLAVE slave_azifilt.slc /濾波后的輔影像文件名AF_OUT_FORMAT cr4 /輸出文件數(shù)據(jù)的格式輸出結(jié)果方位向濾波后的主影像master_azifilt.slc方位向濾波后的輔影像slave_azifilt.slc作用主輔影像方位向濾波FILTRANGE主要參數(shù)RF_METHOD adaptive /距離向濾波的方法RF_FFTLENGTH 256 /Length of window for adaptive methodRF_NLMEAN 15 / odd Take mean over RF_NLMEAN lines to reduce noise for peak estimation. Has to be odd.RF_THRESHOLD 5 / Threshold on SNR of peak estimation to perform range filteringRF_HAMMING 0.75 / alphaRF_OVERSAMPLE 4 /Oversample master and slave with this factor before computing the complex interferogram for peak estimation.RF_WEIGHTCORR OFF /In peak estimation,weight values to bias higher frequencies.RF_OUT_MASTER master_filtrg.slc /距離向濾波后的主影像文件RF_OUT_SLAVE slave_filtrg.slc /距離向濾波后的輔影像文件RF_OUT_FORMAT cr4 /輸出文件數(shù)據(jù)的格式輸出結(jié)果距離向濾波后的主影像文件master_filtrg.slc距離向濾波后的輔影像文件slave_filtrg.slc作用主輔影像距離向濾波INTERFERO主要參數(shù)INT_OUT_CINT cint.raw / 干涉圖INT_MULTILOOK 1 1 / 多視輸出結(jié)果（干涉圖cint.raw ）去除平地和地形的差分干涉圖cint.minrefdem作用主輔影像生成干涉圖4.2 mt_prepmt_prep：按照設(shè)定的幅度差分離差閾值、距離向方位向分塊和重疊范圍等參數(shù)進行處理，選定初始高相干候選點。/cd SMALL BASELINES進入這個文件夾mt_prep 0.6 3 2 50 200 mt_prep在運行過程中還調(diào)用了mt_extract_info、sb_parms_initial.m和mt_extract_candsmt_prep主要參數(shù)da_thresh /幅度差分離差閾值，一般取值0.6，小于此幅度差分離差的選為高相干候選點，默認(rèn)為0.6rg_patches /距離向上分塊，默認(rèn)為1az_patches /方位向上分塊，默認(rèn)為1rg_overlap /距離向上塊間重疊像元數(shù)，默認(rèn)50az_overlap /方位向上塊間重疊像元數(shù)，默認(rèn)200輸出結(jié)果建立PATCH文件夾，主影像行數(shù)和列數(shù)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，幅度定標(biāo)信息calamp.out，平均幅度圖mean_amp.flt，DEM參數(shù)demparm.in，分塊情況patch.list，分塊范圍patch.in和patch_noover.in，PS候選點pacands.1.ij和pscands.1.da；每個候選點的經(jīng)緯度信息pscands.1.ll，每個候選點的高程信息pscands.1.hgt，每個候選點的相位信息pscands.1.ph等作用提取所需信息，幅度校正，建立分塊，選取初始高相干候選點4.3 SB處理運行過程同PS中stamps的使用，可參考上