发明名称 基于边缘和骨架信息的车辆阴影消除方法
摘要 本发明公布了一种基于边缘和骨架信息的车辆阴影消除方法,本发明方法对当前图像进行背景差计算并二值化得到前景二值图,通过对前景二值图进行边缘检测和膨胀得到车辆及阴影区域的边缘图像,通过对背景差图像进行边缘检测和膨胀得到车辆的骨架图像,然后,以车辆区域和阴影区域相对于背景区域的灰度比为判别特征,初步检测阴影区域,并用前景二值图减去初步阴影区域,最后利用骨架图像对车辆区域内的孔洞进行补充并减去边缘图像,得到最终消除阴影区域之后的图像。本发明综合利用了用以区分车辆区域和阴影区域的灰度特征、车辆与阴影区域的边缘信息及车辆骨架的信息,能够准确地检测和消除阴影区域。
申请公布号 CN101739560B 申请公布日期 2012.02.01
申请号 CN200910263082.8 申请日期 2009.12.16
申请人 东南大学 发明人 赵新勇;路小波;朱周;吴晶
分类号 G06K9/40(2006.01)I 主分类号 G06K9/40(2006.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 许方
主权项 1.一种基于边缘和骨架信息的车辆阴影消除方法,其特征在于包括如下步骤:①计算背景差图像并对其二值化交通场景的背景图像为BKG,需要进行车辆阴影消除的当前图像为P,计算当前图像P与背景图像BKG之间的绝对差,得到背景差图像BP,即:BP(x,y)=|P(x,y)-BKG(x,y)|,式中,P(x,y)表示图像P中坐标为(x,y)的像素的灰度,BKG(x,y)表示背景图像BKG中坐标为(x,y)的像素的灰度,BP(x,y)表示背景差图像BP中坐标为(x,y)的像素的灰度,x,y分别表示横坐标和纵坐标,下同;其次,对背景差图像BP进行二值化计算,得到前景二值图OM,<img file="DEST_PATH_FSB00000641006300011.GIF" wi="744" he="164" />式中,THR表示分割阈值;OM(x,y)表示前景二值图OM中坐标为(x,y)的像素的值,若为1表示该像素属于车辆区域,若为0表示该像素属于静止背景区域;②计算初始阴影图像用当前图像P除以背景图像BKG得到图像D,即:<img file="DEST_PATH_FSB00000641006300012.GIF" wi="1049" he="158" />式中,P(x,y)表示图像P中坐标为(x,y)的像素的灰度,BKG(x,y)表示背景图像BKG中坐标为(x,y)的像素的灰度,D(x,y)表示图像D中坐标为(x,y)的像素的值,根据如下判断规则计算初始阴影图像SH,<img file="DEST_PATH_FSB00000641006300013.GIF" wi="722" he="149" />式中,D(x,y)表示图像D中坐标为(x,y)的像素的值,T<sub>1</sub>和T<sub>2</sub>表示取值范围且<img file="DEST_PATH_FSB00000641006300014.GIF" wi="329" he="59" />Sh(x,y)表示初始阴影图像SH中坐标为(x,y)的像素的值,若为1表示该像素属于车辆的阴影区域,若为0则表示该像素不属于车辆的阴影区域;③计算车辆区域的边界利用Sobel边缘算子提取前景二值图OM中车辆区域的初始边界,得到图像<img file="DEST_PATH_FSB00000641006300015.GIF" wi="123" he="54" />即:OM<sub>1m</sub>=OM*G<sub>m</sub>,OM<sub>1n</sub>=OM*G<sub>n</sub>,<img file="DEST_PATH_FSB00000641006300021.GIF" wi="1271" he="137" />式中,G<sub>m</sub>=[-1 -2 -1;0 0 0;1 2 1],G<sub>n</sub>=[-1 01 ;-2 0 2;-1 0 1],S<sub>1</sub>为阈值,OM<sub>1m</sub>是前景二值图OM与G<sub>m</sub>进行卷积运算后得到的图像,OM<sub>1n</sub>是前景二值图OM与G<sub>n</sub>进行卷积运算后得到的图像,之后,利用膨胀算子Se=[1 1 1;1 1 1;1 1 1]对图像<img file="DEST_PATH_FSB00000641006300022.GIF" wi="101" he="51" />进行膨胀,得到车辆区域的边界图像OM1,④计算车辆区域的骨架利用Sobel边缘算子提取背景差图像BP中前景区域的骨架,得到图像<img file="DEST_PATH_FSB00000641006300023.GIF" wi="132" he="52" />即:OM<sub>2m</sub>=BP*G<sub>m</sub>,OM<sub>2n</sub>=BP*G<sub>n</sub>,<img file="DEST_PATH_FSB00000641006300024.GIF" wi="1326" he="150" />式中,G<sub>m</sub>=[-1 -2 -1;0 0 0;1 2 1],G<sub>n</sub>=[-1 0 1;-2 0 2;-1 0 1],S<sub>2</sub>为阈值,OM<sub>2m</sub>是图像BP与G<sub>m</sub>进行卷积运算后得到的图像,OM<sub>2n</sub>是图像BP与G<sub>n</sub>进行卷积运算后得到的图像,之后,利用膨胀算子Se=[1 1 1;1 1 1;1 1 1]对图像<img file="DEST_PATH_FSB00000641006300025.GIF" wi="111" he="52" />进行膨胀,得到车辆区域的骨架图像OM2,⑤计算车辆二值图像首先,用前景二值图OM减去初始阴影图像SH中重叠的区域,得到初始的车辆二值图像MV0,即:<img file="DEST_PATH_FSB00000641006300026.GIF" wi="893" he="144" />然后,用初始的车辆二值图像MV0加上车辆区域的骨架图像OM2,得到初步修正的车辆二值图像MV1,即:MV1(x,y)=MV0(x,y)or OM2(x,y), 最后,用初步修正的车辆图像MV1减去车辆区域的边界图像OM1,得到最终的车辆二值图像MV2,即:<img file="F2009102630828C00031.GIF" wi="950" he="150" />⑥消除阴影利用车辆二值图像MV2获得消除阴影后的车辆图像P,即:<img file="F2009102630828C00032.GIF" wi="846" he="137" />式中,P(x,y)表示图像P中坐标为(x,y)的像素的灰度,BKG(x,y)表示背景图像BKG中坐标为(x,y)的像素的灰度;MV2(x,y)表示车辆二值图像MV2中坐标为(x,y)的像素的值,若值为1表示该像素属于车辆区域,若值为0表示该像素不属于车辆区域。 
地址 210096 江苏省南京市四牌楼2号