opencv判断多边形的邻接关系_基于拓扑关系的BIM室内封闭空间边界建筑要素搜索方法研究...

作 者 信 息

罗　丰1，王丽园1，李　霖2，杨　晶1

(1. 中交第二公路勘察设计研究院有限公司，湖北　武汉　430056；2. 武汉大学 资源与环境科学学院，湖北　武汉　430079)

“

【摘要】BIM作为工程建设领域的热门新兴技术，包含建筑物全生命周期内的各种数字化信息，但其对要素间空间关系的描述限于建筑领域的专业性表达，以及建模过程中的预定义。基于BIM数据的几何-语义特点，提出一种自动化构建BIM建筑要素拓扑关系的方法，并以此为基础搜索室内封闭空间的边界建筑要素，得到以建筑要素为单元的室内空间划分；以整层建筑物空间作为实验数据，通过实验结果验证了算法的准确性，实现BIM室内三维空间分割，为基于BIM的空间分析提供数据支撑。

【关键词】BIM；封闭空间；边界建筑要素；拓扑关系

【中图分类号】TU1 【文献标识码】A 【文章编号】1672-1586(2019)02-0104-06

”

引文格式：罗　丰，王丽园，李　霖，等. 基于拓扑关系的BIM室内封闭空间边界建筑要素搜索方法研究[J].地理信息世界,2019,26(2):104-109.

正文

0　引　言

随着智慧城市的建设，三维建筑物建模在城市规划、建设和管理等领域正成为研究热点。建筑信息模型(BIM)被认为是建筑、工程和管理领域最具有前景的工具，在诸多领域都有所应用，如能源建模、设施管理和建筑设计等。BIM包含精确的几何信息、建筑相关的语义属性数据以及对建筑要素间关系的专业性表达，以支撑建筑物生命周期内的设计、采购、制造和运营等活动。通过BIM技术，不仅能够通过数字化方法构建精细的建筑物模型，还能实现对整个项目周期内各类建筑信息的动态掌控。

构建和维护空间元素的拓扑关系是GIS的一个重要特征，拓扑关系对于空间操作和信息查询具有重要意义，因此基于BIM应用的空间分析研究有其必要性。最新的研究尝试将BIM与GIS相结合来尝试解决不同领域三维对象的描述问题，如将BIM应用于三维地籍，BIM与GIS系统集成，BIM与GIS标准的映射，结合本体解决BIM和GIS的异质性问题。在基于BIM的空间分析研究中，有相关学者基于经典拓扑数据模型，进行建筑要素拓扑关系推理，结合几何-语义特性的空间查询、空间信息提取以及室内导航等研究，但从三维角度基于BIM数据进行拓扑关系构建，以实现空间分析的研究较少。

现有的BIM软件和工具无法解析建筑模型中隐含的几何和拓扑信息，也不具备系统构建建筑要素间拓扑关系的功能。如果在BIM建模过程中未定义各要素间的空间关系，则难以实现基于要素间关系的空间分析。因此，本文提出一种基于拓扑关系的BIM室内封闭空间边界要素搜索方法，通过构建建筑要素间的拓扑关系，得到构成室内封闭空间边界的建筑要素，以建筑要素为单元从三维的角度实现基于建筑要素的室内空间划分。

1　BIM标准

BIM作为工程建设领域近年来的新兴热门技术，已经逐步应用于建筑、交通运输、测绘、制造等多个行业，服务于项目勘察、设计、施工、运营、维护、拆除等各个阶段。工业基础类(IFC)作为国际通用的BIM标准，定义了建筑物要素的类别、属性、几何表达方式、拓扑数据结构等，是BIM在工程应用的主要数据格式。

BIM通过IFC标准化接口定义建筑物内的建筑要素类别，其中表示建筑物空间范围的要素包括IfcBuilding、IfcBuildingStorey、IfcSpace等，具体类别见表1。通过其语义属性来表达该建筑空间的功能特性，如IfcSpace的实例对象中标注其功能属性为廊道、房间等。在IFC标准中，IfcSpace定义为表示建筑物内部空间范围的要素，但仅表示一个三维空间范围，并未指明该空间范围的边界由哪些建筑要素构成，建筑要素之间的关系如何。在实际应用中建筑要素之间的关系具有重要作用，如在室内导航中，门、窗等开口建筑要素的位置和方向决定了整个室内空间的连通方式，通过构建门、窗等要素与墙体的空间关系可以解决此类问题。此外，在BIM实例中IfcSpace不是模型的必备要素类型，即BIM数据可能并未定义IfcSpace对象。

