StaMPS软件操作流程

StaMPS软件操作流

程

文档仅供参考

注意：冒号后面都是解释，如果步骤。

没说更改目录就一直是前一步的目录。每个命令都有一个log对

应

1成像处理（ROLPAC）

1.1SAR影像原始数据准备

建立SLC文件夹，将SAR原数据拷贝到SLC目录里，为每景影

像按照获取时间分别建立yyyymmdd文件夹，每景SAR影像原始

数据对应的数据和头文件分别命名为IMAGERYyymmdd和

SARLEADERyyyymmdd。

1.2主影像成像（对应的是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命令介绍

Step_slc_ersMake_raw.pl命令使用make_raw.plorbit_typeleader_filedate

输入参数轨道类型（ODR）、数据头文件

（SARLEADERyyyymmdd）和影像获取时间

yyyymmdd

输出主要结果<date>.raw、<date>.raw.rsc、

<date>_parse_lines.out、shift.out、

shift.out.rsc

作用SAR原始数据预处理，主要完成1）解析数

文档仅供参考

据产品文件，将其转为ROI_PAC的体部格

式；2）填补原始中缺失的行；3）调整采

样窗口的起始时间；4）提取成像参数

Roi_prep.pl命令使用roi_prep,pldateOrbitType

输入参数影像获H又时间（yyyymmdd）、轨道类型

（ODR）

输出主要结果<date>.slc.rscx<date>.roi.in^<date>.roi.out

作用生成sic参数文件

Roi.pl命令使用roi.pldateno1no

输入参数影像获取时间（yyyymmdd）

输出主要结果<date>.slc

作用生成sic图像

roipac2doris命令使用Roipac2corisdate[s]

输入参数影像获取时间（yyyymmdd）

输出主要结果master.res或slave.res

作用生成Doris处理所需的影像信息

1.3选择主影像兴趣区域重新成像Iyao

观测主影像成像图像，确定研究区域，按照研究区域边界扩展

1000像元的范围来编辑co重新运行step_slc_ers,按照兴

趣区域重新成像。step_slc_ers命令用法同12〃粗裁。直接在SLC

文件夹下生成有这个文件，在里面修改，别忘记去掉#号。

1.4建立主影像

精裁，观察主影像成像图像，编辑master_crop.in文件，再次确定

裁剪区域。运行step_master_setup命令建立主影像。〃从

ROI_PAC_SCR中复制master_crop.in到主影像成像的文件夹中。

还要注意精裁的裁剪范围是在粗裁后的影像中得范围不是在初始

成像后的主影像的行列。

Step_master_setup所需参数文件master_crop.in

生成结果辅影像列表：Make_slcs.list

在与SLC同值录下建立文件夹

将、

INSAR_master_datezmaster_date.slc.rsc

master_crop.slc、dem.dorisin、geocode.dorisn、

文档仅供参考

timing.dorisin^master.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

7SX）

第一次对主影像读取，全部读取，不裁剪

2裁剪主影像

cp$MYSCR/master_crop.in或者从安装文件里考

并编辑

3裁剪后主影像再次读取

step_master_read

4读取各个辅影像

文档仅供参考

返回到SLC目录

make_read

2差分干涉处理（DORIS）

2.1提取精密轨道信息（冒号后面的是此步骤的解释说明，不是让

你这么操作）

在insarmaster文件夹下

step_master_orbit_ODR：主影像精密轨道信息提取。〃没有精轨不

运行这步，直接下步

make_orbits：为每幅辅影像建立相应的文件夹，在每个文件夹内

执行step_orbit命令，读取相应辅影像精密轨道信息。

step_master_orbit_ODR运行Doris中的M_PORBITS步骤:

M_PORBITS

主妻参数

M_ORBD1R$SAR_ODR_DIR/ERS1〃精密轨道文件路径

M_ORB」NTERVAL1/屋历时间间隔(s)

M_ORB_EXTRAT1ME6//Timeinsecondsbeforefirstandlastlineto

outputephemerides.

输出结果提出的精密轨道信息保存在master.res中。

作用提取主影像精密轨道信息

step_orbit执行orbit_ers,dorisin命令卡，运行Doris中S_PORBITS

步骤:

S.PORBITS

正要参数

S_ORBDIR$SAR_ODR_DIR/ERS1//精密轨道文件路径

S_ORB_INTERVAL1/唐历时间间隔(s)

