[发明专利]多光谱遥感影像矢量边缘检测方法

说明书

技术领域

本发明涉及一种多光谱遥感影像矢量边缘检测方法。

背景技术

边缘检测一般借助图像中目标的某种不连续性特征(如亮度、色相、饱和度等)来实现边缘信息地提取，在遥感影像地物目标提取中有着广泛的应用。在多光谱高空间分辨率遥感影像中，由于地物材质、方位、几何形状和光照条件的不同，会存在反射边缘、朝向边缘、遮挡边缘、照明(阴影)边缘，以及镜面(高光)边缘；它们使边缘提取结果中通常存在大量伪边缘，严重影响后续基于边缘的图像分析和理解。

著名的边缘检测Canny算子以其良好的边缘检测效果闻名于世，但该算子仅能应用于单通道灰度图像的边缘信息提取，应用于多通道遥感影像边缘检测时则需要进行必要的扩展。相对单波段(单通道)遥感影像而言，多波段(多通道或多光谱，一般波段数目不小于三个)遥感影像中可用于边缘检测的光谱信息及特征更加丰富，所采用的特征与提取方案需要根据地物类别灵活地确定。

传统边缘检测方法(如Canny算子)应用于多光谱遥感影像时存在的主要问题有：

(1)波段数目和处理策略的限制。传统边缘检测方法一般仅能处理单波段灰度影像(或采用逐波段对多波段影像分别进行处理的策略)，各波段对应地物目标的检测响应程度(一般会受波谱范围差异的影响)及之间的关系特征在检测过程中却很少给予考虑，割裂了多光谱遥感影像各波段数据间的有机联系，造成边缘检测结果的局部缺失，无法针对各类地物进行有效的边缘信息提取。

(2)传统方法对多光谱高空间分辨率遥感影像数据的适宜性问题。多光谱高空间分辨率遥感影像数据具有以下特点，数据量大、地物几何与属性细节信息丰富，同类地物甚至同一地物内部光谱异质性较高、“同物异谱”与“异物同谱”现象普遍，地物目标空间结构格局复杂、边界过渡区繁多，阴影与细小地物干扰严重等。这些问题导致利用传统方法获取的边缘检测结果中通常存在大量伪边缘，甚至有时无法对多类型地物目标边缘同时进行有效提取及分层(类)分析与识别，严重影响基于边缘信息的多光谱遥感影像分析与理解等后续工作的有效开展。

(3)传统边缘检测方法对各波段遥感影像对应的边缘检测结果缺乏相关性分析。多波段高空间分辨率遥感影像是各类地物光谱特征在不同光谱范围内的差异化记录，不同地物在不同波段的光谱响应程度差异明显，显然各波段影像将产生不同的边缘检测结果。因此各波段边缘检测结果缺乏相关性分析将导致边缘信息缺失以及检测精度下降。

发明内容

针对上述问题，本发明的目的在于提供一种多光谱遥感影像矢量边缘检测方法。

为达到上述目的，本发明所述一种多光谱遥感影像矢量边缘检测方法，包括以下步骤：

获取具有n个波段的高空间分辨率遥感影像；

对高空间分辨率遥感影像进行梯度矢量边缘检测获得梯度矢量特征影像；

获取梯度矢量特征影像上的最大变化或不连续变化位置点的角度值以及对应位置点在该角度方向上的变化率；

根据获得的角度值以及变化率利用非极值抑制法确定梯度矢量特征影像的初始边缘点集；

根据初始边缘点集利用双阈值分割处理法进行边缘轮廓检测获得影像边缘轮廓；

对边缘轮廓的各波段在多维色彩空间里叠置合成边缘矢量，并输出边缘检测结果。

优选地，根据获得的角度值以及变化率利用非极值抑制法确定梯度矢量特征影像的初始边缘点集的步骤包括：

逐波段选择n维梯度影像中所有像元点，对梯度值较大像元所形成的屋脊带用非极值抑制法进行细化；

若在该像元方向角度θ上的梯度值是局部最大值，则保留为初步边缘点；

反之，将该像元设置为非边缘点。

优选地，根据初始边缘点集利用双阈值分割处理法进行边缘轮廓检测获得影像边缘轮廓的步骤包括：

选定两个梯度阈值Y_h和Y_s；

逐波段在非极值抑制结果中去除梯度值小于Y_h的像元点，并得到强边缘点集Q；

以Q为基础把边缘点连接成初始轮廓；

对初始轮廓进行搜索，在梯度值介于Y_h与Y_s的非极值抑制结果中寻找可以连接到当前端点的边缘点；

利用递归跟踪方法在介于Y_h与Y_s的梯度值中搜集边缘，直到将Y_h中所有间断相连接。

本发明的有益效果为：