发明名称 一种基于时空多特征的车辆阴影消除方法
摘要 本发明公开了一种基于时空多特征的车辆阴影消除方法。本发明首先获取待处理视频帧图像的背景图像及初始前景区域,然后基于其获取色度、光谱方向和纹理特征的前景掩膜,并进行加权求和得到空域多特征前景概率谱图像,接着对空域多特征前景概率谱图像进行时间滑动窗的时域滤波,得到对应的时域前景概率谱,对时域前景概率谱进行加权处理后得到最终的前景掩膜,完成对视频帧图像的车辆阴影消除处理,本发明相较于传统基于单一特征和多特征级联的阴影去除方法,阴影去除更干净,前景轮廓更完整。在保障了较高的阴影辨识率的同时也维持较高的阴影检测率。
申请公布号 CN106339995A 申请公布日期 2017.01.18
申请号 CN201610764066.7 申请日期 2016.08.30
申请人 电子科技大学 发明人 王正宁;柏祁林;韩明燕;周阳;马姗姗
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 电子科技大学专利中心 51203 代理人 周刘英
主权项 一种基于时空多特征的车辆阴影消除方法,其特征在于包括下列步骤:步骤1:输入待处理视频帧I;步骤2:对当前视频帧I的背景进行建模,得到背景图像B并计算对应的初始前景区域F<sub>0</sub>;步骤3:分别计算色度、光谱方向、纹理三种特征的前景掩膜F<sub>chr</sub>,F<sub>phy</sub>,F<sub>tex</sub>;步骤301:计算色的特征的前景掩膜F<sub>chr</sub>。将视频帧I和背景图像B转换到色彩空间W,色彩空间W包括图像的亮度、色度信息;对初始前景区域F<sub>0</sub>中的任意像素(x,y)计算视频帧I和背景图像B中的亮度比<img file="FDA0001099773090000011.GIF" wi="356" he="135" />其中I<sup>V</sup>(x,y)、B<sup>V</sup>(x,y)表示视频帧I和背景图像B中像素(x,y)的亮度分量;以及计算视频帧I和背景图像B中以任意像素(x,y)为中心、k1为半径的邻域的色度分量差分均值<img file="FDA0001099773090000012.GIF" wi="601" he="127" />其中<img file="FDA0001099773090000013.GIF" wi="418" he="63" />表示视频帧I和背景图像B中像素(x,y)的邻域内所有像素的色度分量累加和;根据亮度比V<sub>ratio</sub>、C<sub>avg</sub>计算色度特征的的阴影掩膜S<sub>chr</sub>以及对应的前景F<sub>chr</sub>:对任意像素(x,y),若α≤V<sub>ratio</sub>≤β,且C<sub>avg</sub>≤τ<sub>C</sub>,则阴影掩膜S<sub>chr</sub>中对应像素值S<sub>chr</sub>(x,y)=1;否则S<sub>chr</sub>(x,y)=0,其中α、β分别表示亮度阈值的下限、上限,τ<sub>C</sub>表示色度阈值;对任意像素(x,y),若S<sub>chr</sub>(x,y)=0且F<sub>0</sub>(x,y)=1,则前景掩膜F<sub>chr</sub>中对应像素值F<sub>chr</sub>(x,y)=1;否则F<sub>chr</sub>(x,y)=0,其中F<sub>0</sub>(x,y)表示初始前景区域F<sub>0</sub>中像素(x,y)的像素值;步骤302:计算光谱方向特征的前景掩膜F<sub>phy</sub>:对于初始前景区域F<sub>0</sub>中的任意前景像素(x,y),计算(x,y)在视频帧I和背景图像B中RGB三个色彩通道的差分向量v(x,y)、以及v(x,y)与B(x,y)的模之比α(x,y)、v(x,y)的G分量v<sup>G</sup>(x,y)与R分量v<sup>R</sup>(x,y)的夹角θ(x,y)、以及B分量与v(x,y)模的夹角<img file="FDA0001099773090000014.GIF" wi="160" he="62" />形成光谱方向特征<img file="FDA0001099773090000018.GIF" wi="763" he="59" />其中<img file="FDA0001099773090000015.GIF" wi="414" he="150" /><img file="FDA0001099773090000016.GIF" wi="566" he="134" /><img file="FDA0001099773090000017.GIF" wi="550" he="143" />对初始前景区域F<sub>0</sub>中任意前景像素的光谱方向特征f(x,y)聚类得到阴影模型M<sub>s</sub>(μ,σ<sup>2</sup>),并计算基于光谱方向的阴影掩膜S<sub>phy</sub>以及前景掩膜F<sub>phy</sub>:对任意像素(x,y),若||f(x,y)‑μ||<σ<sup>2</sup>,则阴影掩膜S<sub>phy</sub>中对应像素值S<sub>phy</sub>(x,y)=1;否则S<sub>phy</sub>(x,y)=0,其中,μ,σ<sup>2</sup>分别表示M<sub>s</sub>(μ,σ<sup>2</sup>)的均值和方差;对任意像素(x,y),若S<sub>phy</sub>(x,y)=0且F<sub>0</sub>(x,y)=1,则前景掩膜F<sub>phy</sub>中对应像素值F<sub>phy</sub>(x,y)=1;否则F<sub>phy</sub>(x,y)=0;步骤303:计算纹理特征的前景掩膜F<sub>tex</sub>。将初始前景区域F<sub>0</sub>转换到色彩空间W;根据视频帧I和背景图像B中像素(x,y)的亮度分量I<sup>V</sup>(x,y)、B<sup>V</sup>(x,y),以及色度分量I<sup>C</sup>(x,y)、B<sup>C</sup>(x,y)计算初选阴影掩膜S<sub>W</sub>:若I<sup>V</sup>(x,y)<B<sup>V</sup>(x,y)且|I<sup>C</sup>(x,y)‑B<sup>C</sup>(x,y)|<T<sub>C</sub>,则初选阴影掩膜S<sub>W</sub>中对应像素值S<sub>W</sub>(x,y)=1;否则S<sub>W</sub>(x,y)=0,其中T<sub>C</sub>表示色度差阈值;计算频帧I和背景图像B的边缘图像I<sub>edge</sub>和B<sub>edge</sub>,并差分得到前景边缘像素集合F<sub>edge</sub>;根据F<sub>edge</sub>优化初选阴影S<sub>W</sub>得到边缘阴影掩膜S<sub>edge</sub>:对任意像素(x,y),若S<sub>hsv</sub>(x,y)=1且<img file="FDA0001099773090000024.GIF" wi="435" he="59" />则边缘阴影掩膜S<sub>edge</sub>中对应像素值S<sub>edge</sub>(x,y)=1;若S<sub>hsv</sub>(x,y)∈Ω<sub>k2</sub>(u,v),则S<sub>edge</sub>(x,y)=0,其中Ω<sub>k2</sub>(u,v)表示以F<sub>edge</sub>(u,v)为中心、k2为半径的邻域,F<sub>edge</sub>(u,v)表示前景边缘像素集合F<sub>edge</sub>中的任意像素(u,v);计算边缘阴影掩膜S<sub>edge</sub>中每个连通域的前景与背景的梯度纹理特征夹角<img file="FDA0001099773090000021.GIF" wi="966" he="166" />其中<img file="FDA0001099773090000025.GIF" wi="159" he="62" />分别表示水平和竖直方向的梯度,上标F和B分别用于区分前景与背景;统计S<sub>edge</sub>中每个连通域中纹理特征与背景相近的像素及其所占比例<img file="FDA0001099773090000022.GIF" wi="502" he="111" />其中n表示连通域像素总数,τ<sub>α</sub>表示相似像素的数量阈值,H(·)表示阶跃函数;计算纹理特征的阴影掩膜S<sub>tex</sub>对应的前景掩膜F<sub>tex</sub>:对任意像素(x,y),若d>T<sub>tex</sub>,则阴影掩膜S<sub>tex</sub>中对应像素值S<sub>tex</sub>(x,y)=1;否则S<sub>tex</sub>(x,y)=0,其中T<sub>tex</sub>表示相似度阈值;对任意像素(x,y),若S<sub>tex</sub>(x,y)=1且F<sub>0</sub>(x,y)=1,则前景掩膜F<sub>tex</sub>中对应像素值F<sub>tes</sub>(x,y)=1;否则F<sub>tes</sub>(x,y)=0;步骤4:基于色度、光谱方向、纹理三种特征的前景掩膜F<sub>chr</sub>,F<sub>phy</sub>,F<sub>tex</sub>的加权求和得到空域前景概率谱图像<img file="FDA0001099773090000023.GIF" wi="602" he="79" />其中ω<sub>c</sub>表示加权权值;步骤5:基于空域前景概率谱图像M<sub>prob</sub>计算时域前景概率普M<sub>prob</sub>(t):基于大小为w的时间滑动窗口,对窗口内的空域前景概率谱进行时域滤波,计算出当前时刻t的空域前景概率谱图像,即时域前景概率普<img file="FDA0001099773090000031.GIF" wi="620" he="127" />其中ω<sub>i</sub>表示(t‑i)时刻空域前景概率谱对应的权值:若i<w‑1,则ω<sub>i</sub>=1/2<sup>i+1</sup>;否则ω<sub>i</sub>=1/2<sup>w</sup>;步骤6:计算最终前景掩膜F<sub>final</sub>,完成对车辆阴影的消除处理:比较时域前景概率普M<sub>prob</sub>(t)与前景置信度阈值T<sub>c</sub>,计算出当前时刻t的空域前景概率谱图像的最终前景掩膜F<sub>final</sub>(t):若M<sub>prob</sub>(t)>T<sub>c</sub>,则F<sub>final</sub>(t)=1;否则F<sub>final</sub>(t)=0。
地址 611731 四川省成都市高新区(西区)西源大道2006号