一种DEM辅助的SAR图像高精度配准方法与流程

本发明涉及微波遥感
技术领域：
，特别涉及一种dem(digitalelevationmodel数字高程模型)辅助的sar(syntheticapertureradar，合成孔径雷达)图像高精度配准方法。
背景技术：
insar(interferometrysyntheticapertureradar，合成孔径雷达干涉测量)是一种对地观测技术，它通过对雷达卫星在不同观测几何下获取的多幅sar复图像进行干涉处理，得到地表的三维信息。sar图像配准是insar处理流程的第一步，sar图像配准的精度对于后续的处理步骤以及最终的insar产品都至关重要，通常要求sar图像配准的精度达到亚像素级。sar图像配准时，一般利用两幅sar图像，可以任意选取其中一幅图像为主图像，另外一幅图像为辅图像。目前最常用的sar图像配准方法是二维函数法(参见文献：《registrationofinsarcompleximagesandinterferogramenhancement》yangqing-youetc,remotesensingtechnologyandapplication，vol.3,no.2,may1999，第123页至第125页)。它首先采用某种优化过程(如相关系数最大化、条纹清晰度最大化等)选取主图像、辅图像中的控制点，再以控制点坐标及控制点处的配准偏移量构造二维形变函数来描述主辅图像各像元之间的偏移，进而得到每个像元处的配准偏移量。该方法操作简单，在地形平坦区域和短基线情况下能够达到sar图像配准的精度要求。但是，由于在长基线情况下，配准偏移量对场景高程起伏敏感，导致真实的像元之间的偏移量起伏与场景地形起伏呈高相关，显然，此时一个二维函数是不足以对真实的配准偏移量进行高精度拟合的。因此在复杂地形区域和长基线情况下，二维函数法会由于拟合精度不够而造成局部配准精度低的现象。技术实现要素：本发明要解决的技术问题是，提供一种dem辅助的sar图像高精度配准方法，有效解决了在长基线和复杂地形区域情况下，现有sar图像二维函数配准法的局部配准精度低现象。本发明的技术方案是：已知对同一场景获得的两幅sar图像，首先选取任意一幅为主图像，则另一幅为辅图像。通过主图像的轨道参数和对应的dem计算主图像各像元对应的高程；然后在主图像上选取适量的像元作为控制点，通过相关系数法与辅图像进行配准，得到控制点在辅图像中对应的像元位置，即可得到控制点处主图像和辅图像之间的方位向配准偏移量、距离向配准偏移量；然后再根据控制点的坐标、控制点像元对应的高程，以及控制点的方位向配准偏移量、距离向配准偏移量，通过最小二乘拟合，在方位向、距离向分别得到一个与像元坐标、像元对应高程相关的多项式；最后将主图像各像元的坐标以及像元对应的高程代入多项式，即可得到主图像像元与辅图像对应像元之间的方位向、距离向配准偏移量，完成sar图像配准。本发明的有益效果是：通过dem和主图像轨道参数计算主图像各像元的高程，并将像元对应高程引入到计算配准偏移量的多项式中，计算得到更高精度的配准偏移量。这样将配准偏移量与sar图像场景对应的地形关联起来，可以更高精度的拟合真实的配准偏移量，解决了在长基线在长基线和复杂地形区域情况下，现有sar图像二维函数配准法的局部配准精度低问题。附图说明图1是本发明的原理流程示意图；图2为主图像和图像对应高程图；图3(a)为利用本发明方法后所得相干系数图，图3(b)为利用二维函数法后所得相干系数图；图4(a)为利用本发明方法后所得局部干涉条纹，图4(b)为利用二维函数法后所得局部干涉条纹；具体实施方式图1是本发明的原理流程示意图。下面结合图1进行详细说明。步骤一：通过主图像的轨道参数和对应的dem计算主图像各像元对应的高程。将dem地理格网上每一点的地理坐标(b,l,h)(经度，纬度，高程)代入到主图像轨道参数和成像参数条件下的距离多普勒模型中进行处理，具体方法参照论文《utilizationofspacebornesardataformapping》(ieeetransactionsongeoscienceandremotesensing，vol.22,no.2,march1984，第107页至第108页)，可求得dem地理格网的各网点对应在主图像中的位置，从而得到主图像各像元(i,j)(i为方位向坐标，j为距离向坐标)对应高程h(i,j)。步骤二：首先在主图像上选取适量的像元作为控制点，通过相关系数法与辅图像进行配准，得到控制点在辅图像中对应的像元位置，即可得到控制点处主图像、辅图像之间的方位向配准偏移量、距离向配准偏移量。具体方法参照论文《registrationofinsarcompleximagesandinterferogramenhancement》(remotesensingtechnologyandapplication，vol.3,no.2,may1999，第123页至第125页)。其中控制点的分布原则应满足场景对应地形的各高程区间区域均分布有控制点，控制点数量具体视图像大小而定，一般不少于100个。步骤三：然后根据控制点的坐标、控制点像元对应的高程，以及控制点的配准偏移量，通过最小二乘拟合，在方位向、距离向分别拟合一个与坐标、高程相关的多项式。由上述处理步骤可知各控制点的坐标(i′,j′)，控制点像元对应的高程h(i′,j′)，控制点对应的方位向配准偏移量a(i′,j′)，距离向配准偏移量r(i′,j′)。则拟合的配准偏移量计算多项式为：其中apq、bpq、cpq、dpq(p＝0,1,...n；q＝0,1,...n)为多项式的待定系数，将所有控制点的坐标(i′,j′)、高程h(i′,j′)、方位向配准偏移量a(i′,j′)、距离向配准偏移量r(i′,j′)代入公式(1)，通过最小二乘法拟合得到这些待定系数。m、n为拟合多项式阶数，通常在m＝2、n＝2阶时可以得到足够的精度和较高的运算效率。步骤四：将主图像各像元的坐标(i,j)以及对应的高程h(i,j)代入多项式(1)，即可得到主图像各像元与对应辅图像各像元之间的方位向、距离向配准偏移量，完成sar图像配准。本发明的实施例所用的实测数据如图2(a)所示，图2(a)是一幅高程起伏较大的地形区域的原始的sar图像，将该图像作为主图像。此sar图像是从现有的palsar系统获得的。辅图像也是从现有的palsar系统获得的同一场景的原始sar图像，只是主图像和辅图像获取时间不同。本实施例的dem数据采用的是全球公开的srtm-cdem数据，dem格网大小为30米×30米。图2(b)为通过步骤一得到的sar图像各像元对应的高程图。图3(a)为利用本发明方法配准后所得的相干系数，图3(b)为利用现有的二维函数法配准后所得的相干系数图。从图3可以看出利用本发明方法配准后所得相干系数整体高于二维函数法配准所得相干系数。尤其在图的上部分，二维函数法所得的相干系数较低，这说明在此区域二维函数法配准误差较大，而本发明的方法所得相干系数在该区域明显高于二维函数法，且保持较高的相干系数，说明本发明方法的配准精度较高。图4(a)为通过本发明的方法配准后得到的局部干涉条纹，对应图3(a)中白色方框区域内的局部干涉条纹。图4(b)为通过现有的二维函数法配准后得到的局部干涉条纹，对应图3(a)中白色方框区域内的局部干涉条纹。该区域对应着整场景中高程较大的区域。从图4可以看出，通过二维函数法配准所得到的干涉条纹存在条纹不完整、条纹不清晰、条纹之间噪声较多等缺点，而通过本方面的方法配准所得到的干涉条纹则完整、平滑清晰、噪声少等优点。由此可以看出本发明方法配准效果较好。本发明方法图3(a)白框区域相干系数均值0.42897二维函数法图3(b)白框区域相干系数均值0.33060表1表1为分别对图3(a)和图3(b)白色方框区域内的相干系数统计均值。可以看出，在使用本发明的方法进行配准后，相干系数有明显的提升。当前第1页12