文章快速检索  
  高级检索
基于双特征的丘陵山区耕地低空遥感图像配准算法
宋飞1,2,3, 杨扬1,2, 杨昆1,2, 张愫1,3, 毕东升1,3     
1. 云南师范大学 信息学院, 昆明 650500;
2. 西部资源环境地理信息技术教育部工程研究中心, 昆明 650500;
3. 云南师范大学 模式识别与人工智能实验室, 昆明 650500
摘要: 针对丘陵山区耕地小型无人机航拍图像(低空遥感图像)中的尺度变化、几何畸变、图像重叠等问题,提出了基于双特征的丘陵山区耕地低空遥感图像配准算法。该算法鉴于丘陵山区耕地背景环境复杂、光照因素等影响,采用尺度不变特征SURF算法提取了遥感图像的特征点,并构建了能够稳健描述航拍图像几何特征的双特征描述子;在此基础上,以高斯混合模型(GMM)为核心,结合2个单一特征差异描述子(基于欧氏距离的全局特征和基于和向量的局部特征)构造的双特征描述子,得到了能够同时通过2种特征进行对应关系评估的双特征有限混合模型(DFMM),并通过再生核希尔伯特空间(RKHS),基于高斯径向基函数(GRBF)对待配准图像进行了全局与局部结构双约束的空间变换更新。为了验证本文算法的可行性及其性能,采用小型无人机航拍的丘陵山区坡耕地多视角遥感图像开展了实验,将本文算法与SIFT、SURF、CPD、AGMReg、GLMDTPS及PRGLS进行了比较。实验结果表明,本文算法不仅在不同坡度的坡耕地航拍图像多视角配准过程中,均具有较好的鲁棒性,也适用于部分复杂地形小型无人机航拍的多视角遥感图像配准。
关键词: 图像配准     小型无人机     双特征     有限混合模型     再生核希尔伯特空间(RKHS)    
Low-altitude remote sensing image registration algorithm based on dual-feature for arable land in hills and mountains
SONG Fei1,2,3, YANG Yang1,2, YANG Kun1,2, ZHANG Su1,3, BI Dongsheng1,3     
1. School of Information Science and Technology, Yunnan Normal University, Kunming 650500, China;
2. The Engineering Research Center of GIS Technology in Western China of Ministry of Education of China, Kunming 650500, China;
3. Laboratory of Pattern Recognition and Artificial Intelligence, Yunnan Normal University, Kunming 650500, China
Received: 2017-10-30; Accepted: 2018-01-12; Published online: 2018-04-10
Foundation item: National Nature Science Foundation of China (41661080); Scientific Research Foundation of Yunnan Provincial Department of Education (2018Y037); Doctoral Scientific Research Foundation of Yunnan Normal University (01000205020503065)
Corresponding author. YANG Yang, E-mail: yyang_ynu@163.com
Abstract: Small unmanned aerial vehicle (UAV) aerial images (low-altitude remote sensing image) of arable land in hills and mountains are confronted with multiple challenges of image processing due to its scale change, geometric distortion and image overlap. To address the problems, the low-altitude remote sensing image registration algorithm based on dual-feature for arable land in hills and mountains was proposed. Due to the complex environment and the influence of light factors in hills and mountains, the feature points of remote sensing images are extracted by using the scale-invariant SURF algorithm. And then the dual-feature descriptor using geometrical structure of aerial images was constructed. On this basis, by taking Gaussian mixture model (GMM) as the core and combining with two single feature difference descriptors (i.e., global distance descriptor based on euclidean distance and local structure descriptor based on sum vectors), the dual-feature finite mixture model (DFMM) was obtained, which can simultaneously evaluate the correspondence between two features. With the reproducing kernel Hilbert space (RKHS), the spatial transformation of the global and local structure of the registration image was carried out based on the Gaussian radial basis function (GRBF). In order to verify the feasibility and performance of the proposed algorithm, experiments were carried out by using UAV images with different viewpoints taken from sloping arable land in hills and mountains. Experimental results show that comparing with SIFT, SURF, CPD, AGMReg, GLMDTPS and PRGLS, our method provides better performances in most cases, and can apply to multi-view remote sensing image registration of other complex terrain by small UAV.
Keywords: image registration     small UAV     dual-feature     finite mixed model     reproducing kernel Hilbert space (RKHS)    

丘陵山区是中国南方主要的地貌结构[1],其地面崎岖不平、岗拗交错,又常年受多雾、多雨气候条件的影响,导致耕地数量减少,水土流失加重,以及土壤环境污染加剧等现象。严重威胁着丘陵山区的农业生产、粮食安全和可持续发展[2-3]。在此严峻背景下,如何及时提取耕地分布信息,对于改善丘陵山区居民生活水平至关重要。

