基于非固定初始面元的无人机影像点云优化算法

杨楠, 邵振峰, 郭丙轩, 彭哲, 黄蕾

杨楠, 邵振峰, 郭丙轩, 彭哲, 黄蕾. 基于非固定初始面元的无人机影像点云优化算法[J]. 武汉大学学报 ( 信息科学版), 2016, 41(8): 1013-1020. DOI: 10.13203/j.whugis20140750
引用本文: 杨楠, 邵振峰, 郭丙轩, 彭哲, 黄蕾. 基于非固定初始面元的无人机影像点云优化算法[J]. 武汉大学学报 ( 信息科学版), 2016, 41(8): 1013-1020. DOI: 10.13203/j.whugis20140750
YANG Nan, SHAO Zhenfeng, GUO Bingxuan, PENG Zhe, HUANG Lei. Point Cloud Optimization for UAV Image Based on Non-fixed Initial Patch[J]. Geomatics and Information Science of Wuhan University, 2016, 41(8): 1013-1020. DOI: 10.13203/j.whugis20140750
Citation: YANG Nan, SHAO Zhenfeng, GUO Bingxuan, PENG Zhe, HUANG Lei. Point Cloud Optimization for UAV Image Based on Non-fixed Initial Patch[J]. Geomatics and Information Science of Wuhan University, 2016, 41(8): 1013-1020. DOI: 10.13203/j.whugis20140750

基于非固定初始面元的无人机影像点云优化算法

基金项目: 

国家重大设备专项 2012YQ16018505

国家科技支撑计划 2013BAH42F03

教育部新世纪优秀人才资助计划 NCET-12-0426

省部产学研结合项目 2012B090500016

详细信息
    作者简介:

    杨楠, 博士生, 研究方向为密集匹配、多视立体匹配与三维重建。candyyanghot@126.com

  • 中图分类号: P231

Point Cloud Optimization for UAV Image Based on Non-fixed Initial Patch

Funds: 

National Key Foundation for Exploring Scientific Instrument, No.2012YQ16018505; the National Key Technology Support Program 2012YQ16018505

the National Key Technology Support Program 2013BAH42F03

the New Century Excellent Talents in University NCET-12-0426

Special Project on the Integration of Industry, Education and Research of Guangdong Province 2012B090500016

