发明名称 基于特征向量定向的椭球形水果尺寸快速检测方法
摘要 本发明公开了一种基于特征向量定向的椭球形水果尺寸快速检测方法。对获取的水果图像进行阈值分割、滤波、边缘提取等操作来得到水果边缘图像。对边缘图像建立直角坐标系;求出边缘坐标的协方差矩阵;进而求得协方差矩阵的特征值与单位特征向量;利用单位特征向量定向水果,使其纵径或横径方向平行于直角坐标系水平轴,再利用计算水果边界的上、下、左、右极值点来完成尺寸检测。本发明通过对水果图像的边界坐标信息求取特征向量来达到快速定向到水果纵径和横径,避免了MER方法通过多次旋转水果图像带来的大量运算,在保证检测精度的同时,提高了检测速度,适于椭球形水果商品化处理过程中的水果尺寸实时检测需要。
申请公布号 CN105184775A 申请公布日期 2015.12.23
申请号 CN201510501000.4 申请日期 2015.08.14
申请人 浙江大学 发明人 饶秀勤;宋晨波;许济海;应义斌
分类号 G06T7/00(2006.01)I 主分类号 G06T7/00(2006.01)I
代理机构 杭州求是专利事务所有限公司 33200 代理人 林怀禹
主权项 一种基于特征向量定向的椭球形水果尺寸快速检测方法,获取水果图像,并经二值分割、滤波和边缘检测后,得到水果的边界E,建立直角坐标系;其特征在于该方法的步骤是:1)将边界E的数据以x<sub>i</sub>,y<sub>i</sub>,i=1,2…N,的形式存放到水果边界信息矩阵M,其中N为水果边界点总数;2)按公式(1)求水果边界信息矩阵M的协方差矩阵C:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>C</mi><mo>=</mo><mfenced open = '[' close = ']'><mtable><mtr><mtd><mfrac><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>&mu;</mi><mi>X</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></mfrac></mtd><mtd><mfrac><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>&mu;</mi><mi>X</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><msub><mi>&mu;</mi><mi>Y</mi></msub><mo>)</mo></mrow></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></mfrac></mtd></mtr><mtr><mtd><mfrac><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>&mu;</mi><mi>X</mi></msub><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><msub><mi>&mu;</mi><mi>Y</mi></msub><mo>)</mo></mrow></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></mfrac></mtd><mtd><mfrac><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>i</mi></msub><mo>-</mo><msub><mi>&mu;</mi><mi>Y</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></mfrac></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000781934000000011.GIF" wi="1891" he="309" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><msub><mi>&mu;</mi><mi>X</mi></msub><mo>=</mo><mfrac><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msub><mi>x</mi><mi>i</mi></msub></mrow><mi>n</mi></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>&mu;</mi><mi>Y</mi></msub><mo>=</mo><mfrac><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></msubsup><msub><mi>y</mi><mi>i</mi></msub></mrow><mi>n</mi></mfrac></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000781934000000012.GIF" wi="1899" he="306" /></maths>式中:μ<sub>X</sub>是边界E上x坐标的均值,求法见公式(2)μ<sub>Y</sub>是边界E上y坐标的均值,求法见公式(2)x<sub>i</sub>是边界E上第i个点的x坐标y<sub>i</sub>是边界E上第i个点的y坐标n是边界E上点的个数3)用公式(3)计算协方差矩阵C的特征值λ<sub>i</sub>,i=1,2|λ<sub>i</sub>E‑C|=0                 (3)用公式(4)求出特征值λ<sub>i</sub>,i=1,2的特征向量V<sub>i</sub>,i=1,2(λ<sub>i</sub>w‑C)V<sub>i</sub>=0,i=1,2              (4)用公式(5)计算单位长度的特征向量V0<sub>i</sub>,<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>V</mi><msub><mn>0</mn><mi>i</mi></msub><mo>=</mo><mfrac><msub><mi>V</mi><mi>i</mi></msub><mrow><mo>|</mo><msub><mi>V</mi><mi>i</mi></msub><mo>|</mo></mrow></mfrac><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000781934000000021.GIF" wi="1892" he="165" /></maths>用单位特征向量V0<sub>i</sub>,i=1,2组成定向矩阵R,即R=[V0<sub>1</sub>V0<sub>2</sub>];4)用公式(6)将水果边界信息矩阵M乘上定向矩阵R得到新的水果边界信息矩阵M’;<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><msubsup><mi>x</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>=</mo><msub><mi>x</mi><mi>i</mi></msub><msub><mi>R</mi><mn>11</mn></msub><mo>+</mo><msub><mi>y</mi><mi>i</mi></msub><msub><mi>R</mi><mn>21</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><msubsup><mi>y</mi><mi>i</mi><mo>&prime;</mo></msubsup><mo>=</mo><msub><mi>x</mi><mi>i</mi></msub><msub><mi>R</mi><mn>12</mn></msub><mo>+</mo><msub><mi>y</mi><mi>i</mi></msub><msub><mi>R</mi><mn>22</mn></msub></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000781934000000022.GIF" wi="1894" he="168" /></maths>式中:x<sub>i</sub>,y<sub>i</sub>为M中第i个元素R<sub>pq</sub>,p=1,2,q=1,2表示R中第p行第q列的元素x′<sub>i</sub>,y′<sub>i</sub>为M’中第i个元素5)在新的水果边界信息矩阵M’分别找出X’坐标和Y’坐标最小值x’<sub>min</sub>,y’<sub>min</sub>和最大值x’<sub>max</sub>,y’<sub>max</sub>;用公式(7)计算水果的纵径和横径:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mfenced open = '{' close = ''><mtable><mtr><mtd><mrow><msub><mi>D</mi><mi>max</mi></msub><mo>=</mo><mi>max</mi><mrow><mo>(</mo><msubsup><mi>x</mi><mi>max</mi><mo>&prime;</mo></msubsup><mo>-</mo><msubsup><mi>x</mi><mi>min</mi><mo>&prime;</mo></msubsup><mo>,</mo><msubsup><mi>y</mi><mi>max</mi><mo>&prime;</mo></msubsup><mo>-</mo><msubsup><mi>y</mi><mi>min</mi><mo>&prime;</mo></msubsup><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>D</mi><mi>min</mi></msub><mo>=</mo><mi>min</mi><mrow><mo>(</mo><msubsup><mi>x</mi><mi>max</mi><mo>&prime;</mo></msubsup><mo>-</mo><msubsup><mi>x</mi><mi>min</mi><mo>&prime;</mo></msubsup><mo>,</mo><msubsup><mi>y</mi><mi>max</mi><mo>&prime;</mo></msubsup><mo>-</mo><msubsup><mi>y</mi><mi>min</mi><mo>&prime;</mo></msubsup><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000781934000000023.GIF" wi="1894" he="191" /></maths>D<sub>max</sub>和D<sub>min</sub>分别代表水果的纵径和横径。
地址 310027 浙江省杭州市西湖区浙大路38号