doris使用手册3.docxVIP

六、input.products批处理PROCESS interfero 计算得到干涉图，输出文件是master_slave.cint，输出格式是complex_real4。这里干涉图的定义为：*表示取共扼，.表示逐点相乘，I表示干涉图，M表示主图像，S表示副图像，R表示相位（振幅为1）。其干涉图的相位主值表示为这里是代表取第四象限的值。表示相位图，I表示干涉图。INT_METHOD OVERSAMPLE在计算干涉条纹数据之前，可以进行过采样，采样因子至少为2。过采样实际上是做内插，采样因子为2即在原数据之间内插2个点。这种采样最好在频域内进行，先对原始数据做FFT变化，然后做过采样，再由逆变换后的数据计算干涉条纹。INT_OUT_CINT filename需要输入前面输出的.cint格式的文件。INT_MULTILOOK 5 1为了对干涉条纹进行显示和评估，在过采样之后进行多视处理。它实际上是一种压缩显示方法，显示的结果与已有图件形态差异较小一些。干涉条纹的多视主力中有采取方位向为20，距离向为4的方案，其比例为5：1。PROCESS comprefpha SAR图像中任一点可视为雷达波束与地球表面的交点。它由地面椭球模型、SAR多普勒方程和卫星与目标点之间的距离方程三个关系式来确定。这里选用WGS84参考椭球模型作为高度为零的参考基准面，去除由此椭球面引入的平地相位。图3-5 去平地效应示意图所以，干涉场景区域的参考基准面上的每一点都就满足椭球方程：式中，a、b分别为椭球长半轴和短半轴。由SAR成像原理可知，目标点P满足距离方程：式中，R=R0+col*spc，R0是图像第一列的卫星的斜距，spc表示距离向分辨率。目标点P也满足多普勒方程：式中，fd是多普勒中心频率。式（3-27）、（3-28）、（3-29）组成一个非线性方程组：式中，波长λ、斜距R已知，多普勒中心频率fd用式（3-31）可以计算得到，其中多普勒系数d0、d1和d2的值从头文件中可以得到。卫星S1的fd=d0+d1(t-t0)+ d2(t-t0)2用牛顿迭代法可以得到目标点P的坐标。在用牛顿迭代法时，坐标初值可以根据中心点的大地坐标，用下式求得：其中N为卯酉圈曲率半径，e为第一偏心率。因而，可计算平地相位：PROCESS subtrrefpha 根据前步计算，去除平地效应。输出结果是master_slave.srp，输出格式为complex_real4。I代表干涉图，.代表逐点相乘，是某点的干涉相位。PROCESS coherence 计算图像的一致性，相干性。输出图像为master_slave.coh，输出格式为real4。一般在图像输出的时候也会做多视处理。COH_OUT_COH Outdata/master_slave.coh COH_MULTILOOK 5 1这里的多视处理要和干涉图的多视处理一样。COH_WINSIZE 16 16 计算相关系数的窗口。这里的复相关系数定义为：E表示取期望值，表示取共轭。相关系数则为：七、input.filter_unwrap批处理PROCESS filtphase 这个步骤是可以选择的，对干涉图进行滤波，减少多余的噪声。这个步骤是在SUBTRREFPHA后进行的，方法是空间卷积运算，其中的模板都是写在后面的CARD中。输入图像是master_slave.srp,输出图像是master_slave.cint.filtered，输出格式为complex_real4。PF_OUT_FILE filenamePF_METHOD spatialconv这个方法是通过FFT来进行卷积运算的。它的1D模板是：而2D模板就是PF_KERNELPF_KERNEL，一般来说这个区域是选的越大越好，2D的模板可以在输入文件里进行特定的规定。不过一般奇数模板才有用。PF_KERNEL 5 1 1 1 1 1对于空间卷积，默认为[1 1 1]的3点移动平均值，实部和虚部是分开来计算的。PROCESS unwrap 对干涉图进行解缠。输出图像为master_slave.uint，输出格式为hgt。在DORIS里用的解缠方法也是从外部调用一个SNAPU软件。UW_OUT_FILE Outdata/master_slave.uint解缠后的文件输出。UW_METHOD SNAPHUUW_OUT_FORMAT REAL4输出格式为浮点型。UW_SNAPHU_LOG snaphu.logUW_SNAPHU_coh