[发明专利]基于非下采样轮廓波变换的遥感图像道路提取方法

说明书

技术领域

本发明属于图像处理技术领域，涉及图像目标检测，特别是一种遥感图像道路提取方法，适用于遥感图像分析和处理。

背景技术

空间技术和信息技术的飞速发展为我们提供了大量的遥感数据，而且从遥感影像中获取目标信息己成为现阶段空间信息更新的重要手段。因此，如何智能解译海量的遥感数据也相应地成为了信息化建设过程中面临的重要问题。作为其中的研究热点，道路的提取研究更是得到了广泛关注，在国民经济生产和军事目标侦察领域具有十分重要的理论和现实意义。遥感图像道路提取是图像目标检测的一个重要内容，对于地图更新、目标识别、影像匹配等方面具有重要的辅助意义，因而是图像目标检测的一个重要内容。

在图像处理领域，形状是目标的一种本质特征，基于形状的目标识别是计算机目标识别中一种重要的技术方法。在各种形状的目标中，线状目标占有很大一部分比例，因而研究线状目标的自动识别具有很重要的理论意义和实用价值。道路是典型的线状目标，因此在处理遥感图像的过程中，对线性地物尤其是对道路的提取占有十分基础和重要的位置。

目前，国内外的专家们提出了许多道路提取的方法，尤其是近年来国内外的研究者们在不同的研究领域提出了许多理论和方法，推动了遥感图像道路网的自动提取技术的巨大进步。线性特征提取方法分为空域的方法和变换域的方法，前者包括基于模板的方法、阈值的方法、形态学的方法以及基于蛇模型的方法等等；后者包括基于小波、小波包的方法，其中Gabor小波在指纹提取、视网膜血管提取等领域应用非常广泛，但是由于Gabor小波基函数的非正交性、倍频性不够显著等缺点，对道路奇异性的描述存在本质的不足。Gabor小波的另一个缺点是参数调整比较复杂，因而限制了其应用。

近几年，多尺度几何理论的发展为图像处理提供了新的思路。多尺度几何变换能够从多方向、多分辨等角度对图像中的轮廓和纹理等方向信息进行稀疏表示，并且对图像中的线状奇异性具有很强的表征能力。非下采样轮廓波以其良好的性能，获得了较为广泛的应用。在非下采样轮廓波变换的基础上进行图像去噪、融合、增强、编码等方法均取得了较好的效果。但将非下采样轮廓波变换很好的用于道路提取还是一个挑战，目前还没有将此变换用于遥感图像道路提取的相关技术。

发明内容

本发明的目的在于克服已有技术的缺陷，提出一种基于非下采样轮廓波变换的遥感图像道路提取方法，以剔除虚假目标，提高道路的连续性，实现对道路的准确定位。

本发明的技术方案是在非下采样轮廓波的多尺度性，多方向性，平移不变性的基础上，通过分析非下采样轮廓波变换对遥感图像中道路和背景像素的系数响应的差异，设计出遥感图像道路检测的方法。其提取的步骤包括：

1)对输入图像X进行包括Frost去噪和自适应直方图均衡化的预处理；

2)对预处理后的图像进行3个尺度非下采样轮廓波变换，每个尺度变换为8个方向，提取第1个尺度和第2个尺度各方向子带的模极大值作为道路的特征向量；

3)对提取的线性特征向量进行模糊C均值聚类；

4)对初始分割结果进行包括非极大值抑制以及基于空间关系的道路后处理，得到遥感图像道路。

上述步骤4)所述的对初始分割结果进行包括非极大值抑制以及基于空间关系的道路后处理，按如下步骤进行：

(4a)用导数算子求出初始分割结果图的梯度大小和方向；

(4b)把梯度分为水平、竖直、45°方向和135°四个方向，将各个方向与邻近像素进行比较，在图像中某个像素的梯度方向上，若该像素的灰度值与其前后两个像素的灰度值相比不是最大的，则判定其为非道路像素，并将该像素的灰度值设定为0，否则判定其为道路像素，并将该像素的灰度值设定为1，得到道路集合ROAD；

(4c)计算道路集合ROAD中每个道路的长度length，剔除掉小于设定阈值的片段，得到最终的道路检测结果road。

本发明与现有技术相比具有如下优点

1)非下采样轮廓波变换在对线状奇异性表示方面优于其它变换域的方法，本发明通过对包含道路的遥感图像进行非下采样轮廓波变换，提取道路的线性特征向量，从而能够准确捕获图像中的高维奇异信息，很好的提取道路；并且非下采样轮廓波变换是一种固定框架、固定参数的变换，它能够有效避免其它方法中需要设置多个参数的问题，使检测的结果更加稳健。

2)本发明采用的模糊C均值聚类由于不需要训练学习，运算速度快，因此提高了检测的效率。

3)本发明由于采用对初始分割进行非极大值抑制，可以使检测出的道路单像素化易于后处理，同时便于将检测道路结果叠加到原始图像上，更好的评价检测的结果。