S_ORB_EXTRATIME6//Timeinsecondsbeforefirstandlastlineto

outputephemerides.

输出结果提出的藉密轨道信息保存在slave.res中。

作用提取辅影像精密轨道信息

文档仅供参考

2.2主辅影像配准

2.2.1主辅影像粗配准

make_coarse：进入每幅辅影像相应的文件夹，在每个文件夹内执

行step_coarse命令，(coarse.dorisin)进行主辅影像一级轨道配

准和二级像素级配准。

Step_coarse命令执行Doris软件的COARSEORB和

COARSECORR两步骤。

COARSEORB

主要参数不需要参数，依靠读取主辅影像的的master.res和slave.res来完成

输出结果轨道配准结果保存在coreg.out文件中

作用基于主辅影像轨道数据配准，精度达30个像元

COARSECORR

主要参数CC_METHODmagffl〃相关系数计算方法，magfft是先使用ffl变换到

频率域，然后计算相关系数。magspace是直接在空间域计算相关系数

CC_NWIN21〃像素级配准划分窗口的数目，大于5.

CC.WINSIZE1024256〃匹配窗口的大小，处理的时候被改变为奇

数C

CC_INITOFForbit〃利用轨道计算的偏移量作为初始偏移量

输出结果计算出Coarse_correla(ion_lranslation_lines与

coarse_correlation_translation_pixels,保存在coreg.out文件中。

作用利用相关系数法实现像素级疝准。

2.2.2主辅影像精配准

make.coreg：精配准，精配准过程中所有基线距小于设定的

bperp_max值的直接与主影像配准，基线距比较大的和离它最近的

三个辅影像配准。然后利用Hooper提出的配准方法，利用中间过

渡的方法实现所有主辅影像的配准。

make_coreg命令运行Doris软件中的FINE和COREGPM步骤，即三

级子像素精配准和估计配准多项式系数。

FINE

文档仅供参考

主要参数FC_METH0Doversample//精配准计算方法，过采样处理

FC_IN_POSfc_pos.in〃ASCII文件，记录主影像上待配准

点点的坐标

FC_W1NSIZE6464//匹配窗口的大小

FC_ACC88//在方位向和距离向的搜索偏移量

FC」NITOFFcoarsecorr//使用像素级配准后的偏移量

FC_OSFACTOR32U过采样的内插倍数

输出结果配准结果保存在INSAR_date文件夹Fcoreg文件夹中的coreg.out和

CPM_DATA.nl.n2文件至

作用计算主辅影像每个配准点的偏移量

COREGPM

主要参数CPMTHRESHOLD0.30〃经配准后生成大量的匹配点，当

相灵系数大于CPMJHRESHOLD时，才参与平差运算。阈值的大小与File

中窗口的大小有夹，小窗口情况下现值趋于1.0,阈值越大越好。

CPM_DEGREE2〃多项式的阶数，2阶即可。

CPMJ'EIGHTbamler!/最小二乘中权矩阵方法选择，

hauler：权矩阵为相关系数的2倍，linear：权矩阵为相关系数的1

倍c

GPMMAXITER8000〃最小二乘运算的最大迭代次数

输出结果计算出方位向和距离向配准多项式系数，输出的结果保存在coreg文

件夹的coreg.ont文件里c

作用计算配准多项式系数

2.3主影像强度图模拟与地形相位计算（平坦地区地形可不执行）

cdINSAR_master_date〃在INSAR_master_date文件夹，修改

timing.dorisin文件参数

step_master_timing：模拟主影像强度图，计算DEM与主影像的

TimingErroro

make_dems：进入到每个辅影像目录，执行step_dem命令，计算

每个干涉正确地形相位。

step_master_timing运行Doris软件的M_SIMAMP和M_TIMING步

骤：

M_SIMAMP

主要参数SAMJN_FORMATI2/i•型的〃DEM格式，主要包括12(整)，R4

(float)和R8(double)

SAM_IN_DEM/sungt/suzhou_test/DEM/suzhou.dem//DEM路径

SAMJN.SIZE36013601〃DEM行数和列数

SAMJN.DELIA0.〃DEM采样间隔

文档仅供参考

SAM_IN_UL39.99583333115.9958333//DEM左上角纬度和精

度，能够在hdr文件里看到，用写

字板打开可看

SAM」N_NODATA-327687-9999〃无值区的取值

SAM_OUT_FILEniaster_sam.ra\v〃模拟主影像幅度图

SAM_OUT_DEMdem_sam.raw〃裁剪的DEM

输出结果模队的主影像幅度图master_san.raw和根据主影像裁剪的DEM

dem_sain_raw或者demcrop.raw:行歹U号在step_master_timing.log

NcolsDEM和NrowsDEM)

