一种遥感图像建筑物提取方法与流程

本发明涉及一种遥感图像处理领域，具体说是一种遥感图像建筑物提取方法。

背景技术：

建筑物是城市主要地理要素之一，是各种城市专题图的重要内容，研究建筑物的提取对综合考察城市地理信息环境具有重要意义。随着高分辨率遥感影像获取技术的快速发展，遥感影像的处理、分析和应用有了更好的数据源，其数字产品则有了更广泛、更深入的应用。计算机图像处理技术、模式识别、人工智能等方面的都取得不同程度的进展，为高效地提取海量影像中的有效信息提供了可能。但建筑物信息的提取要比其他信息如道路、水体的获取难得多，主要原因如下：

(1)数据源主要是二维的遥感影像，大多数情况下缺少直接的三维数据；

(2)不同的遥感影像常因为光谱范围、分辨率、传感器的几何图像以及成像条件等因素的不同而有较大的差异；

(3)不同种类的建筑物其所表现出来的外观和纹理细节等千变万化，表现在遥感图像上差异很大，统一的建筑物模型库难以建立，这使得信息的自动提取变得相当困难；

(4)建筑物所处场景的复杂性，如对比度较低时、房屋相互遮挡、建筑物自身的阴影以及处在其它地物的阴影等，所以想自动地从背景中提取出边界清晰的建筑物较为困难。

技术实现要素：

本发明提供了一种遥感图像建筑物提取方法，可克服目前遥感影像中建筑物提取困难的问题，可以检测遥感影像中具有光谱均一性的建筑物目标，无需人工干预，自动化程度高。

为实现本发明的目标所采用的技术方案方法包括以下步骤：

步骤1：输入遥感图像i；

步骤2：对遥感图像i进行增强，得到增强后的图像ie；

步骤3：通过下式计算图像ie的内部灰度方差特征空间mgv：

式中，(x,y)为坐标点位置，m(x,y)为图像ie中坐标点(x,y)的平均灰度值，l为滑动窗口的边长；

步骤4：将兴趣区域r初始化为图像ie，将聚类种子点集合seedset初始化为空集，提取图像ie的种子点，包括以下步骤：

步骤4.1：计算兴趣区域r的最大灰度值vmax和最小灰度值vmin；

步骤4.2：计算兴趣区域r的灰度直方图h，并用h(i)表示灰度值为i的像素点的数量，其中，i∈[0,k-1]，k为兴趣区域r所包含的灰度等级；

步骤4.3：将兴趣区域r中灰度值等于max(h(i))的像素点添加到候选种子点集合seedset中，max为取最大值函数；

步骤4.4：遍历候选种子点集合seedset中的所有聚类种子点s，用下式计算兴趣区域r中的灰度变分值sd：

式中，vs为候选种子点s的灰度值，gi为兴趣区域r中第i个灰度等级，i∈[0,k-1]；

步骤4.5：判断兴趣区域r的同质特性，当sd大于th时，表示兴趣区域r是同质区域，结束步骤4，进入步骤5，否则进入步骤4.6，其中，th为同质因子；

步骤4.6：分裂兴趣区，当vs-f×sd≤gi≤vs+f×sd时，删除灰度值属于区间[vs-f×sd,vs+f×sd]的像素点，并将剩余的像素点分裂为两个兴趣区：即[vmin,vs-f×sd]和[vs+f×sd,vmax]，其中，f为灰度变分值调整系数；

步骤4.7：将兴趣区域r分别更新为[vmin,vs-f×sd]和[vs+f×sd,vmax]，迭代运行步骤4.1～步骤4.6；

步骤5：将步骤4.3中的候选种子点集合seedset的所有种子点映射到内部灰度方差特征空间mgv中，得到内部灰度方差特征空间中的聚类种子点集合mgvseedset；

步骤6：用聚类种子点集合mgvseedset中的种子点进行聚类操作，得到聚类对象集合co；

步骤7：对聚类对象集合co进行二值化运算，得到二值聚类对象集合bco；

步骤8：对二值聚类对象集合bco进行形态学滤波，得到候选建筑物集合cbs；

步骤9：分别计算候选建筑物集合cbs中每个对象的最小外接矩形，得到矩形建筑物集合rbs；

步骤10：提取建筑物的轮廓，即采用活动轮廓模型分别对矩形建筑物集合rbs中的每个对象所包含的图像ie的区域进行分割，得到建筑物集合bout；

步骤11：输出建筑物集合bout。

步骤2所述的增强方法采用直方图均衡化方法。

步骤6所述的聚类操作采用基于灰度值的最近邻聚类方法。

步骤8所述的形态学滤波包括用于去除小面积干扰区域的形态学开运算及细化操作。

