2. 南京师范大学虚拟地理环境教育部重点实验室, 江苏 南京 210023;
3. 南京师范大学地理环境演变国家重点实验室培育点, 江苏 南京 210023;
4. 上海树融数据科技有限公司, 上海 200082;
5. 自然资源部城市国土资源监测与仿真重点实验室, 广东 深圳 518034
2. Key Laboratory of Virtual Geographic Environment, Ministry of Education, Nanjing Normal University, Nanjing 210023, China;
3. State Key Laboratory Cultivation Base of Geographical Environment Evolution, Nanjing Normal University, Nanjing 210023, China;
4. Shanghai Shurong Data Technology Co. Ltd., Shanghai 200082, China;
5. Key Laboratory of Urban Land Resources Monitoring and Simulation, Ministry of Natural Resources, Shenzhen 518034, China
三维城市模型作为数字城市及智慧城市空间框架数据的核心内容[1-3]，被广泛应用于旅游、室内导航、灾害逃生、城市规划、虚拟战场等领域[4-7]。但目前开放的三维城市模型一般是面向可视化或者具体应用，模型中会存在语义、拓扑和几何不一致性的问题，如建筑物几何模型未闭合、建筑墙面为非平面等[8]，从而导致三维模型并不能满足主题查询、空间数据挖掘以及基于空间关系计算、分析等方面的需求，不利于其在下游领域的深层次应用。对三维空间模型数据的几何、拓扑和语义的一致性检测和修复是三维模型空间分析和下游领域深度应用的前提和保证[9-10]，其中三维空间模型数据的拓扑一致性是保证模型几何、拓扑表达的真实性、正确性的重要方面[11-12]。在CAD领域，三维实体通常以边界表示法(boundary representation，B-rep)表示，并以数据结构强制性要求多面体为流形形体[13-14]。在GIS相关领域，对三维模型的检测主要根据模型建模需求定义几何拓扑一致性规则，再以规则检测模型为主[15-16]。文献[12]针对3D城市模型的拓扑一致性检测给出13条规则，从拓扑和数学公理两个方面验证规则的完整性、正确性、有效性。但该规则未考虑多面体的孔洞问题，无法检测CityGML(city geography markup language)模型LOD3中含有门窗的建筑物，模型的拓扑检测不完整。
本文拟针对CityGML模型数据，提出适合多种LODs数据的语义约束下的拓扑一致性规则，设计相应的检测算法，并采用公开数据对算法进行试验与验证。
1 CityGML拓扑一致性规则 1.1 CityGMLCityGML是开放地理空间信息联盟(open geospatial consortium, OGC)为实现三维城市模型信息共享和可视化而提出的国际标准，基于XML模式实现了虚拟三维城市模型中城市空间对象三维几何、拓扑、语义信息的存储和交换[17-18]。作为一种开放、结合室内外场景的数据模型，CityGML被广泛地应用于三维GIS领域，如室内外一体化导航、无人驾驶、防灾减灾、路径分析等[19-20]。CityGML有如下特点：①按照GML3.1.1实现了建筑物、地形、交通、植被、水体、城市设施、土地利用等三维城市专题模块，其中以建筑物模型的使用最为广泛与复杂；②CityGML中空间对象的几何属性通过几何-拓扑模型进行表示[21]，即更高维度的图元拓扑由低维原语构建，简单几何模型一般有4个维度的基本几何元素，每一个维度对应一个几何-拓扑基元，零维(point)、一维(curve)、二维(surface)、三维(solid)，即点、线、面、体；③针对同一个地理空间对象，CityGML在模型的语义层面上，对空间对象的属性或对象之间的关系进行描述，在空间层面上，对空间对象所在的空间位置和范围等进行描述，因此CityGML的拓扑一体化模型通过对应关系相互联系空间对象的几何语义建模和几何拓扑建模两个方面[22](图 1)；④针对同一空间对象的语义、拓扑与几何信息，CityGML模型将描述的精细程度分为5个细节层次LOD0-LOD4，针对建筑物、建筑物部件以及建筑物附属设施由简到繁的表达(图 2)。
1.2 建筑物拓扑结构信息与语义信息的规则定义
CityGML建筑物模型的一致性规则主要有几何拓扑和几何语义两个方面。其几何要素的空间属性通过GML3的几何-拓扑模型进行表示，GML3的几何模型在ISO19107标准上构建。因此根据ISO19107等通用的国际标准提取建筑物模型的几何检测规则，可以在保证建筑物几何有效性的同时，减少不同模型之间数据模型集成的冲突。CityGML的语义属性代表对象在现实空间中的具体意义，因此可根据不同语义对象在现实空间的属性，定义针对建筑物的拓扑规则。
1.2.1 拓扑规则提取通过分析ISO19107/OGC/GML3[10, 16, 23-24]等标准中几何图元点、线、面、体的定义，根据几何维度不同，分层次对几何模型定义相应拓扑规则。其中，参照CityGML OGC标准，surface不考虑曲面，而是由平面片构成，以polygon表示；curve被强制定义为折线，以line string进行表示；solid定义为用封闭表面shell围起来的边界，体有内外边界，即内外shell之分。具体规则定义如下。
(1) 环规则
规则1：linear ring至少由3点组成，当只有3点时，第1点和最后一点坐标不可相同，记为CL-NumPoints&CL-Close。
规则2：除了首末点，linear ring上其他点的坐标值不同，记为CL-DupPoint。
规则3：linear ring上点不能重复，记为CL-DupRing。
规则4：linear ring具有非自相交性，记为CL-SelfInt。如图 3(a)的solid因为面上内外环相交而无效。
|
| 图 3 4个不符合检测规则的无效solid Fig. 3 Four invalid solids which break the validation rules |
(2) 面规则
规则5：polygon上所有点需在同一个平面上，但允许有一定的平面允许容差，记为CP-Plan-Dist&CP-Plan-Dist-TRI。
规则6：同一个面上内外环方向相反，外环方向为右手系逆时针方向，记为CP-Clockwise。
规则7：polygon上的内外环拓扑关系为外环包含内环，记为CP-InROutside。
规则8：面上内环之间不能相交或重叠，记为CP-NoCrossing。
规则9：面内部必须连通可达，记为CP-Connected，如图 3(b)上表面，一个内环因将一个平面分成两部分而无效。
(3) 体规则
规则10：solid至少由4个有效且位于不同平面上的面组成，记为CS-NUMFaces。
规则11：solid上的每一条边被两个面，且仅被两个面所共享，为流形形体。记为CS-Non ManifoldEdge&CS-NonManifoldVer。如图 3(c)，一条公共边被3个面所共享而无效。
规则12：构成solid的各个面不相交，记为CS-SelfInt。
规则13：构成solid上的每个面都有方向，且面的法向量必须都指向solid外部，记为CS-Face-Orient&CS-Face-Out。
规则14：solid必须有连通性，记为CS-Connected。如图 3(c)，solid两部分不可达。
规则15：内shell必须在外shell内部，且shells之间不能相交重叠，只能相切于点、线，记为CS-Inshell-Outside。如图 3(d)，solid因内部shell之间的交叉而无效。
1.2.2 语义属性规则提取基于几何语义一致性原则，CityGML中的面元素、体元素都有对应的语义属性信息，可从语义属性中分析提取空间实体的拓扑一致性规则，包括空间隐喻规则与非空间隐喻规则。空间隐喻规则如几何对象间的方向约束、拓扑关系约束或者距离约束。如在LOD2、LOD3中语义为wall surface的空间对象相对应的几何表达为polygon，wall surface的语义属性会对该polygon几何元素有法向量方向约束，即可提取相应规则；非空间隐喻规则如门窗与墙体的is-part-of的关系。
(1) 语义属性对空间对象的方向约束需考虑CityGML多层次的特点，模型中不同的LOD层次所含语义对象信息不同，则有不同的方向约束规则。具体见表 1。
| CityGML语义对象 | 所属LOD | 语义规则下的几何特性 |
| ground surface | LOD2、LOD 3、LOD 4 | 法向向下 |
| wall surface | LOD2、LOD 3、LOD 4 | 法向水平 |
| roof surface | LOD2、LOD 3、LOD 4 | 方向不定 |
| window | LOD3、LOD 4 | 法向水平 |
| door | LOD3、LOD 4 | 法向水平 |
| ceiling surface | LOD4 | 法线向下 |
| interior wall surface | LOD4 | 法向水平 |
| floor surface | LOD4 | 法向向上 |
(2) 语义属性对空间对象间的拓扑关系约束。该约束考虑不同语义对象间的拓扑关系，根据CityGML的数据特点，将其支持的三维拓扑关系分为相等(equal)、相离(disjoint)、相交(intersects)、相切(touch)、包含(contains)5种，且根据分层的几何结构规则，仅考虑相同几何基元间的拓扑关系。在LOD1和LOD2中，建筑物的几何信息是单个几何体或者复合面，针对该语义对象考察点与点、线与线、面与面间的拓扑关系；在LOD3中，门窗与墙以面表示，需要考虑门、窗等语义要素在几何上的拓扑关系，wall surface上的点与window/door上的点应在同一平面上，且门窗与墙面内环相切，墙面几何的内环间不能相离、相交或重叠等；在LOD4中建筑物为复合体，包括墙、门窗、室内房间、楼梯等几何体，室内几何对象之间不能相交，且被包含于整个建筑物外围，整个建筑物外围包括roof surface、wall surface和ground surface。部分空间对象的拓扑关系约束形式化表达见表 2。
| 几何要素 | 空间谓词 | 规则描述 | 规则形式化表达 | 所属LOD层 |
| 复合面 | 相切 | door/window与wall间的拓扑关系，即门窗所在几何面应与墙体所在几何面的内环相切 | 
|
LOD3 |
| 相离 | door与window之间相离，不能相交 | Pdoor∩Pwindow=∅ Pwindow∩Pwindow=∅ Pdoor表示门所代表的几何面 Pwindow表示窗所代表的几何面 |
LOD3 | |
| 相等 | 所有的面不重复 | P1≠P2≠…≠Pn | LOD3 | |
| 复合体 | 包含 | 如所有CityGML室内几何对象被整个外部几何所包含 | S∩Sroom=H S表示建筑物不含室内房间的几何实体，Sroom表示建筑内某个房间的几何实体，H表示S中内shell的面集合 |
LOD4 |
| 相离 | 如室内各个room之间不相交 | S1∩S2=∅ S1、S2表示建筑物室内两个房间实体 |
LOD4 | |
| 相切 | 如同楼层中room与上下楼板应相切于面 | dim(Sroom∩Sfloor)=2 Sfloor表示一个floor实体，Sroom表示楼层内一个room |
LOD4 | |
| 相等 | 建筑物不重复 | S1≠S2≠…≠Sn | LOD4 |
(3) 语义属性对空间对象的距离约束。根据GB50352-2019《民用建筑设计统一标准》及GB 50096-2011《住宅设计规范》中对常见室内对象的设计规定，可分析提取针对空间对象层高、净高等距离关系的约束规则[25-26](表 3)。
| 几何一致性检测的语义对象 | 规则描述 | 所属LOD |
| floor surface/floor surfae | 普通住宅层高宜为2.70~2.90 m | LOD4 |
| room/living room | 室内净高不低于2.40 m；厨房、卫生间的室内净高不应低于2.20 m | LOD4 |
| stair way | 两股人流最小宽度不应小于1.10 m，梯段净高不宜小于2.20 m | LOD4 |
(4) 非空间隐喻对空间对象的约束。一些语义关系，如is-part-of等只能在语义层面提取[27]，如每个room必定与一个door关联并与其相切，3楼的stair way必须高于2楼的stair way等常识性规则，则需对stair way所在楼层、floor所在楼层信息，门窗与房间is-part-of的关联关系等非空间隐喻提取拓扑一致性检测规则，具体示例见表 4。
| 拓扑一致性检测的语义对象 | 规则描述 | 所属LOD |
| stair way | 楼层属性大的高程大于楼层属性小的高程 | LOD4 |
| floor | 楼层属性大的高程大于楼层属性小的高程 | LOD4 |
| floor/ground surface与stair way | 两个floor surface之间或者floor surface与ground surface之间必存在stair way | LOD4 |
| room与window/door | room必定有与某个door存在is-part-of的关联关系 | LOD4 |
2 自动检测方法 2.1 分层检测方法
本文以产生式表示法表示规则集，一个产生式对应一个规则，并采用点、线、面、体、多体等不同几何维度的分层检测算法进行拓扑一致性检测。高维度的拓扑检测建立在低维度的几何一致性基础上。
上述规则，最终可转换为诸如平面拟合、距离判断等约束条件，采用最小二乘法、三维布尔运算等算法，进行判断与求解。由于涉及的检测算法较多，考虑到平面检测在模型检测中比较典型，本文仅以此为代表进行介绍。
2.1.1 投影降维检测三维空间数据相较二维空间数据复杂，由直线的几何投影特性可知：三维平面向不垂直于该面的另一个二维面进行投影，投影面上各点的顺序关系不会发生改变，即投影点之间的拓扑关系与相对应空间点之间的拓扑关系一致。因此在拓扑一致性检测的过程中，三维面的自相交检测可通过降维投影到相应的二维平面，通过二维面的自相交检测判断空间面是否相交。但由于投影有积聚性和收缩性，若向垂直于该面的平面进行投影，则所有点会积聚为一条线，投影后各点之间的距离也会发生改变。为了避免投影的积聚性和减少投影收缩可能引起的判断误差，根据空间三维面的单位法向量判定投影方向，选择沿法向量的3个方向轴的最大矢量轴进行投影，若法向量为(0，0，1)，则沿Z轴向xy平面投影。因此点、线、面几何元素的自相交、面上内外环等拓扑一致性通过投影降维进行检测。
2.1.2 拓扑一致性检测算法针对三维平面的检测规则共有5条，面元素的几何一致性检测建立在该面所有环都具有拓扑一致性的基础上，并考虑几何结构和语义对象两种规则约束。针对规则(5)，通过同一面上所有空间点到拟合面的距离可检测构成该面的点是否在给定的精度阈值范围内共面，采用最小二乘法生成拟合面[28]；针对规则(6)、(7)、(8)，先将三维面投影降维到二维面，然后对投影面做环方向检测、包含检测以及内环相交检测；针对规则(9)，可首先将三维面进行三角剖分，通过判断剖分后三角形之间的连通关系，检测面内部的连通性。本文提出另外一种方法，由于面的内部连通性主要由内环与外环、内环之间的关系决定，当其关系满足以下的约束时，面的内部是连通的(各个条件为递进关系，即从上到下，满足条件即可终止判断)：
(1) 面有且只有一个外环，无内环。
(2) 面有一个内环，且内环与外环没有公共顶点。
(3) 面有一个内环，内环与外环之间只有一个公共顶点。
(4) 面有多于一个内环，每个内环都与外环没有公共顶点。
(5) 每个内环与外环之间最多只有一个公共顶点，且内环之间没有公共顶点。
(6) 内环之间最多只有一个公共顶点，且两个相邻内环之间只能通过该公共顶点进行连通。
当CityGML细节层次模型为LOD3及以上，需要检测solid中多个shells的拓扑关系或者多个solids间的拓扑关系，在CGAL库中以Nef-Polyhedral表示对象，对对象进行布尔运算，求其交集、并集、差集，以volume数判断shell之间的空间关系。
2.2 不同LOD中规则的适用性与依赖性不同LOD层级的建筑物适应不同的一致性规则。LOD1是以规则长方体表示建筑物的没有屋顶结构的块状模型，因而在LOD1的检测中，只考虑简单体的拓扑规则，不考虑几何面的语义属性、面上内外环的拓扑关系或者几何体中内外shell的拓扑关系等；LOD2中可以确切地表示建筑物屋顶的形状及外观，是一个带有屋顶结构的简化模型，所有几何面均拥有语义属性信息，因此在这个层级的几何模型中，在LOD1拓扑检测的基础上，增加语义属性对几何结构约束的规则检测；LOD3的建筑模型在LOD2的基础上，增加了墙面、门窗、屋顶形状等详细的几何外观信息以及对“开口”的描述，在LOD3建筑模型的拓扑一致性检测过程中，以LOD2的拓扑检测为基础，增加语义属性中复合面间拓扑关系的规则检测，以检测墙面与窗或者窗与窗之间的几何拓扑关系；LOD4在LOD3的基础上加入了三维空间对象的内部结构信息，如房间、楼梯、家具等，这里不考虑家具等室内设施，因而在LOD3拓扑检测的基础上，增加了内外shell及体与体之间的拓扑一致性检测。由此可见，多层次的语义规则更有助于CityGML模型不同细节层次几何模型的拓扑一致性检测需求。
在同一个CityGML模型中，前文中描述的规则之间具有依赖关系。按照几何要素进行划分。环要素的检测优先于其他所有规则进行检测，其后依次是面要素规则检测，外壳要素规则检测，体要素规则检测，多体要素规则检测；后一个环节的检测依赖于前一个环节规则检测结果，在前一个检测正确或者经过修复后正确之后，后面环节的检测方能开始。
此外，在各个几何要素层的检测过程中，各个子规则也需要遵循一定的依赖性。下面以面检测为例进行说明。
面检测子规则具体见表 5。在保证其前序环规则检测正确后，ID为203的多边形平面性检测是面要素几何一致性检测的前提，ID为201、202、205、206、207的规则之间相互独立，依赖于ID为203、208的规则。其他子规则之间的检测依赖关系见表 6。
| 检测规则ID | 检测规则函数名 | 规则层级 | 检测规则说明 |
| 201 | find(“Self-intersection”) | 面 | 面边界自相交检测 |
| 202 | find(“Duplicate Rings”) | 面 | 内外环重复性检测 |
| 203 | is_face_planar_distance2 plane() | 面 | 面上个点是否共面检测 |
| 204 | is_face_planar_normals() | 面 | 面内各个子面是否法向量方向一致 |
| 205 | find(“Interior is discon-nected”) | 面 | 面内区域是否连通 |
| 206 | find(“Hole lies outside shell”) | 面 | 内外环关系检测 |
| 207 | find(“Holes are nested”) | 面 | 内环嵌套检测 |
| 208 | orientation() | 面 | 外环与内环方向检测 |
| 209 | wallwithWorD() | 多个面 | 门与墙的关系检测 |
| 210 | doorwithWindow() | 多个面 | 门与窗户之间的关系检测 |
| 检测规则ID | 201 | 202 | 203 | 204 | 205 | 206 | 207 | 208 |
| 201 | √ | √ | ||||||
| 202 | √ | √ | ||||||
| 203 | ||||||||
| 204 | √ | √ | √ | √ | √ | √ | √ | |
| 205 | √ | √ | ||||||
| 206 | √ | √ | ||||||
| 207 | √ | √ | ||||||
| 208 | √ | |||||||
| 209 | √ | |||||||
| 210 | √ |
2.3 坐标容差及平面容差
CityGML的几何模型未显示存储solid中的拓扑关系，对solid的每个面采用单独存储的方法。因此在一个正方形中，同一个点将被存储3次，则可能会出现3个坐标表示同一个点。当点数值很大时，在计算机的浮点计算下会导致solid出现未闭合的不一致性，因此需考虑一定的舍入误差。ISO19107及GML中要求CityGML几何模型中表面必须为平面，OGC在OGC CityGML QIE中对几何模型平面性引入平面容差和角度容差[16]。
3 CityGML建筑物模型的拓扑一致性修复一条拓扑一致性规则对应一个类型的拓扑问题。当建筑物模型的拓扑不一致类型较多且数据缺失严重时，较难修复全部类型的错误。这里以常见的平面几何要素、方向一致性、闭合性、重复性等拓扑一致性问题的修复为例进行说明。
3.1 平面几何要素修复平面几何的特性与几何实体的语义属性互相联系，针对平面几何要素的修复，需以语义属性规则进行辅助。以下对常见的地面(ground surface)、墙面(wall surface)及屋顶(roof surface)分别进行说明。
地面的法向量垂直向下，计算该几何面所有点所在最小二乘拟合面，面上所有点到拟合面的对应投影点构成平面Z，设平面Pxy平行于xy水平面并经过平面中最低点，求得Z平面上所有点到该平面的投影点，则对应投影点坐标构成的面即是修复后的几何面，如图 4(a)所示；墙面应与地面共享一条边，以边与法向量构建平面Pyz，原几何面上所有点至该平面的投影点构成的面即为修复面，如图 4(b)所示；屋顶主要分为平屋顶与三角屋顶，平面屋顶修复与地面修复方法类似，如图 4(c)所示；三角屋顶需计算该屋顶所有点所在最小二乘拟合面，所有点沿Z轴至拟合面的投影点构成的面即为屋顶的修复面，如图 4(d)所示。
|
| 图 4 平面几何体的修复 Fig. 4 Fixing errors in plane geometry |
3.2 方向一致性修复
在各个等级规则中，都有对方向性的检测。以solid为例，其外shell法向应指向体外，而内shell方向应指向体的内部。该错误可通过遍历每个面，根据该面的正确朝向，对顺序不一致的边界点索引进行反转，从而修复方向不一致问题。在其他的检测等级中，也可依次处理方向不一致问题。
3.3 闭合性修复闭合性修复分为环与体的情况。环的首末点必须相同，以保证环的闭合性，可通过复制第一点坐标并加入到点序列末尾，对环的闭合性进行修复。
多面体以半边数据结构表示，即几何模型中每条边均存在仅被两个面所共享的特征，如图 5所示，有向边〈V5, V6〉存在于面D中，有向边〈V6, V5〉存在于面C中，则边V5V6仅被C、D两个面共享。针对数据结构表示的特性，可通过存储的半边信息以修复未闭合问题。如图 5所示，边V1V2仅存在于面A中，有向边〈V1, V2〉存在于面A中，则〈V2, V1〉是缺失面B的一边，同理可得有向边〈V1, V4〉、〈V4, V3〉、〈V3, V2〉，即缺失面B由边〈V1, V4〉、〈V4, V3〉、〈V3, V2〉、〈V2, V1〉组成。对于LOD2/LOD3的闭合修复，可通过面B的法向量方向赋予语义信息。若单位法向量向上为屋顶，向下为地面；若法向量为水平方向则需判断面B为墙面还是门窗。有向边〈V1, V2〉在面A中，若其为面A外环中的边，则边〈V2, V1〉为面B的外环，标识B为墙面；否则为门或者窗。对于LOD4的体闭合修复，则为该Solid语义属性，其中针对未闭合面修复的算法如下：
|
| 图 5 建筑物体未闭合 Fig. 5 Solid geometry is not closed |
begin
多面体_polyhedron;
if (_polyhedron->is_closed() == false)
while (_polyhedron->size_of_border_edges()>0)
Halfedge_handle he=++(_polyhedron->border_halfedges_begin());
he->vertex()->point().x()
he->vertex()->point().y()
he->vertex()->point().z()//get the point in the missing polygon
_polyhedron->fill_hole(he);
_polyhedron->normalize_border();
end if
end
3.4 重复性修复以环的重复性修复为例，在环一致性规则中，同一个环的点只允许首末点相同，其他点不允许重复。根据点重复类型不同，修复分为两种情况，一种情况如图 6(a)所示，R=(V0, V1, V2, V3, V3, V0)，其中点V3在点序列中出现连续重复；另一种情况如图 6(b)中的R=(V0, V1, V4, V3, V2, V4, V0)，点V4在序列中不是连续重复。针对(a)图中的重复性，可将重复的第2个点直接删除；针对(b)图中的不一致问题，可将该环分为R=(V0, V1, V4, V0)及R=(V4, V2, V3, V4)两个环，因此多边形被分割成多个多边形。针对环重复、面重复或者shell等不一致问题，则可直接进行删除。
|
| 图 6 环上点的两种重复性情况 Fig. 6 Two repetitive cases of points on the ring |
4 试验与分析
本文研究以CityGML建筑物模型的LOD2数据、LOD3数据及LOD4数据为例，进行试验验证分析。其中LOD2含32个建筑物，LOD3含400个建筑物，数据来源于OGC官网，LOD4建筑物含10个语义对象，由某栋IFC建筑物模型经过数据转换得到的CityGML LOD4数据。统一设置舍入误差为1 mm，平面容差分别为0.01 m，平面角度容差1°。检测系统界面如图 7所示。
|
| 图 7 几何检测试验系统界面 Fig. 7 System interface of geometric validation |
4.1 LOD2建筑物模型的拓扑检测
针对CityGML LOD2试验数据的几何一致性检测结果统计如图 8所示。该模型数据在几何一致性检测中存在环自相交、非平面、非流形边、方向不一致、数据不完整导致几何模型未闭合的不一致问题。如图 8(a)所示，该数据中不一致性问题占比达到一半以上；图 8(b)中给出了不一致问题的类型，其中非平面的几何错误(对应203号面要素规则)相比其他不一致性问题为多。
|
| 图 8 LOD2建筑物几何一致性检测试验结果 Fig. 8 Results of geometric consistency validation on LOD2 model of building |
图 9给出了某些建筑物中出现的非流形问题(图 6(a))及方向不一致问题(图 6(b))共存的情况。其中红色表示roof surface，灰色表示wall surface，蓝色表示ground surface。
|
| 图 9 LOD2数据中多种几何不一致共存的情况 Fig. 9 Multiple geometric inconsistency coexist in LOD2 models |
4.2 LOD3建筑物模型的拓扑错误检测
针对CityGML LOD3试验数据的几何一致性检测结果统计如图 10。该模型数据在几何一致性检测中存在环自相交、环与环相交、非平面、门窗与墙面拓扑关系错误，导致几何模型中内环在外环外部或者门与墙面相交的内外环拓扑关系错误、模型数据不完整导致几何模型未闭合及模型方向不一致等问题。如图 10(a)所示，该数据中不一致性问题占比约33%；图 10(b)中给出了不一致问题的类型，其中非平面的几何错误(对应203号面要素规则)同样最为突出。图 11给出了某些建筑物中因门窗与墙面相离或者相交从而导致数据存在环与环相交的几何不一致问题，如图 11(a)、(b)，其中红色表示roof surface，灰色表示wall surface，土黄色表示窗，棕色表示门。
|
| 图 10 LOD3建筑物几何一致性检测试验结果 Fig. 10 Results of geometric consistency validation on LOD3 models of building |
|
| 图 11 LOD3数据几何不一致可视化 Fig. 11 Visualizing results of inconsistency of LOD3 models |
4.3 建筑物LOD4模型检测试验
CityGML LOD4细节层次的建筑物数据中，语义对象以solid表示，有room、living room、hall way、window、door、floor surface、ground surface、wall surface、roof surface、stair way。该LOD4的建筑物数据存在非流形、方向不一致、room不存在与door的关联关系导致该room不与其他室内空间连通、几何体未闭合及语义对象间的拓扑关系错误等几何不一致等问题。其中低等级的几何错误较少，但方向不一致性错误较多。living room、roof surface、ground surface都出现了该错误，如room与door的is-part-of的错误(图 12(a))、楼板与楼梯的拓扑关系错误(图 12(b))。
|
| 图 12 LOD4几何不一致可视化 Fig. 12 Visualizing results of inconsistency of LOD4 models |
5 结语
本文以CityGML建筑物模型中常见的建筑物为研究对象，根据国际与国内的相关标准，从建筑物几何结构特征和语义属性信息中抽取对建筑物几何结构、拓扑关系、方向关系、距离关系具有约束影响及可转化为约束条件的非空间隐喻的语义信息了构建了建筑物的拓扑一致性规则集。考虑到不同LOD层次和不同语义属性的几何对象，结合拓扑一致性规则对单个语义对象的点、线、面、体等几何要素及多个语义对象间的空间关系，进行模型空间数据拓扑一致性的分层检测，提出了系统的不一致性检测方法，检测同一空间现象/空间实体在表达时存在的拓扑不一致问题，并进行了典型修复算法研究。试验结果表明所提出的规则和算法可检测出存在几何不一致的建筑物模型空间数据，可为后期修复提供有力支持，确保数据集的拓扑一致性，提高空间数据质量，提升三维数据在集成中的互操作性。
| [1] |
朱庆. 三维GIS及其在智慧城市中的应用[J]. 地球信息科学学报, 2014, 16(2): 151-157. ZHU Qing. Full three-dimensional GIS and its key roles in smart city[J]. Journal of Geo-Information Science, 2014, 16(2): 151-157. |
| [2] |
李德仁, 邵振峰, 杨小敏. 从数字城市到智慧城市的理论与实践[J]. 地理空间信息, 2011, 9(6): 1-5. LI Deren, SHAO Zhenfeng, YANG Xiaomin. Theory and practice from digital city to smart city[J]. Geospatial Information, 2011, 9(6): 1-5. |
| [3] |
单杰, 李志鑫, 张文元. 大规模三维城市建模进展[J]. 测绘学报, 2019, 48(12): 1523-1541. SHAN Jie, LI Zhixin, ZHANG Wenyuan. Recent progress in large-scale 3D city modeling[J]. Acta Geodaetica et Cartographica Sinica, 2019, 48(12): 1523-1541. DOI:10.11947/j.AGCS.2019.20190471 |
| [4] |
HUGHES T J R, COTTRELL J A, BAZILEVS Y. Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement[J]. Computer Methods in Applied Mechanics & Engineering, 2005, 194(39-41): 4135-4195. |
| [5] |
ISIKDAG U, UNDERWOOD J, AOUAD G. An investigation into the applicability of building information models in geospatial environment in support of site selection and fire response management processes[J]. Advanced Engineering Informatics, 2008, 22(4): 504-519. DOI:10.1016/j.aei.2008.06.001 |
| [6] |
HAALA N, FRITSCH D, PETER M, et al. Pedestrian navigation and modeling for indoor environments[C]//Proceeding of the 7th International Symposium on Mobile Mapping Technology. Cracow, Poland: [s.n.], 2011.
|
| [7] |
ZHU Jie, YOU Xiong, XIA Qing, et al. Object-oriented battlefield environment simulation process object model based on task-driven[J]. Journal of Geodesy and Geoinformation Science, 2019, 2(3): 31-43. |
| [8] |
BILJECKI F, LEDOUX H, DU Xin, et al. The most common geometric and semantic errors in CityGML datasets[C]//ISPRS Annals of the Photogrammetry, Remote Sensing and Spatial Information Sciences. Athens, Greece: [s.n.], 2016: 13-22.
|
| [9] |
XU D, VAN OOSTEROM P, ZLATANOVA S. A methodology for modelling of 3D spatial constraints[M]. Berne: Springer International Publishing, 2017: 95-117.
|
| [10] |
BOGDAHN J, COORS V. Towards an automated healing of 3D urban models[C]//Proceedings of the 5th International Conference on 3D Geoinformation. Berlin: ISPRS, 2010, 38: 13-17.
|
| [11] |
BILJECKI F, STOTER J, LEDOUX H, et al. Applications of 3D city models: state of the art review[J]. ISPRS International Journal of Geo-Information, 2015, 4(4): 2842-2889. DOI:10.3390/ijgi4042842 |
| [12] |
GRÖGER G, PLVMER L. How to achieve consistency for 3D city models[J]. Geoinformatica, 2011, 15(1): 137-165. DOI:10.1007/s10707-009-0091-6 |
| [13] |
MÄNTYLÄ M. Introduction to solid modeling: tutorial overheads[M]//Introduction to Solid Modeling: Tutorial Overheads. Helsinki: Helsinki University of Technology Press, 1988.
|
| [14] |
FORTUNE S. Polyhedral modelling with multiprecision integer arithmetic[J]. Computer-Aided Design, 1997, 29(2): 123-133. DOI:10.1016/S0010-4485(96)00041-3 |
| [15] |
LEDOUX H, VERBREE E, SI Huang. Geometric validation of GML solids with the constrained delaunay tetrahedralization[C]//Proceedings of the 4th International Workshop on 3D Geo-Information. Ghent, Belgium: [s.n.], 2009: 143-148.
|
| [16] |
LEDOUX H. On the validation of solids represented with the international standards for geographic information[J]. Computer-Aided Civil and Infrastructure Engineering, 2013, 28(9): 693-706. DOI:10.1111/mice.12043 |
| [17] |
GRÖGER G, KOLBE T, NAGEL C, HÄFELE K H. OGC city geography markup language (CityGML) encoding standard[S/OL]. (2012-04-04). http://mediatum.ub.tum.de/?id=1145731.
|
| [18] |
汤圣君, 朱庆, 赵君峤. BIM与GIS数据集成: IFC与CityGML建筑几何语义信息互操作技术[J]. 土木建筑工程信息技术, 2014, 6(4): 11-17. TANG Shengjun, ZHU Qing, ZHAO Junqiao. Towards interoperating of BIM and GIS model: geometric and semantic integration of CityGML and IFC building models[J]. Journal of Information Technology in Civil Engineering and Architecture, 2014, 6(4): 11-17. |
| [19] |
AMIREBRAHIMI S, RAJABIFARD A, SABRI S, et al. Spatial information in support of 3D flood damage assessment of buildings at micro level: a review[C]//Proceedings of ISPRS Annals of the Photogrammetry, Remote Sensing and Spatial Information Sciences. Athens, Greece: ISPRS, 2016: 73-81.
|
| [20] |
MÅNSSON U. BIM & GIS connectivity paves the way for really smart cities[J]. Perspektiv NR, 2015, 14(25): 19-24. |
| [21] |
ALBERTS B. Computer graphics: principles and practice[J]. 3rd ed. Boston, MA: Addison-Wesley Professional. 2013.
|
| [22] |
STADLER A, KOLBE T H. Spatio-semantic coherence in the integration of 3D city models[C]//Proceedings of the 5th International ISPRS Symposium on Spatial Data Quality ISSDQ 2007 in Enschede. Amsterdam: [s.n.], 2007.
|
| [23] |
OGC. OpenGIS geography markup language (GML3) implementation specification version 3.1.1[S/OL]. (2004-02-07). http://xml.coverpages.org/ni2003-02-06-c.html.
|
| [24] |
ISO 19107: 2019 Geographic information-Spatial schema[S]. ISO Technical Committee: ISO/TC 211 Geographic information/Geomatics. https://www.iso.org/standard/66175.html, 2019. 12.
|
| [25] |
中华人民共和国住房和城乡建设部. 民用建筑设计统一标准: GB 50352-2019[S]. 北京: 中国建筑工业出版社, 2019: 4. Ministry of Housing and Urban-Rural Development of the People's Republic of China. Uniform standard for design of civil buildings: GB 50352-2019[S]. Beijing: China Architecture &Building Press, 2019: 4. |
| [26] |
中华人民共和国住房和城乡建设部. 住宅设计规范: GB 50096-2011[S]. 北京: 中国建筑工业出版社, 2012: 12. Ministry of Housing and Urban-rural Development of the People's Republic of China. Design code for residential buildings: GB 50096-2011[S]. Beijing: China Architecture & Building Press, 2012: 12. |
| [27] |
GRÖGER G, KOLBE T H, CZERWINSKI A, et al. OpenGIS city geography markup language (CityGML) encoding standard, version 1.0.0[S/OL]. (2008-01-10). https://www.researchgate.net/publication/258255672_Open_GIS_city_geography_markup_language_City_GML_encoding_standard.
|
| [28] |
李伟, 孙元超, 李宗春, 等. 一种改进的最小二乘平面拟合算法[J]. 测绘科学, 2017, 42(1): 15-19, 100. LI Wei, SUN Yuanchao, LI Zongchun, et al. An improved least-square plane fitting algorithm[J]. Science of Surveying and Mapping, 2017, 42(1): 15-19, 100. |