航天、航空遥感技术是目前提取耕地分布信息的重要方法[4]。但现有大、中型的航天、航空遥感技术存在气象影响因子多、重访周期过长、应急不及时和成本高等缺点。因此,采用小型无人机进行耕地分布信息的遥感获取,具有机动灵活、高时效、高分辨率、能耗低、成本低和受天气状况影响小等优点[5]。在小型无人机获取遥感图像(低空遥感图像)的过程中,由于受飞行器姿态、高度、速度等因素影响,经常会导致获取的遥感图像相对于地面目标实际位置发生挤压、扭曲、伸展、偏移以及图像重叠度不规则等一系列问题。消除这些因素所造成的同一场景(如不同视角的同一区域)不在同一坐标系(相同尺度的空间坐标系)的情形,是顺利实现航拍数据分析的前提和保证。

针对以上问题,遥感图像配准方法可以依据一定的相似性度量,对在不同条件下获取的两幅或多幅遥感图像进行对准、叠加的操作,使得图像在同一坐标系下获得最佳的重合。遥感图像配准方法主要分为2类:基于区域的配准方法和基于特征(例如点、线、面)的配准方法。由于点特征具有光照、旋转不变,可靠性高,且易提取性、计算简单、运算速度快,因此,选择图像中点特征研究的最多,典型的点特征提取算法有Harris算子[6]、SIFT(Scale-Invariant Feature Transform)算子[7]、SURF(Speeded Up Robust Features)算子[8]等。其中,Bay等[8]提出的SURF算子,不仅具有尺度、光照不变性,而且相比于SIFT算子计算量小,处理速度快而受到了领域内学者的青睐。

在基于特征点的图像配准方法研究领域,Brook和Bendor[9]通过在SURF算子中加入拓扑关系提高了遥感图像的配准精度;Zhao等[10]提出了基于SIFT和区域互信息优化的无人机遥感图像方案,提高了配准的精度;Lei等[11]通过在SIFT算法中加入随机抽样一致性和最小二乘法提高了无人机遥感图像自动配准的精度;乔川等[12]针对传统图像配准方法难以对海洋、沙漠和草原等特征不明显区域航空遥感图像进行配准的问题,提出了一种基于地理位置信息的航空遥感图像配准算法;Myronenko和Song[13]提出了可以同时处理刚性和非刚性配准问题的CPD(Coherent Point Drift)算法,可以在牺牲些许配准精度的前提下,通过采用快速高斯变换[14]和矩阵低秩逼近[15]技术减少计算量来提升算法的配准速度。Wang等[16]提出了一种基于非对称高斯模型的非刚性配准算法(AGMReg),在再生核希尔伯特空间(RKHS)决了正则化理论下的优化问题。Yang等[17]提出了一种基于欧氏距离的全局特征描述子和局部距离特征描述子的混合特征描述子的非刚性配准(GLMDTPS)算法,大幅提升了非刚性点阵配准在处理冗余点和噪音等方面的鲁棒性。Ma等[18]提出的PRGLS通过匹配形状上下文(Shape Context,SC)描述子来获取一个二进制对应矩阵,如果分配高斯混合模型(Gaussian Mixture Model,GMM)的成员概率接近1,则表示匹配,从而来改善CPD的配准性能。

本文在综述图像配准方法研究现状的基础上,针对丘陵山区耕地多视角航拍图像中的尺度变化、几何畸变和图像重叠等问题,开展了图像特征提取、匹配等相关算法和技术的研究。主要研究工作以及贡献:首先,针对耕地复杂的空间分布和几何形状,本文构建了能够稳健描述航拍图像几何特征的双特征描述子(Dual Feature Descriptor,DFD),此描述子由基于欧氏距离的全局特征和基于和向量的局部特征构成。在图像配准过程中,通过使用不确定性退火方案合理调节2个特征的强度。其次,采用特征点集之间欧氏距离度量的匹配策略最为常用,但匹配精度不高。为了提高低空遥感图像的匹配精度,本文以GMM为核心,结合构建的双特征描述子,得到了能够同时通过2种特征进行对应关系评估的双特征有限混合模型。最后,为了维护特征点集在进行空间变换更新时的局部与全局结构稳定,构建了基于L2E[19-21]的双约束(局部、全局约束项)能量方程来估计变换函数。

1 方法

首先介绍本文算法的3个贡献:构建双特征描述子、基于双特征有限混合模型、基于L2E的双约束能量方程。其次,详述具体的算法步骤。最后,进行算法分析。假定D维上的2个点集YM×D=[y1, y2, …, ym}TXN×D=[x1, x2, …, xn]T(待配准图像特征点集Y、参考图像特征点集X),其中ymxn分别表示为Y的第m个特征点和X的第n个特征点。

1.1 构建双特征描述子

首先,以特征点ym(注意:X处理类似)为中心提取出K个最近相邻点,把最近相邻点以偏离中心点ym距离升序构成局部片段(见图 1),并将特征点ym与局部片段中各点形成向量求和,得到ym的局部结构描述子(Local Structure Descriptor,LSD)表达式为

图 1 心型特征点集中某中心点与其最近相邻点(5个点)构成的局部片段 Fig. 1 A central point and its nearest neighboring points(five points) construct a small fragment of heart-shaped feature points
(1)

式中:γ为一个权重参数,用来调节μmkLm之间的平衡;μmk为向量。当ymk相对于ym是冗余点时,权重参数γ将起到消除这些冗余点的作用,LmY具体的计算方阵形式为

(2)

式中:IM×M对角矩阵;RY为一个M×M的矩阵;定义rij表示RY中任意一项,yiyj表示Y中任意一点,且

(3)

然后,基于欧氏距离的全局特征GY和基于和向量的局部特征LmY构建具有特征互补的双特征描述子,其表达式为

(4)

式中:DFDY为基于GYLY的代价矩阵(Cost Matrix),通过调节权重参数ε,控制GYLY之间的代价平衡。在迭代开始时,ε值较大,最小化DFDY相当于最小化LY的代价;随着迭代次数增加,利用非确定性退火方案(其中t为当前迭代次数,r为常数)使ε逐渐减小,该过程倾向于最小化GY的代价。

1.2 基于双特征有限混合模型

Myronenko和Song[13]在CPD方法中提出2个特征点集对应关系评估可以视为概率密度估计问题。本节针对特征点集配准问题采用GMM进行建模,结合双特征描述子对混合模型的参数进行赋值,从而改进对应关系评估的效果。双特征有限混合模型的推导过程,具体如下:

1) 构建GMM。将整个参考图像特征点集X拟合成一个GMM,使其高斯模型(Gaussian Model,GM)的中心与待配准图像特征点集Y保持一致,故GMM的概率密度函数(Probability Density Function,PDF)为

(5)

式中:P(m)=1/M为混合系数,描述了第m个高斯模型在混合模型中的权重。

典型地,对于点集配准问题,可以对GMM的所有组成部分采用同一的各向同性的协方差σ2I,且异常点的分布可以假设为均匀分布。令为未知参数集,其中μm为第m个高斯模型ym。于是,混合模型表达式为

(6)

式中:ω为异常点的比率,且大小0≤ω≤1,。然后,基于贝叶斯准则来计算后验概率pmn,即

(7)

式中:dLmn=‖xn-T(ym)‖, (dLmn)2xnT(ym)之间的欧氏距离,T(ym)为变换函数;pmn为待配准图像ym与参考图像xn在当前估计出的T(ym)下的吻合程度; dLjndLmn类似。

2) 构建双特征有限混合模型(DFMM)。若使用式(7)进行对应关系评估,只考虑点集X到点集Y之间的欧氏距离,可能导致配准过程的鲁棒性不好,例如:xaxb的局部结构不同,但xaxbym的欧氏距离一致,导致pmapmb相等。

为了解决以上问题,将双特征描述子式(4)代替方程式(7)中的欧氏距离,重写式(7)的形式得到DFMM,并将其定义如下:

(8)

式中:dDFDmn=dGmn+εdLmn, dDFDmndGmndLmn分别为混合特征描述子、欧氏距离和局部结构描述子;dDFDjndDFDmn类似。

1.3 基于L2E的双约束能量方程

从全局(欧氏距离)的角度结合局部特征(和向量)对点集空间分布和结构进行分析,为点集匹配问题提供了一条新颖而具有潜力的解决途径。但该方法不可避免的会产生误匹配,即异常点。因此,在空间变换更新中,不得不对变换函数T(ym)采取鲁棒的估计。本文构建基于L2E[19-21]的双约束能量方程来估计变换函数T(ym),其表达式为

(9)

式中:L2E为一个最小密度之间L2距离的鲁棒估计子,其特点适用于分析包含冗余点的大规模数据集。ℜ(T)表示基于正则化项的双约束项,对变换函数T(ym)进行平滑性约束。

1) 构建L2E估计子。假设有参数模型P(z|θ),本文目标是对L2距离进行最小化,即

(10)

式中:参数θ的真实值θ0未知,省略了与θ无关的常数项∫p(z|θ0)2dz。基于L2距离,参数θ可以由L2E估计子对模型p(z|θ)通过最小化准则式(11)来估计:

(11)

2) 构建ℜ(T)双约束项。在非刚性图像配准中,由于采用图像特征点集插值的变换函数T(ym)不是唯一的,因此,需要引入基于运动一致性的几何约束以解决这种不确定性问题,运动一致性理论在文献[22]中得到了直观的解释。在运动理论中,有一些重要的现象涉及到一致性。其中,速度一致性是对转换进行平滑的一种特殊的方式。本文使用基于变分原理的Tikhonov正则化框架[23]定义了保持运动一致性的双约束项:全局结构约束项和局部结构约束项,具体描述如下:

全局结构约束项,具体能量方程形式为

(12)

式中:λ>0为控制全局结构约束项的强度;T(ym)为变换函数,相当于在ym原始位置加入位移函数。正则化项‖T(ym)‖G2通过内积来定义‖T(ym)‖G2=〈T(ym), T(ym)〉GG为高斯径向基函数(Gaussian Radial Basis Function,GRBF),方程形式如下:

(13)

式中:β为一个常数,用来控制空间平滑性。

局部结构约束项,核心思想是通过判断某图像在空间变换前后的局部相似性, 来约束其在进行空间变换时的局部形变,方程形式如下:

(14)

