CN113076803B - 一种基于高分辨率遥感影像的建筑物矢量提取方法及系统 - Google Patents

一种基于高分辨率遥感影像的建筑物矢量提取方法及系统 Download PDF

Info

Publication number
CN113076803B
CN113076803B CN202110253141.4A CN202110253141A CN113076803B CN 113076803 B CN113076803 B CN 113076803B CN 202110253141 A CN202110253141 A CN 202110253141A CN 113076803 B CN113076803 B CN 113076803B
Authority
CN
China
Prior art keywords
building
vector
remote sensing
sensing image
vertex
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.)
Active
Application number
CN202110253141.4A
Other languages
English (en)
Other versions
CN113076803A (zh
Inventor
李子明
辛秦川
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sun Yat Sen University
Original Assignee
Sun Yat Sen University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Sun Yat Sen University filed Critical Sun Yat Sen University
Priority to CN202110253141.4A priority Critical patent/CN113076803B/zh
Publication of CN113076803A publication Critical patent/CN113076803A/zh
Application granted granted Critical
Publication of CN113076803B publication Critical patent/CN113076803B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00—Scenes; Scene-specific elements
    • G06V20/10—Terrestrial scenes
    • G06V20/176—Urban or other man-made structures
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/25—Fusion techniques
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/04—Architecture, e.g. interconnection topology
    • G06N3/045—Combinations of networks
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/02—Neural networks
    • G06N3/08—Learning methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/10—Segmentation; Edge detection
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00—Arrangements for image or video recognition or understanding
    • G06V10/40—Extraction of image or video features
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10032—Satellite or aerial image; Remote sensing

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Artificial Intelligence (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Computational Linguistics (AREA)
  • Multimedia (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Biomedical Technology (AREA)
  • Biophysics (AREA)
  • Mathematical Physics (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Image Analysis (AREA)

Abstract

本发明公开了一种基于高分辨率遥感影像的建筑物矢量提取方法及系统，其方法包括：获取目标区域的三波段遥感影像；将三波段遥感影像分别导入训练好的U‑Net网络模型和训练好的Cascade CNN网络模型，对应输出建筑物语义分割图和建筑物角点热力图；对建筑物语义分割图和建筑物角点热力图进行要素提取，构成建筑物矢量顶点集；将三波段遥感影像导入训练好的Cascade R‑CNN网络模型，并利用模型输出的建筑物预测边界框对建筑物矢量顶点集进行更新；基于更新后的建筑物矢量顶点集生成建筑物基本骨架，结合建筑物语义分割图从建筑物基本骨架中筛选出所有相关骨架数据；对所有相关骨架数据进行矢量融合，得到建筑物矢量结果。

Description

一种基于高分辨率遥感影像的建筑物矢量提取方法及系统
技术领域
本发明涉及遥感影像分析技术领域，尤其涉及一种基于高分辨率遥感影像的建筑物矢量提取方法及系统。
背景技术
随着城市化进程快速发展以及遥感技术不断进步，从高分辨率遥感影像中提取出建筑物信息对城市动态监测、规划管理等应用具有重要意义。面对传统人工标注方式所带来的效率问题和准确率问题，加上建筑物的实际大小形状、分布疏密情况等影响因素，使得实现从高分辨率遥感影像中进行自动化建筑物提取的功能成为一个极有挑战的任务。目前针对这一任务，相关技术人员已提出一种基于深度学习的语义分割方法，虽然可提升对建筑物的提取精度，但是最终得到的语义分割结果在执行矢量化处理后往往产生不规则多边形，其与真实地物斑块的边界存在很大的区别，很难直接应用于地理空间分析应用中，即采用此方法的实用效果不佳。
发明内容
本发明的目的在于克服现有技术的不足，本发明提供了一种基于高分辨率遥感影像的建筑物矢量提取方法及系统，可在保证大幅度提升建筑物提取精度的同时，克服现有技术所存在的建筑物边界提取误差。
为了解决上述问题，本发明提出了一种基于高分辨率遥感影像的建筑物矢量提取方法，所述方法包括：
获取目标区域的三波段遥感影像；
将所述三波段遥感影像分别导入训练好的U-Net网络模型和训练好的CascadeCNN网络模型，对应输出建筑物语义分割图和建筑物角点热力图；
对所述建筑物语义分割图和所述建筑物角点热力图进行要素提取，构成建筑物矢量顶点集；
将所述三波段遥感影像导入训练好的Cascade R-CNN网络模型，输出建筑物预测边界框，并利用所述建筑物预测边界框对所述建筑物矢量顶点集进行更新；
基于更新后的建筑物矢量顶点集生成建筑物基本骨架，结合所述建筑物语义分割图从所述建筑物基本骨架中筛选出所有相关骨架数据；
对所述所有相关骨架数据进行矢量融合，得到建筑物矢量结果。
可选的，所述对所述建筑物语义分割图和所述建筑物角点热力图进行要素提取，构成建筑物矢量顶点集包括：
从所述建筑物语义分割图中提取出建筑物斑块与所述三波段遥感影像边缘之间的所有相交点信息，构成虚拟顶点集；
利用非极大抑制算法从所述建筑物角点热力图中提取出建筑物的所有角点信息，构成真实顶点集；
对所述虚拟顶点集和所述真实顶点集进行数据合并，构成建筑物矢量顶点集。
可选的，所述利用所述建筑物预测边界框对所述建筑物矢量顶点集进行更新包括：
获取所述建筑物预测边界框上的最小边界点和最大边界点；
根据所述最小边界点和所述最大边界点建立阈值范围，从所述建筑物矢量顶点集中删除未落在所述阈值范围内的所有矢量顶点数据。
可选的，所述基于更新后的建筑物矢量顶点集生成建筑物基本骨架包括：
基于更新后的建筑物矢量顶点集包含有M个矢量顶点，构建允许涵盖所述M个矢量顶点的初始建筑物三角网；
将所述M个矢量顶点中的每一个矢量顶点插入所述初始建筑物三角网，基于每一个矢量顶点所在的三角形进行空外接圆检测，得到每一个矢量顶点所对应的单个三角网格；
基于所述M个矢量顶点所对应的M个三角网格已更新在所述初始建筑物三角网中，对构成当前初始建筑物三角网边框的所有关联三角形进行删除处理，得到建筑物基本骨架。
可选的，所述基于每一个矢量顶点所在的三角形进行空外接圆检测，得到每一个矢量顶点所对应的单个三角网格包括：
以第i(1≤i≤M)个矢量顶点所在的三角形为中心点，搜索与第i个矢量顶点所在的三角形相邻近的P个三角形，同时获取所述P个三角形中的每一个三角形所对应的一个外接圆；
从所述P个三角形所对应的P个外接圆中筛选出包含有第i个矢量顶点的R个外接圆，并对所述R个外接圆所对应的R个三角形进行删除处理；
将第i个矢量顶点分别与剩下的(P-R)个三角形所包含的全部顶点进行连接，得到第i个矢量顶点所对应的单个三角网格。
可选的，所述结合所述建筑物语义分割图从所述建筑物基本骨架中筛选出所有相关骨架数据包括：
获取所述建筑物基本骨架中所包含的K个骨架数据，并统计出所述K个骨架数据中的每一个骨架数据中所占的建筑物像素个数；
根据第j(1≤j≤K)个骨架数据所构成的三角边框位置，从所述建筑物语义分割图中提取出与第j个骨架数据所构成的三角边框位置相匹配的重合三角形，并统计出所述重合三角形中所占的建筑物像素个数；
基于第j个骨架数据中所占的建筑物像素个数以及与第j个骨架数据匹配的重合三角形中所占的建筑物像素个数，计算第j个骨架数据的关联比例；
通过K次循环操作得到所述K个骨架数据所对应的K个关联比例，从所述K个关联比例中筛选出大于预设值的T个关联比例，并从所述建筑物基本骨架中筛选出所述T个关联比例所对应的T个骨架数据。
可选的，所述对所述所有相关骨架数据进行矢量融合，得到建筑物矢量结果包括：
基于并集处理方式对所述T个骨架数据中的相邻每两个骨架数据进行合并，得到建筑物矢量结果。
另外，本发明实施例还提供了一种基于高分辨率遥感影像的建筑物矢量提取系统，所述系统包括：
图像获取模块，用于获取目标区域的三波段遥感影像；
图像解析模块，用于将所述三波段遥感影像分别导入训练好的U-Net网络模型和训练好的Cascade CNN网络模型，对应输出建筑物语义分割图和建筑物角点热力图；
数据提取模块，用于对所述建筑物语义分割图和所述建筑物角点热力图进行要素提取，构成建筑物矢量顶点集；
数据筛选模块，用于将所述三波段遥感影像导入训练好的Cascade R-CNN网络模型，输出建筑物预测边界框，并利用所述建筑物预测边界框对所述建筑物矢量顶点集进行更新；
骨架建立模块，用于基于更新后的建筑物矢量顶点集生成建筑物基本骨架，结合所述建筑物语义分割图从所述建筑物基本骨架中筛选出所有相关骨架数据；
矢量输出模块，用于对所述所有相关骨架数据进行矢量融合，得到建筑物矢量结果。
另外，本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述中任意一项所述的基于高分辨率遥感影像的建筑物矢量提取方法。
在本发明实施例中，通过结合三种具备不同功能的深度学习网络模型对三波段遥感影像中的相关形态控制点进行提取与筛选，可大幅度提升建筑物提取精度，减少冗余矢量标注；通过三角剖分框架对筛选后的相关形态控制点进行面状要素的构建，使得最终得到的建筑物斑块呈现出规则形状，以克服现有技术所存在的建筑物边界提取误差；本发明实施例可为城市动态监测、规划管理等应用实施提供便利条件。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见的，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其它的附图。
图1是本发明实施例中的基于高分辨率遥感影像的建筑物矢量提取方法的流程示意图；
图2是本发明实施例中的基于高分辨率遥感影像的建筑物矢量提取系统的组成示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
实施例
请参阅图1，图1示出了本发明实施例中的基于高分辨率遥感影像的建筑物矢量提取方法的流程示意图。
如图1所示，一种基于高分辨率遥感影像的建筑物矢量提取方法，所述方法包括：
S101、获取目标区域的三波段遥感影像；
在具体实施过程中，通过搭建高光谱传感器来采集到目标区域的三波段遥感影像，且在所述三波段遥感影像上建立一个平面坐标系，以满足后续对所述三波段遥感影像的分析应用中可统一延用该平面坐标系来表示出各个数据点信息为pi＝(xi，yi)。
S102、将所述三波段遥感影像分别导入训练好的U-Net网络模型和训练好的Cascade CNN网络模型，对应输出建筑物语义分割图和建筑物角点热力图；
其中，所述训练好的U-Net网络模型内部设置有运算精度最高的基于ResNet34(ResNet：Residual Network，残差网络)的U-Net语义分割网络，所述训练好的Cascade CNN(CNN：Convolutional Neural Network，卷积神经网络)网络模型内部设置有运算精度最高的基于VGG19(VGG：Visual Geometry Group，视觉几何组)的Cascade CNN角点提取网络。
S103、对所述建筑物语义分割图和所述建筑物角点热力图进行要素提取，构成建筑物矢量顶点集；
具体实施过程包括：首先从所述建筑物语义分割图中提取出建筑物斑块与所述三波段遥感影像边缘之间的所有相交点信息，构成虚拟顶点集为Pv＝{pv1，pv2，…，pvN}；其次利用非极大抑制算法从所述建筑物角点热力图中提取出建筑物的所有角点信息，构成真实顶点集为Pr＝{pr1，pr2，…，prN′}；最后对所述虚拟顶点集Pv和所述真实顶点集Pr进行数据合并，构成包含有(N+N′)个顶点的建筑物矢量顶点集为P＝Pv∪Pr。
S104、将所述三波段遥感影像导入训练好的Cascade R-CNN网络模型，输出建筑物预测边界框，并利用所述建筑物预测边界框对所述建筑物矢量顶点集进行更新；
具体实施过程包括：首先获取所述建筑物预测边界框上的最小边界点ptopleft＝(xmin，ymin)和最大边界点pbottomright＝(xmax，ymax)；其次根据所述最小边界点和所述最大边界点建立阈值范围为[ptopleft，pbottomright]，从所述建筑物矢量顶点集P中删除未落在所述阈值范围内的所有矢量顶点数据，此时更新后的建筑物矢量顶点集P′内的任意一个矢量顶点pa均需要满足：{pa＝(xa，ya)|xmin≤xa≤xmax且ymin≤ya≤ymax，pa∈P′}。其中，所述训练好的Cascade R-CNN(R-CNN：region-based Convolutional Neural Network，区域卷积神经网络)网络模型内部设置有运算精度最高的基于ResNet101的Cascade R-CNN目标检测网络。
S105、基于更新后的建筑物矢量顶点集生成建筑物基本骨架，结合所述建筑物语义分割图从所述建筑物基本骨架中筛选出所有相关骨架数据；
具体实施过程包括：
(1)基于更新后的建筑物矢量顶点集P′包含有M个矢量顶点，构建允许涵盖所述M个矢量顶点的初始建筑物三角网，进一步表现为：首先设定一个包含所述M个矢量顶点的初始矩形，再连接所述初始矩形的任意一条对角线，构成初始建筑物三角网；
(2)将所述M个矢量顶点中的每一个矢量顶点插入所述初始建筑物三角网，基于每一个矢量顶点所在的三角形进行空外接圆检测，得到每一个矢量顶点所对应的单个三角网格；
进一步的，首先以第i(1≤i≤M)个矢量顶点所在的三角形为中心点，搜索与第i个矢量顶点所在的三角形相邻近的P个三角形，同时获取所述P个三角形中的每一个三角形所对应的一个外接圆；其次，从所述P个三角形所对应的P个外接圆中筛选出包含有第i个矢量顶点的R个外接圆，并对所述R个外接圆所对应的R个三角形进行删除处理；最后，将第i个矢量顶点分别与剩下的(P-R)个三角形所包含的全部顶点进行连接，得到第i个矢量顶点所对应的单个三角网格。
(3)基于所述M个矢量顶点所对应的M个三角网格已更新在所述初始建筑物三角网中，对构成当前初始建筑物三角网边框的所有关联三角形进行删除处理，得到建筑物基本骨架；
其中，对构成当前初始建筑物三角网边框的所有关联三角形进行删除处理实际表现为：将与上述步骤(1)所生成的所述初始矩形的边框相关联的所有三角形进行删除。此外，本发明实施例主要通过依赖Python软件的SciPy库与Bowyer-Watson算法的相结合来完成建筑物基本骨架的构建。
(4)获取所述建筑物基本骨架中所包含的K个骨架数据，并统计出所述K个骨架数据中的每一个骨架数据中所占的建筑物像素个数为
Figure BDA0002960436360000071
其中1≤j≤K；
(5)根据第j个骨架数据所构成的三角边框位置，从所述建筑物语义分割图中提取出与第j个骨架数据所构成的三角边框位置相匹配的重合三角形，并统计出所述重合三角形中所占的建筑物像素个数为
Figure BDA0002960436360000072
(6)基于第j个骨架数据中所占的建筑物像素个数以及与第j个骨架数据匹配的重合三角形中所占的建筑物像素个数，计算第j个骨架数据的关联比例为
Figure BDA0002960436360000073
(7)通过K次循环操作得到所述K个骨架数据所对应的K个关联比例，从所述K个关联比例中筛选出大于预设值的T个关联比例，并从所述建筑物基本骨架中筛选出所述T个关联比例所对应的T个骨架数据。
其中，本发明实施例选定所述预设值为50％，此时可从所述K个关联比例中筛选出大于50％的T个关联比例，即所述T个关联比例所对应的T个骨架数据被判定属于建筑物三角形；相反的，剩下的(K-T)个小于等于50％的关联比例所对应的(K-T)个骨架数据即被判定属于非建筑物三角形。
S106、对所述所有相关骨架数据进行矢量融合，得到建筑物矢量结果。
具体的，基于并集处理方式对通过步骤S105所得到的所述T个骨架数据中的相邻每两个骨架数据进行合并，得到建筑物矢量结果，进一步表现为：基于Python软件的GeoPandas库，通过调用库中的Dissolve函数对所述T个骨架数据进行矢量融合处理，以形成单个矢量多边形输出。
在本发明实施例中，通过结合三种具备不同功能的深度学习网络模型对三波段遥感影像中的相关形态控制点进行提取与筛选，可大幅度提升建筑物提取精度，减少冗余矢量标注；通过三角剖分框架对筛选后的相关形态控制点进行面状要素的构建，使得最终得到的建筑物斑块呈现出规则形状，以克服现有技术所存在的建筑物边界提取误差；本发明实施例可为城市动态监测、规划管理等应用实施提供便利条件。
实施例
请参阅图2，图2示出了本发明实施例中的基于高分辨率遥感影像的建筑物矢量提取系统的组成示意图。
如图2所示，一种基于高分辨率遥感影像的建筑物矢量提取系统，所述系统包括：
图像获取模块201，用于获取目标区域的三波段遥感影像；
在具体实施过程中，通过搭建高光谱传感器来采集到目标区域的三波段遥感影像，且在所述三波段遥感影像上建立一个平面坐标系，以满足后续对所述三波段遥感影像的分析应用中可统一延用该平面坐标系来表示出各个数据点信息为pi＝(xi，yi)。
图像解析模块202，用于将所述三波段遥感影像分别导入训练好的U-Net网络模型和训练好的Cascade CNN网络模型，对应输出建筑物语义分割图和建筑物角点热力图；
其中，所述训练好的U-Net网络模型内部设置有运算精度最高的基于ResNet34(ResNet：Residual Network，残差网络)的U-Net语义分割网络，所述训练好的Cascade CNN(CNN：Convolutional Neural Network，卷积神经网络)网络模型内部设置有运算精度最高的基于VGG19(VGG：Visual Geometry Group，视觉几何组)的Cascade CNN角点提取网络。
数据提取模块203，用于对所述建筑物语义分割图和所述建筑物角点热力图进行要素提取，构成建筑物矢量顶点集；
具体实施过程包括：首先从所述建筑物语义分割图中提取出建筑物斑块与所述三波段遥感影像边缘之间的所有相交点信息，构成虚拟顶点集为Pv＝{pv1，pv2，…，pvN}；其次利用非极大抑制算法从所述建筑物角点热力图中提取出建筑物的所有角点信息，构成真实顶点集为Pr＝{pr1，pr2，…，prN′}；最后对所述虚拟顶点集Pv和所述真实顶点集Pr进行数据合并，构成包含有(N+N′)个顶点的建筑物矢量顶点集为P＝Pv∪Pr。
数据筛选模块204，用于将所述三波段遥感影像导入训练好的Cascade R-CNN网络模型，输出建筑物预测边界框，并利用所述建筑物预测边界框对所述建筑物矢量顶点集进行更新；
具体实施过程包括：首先获取所述建筑物预测边界框上的最小边界点ptopleft＝(xmin，ymin)和最大边界点pbottomright＝(xmax，ymax)；其次根据所述最小边界点和所述最大边界点建立阈值范围为[ptopleft，pbottomright]，从所述建筑物矢量顶点集P中删除未落在所述阈值范围内的所有矢量顶点数据，此时更新后的建筑物矢量顶点集P′内的任意一个矢量顶点pa均需要满足：{pa＝(xa，ya)|xmin≤xa≤xmax且ymin≤ya≤ymax，pa∈P′}。其中，所述训练好的Cascade R-CNN(R-CNN：region-based Convolutional Neural Network，区域卷积神经网络)网络模型内部设置有运算精度最高的基于ResNet101的Cascade R-CNN目标检测网络。
骨架建立模块205，用于基于更新后的建筑物矢量顶点集生成建筑物基本骨架，结合所述建筑物语义分割图从所述建筑物基本骨架中筛选出所有相关骨架数据；
具体实施过程包括：
(1)基于更新后的建筑物矢量顶点集P′包含有M个矢量顶点，构建允许涵盖所述M个矢量顶点的初始建筑物三角网，进一步表现为：首先设定一个包含所述M个矢量顶点的初始矩形，再连接所述初始矩形的任意一条对角线，构成初始建筑物三角网；
(2)将所述M个矢量顶点中的每一个矢量顶点插入所述初始建筑物三角网，基于每一个矢量顶点所在的三角形进行空外接圆检测，得到每一个矢量顶点所对应的单个三角网格；
进一步的，首先以第i(1≤i≤M)个矢量顶点所在的三角形为中心点，搜索与第i个矢量顶点所在的三角形相邻近的P个三角形，同时获取所述P个三角形中的每一个三角形所对应的一个外接圆；其次，从所述P个三角形所对应的P个外接圆中筛选出包含有第i个矢量顶点的R个外接圆，并对所述R个外接圆所对应的R个三角形进行删除处理；最后，将第i个矢量顶点分别与剩下的(P-R)个三角形所包含的全部顶点进行连接，得到第i个矢量顶点所对应的单个三角网格。
(3)基于所述M个矢量顶点所对应的M个三角网格已更新在所述初始建筑物三角网中，对构成当前初始建筑物三角网边框的所有关联三角形进行删除处理，得到建筑物基本骨架；
其中，对构成当前初始建筑物三角网边框的所有关联三角形进行删除处理实际表现为：将与上述步骤(1)所生成的所述初始矩形的边框相关联的所有三角形进行删除。此外，本发明实施例主要通过依赖Python软件的SciPy库与Bowyer-Watson算法的相结合来完成建筑物基本骨架的构建。
(4)获取所述建筑物基本骨架中所包含的K个骨架数据，并统计出所述K个骨架数据中的每一个骨架数据中所占的建筑物像素个数为
Figure BDA0002960436360000101
其中1≤j≤K；
(5)根据第j个骨架数据所构成的三角边框位置，从所述建筑物语义分割图中提取出与第j个骨架数据所构成的三角边框位置相匹配的重合三角形，并统计出所述重合三角形中所占的建筑物像素个数为
Figure BDA0002960436360000102
(6)基于第j个骨架数据中所占的建筑物像素个数以及与第j个骨架数据匹配的重合三角形中所占的建筑物像素个数，计算第j个骨架数据的关联比例为
Figure BDA0002960436360000103
(7)通过K次循环操作得到所述K个骨架数据所对应的K个关联比例，从所述K个关联比例中筛选出大于预设值的T个关联比例，并从所述建筑物基本骨架中筛选出所述T个关联比例所对应的T个骨架数据。
其中，本发明实施例选定所述预设值为50％，此时可从所述K个关联比例中筛选出大于50％的T个关联比例，即所述T个关联比例所对应的T个骨架数据被判定属于建筑物三角形；相反的，剩下的(K-T)个小于等于50％的关联比例所对应的(K-T)个骨架数据即被判定属于非建筑物三角形。
矢量输出模块206，用于对所述所有相关骨架数据进行矢量融合，得到建筑物矢量结果。
具体的，基于并集处理方式对通过所述骨架建立模块205所得到的所述T个骨架数据中的相邻每两个骨架数据进行合并，得到建筑物矢量结果，进一步表现为：基于Python软件的GeoPandas库，通过调用库中的Dissolve函数对所述T个骨架数据进行矢量融合处理，以形成单个矢量多边形输出。
在本发明实施例中，通过结合三种具备不同功能的深度学习网络模型对三波段遥感影像中的相关形态控制点进行提取与筛选，可大幅度提升建筑物提取精度，减少冗余矢量标注；通过三角剖分框架对筛选后的相关形态控制点进行面状要素的构建，使得最终得到的建筑物斑块呈现出规则形状，以克服现有技术所存在的建筑物边界提取误差；本发明实施例可为城市动态监测、规划管理等应用实施提供便利条件。
本发明实施例提供的一种计算机可读存储介质，所述计算机可读存储介质上存储有可执行的计算机程序，该计算机程序被处理器执行时实现上述实施例所提出的基于高分辨率遥感影像的建筑物矢量提取方法。其中，所述计算机可读存储介质包括但不限于任何类型的盘(包括软盘、硬盘、光盘、CD-ROM、和磁光盘)、ROM(Read-Only Memory，只读存储器)、RAM(Random AcceSS Memory，随即存储器)、EPROM(EraSable Programmable Read-Only Memory，可擦写可编程只读存储器)、EEPROM(ElectricallyEraSableProgrammableRead-Only Memory，电可擦可编程只读存储器)、闪存、磁性卡片或光线卡片。也就是说，存储设备包括由设备(如计算机、手机等)以可读的形式存储或传输信息的任何介质，可以是只读存储器、磁盘或光盘等。
以上对本发明实施例所提供的一种基于高分辨率遥感影像的建筑物矢量提取方法及系统进行了详细介绍，本文中采用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。

Claims (9)

1.一种基于高分辨率遥感影像的建筑物矢量提取方法，其特征在于，所述方法包括：
获取目标区域的三波段遥感影像；
将所述三波段遥感影像分别导入训练好的U-Net网络模型和训练好的Cascade CNN网络模型，对应输出建筑物语义分割图和建筑物角点热力图；
对所述建筑物语义分割图和所述建筑物角点热力图进行要素提取，构成建筑物矢量顶点集；
将所述三波段遥感影像导入训练好的Cascade R-CNN网络模型，输出建筑物预测边界框，并利用所述建筑物预测边界框对所述建筑物矢量顶点集进行更新；
基于更新后的建筑物矢量顶点集生成建筑物基本骨架，结合所述建筑物语义分割图从所述建筑物基本骨架中筛选出所有相关骨架数据；
对所述所有相关骨架数据进行矢量融合，得到建筑物矢量结果。
2.根据权利要求1所述的基于高分辨率遥感影像的建筑物矢量提取方法，其特征在于，所述对所述建筑物语义分割图和所述建筑物角点热力图进行要素提取，构成建筑物矢量顶点集包括：
从所述建筑物语义分割图中提取出建筑物斑块与所述三波段遥感影像边缘之间的所有相交点信息，构成虚拟顶点集；
利用非极大抑制算法从所述建筑物角点热力图中提取出建筑物的所有角点信息，构成真实顶点集；
对所述虚拟顶点集和所述真实顶点集进行数据合并，构成建筑物矢量顶点集。
3.根据权利要求1所述的基于高分辨率遥感影像的建筑物矢量提取方法，其特征在于，所述利用所述建筑物预测边界框对所述建筑物矢量顶点集进行更新包括：
获取所述建筑物预测边界框上的最小边界点和最大边界点；
根据所述最小边界点和所述最大边界点建立阈值范围，从所述建筑物矢量顶点集中删除未落在所述阈值范围内的所有矢量顶点数据。
4.根据权利要求1所述的基于高分辨率遥感影像的建筑物矢量提取方法，其特征在于，所述基于更新后的建筑物矢量顶点集生成建筑物基本骨架包括：
基于更新后的建筑物矢量顶点集包含有M个矢量顶点，构建允许涵盖所述M个矢量顶点的初始建筑物三角网；
将所述M个矢量顶点中的每一个矢量顶点插入所述初始建筑物三角网，基于每一个矢量顶点所在的三角形进行空外接圆检测，得到每一个矢量顶点所对应的单个三角网格；
基于所述M个矢量顶点所对应的M个三角网格已更新在所述初始建筑物三角网中，对构成当前初始建筑物三角网边框的所有关联三角形进行删除处理，得到建筑物基本骨架。
5.根据权利要求4所述的基于高分辨率遥感影像的建筑物矢量提取方法，其特征在于，所述基于每一个矢量顶点所在的三角形进行空外接圆检测，得到每一个矢量顶点所对应的单个三角网格包括：
以第i个矢量顶点所在的三角形为中心点，搜索与第i个矢量顶点所在的三角形相邻近的P个三角形，同时获取所述P个三角形中的每一个三角形所对应的一个外接圆，其中，1≤i≤M；
从所述P个三角形所对应的P个外接圆中筛选出包含有第i个矢量顶点的R个外接圆，并对所述R个外接圆所对应的R个三角形进行删除处理；
将第 i个矢量顶点分别与剩下的(P-R)个三角形所包含的全部顶点进行连接，得到第i个矢量顶点所对应的单个三角网格。
6.根据权利要求1所述的基于高分辨率遥感影像的建筑物矢量提取方法，其特征在于，所述结合所述建筑物语义分割图从所述建筑物基本骨架中筛选出所有相关骨架数据包括：
获取所述建筑物基本骨架中所包含的K个骨架数据，并统计出所述K个骨架数据中的每一个骨架数据中所占的建筑物像素个数；
根据第j个骨架数据所构成的三角边框位置，从所述建筑物语义分割图中提取出与第j个骨架数据所构成的三角边框位置相匹配的重合三角形，并统计出所述重合三角形中所占的建筑物像素个数，其中，1≤j≤K；
基于第j个骨架数据中所占的建筑物像素个数以及与第j个骨架数据匹配的重合三角形中所占的建筑物像素个数，计算第j个骨架数据的关联比例；
通过K次循环操作得到所述K个骨架数据所对应的K个关联比例，从所述K个关联比例中筛选出大于预设值的T个关联比例，并从所述建筑物基本骨架中筛选出所述T个关联比例所对应的T个骨架数据。
7.根据权利要求6所述的基于高分辨率遥感影像的建筑物矢量提取方法，其特征在于，所述对所述所有相关骨架数据进行矢量融合，得到建筑物矢量结果包括：
基于并集处理方式对所述T个骨架数据中的相邻每两个骨架数据进行合并，得到建筑物矢量结果。
8.一种基于高分辨率遥感影像的建筑物矢量提取系统，其特征在于，所述系统包括：
图像获取模块，用于获取目标区域的三波段遥感影像；
图像解析模块，用于将所述三波段遥感影像分别导入训练好的U-Net网络模型和训练好的Cascade CNN网络模型，对应输出建筑物语义分割图和建筑物角点热力图；
数据提取模块，用于对所述建筑物语义分割图和所述建筑物角点热力图进行要素提取，构成建筑物矢量顶点集；
数据筛选模块，用于将所述三波段遥感影像导入训练好的Cascade R-CNN网络模型，输出建筑物预测边界框，并利用所述建筑物预测边界框对所述建筑物矢量顶点集进行更新；
骨架建立模块，用于基于更新后的建筑物矢量顶点集生成建筑物基本骨架，结合所述建筑物语义分割图从所述建筑物基本骨架中筛选出所有相关骨架数据；
矢量输出模块，用于对所述所有相关骨架数据进行矢量融合，得到建筑物矢量结果。
9.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至7中任意一项所述的基于高分辨率遥感影像的建筑物矢量提取方法。
CN202110253141.4A 2021-03-03 2021-03-03 一种基于高分辨率遥感影像的建筑物矢量提取方法及系统 Active CN113076803B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110253141.4A CN113076803B (zh) 2021-03-03 2021-03-03 一种基于高分辨率遥感影像的建筑物矢量提取方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110253141.4A CN113076803B (zh) 2021-03-03 2021-03-03 一种基于高分辨率遥感影像的建筑物矢量提取方法及系统

Publications (2)

Publication Number Publication Date
CN113076803A CN113076803A (zh) 2021-07-06
CN113076803B true CN113076803B (zh) 2022-09-30

Family

ID=76612121

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110253141.4A Active CN113076803B (zh) 2021-03-03 2021-03-03 一种基于高分辨率遥感影像的建筑物矢量提取方法及系统

Country Status (1)

Country Link
CN (1) CN113076803B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115619907B (zh) * 2022-12-15 2023-06-09 航天宏图信息技术股份有限公司 一种自适应建筑物直角化方法和装置

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102004915A (zh) * 2010-10-19 2011-04-06 中国科学院遥感应用研究所 基于纹理特征空间的震损建筑遥感快速提取技术
CN102930540A (zh) * 2012-10-26 2013-02-13 中国地质大学（武汉） 城市建筑物轮廓检测的方法及系统
WO2014079477A1 (de) * 2012-11-20 2014-05-30 Siemens Aktiengesellschaft Verfahren zum automatischen erstellen von zwei- oder dreidimensionalen gebäudemodellen
WO2019136623A1 (en) * 2018-01-10 2019-07-18 Nokia Technologies Oy Apparatus and method for semantic segmentation with convolutional neural network
CN110232693A (zh) * 2019-06-12 2019-09-13 桂林电子科技大学 一种结合热力图通道及改进U-Net的图像分割方法
CN111460936A (zh) * 2020-03-18 2020-07-28 中国地质大学（武汉） 基于U-Net网络的遥感影像建筑物提取方法、系统、电子设备
CN111915608A (zh) * 2020-09-11 2020-11-10 北京百度网讯科技有限公司 建筑物提取方法、装置、设备和存储介质

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102004915A (zh) * 2010-10-19 2011-04-06 中国科学院遥感应用研究所 基于纹理特征空间的震损建筑遥感快速提取技术
CN102930540A (zh) * 2012-10-26 2013-02-13 中国地质大学（武汉） 城市建筑物轮廓检测的方法及系统
WO2014079477A1 (de) * 2012-11-20 2014-05-30 Siemens Aktiengesellschaft Verfahren zum automatischen erstellen von zwei- oder dreidimensionalen gebäudemodellen
WO2019136623A1 (en) * 2018-01-10 2019-07-18 Nokia Technologies Oy Apparatus and method for semantic segmentation with convolutional neural network
CN110232693A (zh) * 2019-06-12 2019-09-13 桂林电子科技大学 一种结合热力图通道及改进U-Net的图像分割方法
CN111460936A (zh) * 2020-03-18 2020-07-28 中国地质大学（武汉） 基于U-Net网络的遥感影像建筑物提取方法、系统、电子设备
CN111915608A (zh) * 2020-09-11 2020-11-10 北京百度网讯科技有限公司 建筑物提取方法、装置、设备和存储介质

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Evaluations and comparisons of rule-based and machine-learning-based methods to retrieve satellite-based vegetation phenology using MODIS and USA National Phenology Network data;Qinchuan Xin 等;《INTERNATIONAL JOURNAL OF APPLIED EARTH OBSERVATION AND GEOINFORMATION》;20201231;第93卷;第1-15页 *
Semantic Segmentation of Urban Buildings from VHR Remote Sensing Imagery Using a Deep Convolutional Neural Network;Yaning Yi 等;《remote sensing》;20190831;第11卷(第15期);第1-19页 *
基于全卷积神经网络的建筑物自动提取方法研究;郭亚;《中国优秀硕士学位论文全文数据库 (工程科技Ⅱ辑)》;20201115(第11期);第C028-16页 *
结合像元级和目标级的高分辨率遥感影像建筑物变化检测;张志强 等;《测绘学报》;20180131;第47卷(第1期);第102-112页 *

Also Published As

Publication number Publication date
CN113076803A (zh) 2021-07-06

Similar Documents

Publication Publication Date Title
CN110163198B (zh) 一种表格识别重建方法、装置和存储介质
Chen et al. A methodology for automated segmentation and reconstruction of urban 3-D buildings from ALS point clouds
Chen et al. Topologically aware building rooftop reconstruction from airborne laser scanning point clouds
US10115165B2 (en) Management of tax information based on topographical information
CN109948477B (zh) 一种提取图片中道路网络拓扑点的方法
CN112990086A (zh) 遥感影像建筑物检测方法、装置及计算机可读存储介质
CN114758337B (zh) 一种语义实例重建方法、装置、设备及介质
Khayyal et al. Creation and spatial analysis of 3D city modeling based on GIS data
Tarsha Kurdi et al. Automatic evaluation and improvement of roof segments for modelling missing details using Lidar data
CN103530901A (zh) 一种匹配真实三维环境的动态水域可视化方法及系统
CN113076803B (zh) 一种基于高分辨率遥感影像的建筑物矢量提取方法及系统
CN109657728B (zh) 样例生产方法及模型训练方法
Dragicevic et al. An application of fuzzy logic reasoning for GIS temporal modeling of dynamic processes
CN114708414A (zh) 一种工程三维地形数据切割融合方法与系统
CN112396701A (zh) 卫星图像的处理方法、装置、电子设备和计算机存储介质
CN113420109B (zh) 一种街道界面渗透率的测量方法、计算机及存储介质
CN112215864B (zh) 电子地图的轮廓处理方法、装置及电子设备
CN113870196A (zh) 一种基于锚点切图的图像处理方法、装置、设备和介质
CN117523403A (zh) 一种宅基地图斑变化检测方法、系统、设备及介质
CN115033972B (zh) 一种建筑主体结构批量单体化方法、系统及可读存储介质
CN114004552B (zh) 一种地下水质风险评估方法和系统
CN115861816A (zh) 一种立体低涡识别方法、装置、存储介质及终端
CN116071651A (zh) 一种均压场识别方法、装置、存储介质及终端
CN115984490A (zh) 风场特性的建模分析方法、系统、无人机设备及存储介质
CN115527028A (zh) 地图数据处理方法及装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant