一种基于二维图像信息的三维地形测量外文翻译资料

 2022-05-17 21:47:23

英语原文共 15 页,剩余内容已隐藏,支付完成后下载完整资料


一种基于二维图像信息的三维地形测量

【摘要】

河床地形测量是河流模型实验最重要的任务之一,为了提高测量效率和解决传统方法的河床干扰问题,本研究讨论了两种使用电子图像得到测量信息的技术。一种全新改良的图像校准方法被提出用于加强测量的精准性,其是基于两个相机校正而得的。基于透视投影转换,提出了一个基于复合投影的3D重建技术,它的特点是使用一个近似的最大似然估计方法(AMLE)考虑残差计算转换参数的一阶误差传播。此外,从单幅图像的灰度信息出发,对三维地形进行了理论推导。利用漫射光照模型,假设理想灰度值与地形高程值正相关,推导出一种新的封闭公式,用以解释三维地形高程、灰度值、灰度梯度和太阳方位矢量之间的关系。实验结果表明,这两种方法虽然都不完善,但都有一定的优势。

【关键词】

河床地形测量;数字图像测量技术;径向畸变矫正;投影变换;灰度信息变换;

介绍

河床地形测量是一个在河工模型试验中最复杂的任务。传统的地形测量工具,如电子水平仪,电子中转仪,以及全站仪,很难在河床地形测量之中应用。目前,电阻式测距仪是通常使用的装置。它具有良好的适应性,可以测量多种地形。它适用于不同的含沙量和沙质条件。但是,它的低测量的效率和河床干扰问题是它的主要缺点。

为了满足动态测量的需求,各种各样的地形测量方法已经进化。地面激光扫描,空中光探测和测距(LIDAR),多波束声纳,实时动态(RTK),全球定位系统,和全站仪测量不断被开发以满足增加地貌科学研究的兴趣和工程应用。然而,这样的方法更适合于大规模河床模型,甚至是实际的河流。此外,还有很多软件及硬件资源。在最近几年,随着允许的精确校准方法的发展,非量测相机变得越来越可靠,摄影的自动化变得可访问到广泛的用户群。使用在线结构从运动(SFM)计划等。2013 年 Bouratsis 等。利用自由和开源软件的优势从照片制作中创建的河流环境模型的高分辨率数字高程;2013 有人使用的一对商业摄像机来记录河床的演化,并且由一组计算机视觉的一个计算方法和图像处理算法被用来分析视频和重建的瞬时 3D 床的表面上。詹姆斯和罗布森(2012)集成 SFM 和多视点立体(MVS)算法研究在方法测量地形需要很少的专业知识和实现自动化处理。 阿斯特吕克等在 2012 年研究了立体技术测量沙床抬高斜在区波时间尺度。此方法是非侵入性,所以有很高的精度,时间分辨率也很容易让人抓获。此外,鲁等人在 2008 年提出了一些有关平原和高原地形的专题制图相关的研究图象数据是简单地 2D 数据。然而,许多线索可能在单个图像或多个图像中找到,其中可以通过图像场景进行三维重建。最典型的重建方法包括从校正方法的形状(Zhang 等,1999 年),从质地法的形状(福赛思和庞塞,2002),手动交互方法(Shashua,1997)。也有根据一些方法在多个图像进行重建,例如立体视觉的方法, 运动图像序列的方法, 和光度立体法( 李, 1991; 张 Y.J.,2000;Pollefeys 和 GOOL,2002)。还值得一提的是,三维重建应用广泛的还有地形测量。还有在其他领域的一些显著的应用。 Goacute;mez 等人,(2013 年)进行了新的方法用于重建,从多个 3D T 彩色多普勒 3D T 速度场图像。 Mikes 等,(2013 年)集中在扫描光束偏折表面重建,并提出和分析的新的数学方法从而解决了三维表面形状重建问题。柳(2013)提出了一种新的深度提取方法技术与 计算积分成像小透镜阵列。此外,图像分析技术也已用于从 CCD 空间和时间上检索水面仰角(Benetazzo,2006; Gallego 等,2011)。在一定程度上,这些技术是类似于地形测量。他们用重叠的视频图像恢复地形信息。

摄像机失真校正

