发明名称 一种基于核密度估计的车牌字符分割方法
摘要 本发明涉及一种基于核密度估计的车牌字符分割方法。目前方法在极端情况下存在较多无法二值化好的情况。本发明首先对车牌图像进行归一化、锐化边缘和去除噪声的预处理;其次确定车牌字符区域;然后找出当前核概率密度曲线中,出现概率最大的像素值和该点的核密度半宽并利用这两个参数将图像二值化。最后取出在宽度范围中间得分最高的分割结果作为最终的字符分割的结果。本发明能够准确的确定字符像素的分布,抗环境干扰能力较一般的迭代二值化方法性能大大增强。
申请公布号 CN103207997B 申请公布日期 2017.02.08
申请号 CN201310131732.X 申请日期 2013.04.15
申请人 浙江捷尚视觉科技股份有限公司 发明人 尚凌辉;蒋宗杰;王弘玥;高勇
分类号 G06K9/34(2006.01)I 主分类号 G06K9/34(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 杜军
主权项 一种基于核密度估计的车牌字符分割方法,其特征在于该方法的具体步骤如下:步骤(1)对车牌图像进行归一化、锐化边缘和去除噪声的预处理;步骤(2)利用水平方向梯度投影对预处理后的图像进行车牌上下边缘的分割;计算每个水平块A<sup>i</sup>上的梯度值<img file="FDA0001015721450000011.GIF" wi="101" he="65" />高度值<img file="FDA0001015721450000012.GIF" wi="68" he="60" />及所有水平块梯度均值<img file="FDA0001015721450000013.GIF" wi="93" he="78" />所有水平块的梯度方差<img file="FDA0001015721450000014.GIF" wi="75" he="71" />如果水平块不满足以下两个条件,则判定该水平块不是车牌字符区域,予以删除:条件1)<img file="FDA0001015721450000015.GIF" wi="334" he="151" />条件2)<img file="FDA0001015721450000016.GIF" wi="200" he="62" />其中C<sub>a</sub>的取值范围为[0,1];步骤(3)利用核密度估计迭代计算水平块的核概率密度曲线G(x),找出当前核概率密度曲线中,出现概率最大的像素值<img file="FDA0001015721450000017.GIF" wi="73" he="60" />和该像素值所在点的核密度半宽L<sub>j</sub>;若核密度半宽不满足L<sub>j</sub>&lt;20,则该核概率密度曲线G(x)所表现的像素区域不是属于车牌字符的区域;重复本步骤,找出除<img file="FDA0001015721450000018.GIF" wi="74" he="61" />外的最大像素值<img file="FDA0001015721450000019.GIF" wi="89" he="62" />和核密度半宽L<sub>j+1</sub>,如果半宽满足L<sub>j</sub>&lt;20,则执行步骤(4);步骤(4)将水平块A<sup>i</sup>根据概率最大的像素值<img file="FDA00010157214500000110.GIF" wi="76" he="59" />和该像素值的核密度半宽L<sub>j</sub>做图像二值化,对二值化的图像结果进行垂直投影分割;如字符块不满足以下任何一条件,则将该像素值所表示的核予以剔除,并跳转至步骤(3):a)字符宽度W<sub>j</sub>&lt;4;b)分割得到的字符总数小于5个;步骤(5)设字符宽度范围为W<sub>min</sub>至W<sub>max</sub>,利用该宽度所得到的标准车牌模板和二值化图像垂直投影后的结果进行匹配,取出在[W<sub>min</sub>,W<sub>max</sub>]中间得分最高的分割结果作为最终的字符分割的结果;步骤(3)在计算核密度半宽时,设两个高斯核函数K<sub>1</sub>(u),K<sub>2</sub>(u),分别为正态分布N(0,1),N(0,4);在车牌上像素在3000个以上时,两个核函数K<sub>1</sub>(u),K<sub>2</sub>(u)计算得到的概率密度最小时,取到的半宽L<sub>j</sub>为当前真实的车牌字符像素分布半宽;即计算下式的最小值:<maths num="0001"><math><![CDATA[<mrow><mi>M</mi><mi>I</mi><mi>N</mi><mrow><mo>(</mo><mi>L</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><msup><mi>n</mi><mn>2</mn></msup><mi>L</mi><msqrt><mi>&pi;</mi></msqrt></mrow></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mo>&lsqb;</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mi>e</mi><mrow><mo>-</mo><mfrac><msup><mrow><mo>(</mo><mrow><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>x</mi><mi>j</mi></msub></mrow><mo>)</mo></mrow><mn>2</mn></msup><mrow><mn>4</mn><msup><mi>L</mi><mn>2</mn></msup></mrow></mfrac></mrow></msup><mo>+</mo><mfrac><mn>1</mn><mn>4</mn></mfrac><msup><mi>e</mi><mrow><mo>-</mo><mfrac><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>x</mi><mi>j</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mrow><mn>16</mn><msup><mi>L</mi><mn>2</mn></msup></mrow></mfrac></mrow></msup><mo>-</mo><mfrac><mn>2</mn><msqrt><mn>10</mn></msqrt></mfrac><msup><mi>e</mi><mrow><mo>-</mo><mfrac><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>x</mi><mi>j</mi></msub><mo>)</mo></mrow><mn>2</mn></msup><mrow><mn>10</mn><msup><mi>L</mi><mn>2</mn></msup></mrow></mfrac></mrow></msup><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0001015721450000021.GIF" wi="1444" he="159" /></maths>其中n表示水平块A<sup>i</sup>中的像素总数,L表示可能核密度半宽,x<sub>i</sub>,x<sub>j</sub>表示水平块内的像素值;计算得到是该式最小的核密度半宽L<sub>j</sub>,代入计算核密度函数:<img file="FDA0001015721450000022.GIF" wi="494" he="121" />其中K<sub>σ</sub>是核密度半宽L<sub>j</sub>的高斯核;步骤(5)是对步骤(4)得到的二值化图像区域做垂直投影,利用二值化字符得到的宽度[W<sub>min</sub>,W<sub>max</sub>]和标准车牌模板的比例换算到当前模板车牌的宽度[PlateW<sub>min</sub>,PlateW<sub>max</sub>];用换算的车牌模板计算可能的分割位置[pos<sub>0</sub>,pos<sub>1</sub>,...pos<sub>m</sub>];根据如下公式计算每个可能分割位置的得分:<maths num="0002"><math><![CDATA[<mrow><msub><mi>Score</mi><mi>p</mi></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mn>1</mn><mo>-</mo><mn>1</mn></mrow></munderover><mfrac><mrow><msub><mi>EC</mi><mi>i</mi></msub></mrow><mrow><msub><mi>RC</mi><mi>i</mi></msub></mrow></mfrac><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mn>2</mn><mo>-</mo><mn>1</mn></mrow></munderover><mfrac><mrow><msub><mi>RC</mi><mi>j</mi></msub></mrow><mrow><msub><mi>RI</mi><mi>j</mi></msub></mrow></mfrac><mo>,</mo><mi>p</mi><mo>&Element;</mo><mo>&lsqb;</mo><mn>0</mn><mo>,</mo><mi>m</mi><mo>&rsqb;</mo><mo>;</mo></mrow>]]></math><img file="FDA0001015721450000023.GIF" wi="902" he="143" /></maths>其中EC<sub>i</sub>表示第i个二值化字符面积和标准车牌第i个字符面积相交大小,RC<sub>i</sub>表示第i个字符的面积,RC<sub>j</sub>表示模板第i和i+1个二值化字符之间的面积,RI<sub>j</sub>表示模板第i和i+1个二值化字符之间的面积,Score<sub>p</sub>表示第p个可能分割位置的分数;选择所有可能分割位置中分数的最大值作为最优化的车牌分割位置。
地址 311121 浙江省杭州市余杭区五常街道文一西路998号7幢东楼