基于颜色恒常性特征和视觉注意机制的阴影检测方法技术

本发明专利技术公开了一种基于阴影的多种颜色不变性特征，引入计算机视觉领域的视觉注意机制，提出了一种新的有效的阴影检测方法，具有以下有益效果：1）根据颜色恒常性研究成果，提取了阴影的多种颜色不变性特征，采用视觉注意模型中的特征显著图模型，将多种阴影特征结合在一起，生成阴影显著图作为阴影检测的基础数据，然后基于此数据进行阴影检测；2）采用SLIC算法对图像进行超像素分割，提出基于对象的自适应阈值生长算法，对阴影检测的初步结果进行修正，从而获得更完整的阴影检测结果。

【技术实现步骤摘要】

本专利技术涉及一种根据彩色图像上阴影的颜色恒常性特征和视觉注意机制的阴影检测方法，属于数字图像处理和目标识别
，可用于彩色照片和遥感影像上的地物阴影检测。
技术介绍
地物阴影是自然照片、视频以及遥感影像上不可避免的现象，但阴影的存在会给计算机视觉的任务带来很多困难，往往造成图像分割、边缘检测和对象识别的错误。因而，阴影检测是目前图像处理和计算机视觉领域的研究热点，也是一个难点问题，在阴影检测中，主要分为基于视频的动态图像阴影检测和基于单幅图像的静态图像阴影检测两个大类。由于基于动态图像的阴影检测可以充分利用相邻图像帧之间的相关信息，相对容易，研究文献也相对较多，但基于静态图像的阴影检测难度更大，研究相对较少。目前，对静态图像阴影的检测主要分为四类方法：1)、基于物理模型的方法：该方法主要思想在于先对阴影成像过程进行建模，进而设计检测算法；2)、基于非物理模型的方法：主要是采用数学模型和人类视觉模型对阴影进行建模和检测，如较常用的有Retinex模型、Sigmod模型、Bayes统计模型；3)、基于本征图像的阴影检测方法：该类方法通常先将原始图像变换到光照不敏感空间，如C1C2C3空间、HSV颜色空间，然后在新的空间设计检测算法，该类方法算法简单，但检测结果依赖于本征图像的质量；4)、基于统计学习的阴影检测方法：先对图像阴影进行特征提取，然后利用所提取的特征对分类器进行训练，最后再通过训练好的分类器对图像进行分类从而获得阴影结果。在每类方法中，研究者都提出了不同的思路和具体实现算法，根据研究结果显示，每类方法都有各自的优点和缺陷，因而，到目前为止还没有一个鲁棒而通用的阴影检测方法。鉴于此，研究新的阴影检测理论和方法仍然非常必要。
技术实现思路
为解决现有技术的不足，本专利技术的目的在于提供一种基于阴影的多种颜色不变性特征和视觉注意机制的阴影检测方法，用于改善阴影检测的结果。为了实现上述目标，本专利技术采用如下的技术方案：基于颜色恒常性特征和视觉注意机制的阴影检测方法，其特征在于，包括如下步骤：(1)、对原始图像进行亮度特征提取，记为Lum；(2)、将原始图像转换到C1C2C3颜色不变空间，提取出其中适合于阴影检测的第三个通道，即C3成分，记为C3；(3)、基于赫林拮抗色原理，对原始图像进行STB模型变换，提取其中的蓝黄对比特征图像，记为B_Y；(4)、分别对上述三个特征图像Lum、C3、B_Y进行局部最大值归一化处理，得到归一化的特征图Lum’、C3’、B_Y’；(5)、根据归一化特征图像计算阴影显著图，记为S_Sal，计算公式为：S_Sal＝Lum'*C3'*B_Y'；(6)、采用大津阈值算法OSTU对S_Sal进行自适应阈值二值化，获得初始阴影分割图S0，其中值1表示为阴影，0表示为非阴影；(7)、对原始图像进行SLIC超像素分割，获得分割图像SegImg，将其中每个分割区域分别记为SegImg[k],k＝0,…,n，共获得n个区域；(8)、将初始阴影分割图S0和SLIC分割结果SegImg进行叠加，然后对SegImg分割结果进行标记，具体如下：a、依次取出S0分割图中值为1的像素，假设该像素位于第y行、第x列，则记为S0[y,x]；b、若S0[y,x]位于分割图RegImg的第k个区域RegImg[k]，则将RegImg[k]标记为阴影区(value＝1)，否则将RegImg[k]标记为非阴影区(value＝0)；(9)、对SegImg中的每个分割区域SegImg[k]，计算其颜色均值m[k]和重心位置c[k]，其中m[k]、c[k]分别为对应的颜色向量和坐标向量：m[k]＝[rk,gk,bk]c[k]＝[xk,yk]；其中，rk、gk、bk分别为第k个分割区域的红、绿、蓝颜色均值，xk、yk分别为第k个分割区域的重心的x坐标和y坐标；(10)、根据分割图像RegImg中所有分割区域的重心位置坐标c[k],k＝0,…,n，构建Delaunay三角网，根据Delaunay三角网信息可获得每个区域RegImg[k]的邻接区域；(11)、计算所有具有邻接关系的阴影区域两两之间的颜色距离，如第i个区域和第j个区域之间的颜色距离为：D=(ri-rj)2+(gi-gj)2+(bi-bj)2]]>将所有两两邻接阴影区域的颜色距离值存入列表DIS_L[k],k＝0,…,m；其中ri、gi、bi分别为第i个分割区域对应的红、绿、蓝颜色分量，rj、gj、bj分别为第j个分割区域对应的红、绿、蓝颜色分量；(12)、对RegImg中标记为1的阴影区域逐个进行生长算法；(13)、整个阴影检测算法结束，分割图像RegImg中标记为1的区域为阴影检测结果。优选地，前述步骤(1)中，所述原始图像为包含R,G,B三个波段的彩色图像，具体过程为：先提取原始图的灰度图，然后将灰度图的对数倒数作为lum特征，计算公式为：Gray＝(R+G+B)/3Lum＝1.0/[5+log(Gray+1)]其中，Gray表示3个颜色分量R、G、B的均值，即灰度值；第2个公式中分母中的1是为了避免对0求对数的情况，从而保证数学运算的有效性；5是调节因子。更优选地，前述步骤(2)中，C1、C2、C3的具体计算公式分别为：C1＝arctan(G/max(R,B))C2＝arctan(R/max(G,B))C3＝arctan(B/max(R,G))。再优选地，前述步骤(3)中，具体计算公式为：MaxC＝max(R,G,B)MinC＝min(R,G)B_Y＝(B-MinC)/MaxC其中，MaxC表示3个颜色分量R、G、B中的最大值；MinC表示颜色分量R、G中的较小值。进一步优选地，前述步骤(4)中，假设待归一化的特征图像为ftmap，对其进行局部最大值归一化的伪代码如下：M＝10T＝M*thresholdftmap＝M*Normalize(ftmap)peaks＝FindLocalMax(ftmap)num_peaks＝Number(peaks)if(num_peaks＝＝1):ftmap＝ftmap*M*Melseif(num_peaks>1):locAvg＝Mean(peaks)ftmap＝ftmap*(M-locAvg)*(M-locAvg)其中，threshold通常取值为0.1；Normalize()函数的作本文档来自技高网...