作用模拟主影像幅度图

M_TIMING

主要参数MTE_METHODmagfft〃Magfft是先使用ffl变换到频率域，然

后计算相关系数。Magspace是直接在空间域计算相关系数.Magfft

计算快于magspace。

MTE.NWIN30〃划分窗口的数目

MTE_INITOFF00〃初始偏移量

MTE_WINSIZE40962048//匹配窗口大小

输出结果模拟幅度图与主影像方位向和距离向的TimingError；

Masler_azimuth_timing_error和Master_range_timing_error

模刃,图与主影像的粗配准偏移量Coarse_correlation_transla(ion_lines

与Coarse_coirelation_translation_pixels；

作用计算模拟输度图与主影像的TimingError以及粗配准偏移量。/DEM

偏移量计算。

step.dem运行Doris软件中的COMPREFDEM步骤:

COMPREFDEM

主要参数CRD_METHODtrilinear//求取地形相位方法，有nearest和

trilinear两种。由于是粗DEM,与影像的分辨率不一-致，需要过采样到

同一分辨率下。

CRD_INCLUDE_FEOFF//phasew.r.t.ellipsoid

CRDOUTFILErefdem_11.raw〃由DEM反算得到的地形相位文件

CR[LOLT_DEM_LPdem_radar.raw〃原DEM得到的当前雷达坐标系统下

的就DEM（已栗样）

输出结果地形相位文件refdem_ll.raw和雷达坐标系下的DEM：

dcm_radar.raw

（在各辅影像文件夹中），deni_radar_i.raw（在主影像文件夹中）

这几个文件和主影像行列号一样，floal格式。

作用计算地形相位

2.4辅影像重采样

make_resample：对辅影像按照精配准算出的多项式来重采样。进

入到每个辅影像目录，执行step_resample。（resample.dorisin,自动

加入RS_DBOW10001099915003019）

文档仅供参考

step_resample运行Doris软件中的RESAMPLE步骤:

RESAMPLE

主要参数RS_METHODrcl2p〃重采样的内插核

RS_OUT_FILEslave.res.sic〃重采样后生成的文件名

RS_OUT_FORMATcr4〃输出文件数据类型

输出结果重采样后的辅影像slave_res.slc

作用辅影像重采样

2.5差分干涉处理

makejfgs：主辅影像进行干涉处理.，去除平地相位和地形相位。

进入到每个辅影像目录，运行step_ifg命令。（integram.dorisin）

（make_ifgs_nodem对于平坦地区不用去地形）

step_ifg运行Doris软件中的SUBTRREFDEM、COMPREFPHA和

SUBTRREFPHA步骤，即生成干涉图，利用轨道来计算平

地相位，减去平地相位，减去由粗DEM反算得到的地形相位:

inter

主要参数INTOUTC1NTcint.raw〃生成的干涉图文件，为复数float型

INT_MLLTILOOK11//多视处理系数

输出结果干涉图：cint.raw

作用干涉处理，生成干涉图

COMPREFPHA

主要参数FE_METHODporbits〃采用计算平地相位的方法，精轨

FE_DEGREE3〃位置与平地相位之间的多项式系数

FE_NPOINTS201〃最小二乘法来计算平地相位所用来拟合的点数

输出结果平地相位终端执行过程存于Interferogram.out中

作用根据精轨信息计算平地相位

SUBTRREFPHA

(SUBTRREFDEM)

主要参数SRP_METHODpolynomial〃采用减去平地相位的方法

SRP_OUT_CINTcint.minrefpha.raw〃去除平地效应后的干涉图

输出结果去除平地相位的干涉图：cint.minrefpha.raw（去地形后没有了）

作用去除平地相位影响

SUBTRREFDEM

主要参数SRD_OUT_CINTcint.minrefdem.raw〃去除地形相位后的复干涉图文件

SRD_OFFSET00//基于粗DEM得到的地形相位与原始干涉

合成相位之间的间隔差异。

输出结果去除地形相位后的干涉图文件cint.minrefdem.raw

作用去除地形相位

文档仅供参考

2.6地理编码

step_geo：从雷达坐标系转换到地图坐标系〃只在

INSAR_masterdate文件夹下任意一个辅文件（就是日期那个文件

夹）夹下运行一次就行。

Step_geo运行Doris软件中的GEOCODE步骤:

GEOCODE

主要参数GEO_OUT_I.AM1on.raw〃经度输出文件

GEO_OIT_PHIlat.raw〃纬度输出文件

输出结果像元经纬度信息Ion.raw和lat.raw（就生成一个结果在INSAR文件

夹，就是累积的形变。〉

作用地理编码，参考坐标系转换

3PS处理

3.1准备

注意小基线不是在此步初选点，是在确定小基线对后再初选点，

别弄错了

〃命令是mt_prep0.43250200还是在INSARmasterdate目录入

mt_prep：按照设定的幅度离差阈值、距离向方位向分块和重叠范

围等参数进行处理，选定初始PS候选点。

mt_prep在运行过程中还调用了mt_extract_info和

mt_extract_cands

mt_prep

主要参数da_thresh〃幅度离差阈值，一般取值0.4-0.42,小于比幅

度离差的选为PS候选点，默认为0.4

rgpatches〃距离向上分块，默认为1

az_patches//方位向上分块，默认为1

rg_overlap〃距离向上块间重叠像元数，默认50

az_overlap〃方位向上块间重叠像无数，默认200

输出结果建立PATCH（分几个块，产生几个文件）文件夹，

（直接在工NSARmasterdate目录下）土影像行数和列数

文档仅供参考

width.txt、len.txt,基线分布bperp.1.in,辅影像日期day.1.in,主

影像日期master_day.1.in,干涉对日期ifgday.1.in,文件头信息

heading.1.in,主影像雷达波长lambda.1.in,视角文件

lookangolc.1.in,幅度定标信息calamp.out,DEM参数dcmparm.in

分块情况patch,list,

平均幅度图meanamp.fit,分块范围patch,in和patchnoover.in,

PS候选点pscands.1.ij和pscands.1.da：每个候选点的经纬度信息

pscands.1.11,每个候选点的高程信息pscands.1.hgt,每个候选点的

相位信息pscands.1.ph等

phriame=[,pscands.1.phJ];%foreachPScandida:e,a

floatcomplexvalueforeachifg

ijname=['pscands.1.ij'];%ID#Azimuth#Range#1

lineperPScandidate

bperpname=[*bperp.1.in*];%inmeters11ineper

slaveimage

dayname=[,day.1.in*];%YYYYMMDD,1lineper

slaveimage

ifgdaynamo=[，ifgday.1.in];%YYYYMMDDYYYYMMDD,1

lineperifg

masterdayname=[，masterday.1.in*];%YYYYMMDD

llname=[，pscands.1.IT];%2floatvalues(Ionand

lat)perPScandidate

dariame=[,pscands.1.daJ]：%1floatvalueperPS

candidate

hgtname=[，pscands.1.hgt,];%1floatvalueperPS

candidate

laname=E,lookangle.1.in*];%gridoflookangle

values

headingname=[，heading.1.in*];%satelliteheading

lanbdanamo=[，lambda.1.in*];%wavelength

calnamo=[,calamp.out,];%amplitidocalibrations

widthname=[，width,txt*];%widthofinterferograms

lenname=[，len.txt,];%lengthofinterferograms

作用提取所需信息，幅度校正，建立分块，选取初始PS候选点

3.2stamps处理

〃在INSARmasterdate下先启动matlab

getparm(parmname),无参数时显示所有参数，默认参数从

ps_parms_default.m提取；

setparm(parmname',value),修改某个参数值。

参数结果保存在parms.mat文件里。

文档仅供参考

〃在INSARmasterdate卜，运行stamps

stamps(start-step,end_step,patches_flag,est_gamma_parm)

Patches_flag默认为，y)如果设置为H把所有数据作为一块处

理。

Est_gamma_parm是一个可选参数passedtops_est_gamma_quicko

选y贝iJ利用ps_est_gamma_quick(est_gamma_parm),n利用

ps_est_gamma(est_gamma_parm)。

stainps(l,l)：调用ps_load_initial.m,将PS处理所需各种文件加载

进matlab工作空间

ps_load_inital

所需文件pscands.l.ph,pscands.l.ij,bperp.l.in,day.l.in,ifgday.l.in,

master_day.l.in,pscands.l.ll,pscands.l.da,pscands.l.hgt,

look_angle.l.in,heading.l.in,lambda.1.in，calamp.out,len.txt,

width.txt

输出结果干涉图中每个候选点的基线分布bpl.mat,每个候选点的幅度离差

dal.mat,每个候选点的高度值hgl.mat,每个候选点的视角

laLmat,每个候选点的干涉相位phi.mat,每个候选点的点位、行列

号、经纬度等信息psl.mat

作用将PS处理所需要的各种文件加载到Matlab工作空间，生成相应

的.mat文件

stamps(2,2)：调用ps_est_gamma_quick.m,迭代计算每个干涉图中

每个候选点的噪声相位，计算每个候选点的/值。

ps_est_gamma_quick

主要参数max_topo_err5//MaximumuncorrelatedDEMerror(inin).Pixels

withuncorrclatedDEMerrorgreaterthanthiswillnotbepicked。

filtcr_grid_size5()//Pixelsizeofgrid(inm).Candidatepixelsarc

resampledtoagridwiththisspacingbeforefilteringtodeterminethe

spatiallycorrelatedphase.

filterweightingT-square7/Weightingscheme(PSprobability

squared),theotherpossibilitybeing'SNR'.Candidatepixelsarc

weightedduringresamplingaccordingtothisscheme.

ciap_win32//CLAP(CombinedLow-passandAdaptivePhase)

filterwindowsize。窗口大小单位不是米，是32*32的矩阵

ciap_low_pass_wavclength800//CLAPfilterlow-pass

文档仅供参考

contributioncut-offspatialwavelength(inni).Wavelengthslongerthan

thisarcpassed.反算由频率，因为是低通频域滤波。

clap_alpha1//CLAP。term.TogetherwiththePterm,

determinestherelativecontributionofthelow-passandadaptivephase

elementstotheCLAPfilter.

clap_beta0.3//CLAPPterm

gamma_change_convergence0.005//Thresholdforchangein

gamma

changeinmeanvalueof/.Determineswhenconvergenceisreached

anditerationceases.收敛阈值

Low_coh_thresh31〃就是0.31,较低y的阈值，小于次值时没有ps

点，应用于在P-square方法算weighing

时。

输出结果pml.mat

作用迭代计算每个干涉图中每个候选点的噪声相位，计算每个候选点的

7值

stamps©,3)：调用ps_select.m,基于计算出的每个候选点的7值,

