Doris应用大全

Doris应用大全FFTW的安装FFTW是由MIT 的 Matteo Frigo和 Steven G. Johnson. 开发的，可用来进行快速傅立叶变换。目前最新版本：v3.0.1 相关网址：/打开终端，进入FFTW路径，依次执行：./configure -prefix=pwd -enable-floatmake make install完成安装。1.1 Doris的安装目前，Doris的最新版本为v3.17，下载地址：http:/enterprise.lr.tudelft.nl/doris/software/download.html解压后进入src路径（cd /doris/src），依次执行： ./configuremakemake install完成安装。1.2 SARtools和ENVISAT_TOOLS的安装分别进入各自路径，执行：make make install完成安装1.3 getorb的安装若linux的版本为fedora4,则需首先安装g77利用getorb可以获得ERS/ENVISAT的精确轨道，Doris的 M_PORBITS 和S_PORBITS要调用getorb 目前最新版本v2.3.0 相关网址：http:/www.deos.tudelft.nl/ers/precorbs/tools/getorb_pack.shtml进入getorb路径，依次执行：makemake install完成安装1.4 snaphu的安装Stanford开发的解缠工具，被在Doris中的UNWRAP调用，目前最新版本v1.4.2 相关网址：/sar_group/snaphu/进入snaphu的src路径，依次执行：makemake install完成安装。1.5 GMT的安装2. 软件运行2.1 设置环境变量在fedora默认的bash下，设置如下环境变量：export PAGER=lessexport EDITOR=viexport EDITOR=gedit /为日后方便，可直接在./etc/bashrc中加入上述两行以永久性设置环境变量。2.2 runDoris为方便用户，可用run命令生成初始化文件。2.2.1在终端输入run g在当前路径下生成如下几个文件：Input.m_inital;input.s_intial;input.resample;input.coregistration;ducts;input.filter_unwrap;input.quicklook正常时在终端返回以下信息：*-andylocalhost $ run -grun: *INFO* generation of inputfiles finished.run: *INFO* inputdir: Inputfilesrun: *INFO* outputdir (stout): Outinforun: *INFO* outputdir (matrices): Outdatarun: *INFO* Using executable in path: dorisrun: *INFO* or (if -d option present): doris.debugrun: *TIP* run -e1; run -s1-*2.2.2输入run e1 调用gedit编辑Input.m_inital(主影像初始化文件)Input.m_inital主要实现以下四步：读取SLC文件-m_readfiles数据剪裁及格式转化m_crop过采样（可选）m_ovs确定精确轨道m_porbits在Input.m_inital中要给出SLC文件volume,leader、data file的路径以及Delft轨道数据文件的路径。Input.m_inital文件中若干参数的设定在图像剪裁选项卡中，需要注意：# -# CROP (master)# -M_DBOW 1 5000 1 1000 / database output window若要处理整幅影像时需要将上面的选项注释掉。# -# OVERSAMPLE (master)# -M_OVS_OUT Outdata/master_ovs.raw / 输出文件名称 M_OVS_OUT_FORMAT ci2 /过采样影像格式 ci2 | cr4.M_OVS_FACT_RNG 2 /距离向过采样率（可采用其他值）M_OVS_FACT_AZI 1 / 方位向过采样率（可采用其他值）M_OVS_KERNELSIZE 16 / 过采样所采用插值核的长度 Raffaele Nutricato曾将距离向采样率设为4来进行多时相分析（multitemporal anylsys ）效果不错。对于PS-Insar，两个方向都设为2比较合理Bert Kamps# -# PORBITS (master)# -从DEOS获得的轨道数据为小写的时，需要改成大写。即将odr.*换成ODR.*否则getorb将无法读取（这跟系统有关）。getorb提供了两种确定轨道的方法三次自然样条插值和三次多项式插值。采用三次自然样条插值时在第一行纪录前最后一行纪录后至少有3个数据点；采用三次多项式插值时，时间间隔可取20-30秒，额外时间可取200秒。Doris默认样条插值。c -for polyfit(3）-c M_ORB_EXTRATIME 300 / Time before first linec M_ORB_INTERVAL 30 / 时间间隔c M_ORB_DUMP 0.05 / dump to ascii filec -for spline-M_ORB_EXTRATIME 6 / Time before first lineM_ORB_INTERVAL 1 / 时间间隔。2.2.3在终端输入run s1执行Input.m_inital中设定的步骤。正常时在终端返回以下信息：*-andylocalhost $ run -s1 doris Inputfiles/input.m_initial Outinfo/stdout.input.m_initialPROGRESS: Interpretation inputoptionsfile finished.total cpu: 0 min 0 secPROGRESS: Finished initializationPROGRESS: Start M_READFILES.PROGRESS: readvolume finished.PROGRESS: readleader finished.PROGRESS: readdat (header info) finished.PROGRESS: Finished M_READFILES.total cpu: 0 min 0.01 secPROGRESS: Start M_PORBITS.(getorb: reading arclist /usr/local/doris/delftorbits/ers2/dgme04/arclist)(getorb: orbit file /usr/local/doris/delftorbits/ers2/dgme04/ODR.112)PROGRESS: getorb: program finished ok.PROGRESS: Finished M_PORBITS.total cpu: 0 min 0.01 secPROGRESS: Orbit: interpolation coefficients computed.PROGRESS: Start M_CROP.PROGRESS: Start cropping slc data.PROGRESS: WRITESLC: 0%PROGRESS: WRITESLC: 10%PROGRESS: WRITESLC: 20%PROGRESS: WRITESLC: 30%PROGRESS: WRITESLC: 40%PROGRESS: WRITESLC: 50%PROGRESS: WRITESLC: 60%PROGRESS: WRITESLC: 70%PROGRESS: WRITESLC: 80%PROGRESS: WRITESLC: 90%PROGRESS: WRITESLC: 100%PROGRESS: Finished M_CROP.total cpu: 0 min 2.21 secPROGRESS: calling preview for cropped masterPROGRESS: Start PREVIEW generation.PROGRESS: SUNraster file created of: Outdata/master.raw (see also file: ./master_mag.ras.sh)PROGRESS: Start M_OVS.PROGRESS: OVERSAMPLESLC: 0%(#)Doris software, $Revision: 3.11 $, $Author: kampes $cpxfiddle: INFO: Number of output pixels: 500PROGRESS: OVERSAMPLESLC: 10%PROGRESS: OVERSAMPLESLC: 20%PROGRESS: OVERSAMPLESLC: 30%PROGRESS: OVERSAMPLESLC: 40%PROGRESS: OVERSAMPLESLC: 50%PROGRESS: OVERSAMPLESLC: 60%PROGRESS: OVERSAMPLESLC: 70%PROGRESS: OVERSAMPLESLC: 80%PROGRESS: OVERSAMPLESLC: 90%PROGRESS: OVERSAMPLESLC: 100%PROGRESS: Finished M_OVS.total cpu: 0 min 17.99 secPROGRESS: calling preview for oversampled masterPROGRESS: Start PREVIEW generation.PROGRESS: SUNraster file created of: Outdata/master_ovs.raw (see also file: ./master_ovs_mag.ras.sh)total cpu: 0 min 17.99 sec .Any given program, when running, is obsolete. - WARNING SUMMARY -There were no messages. Redirected output in file: Outinfo/stdout.input.m_initial view with: run -v1 NOTE: before run -s1, you can do run -ql to get a quicklook processing. which will run step 1 and 2 if not done already.-*2.2.4输入run e2调用gedit编辑Input.s_inital(副影像初始化文件)与run e1同2.2.5 输入run s2执行Input.s_inital中设定的步骤与run s1同2.2.6输入run e3调用gedit编辑Input.coregistration(匹配的初始化文件)Input.coregistration主要实现以下四步：基于轨道的粗匹配coarseorb像素级匹配coarsecorr主影像方位向滤波m_filtazi副影像方位相滤波s_filtazi距离向滤波filtrange（可选）首先根据获取主副影像的轨道进行匹配（coarseorb，没有输入参数），精度30像素左右。然后进行像素级匹配（coarsecorr）有两种方法可供选择，一种是在空间域（space domain）进行的匹配（magspace）。要求窗口尺寸为奇数，设为偶数时会自动转化为奇数。这样做只是为了执行起来方便，因为窗口移动的最小尺寸为1个像素。另一种是在频率域进行的匹配（magfft）。我们知道，空间域的卷积运算与频率域的乘法运算是等价的，通过FFT转化到频率域进行相关性的计算，效率较高。Doris缺省方法为magfftInput.coregistration中若干参数的设定：# -# COARSE COREGISTRATION# -CC_METHOD magfft /缺省方法c CC_METHOD magspace / CC_NWIN 21 / 窗口数量（至少为5）CC_WINSIZE 256 256 / 窗口尺寸c CC_INITOFF 0 0 / 如果没有精确轨道数据，则启用该项设定为0 0CC_INITOFF orbit / 利用基于轨道匹配的结果。窗口数和窗口尺寸的设定，没有什么规律，按默认的也可以，但对于所处理的影像而言，往往不是最佳的。注意：这一步完成后，查看一下logfile，如果初始偏移量与估计值的差异超过1(对ERS-1/2而言，1是正常的）则需要将窗口尺寸增大。# -# AZIMUTH FILTERING# -对主影像进行方位向滤波，通过滤波主副影像不相重叠的部分将被滤掉。在SAR处理过程中如果主副影像采用的多普勒质心频率不相同则会产生不相重叠的部分。本步最好在粗匹配和精匹配之间进行。这样粗匹配得到的方位向的偏移量可用来确定多普勒质心频率的多项式，并有助于精匹配获得较好的结果。c AF_METHOD /AF_BLOCKSIZE 1024 / fftlength each column Length of fft per buffer in azimuth direction. In general, the larger the better.AF_OVERLAP 64 / hbs 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. However it has not been studied yet if taking an overlap is requird. Setting this card to 0 is fastest.AF_HAMMING 0.75 / The weighting of the spectrum in azimuth direction. The filtered output spectrum is first de-weighted with the specified hamming flter, then re-weighted with a (newly centered) one.AF_OUT_MASTER Outdata/master.afAF_OUT_SLAVE Outdata/slave.afAF_OUT_FORMAT cr4 /以上参数按默认值即可# -# RANGE FILTERING (master&slave)# -RF_METHOD porbits / based on orbit geometryRF_SLOPE 0.0 / deg and terrain slopeRF_FFTLENGTH 1024 / power of 2RF_HAMMING 0.75 / alphaRF_OUT_MASTER Outdata/master.rf_orbRF_OUT_SLAVE Outdata/slave.rf_orbRF_OUT_FORMAT cr4 /以上参数按默认值即可#- FINE COREGISTRATION# -该步主要在亚像素精度的水平上计算主副影像对应点位的偏移矢量。首先在分布于整幅影像中的大量窗口中，对于给定窗口，在整幅影像范围内通过计算能量影像（magnitude images）的相干系数来估计偏移量（精度约1个像素），其次，对于上一步所计算的最大相关系数的窗口，进行过采样，以在亚像素的水平上确定相干系数最大值。重复进行上述步骤。下面举个例子说明过采样（距离向过采样）的实现方法：影像是按行读入的，过采样时也按行进行。（定义：line对应于方位向，慢时间；pixel对应于距离向，快时间）给定输入信号，及过采样率（假设取为3）首先，生成一个用填充的与输入信号结构一致的复本：然后，将其与插值核（sinc函数）作卷积，生成输出信号：，其中表示新的采样。c FC_METHOD magfft / fastFC_METHOD oversample / best: oversample complex datac FC_METHOD magspace / same as magfftFC_NWIN 601 / number of smaller windowsc FC_IN_POS highcorr.pos / filename with window positionsFC_WINSIZE 64 64 / size of windowsFC_ACC 8 8 / search window, 2nFC_INITOFF coarsecorr / use result of coarseFC_OSFACTOR 32 / oversampling factorFC_PLOT 0.2 NOBG / GMT plot with magnitude backgroundc /c /STOP 窗口数和窗口尺寸当采用默认值效果不理想时，应根据数据的具体情况进行调整，可以考虑将窗口数增大一倍（个人经验）；窗口尺寸建议使用6464FC_ACC为最大相干系数的搜索精度，建议为8，8，这样搜索区域就为8到8,采用FFT方法时必须为2的幂。FC_OSFACTOR为采样因子（采样率），要想达到1/10之一的匹配精度，建议取为32FC_PLOT为GMT的绘图参数，将绘制大于所给参数的偏移矢量。2.2.6输入run s3将执行Input.coregistration中设定的步骤。正常时，在终端返回以下信息：rootlocalhost doris_tailand# run -s3 doris Inputfiles/input.coregistration Outinfo/stdout.input.coregistrationPROGRESS: Interpretation inputoptionsfile finished.total cpu: 0 min 0 secPROGRESS: Finished initializationPROGRESS: Orbit: interpolation coefficients computed.PROGRESS: Orbit: interpolation coefficients computed.PROGRESS: Start COARSE_ORB.PROGRESS: Coarse precise orbits coregistration finished.PROGRESS: Finished COARSE_ORB.total cpu: 0 min 0.04 secPROGRESS: Start COARSE_CORR.PROGRESS: COARSE_CORR: 0%PROGRESS: COARSE_CORR: 10%PROGRESS: COARSE_CORR: 20%PROGRESS: COARSE_CORR: 30%PROGRESS: COARSE_CORR: 40%PROGRESS: COARSE_CORR: 50%PROGRESS: COARSE_CORR: 60%PROGRESS: COARSE_CORR: 70%PROGRESS: COARSE_CORR: 80%PROGRESS: COARSE_CORR: 90%PROGRESS: COARSE_CORR: 100%PROGRESS: Coarse coregistration based on correlation finished.PROGRESS: Finished COARSE_CORR.total cpu: 0 min 29.35 secPROGRESS: Start FILTAZI (master).PROGRESS: FILTAZI: 0%PROGRESS: FILTAZI: 10%PROGRESS: FILTAZI: 20%PROGRESS: FILTAZI: 30%PROGRESS: FILTAZI: 40%PROGRESS: FILTAZI: 50%PROGRESS: FILTAZI: 60%PROGRESS: FILTAZI: 70%PROGRESS: FILTAZI: 80%PROGRESS: FILTAZI: 90%PROGRESS: Finished FILTAZI for master.total cpu: 2 min 18.9 secPROGRESS: Start FILTAZI (slave).PROGRESS: FILTAZI: 0%PROGRESS: FILTAZI: 10%PROGRESS: FILTAZI: 20%PROGRESS: FILTAZI: 30%PROGRESS: FILTAZI: 40%PROGRESS: FILTAZI: 50%PROGRESS: FILTAZI: 60%PROGRESS: FILTAZI: 70%PROGRESS: FILTAZI: 80%PROGRESS: FILTAZI: 90%PROGRESS: Finished FILTAZI slave.total cpu: 4 min 4.84999 secPROGRESS: Start FILTRANGE (porbits).PROGRESS: filtrange: 0%PROGRESS: Finished FILTRANGE.total cpu: 6 min 57.8 secPROGRESS: Start FINE.PROGRESS: FINE: 0%PROGRESS: FINE: 10%PROGRESS: FINE: 20%PROGRESS: FINE: 30%PROGRESS: FINE: 40%PROGRESS: FINE: 50%PROGRESS: FINE: 60%PROGRESS: FINE: 70%PROGRESS: FINE: 80%PROGRESS: FINE: 90%PROGRESS: FINE: 100%PROGRESS: Fine coregistration finished.PROGRESS: Finished FINE.total cpu: 10 min 9 sectotal cpu: 10 min 9.02002 sec .% How would you rate his incompetence? Unmatched . - WARNING SUMMARY -There were no messages.2.2.7输入run e4调用gedit编辑input.resample(重采样的初始化文件)input.resample主要完成以下三个步骤：coregpm计算匹配参数resample副影像重采样filtrange主副影像距离向滤波# -# COREGISTRATION PARAMETERS MODELING# -该步主要根据步骤FINE估计的偏移量，利用最小二乘方法计算匹配模型（该模型应当具有亚像素级的精度）。可以事先对相干系数设一个阈值来滤除一些数据。计算之后，估计模型和和实际偏移量之间的残差将被绘制。由于干涉图对主副影像之间的匹配误差十分敏感，这一步非常关键。需要不断地迭代以获得最佳的变换模型。Doris采用一个CPM_MAXITER选项来设定自动迭代的次数。CPM_THRESHOLD 0.4 / higher threshold-less windowsCPM_DEGREE 2 / 2d-polynomialc CPM_WEIGHT linear / or nonec CPM_WEIGHT quadratic / or noneCPM_WEIGHT bamler / paper R.Bamler, IGARSS 2000/2004CPM_MAXITER 20 / automated removal, max 20 iterationsc CPM_K_ALPHA 1.97 / automated removal crit. valueCPM_PLOT NOBG / plot w/o magnitude backgroundc CPM_DUMP OFF / no creation of filesCPM_THRESHOLD为相干系数的阈值，这一选项决定了有多少数据可以用于估计多项式的系数。该参数的设定取决于FINE中所采用窗口的大小。若采用的为6464的窗口，设为0.2即可。当然最好还是自己试验若干数值，以获得理想的结果。CPM_DEGREE多项式阶数，建议采用2CPM_WEIGHT最小二乘的定权方法，建议选bamler方法CPM_MAXITER迭代次数选项。如果在FINE中定义了600个窗口，则20次即可。其它参数可以采用默认值。# -# RESAMPLING SLAVE# -根据主影像格网和COREGPM步骤中得到的变换模型对副影像重采样（重采样的过程就是利用插值核在空间域重建副影像原始信号的过程）。该环节非常耗时为了评价重采样的质量，可以将采样后的副影像与主影像再匹配。这时偏移矢量应基本为零。Doris提供了大量的插值核供选择，按缺省设置即可。c - Interpolation Kernel -c RS_METHOD tri / linearc RS_METHOD nearest / nearest neighbourc RS_METHOD cc4p / cubic convolutionc RS_METHOD cc6p / good choice; testedc RS_METHOD ts6p / truncated sincc RS_METHOD ts8p / truncated sincc RS_METHOD ts16p / slower but betterc RS_METHOD knab6p / better than cc6 (theory)c RS_METHOD knab8p / KNAB windowc RS_METHOD knab10p / KNAB windowc RS_METHOD knab16p / KNAB windowc RS_METHOD knab4p / KNAB windowRS_METHOD rc6p / best choice (theory 0.9999)c RS_METHOD rc12p / best choice (theory 0.9999)cc - Output file -RS_OUT_FILE Outdata/slave.resampledRS_OUT_FORMAT cr4 /c RS_DBOW 2000 3000 1 5000 / database output windowRS_SHIFTAZI on / keep this on插值在距离向和方位向上是独立的，对于主副影像重叠区重采样是逐点进行的，每次匹配多项式都要重新生成一次，若按格网进行，效率可能会提高Bert Kampes# -# RANGE FILTERING (master&slave)# -该步主要通过对主副影像距离向滤波去除不重叠区域以减少干涉图中的噪声。不重叠区域产生的原因是卫星视角差异。垂直基线越长，重叠区域越小。RF_METHOD adaptive / estimate local fringe freq.c RF_FFTLENGTH 8 / 8 pixels 160m ground rangec RF_FFTLENGTH 32 / 32 pixels 640m ground rangeRF_FFTLENGTH 128 / 128 pixels 2500m ground rangeRF_OVERLAP 50 / half overlapRF_OVERSAMPLE 4 / master/slave before ifg generationRF_NLMEAN 15 / must be odd (15 pixels 60m)RF_THRESHOLD 5 / SNR spectral peak detectionRF_HAMMING 0.75 / alpha of Hamming windowRF_WEIGHTCORR OFFRF_OUT_MASTER Outdata/master.rf_adaptRF_OUT_SLAVE Outdata/slave.rf_adaptc RF_OUT_FORMAT ci2 / reduce output file sizec /c /STOP / mandatoryRF_METHOD为距离向滤波方法，有adaptive和porbits方法，推荐使用adaptiveRF_FFTLENGTH对于自适应滤波方法，此参数为窗口长度。RF_OVERSAMPLE 过采样因子，缺省为4RF_NLMEAN 行数（方位向）将数行平均为1行以减少噪声，该参数必须为奇数，缺省值为15.RF_THRESHOLD信噪比阈值，缺省为5RF_HAMMING汉明滤波权其它参数按缺省值设置。2.2.8输入run s4正常时在终端返回如下信息：rootlocalhost doris_tailand# run -s4 doris Inputfiles/input.resample Outinfo/stdout.input.resamplePROGRESS: Interpretation inputoptionsfile finished.total cpu: 0 min 0 secPROGRESS: Finished initializationPROGRESS: Orbit: interpolation coefficients computed.PROGRESS: Orbit: interpolation coefficients computed.PROGRESS: Start COREGPM.PROGRESS: Start iteration 0PROGRESS: Start iteration 1PROGRESS: Start iteration 2PROGRESS: Start iteration 3PROGRESS: Start iteration 4PROGRESS: Start iteration 5PROGRESS: Start iteration 6PROGRESS: Start iteration 7PROGRESS: Start iteration 8PROGRESS: Start iteration 9PROGRESS: Start iteration 10PROGRESS: Start iteration 11PROGRESS: Start iteration 12PROGRESS: Start iteration 13PROGRESS: Start iteration 14PROGRESS: Start iteration 15PROGRESS: Start iteration 16PROGRESS: Start iteration 17PROGRESS: Start iteration 18PROGRESS: Start iteration 19PROGRESS: Start iteration 20PROGRESS: finished computation of coregistration parameters.PROGRESS: Finished COREGPM.total cpu: 0 min 0.24 secPROGRESS: Start RESAMPLE.P