发明名称 基于无透镜全息成像的淡水藻类粗分类与计数方法
摘要 本发明公开了基于无透镜全息成像的淡水藻类细胞粗分类与计数方法,首先,对利用无透镜全息成像装置获取的淡水藻类全息图像进行二值化处理,并计算所有细胞面积占整个图像面积的比值;然后,利用基于细胞形状特征的方法对细胞分散、数目不多的细胞图像进行简单分类和计数;最后,利用基于归一化互相关匹配方法对细胞集中、数目多的细胞图像进行简单分类和计数。本发明能够更为针对性、准确性地对淡水藻类全息图像中的细胞进行粗分类与计数,而且针对不同细胞分布、不同细胞数目的细胞图像,包含不同的解决方案,从而能实现更好的粗分类与计数效果。同时实现的复杂度较低,图像处理时间短,最终的粗分类与计数的结果也较为理想。
申请公布号 CN106022303A 申请公布日期 2016.10.12
申请号 CN201610392278.7 申请日期 2016.06.06
申请人 南昌航空大学 发明人 陈震;冷健雄;张聪炫;张初华;王官权;江少锋
分类号 G06K9/00(2006.01)I;G06K9/62(2006.01)I;G06T7/00(2006.01)I 主分类号 G06K9/00(2006.01)I
代理机构 南昌市平凡知识产权代理事务所 36122 代理人 张文杰
主权项 基于无透镜全息成像的淡水藻类细胞粗分类与计数方法,其特征在于,包括如下步骤:1)对利用无透镜全息成像装置获取的淡水藻类全息图像进行二值化处理,并计算所有细胞面积占整个图像面积的比值S,当S&lt;=0.5时,转至步骤2);否则,转至步骤3);2)利用基于细胞形状特征的方法对细胞分散、数目不多的细胞图像进行简单分类和计数;步骤如下:a)计算二值化处理后的图像中每个细胞的长、宽、面积和外接矩形的大小数据,并据此计算圆形度、矩形度的形状特征和参数的大小:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>C</mi><mo>=</mo><mfrac><msubsup><mi>P</mi><mi>n</mi><mn>2</mn></msubsup><mrow><mn>4</mn><msub><mi>&pi;A</mi><mi>n</mi></msub></mrow></mfrac><mo>;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001010351510000011.GIF" wi="1086" he="160" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>R</mi><mo>=</mo><mfrac><msub><mi>A</mi><mi>n</mi></msub><msub><mi>A</mi><mi>r</mi></msub></mfrac></mrow>]]></math><img file="FDA0001010351510000012.GIF" wi="137" he="122" /></maths>式中:C表示圆形度;P<sub>n</sub>和A<sub>n</sub>分别为细胞的周长和面积;R表示矩形度;A<sub>n</sub>为面积;A<sub>r</sub>为图形最小外接矩形的面积;b)根据图像实际情况设定不同的分类阈值,用不同颜色对细胞进行分类标记,并得出每一类标记细胞的数目;3)利用基于归一化互相关匹配方法,对细胞集中、数目多的细胞图像进行简单分类和计数;步骤如下:a)创建细胞模版图像库,单个模版图像的大小为9×9像素,且细胞中心位于图像正中间,数目为20;b)计算每个细胞模板图像与待测试图像的所有匹配区域的二维归一化互相关系数:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>R</mi><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>y</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>L</mi><mo>=</mo><mn>1</mn></mrow></msubsup><msubsup><mi>&Sigma;</mi><mrow><mi>x</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>K</mi><mo>=</mo><mn>1</mn></mrow></msubsup><mrow><mo>(</mo><mi>w</mi><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mo>)</mo><mo>-</mo><mover><mi>w</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mrow><mo>(</mo><mi>f</mi><mo>(</mo><mrow><mi>x</mi><mo>+</mo><mi>i</mi><mo>,</mo><mi>y</mi><mo>+</mo><mi>j</mi></mrow><mo>)</mo><mo>-</mo><mover><mi>f</mi><mo>&OverBar;</mo></mover><mo>(</mo><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mo>)</mo><mo>)</mo></mrow></mrow><mrow><msup><mrow><mo>&lsqb;</mo><msubsup><mi>&Sigma;</mi><mrow><mi>y</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>L</mi><mo>=</mo><mn>1</mn></mrow></msubsup><msubsup><mi>&Sigma;</mi><mrow><mi>x</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>K</mi><mo>=</mo><mn>1</mn></mrow></msubsup><msup><mrow><mo>(</mo><mi>w</mi><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mo>)</mo><mo>-</mo><mover><mi>w</mi><mo>&OverBar;</mo></mover><mo>)</mo></mrow><mn>2</mn></msup><mo>&rsqb;</mo></mrow><mfrac><mn>1</mn><mn>2</mn></mfrac></msup><msup><mrow><mo>&lsqb;</mo><msubsup><mi>&Sigma;</mi><mrow><mi>y</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>L</mi><mo>=</mo><mn>1</mn></mrow></msubsup><msubsup><mi>&Sigma;</mi><mrow><mi>x</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>K</mi><mo>=</mo><mn>1</mn></mrow></msubsup><msup><mrow><mo>(</mo><mi>f</mi><mo>(</mo><mrow><mi>x</mi><mo>+</mo><mi>i</mi><mo>,</mo><mi>y</mi><mo>+</mo><mi>j</mi></mrow><mo>)</mo><mo>-</mo><mover><mi>f</mi><mo>&OverBar;</mo></mover><mo>(</mo><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow><mo>)</mo><mo>)</mo></mrow><mn>2</mn></msup><mo>&rsqb;</mo></mrow><mfrac><mn>1</mn><mn>2</mn></mfrac></msup></mrow></mfrac><mo>;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001010351510000013.GIF" wi="1885" he="231" /></maths>式中:<img file="FDA0001010351510000014.GIF" wi="44" he="55" />和<img file="FDA0001010351510000015.GIF" wi="146" he="78" />是图像w和f的平均灰度;R是归一化互相关系数;c)设定粗分类识别阈值为0.65,并寻找超过设定阈值的系数定位,且细胞位置并用不同颜色标记和计数。
地址 330063 江西省南昌市丰和南大道696号