式中:η>0为控制局部结构约束项的强度;Q= LdiagXc-LdiagT(ym)XcX通过对应关系评估得到的模糊对应关系矩阵。

1.4 图像配准过程

针对以上研究,本文提出了基于双特征的丘陵山区耕地低空遥感图像配准的算法——DFMM算法,其流程图如图 2所示。

图 2 图像配准流程图 Fig. 2 Flowchart of image registration

1.4.1 基于SURF的特征点提取

鉴于丘陵山区耕地背景环境复杂、光照因素等影响,本文使用尺度不变特征SURF算法分别提取待配准遥感图像Is和参考图像Ir的特征点集YM×D=[y1, y2, …, ym]TXN×D=[x1, x2, …, xn]TD为特征点集的维度,M=N。具体步骤如下:

1) 对参考图像与待配准图像分别进行SURF特征点检测与定位和特征描述子计算。

2) SURF特征点匹配。本文采用马氏距离进行相似性度量,其定义如下:

(15)

式中:Xi*为参考图像的第i个特征向量;Yj*为待配准图像的第j个特征向量;Σ-1为协方差矩阵;d为2个向量之间的马氏距离。马氏距离越小表示2个特征点越相似,记录马氏距离最小的特征点对。

1.4.2 特征点匹配

1) 对应关系评估

通过式(8)得到一个M×N矩阵P,便是待配准图像特征点集Y和参考图像特征点集X之间的一组模糊对应关系。然后由矩阵P得到对应的目标点阵Xc,其方程为

(16)

2) 空间变换更新

在RKHS[24]中对变换函数T(ym)进行建模。首先,选择一个正定核G(y, ym)(即高斯核)来定义一个RKHS H。然后,采用GRBF,得到变换函数T(ym)最优解表达式为

(17)

式中:系数cm为一个D×1维的待求解向量。于是,在一个无限维H空间中最小化问题归结于求解一个有限的M系数向量cm。然后,采用能量方程式(9)评估变换函数T(ym),具体形式为

(18)

式中:, Um, ·为矩阵U的第m行,ci为系数矩阵Ch×D的第i行。EGlobalELTD分别为全局结构约束项、局部结构约束项。然后,对方程式(14)与系数矩阵C取偏导数并获得如式(19)形式的梯度函数:

(19)

式中:H=exp{d[(UCPX)(UCPX)T]/(2σ2)}为M×1的向量,d(·)为矩阵的对角线;1D所有1的1×D行向量;⊙为Hadamard乘积;⊗为Kronecker卷积;A=[2KRx-2(RY)TRX+2K·(RY)T-2K2I];B=[2(RY)TRY-2KRY-2K·(RY)T+2K2I];R类似于二进制矩阵,只需将“1”的每个条目替换为高斯模糊的向量范数。

1.4.3 图像转换

本节采用“反推法(the backward approach)[25]”来建立TPS(Thin-Plate Spline)变换模型,以避免输出图像发生孔和(或)重叠的现象(由于离散化和四舍五入所引起的配准错误)。从待配准图像得来的转换图像是由目标像素坐标决定的(即与参考图像有着同样的坐标元),并且估计的TPS变换模型参考图像到浮动图像是“反推法”的关键理念,可以将其定义为

(20)

式中:ItIs分别为转换图像和待配准图像。It与参考图像Is具有相同的大小,且T(x, y)为从参考图像到待配准图像的TPS变换估计模型。T(x, y)的方程式为

(21)

式中:YxωYyω为特征点集Yω的坐标,内核函数U(Yiω, Yjω)=‖YiωYjω 2‖lnYiωYjω‖。因此,TPS变化估计模型(ω1, ω2, …, ωn, c1, cx, cy),可由方程式(22)求解得到:

(22)

式中:0为3×3的零矩阵;Yn×2的矩阵;P为一个n×3的矩阵,其中第iP为(1, Yωix, Yωiy);K为一个n×n核矩阵,其中Kij=U(Yiω, Yjω)。在得到TPS变化估计模型(ω1, ω2, …, ωn, c1, cx, cy)后即可得到It。其中,本部分对待配准图像Is的规格网格使用双三次插值(bicubic interpolation)算法。并通过使用反推法避免输出图像发生孔和重叠的现象。

1.5 算法分析

1.5.1 SURF提取丘陵山区耕地低空遥感图像特征点的优势

小型无人机在航拍丘陵山区耕地时,由于受地面崎岖不平和航拍视角等影响,会加剧航拍图像间产生非刚性几何畸变,增加图像配准的复杂度。因此,有效地选择图像特征点集(优于随机选取),会减小配准的计算复杂度。为了精确、鲁棒、快速地获得特征点集,本文开展了提取丘陵山区耕地低空遥感图像特征点集的研究(见图 3),通过实验比较,SUSAN算子提取的轮廓信息会增加匹配的计算复杂度; Harris算子不具有尺度不变性,而且在角点定位方面存在偏差,可能导致匹配不准确;SURF算子对平移、旋转、尺度变化和噪声等具有良好的不变性,对视觉变化、仿射变化也保持一定程度的稳定性,因而能够均匀提取到丘陵山区耕地低空遥感图像的特征点。