图像的非线性失真校正是图像测量及三维重建必不可少的步骤。它是 CCD 照相机的几何模型用于确定三维几何的关系对象表面的信息和相应的图像匹配点。由于光学透镜的结构和理想的线性与小孔成像模型之间的不同,从 CCD 拍摄的图像会产生非线性几何失真。这严重影响了 3D 信息的准确性通过 2D 图像获得的。

数学模型

在设计和制造透镜的过程中不可避免地产生在图像畸变。径向变形,偏心变形和薄:一般效果可以通过数学模型来表示棱镜失真。镜头畸变通常可表示为:

其中, xd 和 yd 是原始失真图像坐标, xu 和 yu 是理想的失真图像坐标。xd 和

yd 是在 x 和 y 方向上的失真,分别描述为:

其中第一项是径向失真,所述第二术语是偏心失真,而第三项是薄棱镜变形;K1,K2,P1,P2,S1 和 S2 为非线性失真参数(翁等人,1992;Ahmed 和法拉格, 2005; Wang 等人,2008)。以前的研究(张 Z.Y.,2000;江等人,2001;个人和 Kovacic,2002; Miks 和诺瓦克,2012)指出,透镜径向失真是多比其他像差更加严重。切向失真和畸变棱镜不必考虑。太多的非线性参数不会提高解的精度, 而是将导致其不稳定。

即最常用的描述径向失真多项式模型(PM)可以被写入为:

其中, rd 和ru 是从各距离扭曲的点( xd , yd )和失真点( xu , yu ),以扭 曲的中心 P; ki(ⅰ= 1,2,...)是径向失真参数。 Devernay 和 Faugeras(1995) 表明所述第一阶径向对称畸变参数,K1,可以实现一定的准确性。 为了实现更高的精度,我们使用了一阶和二阶畸变参数,k1 和 k2,测量失真的图像失真程度。由于图像中心是失真(Tordoff 和 Murray,2000)中心的一个很好的近似,我们随便拿图像作为失真中心 P 的中心( xc , yc )。图像校正公式坐标表示为:

然后,改写后的未失真坐标和扭曲的坐标由下式给出:

校准过程

为了找到畸变参数,我们使用 Devernay 和 Faugeras(1995 年)的直线法和一些新的改进。这种方法使用基本属性(Devernay 和 Faugeras,2001 年)是在空间中的每个直线上的投影是一条直线。如果我们可以发现径向失真图像上的变换,这样在空间中的每个直线被认为是转换后的图像中的直线,那么我们就可以估计的图像的畸变的参数。 通过使用此属性,采用一个反复的过程来估计失真参数 k1 和 k2。

边缘检测

图像边缘具有两个属性:方向和大小。边缘像素沿着边缘方向平缓变化;然而在垂直于边缘的方向存在一个巨大的变化。首先,用一个新的算术 Canny 边缘检测器(穆罕默德等人,2013)被用来获得的幅度和边缘的方向,然后将数据由高斯低通滤波器进行处理。其次,非极大值抑制和滞后阈值被用于的精确位置具有子像素精度的优势。这种方法可以提高检测的边缘的分辨率,保护低强度边缘, 增强抗干扰能力。

提取物扭曲线段

边缘检测之后,我们需要提取扭曲线段,这是在三维空间中最有可能的直线。因为某些部分可以被打破边缘检测,当两端之间的距离小于一个 0.5 像素时,我们加入碎段。我们也可以把最低门槛上段长度定为 2 个像素,理由是较短的片段可以包含比关于失真的有用信息更多的噪音。

测量失真错误

为了得到失真参数,计算测量一个扭曲线段的曲率,通过它多少扭曲来分割我们使用的扭曲线段,通过最小二乘逼近以形成一条直线。在一般情况下,直线的方程式接合端点( x1 , y2 )和( x2 , y2 ),表示为:

为一个线段的任何点(x,y)的,我们使用 作为距离, 从点到直线,然后畸变误差是所有的距离的平方的总和。其结果是,该畸变误差为零,当失真段是一个直线。扭曲段的曲率越大,畸变误差将是较大的。

优化方法

我们结合黄金分割搜索方法 GSSM)和二次插值法(QIM)(龚和王,2009 年)来计算最佳失真参数。 QIM 是准确的,但不非常有效的。所以,我们首先得到一个小的 GSSM 值范围,然后使用 QIM 获取,满足我们的精度要求。

