自动绘制室内平面图的点云向量追踪算法

蔡来良, 宋德云, 胡青峰, 魏峰远, 舒前进

蔡来良, 宋德云, 胡青峰, 魏峰远, 舒前进. 自动绘制室内平面图的点云向量追踪算法[J]. 武汉大学学报 ( 信息科学版), 2021, 46(9): 1404-1411. DOI: 10.13203/j.whugis20190258
引用本文: 蔡来良, 宋德云, 胡青峰, 魏峰远, 舒前进. 自动绘制室内平面图的点云向量追踪算法[J]. 武汉大学学报 ( 信息科学版), 2021, 46(9): 1404-1411. DOI: 10.13203/j.whugis20190258
CAI Lai-liang, SONG De-yun, HU Qing-feng, WEI Feng-yuan, SHU Qian-jin. A Point Cloud Vector Tracing Algorithm for Automatic Drawing of Interior Plan[J]. Geomatics and Information Science of Wuhan University, 2021, 46(9): 1404-1411. DOI: 10.13203/j.whugis20190258
Citation: CAI Lai-liang, SONG De-yun, HU Qing-feng, WEI Feng-yuan, SHU Qian-jin. A Point Cloud Vector Tracing Algorithm for Automatic Drawing of Interior Plan[J]. Geomatics and Information Science of Wuhan University, 2021, 46(9): 1404-1411. DOI: 10.13203/j.whugis20190258

自动绘制室内平面图的点云向量追踪算法

基金项目: 

国家自然科学基金 41701597

国家自然科学基金 U1810203

中国博士后科学基金 2018M642746

详细信息
    作者简介:

    蔡来良, 博士, 副教授, 主要从事开采沉陷、点云处理方面的教学研究工作。cll@hpu.edu.cn

    通讯作者:

    宋德云,硕士。643823940@qq.com

  • 中图分类号: P237

A Point Cloud Vector Tracing Algorithm for Automatic Drawing of Interior Plan

Funds: 

National Natural Science Foundation of China 41701597

National Natural Science Foundation of China U1810203

the Postdoctoral Science Foundation of China 2018M642746