同时考虑幅度离差选取PS像元。

ps_select

主店参数selectmethod'DENSITY'//Otheroption'PERCENT'.

densityrand20//Maximumacceptablespatial

density(perkm2)ofselectedpixelswithrandomphase.点太多的话这个值能够小

些如2

percentrand20//Maximumacceptablepercentageof

selectedpixelswithrandomphase.

输出结果selectl.mat

作用基于计算出的每个候选点的/他,同时考虑幅度图差选取PS像元

stamps(4,4):调用ps_weed.m,去除那些只在一些干涉图中相位稳

定的PS点及其受临近PS点影响而表示为PS点特征的像元。同

时根据噪声标准差去除噪声点。

ps_weed

中耍参数weedalpha8//Smoothingparameterforestimatingphasenoise

distributionforeachpairofneighbouringpixels.Thetimeseriesphaseforeach

pairissmoothedusingaGaussianwndowwithstandarddeviation1/weed

alpha.Theoriginalphaseminusthesmoothedphaseisassumedtobenoise.

weed_standard_dev1.0//Thresholdstandarddeviation.Iftheminimum

standarddeviationisgreaterthantheIhreshold,thepixelisdropped.

输出结果psweed.l.node,psweed2.node,weedl.mat,pm2.mat,ph2.mat,