整个优化过程是:首先,我们使用该范围的金色部分作为初始点,然后执行由 GSSM 迭代,直到有一点,从获得的函数 f(x)值的差相邻的两步骤迭代。然后一个由迭代 QIM 完成。如果没有太大的区别,那么从两种方法求出的 f(x) 的值,另一迭代由 QIM 完成。再次,如果没有求出的 f(x)的值相差无几从由QIM 相邻两步骤迭代,并且自变量之间的差别非常小,迭代继续进行,直至精度要求得到满足。否则,我们改为 GSSM。我们设定精度为11010 。

由于参数数据减少,并且包含许多异常值,有一个从所述最优参数进一步移动的危险。因为这个原因,我们在 K1 第一只优化,直到它给出了一个稳定的解决方案。然后,K2 加入,并最终全面优化失真参数下进行。

2.3 标定实验结果

要验证校准效果,我们进行两个图像校准实验。该原始图像坐标是不可或缺的。 然而,校准后,图像的坐标非整数,且需要估计的整数坐标的新的灰度值。考虑这个问题,我们使用的三次样条插值法灰度级内插。

图像边缘轮廓如图 1 所示,该两个图像都产生桶形失真。通过 Devernay 和Faugeras 实验结果(1995)的方法示于图 2,计算时间是 65.6 秒和 32.9 秒,失真参数是 K1= 7.9415 107 和 K1= 8.6914 107 ,残余误差分别为 24 和 1.8。该通过该方法实验结果如图 3,计算时间是 55.3 秒和 30.2 秒,失真参数 K1=6.9464times;10-7,K2=4.6261times;10-12 和 K1=7.9267times;10-7,K2=6.2357times;10-13,残余误差分别为 1 和 0.8。结果看起来很合理,表明 Devernay ,Faugeras(1995 年)提出的方法是有效的。我们也提高了精度和失真参数的计算效率,因此,我们可以使用校准图像信息做一个更精确的图像调查。

图 1 两个用于校正失真图像图

(a)从我们的CCD 照相机;(二)从http://www.21tx.com/ DC/ 2007/11/05 /11645.html(五月,2013 访问)

图 2 通过 Devernay 和 Faugeras(1995 年)的原始方法使用参数校准来计算图像

(a):经校准图 1(a)中的结果;

(b):校准图 1(b)中的结果。

图 3 通过该方法使用参数校准来计算图像

(a):经校准图 1(a)中的结果;

(b):校准图(b)中结果。

2D 图像间的映射转换和三维空间坐标

透视变换

考虑到校准的图像可以看作地面中央突起,则存在图像与地面之间的透视投影变换的关系。在摄影测量中,常用的数学模型包括直接线性变换(DLT)(Chen等,1994)和空间后方交会的方法(王,2007)。空间后方交会的方法需要较少的控制点,但它具有较高的要求,在实际迭代过程的初始值要有较高的准确度。所以我们使用 DLT 来计算转换参数。

DLT 是一个直接描述关系角度图像和地面之间的解析表达式。它包含 11 个未知参数(L1-L11),并能提高图像的线性错误。让(U,V)是图像坐标,(X,Y,Z)是相应的三维空间坐标,然后直接线性关系可以表示为:

当给定的 n(nge;6)控制点,我们通常通过最小二乘法来计算这些 11 个参数。传统的方法是:

残留方格的最小总和是在传统方法的目标函数。然而,坐标图像之间的统计属性和相应的三维空间中的坐标被致命地忽略。该溶液的精度不足够大,而且它只能用于低精度摄影。在这项研究中,我们使用一个近似最大似然估计法(AMLE),并考虑残差的一阶误差传播(周和邓,2011 年)来计算转换参数。另外,更简单的迭代算法被应用。

变换参数解决方案基于角法

由于每个点带来了两个方程,用 i 来代表点数,让:

那么我们的控制方程:

是 观 测 值 。 假 定 观 测 是 彼 此 独 立 的 ,

全文共10238字,剩余内容已隐藏,支付完成后下载完整资料


资料编号:[12286],资料为PDF文档或Word文档,PDF文档可免费转换为Word

原文和译文剩余内容已隐藏,您需要先支付 30元 才能查看原文和译文全部内容!立即支付

以上是毕业论文外文翻译,课题毕业论文、任务书、文献综述、开题报告、程序设计、图纸设计等资料可联系客服协助查找。