具有孔洞的地下电缆工井模型拓扑重构研究

赵江洪, 董岩, 黄明, 张晓光, 马思宇, 孙铭悦

赵江洪, 董岩, 黄明, 张晓光, 马思宇, 孙铭悦. 具有孔洞的地下电缆工井模型拓扑重构研究[J]. 武汉大学学报 ( 信息科学版), 2019, 44(12): 1849-1858. DOI: 10.13203/j.whugis20180126
引用本文: 赵江洪, 董岩, 黄明, 张晓光, 马思宇, 孙铭悦. 具有孔洞的地下电缆工井模型拓扑重构研究[J]. 武汉大学学报 ( 信息科学版), 2019, 44(12): 1849-1858. DOI: 10.13203/j.whugis20180126
ZHAO Jianghong, DONG Yan, HUANG Ming, ZHANG Xiaoguang, MA Siyu, SUN Mingyue. Topological Reconstruction of Underground Cable Well Model with Holes[J]. Geomatics and Information Science of Wuhan University, 2019, 44(12): 1849-1858. DOI: 10.13203/j.whugis20180126
Citation: ZHAO Jianghong, DONG Yan, HUANG Ming, ZHANG Xiaoguang, MA Siyu, SUN Mingyue. Topological Reconstruction of Underground Cable Well Model with Holes[J]. Geomatics and Information Science of Wuhan University, 2019, 44(12): 1849-1858. DOI: 10.13203/j.whugis20180126

具有孔洞的地下电缆工井模型拓扑重构研究

基金项目: 

国家重点研发计划 2016YFC0802107

国家自然科学基金 41601409

国家自然科学基金 41501495

北京市自然科学基金 8172016

城市空间信息工程北京市重点实验室开放研究基金 2018210

北京建筑大学科学研究基金 00331616056

无人机倾斜摄像及在教学中的应用研究 ZF16095

详细信息
    作者简介:

    赵江洪, 博士, 副教授, 主要研究方向为三维激光雷达数据处理。zhaojiangh@bucea.edu.cn

    通讯作者:

    黄明, 博士, 教授。huangming@bucea.edu.cn

  • 中图分类号: P208

Topological Reconstruction of Underground Cable Well Model with Holes

Funds: 

The National Key Research and Development Program of China 2016YFC0802107

the National Natural Science Foundation of China 41601409

the National Natural Science Foundation of China 41501495

Beijing Natural Science Foundation 8172016

the Open Research Fund of Beijing Key Laboratory of Urban Spatial Information Engineering 2018210

the BUCEA Scientific Research Foundation 00331616056

Research on Tilt Camera of UAV and Its Application in Teaching ZF16095

