CN101359371B - 基于数字高程模型的建筑物自动提取的方法 - Google Patents
基于数字高程模型的建筑物自动提取的方法 Download PDFInfo
- Publication number
- CN101359371B CN101359371B CN2008100411925A CN200810041192A CN101359371B CN 101359371 B CN101359371 B CN 101359371B CN 2008100411925 A CN2008100411925 A CN 2008100411925A CN 200810041192 A CN200810041192 A CN 200810041192A CN 101359371 B CN101359371 B CN 101359371B
- Authority
- CN
- China
- Prior art keywords
- buildings
- point
- partitioning portion
- minimum
- polygonal
- 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.)
- Expired - Fee Related
Links
Landscapes
- Image Processing (AREA)
- Image Analysis (AREA)
- Processing Or Creating Images (AREA)
Abstract
本发明涉及基于数字高程模型的建筑物自动提取的方法,该方法包括以下步骤:(1)预处理;(2)矢量化;(3)噪音去除;(4)边缘简化;(5)自动建筑物提取处理。与现有技术相比,本发明建立的完善的供电系统仿真体系,可以为以后的突发事件应急预案设计、物流系统仿真、自然灾难模拟与仿真等提供前提条件。
Description
技术领域
本发明涉及遥感影像中自动提取建筑物边缘并识别出建筑物的方法,尤其涉及基于数字高程模型(Digital Elevation Model,DEM)的建筑物自动提取的方法。
背景技术
有效地从城市航空遥感影像中提取建筑物信息对城市仿真和城市规划有着举足轻重的作用。然而由于建筑物的形状、大小和周边环境等的多样性,从遥感影像中自动提取建筑物边缘并识别出建筑物是十分困难。所以,找到一个切实有效的从遥感影像中自动提取建筑物的方法成为当前的研究热点。
不同的学者提出了许多不同的方法。Sohn等人提出的基于窗口傅立叶变换和二分树技术并结合建筑物单元形状的知识推理形成建筑物多变形。Stassopoulou等人通过使用一个基于Canny算子的多尺度区域分割与边缘分割相结合的方法对影像进行分割和区域特征提取,并结合其他成像条件在Bayesian网络支持下识别并提取房屋特征。Lin等人建立了一般光照模型下,房屋边缘线与阴影间的解析几何关系,可以用于影像目标地物相关关系的分析。Kirsch等人二维边缘检测算子进行推广也提出了一种三角形网格数据的边缘提取方法。
上述几种方法或是基于边缘提取,或是结合区域分割,其缺点是无法有效地将线、面信息及其上下文相关关系综合起来进行提取。这些没有被综合起来的信息使得最终提取的轮廓信息量过于冗余,为建筑物的提取后进行的重建工作带了非常大的工作量。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷,提供一种高效实用的基于数字高程模型的建筑物自动提取的方法。
本发明的目的可以通过以下技术方案来实现:基于数字高程模型的建筑物自动提取的方法,其特征在于,该方法包括以下步骤:
(1)预处理;
(2)矢量化;
(3)噪音去除;
(4)边缘简化;
(5)自动建筑物提取(Automatic Building Extraction,ABE)处理。
所述的步骤(1)中预处理包括:
a.进行最小值滤波;
b.使用最大值滤波;
c.通过设置高度值可以得到一个空缺部分;
d.结合面向对象分类技术所提供的对图像中分割区域的光谱、纹理和形态特征参数及上下文相关关系的丰富表达方式,得到比较可靠的分割图像。
所述的步骤(2)中矢量化为将经过预处理后得到的分割图像进行矢量化处理,包括以下步骤:
a.去除分割部分内部的像素点;
b.对每个分割部分中相关联的组成部分的多边形轮廓数量进行计算;
c.通过每个分割部分和多边形的轮廓上各点的x,y的坐标值来顺时针确定轮廓线上各点的顺序;
d.为了减少这个矢量图上点的个数,将相邻两点连成的直线上的点省略掉。
与现有技术相比,本发明建立的完善的供电系统仿真体系,可以为以后的突发事件应急预案设计、物流系统仿真、自然灾难模拟与仿真等提供前提条件。
具体实施方式
实施例1
本发明采用如下技术方案:
一、预处理
由于数学形态对图像信号的处理具有直观上的简单性和数学上的严谨性,所以在描述图像信号形态特征上具有独特的优势。而本文预处理中使用最小滤波对图像进行腐蚀可以消除物体边界点,将具有细小连通的两个建筑物分开。而在腐蚀之后进行的膨胀处理则是为了将与物体接触的所有背景点合并到该建筑物中,达到填补图像分割后建筑物中的空洞的最终效果。
因为一般建筑物都是高出它周围的地形平面的其他物体的,所以基于这个事实,我们可以通过处理后的DEM进行预处理。由于建筑物的信息大量的被包含在原始高程图与近似地形表面高程图的差值图像中,因此,我们可以利用灰度数学形态学[5]的方法来对原始高程图进行背景估计和消除。第一步是进行最小值滤波,
Z′=min{Z(x,y)|x,y∈W}
W是结构元素,在这种情况下取的是正方形窗口。对于数学形态学来说,最小值滤波是一种非常特别的腐蚀方法,因为在w(x,y)这个结构元素中z的值始终不变。
接下来使用最大值滤波,
Z″=max{Z′(x,y)|x,y∈W}
最大值滤波是一种非常特别的膨胀方法,
而经过这两步处理步骤以后,通过设置高度值可以得到一个空缺部分,
Z″=ZOW
即为近似的地形表面。
为了能够去除空缺部分中的所有建筑物信息,窗口大小可以根据前面建筑物的最大期望尺寸进行调整,使得结构元素w不完全包含在建筑物的轮廓里。如果近似地形表面正好汇集在建筑物群中的无房区域,我们设置数据,
d=Z-ZOW
即原始数字高程模型的高度值减去近似地形表面。
通过上述的预处理,再结合面向对象分类技术所提供的对图像中分割区域的光谱、纹理和形态特征参数及上下文相关关系的丰富表达方式,就能够得到比较可靠的分割图像。
二、矢量化
这一步要将经过预处理后得到的分割图像进行矢量化处理。由于分割部分内部的像素点会在网格中呈现轮廓,所以首先要去除。接着,我们要对每个分割部分中相关联的组成部分的多边形轮廓数量进行计算。然后通过每个分割部分和多边形的轮廓上各点的x,y的坐标值来顺时针确定轮廓线上各点的顺序。为了减少这个矢量图上点的个数,我们将相邻两点连成的直线上的点省略掉。经过上述追踪法,大体能把建筑物的外轮廓勾勒出来。但是不得不指出的是:这样生成的图形虽然基本恢复了建筑物的形状信息,但是由于受到噪声影响,建筑物的角点等位置会出现变形。基于这个事实,我们在矢量化处理之后将对这些噪音进行去除。
三、噪音去除
矢量化的最后一步是去除直线上各点,下面我们采用一种合并算法来消除离散化噪音。当在多边形迭代过程中三角形的最小高度值比预先设定的域值大,或者当已经达到了预先设定的点数最小值时,这个由多边形中i-1,i,i+1三点组成的相应三角形的高度值di即为多边形的最小高度值。利用这种合并算法可以成功去除多余点。在这种合并算法中,三角形的最小高度值将被看作为域值T。而这个最小高度与输入的数据的Δx,Δy密切相关。如果最小高度值 就将这条近似直线上的点全部去除。根据经验数据,这里设定K=2/3。
四、边缘简化
以上这些步骤提取的建筑物边缘一般比较复杂,不仅包含大量的无用信息,同时与真实边缘也有一定差距。为此,本文提出了将多边形的相邻边缘看作是直角用来进一步实现多边形的重建的边缘简化算法。这样的优化处理可以大大降低最后提取的数据量,并且减少了后续重建过程的工作量。
对于第二阶段的合并过程,我们一般采取将矩形覆盖一个角度的原则。当出现相邻两个拐角点时,我们也首选将矩形直接覆盖在相邻点上。描述长度的取值取决于数据和模型的相互适宜度以及模型的复杂度。但是当出现描述长度还能再缩短的情况时,我们会直接用一个点来代替这两个相邻点。
经过上面两个公式对描述长度DL以及残差平方和Ω进行计算,结果可以非常明显的说明这种边缘简化的优势。
五、ABE处理
综合以上各节,我们这里给出建筑物边缘自动提取的算法ABE。通过对数字高程模型DEM进行预处理,然后对所得图像进行矢量化和噪音去除两步骤的操作,最后对提取到的建筑物轮廓进行边缘简化。这一系列过程处理组成了本文提出的ABE(Automatic Building Extraction)算法,如下所示:
算法:ABE(Image)
对Image进行最小值滤波;
对处理后的图象进行最大值滤波;
newImage←矢量化(newImage);
while(i<n)//噪音去除{
di←由点i-1,i,i+1组成的三角形中的高度值;
if(dmin<2/3*sqrt(Δx2+Δy2))
去除这条近似直线上的点;
}
for(所有点in newImage)//边缘简化{
//p1和p2是相邻两个拐角点
将矩形直接覆盖在相邻点上;
if(描述长度DL还能再缩短)
用一个点来代替这两个相邻点;
}
Claims (3)
1.基于数字高程模型的建筑物自动提取的方法,其特征在于,该方法包括以下步骤:
(1)预处理:由于建筑物的信息大量的被包含在原始高程图与近似地形表面高程图的差值图像中,利用灰度数学形态学的方法来对原始高程图进行背景估计和消除,第一步是进行最小值滤波,
Z′=min{Z(x,y)|x,y∈W}
W是结构元素,在这种情况下取的是正方形窗口,对于数学形态学来说,最小值滤波是一种非常特别的腐蚀方法,因为在w(x,y)这个结构元素中z的值始终不变,
Z′=ZθW
接下来使用最大值滤波,
Z″=max{Z′(x,y)|x,y∈W}
最大值滤波是一种非常特别的膨胀方法,
而经过这两步处理步骤以后,通过设置高度值可以得到一个空缺部分,
Z″=ZOW
即为近似的地形表面,
为了能够去除空缺部分中的所有建筑物信息,窗口大小根据前面建筑物的最大期望尺寸进行调整,使得结构元素w不完全包含在建筑物的轮廓里,如果近似地形表面正好汇集在建筑物群中的无房区域,设置数据,
d=Z-ZOW
即原始数字高程模型的高度值减去近似地形表面,
通过上述的预处理,再结合面向对象分类技术所提供的对图像中分割区域的光谱、纹理和形态特征参数及上下文相关关系的丰富表达方式,就能够得到的分割图像;
(2)矢量化:由于分割部分内部的像素点会在网格中呈现轮廓,所以首先要去除,接着,要对每个分割部分中相关联的组成部分的多边形轮廓数量进行计算,然后通过每个分割部分和多边形的轮廓上各点的x,y的坐标值来顺时针确定轮廓线上各点的顺序,为了减少这个矢量图上点的个数,将相邻两点连成的直线上的点省略掉,经过上述追踪法,把建筑物的外轮廓勾勒出来;
(3)噪音去除:当在多边形迭代过程中三角形的最小高度值比预先设定的域值大,或者当已经达到了预先设定的点数最小值时,这个由多边形中i-1,i,i+1三点组成的相应三角形的高度值di即为多边形的最小高度值,利用这种合并算法可以成功去除多余点,这种合并算法中,三角形的最小高度值将被看作为域值T,而这个最小高度与输入的数据的Δx,Δy密切相关,如果最小高度值dmin<就将这条近似直线上的点全部去除,设定K=2/3;
(4)边缘简化:将多边形的相邻边缘看作是直角用来进一步实现多边形的重建的边缘简化算法,对于第二阶段的合并过程,采取将矩形覆盖一个角度的原则,当出现相邻两个拐角点时,将矩形直接覆盖在相邻点上,描述长度的取值取决于数据和模型的相互适宜度以及模型的复杂度,但是当出现描述长度还能再缩短的情况时,直接用一个点来代替这两个相邻点,
经过上面两个公式对描述长度DL以及残差平方和Ω进行计算;
(5)自动建筑物提取处理。
2.根据权利要求1所述的基于数字高程模型的建筑物自动提取的方法,其特征在于,所述的步骤(1)中预处理包括:
a.进行最小值滤波;
b.使用最大值滤波;
c.通过设置高度值可以得到一个空缺部分;
d.结合面向对象分类技术所提供的对图像中分割区域的光谱、纹理和形态特征参数及上下文相关关系,得到分割图像。
3.根据权利要求1所述的基于数字高程模型的建筑物自动提取的方法,其特征在于,所述的步骤(2)中矢量化为将经过预处理后得到的分割图像进行矢量化处理,包括以下步骤:
a.去除分割部分内部的像素点;
b.对每个分割部分中相关联的组成部分的多边形轮廓数量进行计算;
c.通过每个分割部分和多边形的轮廓上各点的x,y的坐标值来顺时针确定轮廓线上各点的顺序;
d.为了减少这个矢量图上点的个数,将相邻两点连成的直线上的点省略掉。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100411925A CN101359371B (zh) | 2008-07-30 | 2008-07-30 | 基于数字高程模型的建筑物自动提取的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100411925A CN101359371B (zh) | 2008-07-30 | 2008-07-30 | 基于数字高程模型的建筑物自动提取的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101359371A CN101359371A (zh) | 2009-02-04 |
CN101359371B true CN101359371B (zh) | 2011-05-18 |
Family
ID=40331817
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008100411925A Expired - Fee Related CN101359371B (zh) | 2008-07-30 | 2008-07-30 | 基于数字高程模型的建筑物自动提取的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101359371B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101980200B (zh) * | 2010-11-03 | 2013-09-04 | 东莞市高鑫机电科技服务有限公司 | 中国元素工程数据库构建方法、系统及其在工业设计领域的应用 |
CN102938066B (zh) * | 2012-12-07 | 2015-10-28 | 南京大学 | 一种基于多元数据重建建筑物外轮廓多边形的方法 |
CN105550691B (zh) * | 2015-12-29 | 2019-03-19 | 武汉大学 | 基于尺度空间的自适应山谷山脊线提取方法及系统 |
CN106447510A (zh) * | 2016-08-30 | 2017-02-22 | 水利部水土保持监测中心 | 一种梯田信息提取方法和系统 |
CN108154462B (zh) * | 2017-12-11 | 2021-07-13 | 北京控制工程研究所 | 一种白昼环境下星图背景的实时滤波方法 |
CN108280440A (zh) * | 2018-02-09 | 2018-07-13 | 三亚中科遥感研究所 | 一种果林识别方法和系统 |
CN112164142B (zh) * | 2020-10-21 | 2024-07-26 | 江苏科技大学 | 一种基于智能手机的楼盘采光模拟方法 |
JP7466429B2 (ja) * | 2020-10-22 | 2024-04-12 | 株式会社日立製作所 | 計算機システム及び計画評価方法 |
-
2008
- 2008-07-30 CN CN2008100411925A patent/CN101359371B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN101359371A (zh) | 2009-02-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101359371B (zh) | 基于数字高程模型的建筑物自动提取的方法 | |
EP3695384B1 (en) | Point cloud meshing method, apparatus, device and computer storage media | |
CN102004922B (zh) | 基于骨架特征的高分辨率遥感影像飞机提取方法 | |
JP2019514126A5 (zh) | ||
CN107481274B (zh) | 一种三维作物点云的鲁棒性重构方法 | |
CN107784661A (zh) | 基于区域生长法的变电站设备红外图像分类识别方法 | |
Sohn et al. | An implicit regularization for 3D building rooftop modeling using airborne lidar data | |
CA2776238A1 (fr) | Procede et dispositif pour obtenir une image d'un document defroisse a partir d'une image de ce document lorsqu'il est froisse | |
Cao et al. | Sharp feature extraction in point clouds | |
Zakšek et al. | An effective DEM generalization with basic GIS operations | |
Sahu et al. | Vector map generation from aerial imagery using deep learning | |
CN107909625B (zh) | 一种基于等高线的山顶点提取方法 | |
Matkan et al. | Spatial analysis for outlier removal from LiDAR data | |
CN112734675A (zh) | 一种基于金字塔模型和非局部增强密集块的图像去雨方法 | |
CN102663815B (zh) | 一种基于水平集的lod2建筑物模型构建方法 | |
KR101890536B1 (ko) | 건설 시공 단계에서의 3차원 설계 모델 데이터와 레이저 스캔 데이터 간 비교 분석을 위한 차별적 데이터 감소 방법 및 시스템 | |
CN109977965B (zh) | 一种在遥感机场图像中确定检测目标的方法及装置 | |
WO2008021949A2 (en) | Geospatial modeling system for separating foliage data from building data based upon loose and strict tolerance noise filtering operations and related methods | |
CN112561797B (zh) | 花卉浮雕模型构建方法及基于线条图的花卉浮雕重建方法 | |
Neelima et al. | Rank based approach for extracting unit pixel width skeleton | |
Gold et al. | The crust and skeleton–applications in GIS | |
Pang et al. | An effective quad-dominant meshing method for unorganized point clouds | |
Zhang et al. | A Markov Random Field Model for individual tree detection from airborne laser scanning data | |
CN112116561A (zh) | 基于图像处理融合网络权值的电网传输线检测方法及装置 | |
CN118212368B (zh) | 一种基于不完整轮廓点云数据重构三维图像的方法及系统 |
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: 20110518 Termination date: 20170730 |
|
CF01 | Termination of patent right due to non-payment of annual fee |