【技术保护点】
基于颜色恒常性特征和视觉注意机制的阴影检测方法，其特征在于，包括如下步骤：(1)、对原始图像进行亮度特征提取，记为Lum；(2)、将原始图像转换到C1C2C3颜色不变空间，提取出其中适合于阴影检测的第三个通道，即C3成分，记为C3；(3)、基于赫林拮抗色原理，对原始图像进行STB模型变换，提取其中的蓝黄对比特征图像，记为B_Y；(4)、分别对上述三个特征图像Lum、C3、B_Y进行局部最大值归一化处理，得到归一化的特征图Lum’、C3’、B_Y’；(5)、根据归一化特征图像计算阴影显著图，记为S_Sal，计算公式为：S_Sal＝Lum'*C3'*B_Y'；(6)、采用大津阈值算法OSTU对S_Sal进行自适应阈值二值化，获得初始阴影分割图S0，其中值1表示为阴影，0表示为非阴影；(7)、对原始图像进行SLIC超像素分割，获得分割图像SegImg，将其中每个分割区域分别记为SegImg[k],k＝0,…,n，共获得n个区域；(8)、将初始阴影分割图S0和SLIC分割结果SegImg进行叠加，然后对SegImg分割结果进行标记，具体如下：a、依次取出S0分割图中值为1的像素，假设该像素位于第y行、第x列，则记为S0[y,x]；b、若S0[y,x]位于分割图RegImg的第k个区域RegImg[k]，则将RegImg[k]标记为阴影区(value＝1)，否则将RegImg[k]标记为非阴影区(value＝0)；(9)、对SegImg中的每个分割区域SegImg[k]，计算其颜色均值m[k]和重心位置c[k]，其中m[k]、c[k]分别为对应的颜色向量和坐标向量：m[k]＝[rk,gk,bk]c[k]＝[xk,yk]；其中，rk、gk、bk分别为第k个分割区域的红、绿、蓝颜色均值，xk、yk分别为第k个分割区域的重心的x坐标和y坐标；(10)、根据分割图像RegImg中所有分割区域的重心位置坐标c[k],k＝0,…,n，构建Delaunay三角网，根据Delaunay三角网信息可获得每个区域RegImg[k]的邻接区域；(11)、计算所有具有邻接关系的阴影区域两两之间的颜色距离，如第i个区域和第j个区域之间的颜色距离为：D=(ri-rj)2+(gi-gj)2+(bi-bj)2]]>将所有两两邻接阴影区域的颜色距离值存入列表DIS_L[k],k＝0,…,m；其中ri、gi、bi分别为第i个分割区域对应的红、绿、蓝颜色分量，rj、gj、bj分别为第j个分割区域对应的红、绿、蓝颜色分量；(12)、对RegImg中标记为1的阴影区域逐个进行生长算法；(13)、整个阴影检测算法结束，分割图像RegImg中标记为1的区域为阴影检测结果。...