图 3 提取图像特征 Fig. 3 Image feature extraction

1.5.2 特征点匹配性能评估

1) 单一特征与双特征性能比较。

将只使用单一特征(欧氏距离)的CPD算法和使用双特征的DFMM算法的特征点匹配过程作对比(见图 4),在迭代开始时(如:Iter=1), CPD方法最初评估所得的对应点集整体位于质心区域,而使用DFMM算法最初评估,其分布更接近真正的目标点集。随着迭代次数增加(如:Iter=5, Iter=10),CPD算法评估得到的点集依旧整团分布于区域质心周围,相比之下,DFMM算法迭代更新得到的源点集位置更优。

图 4 单一特征与混合特征(双特征)性能比较 Fig. 4 Performance comparison of single feature and mixed feature (dual-feature)

2) L2E的鲁棒性。

参数估计一般采用最大似然估计(Maximum Likelihood Estimation,MLE)来完成。然而,当采用的参考图像与待配准图像近似得不够好,或者提取的特征点集中存在大量冗余点时, MLE将会产生严重的偏差。L2E是一个用于最小化密度之间L2距离的鲁棒估计子,其与相比,更适用于分析包含冗余点的大规模数据集。将L2E和最大似然估计MLE进行对比,评估L2E鲁棒性。

假设一维高斯模型p(z|θ)=Φ(z|u, 1),其中:,本文通过估计正态分布N(0, 1)在含冗余点的正态分布N(5, 1)干扰下所产生的均值,对比L2E与MLE的鲁棒性。例如:定义一个大小为800的数据样本表示内点,其服从正态分布N(0, 1),而5组数量递增的冗余点n=40, 80, 160, 240, 400由正态分布N(5, 1)生成。L2E与MLE的鲁棒性对比如图 5所示, L2E在0.5附近有一个全局最小值,在4附近存在一个局部最小值,并且随着冗余点的增加极值会变得更深。与此相反,MLE的偏差随着冗余点的增加而变大,效果的稳定性变差。

图 5 L2E与MLE随着冗余点数量的改变的鲁棒性对比 Fig. 5 Variation of L2E and MLE with number of redundant point and their robustness comparison

1.5.3 参数设置

1) 退火参数。通过实验分析,式(4)在开始迭代时ε=exp(0)=1随着迭代次数的增加,ε逐渐减小,最后ε=exp(12.5)≈3×10-6。实验中,考虑到效率问题,将式(1)初始化γ参数设置为0.9。

2) 平滑性控制。在GRBF中,β用来控制空间平滑性。将源图像的特征点的空间坐标归一化为[-1.5,1.5],因此β设置为2。

3) 协方差。在式(7)和式(18)中,σ2对非凸性问题处理上作用明显,退火方程σ2=ρσ2,其中ρ=0.9,σ2分别初始化为1和0.05。

4) 正则化参数。式(18)中有2个正则化参数λη,分别代表了全局和局部结构约束的权重大小。在本算法中,将λη分别设置为2。

5) 最近邻近点。在式(1)中,K的设定基于区别局部结构描述子所需的最少的点数量[17]。例如,区别角(包含2个相邻点)和十字(包含4个相邻点),至少需要4个相邻点。通过实验分析,在二维和三维配准中K的默认值设为5。

2 实验

首先,使用丘陵山区耕地的小型无人机航拍多视角遥感图像(低空遥感图像),验证本文提出的基于双特征的丘陵山区耕地低空遥感图像配准算法——DFMM算法;然后,验证本文算法在其他复杂地形的配准效果。

2.1 丘陵山区耕地低空遥感图像配准实验

2.1.1 研究区域及数据

坡耕地是丘陵山区的典型代表[26-28],是分布在山坡上6°~25°之间的地貌类型,被开垦后形成的耕地。根据坡耕地的坡度(地表单元陡缓的程度),本文实验数据集主要划分成2大类:(i)选取60张(30组)坡度为6°~15°(包括15°),不同角度变化的坡耕地小型无人机多视角遥感图像;(ii)选取60张(30组)坡度为15°~25°,不同角度变化的坡耕地小型无人机多视角遥感图像。

本实验所有小型无人机多视角遥感图像来自720云平台,区域分别为云南、贵州、四川、广西、甘肃及湖南丘陵山区的典型坡耕地。每组图像的航拍视角变化为俯角变化(30°~90°)以及水平视角变化(-90°~90°)。其大小范围为640 mm×450 mm ~1 100 mm×850 mm,实验数据集如表 1所示。

表 1 实验数据 Table 1 Experimental data
参数 (i) (ii)
坡度/(°) 6~15 15~25
实验组数 30 30
尺寸/(mm×mm) 640×450~
1 100×850
640×450~
1 100×850
俯角变化/(°) 30~90 30~90
水平视角/(°) -90~90 -90~90
特征点数 282~712 140~600

2.1.2 评估标准

