一种用于小球藻自动计数的图像调色板阈值分割方法技术

本发明专利技术公开了一种用于小球藻自动计数的图像调色板阈值分割方法，以小球藻显微荧光图像各颜色占调色板总颜色的概率的累积分布函数作为变换函数，其中累积分布函数表示为：式中，rj为原图像归一化调色板颜色所占比重，nj为原图像中出现第k级颜色的像素个数，n为图像调色板颜色的总数，Sk为HE（直方图均衡）后图像归一化后的颜色比重，T(rk)为变换关系式，pr(rj)为原图像取第k级颜色的概率，即原图像的调色板直方图，处理后得到的第一图像的调色板分布为k=0，1，.....，N一1，则变换关系为r′k＝round((N-1)sk)，r′k为均衡后所显示的颜色，统计颜色在调色板中的直方图分布情况，找到小球藻颜色在调色板中的阈值分割。

【技术实现步骤摘要】

本专利技术属于水产养殖
，特别涉及。
技术介绍
小球藻，俗称绿藻，作为一种具有保健食品作用的藻类，人工培养的生产已经由来已久。在人工培养中，对于小球藻个数的计数是一项重要的环节。但是，在现有技术中，藻类的个数的计数工作主要依靠人工，劳动强度大，观测时间长，容易引起疲劳误差。目前，还没有一种有效的方法解决上述问题。·
技术实现思路
本专利技术的目的是提供，用于对于小球藻自动计数，该方法利用了藻类荧光激发效应，藻类受一定波长的光波照射后发射出特定范围波长的光的一种现象。本专利技术的技术方案是，，包括以下步骤将图像格式为24位真彩色图像的小球藻显微荧光图像作为原图像，以小球藻显微荧光图像各颜色占调色板总颜色的概率的累积分布函数作为变换函数，将原图像映射为一幅调色板颜色级分布较均匀的图像，其中累积分布函数表示为 KKSk =T(Ijt)=〉」=〉Pr Cfj) j=0J=OO ^ Tj ^ I ；k = O, I, . . . , I (I)式中，r」为原图像归一化调色板颜色所占比重，η」为原图像中出现第k级颜色的像素个数，η为图像调色板颜色的总数，Sk为HE (直方图均衡)后图像归一化后的颜色比重，T(rk)为变换关系式，pJh )为原图像取第k级颜色的概率，即原图像的调色板直方图，假设处理后得到的第一图像的调色板分布为k=0，1，……，N— 1，则变换关系为，T1 k = round((N-1) sk)(2)k=0，l，....，N—l(rOUnd为四舍五人取整运算)，r' k为均衡后所显示的颜色，利用经过调色板均衡化后的第一图像，统计第一图像颜色在调色板中的直方图分布情况，找到小球藻颜色在调色板中的阈值，选取离此阈值正负相差各5个调色板距离的两个点作为区域分割阈值，对第一图像进行分割。本专利技术基于256色图像调色板的阈值分割算法，利用经过调色板均衡化后的图像，统计图像颜色在调色板中的直方图分布情况，找到目标物体颜色在调色板中的阈值，由于原始图像中目标与背景颜色非常接近，调色板均衡化后目标与背景各颜色间的在调色板中距离变大，比较目标藻类在图像中颜色，以此颜色在调色板中的位置设定上下分割阈值，附图说明图1是本专利技术的小球藻自动计数方法流程图。图2是本专利技术实施例中256色小球藻突光图像。图3是本专利技术实施例中调色板直方图均衡后图像。图4是本专利技术实施例中颜色阈值分割后图像。具体实施例方式对于小球藻自动计数实验选用实验室培养的小球藻实验液，小球藻密度为500 600万个/ml。将小球藻培养液压片后,至于奥林巴斯BX51突光显微镜下观察,显微镜目镜放大倍数为10倍，选择物镜放大倍数为40倍。实验利用显微镜自带的DP71彩色显微摄像机拍摄藻类荧光激发图像，相机CCD大小为2/3英寸，1250万有效像素，采集图像尺寸为1360X1024,以tif格式存储。试验中分别利用450 480nm波长的光束和510 550nm波长的光束照射小球藻样品，发现在510 550nm波长的光束照射下小球藻具有明显的红色荧光，其激发光波长处于590nm附近。实验利用红色滤光片滤除其它干扰光线，得到小球藻荧光激发图片。对于小球藻图像处理流程如图1所示，将得到小球藻藻类荧光激发图像，转换成256色图像如图2所示。以提高藻类计数精度为目标，提出一种基于调色板直方图均衡化彩色图像分割方法。在RGB颜色空间中以256色图像调色板中颜色分布坐标为外部激励输入，以相同颜色像素占总图像总像素的比重为均衡化密度，用改进灰度直方图均衡化算法对藻类图像进行增强(图3为增强后图像)，再根据增强后图像颜色直方图特点选取颜色阈值进行图像分割(图4为分割后图像)。对于小球藻荧光激发图像进行藻类计数，由于上述的处理过程能够很好的滤除液体中杂质的影响，因此图像分割后我们认为图像中黑像素点，不论关联像素面积大小皆可认为是一个小球藻个体。按照BMP图像格式从左到右，从下到上扫描图像，如果遇到找到一个黑像素点，假设为A点，把A点做种子点.将其标记值设为1.并向外寻找与其八邻域相连的其他黑色像素点.将找到的每个这样点的标记值都设为I。再分别以每个这样的点为种子点继续寻找.直到找不到未标记过的相连的黑色像索点为止.这样一个连通区域就标记完毕。继续从A点处按照扫描顺序寻找下一个未标记过的黑色像素点做种子点.将其标记值设为2.寻找与其八邻域相连的未标记过的黑色像索点，如此反复，直到整幅图像扫描完毕，最大的区域标记值即为所计的藻类个数。所述的图像处理算法包括灰度均衡算法，该算法的原理如下。小球藻显微荧光图像图像格式为24位真彩色图像，针对24位真彩色图像的处理程序，现在还不能将小球藻个体与背景分割开来，因此必须要将图像进行变换。传统的直方图均衡化方法是以图像各灰度级概率的累积分布函数作为变换函数，将原图像映射为一幅灰度级分布较均匀的图像。其中累积分布函数表示为权利要求1.，其特征在于，包括以下步骤将图像格式为24位真彩色图像的小球藻显微荧光图像作为原图像，以小球藻显微荧光图像各颜色占调色板总颜色的概率的累积分布函数作为变换函数，将原图像映射为一幅调色板颜色级分布较均匀的图像，其中累积分布函数表示为.O ^ Tj ^ I ；k = O, I, . . . , I (I)式中，r」为原图像归一化调色板颜色所占比重，η」为原图像中出现第k级颜色的像素个数，η为图像调色板颜色的总数，Sk为直方图均衡后图像归一化后的颜色比重，T(rk)为变换关系式，Pr (rj)为原图像取第k级颜色的概率，即原图像的调色板直方图，假设处理后得到的第一图像的调色板分布为k=0，1，.....，N — 1，则变换关系为， r' k = round((N-1) sk)(2)k=0,1，. . . .，N — 1，round为四舍五人取整运算，r/ k为均衡后所显示的颜色，利用经过调色板均衡化后的第一图像，统计第一图像颜色在调色板中的直方图分布情况，找到小球藻颜色在调色板中的阈值，选取离此阈值正负相差各5个调色板距离的两个点作为区域分割阈值，对第一图像进行分割。全文摘要本专利技术公开了，以小球藻显微荧光图像各颜色占调色板总颜色的概率的累积分布函数作为变换函数，其中累积分布函数表示为式中，rj为原图像归一化调色板颜色所占比重，nj为原图像中出现第k级颜色的像素个数，n为图像调色板颜色的总数，Sk为HE(直方图均衡)后图像归一化后的颜色比重，T(rk)为变换关系式，pr(rj)为原图像取第k级颜色的概率，即原图像的调色板直方图，处理后得到的第一图像的调色板分布为k=0，1，.....，N一1，则变换关系为r′k＝round((N-1)sk)，r′k为均衡后所显示的颜色，统计颜色在调色板中的直方图分布情况，找到小球藻颜色在调色板中的阈值分割。文档编号G06T7/00GK102999907SQ20121046698公开日2013年3月27日 申请日期2012年11月17日 优先权日2012年11月17日专利技术者刘世晶, 汤涛林, 苗雷 申请人:中国水产科学研究院渔业机械仪器研究所本文档来自技高网...