表1　部分IFC标准建筑要素类别

Tab.1 Partial classification of building elements in IFC standard

关于要素关系的表达，IFC标准以逆关系(Inverse relations)描述要素之间的多种连接关系，如图1所示，标签ConnectedTo ContainedInStructure隐含地表达相邻和包含关系两种空间关系。IFC标准通过相关模块能够表达建筑物内各对象之间的关系，但其基于建筑领域专业性的描述方式，与自然查询语言中的相邻、相离、包含等关系不同。在三维城市模型中，李霖等学者提出一种面向几何基元和语义对象间的拓扑关系分类，以相接、相等、覆盖、包含、相离5种关系描述三维城市对象间的拓扑关系。本文以此为基础，通过提出一种自动化拓扑关系构建方法，对建筑要素间的拓扑关系进行描述，并以此为基础得到建筑物内封闭空间的边界构成，实现三维空间划分。

图1　IFC标准中逆关系示例

Fig.1 Example of inverse relations in IFC standard

2　数据预处理

IFC建筑物的几何构成主要是三角面，在建模时采用三角面的优势是利于三维可视化。本文提出的拓扑关系自动化构建方法，通过公共元素进行对象之间的关联是核心思想，因此算法实现的前提是三角面之间仅在端点相接，否则会导致算法在搜索关联面时不准确。如图2所示，a图中两个平面f1和f2 相接于边界，那么通过关联面搜索可以唯一化公共边；b图中3个三角面相接，面f3和f4、f5分别相接于边界线AB 和BC，此时根据f3的边界线AC 搜索关联面时，面f4和f5无法搜索到。

因此，在构建建筑要素的拓扑关系时，应先对所有的边界面进行预处理。首先，验证建筑要素内几何基元自身的有效性；其次，当出现几何基元之间相交的情况时，打散所有相互之间有公共部分的几何基元。在图2b中，面f3的边界线AC 与面f4的边界线AB 相交，同时与面f5的边界线BC 相交，那么在预处理过程中将面f3的边界线AC 打断成AB 和BC 两部分。最后，记录打散前后新旧几何基元之间的关系，即AC →{AB ,BC }。经过预处理，面f3的边界边由3条变为4条 ，此时通过边界边搜索关联面时，边界线AB 搜索到的关联面为f4，边界线BC 搜索到的关联面为f5。对应到打断之前的几何基元之间的关系，面f3与面f4、f5相接于公共边界。

图2　多边形在边界相交的两种情形

Fig.2 Two cases: polygon-polygon intersect at the boundary

3　封闭空间边界搜索

IFC标准定义的建筑要素种类繁多，本文仅以构成封闭空间的常见边界建筑要素作为对象进行讨论，包括墙体(IfcWall)、立柱(IfcColume)、楼板(IfcSlab)、门(IfcDoor)、窗(IfcWindow)等。其中墙体和立柱为竖向建筑要素，楼板为水平建筑要素，门和窗作为开口附着于竖向建筑要素上，因此在搜索封闭空间边界时，只考虑竖向建筑要素间的拓扑关系。封闭空间边界建筑要素自动搜索方法，其核心在于几何-语义角度的拓扑关系构建，以及基于拓扑关系的边界要素搜索。在几何方面，拓扑关系通过对单个建筑要素的几何基元构建拓扑关系，可验证建筑要素自身的有效性；从语义方面，以单个建筑要素自身拓扑的有效性和完备性为基础，通过布尔运算的交集判断多建筑要素间的拓扑关系；最后结合建筑要素间拓扑关系和二维闭合多边形边界搜索算法，得到BIM室内三维空间划分。整个算法从数据预处理、拓扑关系构建到边界要素搜索是一个连贯的过程，基于现有BIM数据无需人工干预即可实现自动化处理，得到建筑要素间的拓扑关系以及封闭空间边界的组成要素。

3.1 几何基元拓扑关系构建

结合IFC封闭空间边界组成要素和几何数据组织形式，基于边界公共元素搜索最邻近元素，进行单个建筑要素内部的几何基元拓扑关系构建，具体方法包括：

1)提取建筑物封闭空间边界构成要素，以离散面集合表示其几何信息；

2)判断单个边界要素内的离散面是否存在非顶点相交情况，若存在非顶点相交，则通过线面相交、线线相交运算将线、面剖分，确保同一边界要素的几何元素仅在顶点处相交；

