发明名称 一种采用集合距离的交通路标识别方法
摘要 本发明给出了一种采用集合距离的交通路标识别方法。本发明首先建立标准图像模板库。模板库的建立是通过选取能够代表每幅标准图像的标志点,然后分别计算每两个标志点之间的距离和他们所在直线与正x轴间的角度。其次,观测图像特征点的得到。该步利用马尔科夫随机场,极大后验估计准则和ICM算法对观测图像进行去噪得到其边缘。为减少计算量,提高实时性,去掉边缘点中角度相近,距离在一定范围内的点,利用得到的这些点作为标志点计算特征点。尽管如此,观测图像与标准图像特征点数差别还是较大,所以采用可以度量元素值和元素个数的点集距离计算观测图像与标准图像的差别。本发明能有效减少交通路标识别处理时间,同时提高了识别的准确性。
申请公布号 CN106529482A 申请公布日期 2017.03.22
申请号 CN201611000969.4 申请日期 2016.11.14
申请人 叶瀚礼 发明人 叶瀚礼
分类号 G06K9/00(2006.01)I 主分类号 G06K9/00(2006.01)I
代理机构 杭州奥创知识产权代理有限公司 33272 代理人 王佳健
主权项 一种采用集合距离的交通路标识别方法,其特征在于该方法包括以下步骤:步骤1.路标模板库建立与特征点选取;选择具有典型路标的标准图像,建立路标模板库,对每一个路标图像进行图像大小尺寸的标准化;针对不同形状的路标,选取适当数量的标志点,在得到了标准图像的标志点后,特征点的计算采取欧氏距离,计算公式如下:<maths num="0001"><math><![CDATA[<mrow><mi>d</mi><mo>=</mo><msqrt><mrow><msup><mrow><mo>(</mo><msub><mi>x</mi><mn>1</mn></msub><mo>-</mo><msub><mi>x</mi><mn>2</mn></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><msub><mi>y</mi><mn>1</mn></msub><mo>-</mo><msub><mi>y</mi><mn>2</mn></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001152523890000011.GIF" wi="1525" he="103" /></maths>其中(x<sub>1</sub>,y<sub>1</sub>),(x<sub>2</sub>,y<sub>2</sub>)分别是两个标志点的位置,这样便组成了一个特征点;为增加识别的准确性添加每条直线的角度θ,计算公式如下:<img file="FDA0001152523890000012.GIF" wi="1518" he="222" />每个特征点由两个量组成,即距离和该特征点对应的角度;将每幅标准图像表示为一个二维的矩阵std;一幅图像表示如下:<maths num="0002"><math><![CDATA[<mrow><mi>s</mi><mi>t</mi><mi>d</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>d</mi><mn>1</mn></msub></mtd><mtd><msub><mi>d</mi><mn>2</mn></msub></mtd><mtd><mo>...</mo></mtd><mtd><msub><mi>d</mi><mi>m</mi></msub></mtd></mtr><mtr><mtd><msub><mi>&theta;</mi><mn>1</mn></msub></mtd><mtd><msub><mi>&theta;</mi><mn>2</mn></msub></mtd><mtd><mo>...</mo></mtd><mtd><msub><mi>&theta;</mi><mi>m</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001152523890000013.GIF" wi="1488" he="174" /></maths>其中,m为该图像中所选特征点数;通过计算各标志点之间的欧氏距离和标志点之间的角度,把距离和角度组成的矩阵作为标准图像的特征点;这样便形成了标准图像的模板库;步骤2.路标图像去噪处理;对于采集到的路标图像,先对其进行灰度化,使其成为二维图像,这样便可将其看作二维的随机场;用F表示该随机场,记f<sub>s</sub>是F<sub>s</sub>的实现,其中s是像素点;假设像素点s的邻域是N<sub>s</sub>,设S是有限的M*N的点集,S={s=(i,j)|1≤i≤M,1≤j≤N};若满足以下两个条件:P(F<sub>s</sub>=f<sub>s</sub>)&gt;0    (4)P(F<sub>s</sub>=f<sub>s</sub>|F<sub>r</sub>=f<sub>r</sub>,r∈S,r≠s)=P(F<sub>s</sub>=f<sub>s</sub>|F<sub>r</sub>=f<sub>r</sub>,r∈N<sub>s</sub>)    (5)则称F为关于邻域N<sub>s</sub>的马尔科夫随机场;依据马尔科夫随机场与Gibbs场的等价性,采用符合Gibbs分布随机场,即Gibbs场来描述图像;关于邻域N<sub>s</sub>的Gibbs分布表示为:<maths num="0003"><math><![CDATA[<mrow><mi>P</mi><mrow><mo>(</mo><mi>F</mi><mo>=</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>Z</mi></mfrac><msup><mi>e</mi><mrow><mo>-</mo><mi>U</mi><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001152523890000021.GIF" wi="1422" he="136" /></maths>上式中U(f)称为能量函数<maths num="0004"><math><![CDATA[<mrow><mi>U</mi><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><munder><mo>&Sigma;</mo><mrow><mi>c</mi><mo>&Element;</mo><mi>C</mi></mrow></munder><msub><mi>V</mi><mi>c</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001152523890000022.GIF" wi="1342" he="103" /></maths>C表示邻域N<sub>s</sub>中包含的集簇的集合,V<sub>c</sub>(f)称为关于集簇C的势函数,Z为正则化因子;马尔科夫随机场中,相应的集簇势函数定义如下:<img file="FDA0001152523890000023.GIF" wi="1470" he="175" />若观测图像对应的随机场为Y,那么后验概率表示为:<maths num="0005"><math><![CDATA[<mrow><mi>P</mi><mrow><mo>(</mo><mi>F</mi><mo>|</mo><mi>Y</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>P</mi><mrow><mo>(</mo><mi>Y</mi><mo>|</mo><mi>F</mi><mo>)</mo></mrow><mi>P</mi><mrow><mo>(</mo><mi>F</mi><mo>)</mo></mrow></mrow><mrow><mi>P</mi><mrow><mo>(</mo><mi>Y</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001152523890000024.GIF" wi="1453" he="127" /></maths>上式中,P(Y)为常数,P(F)为观测图像的先验概率,P(F|Y)为已知Y时关于原图像的后验概率,应用贝叶斯准则,观测图像的恢复可转化为求解最大后验概率条件下对原图像的估计值;对于给定观测路标图像Y,P(Y=y)是常数,路标图像的观测模型如下式所示:y<sub>i</sub>=f<sub>i</sub>+e<sub>i</sub>    (10)其中y<sub>i</sub>是观测图像在点i处得像素值,f<sub>i</sub>为标准图像在点i处的像素值,e<sub>i</sub>为噪声;则后验概率表示为下式:<maths num="0006"><math><![CDATA[<mrow><mi>P</mi><mrow><mo>(</mo><mi>F</mi><mo>|</mo><mi>Y</mi><mo>)</mo></mrow><mo>&cong;</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>U</mi><mrow><mo>(</mo><mi>y</mi><mo>|</mo><mi>f</mi><mo>)</mo></mrow></mrow></msup><mo>=</mo><msup><mi>e</mi><mrow><mo>-</mo><mrow><mo>(</mo><mi>h</mi><munder><mi>&Sigma;</mi><mi>s</mi></munder><msub><mi>f</mi><mi>s</mi></msub><mo>-</mo><mi>&beta;</mi><munder><mi>&Sigma;</mi><msub><mi>N</mi><mi>s</mi></msub></munder><mi>&delta;</mi><mo>(</mo><msub><mi>f</mi><mi>s</mi></msub><mo>,</mo><msub><mi>f</mi><mi>r</mi></msub><mo>)</mo></mrow><mo>-</mo><mi>&eta;</mi><munder><mi>&Sigma;</mi><mi>s</mi></munder><mi>&delta;</mi><mo>(</mo><mrow><msub><mi>f</mi><mi>s</mi></msub><mo>,</mo><msub><mi>y</mi><mi>s</mi></msub></mrow><mo>)</mo><mo>)</mo></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001152523890000025.GIF" wi="1518" he="91" /></maths><maths num="0007"><math><![CDATA[<mrow><mi>&delta;</mi><mrow><mo>(</mo><msub><mi>f</mi><mi>s</mi></msub><mo>,</mo><msub><mi>f</mi><mi>r</mi></msub><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><mn>1</mn><mo>,</mo><msub><mi>f</mi><mi>s</mi></msub><mo>=</mo><msub><mi>f</mi><mi>r</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mn>1</mn><mo>,</mo><msub><mi>f</mi><mi>s</mi></msub><mo>&NotEqual;</mo><msub><mi>f</mi><mi>r</mi></msub></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001152523890000026.GIF" wi="1382" he="143" /></maths>其中,h,β,η为调节参数,f<sub>s</sub>,f<sub>r</sub>为观测图像s,r像素点的标志值,y<sub>s</sub>为估计图像s像素点的标志值;求解原图像的最大后验概率问题可转化为求解能量函数U(f|y)的最小值;即<img file="FDA0001152523890000031.GIF" wi="1541" he="103" />求解采用ICM算法优化得到;步骤3.点集集合距离识别;由步骤1得到了标准图像的数据,由步骤2得到了观测图像去噪后图像,利用ICM算法得到的标志变量求得了观测图像的边缘点;将角度相近、距离在设定范围内的边缘点合并,以此来减少特征点数,降低OSPA距离的计算复杂度,提高识别的实时性;由采样得到的标志点来计算观测图像的特征点;为了能够得到标准图像与观测图像之间的相似程度,采用了点集集合计算观测图像与标准图像的差别,其定义如下:<maths num="0008"><math><![CDATA[<mrow><mi>O</mi><mi>S</mi><mi>P</mi><mi>A</mi><mrow><mo>(</mo><mi>p</mi><mo>,</mo><mi>c</mi><mo>,</mo><mi>X</mi><mo>,</mo><mover><mi>X</mi><mo>^</mo></mover><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><mrow><mfrac><mn>1</mn><mi>n</mi></mfrac><mrow><mo>(</mo><mrow><munder><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow><mrow><mi>&pi;</mi><mo>&Element;</mo><msub><mi>&Pi;</mi><mi>&pi;</mi></msub></mrow></munder><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mi>d</mi><msup><mrow><mo>(</mo><mrow><mi>c</mi><mo>,</mo><msub><mi>x</mi><mi>i</mi></msub><mo>,</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>&pi;</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow></msub></mrow><mo>)</mo></mrow><mi>p</mi></msup><mo>+</mo><msup><mi>c</mi><mi>p</mi></msup><mrow><mo>(</mo><mrow><mi>m</mi><mo>-</mo><mi>n</mi></mrow><mo>)</mo></mrow></mrow><mo>)</mo></mrow></mrow><mo>)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mi>p</mi></mrow></msup><mo>,</mo><mi>n</mi><mo>&le;</mo><mi>m</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>14</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001152523890000032.GIF" wi="1620" he="126" /></maths>式(14)中的p代表阶数,c代表一个水平参数,用于调节集合势的估计误差比重,X,<img file="FDA0001152523890000033.GIF" wi="46" he="62" />分别代表标准路标图像和观测路标图像特征点矩阵。
地址 443000 湖北省宜昌市西陵区体育场路41号香山铭苑8-1-601