More Information
    Author Bio:

    YANG Nan,PhD candidate, specializes in dense matching, multi-view stereo matching and 3D reconstruction.E-mail:candyyanghot@126.com

  • 摘要: 获取平台的不稳定性容易导致无人机影像的几何变形增大。如何高精度匹配这类影像是当前摄影测量与遥感领域的研究热点之一。针对这个问题,提出基于非固定初始面元的无人机点云优化算法,采用差分代替微分计算地物表面局部正切平面的近似法向量,以此为初值建立初始物方面元进行匹配计算,并用两组数据进行实验验证。结果表明,基于非固定初始面元的无人机影像点云优化算法改进了基于物方面元的最小二乘匹配方法,优于基于面元的多视立体匹配(patch-based multi-view stereo matching and reconstruction,PMVS)中的点云优化方法,提高了点云优化的效率和精度。
    Abstract: UAV platform instability causes large geometric deformation in UAV images and unsatisfactory matching accuracy. To address this problem, Point Cloud Optimization for UAV images based on a non-fixed initial patch algorithm is proposed in this paper. By using the difference to calculate approximate normal vector of local tangent plane instead of the differential, and by using this approximate normal vector to establish the initial patch. Two groups of images in campus of Northwestern University and Yangjiang Area were used to test this method. Experimental results show that the Point Cloud Optimization for UAV Image based on a non-fixed initial patch algorithm improved the Patch-based Least Squares Image Matching method, and was superior to the optimization method in PMVS. This method increased the efficiency and the accuracy of point cloud optimization.
  • 在传感器技术高速发展的信息化测绘时代,利用无人机低空遥感平台快速地获取地表数据已经成为研究的热点[1]。近年来,利用激光扫描系统获取地物三维数据得到了摄影测量与遥感领域的广泛关注,但激光扫描系统作为一种新兴的遥感方式,价格较为昂贵,在短期内难以大范围普及使用。采用影像密集匹配来重建三维空间几何模型的方式成本较低,且包含了激光点云数据所欠缺的光谱信息[2, 3]。运用无人机作为平台的低空遥感具有飞行高度低、飞行成本低廉、响应速度快等优势[4],弥补了卫星遥感和普通航空摄影测量易受天气影响不能实时获取地面影像的缺陷,具有广阔的应用前景。

    无人机平台由于体积小,飞行过程中易受到气流的影响,导致飞行姿态不稳定,影像可能出现严重的变形,大大增加了影像匹配的难度。文献[5]提出的基于物方面元的最小二乘匹配方法将物方面元的概念引入到平差模型,从而减少了未知数的个数,在无人机影像的匹配过程中取得了较好的效果。然而,最小二乘匹配结果对初值的要求极高。本文提出的基于非固定面元的无人机点云优化算法是对物方面元最小二乘匹配方法的改进,旨在提高匹配效率,提高点云精度。

    本文先讨论多视立体匹配(PMVS)的点云优化算法,再分析基于物方面元的最小二乘匹配算法,并提出优化算法。

    在Seitz[6]等人提供的多目立体重构算法的公共测试平台中,Furukawa[7]等人提出的PMVS算法取得了最佳效果。PMVS算法的核心思想是在物方构建近似于物体表面的局部正切平面(patch)。如图 1(b)所示,蓝色符号表示算法中定义的patch,其法向量(蓝色箭头符号)指向影像(image)方向。如图 1(a)所示,patch的定义包括中心c(p),法向量n(p),参考影像R(p),其中,p表示patch。patch的一条边与参考相机(即拍摄参考影像的相机)的x轴平行。patch上有m像素×m像素的格网。

    图  1  PMVS中的局部正切平面
    Figure  1.  Patch in PMVS

    patch中心点就是两张影像上对应特征点前方交会求出的物方模型点。初始patch平面的法向量norm是patch中心点坐标指向参考影像摄影机中心的单位向量。patch优化的过程就是调整法向量norm和patch中心点的过程。

    PMVS定义的patch网格的边平行于参考影像摄影机坐标系的x、y轴。patch网格大小为patch沿着x、y轴的两个方向步进的向量pxaxispyaxis。PMVS在计算的过程中通常取5×5或7×7的格网。

    将求得的初始patch分别投影到参考影像和搜索影像上,根据patch投影到每张影像上的投影区域q(p, I1)、q(p, I2)计算相关系数,进而衡量影像区域间的匹配程度。

    PMVS简化自由度的思想是通过影像与patch间的投影关系将patch的6个自由度(包括中心点坐标c(p)=(X, Y, Z),法向量n(p)=(a, b, c))简化为3个, 即物点沿摄影光线移动的距离PP′和法向量方向角αβ

    设patch法向量n(p)在摄影坐标系中的坐标为(a, b, c),patch法向量在摄影坐标系中的方向角为αβ,则abc的计算公式为:

    (1)

    优化patch的过程是按照空间几何关系不断调整patch的3个自由度,使得patch网格表面投影到影像上的影像窗口相关系数最大。PMVS将patch投影到参考影像与搜索影像所得的影像窗口的相关系数定义为fi, i=1, 2, 3…,i是patch搜索影像集合的序号,不包含参考影像。定义目标函数:

    (2)

    式中,自变量为patch的三个自由度,n为搜索影像个数。优化patch的过程由此转化为求取一个极值点使得函数f(z, α, β)的值最小的过程。在PMVS算法中,使用下山单纯形法(Nelder-Mead)[8]对函数f (z, α, β)进行求解。

    下山单纯形法是一种无须求解函数导数来寻找函数极值的方法。它的核心思想是,首先构造一个初始的包含给定点的单纯形,然后通过映射、收缩和扩展三种方式不断替换函数值最大的点,直到所掌握的单纯形收缩到一个点为止。

    由下山单纯形法的定义可知,在无法估计近似的“谷底”点(函数最小值)时,算法的时间复杂度非常大。而文献[9]提出的基于物方面元的最小二乘匹配方法可以加快点云优化过程中的收敛速度,减少迭代次数。

    最小二乘匹配[10]方法由于具有灵活、可靠、高精度的特点,广泛地应用于特征提取、影像匹配、目标定位等方面[9, 11, 12],并很快推广到多点和多片的情形[13]。直至今日, 最小二乘匹配仍是处理许多重要摄影测量任务的一个强大的解决方案。众多学者针对最小二乘匹配的特点提出改进方法[14, 15]。Baltsavias将共线方程约束引入最小二乘匹配,提出了共线条件约束的多片最小二乘匹配算法[16](multi-photo geometrically constrained matching, MPGC)。MPGC算法不但可以在影像匹配的同时直接得到对应匹配点前方交会的三维坐标,而且可以提高最小二乘匹配的稳定性。文献[5]在传统MPGC算法的基础上引入物方面元的概念,利用影像与地面间的几何关系,重新定义多片最小二乘匹配,大大简化了MPGC算法的自由度。文献[17]提出匹配区域内多中心平行投影的方法,通过构建有理多项式系数(RPC)模型下的投影方程,将基于物方面元的最小二乘匹配算法推广到多中心投影推扫式成像的影像匹配中,取得了较好的匹配效果。

    由于低空遥感影像多采用单中心投影,因此本文着重介绍基于共线条件方程的物方面元最小二乘匹配。如图 2所示,设三维点P (Xc, Yc, Zc)投影到参考影像上的像点为p0c (x0c, y0c),则以点p0c为中心的影像窗口是以点P为中心的某个物方面元f (X, Y, Z)的投影。将物方面元投影到搜索影像上,所得到的影像窗口与参考影像窗口中的像素是一一对应的。令物方面元法向量为(a, b, c), 则物方面元的方程为:

    图  2  物方面元与光束交点
    Figure  2.  Intersection of the Object Patch and the Beam

    (3)

    式中,(a, b, c)可由方向角(α, β)表示(式(2))。

    令(u0, v0, w0)为像点p0(x0, y0)在像空间辅助坐标系中的坐标,(XS0, YS0, ZS0)为投影中心在物空间坐标系中的坐标,λ为投影系数,得三维点坐标(X, Y, Z):

    (4)

    联立式(4)、式(5), 有:

    (5)

    由于面元中心点P (Xc, Yc, Zc)也在物方面元上,令λc为投影系数,(u0c, v0c, w0c)为P在参考影像中的像点p0c(x0c, y0c)(影像窗口中心点)在像空间辅助坐标系中的坐标,则有:

    (6)

    由共线方程[18]可求得物方面元上的点投影到搜索影像上的像点坐标(xi, yi)。根据最小二乘匹配原理建立误差方程:

    (7)

    式中,h0ih1i为第i张搜索影像到参考影像的辐射畸变系数;dh0idh1i为辐射畸变系数的改正值;g0(x0, y0)为参考影像上影像窗口中点的像素亮度值;gi(xi, yi)为像点(xi, yi)的像素亮度值;(∂gi/∂xi, ∂gi/∂yi)为亮度值在xy方向的导数值,由于影像像素间间隔较小,可用差分代替微分来计算。(dxi, dyi)为像点坐标的改正值,联立方程(2)、(5)、(6)、(7),像点坐标(xi, yi)可由物方面元中心点高程坐标Zc与物方面元法向量方向角(α, β)表示:

    (8)

    文献[5]中,Zc的初值由初始匹配点前方交会求得,初始方向角(α, β)固定为(0°, 90°),即水平面元。相比PMVS的点云优化算法,基于物方面元的最小二乘匹配方法同样将几何变形未知数化简为三个,并且基于物方面元的最小二乘匹配方法将改正值作为未知数,可以提高迭代的收敛速度。同时,基于物方面元的最小二乘匹配方法将辐射畸变系数引入误差方程,提高了匹配精度。

    最小二乘匹配方法的结果很大程度上取决于给定的初值是否接近真值,接近真值的初值可以减少迭代次数,提高匹配精度;否则,初值则很可能造成迭代算法不收敛。基于物方面元的最小二乘匹配方法中,初始面元由三维点坐标(Xc, Yc, Zc)、法向量(a, b, c)以及面元尺寸μ确定。通过物方面元与影像间的几何关系,算法将人为设定的参数减少到4个,即初始面元中心点高程坐标Zc、面元法向量方向角(α, β)与影像窗口尺寸μ。其中,初始物方面元的中心点为匹配点前方交会所求得的三维坐标,其高程坐标Zc与真实坐标相近。物方面元的几何意义是地物表面的局部切平面,其法向量方向角(α, β)给定的初值是垂直于水平面的单位向量。但在实际匹配计算过程中,由于地形起伏是不定的,不同位置三维点的局部最小切平面的方向也不是固定的。如果能获得每个物点近似的局部正切平面法向量并以此作为初值,与垂直于水平面的单位向量相比更接近真实地形。将所得到的面元投影到匹配影像上,所得的匹配窗口相对于固定面元变形更小,可以减少迭代次数,提高匹配的效率和精度。

    在单点最小二乘匹配中,每个三维点被视为独立的个体,然而在物空间坐标系中,立体匹配所生成的点云数据可以被拟合成一个近似连续的曲面,即被拍摄地物的表面。因此,地物表面在物空间直角坐标系中可以近似地用一个曲面方程表示[19]

    (9)

    过曲面上一点处所有曲线的切线组成一个平面,即曲面在该点处的切平面,其法向量表示为:

    (10)

    由于点云数据存在一定密度,可用差分代替微分来计算偏导数的近似值。则点云数据中任意一点P (Xi, Yj, Zij)切面元的法向量n=(n1, n2, 1)可近似表示为:

    (11)

    式中,(Xi-1, Yj, Zi-1, j)、(Xi+1, Yj, Zi+1, j)和(Xi, Yj-1, Zi, j-1)、(Xi, Yj+1, Zi, j+1)分别表示点P(Xi, Yj, Zi, j)在物空间坐标xy方向上相邻两点的坐标。

    将所求得的向量n(n1, n2, 1)归一化,得到的单位向量(n1/‖n‖, n2/‖n‖, 1/‖n‖)即为式(3)中的物方面元法向量(a, b, c),由式(1)得法向量的方向角为:

    (12)

    在点云数据优化过程中,每个三维点均采用以上方法计算其近似的法向量方向角(αi, βi), 代替(0°, 90°)作为法向量方向角的初值建立初始物方面元,这样计算得到的初始物方面元更接近该三维点真实的局部切平面,以此作为初值可以有效地减少迭代次数,降低因面元投影变形过大而造成的误匹配,提高影像匹配的效率和精度。算法具体流程如图 3

    图  3  基于非固定初始面元的无人机影像点云优化算法流程
    Figure  3.  Process of UAV Point Cloud Optimization Based on Non-fixed Initial Patch algorithm

    本文的主要改进思路如下。

    (1)最小二乘匹配及其相关算法的匹配精度可以达到子像素级别,但由于其计算的复杂性,这类算法多用于局部特征的定位与优化,且算法的精度和匹配效率依赖于自定义的初值。本文采用差分代替微分计算三维点局部正切平面的近似法向量,以其作为初始法向量建立的初始物方面元更接近实际真值,可以减少最小二乘匹配的迭代次数,降低因投影变形过大造成的误匹配,从而加快收敛速度,提高匹配精度。同时将基于物方面元的最小二乘匹配算法从卫星影像单点匹配推广到无人机点云优化,拓宽了物方面元最小二乘匹配的适用范围。

    (2)影像窗口尺寸μ也是人为定义的参数之一,文献[5, 18]中并没有对影像窗口尺寸的选择做详细的说明。在实际应用中,影像窗口尺寸过大会增加算法的复杂度,同时会降低影像窗口中心点的权重;而尺寸过小导致参与匹配的像素过少,可能无法真实反映影像的纹理,得不到理想的匹配结果。本文通过实验验证影像窗口大小对匹配精度的影响,讨论可以兼顾匹配效率与匹配精度的影像窗口尺寸选择。

    为了验证算法的可行性,实验从迭代次数和匹配精度两方面对实验结果进行分析。实验数据采用阳江地区与西北大学校区两组低空遥感影像数据,影像重叠度大于70%。影像的内、外方位元素已知。影像大小分别为8 206像素×6 078像素、3 888像素×2 592像素。实验采用立体匹配方法生成点云数据,点云数目分别为170 203,55 715。实验中当给定的辐射畸变改正值dh0i,dh1i,高程坐标的改正值Zc与法向量方向角的改正值dα、dβ均小于10-5时认为改正值收敛,迭代成功。而当迭代次数超过1 000次改正值仍没有收敛时认为优化失败并停止迭代。实验数据如图 4所示。

    图  4  实验数据
    Figure  4.  Experimental Data

    由于最小二乘匹配相关算法的匹配精度受影像方位参数精度的影响极大,两组数据无论在匹配成功率和匹配精度上均有较大不同,因此在匹配效率比较实验中,可将影像窗口设为一定大小(如设定为11像素×11像素),分别比较两组数据采用本文方法与基于物方面元的最小二乘匹配方法,当初始法向量方向角(a, b)分别为(0°, 90°)、(5°, 85°)、(25°, 65°)、(35°, 55°)、(45°, 45°)、(65°, 25°)、(85°, 5°)、(90°, 0°)时匹配优化结果的差异。具体见表 1

    表  1  影像数据匹配结果
    Table  1.  Matching Results of Image 1 and Image 2
    影像 面元初始角度 模型范围内概率/% 相关系数大于0.6/% 平均相关系数(优化前/优化后) 高程差ΔZ/m 平均迭代次数/次
    影像1 (0°, 90°) 83.95 67.51 0.70/0.77 1.792 36.46
    (5°, 85°) 84.27 67.88 0.71/0.78 1.482 35.74
    (25°, 65°) 84.30 67.99 0.71/0.78 1.470 35.71
    (35°, 55°) 84.50 68.50 0.70/0.78 1.462 35.26
    (45°, 45°) 84.77 68.95 0.70/0.78 1.446 35.05
    (65°, 25°) 85.55 70.52 0.69/0.79 1.424 34.53
    (85°, 5°) 75.31 54.84 0.51/0.74 1.763 37.63
    (90°, 0°) 72.77 50.31 0.44/0.70 1.941 37.03
    本文方法 88.61 72.11 0.76/0.86 1.360 25.90
    影像2 (0°, 90°) 95.45 81.58 0.76/0.83 0.517 27.92
    (5°, 85°) 95.47 75.08 0.76/0.83 0.526 27.68
    (25°, 65°) 95.62 81.86 0.75/0.83 0.512 27.22
    (35°, 55°) 95.65 81.88 0.75/0.83 0.506 27.04
    (45°, 45°) 95.65 82.01 0.74/0.83 0.506 27.08
    (65°, 25°) 95.70 82.06 0.70/0.83 0.493 26.95
    (85°, 5°) 85.52 58.21 0.44/0.74 0.779 33.91
    (90°, 0°) 85.96 59.60 0.43/0.74 0.770 33.21
    本文方法 98.92 93.05 0.82/0.87 0.426 21.85
    下载: 导出CSV 
    | 显示表格

    表 1中模型范围内概率表示匹配结束后计算出的三维点坐标在三维模型正确范围内的点占全部点数的百分比。从两组实验的匹配结果来看,基于物方面元的最小二乘匹配方法中,当初始物方面元近似垂直于水平面时,迭代的效率与匹配精度明显降低,这是因为近似垂直于水平面的面元在投影过程中会造成较大的变形,进而影响匹配精度。然而无论固定哪一个初始面元方向,匹配效率和精度均不如本文提出的方法。

    由于影像窗口间的相关系数在一定程度上会反映窗口中心点的相关程度,相关系数越接近1,表示影像窗口中心点是匹配点的可能性越大。所以我们选择面元投影到参考影像与搜索影像得到的影像窗口间相关系数大于一定阈值的点作为验证匹配精度的有效数据。在本文实验中,该阈值取值为0.6(同PMVS算法)。在无人机影像三维建模中,三维点的高程坐标误差是最大的,所以本文采用高程精度作为评价匹配精度的标准。精度验证的方法与文献[18]的方法相同,采用空间前方交会所得的高程坐标Z和匹配实验中所得的物方面元中心点高程坐标Zc作比较,以高程差ΔZ=| Zc-Z|作为验证算法精度的指标。从表 1~2中可以看出,在有效数据范围内,本文方法可以显著提高匹配效率和匹配精度。这是因为初始物方面元越接近地物的真实局部切平面,物方面元投影到影像上的影像窗口变形越小,影像窗口中的像素相关性越高,误差方程收敛所需迭代次数越少,匹配失败与误匹配的可能性也就越小。

    表  2  影像窗口尺寸对优化结果影响实验
    Table  2.  Experiments of Image 1 and Image 2 Window Size Effect on Optimization Results
    影像 影像窗口尺寸 模型范围内概率/% 相关系数大于0.6/% 平均相关系数 高程差ΔZ/m
    固定方向物方面元 本文方法 固定方向物方面元 本文方法 固定方向物方面元 本文方法 固定方向物方面元 本文方法
    影像1 7×7 72.23 71.91 49.36 49.46 0.627 0.630 1.812 1.697
    11×11 84.24 88.61 67.78 72.11 0.784 0.863 1.481 1.360
    15×15 88.75 89.11 74.62 74.83 0.807 0.837 1.363 1.261
    21×21 91.90 92.20 81.71 82.10 0.840 0.841 1.348 1.240
    31×31 93.74 94.29 87.32 87.92 0.870 0.870 1.375 1.273
    51×51 94.56 95.83 91.44 92.45 0.899 0.898 1.481 1.376
    影像2 7×7 87.16 87.47 66.74 67.00 0.719 0.752 0.810 0.705
    11×11 95.51 98.92 81.64 93.05 0.831 0.867 0.526 0.426
    15×15 97.64 98.92 87.75 93.41 0.835 0.855 0.458 0.344
    21×21 98.70 98.92 92.33 94.05 0.865 0.877 0.439 0.325
    31×31 99.32 99.47 95.81 96.37 0.892 0.898 0.441 0.328
    51×51 99.50 99.61 97.10 97.21 0.902 0.905 0.492 0.386
    下载: 导出CSV 
    | 显示表格

    最小二乘匹配及其相关算法中,影像窗口尺寸μ代表了参与匹配计算的像素点个数。μ越大,则参与匹配的像素越多,列出的误差方程(μ×μ个)也越多,计算量也越大。而影像窗口过小,参与匹配计算的像素过少,可能无法真实反映影像的纹理,会导致匹配精度的下降甚至匹配失败。实验分别比较影像窗口尺寸为7×7、11×11、15×15、21×21、31×31、51×51时的匹配精度,实验结果见表 2。从表 2的结果可以看到,随着影像窗口尺寸的逐渐变大,得到的有效数据越多,影像窗口间的相关系数越大。然而,当影像窗口的尺寸从21×21增加到31×31时,虽然相关系数有所增长但Z坐标的精度反而有所下降。这是因为参与匹配像素的增加会削弱影像畸变与重采样等因素的影响,并降低噪声点在匹配计算中的权重,提高影像窗口间的相关程度。而当参与匹配计算的像素过多时,影像窗口中心点的权重也会逐渐变小,反而降低了匹配的精确度。因此,从匹配获得的坐标精度考虑,在已选取的实验尺寸中,影像窗口像素个数为21×21时,得到的三维点Z坐标的精度最高,可以同时兼顾匹配精度与匹配效率。

    本文提出基于非固定初始面元的无人机影像点云优化算法, 利用密集点云数据在物方空间坐标系中的几何关系,采用差分代替微分计算地物表面局部正切平面的近似法向量,以此作为初始法向量建立初始物方面元,对面元为不同初值时的匹配效率和匹配精度进行实验验证,并通过实验讨论了匹配计算中影像窗口的尺寸选择。

    从匹配精度与匹配效率对比实验中可以看出,本文算法在匹配效率和三维重建精度上都优于基于物方面元的最小二乘匹配算法。实验还验证了,当影像窗口选取特定值时,本文方法可以同时兼顾优化精度与匹配效率。从最后的精度验证结果可以看出,无论影像窗口尺寸的大小如何选取,本文的算法计算出的三维点云精度均更高。

  • 图  1   PMVS中的局部正切平面

    Figure  1.   Patch in PMVS

    图  2   物方面元与光束交点

    Figure  2.   Intersection of the Object Patch and the Beam

    图  3   基于非固定初始面元的无人机影像点云优化算法流程

    Figure  3.   Process of UAV Point Cloud Optimization Based on Non-fixed Initial Patch algorithm

    图  4   实验数据

    Figure  4.   Experimental Data

    表  1   影像数据匹配结果

    Table  1   Matching Results of Image 1 and Image 2

    影像 面元初始角度 模型范围内概率/% 相关系数大于0.6/% 平均相关系数(优化前/优化后) 高程差ΔZ/m 平均迭代次数/次
    影像1 (0°, 90°) 83.95 67.51 0.70/0.77 1.792 36.46
    (5°, 85°) 84.27 67.88 0.71/0.78 1.482 35.74
    (25°, 65°) 84.30 67.99 0.71/0.78 1.470 35.71
    (35°, 55°) 84.50 68.50 0.70/0.78 1.462 35.26
    (45°, 45°) 84.77 68.95 0.70/0.78 1.446 35.05
    (65°, 25°) 85.55 70.52 0.69/0.79 1.424 34.53
    (85°, 5°) 75.31 54.84 0.51/0.74 1.763 37.63
    (90°, 0°) 72.77 50.31 0.44/0.70 1.941 37.03
    本文方法 88.61 72.11 0.76/0.86 1.360 25.90
    影像2 (0°, 90°) 95.45 81.58 0.76/0.83 0.517 27.92
    (5°, 85°) 95.47 75.08 0.76/0.83 0.526 27.68
    (25°, 65°) 95.62 81.86 0.75/0.83 0.512 27.22
    (35°, 55°) 95.65 81.88 0.75/0.83 0.506 27.04
    (45°, 45°) 95.65 82.01 0.74/0.83 0.506 27.08
    (65°, 25°) 95.70 82.06 0.70/0.83 0.493 26.95
    (85°, 5°) 85.52 58.21 0.44/0.74 0.779 33.91
    (90°, 0°) 85.96 59.60 0.43/0.74 0.770 33.21
    本文方法 98.92 93.05 0.82/0.87 0.426 21.85
    下载: 导出CSV

    表  2   影像窗口尺寸对优化结果影响实验

    Table  2   Experiments of Image 1 and Image 2 Window Size Effect on Optimization Results

    影像 影像窗口尺寸 模型范围内概率/% 相关系数大于0.6/% 平均相关系数 高程差ΔZ/m
    固定方向物方面元 本文方法 固定方向物方面元 本文方法 固定方向物方面元 本文方法 固定方向物方面元 本文方法
    影像1 7×7 72.23 71.91 49.36 49.46 0.627 0.630 1.812 1.697
    11×11 84.24 88.61 67.78 72.11 0.784 0.863 1.481 1.360
    15×15 88.75 89.11 74.62 74.83 0.807 0.837 1.363 1.261
    21×21 91.90 92.20 81.71 82.10 0.840 0.841 1.348 1.240
    31×31 93.74 94.29 87.32 87.92 0.870 0.870 1.375 1.273
    51×51 94.56 95.83 91.44 92.45 0.899 0.898 1.481 1.376
    影像2 7×7 87.16 87.47 66.74 67.00 0.719 0.752 0.810 0.705
    11×11 95.51 98.92 81.64 93.05 0.831 0.867 0.526 0.426
    15×15 97.64 98.92 87.75 93.41 0.835 0.855 0.458 0.344
    21×21 98.70 98.92 92.33 94.05 0.865 0.877 0.439 0.325
    31×31 99.32 99.47 95.81 96.37 0.892 0.898 0.441 0.328
    51×51 99.50 99.61 97.10 97.21 0.902 0.905 0.492 0.386
    下载: 导出CSV
  • [1] 林宗坚, 李德仁, 胥燕婴.对地观测技术最新进展评述[J].测绘科学, 2011, 36(4): 5-8 http://www.cnki.com.cn/Article/CJFDTOTAL-CHKD201104002.htm

    Lin Zongjian, Li Deren, Xu Yanying. General Review on the New Progress of Earth Observations[J]. Science of Surveying and Mapping, 2011, 36(4): 5-8 http://www.cnki.com.cn/Article/CJFDTOTAL-CHKD201104002.htm

    [2]

    Scharstein D, Szeliski R. A Taxonomy and Evaluation of Dense Two-frame Stereo Correspondence Algorithms[J]. International Journal of Computer Vision, 2002, 47(1-3): 7-42 http://cn.bing.com/academic/profile?id=2177253483&encoded=0&v=paper_preview&mkt=zh-cn

    [3]

    Liang T, Heipke C. Automatic Relative Orientation of Aerial Images[J]. Photogrammetric Engineering and Remote Sensing, 1996, 62(1): 47-55 http://cn.bing.com/academic/profile?id=2140022214&encoded=0&v=paper_preview&mkt=zh-cn

    [4] 林宗坚. UAV低空航测技术研究[J].测绘科学, 2011, 36(1): 5-9 http://www.cnki.com.cn/Article/CJFDTOTAL-CHKD201101002.htm

    Lin Zongjian. UAV Borne Low Altitude Photogrammetry System[J]. Science of Surveying and Mapping, 2011, 36(1): 5-9 http://www.cnki.com.cn/Article/CJFDTOTAL-CHKD201101002.htm

    [5] 江万寿.航空影像多视匹配与规则建筑物自动提取方法与研究[D].武汉:武汉大学, 2004 http://d.wanfangdata.com.cn/ExternalResource-lngcjsdxxb201006014%5e11.aspx

    Jiang Wanshou. Multiple Aerial Image Matching and Automatic Building Detection [D]. Wuhan: Wuhan University, 2004 http://d.wanfangdata.com.cn/ExternalResource-lngcjsdxxb201006014%5e11.aspx

    [6]

    Seitz S M, Curless B, Diebel J, et al. A Comparison and Evaluation of Multi-view Stereo Reconstruction Algorithms[C]. 2006 IEEE Computer Society Conference on Computer Vision and Pattern Recognition, New York, USA, 2006

    [7]

    Furukawa Y, Ponce J. Accurate, Dense, and Robust Multiview Stereopsis[J]. Pattern Analysis and Machine Intelligence, IEEE Transactions on, 2010, 32(8): 1 362-1 376 doi: 10.1109/TPAMI.2009.161

    [8]

    Nelder J A, Mead R. A Simplex Method for Function Minimization[J]. The Computer Journal, 1965, 7(4): 308-313 doi: 10.1093/comjnl/7.4.308

    [9] 林宗坚, 宣文玲, 王艳.高分辨率影像定位的一种新方法[J].武汉大学学报.信息科学版, 2004, 29(4): 363-366 http://www.cnki.com.cn/Article/CJFDTOTAL-WHCH200404019.htm

    Lin Zongjian, Xuan Wenling, Wang Yan. A Novel Method for High Resolution Image Positioning[J]. Geomatics and Information Science of Wuhan University, 2004, 29(4): 363-366 http://www.cnki.com.cn/Article/CJFDTOTAL-WHCH200404019.htm

    [10]

    Ackermann F. High Precision Digital Image Correlation[C]. 39th Photogrammetry Week, Stuttgart, Germany, 1983

    [11] 张祖勋, 张剑清, 胡翔云.基于物方空间几何约束最小二乘匹配的建筑物半自动提取方法[J].武汉大学学报信息科学版, 2001, 26(4): 290-295 http://www.cnki.com.cn/Article/CJFDTOTAL-WHCH200104001.htm

    Zhang Zuxun, Zhang Jianqing, Hu Xiangyun. Semiautomatic Building Extraction by Least Squares Matching Based on GeoMetrical Constraints in Object Space [J]. Geomatics and Information Science of Wuhan University, 2001, 26(4): 290-295 http://www.cnki.com.cn/Article/CJFDTOTAL-WHCH200104001.htm

    [12] 林祥国, 张继贤, 李海涛, 等.基于T型模板匹配半自动提取高分辨率遥感影像带状道路[J].武汉大学学报信息科学版, 2009, 34(3): 293-296 http://www.cnki.com.cn/Article/CJFDTOTAL-WHCH200903010.htm

    Lin Xiangguo, Zhang Jixian, Li Haitao, et al. Semi-automatic Extraction of Ribbon Road from High Resolution Remotely Sensed Imagery by a T-shaped Template Matching[J]. Geomatics and Information Science of Wuhan University, 2009, 34(3): 293-296 http://www.cnki.com.cn/Article/CJFDTOTAL-WHCH200903010.htm

    [13]

    Rosenholm D. Multi-point Matching Using the Least Squares Technique for Evaluation of Three-dimensional Models[J]. Photogrammetric Engineering and Remote Sensing, 1987, 53(6): 621-626 http://d.wanfangdata.com.cn/ExternalResource-Y1030332%5e41.aspx

    [14]

    Wu J, Chang C, Wang C. Least Squares Image Matching and Its Accuracy Improvements[J]. Journal of Photogrammetry and Remote Sensing, 2007, 12(3): 217-224

    [15]

    Gruen A, Akca D. Least Squares 3D Surface and Curve Matching[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2005, 59(3): 151-174 doi: 10.1016/j.isprsjprs.2005.02.006

    [16]

    Baltsavias E P. Multiphoto Geometrically Constrained Matching[D]. Zurich: Institute of Geodesy and Photogrammetry, ETH, 1991

    [17] 张过, 陈钽, 潘红播, 等.基于有理多项式系数模型的物方面元最小二乘匹配[J].测绘学报, 2011, 40(5): 592-597 http://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201105012.htm

    Zhang Guo, Chen Tan, Pan Hongbo, et al. Patch-based Least Squares Image Matching Based on Rational Polynomial Coefficients Model[J]. Acta Geodaetica et Cartographica Sinica, 2010, 39(3): 592-597 http://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201105012.htm

    [18] 张剑清, 潘励, 王树根.摄影测量学[M].武汉大学出版社, 2003

    Zhang Jianqing, Pan Li, Wang Shugen. The Photogrammetry[M]. Wuhan: Wuhan University Press, 2003

    [19] 同济大学数学系.高等数学[M]. 6版.北京:高等教育出版社, 2007

    Mathematics Department of Tongji University. Higher Mathematics[M]. 6th ed. Beijing: Higher Education Press, 2007

  • 期刊类型引用(5)

    1. 花奋奋,王萌萌,于洋,卢丽君,黄国满. 融合枝切法和最小费用流的多级相位解缠方法. 测绘科学. 2022(12): 138-146+183 . 百度学术
    2. 唐桂彬,周波. 基于遥感GIS的地震应急信息质量控制系统设计. 计算机测量与控制. 2021(10): 98-102 . 百度学术
    3. 景贵飞. 遥感可信度及可信系统互操作分析. 遥感信息. 2020(02): 1-7 . 百度学术
    4. 黄国满,程春泉,赵争,卢丽君,杨书成. 机载SAR遥感测图技术及应用. 测绘科学. 2019(06): 105-113 . 百度学术
    5. 金鼎坚,王建超,吴芳,高子弘,韩亚超,李奇. 航空遥感技术及其在地质调查中的应用. 国土资源遥感. 2019(04): 1-10 . 百度学术

    其他类型引用(1)

图(4)  /  表(2)
计量
  • 文章访问数:  1793
  • HTML全文浏览量:  59
  • PDF下载量:  456
  • 被引次数: 6
出版历程
  • 收稿日期:  2015-11-02
  • 发布日期:  2016-08-04

目录

/

返回文章
返回