[发明专利]卫星影像中建筑物水平矢量轮廓自动批量提取方法

说明书

技术领域

本发明涉及一种从单幅建筑群卫星影像中自动批量提取建筑物水平矢量轮廓的方法，特别是对俯视图为直线结构多边形的建筑物的水平矢量轮廓的自动批量提取。

背景技术

利用单幅卫星影像实现三维场景虚拟重建是一个很活跃的研究课题，它主要应用在城市建设规划、军事场景仿真、资源管理、抗震救灾模拟等方面。在真实场景的三维虚拟重建中，绝大部分为结构简单、俯视图为直线段结构多边形的普通建筑物，如何实现这类大量存在的普通建筑物的快速建模是高效重建三维建筑群的关键，而如何实现从建筑群影像中自动批量地提取建筑物的水平轮廓是高效重建三维建筑群的基础，并且决定着未来重建的虚拟场景中建筑群的分布格局与真实场景的符合程度。

当前的研究中，针对不同的应用目的，人们提出了各种各样的建筑物图像边缘提取算法，比较常见的有snake模型法、水平集曲线演化法、区域生长法等。以上方法的研究对象均为灰度图像，不适用于彩色图像，因此无法利用彩色图像中丰富的颜色信息。除此之外，snake模型法对初始位置敏感，需要依赖其他机制将初始轮廓放置在感兴趣的图像特征附近，否则轮廓提取会失败，目前大多采用手工选取的办法来设置初始边界，不仅非常繁琐，而且给轮廓线自动生成造成了困难。而较为常用的区域生长法存在着两点不足：第一，是初始种子点的选取问题。当前方法对种子点的选取大多数是手动选取，需要大量的人工干涉，耗时费力，效率很低。第二，是生长阈值的选取问题。生长阈值过大将会出现过分割，即生成的目标区域面积往往比实际面积大；而生长阈值过小又会导致分割不足，即目标区域生长不完整。因此不同建筑物需要选取不同的生长阈值，这一工作需要人工逐一根据经验完成。另外，snake模型法与区域生长法针对的是单个目标区域，也即每次只能生成一个目标区域，导致了方法的低效率。

以上方法提取的边缘轮廓，均由像素组成，这里称之为像素轮廓。本专利研究建筑物水平轮廓提取的目的，是为后续建筑群的三维虚拟重建提供建模数据，需要对提取出的建筑物像素轮廓进一步作直线矢量化处理，这里将直线矢量化处理后的轮廓称为矢量轮廓。目前的轮廓矢量化方法是将边缘像素提取出来并组成列表，然后对该列表进行直线段拟合，其中设置一个距离阈值，定义为偏离直线的最大距离，当超过该距离阈值后，直线段就会按比例一分为二。这种方法的结果受距离阈值选取的影响较大，为了得到较好的效果，影像中的每个建筑物都需要对应不同的距离阈值，以往设置距离阈值的方法往往是根据经验手工设置，操作效率低。

发明内容

为了克服现有技术的不足，本发明提供一种从单幅建筑群卫星影像中批量快速提取建筑物水平像素轮廓的自动化方法；进一步提出一种新的无需距离阈值及其他阈值选取的直线矢量化自动处理方法，对提取到的建筑物水平像素轮廓进行直线矢量化处理，得到满足后续建筑物三维虚拟重建所需的水平矢量轮廓。该方法特别适用于卫星影像中俯视图为直线段结构多边形的普通建筑物的水平矢量轮廓的批量快速提取，对航拍建筑群影像中的这类建筑群同样适用。

本发明解决其技术问题所采用的技术方案是：首先利用K-均值算法进行图像分类，得到建筑物的主干部分，取各建筑物区域质心解决初始种子点的选取问题。生长出所有种子点所在区域后，借助图像边缘信息将建筑物区域与周围区域分开，并根据矩形度、带状指数等特征参数去除掉非建筑物区域，从而实现建筑物水平像素轮廓的自动提取。然后利用Hough变换、块处理等技术对水平像素轮廓进行直线矢量化处理，最终批量得到所有建筑物的水平矢量轮廓。

本发明中，图像坐标系的原点位于图像左上角，x轴正方向竖直向下，y轴正方向水平向右。具体步骤如下：

步骤一、采用高斯平滑滤波器对建筑群卫星影像进行平滑处理，所述高斯平滑滤波器的方差σ＝0.8，窗口大小为7*7像素；然后在Hsv色彩空间内，保持每个像素的色调不变，按照步进值Δs调整饱和度，按照步进值Δv调整亮度，对卫星影像进行图像增强处理，得到源图像I_orig，其中Δs为0.06-0.08，Δv为0.05-0.07；

步骤二、轮廓自动提取，包括以下步骤：

1.不同建筑物种子点的自动生成，包括以下步骤：

1）在Lab色彩空间利用K-均值算法，在代表颜色信息的a、b色彩子空间上对源图像I_orig进行建筑物区域分割，K-均值算法执行中需要的分类类别数等于图像I_orig在a、b色彩子空间的二维直方图中的峰值个数，将分割出建筑物的图像存储为I_seg；