专利名称：遥感图像的非线性配准方法
技术领域：
本发明涉及遥感图像的非线性配准方法，采用局部三角网格非线性变换配准方法，对具有几何失真地理图像及非刚性畸变图像实现配准、定位与几何修正，属于计算机信息图像处理技术领域。
配准，即对同一个景物在不同时间、用不同探测器、从不同视角获得的图像，利用图像中公有的景物，通过比较和匹配，找出图像之间的相对位置关系。更准确的说，图像配准的目标就是找到把一幅图像中的点映射到另一幅图像中对应点的最佳变换。由于图像在不同时间、用不同探测器、从不同视角获得时成像条件不同，即使是包含了同一个物体，在图像中物体所表现出来的光学特性(灰度值，颜色值等)、几何特性(外形，大小等)及空间位置(图像中位置，方向等)都会有很大的变化。加之噪声、干扰物体等因素的存在使得图像有很大的差异。
图像配准中的变换，一般是指空间的几何变换。典型的几何变换包括平移、旋转、比例尺变换和非刚性变换(膨胀、收缩、切变、拉变、扭变、等等)。
目前的配准方法仅限于平移、旋转、比例尺变换，对非刚性变换还没有有效的方法。配准时，认为标准图与待配图之间不存在几何失真，只有因几何位置差异和尺度不同所引起的图像不重合，比如，同一图像获取设备在不同时段所拍摄的同一地理区域的图像。由于图像的整体配准只需要在坐标旋转、平移和尺度缩放中进行，配准速度较快。实际中，标准图与待配图之间不仅存在几何位置差异和尺度不同，并且存在一定的几何失真。不同图像获取设备对同一地理区域所拍摄的图像往往存在这种现象。对不同获取设备所获取的同一地区的图像，即使经过完美的定标和几何校正，也常常很难符合线性配准所需要的条件。
为实现这样的目的，本发明采用局部三角匹配网格非线性变换方法，待配图像的配准与几何校正在每个以网点为顶点的三角形局部图像区域进行。为实现快速的配准，采用人机交互选取适配区域，自动寻找同名特征点。利用基准图像与待配图像的若干子图像(即若干局部区域)的配准，实现基准图像与待配图像间的整体配准。利用子图像同名点在整体图像形成的网格和网格内图像的变换，同时实现待配图像的位置配准与几何校正。
本发明的方法包括如下主要步骤1、比例尺变换当基准图像与待配图像具有不同的比例尺时，图像必须首先进行比例尺变换，使两图像具有同样的比例尺，即基准图像与待配图像的每个象素在地面上对应几何尺度相等。比例尺变换采用插值法或重取样法。2、适配区的选择由于在实际应用中，遥感遥测地理图像数据巨大，对整幅图像逐点进行配准计算量极大。采用图像的局部配准方法，人机交互选择一定大小的适配区，可以大大地减少运算时间。
在待配图像上选取N块区域，大小为M*M，块分布应较为均匀。选取块规则如下若第i行选取j块区域，则第i+1行选取j+1或j-1块区域。
在基准图像上分别选取这N块大致所在区域的块，可重叠，块大小为R*R(R＞＞M)。3、灰度归一化处理不同获取设备、不同时间、不同气象情况下所获取的遥感图像，在平均灰度、对比度等方面存在很大的差异。采用灰度归一化处理则是使基准图像与待配图像在平均灰度、对比度方面基本一致，以达到易于配准的目的。4、局部配准采用基于灰度信息的局部配准方法。该方法配准速度快，精度高，对图像的小角度旋转不敏感，适合图像局部区域的配准。将待配图像的适配区域中心、基准图像上配准区域中心作为同名特征点记录下来。5、变换网格生成采用两次网格生成法，得到最终三角网格。先将待配图像上适配区域中心连接成三角网格，然后以每个三角网格中心作为新的网格顶点，连成最终三角网格。6、待配图像的局部非线性变换待配图像的配准与几何校正在每个三角网格内局部区域进行。采用三角网格基函数插值法，对网格内每点进行加权插值变换，获得配准和纠正后图像。
由于传感器、辐射度、大气折射等各个因素的影响，原始遥感图像会产生几何失真及非刚性畸变误差，不能满足量测和定位要求，本发明采用局部三角匹配网格非线性变换方法，对具有几何失真地理图像及非刚性畸变图像实现配准、定位与几何修正。通过比例尺变换、灰度归一化等预处理使得图像在对比度方面达到一致，易于配准。通过采用人机交互选取适配区域，局部配准方法及自动寻找同名特征点实现快速配准。两次三角网格生成，确定网格点的基函数，采用三角网格基函数插值法，对网格内每点进行加权插值变换，可以对图像的几何失真及非刚性畸变有效实现配准和纠正，基函数插值法使得生成图像达到整体连续。
图2为在待配图像上较均匀选取八块区域的示意。
图3为在基准图上相应位置选取八块较大区域的示意。
其中，图3(a)为待配图像取块，图3(b)为基准图像上取块，编号相同为对应块。
图4为第一次划分三角网格的实例。
图5为第二次划分三角网格的实例。
图6为三角网格上基函数插值变换示意图。
本发明所要求的输入数据是下述几种卫星与航空遥感成像的图像数据，格式为标准的BSQ与TIFF1、IKONOS、Landsat、SPOT、FY-1C等不同空间分辨率和重复周期的卫星多通道遥感器的卫星图像；2、OMIS和PHI高光谱航空遥感影像；3、3K×2K彩色CCD相机成像；4、6K全色CCD高空间分辨率航空遥感影像。输出为TIFF图像格式。
本发明采用如
图1所示的遥感图像的非线性配准方案，三角网格生成及非线性变换进行配准与纠正的具体实施步骤如下1、比例尺变换输入基准图像与待配图像，一般具有不同的比例尺，将低分辨率图像进行比例尺变换，使用双线性插值法，使两图像具有同样的比例尺，即基准图像与待配图像的每个象素在地面上对应几何尺度相等。2、适配区交互选择采用人机交互方法选择一定大小的适配区。在待配图像上选取八块区域，大小为M*M，块分布应较为均匀。选取块规则如下若第i行选取j块区域，则第i+1行选取j+1或j-1块区域。如图2所示。
对应待配图像上选取的八块区域，在基准图像上相应位置分别选取八块较大的区域，可互相重叠，块大小为R*R(R＞＞M)。如图3所示。3、灰度归一化处理不同获取设备、不同时间、不同气象情况下所获取的遥感图像，在平均灰度、对比度等方面存在很大的差异。采用灰度归一化处理则是使基准图像与待配图像在平均灰度、对比度方面基本一致，以达到易于配准的目的。设图像灰度函数是I(x，y)，灰度归一化处理如下式所示I(x,y)=I(x,y)-min(I)max(I)-min(I)×(IQ-1)]]>其中min(I)和max(I)表示图像中最小和最大灰度值，IQ表示灰度量化级数。4、局部配准采用基于灰度信息的局部配准方法，确定待配图像上每块在基准图像上的匹配块，并记录下中心位置，作为匹配同名点。5、变换网格生成二次生成网格法，步骤如下a、第一次网格划分。在待配图像上将上下相邻三个中心点连接起来，确定一个三角网格，规则如下同一行三点不连成三角网格，相隔两行的三点不连成三角网格，图像的顶点和边缘中点也作为网格点。如图4所示。
b、设待配图像三角网格顶点(x1，y1)、(x2，y2)、(x3，y3)与基准图像上对应点(X1，Y1)、(X2，Y2)、(X3，Y3)间存在如下变换X1Y1=abcdx1y1+txtyi=1,2,3]]>只要(x1，y1)、(x2，y2)、(x3，y3)三点不共线，变换参数a，b，c，d，tx，ty唯一存在。设三角网格中心为(Ox，Oy)，将变换作为网格中心基函数f(a，b，c，d，tx，ty)。A1～A22为三角网格中心点。
c、第二次网格划分。以A1～A22为中心点，重新生成三角网格，规则同前将上下相邻三个中心点连接起来，同一行三点不连成三角网格，相隔两行的三点不连成三角网格，图像的顶点和边缘中点也作为网格点。如图5所示。6、待配图像的局部非线性变换待配图像的配准与几何校正在每个三角网格内局部区域进行。采用三角网格基函数插值法，对网格内每点进行加权插值变换，获得配准和纠正后图像，输出为tiff格式。三角网格基函数插值法如下设三角网格顶点为Ai，Aj，Ak，三点基函数分别为fi，fj，fk，三角网格内任一点p(x，y)，Li，Ll，Lk分别表示p点到Ai，Aj，Ak的距离，p(x，y)为p(x，y)变换后点，则插值为p′(x,y)=1(Li+Lj+Lk)(Lifi+Ljfj+Lkfk).]]>如图6所示。
权利要求
1.一种遥感图像的非线性配准方法，其特征在于包括如下步骤1)采用插值法或重取样法进行比例尺变换，使基准图像与待配图像的每个象素在地面上对应几何尺度相等；2)采用图像的局部配准方法进行适配区的选择，在待配图像上选取布应较为均匀的若干块区域，在基准图像上选取对应的较大的区域；3)进行灰度归一化处理使基准图像与待配图像在平均灰度、对比度方面基本一致；4)采用基于灰度信息的局部配准方法，将待配图像的适配区域中心、基准图像上配准区域中心作为同名特征点记录下来；5)采用两次网格生成法，先将待配图像上适配区域中心连接成三角网格，然后以每个三角网格中心作为新的网格顶点，连成最终三角网格；6)待配图像的配准与几何校正在每个三角网格内局部区域进行，采用三角网格基函数插值法，对网格内每点进行加权插值变换，获得配准和纠正后图像。
2.如权利要求1所说的遥感图像的非线性配准方法，其特征在于进行网格生成时，同一行三点不连成三角网格，相隔两行的三点不连成三角网格，图像的顶点和边缘中点也作为网格点。
全文摘要
一种遥感图像的非线性配准方法，采用局部三角匹配网格非线性变换方法，对具有几何失真地理图像及非刚性畸变图像实现配准、定位与几何修正。待配图像的配准与几何校正在每个以网点为顶点的三角形局部图像区域进行。为实现快速的配准，采用人机交互选取适配区域，自动寻找同名特征点。利用基准图像与待配图像的若干子图像的配准，实现基准图像与待配图像间的整体配准。利用子图像同名点在整体图像形成的网格和网格内图像的变换，同时实现待配图像的位置配准与几何校正。
文档编号G06T7/60GK1426021SQ0215501
公开日2003年6月25日 申请日期2002年12月19日 优先权日2002年12月19日
发明者庹红娅, 刘允才 申请人:上海交通大学