在评估检验时,首先,手动的在待配准图像和参考图像上确定了至少20个地面真值。为了减小评估过程中产生的误差,所有地面真值选择的都是合理分布于容易识别的区域。然后,使用均方根误差(Root of Mean Square Error,RMSE)、平均绝对误差(Mean Absolute Error,MAE)来量化配准精度。统计中的相关配方和定义如下:

(23)

均方根误差能够很好地反映出选定地标与实际位置对应点之间的精密度,其中m′为所选地标的总数,yti为对应于xti的地标。

(24)

平均绝对误差通过计算选定地标与实际位置对应点之间的平均绝对差异,来检验算法的效果。

2.1.3 实验结果及分析

为了证明本文算法的可行性及其性能,将DFMM与SIFT[7]、SURF[8]、CPD[13]、AGMReg[16]、GLMDTPS[17]及PRGLS[18]算法进行了比较。结果表明,本文算法在进行坡耕地的小型无人机多视角遥感图像(低空遥感图像)配准过程中,具有较好的鲁棒性。在图 6图 7中分别列举出了5组坡度为6°~15°、15°~25°坡耕地的小型无人机多视角遥感图像(低空遥感图像)的匹配结果。图 6图 7前2行表示两幅取自不同视角的无人机航拍图像(待配准图像和参考图像)。每种配准方法用2行进行展示,首行表示特征点的配准结果,第2行中表示给定一个5×5的棋盘格,交替显示转换图像和参考图像。数据集(i)和(ii)的RMSE、MAE评估结果如表 2所示。

图 6 5组数据集(i)的图像配准结果示例 Fig. 6 Examples of image registration results for five sets of data sets (i)
图 7 5组数据集(ii)的图像配准结果示例 Fig. 7 Examples of image registration results for five sets of data sets (ii)
表 2 使用RMSE和MAE进行数据集(i)和(ii)实验的评估比较结果 Table 2 Experimental result of exprimental database (i) and (ii) evaluation and comparison using RMSE and MAE
误差 数据集 SURF SIFT CPD AGMReg GLMDTPS PRGLS DFMM
RMSE (i) 7.283 7 12.528 7 5.436 5 5.112 3 4.919 1 3.716 3 1.321 1
(ii) 6.062 7 11.556 6 3.257 6 3.190 4 3.106 5 2.483 2 1.012 7
MAE (i) 4.234 4 7.288 9 3.110 2 3.044 0 3.058 1 2.038 0 1.074 9
(ii) 5.441 1 6.008 0 2.239 4 2.170 2 2.103 8 1.504 5 0.645 2

结果分析:首先,分析以上5种算法的主要特征,SURF、SIFT特征点集的提取是基于强度信息,或者更确切地说,是尺度空间极值;CPD是采用欧氏距离,利用GMM建立对应关系;GLMDTPS主要利用基于欧氏距离的全局特征描述子和局部距离特征描述子的双特征描述子建立对应关系;AGMReg是以非对称高斯模型为核心,使用软分配技术来恢复对应关系;PRGLS基于全局特征和基于SC描述子来获取一个二进制对应矩阵,来完成图像相关匹配。DFMM主要以GMM为核心,结合2个单一特征差异描述子(基于欧氏距离的全局特征和基于和向量的局部特征)构建的双特征描述子(式(4)), 得到了改进后的GMM,使其能够同时通过2种特征进行对应关系评估(式(8))。

总结这5种算法中使用的主要特征包括:①强度信息差异;②全局几何结构差异;③局部几何结构差异;④全局几何约束;⑤局部几何约束。从表 2中的实验结果容易看到CPD、AGMReg、及PRGLS配准结果优于SURF、SIFT,其原因在于SURF、SIFT只采用特征①提取对应的特征点完成图像配准,而CPD、AGMReg、GLMDTPS及PRGLS分别还采用了特征②④、②④、②③及②③④,能够在一定程度上提高图像配准的精度。本文DFMM算法配准结果优于其他配准结果,其主要原因是采用了特征①②③④⑤,不仅能够提取对图像尺度和旋转变化具有不变性、对光照变化和图像变形具有较强适应性的特征点,而且还能有效地利用了图像特征点的局部信息和全局信息完成图像之间的配准。

2.2 其他复杂地形低空遥感图像配准实验

复杂地形地质结构特殊,时常发生滑坡、山体崩塌等自然地质灾害, 因此实时监测复杂地形非常必要。由于小型无人机具有机动灵活、高时效、高分辨率、能耗低、成本低、受天气状况影响小等优点,可以广泛地运用于复杂地形的监测。本节通过实验,验证了DFMM算法在不同复杂地形小型无人机航拍的多视角遥感图像(低空遥感图像)中的配准效果,所有数据来自720云平台。实验结果表明(表 3),DFMM算法具有较好的配准效果。在图 8分别列举出了4组不同复杂地形的对应关系(图 8(a))、匹配结果(图 8(b)),图 8(a)第1行表示不同视角的无人机航拍图像(待配准图像和参考图像),第2行表示匹配的对应关系;图 8(b)第1行表示特征点的配准结果,第2行表示给定一个5×5的棋盘格,交替显示转换图像和参考图像。