More Information
    Author Bio:

    ZHAO Jianghong, PhD, associate professor. His main research direction is three-dimensional LiDAR data processing. E-mail:zhaojiangh@bucea.edu.cn

    Corresponding author:

    HUANG Ming, PhD, professor. E-mail:huangming@bucea.edu.cn

  • 摘要: 城市地下电缆工井是发挥城市功能和确保城市快速协调发展的重要基础设施,但地下管线结构复杂、信息量大,对其进行三维显示、管理与分析具有一定的难度。激光雷达技术出现后,因其快速、不接触、高密度、高精度等特征,被逐步应用于建筑物逆向重建,基于激光雷达点云数据的地下电缆工井三维重建也逐渐成为主流方法。提出了一种针对具有孔洞的地下工井的三维模型拓扑重构算法,以地下工井三维激光点云为数据源,建立具有孔洞的地下工井三维模型,可解决直接采用点云进行构网时模型不完整及速度较慢的问题。三维模型的构建采用与常用的管线化地下管网结构不同的CSG-BREP(constructive solid geometry-boundaryrepresentation)拓扑模型结构,整体上将工井分为墙面、电缆、井筒等,此模型可更加细致地表达地下电缆工井模型的内部拓扑结构。实验证明,该算法可以对多种类型的墙壁进行拓扑重构,在表达模型拓扑关系时突破了传统管网中单一的管孔与管线的点线关系及管线间的线线关系,增加管孔与墙壁的点面关系以及管线与墙壁的线面关系,更加细致准确地对模型进行描述。
    Abstract: Urban underground cable well, as an important infrastructure, contributes to the performance of urban functions and the coordinated urban development. However, the underground pipeline's complex structure and huge amount of data present great difficulty in displaying, management and analysis. With the development of LiDAR technology characterized by high-speed, non-contact, high-density, and high-precision, point cloud data are gradually introduced into the field of building reverse reconstruction. The threedimensional reconstruction of underground cable wells based on point cloud data has gradually gathered momentum. An algorithm for topological reconstruction of underground cable well model with holes is proposed, which can overcome the problem of model incompleteness and modeling inefficiency when point cloud as data source is used to build triangular network directly. A new hybrid model-CSG-BREP (constructive solid geometry-boundary representation) topological model is adopted to construct the 3D model, which different from the commonly used pipe-lined underground pipe network structure. In this new model, the wells are divided into walls, cables, wellbore, et al, this model can more carefully express the internal topological structure of the underground cable well. Experimental results show that the proposed algorithm can reconstruct many types of walls topologically, and the point-surface relationship between pipe holes and walls and the line-surface relationship between pipes and walls are added into the point-line and line-surface relationship of traditional pipe networks, which allows for a more accurate and detailed description.Through the test of the cable management project, the model shows good results.
  • 每个地域往往具有特定的建筑风格特征,精细化探究建筑风格的地理分布有利于历史古建筑的保护、特色旅游资源的开发以及城市建筑区的科学规划。然而城市建筑数量多、区域广,人工采集和调查不满足大范围作业的需求。近年来,以谷歌为代表的互联网公司推出城市街景影像(street view images,SVI)。街景影像分辨率高,含有全方位的城市街道景观以及精确的位置和姿态信息,为探究城市建筑风格的地理分布提供了可能。

    街景建筑区域的风格识别与地图匹配是生成城市建筑风格地图的技术难点。一方面,建筑风格特征复杂多样,图像分类困难。相比传统分类模型[1-4],卷积神经网络[5](convolutional nueral network,CNN)在建筑风格分类上的精度更高。另一方面,建筑图像的跨视角匹配成为建筑地理定位和土地利用分类[6]等领域的重要手段,现有研究主要从单一建筑的空间几何特征考量[7-11],不适合处理场景复杂、建筑众多的全景像片。因此,本文借助深度学习实现街景建筑区域的风格分类,并建立建筑影像与建筑俯视轮廓的映射关系,精细化探究大范围城市建筑风格地图的生成方法。

    更快的基于区域CNN(faster region-based CNN,Faster R-CNN)目标检测模型[12]可以提取各类风格的建筑目标区域。在此基础上,采用双像建筑区域点位映射方法,通过匹配相邻两张街景影像的同名建筑区域来估算每个单体建筑的实地位置,建立与建筑轮廓的点位映射关系。针对无同名建筑区域的建筑影像,本文提出了单像建筑区域方位映射方法,通过街景与地图中建筑物方位角的对应关系实现单体建筑影像与对应建筑轮廓的匹配。对成功匹配的建筑轮廓获取其建筑风格属性,生成城市建筑风格地图。本文方法的流程见图 1图 1中,PASCAL VOC是国际顶级的计算机视觉竞赛提供的标准格式的目标检测数据集。

    图  1  建筑风格地图生成方法流程图
    Figure  1.  Flowchart of Architectural Style Map Production Method

    同名建筑区域是指相邻两张街景中代表同一建筑物的影像区域,是进行双像建筑区域点位映射的必要条件(见图 2)。本文首先通过建筑影像的纹理和颜色特征剔除相似度小的单体建筑区域影像对,然后对少量的候选影像对进行核线约束下的加速稳健特征[13](speeded up robust features,SURF)匹配,得到最佳同名建筑区域及核线距离最近的一个同名点对,用于前方交会定位,匹配流程如图 3所示。

    图  2  同名建筑区域示例
    Figure  2.  An Example of Two Corresponding Architectural Images
    图  3  同名建筑区域匹配流程图
    Figure  3.  Flowchart of Matching Two Corresponding Architectural Images

    街景单体建筑影像的纹理特征采用不同方向上的灰度共生矩阵[14](gary level co-occurrence matrix,GLCM)的熵、二阶矩、对比度和逆差矩表示,本文选取0°、45°、90°和135°这4个方向的灰度共生矩阵,构成16维的纹理特征向量。色调、饱和度、亮度(hue,saturation,value,HSV)颜色空间的低阶矩能很好地反映图像明暗程度、颜色分布范围和分布对称性,共生成9维的颜色特征向量。将两个特征向量融合,得到代表建筑区域的25维全局特征。

    两个单体建筑区域特征向量的差向量可作为判定相似程度的特征向量。首先,将差向量作为后验概率的支持向量机[15](posteriori probability support vector machine,PPSVM)的输入,计算同名建筑区域的概率值,将概率值大于阈值的同类风格建筑区域对作为候选匹配对。然后,在候选匹配对中,以左片的单体建筑区域为基准,分别与右片的每个单体建筑区域进行球面核线约束[16]下的SURF匹配,选取同名点核线距离均值最小的右片建筑区域作为最佳匹配对。最后,在每个最佳匹配对中,取核线距离最小的一对同名点,通过前方交会得到建筑物的坐标信息。

    街景建筑影像形变较大,存在同名点匹配困难、误匹配较多的问题。对此,本文通过检验左、右片单体建筑影像区域中心点的核线距离是否小于阈值来辅助判别是否为同名建筑区域。

    建立街景影像中建筑区域与地图中建筑轮廓的对应关系,是生成建筑风格地图的关键。本文提出了双像建筑区域点位映射和单像建筑区域方位映射方法,流程如图 4所示。

    图  4  街景建筑区域映射方法流程图
    Figure  4.  Flowchart of Building Outline Mapping Methods Based on Architectural Area in Street View

    双像建筑区域点位映射需要满足两个条件:(1)相邻两张街景影像存在同名单体建筑区域;(2)两张影像可以完成相对定向。图 5显示了同名建筑区域点位映射的过程。图 5中,$ {P}_{1} $和$ {P}_{2} $为同名建筑区域的最佳同名点对;$ {S}_{1} $和$ {S}_{2} $是相邻两张街景影像的摄影中心,分别构成球形全景坐标系S1-X1Y1Z1S2-X2Y2Z2;$ P $是$ {P}_{1} $和$ {P}_{2} $对应的前方交会点。

    图  5  双像建筑区域点位映射示意图
    Figure  5.  Location Mapping Method Based on a Panoramic Image Pair

    完成两张街景的相对定向需要将平面街景影像转换成球形全景影像[17-18]。首先,采用随机抽样一致性(random sample consensus,RANSAC)的8点法[16, 19],由式(1)解算相邻两张街景的旋转矩阵和平移向量。然后,将同名点对进行前方交会,利用街景影像的坐标和姿态参数完成绝对定向[17],得到对应单体建筑的实地坐标。

    $$ \left\{\begin{array}{l}{\mathit{P}}_{2}\cdot \underset{3\times 3}{\mathit{E}}\cdot {\mathit{P}}_{1}=0\\ \mathit{E}={\left[\mathit{T}\right]}_{\times }\cdot {\mathit{R}}_{s}\end{array}\right. $$ (1)

    式中,$ \mathit{E} $为本质矩阵;$ {\mathit{P}}_{1} $和$ {\mathit{P}}_{2} $表示至少8对同名点球面坐标的矩阵;$ {\left[\mathit{T}\right]}_{\times } $表示平移向量$ \mathit{T} $的反对称矩阵;$ {\mathit{R}}_{s} $为像对的旋转矩阵。

    由于无法进行严格的绝对定向,估算的单体建筑位置点与建筑轮廓存在偏差。首先,通过限制前方交会距离小于150 m和交会角大于10°提高定位精度[20-21]。然后,采用空间近邻分析,使每个位置点匹配距离最近的单体建筑轮廓,且限制匹配距离小于10 m,提高单体建筑的匹配准确率。

    当球形全景像空间坐标系旋转至与地图坐标系平行时,单体建筑区域的方位范围与街景采集点到对应建筑轮廓的方位范围相同。图 6为建筑区域方位映射示意图,$ {S}_{1} $为街景影像的摄影中心,构成球形全景坐标系S1-X1Y1Z1

    图  6  单像建筑区域方位映射示意图
    Figure  6.  Diagram of Azimuth Mapping Method Based on a Panoramic Image

    在实际场景中,检测出的单体建筑区域范围存在偏差,地图中对应的建筑俯视轮廓也可能存在位置偏移,导致两个方位范围不严格相同。定义两个方位范围的交并比(intersection over union,IoU),用于描述单体建筑区域与建筑俯视轮廓的方位相似度(见图 7)。

    图  7  两个方位范围交并比的定义
    Figure  7.  Definition of Two Azimuth Coverage ?s IoU

    交并比$ \alpha $越接近1,则相似度越大;越接近0,则相似度越小。计算公式为:

    $$ \alpha =\frac{I\left[\right({a}_{i}, {b}_{i}), \mathrm{ }({c}_{j}, {d}_{j}\left)\right]}{U\left[\right({a}_{i}, {b}_{i}), \mathrm{ }({c}_{j}, {d}_{j}\left)\right]}\begin{array}{c}\end{array}(i, j=\mathrm{1, 2}, 3\cdots ) $$ (2)

    式中,$ ({a}_{i}, {b}_{i}) $表示第$ i $个建筑区域的方位角范围;$ ({c}_{j}, {d}_{j}) $表示街景采集点到第$ j $个建筑轮廓的方位角范围;$ I $表示两个方位范围的重合角度;$ U $表示两个方位范围的合并角度。

    所有单体建筑区域的方位范围构成集合$ {A}_{1}=\left\{\right({a}_{i}, {b}_{i}\left)\mathrm{ }\right|i=\mathrm{1, 2}, 3\cdots \} $。采集点到单个建筑轮廓的方位范围可以由采集点到建筑轮廓的两个边界方位角顺时针组成(见图 8),获取缓冲区内所有单体建筑轮廓的方位角范围,构成集合$ {A}_{2}=\left\{\right({c}_{j}, {d}_{j}\left)\mathrm{ }\right|j\in B\} $,$ B $表示缓冲区内的单体建筑轮廓集合。利用式(2)选取IoU最大的匹配对,完成单像建筑区域方位映射。

    图  8  点到建筑俯视轮廓的方位范围示意图
    Figure  8.  Diagram of Azimuth Coverage from One Position to Building Outline

    图 9描述了单个建筑轮廓匹配多个不同建筑区域影像的情形,其中$ {S}_{i}(i=\mathrm{1, 2}\dots 6) $代表各街景采集点位置,$ {P}_{1} $和$ {P}_{2} $代表两个建筑区域的映射点位置。单个建筑轮廓可能映射多类建筑风格(简称多映射问题),故需要确定每个建筑轮廓的唯一风格属性。

    图  9  单体建筑轮廓匹配多个建筑区域影像的示意图
    Figure  9.  Diagram of One Building Outline Matching Mutiple Architectural Images

    对于双像建筑区域点位映射结果,每个建筑轮廓可能匹配到多个建筑区域位置点。针对一个单体建筑轮廓,将匹配到的每类建筑风格的位置点各构成一个集合,考虑每个集合中所有位置点的数量$ {m}_{1} $、Faster R-CNN分类概率均值$ {m}_{2} $、到建筑轮廓的距离均值$ {m}_{3} $以及用于交会的同名点对的核线距离均值$ {m}_{4} $,用这4个参数构成每类建筑风格的评价指标,且权重相同。当一类建筑风格的$ {m}_{1} $和$ {m}_{2} $越大,$ {m}_{3} $和$ {m}_{4} $越小,该类风格的得分越高。将指标归一化、正向化后,构成的标准化矩阵$ \mathit{M} $为:

    $$ \mathit{M}=\left[\begin{array}{cccc}{m}_{\mathrm{1, 1}}& {m}_{\mathrm{1, 2}}& {m}_{\mathrm{1, 3}}& {m}_{\mathrm{1, 4}}\\ {m}_{\mathrm{2, 1}}& {m}_{\mathrm{2, 2}}& {m}_{\mathrm{2, 3}}& {m}_{\mathrm{2, 4}}\\ ⋮& ⋮& ⋮& ⋮\\ {m}_{\mathrm{19, 1}}& {m}_{\mathrm{19, 2}}& {m}_{\mathrm{19, 3}}& {m}_{\mathrm{19, 4}}\end{array}\right] $$ (3)

    式中,$ {m}_{i, j} $表示第$ i $类建筑风格的第$ j $个指标的标准化值。

    结合优劣解距离方法[22](technique for order preference by similarity to an ideal solution,TOPSIS),第$ i $类建筑风格的得分计算式如式(4)所示,取得分最高的建筑风格为该单体建筑轮廓的风格属性$ {s}_{i} $。

    $$ {s}_{i}=\frac{\sqrt{\sum _{j=1}^{4}({m}_{j}^{+}-{m}_{i, j}{)}^{2}}}{\sqrt{\sum _{j=1}^{4}({m}_{j}^{+}-{m}_{i, j}{)}^{2}}+\sqrt{\sum _{j=1}^{4}({m}_{j}^{-}-{m}_{i, j}{)}^{2}}} $$ (4)

    式中,$ {\mathit{m}}^{+} $为4维向量,且$ {\mathit{m}}^{+} $=[max($ {m}_{\mathrm{1, 1}} $,$ {m}_{\mathrm{2, 1}}\cdots {m}_{\mathrm{19, 1}} $),max($ {m}_{\mathrm{1, 2}} $,$ {m}_{\mathrm{2, 2}}\cdots {m}_{\mathrm{19, 2}} $)$ \cdots $];$ {\mathit{m}}^{-} $也为4维向量,且$ {\mathit{m}}^{-} $=[min($ {m}_{\mathrm{1, 1}} $,$ {m}_{\mathrm{2, 1}}\cdots {m}_{\mathrm{19, 1}} $),min($ {m}_{\mathrm{1, 2}} $,$ {m}_{\mathrm{2, 2}}\cdots {m}_{\mathrm{19, 2}} $)$ \cdots $]。

    针对单像建筑区域方位映射结果,考虑每个建筑区域影像与该建筑俯视轮廓的方位交并比和建筑区域影像的风格分类概率,利用式(4)确定建筑轮廓的唯一风格属性。

    本文对北京、上海、武汉、西安4座城市进行建筑风格地图生成实验,实验区域选取各类建筑风格分布密集的中心城区,并获取了2017年间共73 420张百度街景影像,以及共含约26.29万个建筑俯视轮廓的矢量地图数据(.shapfile格式)。每张街景影像分辨率为2 048×1 024像素,相邻两张影像的实地距离为8~20 m,两个街景影像对的实地距离为50~150 m,每张街景影像带有地理信息标签,包括影像拍摄地的经纬度、正北方向在影像中的方位角、影像的拍摄姿态信息和相邻街景的唯一标识符。

    为了训练检测街景建筑风格的Faster R-CNN网络,收集了4座城市中含有19类建筑风格的65个代表性景点或街区,如图 10所示。人工标注生成7 296张影像的PASCAL VOC样本集,训练样本和测试样本的数量比例为4∶1。图 11显示了各类建筑风格的原始标定数量,为了使训练样本数量均衡,将12类数量较少的建筑区域影像水平翻转,使样本数量增加一倍。

    图  10  实验区域及代表性街区或景点的位置分布
    Figure  10.  Experimental Region and Location Distribution of Representative Blocks or Scenic Spots
    图  11  各类建筑风格的原始标定数量
    Figure  11.  Number of Calibration of Different Architectural Styles

    为训练匹配同名建筑区域的支持向量机模型,人工截取200对同名建筑区域影像,每对同名建筑区域影像作为一个正样本,共获得200个正样本数据。同时,任意不同名建筑区域影像两两配对,共组成200个负样本。训练集和测试集的数量比例为3∶1。实验所采用的计算硬件配置为Intel 8 GB CPU,NVIDIA 1660Ti 6 GB GPU,算法由Python 3.6实现,深度学习框架采用TensorFlow-GPU 1.13,机器学习框架利用scikit-learn库,图像处理通过OpenCV库完成。

    Faster R-CNN网络的参数设置为:权重衰减为0.000 5,学习率为0.001,动量参数为0.9,学习率衰减参数为0.1,区域建议网络层生成的候选检测框的数量为300,IoU为0.7,网络迭代次数为50 000次。采用测试集上的平均精度值(average precision,AP)作为Faster R-CNN的精度指标,结果见表 1,检测结果的总体精度(mean average precision,mAP)为各建筑风格类别检测的AP值的算术平均值。IoU设定为0.7时的测试集检测结果示例见图 12图 12中记录了建筑风格类型和检测概率值。

    表  1  各类风格建筑区域的检测精度表
    Table  1.  Detection Precision of Architectural Area of Different Styles on Test Set
    建筑区域的风格类别 AP/%
    战国时期楚国建筑风格 69.23
    汉代建筑风格 57.17
    唐代建筑风格 89.43
    宋代建筑风格 86.73
    元代建筑风格 73.61
    明代建筑风格 78.53
    清代建筑风格 79.59
    京派民居 88.53
    苏派民居 72.19
    徽派民居 84.37
    民国民居 78.61
    现代建筑风格 84.13
    古希腊建筑风格 76.97
    古罗马建筑风格 73.37
    哥特式建筑风格 63.17
    法国古典风格 43.85
    巴洛克建筑风格 48.29
    拜占庭建筑风格 89.20
    其他西式风格 65.40
    mAP 73.81
    下载: 导出CSV 
    | 显示表格
    图  12  测试集的建筑区域检测结果示例
    Figure  12.  Selected Examples of Architectural Area Detection Results on Test Set

    支持向量机选取高斯核函数。通过网格搜索确定模型最佳参数,经过十折交叉验证得到测试集上的AP为92%。最佳参数值包括核函数次数为1.0,Gamma为0.001,最大迭代次数为1 000。为真实反映匹配方法的精度,随机选取100个相邻街景像对,利用Faster R-CNN检测建筑区域,得到808个建筑区域匹配对,经人工核验,其中有206对真实同名建筑区域,如图 13所示。利用相同的检验数据,将本文方法与仅采用附加核线约束的SURF匹配结果的精确率、召回率和F1分数[23]进行对比,结果如表 2所示。

    图  13  同名建筑区域匹配结果的混淆矩阵
    Figure  13.  Confusion Matrix of Matching Results with the Same Name Architectural Area
    表  2  两种同名建筑区域匹配方法的精度对比结果
    Table  2.  Accuracy Comparison of Two Matching Methods
    同名建筑区域匹配方法 精确率/% 召回率/% F1分数 耗时/ms
    本文方法 86.1 90.3 0.882 180.1
    核线约束下的SURF匹配 66.7 94.2 0.781 241.6
    下载: 导出CSV 
    | 显示表格

    表 2可知,本文方法的F1分数相比SURF算法提高10.1%,匹配一个建筑区域对的平均耗时减少25.4%。召回率下降说明PPSVM在筛选候选匹配对时,存在将同名建筑区域误判成非同名建筑区域的情况,而核线约束下的SURF匹配能通过定量的核线距离关系辅助判别两建筑区域是否同名,在一定程度上减少对影像特征的依赖,但核线约束不具有唯一性,精确率降低。因此,融合影像特征和核线几何约束可以提高建筑区域匹配的精确率和时间效率。图 14展示了相邻两张街景进行同名建筑区域匹配的过程,左、右两张街景记录了建筑风格类型和检测概率值。

    图  14  相邻两张街景的同名建筑区域匹配过程
    Figure  14.  Matching Process of Two Corresponding Architectural Areas in an Image Pair

    为检验两种位置映射方法的准确率,实验选取150个不同的单体建筑区域影像,分别通过本文的两种方法映射至建筑俯视轮廓。结合街景影像,人工统计正确映射的单体建筑数量,得到表 3的结果。从表 3可以看出,对于单像方位映射方法,单体建筑区域与对应俯视轮廓的IoU均值为0.682,匹配正确率达80.3%,单个建筑区域完成方位映射的平均时间为0.017 s;双像点位映射得到的单体建筑位置点与对应俯视轮廓的匹配距离均值为5.498 m,匹配正确率为85.1%,一对同名建筑区域完成点位映射的平均时间为2.741 s。两种方法的时间效率对比如图 15所示。

    表  3  两种位置映射方法准确率的对比结果
    Table  3.  Accuracy Results of Two Mapping Methods
    位置映射方法 正确映射/个 错误映射/个 映射失败/个
    单像方位映射 106 26 18
    双像点位映射 63 11 76
    下载: 导出CSV 
    | 显示表格
    图  15  两种映射方法的耗时对比
    Figure  15.  Time Consumption Comparison of Two Mapping Methods

    表 3图 15可以看出,双像点位映射方法的正确率较高,但满足映射条件的街景像对的比例仅为49.33%,丢失了较多不同风格的建筑区域影像;同时,双像点位映射方法涉及同名建筑区域匹配和相对定向,计算耗时长,不适合处理海量街景数据。单像方位映射方法在保证一定精度的同时,成功率高,耗时短,但不适合处理过小的建筑区域。因此,在实际处理大量街景数据时,一般采用单像方位映射匹配建筑轮廓。当检测到过小的单体建筑区域时,可采用双像点位映射进行辅助定位匹配。

    表 4显示了在同一个区域两种映射方法导致多映射的单体建筑数量。在两种映射方法下,存在多映射的建筑数量占比分别为42.9%和23.4%。可见,单像方位映射出现多映射的比例高于双像点位映射,主要因为单像方位映射的成功率高,更多单体建筑影像匹配到建筑轮廓,导致出现较多的多匹配现象,这也进一步说明了建筑风格综合判定方法的必要性。

    表  4  存在多映射问题的单体建筑数量
    Table  4.  Number of Buildings with Multiple Mapping Problems
    位置映射方法 成功映射的建筑/个 存在多映射的建筑/个
    单像方位映射 13 522 5 805
    双像点位映射 9 595 2 245
    下载: 导出CSV 
    | 显示表格

    图 16~19为实验生成的北京、西安、上海、武汉4个城市的建筑风格地图。为检验地图的准确率,本文在生成结果中选取了331个单体建筑轮廓,人工标注其建筑风格特征作为先验信息。将实验结果与先验类别进行验证,得到图 20所示的建筑风格地图分类精度图。

    图  16  北京市建筑风格地图
    Figure  16.  Architectural Style Map of Beijing
    图  17  西安市建筑风格地图
    Figure  17.  Architectural Style Map of Xi'an
    图  18  上海市建筑风格地图
    Figure  18.  Architectural Style Map of Shanghai
    图  19  武汉市建筑风格地图
    Figure  19.  Architectural Style Map of Wuhan
    图  20  建筑风格地图分类精度条形图
    Figure  20.  Bar Chart of Classification Results of Architectural Styles

    由于实验区域中没有用于验证的元代、宋代以及法国古典风格的先验建筑轮廓信息,在忽略这3类建筑风格的情况下,建筑轮廓分类精确率的平均值为55.1%,召回率的平均值为76.4%,F1得分均值为0.601。结果表明,单体建筑轮廓的分类平均精确率不高,但能基本反映建筑风格的地理分布。此外,楚国风格、汉代风格、明代风格、民国民居风格、拜占庭风格的精确率远低于平均值。结合Faster R-CNN模型测试结果,分类精确率低的建筑风格对应的训练样本数量少。从建筑风格特征上考量,一方面,同一类建筑风格在不同地域有着内在差异,例如上海和武汉的民国民居在建筑立面外观和建筑结构上存在不同。另一方面,不同类的建筑风格有相似性,在建筑结构或构件上有相同或相似的元素,例如民国民居往往会采用欧式建筑立面,古典主义的欧式建筑会采用相同的窗形或柱式。上述两个因素会导致建筑风格的影像分类较困难,影响建筑风格地图的分类准确率。

    本文以城市街景影像为研究对象,在完成19类建筑风格影像提取的基础上,结合球形全景影像的空间几何与图像特征,探究同名建筑区域的匹配方法,建立单、双像建筑区域的位置映射关系,实现街景建筑区域影像与单体建筑俯视轮廓的匹配,实际生成了4座城市中心区域的建筑风格地图。在未来仍需进一步提高风格分类精度和地图的精细度。

  • 图  1   地下工井模型拓扑重构算法流程

    Figure  1.   Algorithm Flowchart of the Topological Reconstruction of Underground Cable Well Model

    图  2   顶面点云数据及边界线关键点

    Figure  2.   Top Point Cloud Data and Boundary Key Points

    图  3   工井模型及不同类型的墙壁模型

    Figure  3.   Working Well Model and Different Types of Wall Models

    图  4   平面墙壁关键点

    Figure  4.   Flat Wall Key Points

    图  5   曲面墙壁关键点

    Figure  5.   Curve Surface Wall Key Points

    图  6   模型拓扑重构流程

    Figure  6.   Flowchart of the Topological Reconstruction of the Model

    图  7   相交三角形交线示意图

    Figure  7.   The Intersection of Two Triangles

    图  8   墙面拓扑重构

    Figure  8.   Topological Reconstruction of Wall Surface

    图  9   地下电缆工井模型拓扑结构

    Figure  9.   Topological Structure of Underground Cable Working Well Model

    图  10   地下管网层次关系图

    Figure  10.   Underground Pipe Network Class Diagram

    图  11   软件建模效果图

    Figure  11.   Modeling Results of Software

    图  12   工井模型(部分)以及通过孔洞连接的电缆

    Figure  12.   Well Model (Part) and the Cable Connected Through the Holes

    图  13   工井及电缆拓扑关系查询

    Figure  13.   Topology Relation Query of Well and Cable

    图  14   对具有孔洞的面点云进行多种方法墙壁模型重建

    Figure  14.   Different Methods for the Reconstruction of Wall Models with Holes Depands on Point Cloud Data

  • [1] 王舒, 宁芊.地下管线空间数据模型及三维可视化[J].软件导刊, 2015, 14(2):78-80 doi: 10.3969/j.issn.1672-7800.2015.02.029

    Wang Shu, Ning Qian. Underground Pipeline Spatial Data Model and 3D Visualization[J]. Software Guide, 2015, 14(2):78-80 doi: 10.3969/j.issn.1672-7800.2015.02.029

    [2] 危双丰, 刘明蕾, 赵江洪, 等.利用点云检测室内导航元素的方法综述[J].武汉大学学报·信息科学版, 2018, 43(12): 2003-2011 http://ch.whu.edu.cn/CN/abstract/abstract6284.shtml

    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): 2003-2011 http://ch.whu.edu.cn/CN/abstract/abstract6284.shtml

    [3] 卢丹丹, 谭仁春, 郭明武, 等.城市地下管线三维建模关键技术研究[J].测绘通报, 2017(5):117-119 http://d.old.wanfangdata.com.cn/Periodical/chtb201705027

    Lu Dandan, Tan Renchun, Guo Mingwu, et al. Research on the Key Technology of Urban Underground Pipeline Three-dimensional Modeling[J]. Bulletin of Surveying and Mapping, 2017(5):117-119 http://d.old.wanfangdata.com.cn/Periodical/chtb201705027

    [4] 毕天平, 孙立双, 钱施光.城市地下管网三维整体自动建模方法[J].地下空间与工程学报, 2013, 9(s1):1473-1476 http://d.old.wanfangdata.com.cn/Periodical/dxkj2013z1002

    Bi Tianping, Sun Lishuang, Qian Shiguang. Automatic 3D Modeling Method for Urban Underground Pipe Network[J]. Chinese Journal of Underground Space and Engineering, 2013, 9(s1):1473-1476 http://d.old.wanfangdata.com.cn/Periodical/dxkj2013z1002

    [5] 李清泉, 严勇, 杨必胜, 等.地下管线的三维可视化研究[J].武汉大学学报·信息科学版, 2003, 28(3):277-282 http://ch.whu.edu.cn/CN/abstract/abstract4765.shtml

    Li Qingquan, Yan Yong, Yang Bisheng, et al. Research on 3D Visualization of Underground Pipeline[J]. Geomatics and Information Science of Wuhan University, 2003, 28(3):277-282 http://ch.whu.edu.cn/CN/abstract/abstract4765.shtml

    [6] 周京春.地下管网三维空间数据模型及自动化精细建模方法研究[D].武汉: 武汉大学, 2016 http://cdmd.cnki.com.cn/Article/CDMD-10486-1016113516.htm

    Zhou Jingchun. Research on 3D Spatial Data Models and Refined 3D Automatic Modeling Methods for Underground Pipeline[D]. Wuhan: Wuhan University, 2016 http://cdmd.cnki.com.cn/Article/CDMD-10486-1016113516.htm

    [7] 刘军, 钱海峰, 孙永新.基于Skyline的三维综合地下管线应用与研究[J].城市勘测, 2011(4):43-45 doi: 10.3969/j.issn.1672-8262.2011.04.013

    Liu Jun, Qian Haifeng, Sun Yongxin. Application of Underground Pipeline Three-dimensional Modeling Based on Skyline[J]. Urban Geotechnical Investigation and Surveying, 2011(4):43-45 doi: 10.3969/j.issn.1672-8262.2011.04.013

    [8] 吴思, 杨艳梅, 王明洋, 等.一种真三维地下管线井室自动建模方法[J].测绘科学技术学报, 2016, 33(4):400-404 doi: 10.3969/j.issn.1673-6338.2016.04.014

    Wu Si, Yang Yanmei, Wang Mingyang, et al. An Automatic Modeling Method of True Three-Dimensional Well Chamber of Underground Pipeline[J]. Journal of Geomatics Science and Technology, 2016, 33(4):400-404 doi: 10.3969/j.issn.1673-6338.2016.04.014

    [9] 钟远根, 戴相喜, 李颖捷, 等.三维地下管线建模及系统实现研究[J].现代测绘, 2014, 37(1):25-27 doi: 10.3969/j.issn.1672-4097.2014.01.008

    Zhong Yuangen, Dai Xiangxi, Li Yingjie, et al. Study on Three-Dimensional Underground Pipeline Modeling and System Implementation[J]. Modern Surveying and Mapping, 2014, 37(1):25-27 doi: 10.3969/j.issn.1672-4097.2014.01.008

    [10]

    Edelsbrunner H. Surface Reconstruction by Wrapping Finite Sets in Space[M]//Discrete and Computational Geometry. Berlin, Heidelberg: Springer, 2003: 379-404

    [11]

    Boissonnat J D, Cazals F. Smooth Surface Reconstruction via Natural Neighbour Interpolation of Distance Functions[J]. Computational Geometry Theory and Applications, 2002, 22(1):185-203 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=c9ae433dee46793b470fa75878a43437

    [12]

    Vosselman G. Fusion of Laser Scanning Data, Maps, and Aerial Photographs for Building Reconstruction[C]. IEEE Geoscience and Remote Sensing Symposium, Toronto, Canada, 2002

    [13] 贺彪, 李霖, 郭仁忠, 等.顾及外拓扑的异构建筑三维拓扑重建[J].武汉大学学报·信息科学版, 2011, 36(5):579-583 http://ch.whu.edu.cn/CN/abstract/abstract545.shtml

    He Biao, Li Lin, Guo Renzhong, et al. 3D Topological Reconstruction of Heterogeneous Building Considering Exterior Topology[J]. Geomatics and Information Science of Wuhan University, 2011, 36(5):579-583 http://ch.whu.edu.cn/CN/abstract/abstract545.shtml

    [14] 朱冒冒, 冯祖军, 郭建华, 等.三维栅格法在数据点云精简中的应用[J].上海大学学报(自然科学版), 2005, 11(3):242-246 doi: 10.3969/j.issn.1007-2861.2005.03.006

    Zhu Maomao, Feng Zujun, Guo Jianhua, et al. Data Point Reduction Using a 3-Dimensional Cube Algorithm[J]. Journal of Shanghai University (Natural Science), 2005, 11(3):242-246 doi: 10.3969/j.issn.1007-2861.2005.03.006

    [15] 周波, 陈银刚, 顾泽元.基于八叉树网格的点云数据精简方法研究[J].现代制造工程, 2008(3):64-67 doi: 10.3969/j.issn.1671-3133.2008.03.019

    Zhou Bo, Chen Yingang, Gu Zeyuan. Data Point Reduction on Octree Cube Algorithm[J]. Modern Manufacturing Engineering, 2008(3):64-67 doi: 10.3969/j.issn.1671-3133.2008.03.019

    [16]

    Schnabel R, Wahl R, Klein R. Efficient RANSAC for Point-Cloud Shape Detection[J]. Computer Graphics Forum, 2007, 26(2):214-226 doi: 10.1111/j.1467-8659.2007.01016.x

    [17]

    Schnabel R, Wahl R, Klein R. Shape Detection in Point Clouds[OL]. https://cg.cs.uni-bonn.de/aigaion2root/attachments/cg-2006-2.pdf, 2006

    [18]

    Schnabel R, Wahl R, Wessel R, et al. Shape Recognition in 3D Point Clouds[J]. Václav Skala-Union Agency, 2008, 272(1):65-72 http://cn.bing.com/academic/profile?id=2305f57670d8f2369dc1a37532f0344f&encoded=0&v=paper_preview&mkt=zh-cn

    [19] 刘如飞, 卢秀山, 岳国伟, 等.一种车载激光点云数据中道路自动提取方法[J].武汉大学学报·信息科学版, 2017, 42(2):250-256 http://ch.whu.edu.cn/CN/abstract/abstract5669.shtml

    Liu Rufei, Lu Xiushan, Yue Guowei, et al. An Automatic Extraction Method of Road from Vehicle-Borne Laser Scanning Point Clouds[J]. Geomatics and Information Science of Wuhan University, 2017, 42(2):250-256 http://ch.whu.edu.cn/CN/abstract/abstract5669.shtml

    [20] 于庆, 王结臣, 陈焱明.一种三角网边界提取的方法[J].测绘科学, 2009(s2):82-83 http://www.cnki.com.cn/Article/CJFDTotal-CHKD2009S2031.htm

    Yu Qing, Wang Jiechen, Chen Yanming. A Method of Boundary Extraction for Trianglemesh [J]. Science of Surveying and Mapping, 2009(s2):82-83 http://www.cnki.com.cn/Article/CJFDTotal-CHKD2009S2031.htm

    [21] 谷金良. B样条边界面法及边界积分方程中的等几何方法研究[D].长沙: 湖南大学, 2012

    Gu Jinliang. Research of the B-spline Boundary Face Method and the Isogeometric Analysis in the Boundary Integral Equation[D].Changsha: Hunan University, 2012

    [22] 熊运阳. CAD/CAE中样条曲线曲面的研究[D].杭州: 浙江大学, 2014 http://cdmd.cnki.com.cn/Article/CDMD-10335-1014172903.htm

    Xiong Yunyang. Study on Spline Curves and Surfaces in CAD/CAE[D].Hangzhou: Zhejiang University, 2014 http://cdmd.cnki.com.cn/Article/CDMD-10335-1014172903.htm

    [23]

    Hooton J, Jones M H, Shur J, et al. A Method for the Selection of Algorithms for Form Characterization of Nominally Spherical Surfaces[J]. Precision Engineering, 1999(25):39-56 http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=471e0e774aea0100602ff3207fab5d4a

    [24] 朱春晓, 黄明, 倪春迪.三维CSG-BREP拓扑模型的研究[J].测绘工程, 2017, 26(8):20-23 http://d.old.wanfangdata.com.cn/Periodical/chgc201708005

    Zhu Chunxiao, Huang Ming, Ni Chundi. Research on the Topological Model of Three-Dimisional CSG-BREP[J]. Engineering of Surveying and Mapping, 2017, 26(8):20-23 http://d.old.wanfangdata.com.cn/Periodical/chgc201708005

    [25] 陈波.基于构造立体几何(CSG)法的三维重建技术的研究[D].苏州: 苏州大学, 2013 http://cdmd.cnki.com.cn/Article/CDMD-10285-1013231577.htm

    Chen Bo. Research on 3D Solid Reconstruction Based on Constructive Solid Geometry Method[D]. Suzhou: Soochow University, 2013 http://cdmd.cnki.com.cn/Article/CDMD-10285-1013231577.htm

    [26] 王寒冰.基于特征的BREP→CSG模型转换方法及其应用[D].合肥: 合肥工业大学, 2013

    Wang Hanbing. Feature Based BREP→CSG Model Conversion Method and Its Application[D]. Hefei: Hefei University of Technology, 2013

    [27] 赵江洪, 马思宇, 黄明, 等.地下工井电缆网络模型及拓扑关系构建[J].测绘科学, 2019, 44(3):142-148 http://d.old.wanfangdata.com.cn/Periodical/chkx201903023

    Zhao Jianghong, Ma Siyu, Huang Ming, et al. The Construction of Network Model of Underground Cable Working Well and Its Topological Relationship[J]. Science of Surveying and Mapping, 2019, 44(3):142-148 http://d.old.wanfangdata.com.cn/Periodical/chkx201903023

    [28] 贾军辉, 黄明, 刘祥磊.基于三维狄洛尼三角网的曲面重建算法[J].测绘学报, 2018, 47(2):281-290 http://d.old.wanfangdata.com.cn/Periodical/chxb201802017

    Jia Junhui, Huang Ming, Liu Xianglei.Surface Reconstruction Algorithm Based on 3D Delaunay Triangulation [J]. Acta Geodaetica et Cartographica Sinica, 2018, 47(2):281-290 http://d.old.wanfangdata.com.cn/Periodical/chxb201802017

    [29]

    Gilbert B.3D Boolean/CSG Library[OL]. https://github.com/gilbo/cork, 2016

    [30]

    Gilbert B. Putting Holes in Holey Geometry: Topology Change for Arbitrary Surfaces[OL]. https://github.com/gilbo/toptop, 2013

  • 期刊类型引用(1)

    1. Jianjun ZHU,Leyang WANG,Jun HU,Bofeng LI,Haiqiang FU,Yibin YAO. Recent Advances in the Geodesy Data Processing. Journal of Geodesy and Geoinformation Science. 2023(03): 33-45 . 必应学术

    其他类型引用(2)

图(14)
计量
  • 文章访问数:  1199
  • HTML全文浏览量:  204
  • PDF下载量:  110
  • 被引次数: 3
出版历程
  • 收稿日期:  2019-05-27
  • 发布日期:  2019-12-04

目录

/

返回文章
返回