专利名称：基于非下采样轮廓波变换的遥感图像道路增强方法
技术领域：
本发明属于图像处理技术领域，涉及图像增强，特别是一种遥感图像道路增强方
法，适用于遥感图像分析和处理。
背景技术：
空间技术和信息技术的飞速发展为我们提供了大量的遥感数据，而且从遥感影像
中获取目标信息已成为现阶段空间信息更新的重要手段。因此，如何智能解译海量的遥感
数据也相应地成为了信息化建设过程中面临的重要问题。作为其中的研究热点，道路的检
测研究更是得到了广泛关注，在国民经济生产和军事目标侦察领域具有十分重要的理论和
现实意义。遥感图像道路检测是图像目标检测的一个重要内容，对于地图更新、目标识别、
影像匹配等方面具有重要的辅助意义，因而是图像目标检测的一个重要内容。 由于遥感图像的获取过程受多种因素的影B向，因此图像的视觉效果和清晰度等方
面很容易变差，影响道路等有效信息的检测。遥感图像增强是一种能改善图像质量、丰富图像信息量、加强图像判读和识别效果的图像处理方法。
已有的遥感图像增强技术主要有基于空间域和基于变换域的方法两大类。空间域
的方法主要包括灰度变换、空域滤波和直方图处理、形态学处理等方法。基于变换域的方法
是将图像由时域变换到变换域，再通过修正变换域内的系数达到图像增强的目的。 近几年，多尺度几何理论的发展为图像处理提供了新的思路。多尺度几何变换能
够从多方向、多分辨等角度对图像中的轮廓和纹理等方向信息进行稀疏表示，并且对图像
中的线状奇异性具有很强的表征能力。非下采样轮廓波以其良好的性能，获得了较为广泛
的应用。在非下采样轮廓波变换的基础上进行图像去噪、融合、增强、编码等方法均取得了
较好的效果。 但目前常见的遥感图像的增强方法都仅仅利用空域或者频域的信息，增强后的道路失真会比较大，影响后续道路的检测工作，致使道路的定位不准确，道路目标识别不够准确。例如仅利用形态学增强的方法会使原有的道路宽度变宽，致使道路中心线的定位会发生偏移；而仅利用频域，例如小波域、轮廓波域，则会丢失掉道路的一些细节信息，致使道路断裂较多，道路不够完整。
发明内容
本发明的目的在于克服上述已有技术的缺陷，提出一种基于非下采样轮廓波变换
的遥感图像道路增强方法，以减小道路增强时的失真，实现对道路的准确定位。 本发明的技术方案是在非下采样轮廓波的多尺度性，多方向性，平移不变性的基
础上，通过分析非下采样轮廓波变换对遥感图像中道路和背景像素的系数响应的差异，设
计出遥感图像道路增强的方法。其步骤包括 1)对输入图像X进行非下采样的轮廓波变换，其中每层变换的方向个数由低到高排列分别分为4，4，8，并设轮廓波变换后的第d层的任意一个方向的系数为Xd(m， n);
2)根据各层变换的方向个数，设定相应的结构元素sei(i = 1，2，…，r)的数目；
3)将步骤l)的Xd(m，n)采用与之相近方向的结构元素进行方向增强，得到第d层的任意一个方向的增强后的系数Id(m，n); 4)对任意一个方向的增强系数Id(m， n)中的噪声和背景进行抑制，通过计算输入图像X中的像素点(m，n)的方向D(m，n)判断道路像素，如果D(m，n) 二8;则判Id(m，n)为背景或者噪声，否则将其判为道路像素； 5)重复步骤3)到步骤5)，对每一层的各个方向上的系数进行增强；
6)对所有的增强系数进行轮廓波反变换，得到图像的增强结果。
所述的对xd(m， n)进行增强，其步骤是 首先，采用与Xd(m，n)的方向最相近的两个结构元素Sei和se」分别对Xd(m，n)进行形态学开运算，得到Xd(m， n)在与i和j两个不同方向上开运算后的图像yd,i(m， n)和yd, j (m， n);然后，将Xd(m，n)分别与yd,i(m，n)和yd,j(m，n)进行差运算，得到xd(m，n)与yd,i(m，n)的差值图像Zd, i (m， n)和xd(m， n)与yd,」(m， n)的差值图像Zd,」(m， n);
最后，将差值图像Zd,i(m，n)和Zd,j(m，n)通过下式运算，得到增强后的系数
Id(m， n) = xd(m， n)+Zd(m， n)，其中Zd(m，n) = K* (Zd, i (m， n)+Zd,」(m， n)) ， K是常数，取经验值3 5。
本发明与现有技术相比具有如下优点 1)本发明由于采用非下采样轮廓波变换，在每个变换域方向上，均采用与之方向最相近的两个结构元素对其进行形态学开运算，因而能够更好的将图像的非下采样轮廓波变换域特征和空域特征更好的结合起来，与单纯的采用空域特征进行增强相比较，本发明能在增强道路的同时，很好地保持道路的原有宽度，使后续道路检测以及定位能够更加精确。 2)本发明由于采用图像中每个像素点的方向来区分道路和背景，能够很好的对道路进行增强，有效的抑制噪声和背景，达到较好的效果； 3)本发明流程简单，且不需要复杂的参数设置，实现过程简单快速，复杂度低。
图1是本发明的实现流程框图； 图2是本发明的9X9窗口示意图； 图3是本发明输入的两幅原始遥感图像； 图4是对图2进行轮廓波域增强的结果图； 图5是对图2进行空域增强结果图； 图6是本发明与现有空域增强的效果对比图； 图7是本发明增强前后的道路结果对比图。
具体实施例方式参照图l，本发明的实施如下 步骤l，对输入图像X进行非下采样的轮廓波变换。
对输入图像X进行非下采样的轮廓波变换，其中每层变换的方向个数由低到高排
列分别分为4，4，8，并设轮廓波变换后的第d层的任意一个方向的系数为xd(m， n)。 步骤2，根据各层变换的方向个数，设定相应的结构元素sei(i = 1，2，…，r)的数目。 对于轮廓波变换的8个方向，设定8个[se。， se!， se2， se3， se4， se5， se6， se7]结构
元素，对于轮廓波变换的4个方向，对应取结构元素[se。， se2， se4， see]，每个结构元素sei
均为一个大小为9X9的矩阵，它们将平面大致分为8个方向。 步骤3，对变换后的轮廓波系数xd(m， n)进行形态学方向增强。 3a)采用与xd(m， n)的方向最相近的两个结构元素和se」分别对xd(m， n)进
行形态学开运算，得到Xd(m， n)在与i和j两个不同方向上开运算后的图像yd,i(m， n)和
yd, j (m， n): yd,i(m， n) = xd(m， n) o sei， yd,j(m， n) = xd(m， n) o sej，其中o表示开运算符号；
3b)将xd(m， n)与yd,丄(m， n)进行差运算，得到xd(m， n)与yd,丄(m， n)的差值图像Zd,i(m，n);将Xd(m，n)与yd,j(m，n)进行差运算，得到的差值图像Zd,」(m， n):
Zd, iOn'rO = xd (m, n)-yd, i (m, n)，
Zd, j (m， n) = xd (m， n) -yd, j (m， n); 3c)将差值图像Zd,i(m， n)与Zd,j(m， n)进行和运算，得到增强后的系数
Id0ii， n) = xd0ii， n)+Zd(m， n)其中Zd(m，n) = K* (Zd, i (m， n)+Zd, j (m， n)) ， K是常数，取经验值3 5。
步骤4，对增强系数Id(m， n)中的噪声和背景进行抑制。 4a)以该像素点为中心取一个9X9的窗口，在结构元素sei所代表的8个方向上，分别计算该像素两侧约9个像素的灰度平均值M工和M2 ; 4b)根据灰度平均值M工和M2分别计算8个方向的像素对比度Vi(i = 0，1，…，7):
Vi = max(Vi,p Vi,2) 其中= (M/M》i， Vi,2 = (M2/M》j，即每个方向的像素的对比度是M:和M2比值M/M2及其倒数中的较大值； 4c)计算8个方向中对比度Vi(i = O，l，…，7)的最大值Ci = maX(Vi)， (i =0，
1， ".，7); 4d)设定阈值T，若Ci《T，则判该像素点为无方向，此时设定D(m，n) = 8 ;若Q >T，则结构元素sei代表的方向就是该像素点的方向，此时设定D(m， n) =i(i=0，l，，7)。 步骤5，对任意一个方向的增强系数Id(m， n)中的噪声和背景进行处理。[OO54] 计算输入图像X中像素点(m，n)的方向D(m，n)判断道路像素，如果D(m，n) = 8 ;则判Id(m，n)为背景或者是噪声，将Id(m，n)的值变为0，否则将像素(m，n)判为道路像素，保持Id(m，n)的值不变。 步骤6，重复步骤3到步骤6，对每一层的8个方向上的系数均进行方向增强，并对
噪声和背景进行处理，得到输入图像每层每个方向的处理后的增强系数。 步骤7，对所有处理后的增强系数进行轮廓波反变换，得到遥感图像的最终增强结果。
实验结果与分析 图3是输入的两幅原始遥感图像，白色方框标示出了待评价图像的局部区域。
图4是采用本发明的方法进行增强的结果，图3(a)和图3(b)对应的增强结果分 别是图4(a)和图4(b)。 图5是采用现有的空域方法进行增强的结果，图3(a)和图3(b)对应的增强结果 分别是图5(a)和图5(b)。 图6(a)是图3(a)的局部放大图，图6(b)是图4(a)的局部放大图，图6(c)是图 5(a)的局部放大图。将6(a)、6(b)和6 (c)进行对比从图6(b)的白色方框区域中，可以 很明显的看出本发明的方法可以很好的将原始图像中的道路很好的增强，使道路更加明显 区别与背景，同时不会明显改变道路的原有宽度，而图6(c)将道路的宽度明显增大。
图7(a)是采用图3(a)作为输入图像进行道路检测的结果，图7 (c)是采用图4(a) 作为输入图像进行道路检测结果，图7(e)是采用图5(a)作为输入图像进行道路检测结果， 图7(b)是图7(a)叠加到图3(a)上的结果，图7(d)是图7(c)叠加到图3(a)上的结果，图 7(f)是图7(e)叠加到图3(a)上的结果。 从图7(a)和图7(d)中可以很明显的看出，采用未增强的图像，在道路的检测结果 中，道路线存在有明显的断裂，不能形成明显的道路网，且有很多分支道路的检测结果和原 始图像中的道路线不一致，有较多的漏检和误检。 从图7(b)和图7(e)中可以看出，本发明的方法检测的道路较完整，且漏检和误检 较少，检测出的道路也基本和实际的道路线重合。 从7(c)和图7(f)中可以看出，检测道路与实际的道路线存在一定的偏移，并且丢 失掉了道路的一些细节信息。 综上所述，本发明将遥感图像的空域特征和非下采样轮廓波变换域特征有效的结 合起来，使得道路的增强结果更加明显；且不需要复杂的参数设置，实现过程简单快速，复 杂度低，与空域增强相比较，本发明方法在增强道路的同时，能很好的保持道路原有的宽 度，使后续道路的检测能够更加准确。
权利要求
一种基于非下采样轮廓波变换的遥感图像道路增强方法，包括如下步骤1)对输入图像X进行非下采样的轮廓波变换，其中每层变换的方向个数由低到高排列分别分为4，4，8，并设轮廓波变换后的第d层的任意一个方向的系数为xd(m，n)；2)根据各层变换的方向个数，设定相应的结构元素sei(i＝1，2，…，r)的数目；3)将步骤1)的xd(m，n)采用与之相近方向的结构元素进行方向增强，得到第d层的任意一个方向的增强后的系数Id(m，n)；4)对任意一个方向的增强系数Id(m，n)中的噪声和背景进行抑制，通过计算输入图像X中的像素点(m，n)的方向D(m，n)判断道路像素，如果D(m，n)＝8；则判Id(m，n)为背景或者噪声，否则将其判为道路像素；5)重复步骤3)到步骤5)，对每一层的各个方向上的系数进行增强；6)对所有的增强系数进行轮廓波反变换，得到图像的增强结果。
2. 根据权利要求l的遥感图像道路增强方法，其中步骤2)所述的结构元素sei(i = 1， 2，…，r)为一个大小为9X9的矩阵，它们将平面大致分为8个方向，每一个结构元素代表 其中的一个方向。
3. 根据权利要求1的遥感图像道路增强方法，其中步骤3)所述的对Xd(m， n)进行增 强，按如下步骤进行3a)采用与Xd(m，n)的方向最相近的两个结构元素Sei和Sei分别对Xd(m，n)进行形态 学开运算，得到Xd(m， n)在与i和j两个不同方向上开运算后的图像yd,i(m， n)和yd,j(m， n);3b)将xd (m， n)分别与yd, t (m， n)和yd,」(m， n)进行差运算，得到xd (m， n)与yd, t (m， n) 的差值图像Zd, i (m， n)和xd (m， n)与yd,」(m， n)的差值图像Zd,」(m， n);3c)将差值图像Zd,i(m， n)和Zd,j(m， n)通过下式运算，得到增强后的系数 Id (m， n) = xd (m， n) +Zd (m， n)，其中Zd (m， n) = K* (Zd, i (m， n) +Zd,」(m， n)) ， K是常数，取经验值3 5 。
4. 根据权利要求l的遥感图像道路增强方法，其中步骤4)所述的计算输入图像X中的 像素点(m， n)的方向D(m， n)，按如下步骤计算:4a)以该像素点为中心取一个9X9的窗口，在结构元素所代表的8个方向上，分别计算 该像素两侧约9个像素的灰度平均值M工和M2 ;4b)根据灰度平均值M工和M2分别计算8个方向的像素对比度Vi(i =0，1，…，7):<formula>formula see original document page 2</formula>4c)计算所有方向中对比度Vi(i = O，l，…，7)的最大值Ci = maX(Vi)， (i = O，l，…，7);4d)设定阈值T，若Ci《T，则判该像素点为无方向，设定D(m， n) = 8 ;若& > T，则结 构元素sei代表的方向就是该像素点的方向，设定D(m，n) =i(i=0，l，…，7)。
全文摘要
本发明公开了一种遥感图像道路增强的方法，主要解决现有技术增强后的道路失真大，道路目标检测不准确的问题。其实现过程是首先对遥感图像进行3层非下采样的轮廓波变换，其中每层变换的方向个数由低到高排列分别分为4，4，8，再根据各层变换的方向个数，设定相应的结构元素；然后对变换后的系数采用与之相近方向的结构元素进行方向性增强；最后计算图像中每个像素点的方向，得到图像的方向矩阵，通过方向矩阵对增强系数中的噪声和背景进行处理，再对处理后的增强系数进行轮廓波反变换，得到图像的增强结果。本发明能在增强道路的同时，很好的保持道路原有的宽度，使后续道路的检测能够更加准确，可用于遥感图像道路的分析和处理。
文档编号G06T5/00GK101739667SQ20091021932
公开日2010年6月16日 申请日期2009年12月4日 优先权日2009年12月4日
发明者侯彪, 冯颖涛, 张小华, 杨淑媛, 焦李成, 王桂婷, 王爽, 钟桦 申请人:西安电子科技大学