航空影像与地图的配准纠正.docVIP

航空影像与地图的配准纠正   1 引言 航空影像由于具有较大的比例尺(，因此包含十分丰富的地面信息。对地面信息的分析、研究和利用必须首先纠正由于像片倾斜和地形起伏等引起的变形，即影像纠正。在航空摄影测量进入“数字化”时代后，航空影像的纠正多采用数字微分纠正的方法[1]。数字微分纠正具有较高的几何精度，但必须首先生成该影像范围的数字地面高程(DEM)。在缺乏数字地面高程的情况下，航空影像的纠正将变得复杂和困难。 为了解决这一问题，有效的途径是利用地形图提供的地形信息。本文提出了将航空影像纠正至地形图的新方法，由于地图是建立在大地坐标系下的，因此纠正的航空影像也将是正射的。为了将航空影像准确地纠正至地形图上，本文使用了图形与图像叠加显示的技术，将地形图一矩形区与影像匹配选择控制点，这一方法不仅加快了选择控制点的速度，更重要的是避免了人工量算的错误。在选择四个控制点后，即将第一次的纠正影像与地图相叠加，以后在叠加的图像中选择控制点，并进行迭代纠正，最后由人根据叠加图像的吻合程度终止纠正过程。 影像与地图的配准 影像与地图的叠加配准是本文纠正试验方法的基础。航空影像的纠正首先要在影像和地形图上选择相应的控制点，然后通过建立的数学模型对影像进行纠正，而控制点的多少与精度直接影响着航空影像纠正的质量。目前，在没有数字地面高程的情况下，大多采用的是人工在地图上选择控制点，其缺点是，由于未纠正影像与地图之间存在比例、旋转和地形差变形等差异，使得选点精度受到一定的影响，并在没有明显点处无法选点或所选点精度不足。本文使用了数字影像与数字地图叠加显示的技术，有效地克服了上述缺点，大大提高了控制点的选点精度。 PC显示模式的设定 航空影像多为8BIT灰度图像，彩色图像也可分为多个8BIT图像表示。因此，本文仅讨论灰度图像的情况。8BIT图像的显示在PC机上以256色显示为最好，在此模式下，既可真实地显示图像，又可实现图像的快速显示[2]。 在256色显示模式下，图形和图像的显示采用的是索引彩色的方式。对于256个索引值的每一个，都对应着一个彩色量，彩色量由绿色亮度、蓝色亮度和红色亮度的总和来表示[3]。只有绿色、蓝色和红色的亮度都相等时，才能表达出灰色。在256色模式下，绿色、蓝色和红色的亮度使用6BIT索引值来表达，所以它们的亮度范围为0～63。因此在PC机256色显示模式下，只能表达出64级灰度。灰度图像的显示只需6个信息位。 地图与图像的叠加显示 在使用地图与图像叠加匹配时,由于只使用地图的空间位置信息,因此只需对地图进行二值化扫描输入,也就是说数字地图只需1个信息位来表达。根据图像只需6个信息位的特点,我们将8位灰度图像进行如下变换: X X X X X X X X 置0 X X X X X X 0 0 7 6 5 4 3 2 1 0 高位 低位 图1 图像灰度变换示意图 Fig.1 Transformation of image gray level 这一操作实际上是对灰度值以4为倍数求整。空出的第0位可将数字地图写入,这样即实现了地图和图像的叠加,又节省了存贮空间。对于以上的数据结构,本文设计了以下的调色板技术来实现图像的单独显示或叠加显示,以便于正确地匹配选点。 当需要关闭地图层时,可将调色板设置中的CASE 1语句SetPalette(i,64,0,0)改为SetPalette(i,i/4,i/4,i/4)。 3 航空影像的纠正 航空影像的纠正可分为两个步骤:初纠正和局部调整。首先,在以上讨论的基础上,准确地选择四个控制点,然后利用这四点,将地图进行第一次纠正,即初纠正。然后将纠正后地图叠加至航空影像上,在叠加的基础上和图像漫游的支持下,可以容易地发现残差点,并准确地调整至正确位置,算出残差值,然后进行局部纠正,若仍有残差,继续匹配选点并纠正,直至满足为止。最后,根据以上记录的控制点,对航空影像进行纠正。具体纠正流程如图2所示(见下页)。 初纠正 初纠正采用了以下的方程式: X=a0+a1x+a2y+a3xy (1) Y=b0+b1x+b2y+b3xy 在经过地图和航空影像的叠加并选择了四个控制点后,可解算出(1)式的系数,利用(1)式将地图纠正至航空影像上,并进行整幅图像的叠加,然后对叠加的图像进行漫游观察来发现残差点。图3a, 3b 给出了初纠正后地图与影像叠加显示和准确匹配的结果。 图3a(左) 地图与航空影像的叠加 Fig.3a(Left) Overlay of map and aerial image 图3b(右) 地图与航空影像的配准 Fig.3b