3)以离散面集合为对象，通过自动搜索封闭几何体算法构建每个边界要素内部几何基元的拓扑关系，得到点、线、面的拓扑关系集合，包括以下步骤：

①寻找离散面集合的最外离散面。利用计算几何算法库(CGAL)的三维凸壳算法得到离散面集合的凸壳，计算离散面集合中所有离散面到凸壳几何中心的距离，找到距离最远的离散面f最外；

②根据f最外的边界线搜索集合中的所有面，从f最外边界上的任意一条边e出发，搜索共用边e 的其他所有邻接面；

③以步骤②搜索到的第一个面为起始面，计算具有公共边的两个面之间的夹角，利用最大夹角或最小夹角原则寻找最邻近面；

④循环利用新搜索到的最邻近面重复步骤③，直至没有新的面被搜索到。整个离散面拓扑关系构建算法的终止条件是所有的离散面都被搜索过一次，并且公共面被搜索过两次。

3.2 建筑要素拓扑关系构建

构建建筑要素间的拓扑关系，实质就是对不同建筑要素的几何基元进行处理。以两个建筑要素为例，建筑要素间拓扑关系构建的方法流程归纳如下：

1)进行空间运算，得到两个建筑要素的交集。对两个建筑要素的几何基元(体集合或者面集合)进行布尔交运算。如果交集为空，那么建筑要素间的拓扑关系为相离，如果交集不为空，那么建筑要素间的拓扑关系为相接。

2)若要素间的拓扑关系为相接，则进行几何体集合与几何体集合的拓扑关系构建：

①提取几何体的边界点，对处于容差范围内的点进行处理，保证数据的正确性和一致性；

②去除重复的几何元素，使相同的几何元素唯一化。如建筑要素A包含边界面f，建筑要素B同时也包含边界面f，那么公共边界面f应该唯一化；

③进行线线判交和面面判交，根据交集打断具有公共部分的几何元素，并记录打断操作前后新旧几何元素的关系；

④根据打断产生的新旧几何元素更新原始几何元素所在的建筑要素。

3.3 封闭空间边界搜索

基于BIM的几何结构组织特点，如果直接以三维方式搜索封闭空间边界，需要进行大量而且复杂的几何运算。因此，我们采用将竖向边界要素(墙体、立柱等)投影到二维平面，以二维线段搜索封闭多边形的方法来简化这个问题，并结合竖向建筑要素与楼板之间的空间关系得到封闭空间。在封闭空间搜索算法中，以二维多边形表示三维建筑要素，二维多边形的中轴线对应三维空间的墙体。因此，二维中轴线之间的拓扑关系在搜索封闭多边形中起着重要作用。以几何-建筑要素拓扑关系构建的有效性为前提，封闭空间边界搜索算法步骤归纳如下：

1)将竖向三维建筑要素的边界点投影到水平面xoy，连接各投影边界点，得到二维投影多边形P投影；

2)根据P投影的几何形态提取中轴线，若中轴线包含多条线段，则以线段集合表示该投影多边形的中轴线，且应确保这些线段仅相交于端点，记录线段集合与投影多边形P投影的对应关系；

3)根据建筑要素拓扑关系检验其投影多边形中轴线是否相交。若不相交，则分别延长中轴线的起始点直至相交于一点；若相交于非端点，则在交点处对线段进行打断处理，确保所有二维线段只在端点处相交，打断之后记录新旧线段的关系；

4)以任意一条中轴线线段为起始线段，以线段的起点为原点，按照逆时针(或顺时针)方向，根据与起始线段的最大角(或最小角)搜索最邻近的线段；

5)将新搜索到的线段作为起始线段重复前一步骤，直至没有新的关联线段被搜索到，此时封闭多边形搜索完成；

6)根据中轴线、多边形、三维建筑要素之间的对应关系得到各建筑要素与室三维封闭空间的构成关系，二维离散线段围成的封闭多边形对应建筑要素围成的封闭空间。

4 案例分析

本文的IFC实验数据是一栋四层建筑物顶楼，如图3所示，直观上可以分辨出卧室、卫生间、大厅、阳台等，该楼层包含多种建筑要素，包括墙体、窗户、门、栏杆、家具等。栏杆和家具几何结构复杂，并且不属于封闭空间的边界要素，因此我们在实验中将这类建筑要素删除。此外，实验数据的竖向建筑要素不包含立柱，经过几何信息提取后的数据如图4所示。实验的IFC模型几何数据都是三角面，因此在构建拓扑关系之前要检验是否存在非顶点处相交的情况。在确保几何-拓扑有效性的前提下，对实验数据进行拓扑关系构建。