【技术保护点】
一种用于小球藻自动计数的图像调色板阈值分割方法，其特征在于，包括以下步骤：将图像格式为24位真彩色图像的小球藻显微荧光图像作为原图像，以小球藻显微荧光图像各颜色占调色板总颜色的概率的累积分布函数作为变换函数，将原图像映射为一幅调色板颜色级分布较均匀的图像，其中累积分布函数表示为：SK=T(rk)=Σj=0Knjn=Σj=0Kpr(rj)0≤rj≤1；k＝0，1，...，I????（1）式中，rj为原图像归一化调色板颜色所占比重，nj为原图像中出现第k级颜色的像素个数，n为图像调色板颜色的总数，Sk为直方图均衡后图像归一化后的颜色比重，T(rk)为变换关系式，pr(rj)为原图像取第k级颜色的概率，即原图像的调色板直方图，假设处理后得到的第一图像的调色板分布为k=0，1，.....，N一1，则变换关系为,r′k＝round((N?1)sk)????????（2）k=0，1，....，N一1，round为四舍五人取整运算，r′k为均衡后所显示的颜色，利用经过调色板均衡化后的第一图像，统计第一图像颜色在调色板中的直方图分布情况，找到小球藻颜色在调色板中的阈值，选取离此阈值正负相差各5个调色板距离的两个点作为区域分割阈值，对第一图像进行分割。...

【技术特征摘要】

【专利技术属性】
技术研发人员：刘世晶，汤涛林，苗雷，
申请(专利权)人：中国水产科学研究院渔业机械仪器研究所，
类型：发明
国别省市：