图 8 4组图像配准结果示例 Fig. 8 Examples of four sets of image registration results
表 3 使用RMSE和MAE进行复杂地形低空遥感图像实验的评估实验结果 Table 3 Experimental results of low-altitude remote sensing images with complex terrain evaluation using RMSE and MAE
误差 云南昆明老街 长城 大坝 铁路 均值
RMSE 1.874 1 1.231 8 1.042 3 1.453 6 1.400 5
MAE 1.321 6 1.875 7 1.031 0 1.234 2 1.365 6

3 结论

采用小型无人机提取丘陵山区耕地分布信息,具有机动灵活、高时效、高分辨率、能耗低、成本低、受天气状况影响小等优点,但地面的起伏变化和航拍视角变化会加剧图像间产生非刚性几何畸变及重叠度不规则等问题。针对以上问题,本文提出了一种基于双特征的丘陵山区耕地低空遥感图像配准算法——DFMM算法。

实验验证表明,该算法可以有效用于解决丘陵山区耕地小型无人机低空遥感图像间存在的非刚性几何畸变及重叠度不规则等问题。DFMM算法不仅能够稳健描述航拍图像几何特征,而且配准结果的鲁棒性和精度有明显提高。这种基于双特征的低空遥感配准算法,也适用于部分复杂地形小型无人机航拍的多视角遥感图像配准(低空遥感图像),将会进一步研究能否通用到所有复杂地形。