步骤10所述的活动轮廓模型可以采用snake模型，并将其初始轮廓曲线设定为候选建筑物集合cbs中每个对象的最小外接矩形。

步骤11所述的输出包括显示建筑物集合bout中每个建筑物的轮廓信息。

本发明的有益效果是：解决了遥感影像中建筑物提取准确性不高的问题，达到完全自动化的效果。可以用于遥感影像制图、地理信息系统的数据获取和自动更新。

附图说明

图1是本发明的总体处理流程图；

图2是本发明的提取候选种子点的处理方法流程图。

具体实施方式

下面结合附图详细描述本发明的具体实施方式。

图1是本发明的总体处理流程图。如图1所示，101是输入遥感图像步骤，102是图像增强步骤，103是计算内部灰度方差特征空间步骤，104是提取候选种子点步骤，105是映射聚类种子点集合步骤，106是聚类步骤，107是二值化步骤，108是形态学滤波步骤，109是求最小外接矩形步骤，110是分割步骤，111是输出建筑物集合步骤。

步骤101：输入遥感图像i；

步骤102：利用直方图均衡化对遥感图像i进行增强，得到增强后的图像ie；

步骤103：通过下式计算图像ie的内部灰度方差特征空间mgv：

式中，(x,y)为坐标点位置，m(x,y)为图像ie中坐标点(x,y)的平均灰度值，l为滑动窗口的边长，将l设置为5；

步骤104：将兴趣区域r初始化为图像ie，将聚类种子点集合seedset初始化为空集，提取图像ie的种子点，并添加到候选种子点集合seedset中；

步骤105：将候选种子点集合seedset的所有种子点映射到内部灰度方差特征空间mgv中，得到内部灰度方差特征空间中的聚类种子点集合mgvseedset；

步骤106：用聚类种子点集合mgvseedset中的种子点进行聚类操作，其中聚类方法采用基于灰度值的最近邻聚类方法，得到聚类对象集合co；

步骤107：对聚类对象集合co进行二值化运算，得到二值聚类对象集合bco；

步骤108：对二值聚类对象集合bco进行形态学滤波，包括用于去除小面积干扰区域的形态学开运算及细化操作，得到候选建筑物集合cbs；

步骤109：分别计算候选建筑物集合cbs中每个对象的最小外接矩形，得到矩形建筑物集合rbs；

步骤110：提取建筑物的轮廓，即采用活动轮廓模型分别对矩形建筑物集合rbs中的每个对象所包含的图像ie的区域进行分割，活动轮廓模型采用snake模型，并将其初始轮廓曲线设定为候选建筑物集合cbs中每个对象的最小外接矩形，得到建筑物集合bout；

步骤111：输出建筑物集合bout，并显示建筑物集合bout中每个建筑物的轮廓信息。

图2是本发明的提取候选种子点的处理方法流程图。如图2所示，in4是输入兴趣区步骤，401是计算最大灰度值和最小灰度值步骤，402是计算灰度直方图步骤，403是提取候选种子点步骤，404是计算灰度变分值步骤，405是判断兴趣区是否为同质区域步骤，406是分裂兴趣区步骤，407是更新兴趣区步骤，out4是输出候选种子点步骤。

步骤in4：输入兴趣区r；

步骤401：计算兴趣区域r的最大灰度值vmax和最小灰度值vmin；

步骤402：计算兴趣区域r的灰度直方图h，并用h(i)表示灰度值为i的像素点的数量，其中，i∈[0,k-1]，k为兴趣区域r所包含的灰度等级；

步骤403：将兴趣区域r中灰度值等于max(h(i))的像素点添加到候选种子点集合seedset中，max为取最大值函数；

步骤404：遍历候选种子点集合seedset中的所有聚类种子点s，用下式计算兴趣区域r中的灰度变分值sd：

式中，vs为候选种子点s的灰度值，gi为兴趣区域r中第i个灰度等级，i∈[0,k-1]；

步骤405：判断兴趣区域r的同质特性，当sd大于th时，表示兴趣区域r是同质区域，进入步骤out4，否则进入步骤406，其中，th为同质因子；

步骤406：分裂兴趣区，当vs-f×sd≤gi≤vs+f×sd时，删除灰度值属于区间[vs-f×sd,vs+f×sd]的像素点，并将剩余的像素点分裂为两个兴趣区：即[vmin,vs-f×sd]和[vs+f×sd,vmax]，其中，f为灰度变分值调整系数；

步骤407：将兴趣区域r分别更新为[vmin,vs-f×sd]和[vs+f×sd,vmax]，迭代运行步骤401～步骤406；

步骤out4：输出候选种子点。