发明名称 基于形状特征不变子空间的交通标志识别方法
摘要 本发明涉及一种基于形状特征不变子空间的交通标志识别方法,以交通标志的二值图像为特征提取对象,将主成分分析方法和线性判别分析方法相结合,先利用主成分分析方法对交通标志图像进行特征提取,得到具有最佳描述效果的特征矩阵;再在该矩阵上用线性判别分析方法进行二次特征提取,得到具有最佳分类效果的特征矩阵,所以本方法所提取的特征既具有最佳描述性也具有最佳区分性;最后采用最小距离分类法对所提取特征进行识别,通过实验验证可以将交通标志准确的识别出来。
申请公布号 CN103761538A 申请公布日期 2014.04.30
申请号 CN201410012899.9 申请日期 2014.01.10
申请人 沈阳工业大学 发明人 张志佳;何纯静;李雅红;崔世昊
分类号 G06K9/64(2006.01)I;G06K9/46(2006.01)I 主分类号 G06K9/64(2006.01)I
代理机构 沈阳智龙专利事务所(普通合伙) 21115 代理人 宋铁军;周楠
主权项 1.一种基于形状特征不变子空间的交通标志识别方法,其特征在于:步骤如下:(1)建立交通标志图像数据库:针对交通标志主标志进行识别,分别建立指示交通标志图像数据库、警告交通标志图像数据库和禁令交通标志图像数据库,分为以下几步:1)下载标准交通标志的彩色图像,建立相应的彩色图像数据库;2)利用公式gray=0.299R+0.587G+0.114B将交通标志图像转化为灰度图像,建立相应的灰度图像数据库;3)对图像加入噪声;4)对灰度图像进行二值化,建立交通标志的形状数据库;(2)计算主成分分析(PCA)特征子空间:假设交通标志图像集合为:X={f<sub>1</sub>(x,y),f<sub>2</sub>(x,y),...,f<sub>N</sub>(x,y)},每一幅图像f<sub>i</sub>(x,y)都能够按照行顺序展开成一个M维的向量<img file="FDA0000455544350000011.GIF" wi="496" he="110" />M为图像像素个数,它的协方差矩阵定义为:C=XX<sup>T</sup>=E{(X-u)(X-u)<sup>T</sup>},其中u=E{X};对C进行奇异值分解:<img file="FDA0000455544350000012.GIF" wi="1701" he="331" />矩阵U列向量u<sub>1</sub>,u<sub>2</sub>,...,u<sub>M</sub>是单位正交的,构成主成分空间的基,它们被称之为主成分;KL变换把协方差矩阵对角化了,即通过KL变换,消除了原有向量<img file="FDA0000455544350000013.GIF" wi="47" he="59" />的各分量之间的相关性,从而达到降低特征空间维数的目的;设图像向量<img file="FDA0000455544350000014.GIF" wi="54" he="72" />在主成分空间中坐标为a<sub>i</sub>=[a<sub>1i</sub>,a<sub>2i</sub>,...,a<sub>Mi</sub>]<sup>T</sup>∈R<sup>M</sup>,则有<maths num="0001"><![CDATA[<math><mrow><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>=</mo><msub><mi>Ua</mi><mi>i</mi></msub><mo>=</mo><mo>[</mo><msub><mi>u</mi><mn>1</mn></msub><mo>,</mo><msub><mi>u</mi><mn>2</mn></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>u</mi><mi>M</mi></msub><mo>]</mo><mo>&CenterDot;</mo><mo>[</mo><msup><mrow><mo>[</mo><msub><mi>a</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mo>,</mo><msub><mi>a</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>a</mi><mi>Mi</mi></msub><mo>]</mo></mrow><mi>T</mi></msup><mo>]</mo><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>a</mi><mi>ki</mi></msub><mo>&CenterDot;</mo><msub><mi>u</mi><mi>k</mi></msub><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>因为U为单位正交矩阵,UU<sup>T</sup>=U<sup>T</sup>U=E,则<img file="FDA0000455544350000021.GIF" wi="239" he="90" /><maths num="0002"><![CDATA[<math><mrow><msup><mrow><mo>[</mo><msub><mi>a</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mo>,</mo><msub><mi>a</mi><mrow><mn>2</mn><mi>i</mi></mrow></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>a</mi><mi>Mi</mi></msub><mo>]</mo></mrow><mi>T</mi></msup><mo>=</mo><msup><mrow><mo>[</mo><msub><mi>u</mi><mn>1</mn></msub><mo>,</mo><msub><mi>u</mi><mn>2</mn></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>u</mi><mi>M</mi></msub><mo>]</mo></mrow><mi>T</mi></msup><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>=</mo><mo>[</mo><msubsup><mi>u</mi><mn>1</mn><mi>T</mi></msubsup><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>,</mo><msubsup><mi>u</mi><mn>2</mn><mi>T</mi></msubsup><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msubsup><mi>u</mi><mi>M</mi><mi>T</mi></msubsup><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>]</mo><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msub><mi>a</mi><mi>ki</mi></msub><mo>=</mo><msubsup><mi>u</mi><mi>k</mi><mi>T</mi></msubsup><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>,</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>代入上式得到:<maths num="0004"><![CDATA[<math><mrow><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><mo>{</mo><msubsup><mi>u</mi><mi>k</mi><mi>T</mi></msubsup><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><msub><mi>u</mi><mi>k</mi></msub><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>它物理意义表示为:使用主成分向量的加权累加和来拟合输入图像<img file="FDA0000455544350000025.GIF" wi="70" he="73" />假设只用前d个有限项来估计<img file="FDA0000455544350000026.GIF" wi="77" he="83" /><maths num="0005"><![CDATA[<math><mrow><msub><mover><mi>x</mi><mo>^</mo></mover><mi>i</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>d</mi></munderover><msubsup><mi>u</mi><mi>k</mi><mi>T</mi></msubsup><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mi>i</mi></msub><msub><mi>u</mi><mi>k</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>这样就完成了对原图像的降维,由于KL变换所具有的优良特性,认为降维后的特征空间足以代表原图像的本质特征;取前d个特征值所对应的特征向量作为PCA特征子空间的基,将所有交通标志图像向该子空间投影,就得到了能够表示原图像集的特征子空间;(3)计算线性判别分析(LDA)特征子空间:经典的线性判别分析中使用的是Fisher准则函数,所以线性判决分析又被称为Fisher线性判别分析(Fisher Linear Discriminant Analysis/FLDA);Fisher准则函数是这样定义的:<maths num="0006"><![CDATA[<math><mrow><mi>J</mi><mrow><mo>(</mo><mi>w</mi><mo>)</mo></mrow><mo>=</mo><mi>arg</mi><munder><mi>max</mi><mi>w</mi></munder><mfrac><mrow><mo>|</mo><msup><mi>w</mi><mi>T</mi></msup><msub><mi>S</mi><mi>B</mi></msub><mi>w</mi><mo>|</mo></mrow><mrow><mo>|</mo><msup><mi>w</mi><mi>T</mi></msup><msub><mi>S</mi><mi>w</mi></msub><mi>w</mi><mo>|</mo></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,类间离散度S<sub>B</sub>和类内离散度S<sub>w</sub>定义为:<maths num="0007"><![CDATA[<math><mrow><msub><mi>S</mi><mi>B</mi></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mi>C</mi><mrow><mo>(</mo><mi>C</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>C</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>C</mi></munderover><mrow><mo>(</mo><msub><mi>u</mi><mi>i</mi></msub><mo>-</mo><msub><mi>u</mi><mi>j</mi></msub><mo>)</mo></mrow><msup><mrow><mo>(</mo><msub><mi>u</mi><mi>i</mi></msub><mo>-</mo><msub><mi>u</mi><mi>j</mi></msub><mo>)</mo></mrow><mi>T</mi></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><msub><mi>S</mi><mi>B</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi>C</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>C</mi></munderover><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>i</mi></msub></munderover><mrow><mo>(</mo><msubsup><mi>x</mi><mi>j</mi><mi>i</mi></msubsup><mo>-</mo><msub><mi>u</mi><mi>i</mi></msub><mo>)</mo></mrow><msup><mrow><mo>(</mo><msubsup><mi>x</mi><mi>j</mi><mi>i</mi></msubsup><mo>-</mo><msub><mi>u</mi><mi>i</mi></msub><mo>)</mo></mrow><mi>T</mi></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow></math>]]></maths>C为总的类数,N<sub>i</sub>表示第i类的样本数,<img file="FDA00004555443500000211.GIF" wi="58" he="82" />为第i类中的第j个样本,<img file="FDA00004555443500000212.GIF" wi="215" he="146" />为第i类中的样本均值;在数学上,求解Fisher准则函数的最优解就等于求解<img file="FDA0000455544350000031.GIF" wi="134" he="78" />的特征值问题;求使<img file="FDA0000455544350000032.GIF" wi="126" he="78" />取得最大值时的特征值,所对应的特征向量就是LDA子空间的基;将PCA特征子空间中的特征集向LDA子空间投影,就得到最终的代表原交通标志图像的PLA特征空间;(4)采用最小距离分类法验证所特征提取特征的识别效果。
地址 110870 辽宁省沈阳市经济技术开发区沈辽西路111号