p$2.mat/hgt2.matzIa2.mat,bp2.mat

作用去除那些只在一些干涉图中相位稳定的PS点及其受临近PS点影响

而表示为PS点特征的像元

文档仅供参考

stamps（5,5）:调用ps_correcl_phase.m和ps_calc_ifg_std.m（计算每一

个干涉图噪声标准差，作为权重），对选取的PS像元的缠绕相位

去除spatiall-uncorrelatedlookangel（DEM）error.还要调用

ps_merge_patches

ps_correct_phase

ps_calc_ifg_std

主族参薮merge_resample_size0//Coarserposting(inm)toresampleto.If

setto0,noresamplingisapplied.

输出结果rc2.mat(去除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'//Prefilterphasebeforeunwrappingto

reducenoise.

Unwrap_patch_phase,nz//UsethepatchphasefromStep3as

prefilteredphase.Ifsetto'n'(recommended),PSphaseisfilteredusinga

Goldsteinadaptivephasefilter.

Unwrap_grid_size100//Resamplinggridspacing.

Unwrap_gold_n_win32//WindowsizeforGoldsteinfilter.

Unwrap_time_win180HSmoothingwindow(indays)forestimating

phasenoisedistributionforeachpairofneighbouringpixels.

unwrapgoldalpha0.8//ValueofaforGoldsteinfilter.

输出结果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.m

ps_smooth_scla.m

ps_scn_filt.m

主旅参数R2calc_index'all'//IndextointerferogramstobeusedintheSCLA

estimation.

Scla_deramp'y'//Ifsetto'yz,aphaserampisestimatedforeach

interferogram.Otheroptionis'nz.

输出结果scla2.mat（SCLA和C主影像大气和轨道误差也在里面），

scla_smooth2.mat（滤波后的,scn2.mat（滤波后的各项误差相

文档仅供参考

位）

作用估计spatially-correlatedlookangleerror,同时估计出主影像大气影

响相位和轨道误差相位（AOE）

最后输出显示

ps_plot（'v',4）

4小基线处理

如果没有运行前面的PS处理，则需要先在INSARmasterdate目录

下读取基线信息加载到matlab工作空间，命令如下

mt_extract_info

matlab

»ps_load_info

4.1确定小基线对

在matlab里sb_find:寻找小基线干涉对

〃sb_find(rho_min,ddiff_max,bdiff_max)运行后显示基线图。也能够

使用命令plot_sb_baselines调出基线图

sb_find

主屡参数

RHO_MIN〃最小相关系数（默认0.50）能够减小这个参数获得更多

小基线对设置成0.1或0.2

DDIFF_MAX〃完全失相关的时间间隔（默认1500天）

BDIFF_MAX〃临界基线值（默认1070）

输出结果small_baselines.list

以相灵系数为条件寻找小基线干涉对

作用

〃终端在INSARmasterdate目录下运行

make_small_baselines：

方位向、距离向滤波，按照选择的小基线对，生成差分干涉图。

文档仅供参考

运行Doris软件中的m_filtazi,s.filtazi,FILTRANGE与

INTERFERO步骤，及其Doris的cpxsum与cpxfiddle命令。并在

INSARmasterdate下生成SMALLBASELINES文件夹

将原来配准后的文件slave_res.slc考入到相应的干涉对组成的文件

夹内分别进行主辅影像方位向滤波与距离向滤波，生成干涉图。

利用PS处理中每个干涉对中的干涉图和去除地形相位的差分干涉

图，结合现有的干涉图，用Doris命令cpxsum处理，生成新的差

分干涉图。当前，对此过程的执行过程还不是很清楚。

m_filtazi,s_filtazi

主要参数AF_BLOCKSIZE4096//fftlengtheachcolumnFFT变换

中龌次距离向中的行数，一般是越大越好。

AF_OVERLAP64//Halfoftheoverlapbetween

consecutivebufferesinazimuthdirection.Partiallythe

samedataisusedtoestimatethespectrum,whichmight

havecertainadvantages.

cAF_HAMMING0.75〃光谱灌波的权值。滤波器的alpha系

飞

AF_OIT_MASTERmaster_azifilt.sic〃滤波后的主影像文件名

AF_OIT_SLAVEslave_azifik.sic〃滤波后的辅影像文件名

AFOUTb'ORMATcr4〃输出文件数据的格式

输出结果方位向滤波后的主影像masterazitilt.sic

方位向滤波后的辅影像slaveazifilt.sic

作用主轴影像方位向滤波

FILTRANGE

主要参数RF_METHODadaptive//距离向滤波的方法

RF_FFTLENGTH256//Lengthofwindowforadaptive

method

RF_NLMEAN15//oddTakemeanoverRF_NLMEANlines

toreducenoiseforpeakestimation.Hastobeodd.

RF_THRESHOLD5//ThresholdonSNRofpeakestimation

toperformrangefiltering

RF_HAMMING0.75//alpha

RFOVERSAMPLE4//Oversamplemasterandslavewith

thisfactorbeforecomputingthecomplexinterferogramfor