More Information
    Author Bio:

    CAI Lai-liang: CAI Lailiang, PhD, associate professor, specializes in mining subsidence and point cloud processing. E-mail: cll@hpu.edu.cn

    Corresponding author:

    SONG De-yun: SONG Deyun, master. E-mail: 643823940@qq.com

  • 摘要: 根据建筑物室内墙壁的空间姿态特征, 建立了一种自动绘制室内平面图的点云向量追踪算法。首先截取一定厚度的室内墙体点云并将其投影至水平面, 在投影平面上建立正方形格网, 对投影后的平面点进行分割管理, 并通过八邻域算法对网格内点云进行聚类。然后采用网格重心法对聚类后的点云数据进行抽稀, 并根据邻域内点间距与连线向量夹角大小对抽稀后的点进行追踪排序, 建立追踪方向的向量序列, 结合向量序列中相邻值夹角的突变情况完成不同墙面的点云分割。最后采用最小二乘算法对分割得到的相同墙面点云进行直线拟合, 求取相邻墙体直线的交点, 建立房屋墙体平面投影线段及其空间连接关系, 依序输出墙体投影线段, 完成房屋平面图的绘制并导出DXF格式数据交换文件。通过对某小区建筑物室内扫描数据的分析, 对所提算法进行验证, 结果表明所提算法可准确快速地对室内三维激光扫描点云进行分析处理, 并完成室内建筑物平面图的绘制。
    Abstract:
      Objectives  According to the spatial attitude characteristics of the indoor walls, a point cloud vector tracking algorithm for automatically drawing the indoor plan was proposed.
      Methods  Firstly, the indoor wall point cloud with a certain thickness was intercepted and projected onto the horizontal plane, and square grids were established on the projection plane. The plane points after projection were segmented and managed, and the point cloud in every grid was clustered by the eight neighborhood algorithm. Then the grid center of gravity method was used to generate sparse points, which were tracked and sorted according to the distance between points in the neighborhood and the angle of connection vector. The vector sequence of tracking direction was established, and the point cloud segmentation of different walls was completed combined with mutation of angle between adjacent values in the vector sequence. Finally, the least squares algorithm was used to fit the line of the same wall point cloud, and the intersection point of the adjacent wall line was obtained.
      Results  The plane projection line segment of the house wall and its spatial connection relationship were established. The wall projection line segment was output in sequence. The drawing of the house plan was completed and the DXF format data exchange file was exported.
      Conclusions  The proposed algorithm was verified by the analysis of indoor scanning data of a residential building. The results show that the proposed algorithm can accurately and quickly analyzed and processed the indoor laser scanning point cloud, and completed the drawing of interior plan.
  • 土壤水分是影响生态、水文和气候等环境过程的重要参数,土壤水分的准确获取和长期监测对于作物生长、旱情监测与预报等具有积极作用。重量法和时域反射法等传统的土壤水分测量方法需要定时定点测量,无法进行大面积长时效的实时观测。遥感技术和理论的不断发展为大范围土壤水分的获取提供了有效途径[1]

    微波遥感具有全天时、全天候工作的特点,且具有一定的穿透能力,已成为土壤水分反演的重要手段。被动微波遥感起步较早、算法成熟,但空间分辨率较低,在较小尺度研究中具有一定的局限性。主动微波遥感中的极化合成孔径雷达(polarimetric synthetic aperture radar,PolSAR)空间分辨率较高,同时可提供对土壤水分反演有利的极化特征参数,国内外学者已对此进行了大量研究。文献[2]提出了一种双分量模型(表面散射和体散射)的极化分解方法,忽略了二面角散射分量,在农田土壤水分反演中取得了较好的效果。文献[3]基于Radarsat-2数据,利用多种分解方法提取多种极化特征参数,对美国黑云杉林的土壤水分进行反演,效果良好。上述研究表明,PolSAR极化特征参数在土壤水分反演中潜力巨大。然而,土壤水分和极化特征参数之间并非简单的线性关系,通常为了提高反演精度,往往考虑较多参数。简单的线性回归分析方法无法准确、全面地衡量大量特征参数在土壤水分反演中的作用,而机器学习方法对于参数的数量和类型没有限制,能够较好地弥补简单回归分析方法在土壤水分反演中的不足。因此,当考虑多极化特征参数时,常常与机器学习方法相结合。

    随机森林(random forest,RF)[4]和支持向量回归机(support vector regression,SVR)[5]是目前较为常用的机器学习方法。文献[6]将后向散射系数作为输入变量、土壤水分作为输出变量,对RF和SVR等12种机器学习方法的性能进行了综合评估,认为RF性能最佳,适用于土壤水分反演。文献[7]基于Radarsat-2数据提取多种极化特征参数并进行特征筛选,评价了RF在土壤水分反演中的适用性。结果表明,RF能较好地模拟冬小麦各生长阶段的土壤水分含量,其性能优于人工神经网络(artifical neural network,ANN)和SVR。综上,机器学习与PolSAR数据相结合能够较好地实现土壤水分模拟。然而,在植被覆盖地区,土壤水分的反演精度常常受到植被覆盖和地表粗糙度的影响,仅考虑极化特征具有一定局限性。

    因此,不少学者协同利用PolSAR数据和光学数据,以提升土壤水分反演精度。文献[8]基于Sentinel-1/2数据,结合水云模型对冬小麦区域的土壤水分进行反演,结果表明,融合植被指数(fusion vegetation index,FVI)的反演效果最优,且能在一定程度上去除冬小麦对后向散射系数的影响。文献[9]基于Sentinel-1/2数据提取后向散射系数和植被指数,利用SVR和广义回归神经网络(generalized regression neural network,GRNN)反演了生长初期冬小麦区的土壤水分,认为加入植被指数的SVR模型效果最优。文献[10]基于Sentinel-1和Landsat 8影像,建立SVR模型,将后向散射系数、归一化植被指数(normalized difference vegetation index,NDVI)、地形因子等作为特征变量,在洪泛平原的水分反演中取得了较好效果,并认为NDVI能够在一定程度上减少植被覆盖的影响。上述研究表明,协同利用PolSAR和光学数据能够在一定程度上减少反演的不确定性,使反演结果更接近真实情况。然而,已有研究多基于C波段SAR数据,对于L波段在土壤水分反演中的效果有待深入[11]

    本文以内蒙古额济纳绿洲为研究区,利用L波段ALOS-2 PALSAR-2影像和Landsat 8影像提取雷达、光学特征参数,通过平均精度下降(mean decrease accuracy,MDA)评估各参数重要性,建立并对比雷达、光学、雷达-光学协同的多参数组合方案的随机森林土壤水分反演模型,对额济纳绿洲区进行土壤水分反演。

    额济纳绿洲位于内蒙古自治区西部,气候干旱,降水稀少,年均气温8.3 ℃,年均降水量约34 mm,潜在蒸发量约3 700~4 000 mm,属温带大陆性气候。2019⁃09⁃24-2019⁃09⁃29,在研究区共布设土壤采样点111个,采集表层0~10 cm的土壤样品,在实验室通过烘干称重法[12]获取各采样点的土壤水分含量。土壤采样点覆盖林地、草地、耕地和裸地等地物类型,实测土壤水分含量为0.02%~15.34%。研究区位置和采样点分布如图1所示。图1采用2019⁃09⁃24的Landsat 8影像,波段组合为7、5、4。

    图  1  研究区位置及采样点分布图
    Figure  1.  Location of the Study Area and Distribution of the Sampling Points

    2014年日本宇航局发射ALOS-2卫星,搭载L波段全极化合成孔径雷达传感器PALSAR-2,卫星轨道高度为628 km,重访周期为14 d。本文采用2019-09-29的ALOS-2 PALSAR-2升轨影像,包含水平极化(horizontal polarization,HH)、垂直极化(vertical polarization,VV)、水平-垂直极化(horizontal-vertical polarization,HV)、垂直-水平极化(vertical-horizontal polarization,VH)这4种极化方式,视向为右视,产品等级为Level 1.1级,中心入射角为33.87°,空间分辨率为6 m。利用SARScape软件首先对ALOS-2影像进行多视、滤波和地理编码,提取后向散射系数σVVσVHσHVσHH;然后利用PolSARpro软件对ALOS-2影像进行滤波处理,通过极化分解提取极化特征参数;最后通过距离多普勒地形校正和重采样,得到与光学影像相匹配的后向散射系数和极化特征参数图像。

    Landsat 8卫星于2013年由美国成功发射,卫星上携带运行陆地成像仪(operational land imager,OLI)和热红外传感器(thermal infrared sensor,TIRS)2种载荷。其中,OLI包括9个波段,除全色波段的空间分辨率为15 m外,其余波段的空间分辨率均为30 m;TIRS包括2个波段,空间分辨率为100 m,卫星重访周期为16 d。本文采用2019-09-24的Landsat 8影像,云量为0.02%。为了消除光照和大气等因素的影响,首先对Landsat 8影像进行辐射定标和大气校正,然后提取地表温度(land surface temperature,LST)和绿叶指数(green leaf index,GLI)等光学特征参数。

    极化分解是PolSAR影像目标特征描述与参数提取的主要方法[13],其目的在于将复杂的散射机制表示为几种具有对应散射矩阵的单一散射分量。根据分解机制和针对目标的不同,已有学者提出了Freeman-Durden分解[14](FD)、Cloude-Pottier分解[15](CP)、Yamaguchi分解[16](YM)、An &Yang分解[17](AY)、H/A/α分解[18]、van Zyl分解[19](VZ)和Bhattacharya & Frery分解[20](BF)等多种极化分解方法。基于ALOS-2 PALSAR影像和上述极化分解方法,提取36个雷达特征参数,具体见表1,其中σVVσHH为同极化后向散射系数,σHVσVH为交叉极化后向散射系数,下标Vol、Odd、Surf、Dbl和Hlx分别代表体散射、奇次、表面、二面角和螺旋体散射分量。

    表  1  雷达特征参数
    Table  1.  Radar Characteristic Parameters
    参数类型处理方法参数个数雷达特征参数
    后向散射系数标准强度和相位处理4σVVσVHσHVσHH
    极化特征参数Freeman-Durden分解[14]3FDVol、FDOdd、FDDbl
    Cloude-Pottier分解[15]3CPDbl、CPSurf、CPVol
    Yamaguchi分解[16]3YM3Vol、YM3Odd、YM3Dbl
    4YM4Vol、YM4Odd、YM4Dbl、YM4Hlx
    An &Yang分解[17]3AY3Vol、AY3Odd、AY3Dbl
    4AY4Vol、AY4Odd、AY4Dbl、AY4Hlx
    H/A/α分解[18]5Alpha、Anisotropy、Entropy、Lambda、RVI
    van Zyl分解[19]3VZVol、VZOdd、VZDbl
    Bhattacharya & Frery分解[20]4BFVol、BFOdd、BFDbl、BFHlx
    下载: 导出CSV 
    | 显示表格

    基于Landsat 8 OLI/TIRS数据,提取11个光学特征参数,其中植被指数通过波段反射率计算得到,LST通过ENVI扩展工具采用大气校正法提取。各参数计算公式参见文献[21-31],包含NDVI、综合光谱响应指数(combined spectral response index,COSRI)、增强型植被指数(enhanced vegetation index,EVI)、GLI、全球植被水分指数(global vegetation moisture index,GVMI)、红外百分比植被指数(infrared percentage vegetation index,IPVI)、植被供水指数(vegetation supply water index,VSWI)、土壤改良植被指数

    (modified soil adjusted vegetation index,MSAVI)、温度植被干旱指数(temperature vegetation

    drought index,TVDI)、比值植被指数(ratio vegetation index,RVI)和LST。

    随机森林是一种基于决策树和Bagging的集成学习方法。决策树是依据一定的划分规则对变量空间进行拆分,并根据相似实例对数据集分组的树形结构。Bagging是一种从原样本集中随机抽取训练样本的技术,这意味着某些样本可能在训练中多次使用,而一些样本可能永远都不会被使用[32],未被使用的样本称为袋外(out of bag,OOB)数据,可用来评估RF的预测效果。MDA是RF中评估特征变量重要性的参数,通过随机打乱OOB数据中某一特征变量值,计算OOB误差,根据其下降程度评估特征变量重要性。通常,OOB数据误差越大,参数越重要。

    本文依据重要性评分结果选取最优特征参数方案建立RF模型,包括2个关键参数,即决策树数量N和节点分裂时输入的特征变量数M。采用N={400,800,1 200,1 600,2 000,2 400},M={1,2,3,4,5,6}的网格搜索,获得最优参数为N=1 600,M=3。模型精度评定采用K折交叉验证,K=10,精度评定指标采用决定系数R2和均方根误差(root mean square error,RMSE)。

    由于特征参数较多,因此本文采用MDA分别对雷达特征参数和光学特征参数的重要性进行评分,结果如图2所示。

    图  2  重要性评分
    Figure  2.  Importance Scores

    图2(a)可以看出,从散射方式看,奇次散射、表面散射和体散射在雷达特征参数中的重要性较高,而二面角散射和螺旋体散射的重要性相对较低。原因可能是由于研究区内部虽有胡杨、柽柳等植被覆盖,但大部分地物还是以低矮草地、戈壁、沙漠为主,因此主要散射机制仍为奇次散射、表面散射和体散射。(1)在奇次散射和表面散射分量中,FD、YM4、AY4、AY3和YM3分解的重要性程度相对较高,在重要性排名中占据前5;而BF、VZ、CP分解相对其他分解方式的重要性程度有所下降,分别占据7~9名。(2)在体散射分量中,YM3、AY3和FD分解的重要性程度相对较高,分别排名12~14位,CP分解的重要性偏低,对土壤水分含量变化不太敏感。(3)二面角散射分量的重要性程度大部分低于体散射分量,但也存在例外,如VZ、AY3分解的二面角散射分量分别优于YM4、AY4的体散射分量,重要性排名分别为15、17位。(4)螺旋体散射分量的重要性偏低,原因是由于该分量的引入是为了解决城区的人造地物不满足反射对称问题[33],研究区内人造地物相对较少,因而其重要性偏低。(5)在H/A/α分解中,Lambda对土壤水分反演的贡献率较高,重要性排名第6,仅次于表面散射分量;RVI相对Lambda重要性程度有所下降,主要原因是RVI能够量化区域植被影响差异[7],但研究区植被较为稀疏,因此重要性偏低;其他参数的重要性排序依次为Entropy>Anisotropy>Alpha。(6)在4个后向散射系数中,同极化后向散射系数的重要性高于交叉极化的后向散射系数,说明同极化后向散射系数相较交叉极化后向散射系数更适合于稀疏植被,与文献[7]在冬小麦返青期的研究较为一致。

    图2(b)可见,光学特征参数的重要性排序依次为VSWI>COSRI>IPVI>LST>RVI>TVDI>MSAVI>NDVI>GVMI>GLI>EVI。VSWI重要性最高,EVI最低,这说明VSWI对于0~10 cm的土壤水分敏感度较高,能够较好地反映土壤水分的变化情况;EVI是为改善高植被区大气校正和土壤背景的影响而提出的[34],研究区多为低矮植被和裸地,因此EVI重要性偏低。MSAVI、NDVI、GVMI和GLI相比EVI重要性程度有所提升,但总体依然较低,可能是GVMI不太适合小尺度区域的干旱监测,NDVI和GLI适用于中等植被覆盖区,在稀疏植被和较高植被覆盖区,其土壤水分敏感性有所降低[35]。相比NDVI,MSAVI能在一定程度上改善土壤背景影响,因此重要性略高于NDVI。其余光学特征参数均对土壤水分反演具有较高的重要性。

    在雷达特征参数中,剔除重要性排名靠后的6个参数,即对土壤水分含量变化敏感度较低的参数,将剩余的30个参数按重要性排序,从VIMR4(VIM表示重要性,下标R表示雷达特征参数,4表示前4个参数)开始,依次放入2个特征变量,分别计算R2和RMSE,选出最优特征组合方案,建立RF模型,各模型精度如图3所示。

    图  3  雷达特征参数组合模型精度对比
    Figure  3.  Accuracy Comparison of Radar Characteristic Parameter Combination Models

    图3可见,基于雷达特征参数组合的14种模型训练集R2普遍高于0.8,RMSE普遍低于1.5%,验证集R2普遍高于0.5,RMSE普遍低于2.5%。在所有模型中,VIMR30的R2和RMSE并非最优,说明模型精度与输入特征数量并非简单的正相关关系,特征选择是提高模型精度的重要环节。从VIMR4至VIMR12,随着输入特征增加,精度显著提升,R2从0.54提升至0.64,RMSE从2.48%降低至2.26%;从VIMR12至VIMR18,虽然精度有所提升,但提升速度相比之前有所减缓,R2由0.64提升至0.67,RMSE由2.26%降低至2.16%;从VIMR18之后,R2和RMSE均呈现不同程度的下降,推测是前后因子间存在强相关性造成了数据冗余,使模型在训练时可能存在轻微过拟合现象,尽管如此,VIMR18的R2和RMSE分别达到0.67和2.16%。C波段极化雷达数据在植被覆盖地表土壤水分研究中应用较多,笔者团队前期对研究区Radarsat-2 C波段数据的研究表明,总体样本中各参数组合方案的RF模型验证集R2均低于0.6,RMSE在3.5%~5%之间[36]。本文基于L波段的参数方案相比C波段在精度上均有提升,尤其是VIMR18效果更佳,说明相比C波段而言,L波段全极化SAR对土壤水分敏感性更高。文献[37-38]对比了不同NDVI状态下ALOS-2 PALSAR-2 L波段雷达数据和Sentinel-1 C波段雷达数据在法国南部、突尼斯中部地区土壤水分反演中的敏感性,认为不同NDVI状态下L波段对土壤水分的敏感程度优于C波段,与本文的研究结果相吻合。因此,选取VIMR18作为最优雷达特征参数组合模型反演土壤水分。

    在光学特征参数中,按重要性排名从VIMO4(下标O代表光学特征参数)开始,每次放入1个特征参数,将其作为RF模型的输入变量,以土壤水分含量作为输出变量,建立RF模型,以R2和RMSE评估模型精度,结果如图4所示。

    图  4  光学特征参数组合模型精度对比
    Figure  4.  Accuracy Comparison of Optical Characteristic Parameter Combination Models

    图4可以看出,基于光学特征参数组合的8种模型训练集R2均在0.89左右,RMSE在1.25%左右,验证集R2在0.5左右,RMSE在2.47%左右。从VIMO4至VIMO11,验证集R2呈现先缓慢增长、后逐渐降低、再增长、又降低的趋势;RMSE则呈现先降低、后增长、再降低、又增长的趋势;从VIMO6至VIMO9,模型精度的第一次降低可能是由于因子间的强相关性导致模型训练集的过度拟合所致;从VIMO9至VIMO11,验证集R2呈现先增长、后下降的趋势,但对比VIMO9和VIMO11两组方案可知,VIMO11虽然R2有所提升,但RMSE也呈现上升趋势,可能是由于采样点土壤水分低值数量较多、高值数量过少,导致出现低值点高估、高值点低估现象,最终导致验证集RMSE有所上升。由于在光学特征参数组合方案中综合使用R2和RMSE无法确定最优组合方案,因此,选取验证集R2位居前3位的VIMO6、VIMO7、VIMO10光学特征参数组合方案进行土壤水分反演。

    将上述优选的雷达和光学特征参数组合方案进一步组合,得到3种协同反演方案,将3种方案的反演值分别与验证点的土壤水分含量实测值进行对比分析,计算R2和RMSE,结果见表2

    表  2  雷达-光学参数协同反演模型精度对比
    Table  2.  Accuracy Comparison of Radar-Optical Parameter Integrated Inversion Model
    方案组合类型R2RMSE/%
    方案1VIMR18+ VIMO60.721.99
    方案2VIMR18+VIMO70.712.00
    方案3VIMR18+VIMO100.702.02
    下载: 导出CSV 
    | 显示表格

    表2可知,3种协同反演模型R2均不低于0.7,RMSE均低于2.05%。其中方案1的R2为0.72,RMSE为1.99%,在3种协同模型中精度最高。与仅使用雷达特征参数的最佳模型VIMR18相比,R2提高7.46%,RMSE降低8.54%,表明光学特征参数能够在一定程度上消除植被覆盖的影响,弥补雷达数据在土壤水分反演中的不足。与仅使用光学特征参数的模型VIMO6相比,R2提升38.4%,RMSE降低22.6%,说明相较于单一数据源而言,多源数据融合能够在一定程度上提升土壤水分反演精度,使预测结果更接近真实情况。相比方案1,方案2和方案3总体精度虽然有所下降,但仍然较为理想,进一步说明融合多源数据的机器学习方法在土壤水分反演中的应用潜力,其适用范围相较单一数据源更为广泛。

    图5为基于雷达、光学以及雷达-光学特征参数协同反演模型的土壤水分反演结果。由图5可以看出,仅使用光学特征参数的模型(图5(b))存在水分整体高估的现象,特别是在裸地和部分低植被覆盖区,这些区域易受土壤背景和天气的影响[35],从而导致预测精度较低。相比图5(b),基于雷达特征参数(图5(a))和雷达-光学特征参数协同(图5(c))的模型对于裸地和低植被覆盖区的估算效果较好,与实测水分含量基本一致,表明雷达特征参数在裸地和低植被区土壤水分反演中的应用潜力。已有研究表明,对于植被覆盖较高地区而言,仅使用雷达特征参数难以完全消除植被覆盖和地表粗糙度的影响,由图5(a)可以看出,东部植被覆盖区土壤水分含量相比实测值有所低估,加入光学特征参数后(图5(c)),低估现象有所改善,表明多源数据融合能在一定程度上克服因影响因子考虑不足导致的精度受限问题,这与文献[39]的研究结果较为吻合。但由于研究区土壤水分含量整体较低,导致土壤水分高值采样点较少,高值低估现象仍然存在。总体来说,图5能反映研究区土壤水分含量的整体趋势,其中雷达-光学协同模型反演效果最佳,能够更好地反映研究区土壤水分的空间分布格局。

    图  5  多种模型土壤水分反演结果
    Figure  5.  Soil Moisture Inversion Results of Various Models

    本文基于L波段ALOS-2 PALSAR-2和Landsat 8影像提取雷达特征参数和光学特征参数,通过MDA优选特征参数,建立基于雷达、光学以及雷达-光学特征参数协同的RF土壤水分反演模型,定量评估模型精度,反演绿洲土壤水分。结论如下:(1)相比C波段而言,L波段对干旱荒漠绿洲地区的土壤水分含量敏感性更高;(2)雷达特征参数中,表面散射分量重要性最高,体散射分量次之,二面角和螺旋体散射分量重要性较低,4个后向散射系数中,同极化后向散射系数的重要性优于交叉极化后向散射系数,光学特征参数中,VSWI重要性最高,EVI最低;(3)雷达特征参数方案最优模型R2、RMSE分别为0.67、2.16%。光学特征参数方案模型精度普遍较低且精度相当,R2、RMSE分别为0.5、2.47%;(4)相比单一数据源而言,协同利用雷达和光学特征参数,模型适用范围更广,预测精度更高。最优协同反演模型的R2、RMSE分别为0.72、1.99%。相比仅使用雷达特征参数模型VIMR18,R2提升了7.46%,RMSE降低了8.54%,相比仅使用光学特征参数模型VIMO6,R2提升了38.4%,RMSE降低了22.6%。

    本文研究表明,雷达-光学协同的随机森林模型在干旱荒漠绿洲区土壤水分反演中具有良好的适用性,能够较好地反映区域土壤水分空间分布格局。后期将采用更多的机器学习方法,探索协同利用L波段全极化雷达数据和光学数据在较深层次土壤水分含量反演中的效果和适用性。

  • 图  1   建筑物室内平面图绘制流程图

    Figure  1.   Flowchart of Building Interior Plan Drawing

    图  2   点云投影示意图

    Figure  2.   Diagrams of Point Cloud Projection

    图  3   基于八邻域的点云聚类

    Figure  3.   Point Cloud Clustering Based on Eight Neighborhood

    图  4   点云抽稀示意图

    Figure  4.   Diagram of Point Cloud Extraction

    图  5   点云排序流程图

    Figure  5.   Flowchart of Point Cloud Sorting

    图  6   点云排序示意图

    Figure  6.   Diagram of Point Cloud Sorting

    图  7   点云分割流程图

    Figure  7.   Flowchart of Point Cloud Segmentation

    图  8   点云分割并求交点

    Figure  8.   Segmentation and Intersection of Point Clouds

    图  9   窗户拐角示意图

    Figure  9.   Diagram of Tiny Corner

    图  10   室内平面图绘制结果

    Figure  10.   Results of Interior Plan Drawing

    图  11   门与窗户识别流程图

    Figure  11.   Flowchart of Door and Window Identification

    图  12   室内点云的截面与投影

    Figure  12.   Cross-Section and Projection of Indoor Point Cloud

    图  13   点云聚类单元

    Figure  13.   Point Cloud Clustering Units

    图  14   不同算法对墙体点云的处理结果对比

    Figure  14.   Processing Results of Wall Point Clouds by Different Algorithms

    图  15   向量追踪算法室内平面图绘制结果

    Figure  15.   Result of Point Cloud Vector Tracing Algorithm for Drawing Interior Plan

    表  1   不同算法绘图效果对比

    Table  1   Results Comparison of Different Algorithms

    算法 精度/% 召回率/% F1综合评价系数
    IEPF算法 95 83 0.88
    HT算法 71 100 0.83
    向量追踪算法 100 100 1
    下载: 导出CSV

    表  2   不同长度的直线识别准确率

    Table  2   Line Recognition Accuracy of Different Lengths

    墙壁长度/cm 墙壁数量 正确识别数量 正确率/%
    (5, 15] 18 16 89
    (15, 25] 19 19 100
    (25, 35] 25 25 100
    > 35 28 28 100
    下载: 导出CSV

    表  3   拐角点位误差统计

    Table  3   Statistics of Corner Point Position Errors

    点位误差/mm 拐角数量
    (0, 1.5] 13
    (1.5, 3.0] 26
    (3.0, 4.5] 29
    (4.5, 10] 10
    下载: 导出CSV
  • [1]

    Weerakoon T, Ishii K, Nassiraei A A F. Geometric Feature Extraction from 2D Laser Range Data for Mobile Robot Navigation[C]//2015 IEEE 10th International Conference on Industrial and Information Systems(ICⅡS), Peradeniya, Sri Lanka, 2015

    [2] 危双丰, 刘明蕾, 赵江洪, 等. 利用点云检测室内导航元素的方法综述[J]. 武汉大学学报·信息科学版, 2018, 43(12): 2 003- 2 011 doi: 10.13203/j.whugis20180144

    Wei Shuangfeng, Liu Minglei, Zhao Jianghong, et al. A Survey of Methods for Detecting Indoor Navigation Elements from Point Clouds[J]. Geomatics and Information Science of Wuhan University, 2018, 43(12): 2 003-2 011 doi: 10.13203/j.whugis20180144

    [3] Nguyen Tien Thanh, 刘修国, 王红平, 等. 基于激光扫描技术的三维模型重建[J]. 激光与光电子学进展, 2011, 48 (8): 112-117 https://www.cnki.com.cn/Article/CJFDTOTAL-JGDJ201108016.htm

    Nguyen Tien Thanh, Liu Xiuguo, Wang Hongping, et al. 3D Model Reconstruction Based on Laser Scanning Technique[J]. Laser & Optoelectronics Progress, 2011, 48(8): 112-117 https://www.cnki.com.cn/Article/CJFDTOTAL-JGDJ201108016.htm

    [4] 王鹏, 刘如飞, 马新江, 等. 一种车载激光点云中杆目标自动提取方法[J]. 武汉大学学报·信息科学版, 2020, 45 (7): 1 035-1 042 doi: 10.13203/j.whugis20170421

    Wang Peng, Liu Rufei, Ma Xinjiang, et al. An Automatic Extraction Method for Pole-Like Objects from Vehicle-Borne Laser Point Cloud[J]. Geomatics and Information Science of Wuhan University, 2020, 45(7): 1 035-1 042 doi: 10.13203/j.whugis20170421

    [5] 朱庆, 李世明, 胡翰, 等. 面向三维城市建模的多点云数据融合方法综述[J]. 武汉大学学报·信息科学版, 2018, 43 (12): 1 962-1 971 doi: 10.13203/j.whugis20180109

    Zhu Qing, Li Shiming, Hu Han, et al. Multiple Point Clouds Data Fusion Method for 3D City Modeling[J]. Geomatics and Information Science of Wuhan University, 2018, 43(12): 1 962-1 971 doi: 10.13203/j.whugis20180109

    [6]

    Garulli A, Giannitrapani A, Rossi A, et al. Mobile Robot SLAM for Line-Based Environment Representation[C]//44th IEEE Conference on Decision and Control, Seville, Spain, 2005

    [7] 满增光, 叶文华, 楼佩煌, 等. 基于自适应阈值的距离图像线段特征提取[J]. 深圳大学学报(理工版), 2011, 28 (6): 483-488 doi: 10.3969/j.issn.1000-2618.2011.06.003

    Man Zengguang, Ye Wenhua, Lou Peihuang, et al. Line Feature Extraction Based on Adaptive Threshold in Range Images[J]. Journal of Shenzhen University Science and Engineering, 2011, 28(6): 483-488 doi: 10.3969/j.issn.1000-2618.2011.06.003

    [8]

    Borges G A, Aldon M J. Line Extraction in 2D Range Images for Mobile Robotics[J]. Journal of Intelligent and Robotic Systems, 2004, 40 (3): 267-297 doi: 10.1023/B:JINT.0000038945.55712.65

    [9]

    Vandorpe J, Van Brussel H, Xu H. Exact Dynamic Map Building for a Mobile Robot Using Geometrical Primitives Produced by a 2D Range Finder[C]//IEEE International Conference on Robotics and Automation, Minneapolis, MN, USA, 1996

    [10]

    Borges G A, Aldon M J. A Split-and-Merge Segmentation Algorithm for Line Extraction in 2D Range Images[C]//15th International Conference on Pattern Recognition (ICPR-2000), Barcelona, Spain, 2000

    [11]

    Nguyen V, Gächter S, Martinelli A, et al. A Comparison of Line Extraction Algorithms Using 2D Range Data for Indoor Mobile Robotics[J]. Autonomous Robots, 2007, 23 (2): 97-111 doi: 10.1007/s10514-007-9034-y

    [12]

    Choi Y H, Lee T K, Oh S Y. A Line Feature Based SLAM with Low Grade Range Sensors Using Geometric Constraints and Active Exploration for Mobile Robot[J]. Autonomous Robots, 2008, 24 (1): 13-27 doi: 10.1007/s10514-007-9050-y

    [13] 田朋举, 梅文胜, 花向红, 等. 一种基于点云数据快速生成建筑物断面图的方法[J]. 测绘地理信息, 2020, 45 (4): 88-92 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXG202004019.htm

    Tian Pengju, Mei Wensheng, Hua Xianghong, et al. A Method for Generating Building Sectional Drawing Using 3D Laser Point Cloud Data[J]. Journal of Geomatics, 2020, 45(4): 88-92 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXG202004019.htm

    [14]

    He X, Cai Z H. Feature Extraction from 2D Laser Range Data for Indoor Navigation of Aerial Robot [C]//2013 Chinese Automation Congress, Changsha, China, 2013

    [15] 湛逸飞, 杨必胜. 基于LiDAR点云法向量稳健估计的室内平面图自动生成方法[J]. 地理信息世界, 2017, 24 (1): 92-96 doi: 10.3969/j.issn.1672-1586.2017.01.018

    Zhan Yifei, Yang Bisheng. An Automatic Indoor Plan Generation Method Based on Robust LiDAR Point Cloud Normal Estimation[J]. Geomatics World, 2017, 24 (1): 92-96 doi: 10.3969/j.issn.1672-1586.2017.01.018

    [16] 郑德华. ICP算法及其在建筑物扫描点云数据配准中的应用[J]. 测绘科学, 2007, 32 (2): 31-32 https://www.cnki.com.cn/Article/CJFDTOTAL-CHKD200702008.htm

    Zheng Dehua. ICP Algorithm and Application in the Data Registration of Building Scanning Point Cloud [J]. Science of Surveying and Mapping, 2007, 32 (2): 31-32 https://www.cnki.com.cn/Article/CJFDTOTAL-CHKD200702008.htm

    [17] 何君. 基于局部拟合平面投影搜索最近点的ICP配准[J]. 测绘地理信息, 2019, 44 (4): 86-89 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXG201904021.htm

    He Jun. ICP Registration Algorithm Based on Closet Point Searching in Local Plane Projection[J]. Journal of Geomatics, 2019, 44(4): 86-89 https://www.cnki.com.cn/Article/CJFDTOTAL-CHXG201904021.htm

    [18] 屈亚运, 程英蕾, 邱浪波. 改进的八邻域搜索提取建筑物立体特征方法[J]. 空军工程大学学报(自然科学版), 2015, 16 (4): 66-69 https://www.cnki.com.cn/Article/CJFDTOTAL-KJGC201504016.htm

    Qu Yayun, Cheng Yinglei, Qiu Langbo. An Improved Method of Extracting 3D Features of Buildings Based on Eight Neighborhood Search Method[J]. Journal of Air Force Engineering University(Natural Science Edition), 2015, 16(4): 66-69 https://www.cnki.com.cn/Article/CJFDTOTAL-KJGC201504016.htm

    [19] 衡星辰, 罗俊颉, 郭俊文, 等. 八邻域网格聚类的多样性XML文档近似查询算法[J]. 西安交通大学学报, 2007, 41(8): 907-911 https://www.cnki.com.cn/Article/CJFDTOTAL-XAJT200708009.htm

    Heng Xingchen, Luo Junjie, Guo Junwen, et al. Approximate Query Algorithm Based on EightNeighbor Grid Clustering for Heterogeneous XML Documents[J]. Journal of Xi'an Jiaotong University, 2007, 41(8): 907-911 https://www.cnki.com.cn/Article/CJFDTOTAL-XAJT200708009.htm

    [20] 彭晨, 余柏蒗, 吴宾, 等. 基于移动激光扫描点云特征图像和SVM的建筑物立面半自动提取方法[J]. 地球信息科学学报, 2016, 18 (7): 878-885 https://www.cnki.com.cn/Article/CJFDTOTAL-DQXX201607004.htm

    Peng Chen, Yu Bailang, Wu Bin, et al. A Method for Semiautomated Segmentation of Building Facade from Mobile Laser Scanning Point Cloud Based on Feature Images and SVM[J]. Journal of Geo-Information Science, 2016, 18(7): 878-885 https://www.cnki.com.cn/Article/CJFDTOTAL-DQXX201607004.htm

图(15)  /  表(3)
计量
  • 文章访问数:  1427
  • HTML全文浏览量:  415
  • PDF下载量:  97
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-10-25
  • 发布日期:  2021-09-17

目录

/

返回文章
返回