14.加入光学畸变差改正的二维直接线性变换方法用于建筑立面测绘的研究.doc

加入光学畸变差改正的二维直接线性变换方法 用于建筑立面测绘的研究 王 涛1，闫 利2，尹志永1 (1.天津市测绘院，天津 300381；2.武汉大学测绘学院，湖北 武汉 430079) 摘要：普通数码相机拍摄的建筑立面相片采用二维直接线性变换方法经纠正后可作为绘制建筑立面图的底图，但其精度不高。本文研究在二维直接线性变换模型中加入相机光学畸变差改正以改进成像模型的严格性，进而提高结果的精度，采用的是最小二乘法求解近似的光学畸变差改正项。本文通过试验得到了较高几何精度的纠正影像，证明了该方法的正确性。 关键词：二维直接线性变换；光学畸变差改正；建筑立面测绘；最小二乘法；数码相机 一、引言 随着我国城市规划建设的发展，城市建筑立面测绘的社会需求越来越大，其主要应用方面有三维城市、三维社区系统建设的建筑侧面纹理数据采集，规划审批时的建筑日照分析，等等。建筑立面测绘常用方法有地面常规测量方法和近景摄影测量方法。地面常规测量方法所得结果的几何精度较高，但外业采集工作量大,相应的制作周期也长。采用近景摄影测量方法具有时间短、工作效率高的特点，当采用量测相机也可获得较高精度的结果。 当受设备限制而无量测相机可用，且对建筑物立面图几何精度要求不是特别高时，宜采用直接线性变换方法，相应的设备可选择普通数码相机。直接线性变换（Direct Linear Transformation）解法，简称DLT，因无需相机的内方位元素和摄影时外方位元素的初始近似值，适用于非量测相机（包括普通数码相机）所摄影像的摄影测量处理。 理论上，采用DLT方法处理结果的几何精度不高，若想提高建筑立面图的几何精度，须从提高直接线性变换模型的严格性着手，对成像时的几何变形加以改正，消弱各项误差的影响。对于普通数码相机来说，DLT模型可改正成像时的线性变形，而相机的光学畸变这种非线性变形是普遍存在的，因而可考虑通过改正光学畸变来提高成像模型的精度。本文方法即是基于此思路。 二、加入光学畸变改正的二维直接线性变换方法 1. 二维直接线性变换 二维直接线性变换的关系式为： (1) 建筑物立面可视作二???对象，故建筑立面图测绘可采用二维DLT解法，其计算过程为： 1) 依据最小二乘法，利用4个以上控制点建立误差方程，求解8个系数。 2) 将未知点的像点坐标代入二维DLT关系式，计算其二维平面坐标。 3) 根据计算结果，进行影像重采样得到几何纠正的影像底图，用于绘制建筑立面图。 二维DLT的误差方程式为： (2) 式中， (3) 该解算采用迭代方法，8个系数的初始值可设置为一较小常数，本文试验中其初始值均设为1.0，收敛条件为相邻两次运算差值小于0.1。 2.加入近似光学畸变改正的整体模型 利用上述解算得到的8个系数即可建立二维DLT模型。二维DLT模型可改正数码相机成像过程中的线性误差，如方向、方向单位长度不一致的线性误差，以及坐标轴不垂直所引起的线性误差。但二维DLT模型并未考虑非线性变形的影响，而对普通数码相机来说，其物镜光学畸变差对结果精度的影响须予以考虑。普通数码相机光学畸变差改正模型为: (4) 其中, 为光学畸变差改正项，为未知的像主点坐标，为待定的畸变系数。对普通数码相机来说，借助精密的控制场通过检校可测定此三项未知量，进而得到光学畸变改正模型以提高成像几何精度。但对于科研单位以外人员和生产项目，相机检校有所限制，因此本文采用如下方法来求解其近似值。 在已建立起的二维DLT模型基础上，不妨设 (5) 可见,均为常数，则加入光学畸变差改正的二维DLT模型的观测值方程式为： (6) 易见此关系式是观测值与未知量的非线性关系式，则宜采用泰勒级数展开，再运用改正数累加的方法求解。 采用泰勒级数展开的观测值方程式为： (7)近似值、的计算式为： (8)其中，为未知量的初始值。 运用隐函数偏导求解法得到式(7)中6系数的计算公式为： (9) 其中，为未知量的当前值。 相应的误差方程式为: (10) 依据间接平差原理，求解出未知量的改正数，采用迭代方法反复趋近，以改正数3元素均小于限值为收敛条件，通过改正数累加得到未知量的正解。 至此，建立了加入近似光学畸变改正的二维DLT模型，不妨