【技术特征摘要】
1.基于颜色恒常性特征和视觉注意机制的阴影检测方法，其特征在于，包括如下步骤：
(1)、对原始图像进行亮度特征提取，记为Lum；
(2)、将原始图像转换到C1C2C3颜色不变空间，提取出其中适合于阴影检测的第三个通
道，即C3成分，记为C3；
(3)、基于赫林拮抗色原理，对原始图像进行STB模型变换，提取其中的蓝黄对比特征图
像，记为B_Y；
(4)、分别对上述三个特征图像Lum、C3、B_Y进行局部最大值归一化处理，得到归一化的
特征图Lum’、C3’、B_Y’；
(5)、根据归一化特征图像计算阴影显著图，记为S_Sal，计算公式为：
S_Sal＝Lum'*C3'*B_Y'；
(6)、采用大津阈值算法OSTU对S_Sal进行自适应阈值二值化，获得初始阴影分割图S0，
其中值1表示为阴影，0表示为非阴影；
(7)、对原始图像进行SLIC超像素分割，获得分割图像SegImg，将其中每个分割区域分
别记为SegImg[k],k＝0,…,n，共获得n个区域；
(8)、将初始阴影分割图S0和SLIC分割结果SegImg进行叠加，然后对SegImg分割结果进
行标记，具体如下：
a、依次取出S0分割图中值为1的像素，假设该像素位于第y行、第x列，则记为S0[y,x]；
b、若S0[y,x]位于分割图RegImg的第k个区域RegImg[k]，则将RegImg[k]标记为阴影区
(value＝1)，否则将RegImg[k]标记为非阴影区(value＝0)；
(9)、对SegImg中的每个分割区域SegImg[k]，计算其颜色均值m[k]和重心位置c[k]，其
中m[k]、c[k]分别为对应的颜色向量和坐标向量：
m[k]＝[rk,gk,bk]
c[k]＝[xk,yk]；
其中，rk、gk、bk分别为第k个分割区域的红、绿、蓝颜色均值，xk、yk分别为第k个分割区域
的重心的x坐标和y坐标；
(10)、根据分割图像RegImg中所有分割区域的重心位置坐标c[k],k＝0,…,n，构建
Delaunay三角网，根据Delaunay三角网信息可获得每个区域RegImg[k]的邻接区域；
(11)、计算所有具有邻接关系的阴影区域两两之间的颜色距离，如第i个区域和第j个
区域之间的颜色距离为：
D=(ri-rj)2+(gi-gj)2+(bi-bj)2]]>将所有两两邻接阴影区域的颜色距离值存入列表DIS_L[k],k＝0,…,m；其中ri、gi、bi分别为第i个分割区域对应的红、绿、蓝颜色分量，rj、gj、bj分别为第j个分割区域对应的红、
绿、蓝颜色分量；
(12)、对RegImg中标记为1的阴影区域逐个进行生长算法；
(13)、整个阴影检测算法结束，分割图像RegImg中标记为1的区域为阴影检测结果。
2.根据权利要求1所述的基于颜色恒常性特征和视觉注意机制的阴影检测方法，其特
征在于，步骤(1)中，所述原始图像为包含R,G,B三个波段的彩色图像，具体过程为：先提取
原始图的灰度图，然后将灰度图的对数倒数作为lum特征，计算公式为：
Gray＝(R+G+B)/3
Lum＝1.0/[5+log(Gray+1)]...

【专利技术属性】
技术研发人员：陈仁喜，
申请(专利权)人：河海大学，
类型：发明
国别省市：江苏;32