论二维几何纠正与三维数字微分纠正异同.docVIP

论二维几何纠正与三维数字微分纠正异同 论二维几何纠正与三维数字微分纠正异同 郑天赐 ① (福建省地质测绘院 , 福州 , 350011) 方子岩 (东南大学 , 南京 , 210096 ) 摘 要 用多项式进行遥感影像的二维纠正回避构像的空间几何过程 , 简单易行 , 适用于 各种传感器影像 , 普遍用于 2 幅图像之间的配准 。二维纠正未考虑地形引起的变形影响 , 是一 种近似解法 。结合 D EM 的共线方程三维数字微分纠正 , 是对构像空间几何形态直接描述 , 对传 感器成像时的位置和姿态进行模拟 。当外方位元素解算精度较高时 , 可获得理想的正射影像图 。 关键词 数字图像 几何纠正 多项式 共线方程 D EM 几何位置变换 遥感图像作为空间数据 , 具有空间地理位置的信息 。在应用遥感图像之前 , 必须将其投 影到所处的地理坐标系中 。因此 , 遥感图像的几何处理是遥感信息处理过程中的一个重要环 节 。随着遥感技术的发展 , 来自不同空间分辨率 , 不同光谱分辨率和不同时相的多源遥感数 据 , 形成了空间对地貌观测的影像金字塔 。当处理 、分析和综合利用这些多尺度的遥感数 据 , 进行多源遥感信息的表示 、融合及混合像元的分解时 , 必须保证各不同数据之间几何的 一致性 , 需要进行影像间的几何配准 。同时高分辨率影像的出现对几何处理提出了更高的要 求 。 1 直接法与间接法纠正 像元几何位置变换是按选定的纠正变换函数把原始的数字影像逐个像元地变换到输出影 像相应的位置上去 , 变换方法分直接纠正和间接纠正 (或称正解法与反解法 、直接成图法与 间接成图法) , 2 种方法除考虑影像坐标的出发点不一样外 , 纠正后影像像元的灰度赋值也 不一样 。 直接法方案 , 是从原始图像阵列 (畸变空间) 出发 , 按行列的顺序依次对每个原始像元 点位求其在地面坐标系或输出图像坐标系 (标准空间) 中的正确位置 。 X = Fx ( x , y) ( 1) Y = Fy ( x , y) 式中 : Fx 和 F y 为直接纠正变换函数 。 用直接纠正变换函数计算确定新像元在标准空间中的位置之后 , 该位置处还没有相应的 ① 收稿日期 : 2008204214 作者简介 : 郑天赐 ( 19522) , 男 , 高级工程师 , 测绘专业 。 灰度值 , 应该把原始像元的灰度值复制到 灰度值 , 应该把原始像元的灰度值复制到由上式计算得到的输出图像中的相应点位 ( 像元) 上去 。这种纠正后像元灰度获取办法称为灰度重配置 ( 图 1) 。 图 1 直接成图法 Fig1 1 D iagra m sho wing the method of the direct image f ormat ion 从图 1 可以看出 , 纠正后输出图像的像元其位置已不按阵列排列 , 是不规则像元影像 。 而且有的像元因位置原因灰度会重叠 , 有的又散开了 。可见 , 直接法的优点是灰度重配置简 单快速 。缺点是图像质量变差 。因此 , 在实践中常常采用间接法进行纠正 。 间接法纠正方案 , 是从空白的输出图像阵列 (标准空间) 出发 , 亦按行列的顺序依次对 每个输出像元点位反求其原始图像坐标 (畸变空间) 中的位置 。 x = GX ( X , Y) ( 2) y = GY ( X , Y ) 式中 : GX 和 GY 是间接纠正变换函数 。 用标准空间中的 “空”像元坐标 ( X , Y ) 按间接纠正变换函数计算得到相应的原始图 像 ( 畸变空间) 坐标 ( x , y) , 不一定刚好落在原始图像阵列像元中心 , 因而没有明确的灰 度值 。解决的方法是进行灰度重采样 , 例如采用双线性内插法 ( 图 2) 。 直接法和间接法这 2 种方案本质上并无差别 , 主要不同仅在于所用的纠正变换函数不 同 , 互为逆变换 。其次 , 纠正后像元获得的灰度值的方法 , 对于直接法方案 , 称为灰度重配 置 , 而对间接法方案 , 称为灰度重采样 。由于直接法纠正方案要进行像元的重新排列 , 要求 存储空间大 1 倍 , 计算时间也长 , 所以在实践中通常使用的方案是间接法方案 。 图 2 间接成图法 图 2 间接成图法 ( 图中只用了 5 个像元例子) Fig1 2 Diagra m sho wing the method of the indirect image f ormat ion 2