参考文献
[1] 易中懿, 曹光乔, 张宗毅. 我国南方丘陵山区农业机械化宏观影响因素分析[J]. 农机化研究, 2010, 32 (8): 229–233.
YI Z Y, CAO G Q, ZHANG Z Y. Analysis of macro factors influencing agricultural mechanization in hilly and mountainous areas of China[J]. Journal of Agricultural Mechanization, 2010, 32 (8): 229–233. DOI:10.3969/j.issn.1003-188X.2010.08.059 (in Chinese)
[2] 祝功武. 南方丘陵山区耕地现状, 潜力与开发对策:以德庆县为例[J]. 地理科学, 1998, 18 (1): 10–15.
ZHU G W. Status quo hilly area south of arable land, potential and development strategies:A case study in Deqing county[J]. Geographical Sciences, 1998, 18 (1): 10–15. (in Chinese)
[3] 吴兆娟, 魏朝富, 丁声源. 丘陵山区地块尺度耕地社会保障功能价值研究[J]. 资源科学, 2013, 35 (4): 95–103.
WU Z J, WEI Z F, DING S Y. Study on the social security function value of cultivated land in the hilly area of the mountainous region[J]. Resources Science, 2013, 35 (4): 95–103. (in Chinese)
[4] 杨邦杰, 裴志远, 周清波, 等. 我国农情遥感监测关键技术研究进展[J]. 农业工程学, 2002, 18 (3): 191–194.
YANG B J, PEI Z Y, ZHOU Q B, et al. Study on key technology of remote sensing monitoring in China[J]. Journal of Agricultural Engineering, 2002, 18 (3): 191–194. DOI:10.3321/j.issn:1002-6819.2002.03.045 (in Chinese)
[5] 汪沛, 罗锡文, 周志艳, 等. 基于微小型无人机的遥感信息获取关键技术综述[J]. 农业工程学报, 2014, 30 (18): 1–12.
WANG P, LUO X W, ZHOU Z Y, et al. Study on key technology of remote sensing information acquisition based on micro-UAV[J]. Journal of Agricultural Engineering, 2014, 30 (18): 1–12. DOI:10.3969/j.issn.1002-6819.2014.18.001 (in Chinese)
[6] HARRIS C. A combined corner and edge detector[J]. Alvey Vision Conference, 1988 (3): 147–151.
[7] LOWE D G. Distinctive image features from scale-invariant keypoints[J]. International Journal of Computer Vision, 2004, 60 (2): 91–110. DOI:10.1023/B:VISI.0000029664.99615.94
[8] BAY H, ESS A, TUYTELAARS T, et al. Speeded-up robust features[J]. Computer Vision and Image Understanding, 2008, 110 (3): 404–417.
[9] BROOK A, BENDOR E. Automatic registration of airborne and spaceborne images by topology map matching with SURF processor algorithm[J]. Remote Sensing, 2011, 3 (1): 65–82. DOI:10.3390/rs3010065
[10] ZHAO L Y, LYU B Y, LI X R, et al. Multi-source remote sensing image registration based on SIFT and optimization of regional mutual information[J]. Acta Physica Sinica, 2015, 64 (12): 124204.
[11] LEI T, LI L, KAN G, et al.Automatic registration of unmanned aerial vehicle remote sensing images based on an improved SIFT algorithm[C]//Eighth International Conference on Digital Image Processing.Chengdu: ICDIP, 2016: 100333R.
[12] 乔川, 丁亚林, 许永森, 等. 基于地理位置信息的图像配准方法及精度分析[J]. 光学学报, 2017, 37 (8): 318–328.
QIAO C, DING Y L, XU Y S, et al. Image registration method and precision analysis based on geographic location information[J]. Acta Optica Sinica, 2017, 37 (8): 318–328. (in Chinese)
[13] MYRONENKO A, SONG X. Point set registration:Coherent point drift[J]. IEEE Transactions on Pattern Analysis and Machine Intelligence, 2010, 32 (12): 2262–2275. DOI:10.1109/TPAMI.2010.46
[14] GREENGARD L, STRAIN J. The fast Gauss transform[J]. SIAM Journal of Scientific and Statistical Computing, 1991, 12 (1): 79–94. DOI:10.1137/0912004
[15] MARKOVSKY I. Structured low-rank approximation and its applications[J]. Automatica, 2008, 44 (4): 891–909. DOI:10.1016/j.automatica.2007.09.011
[16] WANG G, WANG Z, CHEN Y, et al. Robust point matching method for multimodal retinal image registration[J]. Biomedical Signal Processing and Control, 2015, 19 : 68–76. DOI:10.1016/j.bspc.2015.03.004
[17] YANG Y, ONG S H, FOONG K W C. A robust global and local mixture distance based non-rigid point set registration[J]. Pattern Recognition, 2015, 48 (1): 156–173. DOI:10.1016/j.patcog.2014.06.017
[18] MA J, ZHAO J, YUILLE A L. Non-rigid point set registration by preserving global and local structures[J]. IEEE Transactions on Image Processing, 2016, 25 (1): 53–64. DOI:10.1109/TIP.2015.2467217
[19] SCOTT D W. Parametric statistical modeling by minimum integrated square error[J]. Technometrics, 2001, 43 (3): 274–285. DOI:10.1198/004017001316975880
[20] SEVERINI T A. Robust and efficient estimation by minimising a density power divergence[J]. Biometrika, 1998, 85 (3): 549–559. DOI:10.1093/biomet/85.3.549
[21] MA J Y, ZHAO J W, TIAN J, et al.Robust estimation of nonrigid transformation for point set registration[C]//IEEE Conference on Computer Vision and Recognition. Piscataway, NJ: IEEE Press, 2013: 2147-2154.
[22] YUILLE A L, GRZYWACZ N M. A Mathematical analysis of the motion coherence theory[J]. International Journal of Computer Vision, 1989, 3 (2): 155–175. DOI:10.1007/BF00126430
[23] GROETSCH C W. The theory of Tikhonov regularization for fredholm equations of the first kind[M]. Boston: Pitman Advanced Publishing Program, 1984: 98-100.
[24] ARONSZAJN N. Theory of reproducing Kernels[J]. Transactions of the American Mathematical Society, 1950, 68 (3): 337–404. DOI:10.1090/S0002-9947-1950-0051437-7
[25] ZITOVA B, FLUSSER J. Image registration methods:A survey[J]. Image and Vision Computing, 2003, 21 (11): 977–1000. DOI:10.1016/S0262-8856(03)00137-9
[26] 范玉芳, 罗友进, 魏朝富. 西南丘陵山区坡耕地水平梯田工程设计分析[J]. 山地学报, 2010, 28 (5): 560–565.
FAN Y F, LUO Y J, WEI C F. Design and analysis of horizontal terraces of sloping land in hilly and hilly areas of southwest China[J]. Acta Terrestrial Journal, 2010, 28 (5): 560–565. DOI:10.3969/j.issn.1008-2786.2010.05.008 (in Chinese)
[27] 马芊红, 张光辉, 耿韧, 等. 我国水蚀区坡耕地土壤重金属空间分布及其污染评价[J]. 水土保持研究, 2017, 24 (2): 112–118.
MA Q H, ZHANG G H, GENG R, et al. Spatial distribution and pollution assessment of heavy metals in slope farmland of water erosion area in China[J]. Journal of Soil and Water Conservation, 2017, 24 (2): 112–118. (in Chinese)
[28] 陈伟, 朱继平, 陈小兵, 等. 我国坡耕地农业机械化发展现状[J]. 农机化研究, 2017, 39 (5): 1–5.
ZHU W, ZHU J P, CHEN X B, et al. China's sloping farmland agricultural mechanization development status[J]. Agricultural Mechanization Research, 2017, 39 (5): 1–5. DOI:10.3969/j.issn.1003-188X.2017.05.001 (in Chinese)
http://dx.doi.org/10.13700/j.bh.1001-5965.2017.0674
北京航空航天大学主办。
0

文章信息

宋飞, 杨扬, 杨昆, 张愫, 毕东升
SONG Fei, YANG Yang, YANG Kun, ZHANG Su, BI Dongsheng
基于双特征的丘陵山区耕地低空遥感图像配准算法
Low-altitude remote sensing image registration algorithm based on dual-feature for arable land in hills and mountains
北京航空航天大学学报, 2018, 44(9): 1952-1963
Journal of Beijing University of Aeronautics and Astronsutics, 2018, 44(9): 1952-1963
http://dx.doi.org/10.13700/j.bh.1001-5965.2017.0674

文章历史

收稿日期: 2017-10-30
录用日期: 2018-01-12
网络出版时间: 2018-04-10

相关文章

工作空间