CN107292898A - 一种基于hsv的车牌阴影检测和去除方法 - Google Patents
一种基于hsv的车牌阴影检测和去除方法 Download PDFInfo
- Publication number
- CN107292898A CN107292898A CN201710309363.7A CN201710309363A CN107292898A CN 107292898 A CN107292898 A CN 107292898A CN 201710309363 A CN201710309363 A CN 201710309363A CN 107292898 A CN107292898 A CN 107292898A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- image
- shadow
- coordinate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 20
- 238000005516 engineering process Methods 0.000 title abstract description 5
- 238000000034 method Methods 0.000 claims abstract description 32
- 238000007781 pre-processing Methods 0.000 claims description 9
- 238000003708 edge detection Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 abstract description 13
- 230000000295 complement effect Effects 0.000 abstract 1
- 238000012216 screening Methods 0.000 abstract 1
- 238000013459 approach Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000009191 jumping Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000013135 deep learning Methods 0.000 description 1
- 238000012850 discrimination method Methods 0.000 description 1
- 238000003709 image segmentation Methods 0.000 description 1
- 230000001502 supplementing effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/13—Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/155—Segmentation; Edge detection involving morphological operators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/60—Type of objects
- G06V20/62—Text, e.g. of license plates, overlay texts or captions on TV images
- G06V20/63—Scene text, e.g. street names
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/60—Type of objects
- G06V20/62—Text, e.g. of license plates, overlay texts or captions on TV images
- G06V20/625—License plates
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Multimedia (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于HSV的车牌阴影检测和去除方法。它是一种针对处理车牌阴影的阴影检测和去除方法,它首先通过处理车牌S、V通道确定车牌背景阴影区域,然后通过阴影边缘点集提取、筛选、连接和补充操作,得到完整阴影边缘,最后根据阴影区域和非阴影区域的平均灰度比,对H、S、V通道进行阴影去除,同时还利用基于筛选的中值滤波方法对阴影边缘区域进行修正,保证图像过度自然。
Description
技术领域
本发明涉及智能交通技术领域,具体涉及一种基于HSV的车牌阴影检测和去除方法。
背景技术
车牌识别系统是智能交通系统的重要组成部分,车牌预处理是车牌识别系统的重要环节,预处理的质量高低直接影响最终识别结果。而当下,车牌存在阴影是拉低预处理质量最主要的因素,因此去除车牌阴影,提高车牌的可识别率是车牌识别系统亟待解决的问题。
由于图像阴影严重影响目标检测、图像切割以及图像识别等算法的处理效果,因此阴影检测和去除是近几年较热门的一个研究点。阴影检测的方法大致可分为以下几种:基于分类的方法、基于边缘检测的方法和基于自学习特征的方法。基于分类的方法根据不同特征设计分类方法或分类器,以实现区域级或像素级的阴影分类,主要用到的特征包括亮度、纹理、颜色以及通道特性;基于边缘检测的方法利用阴影边缘区域的颜色、梯度、亮度变化特征进行阴影边缘检测,同时通过阴影区域判别方法确定阴影区域;基于自学习特征的方法通过引入深度学习算法自动学习阴影特征用于阴影检测,这种方法无需任何前提假设且鲁棒性较好。目前已有的阴影去除方法步骤大多类似,首先计算阴影区域和非阴影区域的亮度比,然后以重新点亮阴影区域的方法实现阴影去除,计算亮度比的方式可分为基于多通道的方法和基于亮度模型的方法。
目前大多阴影检测和去除的算法主要研究复杂场景下的阴影问题,算法的鲁棒性较好,但复杂性过高,不适用于处理颜色和结构较为简单的车牌图像,因此本文通过分析车牌的结构和通道特性,提出了一种基于HSV的阴影检测和去除方法。
发明内容
针对现有技术中存在的上述问题,本发明的目的在于提供一种基于HSV的车牌阴影检测和去除方法。它可以处理蓝底白字和黄底黑字且由于车辆自遮挡造成的阴影车牌,在该类阴影车牌中,阴影区域和非阴影区域上下分布,且阴影下边界曲线较为平滑。
所述的一种基于HSV的车牌阴影检测和去除方法,其特征在于包括如下步骤:
步骤1:对获取的车牌阴影图像D进行图像预处理,得到二值化图像Db,其中,车牌图像的宽度为W,高度为H,单位为像素,图像预处理过程具体如下:
步骤1.1:对图像D的S通道进行OTSU二值化和膨胀操作,得到二值化图D2;
步骤1.2:对图像D的V通道进行中值滤波操作,得到灰度图D3;
步骤1.3:根据公式(1)和(2)获得字符区域的灰度值集合E1和背景区域的灰度值集合E2,并利用OTSU计算E2的最佳阈值,设为thE2,然后根据公式(3)得到图像D3的二值化图像D4;
E1={PD3(x,y)|PD2(x,y)=255,0≤x<W,0≤y<H} (1)
E2={PD3(x,y)|PD2(x,y)=0,0≤x<W,0≤y<H} (2)
其中,表示图像D2中坐标(x,y)处的像素的灰度值,表示图像D3中坐标(x,y)处的像素的灰度值,表示图像D4中坐标(x,y)处的像素的灰度值;
步骤1.4:去除图像D4中轮廓面积小于A的小团块,得到用于阴影边缘检测的二值化图像Db,其中A表示预先设置的轮廓面积最小阈值;
步骤2:对图像Db进行阴影边缘提取,确定车牌阴影边缘点坐标集B,具体为:
步骤2.1:根据图像Db构建集合B1={(xi,yi)|i=0,1,…,Nb1-1},其中,Nb1表示集合B1中元素的个数,(xi,yi)表示像素的坐标且须满足式(4);
且
其中,表示图像Db中坐标(xi,yi)处的像素的灰度值;
步骤2.2:从坐标集B1中剃除那些满足255∈N8D2((xi,yi))的(xi,yi),其中N8D2((xi,yi))表示图像D2中坐标(xi,yi)处的像素的8邻域的灰度值集合;
步骤2.3:对B1中的元素进行聚类,得到满足式(5)~(7)的Nh个坐标子集的集合H={Hj|j=0,1,…,Nh-1},其中第j个子集为 表示子集Hj中元素的个数;
且当p≠q时Hp∩Hq=φ (5)
对于任意存在使
对于任意不存在使且i≠j(7)其中,表示坐标的8邻域的坐标集合;
步骤2.4:将集合H中的子集按子集第一个元素的横坐标从小到大进行排序,然后从H中找到元素个数最多的子集并将其标记为正确的,最后根据如下步骤从坐标集B1中剔除那些位置跳变的坐标点:
步骤2.4.1:从H中找到被标记为正确的子集,设为Hk,若k>0且Hk-1未标记过,令Hc=Hk-1,转步骤2.4.2;否则,若k<Nh-1且Hk+1为未标记过,令Hc=Hk+1,转步骤2.4.2;直至H中不存在满足上述两个条件之一的Hk转步骤2.5;
步骤2.4.2:找到Hc与Hk两个子集中横坐标距离最近的坐标,分别设为和若满足且则判定子集Hc中都为跳变坐标点,将包含在子集Hc中的元素从集合B1中剔除,同时从H集合中剔除Hc,转步骤2.4.1;否则将Hc标记为正确的,转步骤2.4.1;
步骤2.5:顺序遍历坐标集B1中下标满足式(8)的相邻坐标点对每组这样的相邻三个坐标点用B样条二次曲线法进行连接,得到一条连接集合B1中所有坐标点的连续曲线,并将这条曲线经过的所有坐标重新保存为坐标集B2={(xi,yi)|i=0,1,…,Nb2-1},其中Nb2表示集合B2中元素的个数;
步骤2.6:若B2中或则不再执行下面步骤,阴影检测失败;否则,根据式(9)和(10)依次确定0≤ai<x0以及的残缺边缘坐标点(ai,bi)和(ak,bk),并将这些坐标点依次加入到新的坐标集合B3中;
步骤2.7:合并集合B2和B3得到边缘点坐标集B={(xi,yi)|i=0,1,…,Nb-1},其中Nb表示集合B中元素的个数,并将集合B按yi从小到大进行排序;
步骤3:对车牌图像D进行分通道阴影去除,得到去除阴影后的车牌彩图DR,具体为:
步骤3.1:记车牌图像D的V通道为图像Dv,对Dv的阴影区域进行灰度值调整,首先计算图像Dv中阴影区域和非阴影区域的平均灰度比值,然后根据平均灰度比值对阴影区域进行灰度值调整,并对阴影边缘进行灰度值修正,具体为:
步骤3.1.1:得到Dv中所有坐标(x,y)满足式(11)的像素的灰度值集合M1={pi|i=0,1,…,Nm1-1},其中Nm1表示集合M1中元素的个数;将集合M1按照pi从小到大进行排序,并根据式(12)计算M1对应的平均灰度值pvb;
存在(xi,yi)∈B,使x=xi且且PD2(x,y)=0(11)
步骤3.1.2:对Dv中所有坐标(x,y)满足式(13)的像素进行聚类,得到满足式(14)~(15)的Nu个子集其中,表示像素在Dv中的坐标,表示像素的灰度值,表示集合Ui中元素的个数;将每个Ui中的元素按照从小到大进行排序,并根据式(15)得到各个集合Ui对应的平均灰度值
存在(xi,yi)∈B,使x=xi且y<yi且PD2(x,y)=0 (13)
存在(xz,yz),(xk,yk)∈B,
使
不存在(xz,yz),(xk,yk)∈B,
使
其中,γ表示预先设置的像素个数的最小阈值;
步骤3.1.3:顺序遍历i=0,1,…,Nu-1,若则找到某个并令其中,不等于-1且离距离最近,此处距离表示为|t-i|;
步骤3.1.4:依次对Dv中坐标(x,y)满足式(17)的像素的灰度值根据式(18)进行调整;
存在(xi,yi)∈B,使x=xi且y<yi (17)
(xk,yk)∈B且xk=x (18)
步骤3.1.5:构建集合O={(ni,mi)|i=0,1,…,No-1},其中,No表示满足式(19)的坐标数量,(ni,mi)表示图像Dv中某个像素的坐标且满足式(19);将集合O中的元素按照di从大到小进行排序,其中di=|mi-yk|,yk满足(xk,yk)∈B且xk=ni;
存在(xk,yk)∈B,使xk=ni且
从D中得到原始V通道图像,记为Dv2,依次遍历坐标集O中元素(ni,mi),并对Dv中坐标(ni,mi)处的像素的灰度值按照步骤(a)–(c)进行调整:
(a)遍历图像Dv2和Dv中的所有像素并构建集合L={(xj,yj,qj,pj)|j=0,1,…,NL-1},其中,qj为Dv2中坐标(xj,yj)处的像素的灰度值,pj为Dv中坐标(xj,yj)处的像素的灰度值,NL表示集合L中元素的个数,(xj,yj)须满足式(20);
其中,分别表示图像Ds中坐标(xj,yj)和(mj,nj)处像素的灰度值,thS表示图像Ds的OTSU最佳阈值,Ds表示车牌图像D的S通道图像;
(b)将集合L按照从小到大进行排序,其中表示Dv2中坐标(mj,nj)处像素的灰度值;保留L前一半的元素为新的集合并将集合L2按照pj从小到大进行排序;
(c)从L2中找到第个元素并将Dv中(ni,mi)处像素的灰度值设置为同时从集合O中剔除(ni,mi);
步骤3.2:对图像Ds进行灰度值调整,首先计算Ds中阴影区域和非阴影区域的平均灰度比值,然后根据平均灰度比值对Ds中的阴影区进行灰度值调整,具体为:
步骤3.2.1:得到Ds中所有坐标(x,y)满足式(11)的像素的灰度值集合M2={pi|i=0,1,…,Nm2-1},其中Nm2表示集合M2中元素的个数;将集合M2按照pi从小到大进行排序,并根据式(21)计算M2对应的平均灰度值psb;
步骤3.2.2:得到Ds中所有坐标(x,y)满足式(22)的像素的灰度值集合M3={pi|i=0,1,…,Nm3-1},其中Nm3表示集合M3中元素的个数;将集合M3按照pi从小到大进行排序,并根据式(23)计算M3对应的平均灰度值psa;
存在(xi,yi)∈B,使x=xi且且PD2(x,y)=0 (22)
步骤3.2.3:依次对Ds中坐标(x,y)满足式(17)的像素的灰度值根据式(24)进行调整;
步骤3.3:记车牌图像D的H通道图像为Dh,对Dh进行灰度值调整,首先计算Dh中阴影区域和非阴影区域的平均灰度比值,然后根据平均灰度比值对Dh中的阴影区进行灰度值调整,具体为:
步骤3.3.1:得到Dh中所有坐标(x,y)满足式(11)的像素的灰度值集合M4={pi|i=0,1,…,Nm4-1},其中Nm4表示集合M4中元素的个数;将集合M4按照pi从小到大进行排序,并根据式(25)计算M4对应的平均灰度值phb;
步骤3.3.2:得到Dh中所有坐标(x,y)满足式(22)的像素的灰度值集合M5={pi|i=0,1,…,Nm5-1},其中Nm5表示集合M5中元素的个数;将集合M5按照pi从小到大进行排序,并根据式(26)计算M5对应的平均灰度值pha;
步骤3.3.3:依次对Dh中坐标(x,y)满足式(17)的像素的灰度值根据式(27)进行调整;
步骤3.4:对调整后的图像Dh,Ds,Dv进行整合,得到去除阴影后的车牌彩图DR。
通过采用上述技术,与现有技术相比,本发明的有益效果如下:
本发明通过采用上述限定技术得到的一种针对处理车牌阴影的阴影检测和去除方法,它利用车牌结构和通道特性实现了车牌背景阴影的精准定位,同时基于HSV通道的阴影区域方法对阴影车牌图像的修复效果较好,修复得到的车牌彩图显得十分自然。
附图说明
图1为原阴影车牌彩图的灰度图;
图2为车牌S通道图像;
图3为根据步骤1.1处理得到的图像D2;
图4为车牌V通道图像;
图5为根据步骤1.2处理得到的灰度图D3;
图6为车牌V通道背景区域图像;
图7为车牌V通道字符区域图像;
图8为根据步骤1.3处理得到的二值化图像D4;
图9为根据步骤2.1处理得到的坐标集B1;
图10为根据步骤2.2-2.4处理得到的坐标集B1;
图11为根据步骤2.5处理得到的坐标集B2;
图12为根据步骤2.6-2.7处理得到的坐标集B;
图13为根据步骤3.1.1-3.1.4处理得到的图像Dv;
图14为根据步骤3.1.5处理得到的图像Dv
图15为步骤3.2处理得到的图像Ds;
图16为车牌H通道图像;
图17为根据步骤3.3处理得到的图像Dh;
图18为根据步骤3.4处理得到的图像DR。
具体实施方式
下面结合具体实例,对本发明实施的过程进行详细的说明。
如图1-18所示,本发明的一种基于HSV的车牌阴影检测和去除方法,包括如下步骤:
步骤1:对获取的车牌阴影图像D进行图像预处理,得到二值化图像Db,其中,车牌图像的宽度为W,高度为H,单位为像素,图像预处理过程具体如下:
本实例中,车牌原阴影图像D经灰度处理后如图1所示,图像宽度204像素,高度54像素;
步骤1.1:对图像D的S通道进行OTSU二值化和膨胀操作,得到二值化图D2;
阴影车牌的S通道如图2所示,根据步骤1.1对S通道进行二值化和膨胀操作,得到的二值化图像D2如图3所示;
步骤1.2:对图像D的V通道进行中值滤波操作,得到灰度图D3;
阴影车牌的V通道如图4所示,根据步骤1.2对V通道进行中值滤波,得到的灰度图D3如图5所示;
步骤1.3:根据公式(1)和(2)获得字符区域的灰度值集合E1和背景区域的灰度值集合E2,并利用OTSU计算E2的最佳阈值,设为thE2,然后根据公式(3)得到图像D3的二值化图像D4;
E1={PD3(x,y)|PD2(x,y)=255,0≤x<W,0≤y<H} (1)
E2={PD3(x,y)|PD2(x,y)=0,0≤x<W,0≤y<H} (2)
其中,表示图像D2中坐标(x,y)处的像素的灰度值,表示图像D3中坐标(x,y)处的像素的灰度值,表示图像D4中坐标(x,y)处的像素的灰度值;
灰度图D3中的背景区域和字符区域图像如图6和7所示,根据步骤1.3所述二值化方法得到的二值化图像D4如图8所示;
步骤1.4:去除图像D4中轮廓面积小于A的小团块,得到用于阴影检测的二值化图像Db,其中A表示预先设置的轮廓面积最小阈值;
在本实例中,轮廓面积最小阈值A取20,由于图像D4中不存在小团块,因此根据步骤1.4处理图像D4后图像没有发生改变,图像Db如图8所示;
步骤2:对图像Db进行阴影边缘提取,确定车牌阴影边缘点坐标集B,具体为:
步骤2.1:根据图像Db构建集合B1={(xi,yi)|i=0,1,…,Nb1-1},其中,Nb1表示集合B1中元素的个数,(xi,yi)表示像素的坐标且须满足式(4)。
且
其中,表示图像Db中坐标(xi,yi)处的像素的灰度值;
根据步骤2.1得到坐标集B1如图9所示,其中白点表示集合B1中的坐标点。
步骤2.2:从坐标集B1中剃除那些满足255∈N8D2((xi,yi))的(xi,yi),其中N8D2((xi,yi))表示图像D2中坐标(xi,yi)处的像素的8邻域的灰度值集合;
步骤2.3:对B1中的元素进行聚类,得到满足式(5)~(7)的Nh个坐标子集的集合H={Hj|j=0,1,…,Nh-1},其中第j个子集为 表示子集Hj中元素的个数;
且当p≠q时Hp∩Hq=φ (5)
对于任意存在使
对于任意不存在使且i≠j(7)其中,表示坐标的8邻域的坐标集合;
步骤2.4:将集合H中的子集按子集第一个元素的横坐标从小到大进行排序,然后从H中找到元素个数最多的子集并将其标记为正确的,最后根据如下步骤从坐标集B1中剔除那些位置跳变的坐标点:
步骤2.4.1:从H中找到被标记为正确的子集,设为Hk,若k>0且Hk-1未标记过,令Hc=Hk-1,转步骤2.4.2;否则,若k<Nh-1且Hk+1为未标记过,令Hc=Hk+1,转步骤2.4.2;直至H中不存在满足上述两个条件之一的Hk转步骤2.5。
步骤2.4.2:找到Hc与Hk两个子集中横坐标距离最近的坐标,分别设为和若满足且则判定子集Hc中都为跳变坐标点,将包含在子集Hc中的元素从集合B1中剔除,同时从H集合中剔除Hc,转步骤2.4.1;否则将Hc标记为正确的,转步骤2.4.1。
根据步骤2.2-2.4剔除坐标集B1中的字符边缘点坐标集和跳变曲线点坐标集,得到的坐标集B1如图10所示,其中白点表示集合B1中的坐标点;
步骤2.5:顺序遍历坐标集B1中下标满足式(8)的相邻坐标点对每组这样的相邻三个坐标点用B样条二次曲线法进行连接,得到一条连接集合B1中所有坐标点的连续曲线,并将这条曲线经过的所有坐标重新保存为坐标集B2={(xi,yi)|i=0,1,…,Nb2-1},其中Nb2表示集合B2中元素的个数;
根据步骤2.5进行曲线拟合,得到边缘点坐标集B2如图11所示,其中白点表示集合B1中的坐标点;
步骤2.6:若B2中或则不再执行下面步骤,阴影检测失败;否则,根据式(9)和(10)依次确定0≤ai<x0以及的残缺边缘坐标点(ai,bi)和(ak,bk),并将这些坐标点依次加入到新的坐标集合B3中。
步骤2.7:合并集合B2和B3得到边缘点坐标集B={(xi,yi)|i=0,1,…,Nb-1},其中Nb表示集合B中元素的个数,并将集合B按yi从小到大进行排序;
根据步骤2.6-2.7对残缺边缘点坐标进行补充,得到边缘点坐标集B如图12所示,其中白点表示集合B1中的坐标点;
步骤3:对阴影车牌D进行分通道阴影去除,得到去除阴影后的车牌彩图DR,具体为:
步骤3.1:记车牌图像D的V通道为图像Dv,对Dv的阴影区域进行灰度值调整,首先计算图像Dv中阴影区域和非阴影区域的平均灰度比值,然后根据平均灰度比值对阴影区域进行灰度值调整,并对阴影边缘进行灰度值修正,具体为:
步骤3.1.1:得到Dv中所有坐标(x,y)满足式(11)的像素的灰度值集合M1={pi|i=0,1,…,Nm1-1},其中Nm1表示集合M1中元素的个数;将集合M1按照pi从小到大进行排序,并根据式(12)计算M1对应的平均灰度值pvb;
存在(xi,yi)∈B,使x=xi且且PD2(x,y)=0 (11)
步骤3.1.2:对Dv中所有坐标(x,y)满足式(13)的像素进行聚类,得到满足式(14)~(15)的Nu个子集其中,表示像素在Dv中的坐标,表示像素的灰度值,表示集合Ui中元素的个数;将每个Ui中的元素按照从小到大进行排序,并根据式(15)得到各个集合Ui对应的平均灰度值
存在(xi,yi)∈B,使x=xi且y<yi且PD2(x,y)=0 (13)
存在(xz,yz),(xk,yk)∈B,
使
不存在(xz,yz),(xk,yk)∈B,
使
其中,γ表示预先设置的像素个数的最小阈值;
步骤3.1.3:顺序遍历i=0,1,…,Nu-1,若则找到某个并令其中,不等于-1且离距离最近,此处距离表示为|t-i|;
步骤3.1.4:依次对Dv中坐标(x,y)满足式(17)的像素的灰度值根据式(18)进行调整;
存在(xi,yi)∈B,使x=xi且y<yi (17)
(xk,yk)∈B且xk=x (18)
在本实例中,预先设置的像素个数的最小阈值γ为0.02×W×H,阴影车牌原V通道图像如图4所示,根据步骤3.1.1-3.1.4对V通道中的阴影区域进行像素调整,得到的图像Dv如图13所示;
步骤3.1.5:构建集合O={(ni,mi)|i=0,1,…,No-1},其中,No表示满足式(19)的坐标数量,(ni,mi)表示图像Dv中某个像素的坐标且满足式(19);将集合O中的元素按照di从大到小进行排序,其中di=|mi-yk|,yk满足(xk,yk)∈B且xk=ni;
存在(xk,yk)∈B,使xk=ni且
从D中得到原始V通道图像,记为Dv2,依次遍历坐标集O中元素(ni,mi),并对Dv中坐标(ni,mi)处的像素的灰度值按照步骤(a)–(c)进行调整:
(a)遍历图像Dv2和Dv中的所有像素并构建集合L={(xj,yj,qj,pj)|j=0,1,…,NL-1},其中,qj为Dv2中坐标(xj,yj)处的像素的灰度值,pj为Dv中坐标(xj,yj)处的像素的灰度值,NL表示集合L中元素的个数,(xj,yj)须满足式(20)。
其中,分别表示图像Ds中坐标(xj,yj)和(mj,nj)处像素的灰度值,thS表示图像Ds的OTSU最佳阈值,Ds表示车牌图像D的S通道图像。
(b)将集合L按照从小到大进行排序,其中表示Dv2中坐标(mj,nj)处像素的灰度值;保留L前一半的元素为新的集合并将集合L2按照pj从小到大进行排序。
(c)从L2中找到第个元素并将Dv中(ni,mi)处像素的灰度值设置为同时从集合O中剔除(ni,mi)。
根据步骤3.1.5对Dv中的阴影边缘区域的像素的灰度值进行修正,得到的图像Dv如图14所示;
步骤3.2:对图像Ds进行灰度值调整,首先计算Ds中阴影区域和非阴影区域的平均灰度比值,然后根据平均灰度比值对Ds中的阴影区进行灰度值调整,具体为:
步骤3.2.1:得到Ds中所有坐标(x,y)满足式(11)的像素的灰度值集合M2={pi|i=0,1,…,Nm2-1},其中Nm2表示集合M2中元素的个数;将集合M2按照pi从小到大进行排序,并根据式(21)计算M2对应的平均灰度值psb;
步骤3.2.2:得到Ds中所有坐标(x,y)满足式(22)的像素的灰度值集合M3={pi|i=0,1,…,Nm3-1},其中Nm3表示集合M3中元素的个数;将集合M3按照pi从小到大进行排序,并根据式(23)计算M3对应的平均灰度值psa;
存在(xi,yi)∈B,使x=xi且且PD2(x,y)=0 (22)
步骤3.2.3:依次对Ds中坐标(x,y)满足式(17)的像素的灰度值根据式(24)进行调整;
阴影车牌原S通道图像如图2所示,根据步骤3.2对S通道中阴影区域的像素进行灰度值调整,得到的图像Ds如图15所示;
步骤3.3:记车牌图像D的H通道图像为Dh,对Dh进行灰度值调整,首先计算Dh中阴影区域和非阴影区域的平均灰度比值,然后根据平均灰度比值对Dh中的阴影区进行灰度值调整,具体为:
步骤3.3.1:得到Dh中所有坐标(x,y)满足式(11)的像素的灰度值集合M4={pi|i=0,1,…,Nm4-1},其中Nm4表示集合M4中元素的个数;将集合M4按照pi从小到大进行排序,并根据式(25)计算M4对应的平均灰度值phb;
步骤3.3.2:得到Dh中所有坐标(x,y)满足式(22)的像素的灰度值集合M5={pi|i=0,1,…,Nm5-1},其中Nm5表示集合M5中元素的个数;将集合M5按照pi从小到大进行排序,并根据式(26)计算M5对应的平均灰度值pha;
步骤3.3.3:依次对Dh中坐标(x,y)满足式(17)的像素的灰度值根据式(27)进行调整;
阴影车牌原H通道图像如图16所示,根据步骤3.3对H通道中阴影区域的像素进行灰度值调整,得到的图像Dh如图17所示;
步骤3.4:对调整后的图像Dh,Ds,Dv进行整合,得到去除阴影后的车牌彩图DR;
根据步骤3.4得到去除阴影后的车牌彩图DR经灰度处理后的图如图18所示;
本说明书具体实施方法中所列举的处理对象,仅用于说明本发明的实现过程,本发明所能处理的对象情况并非仅限于所举实例。
Claims (2)
1.一种基于HSV的车牌阴影检测和去除方法,其特征在于包括如下步骤:
步骤1:对获取的车牌阴影图像D进行图像预处理,得到二值化图像Db,其中,车牌图像的宽度为W,高度为H,单位为像素;
步骤2:对图像Db进行阴影边缘提取,确定车牌阴影边缘点坐标集B,具体为:
步骤2.1:根据图像Db构建集合B1={(xi,yi)|i=0,1,…,Nb1-1},其中,Nb1表示集合B1中元素的个数,(xi,yi)表示像素的坐标且须满足式(4);
且
其中,表示图像Db中坐标(xi,yi)处的像素的灰度值;
步骤2.2:从坐标集B1中剃除那些满足255∈N8D2((xi,yi))的(xi,yi),其中N8D2((xi,yi))表示图像D2中坐标(xi,yi)处的像素的8邻域的灰度值集合;
步骤2.3:对B1中的元素进行聚类,得到满足式(5)~(7)的Nh个坐标子集的集合H={Hj|j=0,1,…,Nh-1},其中第j个子集为 表示子集Hj中元素的个数;
且当p≠q时Hp∩Hq=φ (5)
对于任意存在使
对于任意不存在使且i≠j(7)其中,表示坐标的8邻域的坐标集合;
步骤2.4:将集合H中的子集按子集第一个元素的横坐标从小到大进行排序,然后从H中找到元素个数最多的子集并将其标记为正确的,最后根据如下步骤从坐标集B1中剔除那些位置跳变的坐标点:
步骤2.4.1:从H中找到被标记为正确的子集,设为Hk,若k>0且Hk-1未标记过,令Hc=Hk-1,转步骤2.4.2;否则,若k<Nh-1且Hk+1为未标记过,令Hc=Hk+1,转步骤2.4.2;直至H中不存在满足上述两个条件之一的Hk转步骤2.5;
步骤2.4.2:找到Hc与Hk两个子集中横坐标距离最近的坐标,分别设为和若满足且则判定子集Hc中都为跳变坐标点,将包含在子集Hc中的元素从集合B1中剔除,同时从H集合中剔除Hc,转步骤2.4.1;否则将Hc标记为正确的,转步骤2.4.1;
步骤2.5:顺序遍历坐标集B1中下标满足式(8)的相邻坐标点对每组这样的相邻三个坐标点用B样条二次曲线法进行连接,得到一条连接集合B1中所有坐标点的连续曲线,并将这条曲线经过的所有坐标重新保存为坐标集B2={(xi,yi)|i=0,1,…,Nb2-1},其中Nb2表示集合B2中元素的个数;
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>q</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mi>max</mi>
<mrow>
<mo>(</mo>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mi>i</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>q</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mi>min</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>b</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>b</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<mn>1</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
1
步骤2.6:若B2中或则不再执行下面步骤,阴影检测失败;否则,根据式(9)和(10)依次确定0≤ai<x0以及的残缺边缘坐标点(ai,bi)和(ak,bk),并将这些坐标点依次加入到新的坐标集合B3中;
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mi>i</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>y</mi>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mn>0</mn>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mi>min</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>max</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>b</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<mn>0</mn>
</mrow>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>H</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>0</mn>
<mo>,</mo>
<mn>1</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mn>0</mn>
</msub>
<mo>-</mo>
<mn>1</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<mi>k</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>b</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<msub>
<mi>y</mi>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>b</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<mn>1</mn>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mrow>
<mi>k</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>b</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>-</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>y</mi>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>b</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<mn>1</mn>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mrow>
<mi>W</mi>
<mo>-</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>b</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>b</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>b</mi>
<mi>k</mi>
</msub>
<mo>=</mo>
<mi>min</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>max</mi>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>b</mi>
<mi>k</mi>
</msub>
<mo>,</mo>
<mn>0</mn>
</mrow>
<mo>)</mo>
</mrow>
<mo>,</mo>
<mi>H</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>k</mi>
<mo>=</mo>
<msub>
<mi>x</mi>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>b</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>b</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>W</mi>
<mo>-</mo>
<mn>1</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤2.7:合并集合B2和B3得到边缘点坐标集B={(xi,yi)|i=0,1,…,Nb-1},其中Nb表示集合B中元素的个数,并将集合B按yi从小到大进行排序;
步骤3:对车牌图像D进行分通道阴影去除,得到去除阴影后的车牌彩图DR,具体为:
步骤3.1:记车牌图像D的V通道为图像Dv,对Dv的阴影区域进行灰度值调整,首先计算图像Dv中阴影区域和非阴影区域的平均灰度比值,然后根据平均灰度比值对阴影区域进行灰度值调整,并对阴影边缘进行灰度值修正,具体为:
步骤3.1.1:得到Dv中所有坐标(x,y)满足式(11)的像素的灰度值集合M1={pi|i=0,1,…,Nm1-1},其中Nm1表示集合M1中元素的个数;将集合M1按照pi从小到大进行排序,并根据式(12)计算M1对应的平均灰度值pvb;
存在(xi,yi)∈B,使x=xi且且PD2(x,y)=0 (11)
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>v</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>m</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>/</mo>
<mn>4</mn>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>m</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>m</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>/</mo>
<mn>4</mn>
</mrow>
</munderover>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>m</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>/</mo>
<mn>2</mn>
</mrow>
</mfrac>
<mo>;</mo>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
<mo>&Element;</mo>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤3.1.2:对Dv中所有坐标(x,y)满足式(13)的像素进行聚类,得到满足式(14)~(15)的Nu个子集其中,表示像素在Dv中的坐标,表示像素的灰度值,表示集合Ui中元素的个数;将每个Ui中的元素按照从小到大进行排序,并根据式(15)得到各个集合Ui对应的平均灰度值
存在(xi,yi)∈B,使x=xi且y<yi且PD2(x,y)=0 (13)
存在(xz,yz),(xk,yk)∈B,
使
不存在(xz,yz),(xk,yk)∈B,
使
<mrow>
<msubsup>
<mi>p</mi>
<mrow>
<mi>v</mi>
<mi>a</mi>
</mrow>
<mi>i</mi>
</msubsup>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>j</mi>
<mo>=</mo>
<msubsup>
<mi>N</mi>
<mi>u</mi>
<mi>i</mi>
</msubsup>
<mo>/</mo>
<mn>6</mn>
</mrow>
<mrow>
<msubsup>
<mi>N</mi>
<mi>u</mi>
<mi>i</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>N</mi>
<mi>u</mi>
<mi>i</mi>
</msubsup>
<mo>/</mo>
<mn>6</mn>
</mrow>
</munderover>
<msubsup>
<mi>p</mi>
<mi>j</mi>
<mi>i</mi>
</msubsup>
</mrow>
<mrow>
<msubsup>
<mi>N</mi>
<mi>u</mi>
<mi>i</mi>
</msubsup>
<mo>/</mo>
<mn>3</mn>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mfrac>
<msubsup>
<mi>N</mi>
<mi>u</mi>
<mi>i</mi>
</msubsup>
<mn>3</mn>
</mfrac>
<mo>></mo>
<mi>&gamma;</mi>
<mo>,</mo>
<msubsup>
<mi>u</mi>
<mi>j</mi>
<mi>i</mi>
</msubsup>
<mo>&Element;</mo>
<msub>
<mi>U</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
<mo>;</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mfrac>
<msubsup>
<mi>N</mi>
<mi>u</mi>
<mi>i</mi>
</msubsup>
<mn>3</mn>
</mfrac>
<mo>&le;</mo>
<mi>&gamma;</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>16</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,γ表示预先设置的像素个数的最小阈值;
步骤3.1.3:顺序遍历若则找到某个并令其中,不等于-1且离距离最近,此处距离表示为|t-i|;
步骤3.1.4:依次对Dv中坐标(x,y)满足式(17)的像素的灰度值根据式(18)进行调整;
存在(xi,yi)∈B,使x=xi且y<yi (17)
(xk,yk)∈B且xk=x (18)
步骤3.1.5:构建集合O={(ni,mi)|i=0,1,…,No-1},其中,No表示满足式(19)的坐标数量,(ni,mi)表示图像Dv中某个像素的坐标且满足式(19);将集合O中的元素按照di从大到小进行排序,其中di=|mi-yk|,yk满足(xk,yk)∈B且xk=ni;
存在(xk,yk)∈B,使xk=ni且
从D中得到原始V通道图像,记为Dv2,依次遍历坐标集O中元素(ni,mi),并对Dv中坐标(ni,mi)处的像素的灰度值按照步骤(a)–(c)进行调整:
(a)遍历图像Dv2和Dv中的所有像素并构建集合L={(xj,yj,qj,pj)|j=0,1,…,NL-1},其中,qj为Dv2中坐标(xj,yj)处的像素的灰度值,pj为Dv中坐标(xj,yj)处的像素的灰度值,NL表示集合L中元素的个数,(xj,yj)须满足式(20);
其中,分别表示图像Ds中坐标(xj,yj)和(mj,nj)处像素的灰度值,thS表示图像Ds的OTSU最佳阈值,Ds表示车牌图像D的S通道图像;
(b)将集合L按照从小到大进行排序,其中表示中坐标(mj,nj)处像素的灰度值;保留L前一半的元素为新的集合并将集合L2按照pj从小到大进行排序;
(c)从L2中找到第个元素并将Dv中(ni,mi)处像素的灰度值设置为同时从集合O中剔除(ni,mi);
步骤3.2:对图像Ds进行灰度值调整,首先计算Ds中阴影区域和非阴影区域的平均灰度比值,然后根据平均灰度比值对Ds中的阴影区进行灰度值调整,具体为:
步骤3.2.1:得到Ds中所有坐标(x,y)满足式(11)的像素的灰度值集合M2={pi|i=0,1,…,Nm2-1},其中Nm2表示集合M2中元素的个数;将集合M2按照pi从小到大进行排序,并根据式(21)计算M2对应的平均灰度值psb;
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>s</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>m</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>/</mo>
<mn>4</mn>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>m</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>m</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>/</mo>
<mn>4</mn>
</mrow>
</munderover>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>m</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>/</mo>
<mn>2</mn>
</mrow>
</mfrac>
<mo>;</mo>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
<mo>&Element;</mo>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>21</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤3.2.2:得到Ds中所有坐标(x,y)满足式(22)的像素的灰度值集合M3={pi|i=0,1,…,Nm3-1},其中Nm3表示集合M3中元素的个数;将集合M3按照pi从小到大进行排序,并根据式(23)计算M3对应的平均灰度值psa;
存在(xi,yi)∈B,使x=xi且且PD2(x,y)=0(22)
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>m</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>/</mo>
<mn>4</mn>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>m</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>m</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>/</mo>
<mn>4</mn>
</mrow>
</munderover>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>m</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>/</mo>
<mn>2</mn>
</mrow>
</mfrac>
<mo>;</mo>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
<mo>&Element;</mo>
<msub>
<mi>M</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>23</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤3.2.3:依次对Ds中坐标(x,y)满足式(17)的像素的灰度值根据式(24)进行调整;
<mrow>
<msub>
<mi>P</mi>
<msub>
<mi>D</mi>
<mi>S</mi>
</msub>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>P</mi>
<msub>
<mi>D</mi>
<mi>S</mi>
</msub>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>p</mi>
<mrow>
<mi>s</mi>
<mi>b</mi>
</mrow>
</msub>
<msub>
<mi>p</mi>
<mrow>
<mi>s</mi>
<mi>a</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>24</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤3.3:记车牌图像D的H通道图像为Dh,对Dh进行灰度值调整,首先计算Dh中阴影区域和非阴影区域的平均灰度比值,然后根据平均灰度比值对Dh中的阴影区进行灰度值调整,具体为:
步骤3.3.1:得到Dh中所有坐标(x,y)满足式(11)的像素的灰度值集合M4={pi|i=0,1,…,Nm4-1},其中Nm4表示集合M4中元素的个数;将集合M4按照pi从小到大进行排序,并根据式(25)计算M4对应的平均灰度值phb;
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>h</mi>
<mi>b</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>m</mi>
<mn>4</mn>
</mrow>
</msub>
<mo>/</mo>
<mn>4</mn>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>m</mi>
<mn>4</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>m</mi>
<mn>4</mn>
</mrow>
</msub>
<mo>/</mo>
<mn>4</mn>
</mrow>
</munderover>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>m</mi>
<mn>4</mn>
</mrow>
</msub>
<mo>/</mo>
<mn>2</mn>
</mrow>
</mfrac>
<mo>;</mo>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
<mo>&Element;</mo>
<msub>
<mi>M</mi>
<mn>4</mn>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>25</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤3.3.2:得到Dh中所有坐标(x,y)满足式(22)的像素的灰度值集合M5={pi|i=0,1,…,Nm5-1},其中Nm5表示集合M5中元素的个数;将集合M5按照pi从小到大进行排序,并根据式(26)计算M5对应的平均灰度值pha;
<mrow>
<msub>
<mi>p</mi>
<mrow>
<mi>h</mi>
<mi>a</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>m</mi>
<mn>5</mn>
</mrow>
</msub>
<mo>/</mo>
<mn>4</mn>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>m</mi>
<mn>5</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>N</mi>
<mrow>
<mi>m</mi>
<mn>5</mn>
</mrow>
</msub>
<mo>/</mo>
<mn>4</mn>
</mrow>
</munderover>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mrow>
<mi>m</mi>
<mn>5</mn>
</mrow>
</msub>
<mo>/</mo>
<mn>2</mn>
</mrow>
</mfrac>
<mo>;</mo>
<msub>
<mi>p</mi>
<mi>i</mi>
</msub>
<mo>&Element;</mo>
<msub>
<mi>M</mi>
<mn>5</mn>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>26</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤3.3.3:依次对Dh中坐标(x,y)满足式(17)的像素的灰度值根据式(27)进行调整;
<mrow>
<msub>
<mi>P</mi>
<msub>
<mi>D</mi>
<mi>h</mi>
</msub>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>P</mi>
<msub>
<mi>D</mi>
<mi>h</mi>
</msub>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>p</mi>
<mrow>
<mi>h</mi>
<mi>b</mi>
</mrow>
</msub>
<msub>
<mi>p</mi>
<mrow>
<mi>h</mi>
<mi>a</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>27</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤3.4:对调整后的图像Dh,Ds,Dv进行整合,得到去除阴影后的车牌彩图DR。
2.根据权利要求1所述的一种基于HSV的车牌阴影检测和去除方法,其特征在于步骤1)中的对获取的车牌阴影图像D进行图像预处理的具体过程如下:
步骤1.1:对图像D的S通道进行OTSU二值化和膨胀操作,得到二值化图D2;
步骤1.2:对图像D的V通道进行中值滤波操作,得到灰度图D3;
步骤1.3:根据公式(1)和(2)获得字符区域的灰度值集合E1和背景区域的灰度值集合E2,并利用OTSU计算E2的最佳阈值,设为thE2,然后根据公式(3)得到图像D3的二值化图像D4;
E1={PD3(x,y)|PD2(x,y)=255,0≤x<W,0≤y<H} (1)
E2={PD3(x,y)|PD2(x,y)=0,0≤x<W,0≤y<H} (2)
其中,表示图像D2中坐标(x,y)处的像素的灰度值,表示图像D3中坐标(x,y)处的像素的灰度值,表示图像D4中坐标(x,y)处的像素的灰度值;
步骤1.4:去除图像D4中轮廓面积小于A的小团块,得到用于阴影边缘检测的二值化图像Db,其中A表示预先设置的轮廓面积最小阈值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710309363.7A CN107292898B (zh) | 2017-05-04 | 2017-05-04 | 一种基于hsv的车牌阴影检测和去除方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710309363.7A CN107292898B (zh) | 2017-05-04 | 2017-05-04 | 一种基于hsv的车牌阴影检测和去除方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107292898A true CN107292898A (zh) | 2017-10-24 |
CN107292898B CN107292898B (zh) | 2019-09-10 |
Family
ID=60094347
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710309363.7A Active CN107292898B (zh) | 2017-05-04 | 2017-05-04 | 一种基于hsv的车牌阴影检测和去除方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107292898B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109165643A (zh) * | 2018-08-21 | 2019-01-08 | 浙江工业大学 | 一种基于深度学习的车牌识别方法 |
CN113160311A (zh) * | 2020-01-07 | 2021-07-23 | 北京赛博联物科技有限公司 | 太阳能电池板调节方法及调节装置、追日系统和存储设备 |
CN113487502A (zh) * | 2021-06-30 | 2021-10-08 | 中南大学 | 一种坑洼图像的阴影去除方法 |
CN113689710A (zh) * | 2021-09-01 | 2021-11-23 | 深圳市博思凯电子有限公司 | 一种用在智慧停车无人值守的边缘计算方法及系统 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101789080A (zh) * | 2010-01-21 | 2010-07-28 | 上海交通大学 | 车辆车牌实时定位字符分割的检测方法 |
CN102521587A (zh) * | 2011-11-24 | 2012-06-27 | 华中科技大学 | 一种车牌定位方法 |
CN103268493A (zh) * | 2013-05-29 | 2013-08-28 | 上海电机学院 | Rgb格式的汽车牌照图像定位方法 |
CN103488978A (zh) * | 2013-09-26 | 2014-01-01 | 浙江工业大学 | 一种基于灰度跳变和字符投影间隔模式的车牌定位方法 |
CN103632351A (zh) * | 2013-12-16 | 2014-03-12 | 武汉大学 | 一种基于亮度基准漂移的全天候交通图像增强方法 |
CN104036262A (zh) * | 2014-06-30 | 2014-09-10 | 南京富士通南大软件技术有限公司 | 一种lpr车牌筛选识别的方法和系统 |
CN104134079A (zh) * | 2014-07-31 | 2014-11-05 | 中国科学院自动化研究所 | 一种基于极值区域和极限学习机的车牌识别方法 |
CN105117726A (zh) * | 2015-08-07 | 2015-12-02 | 南京富士通南大软件技术有限公司 | 基于多特征区域累积的车牌定位方法 |
CN105488797A (zh) * | 2015-11-25 | 2016-04-13 | 安徽创世科技有限公司 | 一种hsv空间的车牌定位方法 |
CN105913058A (zh) * | 2016-04-05 | 2016-08-31 | 浙江工业大学 | 一种多车牌定位方法 |
CN106203433A (zh) * | 2016-07-13 | 2016-12-07 | 西安电子科技大学 | 一种车辆监控图像中车牌位置自动提取并透视校正的方法 |
CN106529533A (zh) * | 2016-11-07 | 2017-03-22 | 湖南源信光电科技有限公司 | 一种基于多尺度分析和匹配排序的复杂天气车牌定位方法 |
-
2017
- 2017-05-04 CN CN201710309363.7A patent/CN107292898B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101789080A (zh) * | 2010-01-21 | 2010-07-28 | 上海交通大学 | 车辆车牌实时定位字符分割的检测方法 |
CN102521587A (zh) * | 2011-11-24 | 2012-06-27 | 华中科技大学 | 一种车牌定位方法 |
CN103268493A (zh) * | 2013-05-29 | 2013-08-28 | 上海电机学院 | Rgb格式的汽车牌照图像定位方法 |
CN103488978A (zh) * | 2013-09-26 | 2014-01-01 | 浙江工业大学 | 一种基于灰度跳变和字符投影间隔模式的车牌定位方法 |
CN103632351A (zh) * | 2013-12-16 | 2014-03-12 | 武汉大学 | 一种基于亮度基准漂移的全天候交通图像增强方法 |
CN104036262A (zh) * | 2014-06-30 | 2014-09-10 | 南京富士通南大软件技术有限公司 | 一种lpr车牌筛选识别的方法和系统 |
CN104134079A (zh) * | 2014-07-31 | 2014-11-05 | 中国科学院自动化研究所 | 一种基于极值区域和极限学习机的车牌识别方法 |
CN105117726A (zh) * | 2015-08-07 | 2015-12-02 | 南京富士通南大软件技术有限公司 | 基于多特征区域累积的车牌定位方法 |
CN105488797A (zh) * | 2015-11-25 | 2016-04-13 | 安徽创世科技有限公司 | 一种hsv空间的车牌定位方法 |
CN105913058A (zh) * | 2016-04-05 | 2016-08-31 | 浙江工业大学 | 一种多车牌定位方法 |
CN106203433A (zh) * | 2016-07-13 | 2016-12-07 | 西安电子科技大学 | 一种车辆监控图像中车牌位置自动提取并透视校正的方法 |
CN106529533A (zh) * | 2016-11-07 | 2017-03-22 | 湖南源信光电科技有限公司 | 一种基于多尺度分析和匹配排序的复杂天气车牌定位方法 |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109165643A (zh) * | 2018-08-21 | 2019-01-08 | 浙江工业大学 | 一种基于深度学习的车牌识别方法 |
CN109165643B (zh) * | 2018-08-21 | 2021-11-16 | 浙江工业大学 | 一种基于深度学习的车牌识别方法 |
CN113160311A (zh) * | 2020-01-07 | 2021-07-23 | 北京赛博联物科技有限公司 | 太阳能电池板调节方法及调节装置、追日系统和存储设备 |
CN113160311B (zh) * | 2020-01-07 | 2024-09-17 | 北京赛博联物科技有限公司 | 太阳能电池板调节方法及调节装置、追日系统和存储设备 |
CN113487502A (zh) * | 2021-06-30 | 2021-10-08 | 中南大学 | 一种坑洼图像的阴影去除方法 |
CN113487502B (zh) * | 2021-06-30 | 2022-05-03 | 中南大学 | 一种坑洼图像的阴影去除方法 |
CN113689710A (zh) * | 2021-09-01 | 2021-11-23 | 深圳市博思凯电子有限公司 | 一种用在智慧停车无人值守的边缘计算方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN107292898B (zh) | 2019-09-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106651872B (zh) | 基于Prewitt算子的路面裂缝识别方法及系统 | |
CN105809138B (zh) | 一种基于分块识别的道路警告标志检测与识别方法 | |
CN107230202B (zh) | 路面病害图像的自动识别方法和系统 | |
CN107610114B (zh) | 基于支持向量机的光学卫星遥感影像云雪雾检测方法 | |
CN105373794B (zh) | 一种车牌识别方法 | |
CN102375982B (zh) | 一种融合多字符特征的车牌定位方法 | |
CN107123088B (zh) | 一种自动更换证件照背景颜色的方法 | |
CN107292898B (zh) | 一种基于hsv的车牌阴影检测和去除方法 | |
CN103886589B (zh) | 面向目标的自动化高精度边缘提取方法 | |
CN105809121A (zh) | 多特征协同的交通标志检测与识别方法 | |
CN107103317A (zh) | 基于图像融合和盲反卷积的模糊车牌图像识别算法 | |
CN110706235B (zh) | 一种基于双级级联分割的远红外行人检测方法 | |
CN109086687A (zh) | 基于pca降维的hog-mblbp融合特征的交通标志识别方法 | |
CN102306384B (zh) | 一种基于单幅图像的颜色恒常性处理方法 | |
CN110866879B (zh) | 一种基于多密度雨纹感知的图像去雨方法 | |
CN106529532A (zh) | 一种基于积分特征通道与灰度投影的车牌识别系统 | |
CN104766096A (zh) | 一种基于多尺度全局特征和局部特征的图像分类方法 | |
CN109948625A (zh) | 文本图像清晰度评估方法及系统、计算机可读存储介质 | |
CN102915544A (zh) | 基于纹理检测和颜色分割的视频图像运动目标提取方法 | |
CN107154044B (zh) | 一种中餐食物图像的分割方法 | |
CN103824081A (zh) | 一种室外恶劣光照条件下的快速鲁棒交通标志检测方法 | |
CN113673541B (zh) | 一种用于目标检测的图像样本生成方法及应用 | |
CN105678218A (zh) | 一种移动物体分类的方法 | |
CN109886168B (zh) | 一种基于层阶的地面交通标志识别方法 | |
CN104951795A (zh) | 图像分类识别判断方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |