发明名称 一种基于机器视觉的免疫层析结果判读方法
摘要 本发明公开一种基于机器视觉的免疫层析结果判读方法,其特征步骤是:一:采集免疫层析测试区显色特征区域及其周边非特征区的原始图像I,并获取免疫层析测试区显色特征区域的信息B;二:以非特征区的背景,采用最大似然法拟合出免疫层析测试区显色特征区域可能的背景图像信息C;三:根据背景差分方法,B与C差分,得到免疫层析测试区显色特征区域中前景的真实颜色信息D;四:根据息D统计出免疫层析测试区显色特征区域包含颜色信息熵的大小。其显著效果是:根据数字图像处理技术和仿人视觉技术,结合背景建模技术,剔除掉图像中的毛刺、瑕疵和噪声,精确提取出免疫层析测试区显色区域的颜色信息熵,以量化形式反应出显色的深浅程度。
申请公布号 CN102384972B 申请公布日期 2013.09.11
申请号 CN201010512048.2 申请日期 2010.10.19
申请人 刘江;李洲 发明人 刘江;李洲;张路
分类号 G01N33/558(2006.01)I 主分类号 G01N33/558(2006.01)I
代理机构 重庆市前沿专利事务所(普通合伙) 50211 代理人 余锦曦
主权项 1.一种基于机器视觉的免疫层析结果判读方法,其特征在于按以下步骤进行:步骤一:采集免疫层析测试区显色特征区域及其周边非特征区的原始图像I,并获取免疫层析测试区显色特征区域的信息B;步骤二:以非特征区的背景,采用最大似然法,拟合出免疫层析测试区显色特征区域可能的背景图像信息C;步骤三:根据背景差分方法,用所述免疫层析测试区显色特征区域的图像信息B与背景图像信息C进行差分,得到免疫层析测试区显色特征区域中前景的真实颜色信息D;步骤四:然后根据真实颜色信息D统计出免疫层析测试区显色特征区域包含信息熵的大小;所述原始图像I为:I={f(x,y)|(0<x<=width)&amp;(0<y<=height)}其中,f(x,y)代表图像中在(x,y)像素点处的RGB三维向量,width为图像的宽度,height为图像的高度;所述免疫层析测试区显色特征区域的信息B为:B={B<sup>i</sup>|i∈1…n}<maths num="0001"><![CDATA[<math><mrow><msup><mi>B</mi><mi>i</mi></msup><mo>=</mo><mo>{</mo><msup><mi>B</mi><mi>i</mi></msup><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>|</mo><msubsup><mi>x</mi><mi>left</mi><mi>i</mi></msubsup><mo>&lt;</mo><mi>x</mi><mo>&lt;</mo><msubsup><mi>x</mi><mi>right</mi><mi>i</mi></msubsup><mo>}</mo></mrow></math>]]></maths>B<sup>i</sup>(x)=(B<sup>i</sup>(x,1),B<sup>i</sup>(x,2),B<sup>i</sup>(x,3))B<sup>i</sup>(x,m)={f(x,y,m)|0<y<=height}其中i为免疫层析测试区显色特征区域的个数,<img file="FDA00003460406900021.GIF" wi="155" he="133" />为第i个免疫层析测试区显色特征区域的左边缘,<img file="FDA00003460406900022.GIF" wi="144" he="112" />为第i个免疫层析测试区显色特征区域的右边缘;将图像I的向x轴的投影设为E(x),即使关于x轴方向的height×3二维向量,其表达式如下:E(x)=(E<sup>1</sup>(x)',E<sup>2</sup>(x)',E<sup>3</sup>(x)')<maths num="0002"><![CDATA[<math><mrow><msup><mi>E</mi><mi>m</mi></msup><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>height</mi></mfrac><mo>&times;</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>height</mi></munderover><mrow><mo>(</mo><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>m</mi><mo>)</mo></mrow><mo>)</mo></mrow></mrow></math>]]></maths>其中m的取值为1、2、3分别代表着R、G、B维度,f(x,j,m)代表着在图像I中(x,j)点处的第m分量的颜色值;所述免疫层析测试区显色特征区域建模过程如下:第i个免疫层析测试区显色特征区域第m个颜色维度的左、右边缘非特征区采样点集合为:<maths num="0003"><![CDATA[<math><mrow><msubsup><mi>H</mi><mi>left</mi><mi>i</mi></msubsup><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><mo>{</mo><mrow><mo>(</mo><mi>x</mi><mo>,</mo><msup><mi>E</mi><mi>m</mi></msup><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>|</mo><mn>0</mn><mo>&lt;</mo><msubsup><mi>x</mi><mi>left</mi><mi>i</mi></msubsup><mo>-</mo><mi>&Delta;</mi><mo>&lt;</mo><mi>x</mi><mo>&lt;</mo><msubsup><mi>x</mi><mi>left</mi><mi>i</mi></msubsup><mo>}</mo></mrow></math>]]></maths><maths num="0004"><![CDATA[<math><mrow><msubsup><mi>H</mi><mi>right</mi><mi>i</mi></msubsup><mrow><mo>(</mo><mi>m</mi><mo>)</mo></mrow><mo>=</mo><mo>{</mo><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>E</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>|</mo><msubsup><mi>x</mi><mi>right</mi><mi>i</mi></msubsup><mo>&lt;</mo><mi>x</mi><mo>&lt;</mo><msubsup><mi>x</mi><mi>right</mi><mi>i</mi></msubsup><mo>+</mo><mi>&Delta;</mi><mo>&lt;</mo><mi>width</mi><mo>}</mo></mrow></math>]]></maths>Δ值决定拟合采样点的多少;拟合背景所选取的曲线形式,则需要根据图像采集光源的具体情况进行分析,如采用单光源,则图像的背景程抛物线形式分布,即第i个免疫层析测试区显色特征区域第m个颜色维度关于x轴的函数表达式为:c<sup>i</sup>(x,m)=k×(x-x0)<sup>2</sup>+a   <maths num="0005"><![CDATA[<math><mrow><mi>x</mi><mo>&Element;</mo><mrow><mo>(</mo><msubsup><mi>x</mi><mi>left</mi><mi>i</mi></msubsup><mo>-</mo><mi>&Delta;</mi><mo>,</mo><msubsup><mi>x</mi><mi>right</mi><mi>i</mi></msubsup><mo>+</mo><mi>&Delta;</mi><mo>)</mo></mrow></mrow></math>]]></maths>其中k、x0、a为抛物线的形状控制参数;假如采用平行光源,则可通过线性拟合便可进行背景建模,函数形式为:c<sup>i</sup>(x,m)=k×x+a在确定函数的形式后,由于集合<img file="FDA00003460406900032.GIF" wi="266" he="110" /><img file="FDA00003460406900033.GIF" wi="238" he="102" />均满足函数形式,即采用常用的最小二乘法或者是最大似然法,可得到相关位置参数,然后便拟合出免疫层析测试区显色特征区域的背景函数;所述免疫层析测试区显色特征区域中前景的真实颜色信息D获取如下:获取第i个免疫层析测试区显色特征区域x轴上第m个颜色维度的真实颜色信息D<sup>i</sup>(x,m):D<sup>i</sup>(x,m)=B<sup>i</sup>(x,m)-c<sup>i</sup>(x,m)D<sup>i</sup>(x,y,m)=f(x,y,m)-c<sup>i</sup>(x,m)D<sup>i</sup>(x)=(D<sup>i</sup>(x,1),D<sup>i</sup>(x,2),D<sup>i</sup>(x,3))<maths num="0006"><![CDATA[<math><mrow><msup><mi>D</mi><mi>i</mi></msup><mo>=</mo><mo>{</mo><msup><mi>D</mi><mi>i</mi></msup><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>|</mo><msubsup><mi>x</mi><mi>left</mi><mi>i</mi></msubsup><mo>&lt;</mo><mi>x</mi><mo>&lt;</mo><msubsup><mi>x</mi><mi>rightt</mi><mi>i</mi></msubsup><mo>}</mo></mrow></math>]]></maths>D={D<sup>i</sup>|i=1…n}其中,D为图像上所有免疫层析测试区显色特征区域真实的颜色矩阵;则总体的颜色信息熵Entropy<sup>i</sup>包括R、G、B分别得颜色信息熵Entropy<sub>R</sub>、Entropy<sub>G</sub>、Entropy<sub>B</sub>,其分别得表达式如下所示:<maths num="0007"><![CDATA[<math><mrow><msubsup><mi>Entropy</mi><mi>R</mi><mi>i</mi></msubsup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>x</mi><mo>=</mo><msubsup><mi>x</mi><mi>left</mi><mi>i</mi></msubsup></mrow><msubsup><mi>x</mi><mi>right</mi><mi>i</mi></msubsup></munderover><munderover><mi>&Sigma;</mi><mrow><mi>y</mi><mo>=</mo><mn>0</mn></mrow><mi>height</mi></munderover><msup><mi>D</mi><mi>i</mi></msup><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><msubsup><mi>Entropy</mi><mi>G</mi><mi>i</mi></msubsup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>x</mi><mo>=</mo><msubsup><mi>x</mi><mi>left</mi><mi>i</mi></msubsup></mrow><msubsup><mi>x</mi><mi>right</mi><mi>i</mi></msubsup></munderover><munderover><mi>&Sigma;</mi><mrow><mi>y</mi><mo>=</mo><mn>0</mn></mrow><mi>height</mi></munderover><msup><mi>D</mi><mi>i</mi></msup><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0009"><![CDATA[<math><mrow><msubsup><mi>Entropy</mi><mi>B</mi><mi>i</mi></msubsup><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>x</mi><mo>=</mo><msubsup><mi>x</mi><mi>left</mi><mi>i</mi></msubsup></mrow><msubsup><mi>x</mi><mi>right</mi><mi>i</mi></msubsup></munderover><munderover><mi>&Sigma;</mi><mrow><mi>y</mi><mo>=</mo><mn>0</mn></mrow><mi>height</mi></munderover><msup><mi>D</mi><mi>i</mi></msup><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0010"><![CDATA[<math><mrow><msup><mi>Entropy</mi><mi>i</mi></msup><mo>=</mo><mi>&lambda;</mi><mn>1</mn><mo>&times;</mo><msubsup><mi>Entropy</mi><mi>R</mi><mi>i</mi></msubsup><mo>+</mo><mi>&lambda;</mi><mn>2</mn><mo>&times;</mo><msubsup><mi>Entropy</mi><mi>G</mi><mi>i</mi></msubsup><mo>+</mo><mi>&lambda;</mi><mn>3</mn><mo>&times;</mo><msubsup><mi>Entropy</mi><mi>B</mi><mi>i</mi></msubsup></mrow></math>]]></maths>上式中λ1、λ2、λ3分别为对总体颜色信息熵的作用,选取λ1=-2、λ2=1、λ3=1。
地址 401120 重庆市渝北区紫荆路1号7幢2单元4-1