图3　IFC实验数据三维可视化

Fig.3 Visualization of experimental IFC data

图4　实验数据建筑要素提取结果

Fig.4 The extraction result of building elements in experiment data

需要注意的是，实验数据中出现一种跨度多个房间的长墙体，如图5所示，它同时作为多个封闭空间的边界，并且每个封闭空间的边界墙体只是它的子集。在封闭空间边界搜索时应将二维投影多边形的中轴线在交点处剖分，以线段集合表示该长墙体的中轴线。待二维闭合多边形搜索完毕后，再根据对应关系得到封闭空间的边界要素。本实验通过对单个建筑要素进行拓扑关系构建，得到点、线、面间的拓扑关系集合，结合布尔运算结果判断要素间的拓扑关系。根据实验结果得到墙体间的拓扑关系为相接和相离两种，门、窗户与对应墙体的拓扑关系也为相接。

图5　跨度多个房间的长墙体示例

Fig.5 Instance of long wall which spans multiple rooms

对实验数据进行封闭空间边界搜索，整个楼层墙体的二维化结果如图6所示。整个楼层空间由18个竖直规则墙体围成，各线段表示三维建筑要素在二维平面投影多边形的中轴线，红色线段表示原始IFC数据中跨度多个封闭空间的长墙体，黑色线段表示只作为完整对象参与构建封闭空间边界的墙体，黑色节点表示投影多边形中轴线的交点。图5所示的长墙体W 1在交点处剖分成3个子线段l 1、l 2和l 3，即W 1←→{l 1,l 2,l 3}。根据搜索算法得到l 2，l 4，l 5，l 6，l 7围成二维封闭多边形，由对应关系可知W 1、W 2、W 3和W 4构成一个封闭空间。如图6所示，案例数据由于原始几何建模方式的原因，搜索出的封闭空间大多数都包含作为多个封闭空间边界的长墙体。本实验数据通过自动搜索得到8个封闭空间，每个封闭空间的边界墙组成情况见表2，其中封闭空间No.3的边界墙体要素数量多达12个，几何形态复杂。

图6　IFC实验数据搜索二维封闭多边形结果

Fig.6 Result of searching two-dimensional closed polygon based on experiment IFC data

表2　楼层IFC数据封闭空间搜索结果

Tab.2 Results of searching boundary building elements in a one-storey IFC data

5　结束语

本文基于BIM几何-语义数据特点，提出一种通过构建建筑要素间拓扑关系，实现自动化搜索BIM数据室内封闭空间的边界建筑要素的方法，从而对整个BIM内部空间进行三维划分。通过对整层建筑物数据开展实验，验证了算法的有效性和完备性，可为基于室内空间划分的分析应用提供数据支撑。本研究提取的建筑要素类别有限且几何结构规则，并且楼层间的空间关联尚未建立，在后续的研究中，将针对这些问题开展进一步研究以提高方法的适用性。

本期回顾

理论研究

 · 基于地理空间信息的文化遗产可持续发展指标建设  ·基于Landsat 8的成都市城市热岛时空格局变化研究  ·基于视频场景多因素的人群运动状态分析  ·岩体边坡影像的控制测量方法研究  ·基于套件技术的一站式GIS系统框架结构研究  · JPEG 2000编码参数对遥感影像压缩质量的影响  · 基于Mann-Kendall法的长三角城市集聚效应空间结构研究  · 3S技术在我国生态环境动态演变研究中的应用进展  · 基于异常误差补偿的多期水准网动态平差  · 基于立体网格的放射性污染物扩散过程模拟与表达   创新应用  · 基于感兴趣点和滴滴数据的打车需求分析  · 新世纪版《国家普通地图集》的内容策划与结构设计  · 基于支持向量机的新闻事件类型识别  · 国土资源执法监察动态巡查系统设计与实现  · 基于AutoCAD的网络地图服务系统   技术探讨  · 一种基于地理格网的三维模型数据统一表达方法  · 基于用户历史瓦片浏览记录的兴趣点智能搜索方法优化研究  · 稀无控制点无人机遥感影像几何定位方法研究  · 贵州省市州域经济差异的时空演变

邮箱变更声明

·《地理信息世界》邮箱变更声明

网站开通公告

·关于开通《地理信息世界》网站的公告

诚聘特约审稿专家

·诚聘|《地理信息世界》诚聘特约审稿专家

专题组稿

·约稿函|《地理信息世界》关于开辟“博士综述论坛”专栏的约稿函