CN102663815A - 一种基于水平集的lod2建筑物模型构建方法 - Google Patents
一种基于水平集的lod2建筑物模型构建方法 Download PDFInfo
- Publication number
- CN102663815A CN102663815A CN2012100906286A CN201210090628A CN102663815A CN 102663815 A CN102663815 A CN 102663815A CN 2012100906286 A CN2012100906286 A CN 2012100906286A CN 201210090628 A CN201210090628 A CN 201210090628A CN 102663815 A CN102663815 A CN 102663815A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- phi
- building
- msup
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000010276 construction Methods 0.000 claims abstract description 18
- 230000011218 segmentation Effects 0.000 claims abstract description 12
- 239000012634 fragment Substances 0.000 claims abstract description 7
- 238000012545 processing Methods 0.000 claims abstract description 6
- 239000013598 vector Substances 0.000 claims description 10
- 230000000873 masking effect Effects 0.000 claims description 5
- 238000010606 normalization Methods 0.000 claims description 3
- 238000011478 gradient descent method Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Landscapes
- Image Processing (AREA)
- Processing Or Creating Images (AREA)
Abstract
一种基于水平集的LOD2建筑物模型构建方法,属于应用水平集算法对DSM数据的分割处理领域。为了解决现有的基于二维图像的建筑物模型构建方法中存在的顶面结构不够细致而导致构建精度不高的问题。具体过程为:提取建筑物轮廓掩膜Ωm,选取DSM数据,配入到统一坐标系下;获取建筑物顶面数据T;获取建筑物顶面数据T的特征空间;进行多相水平集分割,获取子区域;提取每个子区域的点集,检测每个片元的边界点,获取建筑物每个基元的角点的图像坐标;建立建筑物顶面数据T的拓扑结构;根据航拍可见光图像提取建筑物表面的纹理数据,与建筑物不同基元相对应,完成LOD2建筑物模型的构建。用于LOD2级别的大规模建筑物的三维构建任务。
Description
技术领域
本发明属于应用水平集算法对DSM数据的分割处理领域,具体涉及一种基于水平集的LOD2建筑物模型构建方法。
背景技术
随着对数字城市需求的日益提高,城市和郊区人造地物的高精度模型化能够使人们参与到城市的规划设计与建设,对经济的发展和避难措施的策划提出参考性的建议。游戏娱乐、汽车导航、电子商务等领域的飞速发展促进了人们对数字城市的迫切需求。三维建筑物模型通常采用多尺度的表达方式,以满足不同的应用需求。根据CityGML的定义,LOD1(Level Of Detail 1)是指仅有平面屋顶的建筑物模型,LOD1是第一级尺度,LOD2(LevelOf Detail 2)建筑物模型具有不同类型的建筑物屋顶和真实的表面纹理。目前,许多研究集中在LOD1建筑物模型的构建方法上,主要研究方向是对具有精细几何结构以及表面纹理的建筑物模型的构建方法。
在过去的几十年里,三维建筑物模型的构建技术可以使用多种数据源,尤其是DSM(Digital Surface Model,数字表面模型)的应用愈加广泛,通过多角度航拍图像或LiDAR(Light Detection And Ranging)数据处理所得的高精度DSM数据,其分辨率高达0.3m,可提供丰富的三维信息,从而构建出精细的屋顶几何结构。但是现有技术还没有从DSM数据中提取复杂建筑物模型的顶面基元以及拓扑结构、实现较高精度的LOD2建筑物模型构建的方法。
发明内容
本发明是为了解决现有的基于二维图像的建筑物模型构建方法中存在的顶面结构不够细致而导致构建精度不高的问题,提出的一种基于水平集的LOD2建筑物模型构建方法。
本发明的一种基于水平集的LOD2建筑物模型构建方法,其具体过程为:
步骤1、提取建筑物轮廓作为建筑物轮廓掩膜Ωm,同时根据相对应的建筑物区域的地理坐标选取DSM数据,通过自动配准或人工配准的方式,将建筑物轮廓掩膜Ωm和DSM数据配入到统一的坐标系下;
步骤2、根据步骤1选取的建筑物轮廓掩膜Ωm和DSM数据获取建筑物顶面数据T;
步骤3、根据步骤2获得的建筑物顶面数据T获取建筑物顶面数据T的特征空间;
步骤4、根据步骤3获取的特征空间将建筑物区域进行多相水平集分割,获取子区域;
步骤5、根据多相水平集的定义分别提取步骤4获得的每个子区域的点集,然后检测建筑物每个片元的边界点,获取建筑物每个基元的角点的图像坐标;
步骤6、根据步骤5获得的的建筑物角点的图像坐标,与步骤1中选取的DSM数据中的地理坐标相对应,建立建筑物顶面数据T的拓扑结构;
步骤7、根据步骤6获取的建筑物顶面数据T的拓扑结构,同时根据相应的建筑物区域的航拍可见光图像提取建筑物表面的纹理数据,与建筑物的不同基元相对应,完成LOD2建筑物模型的构建。
本发明基于多相多通道水平集框架,将局部区域信息引入该框架,以DSM数据点的法向量作为特征向量进行水平集分割,分离出建筑物顶面的不同基元,从而实现顶面几何结构的构建。
本发明的优势在于能够从DSM数据中提取复杂建筑物模型的顶面基元以及拓扑结构,实现较高精度的LOD2建筑物模型构建。
本发明可用于LOD2级别的大规模建筑物的三维构建任务。
附图说明
图1是本发明的建筑物模型构建方法流程示意图;图2是多想水平集分割过程的流程示意图;图3是本发明所采用的多相水平集方法的示意图;图4和图5是本发明在水平集框架中引入局部化区域信息方法的示意图。
具体实施方式
具体实施方式一、结合图1说明本实施方式,一种基于水平集的LOD2建筑物模型构建方法,其具体过程为:
步骤1、提取建筑物轮廓作为建筑物轮廓掩膜Ωm,同时根据相对应的建筑物区域的地理坐标选取DSM数据,通过自动配准或人工配准的方式,将建筑物轮廓掩膜Ωm和DSM数据配入到统一的坐标系下;
步骤2、根据步骤1选取的建筑物轮廓掩膜Ωm和DSM数据获取建筑物顶面数据T;
步骤3、根据步骤2获得的建筑物顶面数据T获取建筑物顶面数据T的特征空间;
步骤4、根据步骤3获取的特征空间将建筑物区域进行多相水平集分割,获取子区域;
步骤5、根据多相水平集的定义分别提取步骤4获得的每个子区域的点集,然后检测建筑物每个片元的边界点,获取建筑物每个基元的角点的图像坐标;
步骤6、根据步骤5获得的的建筑物角点的图像坐标,与步骤1中选取的DSM数据中的地理坐标相对应,建立建筑物顶面数据T的拓扑结构;
步骤7、根据步骤6获取的建筑物顶面数据T的拓扑结构,同时根据相应的建筑物区域的航拍可见光图像提取建筑物表面的纹理数据,与建筑物的不同基元相对应,完成LOD2建筑物模型的构建。
具体实施方式二、本实施方式与具体实施方式一的区别在于,本实施方式是对具体实施方式一中步骤2的进一步限定,所述步骤2中获取建筑物顶面数据T的方法为:用建筑物轮廓掩膜Ωm与配准后的DSM数据进行对应点相乘,然后用获得的结果去除建筑物轮廓外点集,获得建筑物的顶面数据T。
具体实施方式三、本实施方式与具体实施方式一的区别在于,本实施方式是对具体实施方式一中步骤3的进一步限定,所述步骤3中获取建筑物顶面数据T的特征空间的具体过程为:对建筑物顶面数据T应用Delaunay三角形剖分算法计算每一个三角形的法向量,然后对法向量进行归一化处理,将建筑物顶面数据T中的每一个点pi的法向量Ni用共用点pi的相邻三角形的法向量的平均值来表示,结果作为建筑物顶面数据T的特征空间。
具体实施方式四、本实施方式与具体实施方式一的区别在于,本实施方式是对具体实施方式一中步骤4的进一步限定,所述步骤4中将建筑物区域进行多相水平集分割的具体过程为:
步骤4-1、在待分割数据I的定义域Ω内,以两条等间距、等半径的圆圈为初始曲线C1和C2,把定义域Ω划分为四个区域,分别将曲线C1和C2转化成水平集方程φ1和φ2的形式,令X表示任意点(x,y),用符号距离函数表示一个点X的水平集方程φ1(X)和φ2(X)为:
其中,D(X,C1)表示点X与曲线C1之间的Euclidean距离,
其中,D(X,C2)表示点X与曲线C2之间的Euclidean距离;
步骤4-2、分别计算两个水平集方程φ1(X)和φ2(X)曲线的长度分量和面积分量,长度分量分别为:
length(C1)=∫Ω|▽H(φ1(X))|dX=∫Ωδ(φ1(X))|▽φ1(X)|dX
length(C2)=∫Ω|▽H(φ2(X))|dX=∫Ωδ(φ2(X))|▽φ2(X)|dX
面积分量分别为:
area(inside(C1))=∫ΩH(φ1(X))dX
area(inside(C2))=∫ΩH(φ2(X))dX
其中,
φ为φ1或φ2;
步骤4-3、在多相水平集框架中引入局部化区域函数B(X,Y),设置局部化区域函数B(X,Y)的半径参数r,以局部化区域参数B(X,Y)为局部区域掩膜,计算以X为中心、r为半径的圆圈内的特征空间的平均值:
I(Y)表示待分割数据;
步骤4-4、设置权重因子μ1,μ2,ν1,ν2,λ1,λ2,λ3,λ4,定义多相水平集能量泛函:
Etotal=μ1·E1+μ2·E2+ν1·E3+ν2·E4+λ1·E5+λ2·E6+λ3·E7+λ4·E8;
步骤4-5、以局部化区域函数B(X,Y)为局部区域掩膜,计算以X为中心、r为半径的圆圈内的能量分量,
水平集方程φ1的能量泛函表示为:
水平集方程φ2的能量泛函表示为:
步骤4-6、根据梯度下降法,分别计算步骤4-5获得的能量泛函Etotal(φ1)和Etotal(φ2)的演化方程:
步骤4-7、判断能量泛函是否达到了最小值,如果是则完成将待分割区域进行多相水平集分割的过程,将待分割区域分割成四个子区域,如果否,则返回步骤4-2,迭代更新水平集方程φ1和φ2。
本实施方式中,步骤4-4中设置权重因子的经验值为μ1=μ2=0.2,ν1=ν2=λ1=λ2=λ3=λ4=1。步骤4-5的作用是减少冗余信息的影响,并加速计算的速度。设置局部化区域函数B(X,Y)的半径参数r的经验值是r=5。
具体实施方式五、本实施方式与具体实施方式四的区别在于,本实施方式是对具体实施方式四中步骤4-7的进一步限定,所述步骤4-7中判断能量泛函是否达到了最小值的依据为:通过设定阈值或设定迭代次数。
具体实施方式六、本实施方式与具体实施方式一的区别在于,本实施方式是对具体实施方式一中步骤5的进一步限定,所述步骤5中获取建筑物每个基元的角点的图像坐标具体方法为:
步骤5-1、令水平集方程φ1=0,φ2=0,获得建筑物顶面几何结构边缘的二值图;
Claims (6)
1.一种基于水平集的LOD2建筑物模型构建方法,其特征在于:该方法的具体过程为:
步骤1、提取待构建模型的建筑物轮廓作为该建筑物轮廓掩膜Ωm,同时根据相对应的建筑物区域的地理坐标选取DSM数据,通过自动配准或人工配准的方式,将建筑物轮廓掩膜Ωm和DSM数据配入到统一的坐标系下;
步骤2、根据步骤1选取的建筑物轮廓掩膜Ωm和DSM数据获取建筑物顶面数据T;
步骤3、根据步骤2获得的建筑物顶面数据T获取建筑物顶面数据T的特征空间;步骤4、根据步骤3获取的特征空间将建筑物区域进行多相水平集分割,获取子区域;
步骤5、根据多相水平集的定义分别提取步骤4获得的每个子区域的点集,然后检测建筑物每个片元的边界点,获取建筑物每个基元的角点的图像坐标;
步骤6、根据步骤5获得的建筑物角点的图像坐标,与步骤1中选取的DSM数据中的地理坐标相对应,建立建筑物顶面数据T的拓扑结构;
步骤7、根据步骤6获取的建筑物顶面数据T的拓扑结构,同时根据相应的建筑物区域的航拍可见光图像提取建筑物表面的纹理数据,与建筑物的不同基元相对应,完成LOD2建筑物模型的构建。
2.根据权利要求1所述的一种基于水平集的LOD2建筑物模型构建方法,其特征在于:所述步骤2中获取建筑物顶面数据T的方法为:用建筑物轮廓掩膜Ωm与配准后的DSM数据进行对应点相乘,然后用获得的结果去除建筑物轮廓外点集,获得建筑物的顶面数据T。
3.根据权利要求1所述的一种基于水平集的LOD2建筑物模型构建方法,其特征在于:所述步骤3中获取建筑物顶面数据T的特征空间的具体过程为:对建筑物顶面数据T应用Delaunay三角形剖分算法计算每一个三角形的法向量,然后对所有法向量进行归一化处理,将建筑物顶面数据T中的每一个点pi的法向量Ni用共用点pi的相邻三角形的法向量的平均值来表示,结果作为建筑物顶面数据T的特征空间。
4.根据权利要求1所述的一种基于水平集的LOD2建筑物模型构建方法,其特征在于:所述步骤4中将建筑物区域进行多相水平集分割的具体过程为:
步骤4-1、在待分割数据I的定义域Ω内,以两条等间距、等半径的圆圈为初始曲线C1和C2,把定义域Ω划分为四个区域,分别将曲线C1和C2转化成水平集方程φ1和φ2的形式,令X表示任意点(x,y),用符号距离函数表示一个点X的水平集方程φ1(X)和φ2(X)为:
其中,D(X,C1)表示点X与曲线C1之间的Euclidean距离,
其中,D(X,C2)表示点X与曲线C2之间的Euclidean距离;
步骤4-2、分别计算两个水平集方程φ1(X)和φ2(X)曲线的长度分量和面积分量,长度分量分别为:
length(C1)=∫Ω|▽H(φ1(X))|dX=∫Ωδ(φ1(X))|▽φ1(X)|dX
length(C2)=∫Ω|▽H(φ2(X))|dX=∫Ωδ(φ2(X))|▽φ2(X)|dX
面积分量分别为:
area(inside(C1))=∫ΩH(φ1(X))dX
area(inside(C2))=∫ΩH(φ2(X))dX
其中,
φ为φ1或φ2;
步骤4-3、在多相水平集框架中引入局部化区域函数B(X,Y),设置局部化区域函数B(X,Y)的半径参数r,以局部化区域参数B(X,Y)为局部区域掩膜,计算以X为中心、r为半径的圆圈内的特征空间的平均值:
I(Y)表示待分割数据;
步骤4-4、设置权重因子μ1,μ2,ν1,ν2,λ1,λ2,λ3,λ4,定义多相水平集能量泛函:
Etotal=μ1·E1+μ2·E2+ν1·E3+ν2·E4+λ1·E5+λ2·E6+λ3·E7+λ4·E8;
步骤4-5、以局部化区域函数B(X,Y)为局部区域掩膜,计算以X为中心、r为半径的圆圈内的能量分量,
水平集方程φ1的能量泛函表示为:
水平集方程φ2的能量泛函表示为:
步骤4-6、根据梯度下降法,分别计算步骤4-5获得的能量泛函Etotal(φ1)和Etotal(φ2)的演化方程:
步骤4-7、判断能量泛函是否达到了最小值,如果是则完成将待分割区域进行多相水平集分割的过程,将待分割区域分割成四个子区域,如果否,则返回步骤4-2,迭代更新水平集方程φ1和φ2。
5.根据权利要求4所述的一种基于水平集的LOD2建筑物模型构建方法,其特征在于:所述步骤4-7中判断能量泛函是否达到了最小值的依据为:通过设定阈值或设定迭代次数。
6.根据权利要求1所述的一种基于水平集的LOD2建筑物模型构建方法,其特征在于:所述步骤5中获取建筑物每个基元的角点的图像坐标具体方法为:
步骤5-1、令水平集方程φ1=0,φ2=0,获得建筑物顶面几何结构边缘的二值图;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210090628.6A CN102663815B (zh) | 2012-03-30 | 2012-03-30 | 一种基于水平集的lod2建筑物模型构建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210090628.6A CN102663815B (zh) | 2012-03-30 | 2012-03-30 | 一种基于水平集的lod2建筑物模型构建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102663815A true CN102663815A (zh) | 2012-09-12 |
CN102663815B CN102663815B (zh) | 2014-11-26 |
Family
ID=46773292
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210090628.6A Expired - Fee Related CN102663815B (zh) | 2012-03-30 | 2012-03-30 | 一种基于水平集的lod2建筑物模型构建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102663815B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105931294A (zh) * | 2016-04-19 | 2016-09-07 | 西南交通大学 | 一种将bim实体模型转换为多细节层次gis标准化模型的方法 |
CN106600696A (zh) * | 2016-12-01 | 2017-04-26 | 华中科技大学 | 一种结构及其多类型边界的水平集隐式表达方法 |
CN111161411A (zh) * | 2020-04-08 | 2020-05-15 | 速度时空信息科技股份有限公司 | 一种基于八叉树的三维建筑物模型lod方法 |
WO2023104598A1 (fr) | 2021-12-10 | 2023-06-15 | Akidaia | Procédé de contrôle d'accès à une zone à sécuriser et procédé d'initialisation associé |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090110267A1 (en) * | 2007-09-21 | 2009-04-30 | The Regents Of The University Of California | Automated texture mapping system for 3D models |
CN101655993A (zh) * | 2008-08-20 | 2010-02-24 | 武汉大学 | 复杂建筑物三维模型多分辨率建模方法 |
-
2012
- 2012-03-30 CN CN201210090628.6A patent/CN102663815B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090110267A1 (en) * | 2007-09-21 | 2009-04-30 | The Regents Of The University Of California | Automated texture mapping system for 3D models |
CN101655993A (zh) * | 2008-08-20 | 2010-02-24 | 武汉大学 | 复杂建筑物三维模型多分辨率建模方法 |
Non-Patent Citations (2)
Title |
---|
BING JIA ET AL.: "A NOVEL LEVEL SET FRAMEWORK FOR LOD2 BUILDING MODELING", 《2012 19TH IEEE INTERNATIONAL CONFERENCE ON IMAGE PROCESSING(ICIP)》, 3 October 2012 (2012-10-03), pages 1781 - 1784 * |
贾冰: "立体目标的模型构建与变化检测技术", 《万方学位论文》, 25 December 2012 (2012-12-25), pages 30 - 39 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105931294A (zh) * | 2016-04-19 | 2016-09-07 | 西南交通大学 | 一种将bim实体模型转换为多细节层次gis标准化模型的方法 |
CN105931294B (zh) * | 2016-04-19 | 2019-02-12 | 西南交通大学 | 一种将bim实体模型转换为多细节层次gis标准化模型的方法 |
CN106600696A (zh) * | 2016-12-01 | 2017-04-26 | 华中科技大学 | 一种结构及其多类型边界的水平集隐式表达方法 |
CN111161411A (zh) * | 2020-04-08 | 2020-05-15 | 速度时空信息科技股份有限公司 | 一种基于八叉树的三维建筑物模型lod方法 |
WO2023104598A1 (fr) | 2021-12-10 | 2023-06-15 | Akidaia | Procédé de contrôle d'accès à une zone à sécuriser et procédé d'initialisation associé |
FR3130481A1 (fr) | 2021-12-10 | 2023-06-16 | Akidaia | Procédé de contrôle d’accès à une zone à sécuriser et procédé d’initialisation associé |
Also Published As
Publication number | Publication date |
---|---|
CN102663815B (zh) | 2014-11-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chen et al. | Automatic building information model reconstruction in high-density urban areas: Augmenting multi-source data with architectural knowledge | |
CN112927370B (zh) | 三维建筑物模型构建方法、装置、电子设备及存储介质 | |
CN107341453B (zh) | 一种车道线提取方法及装置 | |
Jung et al. | Automated 3D wireframe modeling of indoor structures from point clouds using constrained least-squares adjustment for as-built BIM | |
Matei et al. | Building segmentation for densely built urban regions using aerial lidar data | |
CN103366631B (zh) | 一种制作室内地图的方法和制作室内地图的装置 | |
Ebrahimpour et al. | Vanishing point detection in corridors: using Hough transform and K-means clustering | |
CN103727930B (zh) | 一种基于边缘匹配的激光测距仪与相机相对位姿标定方法 | |
CN106600622A (zh) | 一种基于超体素的点云数据分割方法 | |
CN104484668A (zh) | 一种无人机多重叠遥感影像的建筑物轮廓线提取方法 | |
CN103761739B (zh) | 一种基于半全局能量优化的影像配准方法 | |
CN103714541A (zh) | 一种利用山体轮廓区域约束识别定位建筑物的方法 | |
CN104376595A (zh) | 一种基于机载LiDAR和GIS协同的三维道路生成方法 | |
CN103514604B (zh) | 电子散斑干涉条纹图像骨架线提取方法 | |
CN102663815B (zh) | 一种基于水平集的lod2建筑物模型构建方法 | |
US20210192840A1 (en) | Maps comprising covariances in multi-resolution voxels | |
WO2022104260A1 (en) | Data normalization of aerial images | |
CN115063465B (zh) | 一种基于激光雷达的无人车行驶路况建模方法 | |
CN108875127B (zh) | 一种计算机气象软件中基于风场数据的槽线修正方法 | |
Li et al. | A deep learning-based indoor acceptance system for assessment on flatness and verticality quality of concrete surfaces | |
Wang et al. | A method for detecting windows from mobile LiDAR data | |
Paiva et al. | Historical building point cloud segmentation combining hierarchical watershed transform and curvature analysis | |
CN116310355A (zh) | 一种复杂结构件的激光点云去噪和缺陷检测方法 | |
Huang et al. | Automatic CAD model reconstruction from multiple point clouds for reverse engineering | |
WO2023223284A1 (en) | System and method for triggering data transfer using progress tracking |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20141126 Termination date: 20150330 |
|
EXPY | Termination of patent right or utility model |