本发明涉及图像处理技术领域，尤其涉及一种固点噪声定位、精子显微环境洁净度评价方法。
背景技术：
精液检测是对男性生育能力评估的重要依据，计算机辅助精子质量分析系统因具有客观性、准确性等优点而成为精子质量分析的发展趋势。
在20倍精子显微镜观测野内，静态精子、杂质、灰尘、气泡、镜头污染物等在显微镜放大作用下，形成了位置固定、灰度较暗的颗粒影像，我们将这种位置固定的颗粒统称为固点噪声。其中，由静态精子形成的颗粒噪声称之为精子固点噪声；由杂质、灰尘、气泡及污物形成的颗粒噪声称为杂质固点噪声。杂质固点噪声的多少可在一定程度上可反映精液样本制作过程操作规范性以及显微镜光学系统洁净度。临床进行精子存活率计算时，首先必须进行静态精子的提取，而精子固点噪声是由静态精子形成的，所以静态精子提取的实质就是将精子固点噪声从杂质固点噪声中分离出来。显微环境的洁净度越低，两种固点噪声分离的难度越大，对应精子存活率指标的准确性越低，因此显微环境洁净度是精子质量分析系统的质量控制指标之一。
目前静态精子的提取大多基于最大类间差等相关分割算法进行，这种提取方法存在如下不足：
(1)精子固点噪声与杂质固点噪声无法有效分离，提取的精子固点噪声夹杂部分杂质固点噪声；
(2)对无明显双峰分布的视频图像，特别是静态精子稀疏的图像，最大类间差等传统分割算法根本无法将精子固点噪声从图像中分离出来；
(3)由于动态精子轨迹的灰度与部分杂质固点噪声相似，造成两者难以有效分离。
目前，由于缺乏有效的固点噪声提取算法，动态精子轨迹、精子固点噪声、杂质固点噪声无法进行有效分离，精子图像获取过程中显微环境洁度难以进行量化计算，致使目前计算机精子质量分析系统缺乏显微环境洁净度质控指标。
因此，研究有效的固点噪声定位算法，并在此基础上建立显微环境洁净度质控指标迫在眉睫。
技术实现要素：
本发明的目的就是为了解决上述问题，提供一种固点噪声定位、精子显微环境洁净度评价方法，能够有效分离静态精子固点噪声及杂质固点噪声，在此基础上给出了一种实用的精子显微环境洁净度评价方法。
为了实现上述目的，本发明采用如下技术方案：
一种固点噪声定位方法，包括以下步骤：
步骤一、对精子视频进行建模，获得含有固点噪声的背景图像；
步骤二、利用中智灰度和中智梯度两个特征的融合信息将固点噪声从背景像素中分离出来；
步骤三、通过调节噪声因子调节像素与理想背景像素的相似性，控制提取固点噪声的颗粒大小，将精子固点噪声和杂质固点噪声进行分离。
所述步骤一中，采用均值法对精子视频进行建模。
采用均值法对精子视频进行建模的具体建模公式为：
所述步骤二中，利用中智灰度进行分离时用到的具体函数为：
确定性
不确定性
错误性fg(x,y)＝1-tg(x,y)；
其中，g(x,y)表示背景图像中像素的灰度；gmin表示背景图像灰度最小值；gmax表示背景图像灰度最大值；gm为像素均值；dif(x,y)表示像素与图像均值gm的灰度绝对差；difmax和difmin表示dif(x,y)的最大值和最小值。
固点噪声在显微镜下的灰度越暗，像素灰度越低，对应的不确定性ig(x,y)越强，像素为噪声的可能性越大；如果像素灰度高于均值gm，则其为固点噪声的可能性小，为背景像素的可能性大，ig(x,y)降为0。
所述步骤二中，利用中智梯度进行分离时用到的具体函数为：
确定性：
不确定性：
错误性：fg(x,y)＝1-tg(x,y)；
其中，f(x,y)表示背景图像中值滤波后像素的灰度值；fmin表示滤波图像的灰度最小值；fmax表示滤波图像灰度最大值；i′g(x,y)为像素归一化前的不确定性；gx为像素在水平方向的梯度；gy为像素在垂直平方向的梯度；gmax为梯度最大值；gm为梯度均值；λ为噪声调节因子。
像素梯度越大，对应的中智不确定性越高，为噪声的可能性越大；通过改变i′g(x,y)中λ的大小，可改变像素梯度基准，进而改变噪声不确定性的大小，从而完成不同噪声颗粒的过滤筛选。
所述步骤二中，利用灰度与梯度融合信息将固点噪声从背景中提取出来的具体方法包括：
在中智域内，理想背景像素的灰度确定性tg(x,y)为1、不确定性ig(x,y)为0、错误性fg(x,y)为0，像素表示为：
在中智域内，理想背景像素的梯度确定性tg(x,y)为1、不确定性ig(x,y)为0、错误性fg(x,y)为0，像素表示为：
如果像素为噪声，则其与理想背景像素的相似性较低；因此像素与理想背景像素的相似性大小反映像素的噪声特性：相似性越低，为噪声的可能性越大；反之，则为背景的可能性越大；
像素与理想背景像素的灰度相似性sg(g(x,y),g*)表示为：
像素与理想背景像素的梯度相似性sg(g(x,y),g*)表示为：
基于灰度与梯度双特征融合的相似性snoise(x,y)表示为：
其中，w1和w2分别为灰度与梯度的权重系数，且
将snoise(x,y)与最大类间差法相结合，将固点噪声从背景中提取出来。
所述步骤三中还包括对筛选出的精子固点噪声与设定精子半径参考阈值进行比较，将小于设定精子半径参考阈值的杂质固点噪声剔除。
采用所述一种固点噪声定位方法的精子显微环境洁净度评价方法，包括：
获取杂质固点噪声的数量nns和静态精子的数量ns，则非精比rnso为：rnso＝nns/ns；
非精比rnso越大，说明视频序列显微图像中杂质或污物越多，表明精子样本本身、样本制作流程及显微图像获取各环节的洁净度越低。
本发明的有益效果：
基于噪声调节因子的固点噪声提取算法可将精子固点噪声和杂质固点噪声分别提取出来，能够计算非精比指标。
能够对显微环境洁净度进行衡量。
附图说明
图1(a)为实验用显微视频背景模型灰度图；
图1(b)为λ取1时对应的固点噪声；
图1(c)为λ取6时对应的固点噪声；
图1(d)为λ取10时对应的固点噪声；
图2为表1中标准实验序列的精子固点噪声与杂质固点噪声分布图；
图3为表1中130号序列的精子固点噪声与杂质固点噪声分布图。
具体实施方式
下面结合附图与实施例对本发明作进一步说明。
首先建立具有噪声调节功能的固点噪声提取算法；然后进行杂质固点噪声和精子固点噪声的分离；最后计算精子与杂质固点噪声数量之比-非精比，并将此参数作为显微环境洁净度评价指标。具体实施步骤如下：
(1)背景建模：采用均值法对精子视频进行建模，获得含有固点噪声的背景图像gb。均值建模公式如下：
其中：gi(x,y)为视频序列的第i图像，n是视频序列图像的帧数，取值25。
(2)固点噪声定位算法的建立：本发明提出具有噪声调节功能的固点噪声定位算法，该算法利用中智灰度和中智梯度两个特征的融合信息将固点噪声从背景像素中分离出来。
①在中智域内，像素基于灰度特征的确定性tg(x,y)、不确定性ig(x,y)、错误性fg(x,y)量化函数如下：
fg(x,y)＝1-tg(x,y)
其中，g(x,y)表示背景图像中像素的灰度；gmin表示背景图像灰度最小值；gmax表示背景图像灰度最大值；gm为像素均值；dif(x,y)表示像素与图像均值gm的灰度绝对差。difmax和difmin表示dif(x,y)的最大值和最小值。
固点噪声在显微镜下的灰度较暗，所以像素灰度越低，对应的不确定性ig(x,y)越强，说明像素为噪声的可能性越大。如果像素灰度高于均值gm，则其为固点噪声的可能性小，为背景像素的可能性越大，其ig(x,y)降为0。
②在中智域内，利用梯度特征对背景像素噪声特性进行衡量，噪声确定性tg(x,y)、不确定性ig(x,y)、错误性fg(x,y)量化函数如下：
fg(x,y)＝1-tg(x,y)
其中，f(x,y)表示背景图像中值滤波后像素的灰度值；fmin表示滤波图像的灰度最小值；fmax表示滤波图像灰度最大值；i′g(x,y)为像素归一化前的不确定性；gx为像素在水平方向的梯度；gy为像素在垂直平方向的梯度；gmax为梯度最大值；gm为梯度均值；λ为噪声调节因子。
如果一个像素为背景像素，则其灰度变化较小，只有灰度变化较大的地方才有噪声出现的可能。所以像素梯度越大，对应的不确定性越高，为噪声的可能性越大。通过改变i′g(x,y)中λ的大小，可以改变像素梯度基准，进而改变噪声不确定性的大小。所以，利用λ这个噪声调节因子，可调节像素的噪声特性，进而完成不同噪声颗粒的过滤筛选。
③基于双特征融合的相似性度量：
在中智域内，理想背景像素的灰度确定性tg(x,y)为1、不确定性ig(x,y)为0、错误性fg(x,y)为0，该像素可表示为：
在中智域内，理想背景像素的梯度确定性tg(x,y)为1、不确定性ig(x,y)为0、错误性fg(x,y)为0，该像素可表示为：
如果像素为噪声，则其与理想背景像素的相似性较低。因此像素与理想背景像素的相似性大小可以反映像素的噪声特性：相似性越低，为噪声的可能性越大；反之，则为背景的可能性越大。
像素与理想背景像素的灰度相似性sg(g(x,y),g*)可用下式衡量：
像素与理想背景像素的梯度相似性sg(g(x,y),g*)可用下式计算：
基于灰度与梯度双特征融合的相似性snoise(x,y)用下式衡量：
将snoise(x,y)与最大类间差法相结合，便可将固点噪声从背景中提取出来。
(3)精子与杂质的分离：通过调节噪声因子λ，可调节像素与理想背景像素的相似性，进而可控制提取固点噪声的颗粒大小，从而将精子固点噪声和杂质固点噪声进行分离。精子样本千差万别，但人类精子的大小在一定范围内变化。一般情况下显微镜下精子细胞影像比杂质颗粒、灰尘及气泡等的影像大。对筛选出的精子固点噪声与设定精子半径参考阈值比较，可进一步剔除精子固点噪声中残存的杂质固点噪声，进而较准确地将静态精子提取出来。
(4)显微环境洁净度量化参数非精比的计算：nns、ns分别为杂质固点噪声和静态精子的数量，则非精比rnso定义为：rnso＝nns/ns。
(5)显微环境洁净度评价结论：非精比rnso越大，说明视频序列显微图像中杂质或污物越多，表明精子样本本身、样本制作流程及显微图像获取各环节的洁净度越低。
具体的：
第一步，计算固点噪声背景模型，实验用显微视频背景模型灰度图如图1(a)所示。图中已经将精子、杂质、动态精子轨迹标记出来。
第二步，改变噪声调节因子λ的大小，λ取不同值时对应的固点噪声如图1(b)-图1(d)所示。λ＝1时，提取的固点噪声像素最多，但动态精子运动轨迹被提取出来；当λ＝6时，动态精子轨迹基本被过滤，剩余的为精子与杂质的混合颗粒；λ＝10时，杂质等较小颗粒基本被过滤，剩下的主要为静态精子等颗粒较大的噪声。按精子大小半径进一步过过滤，便可将精子提取出来。从精子与杂质的混合物中将精子检出，剩下的即为杂质固点噪声。
第三步，计算显微图像的非精比指标。7个实验用精子显微图像的非精比rnso如表1所示。以标准实验样本和130号序列的图像为例来说明rnso越大，显微图像中杂质固点噪声越多。
标准实验样本和130号序列中精子与杂质固点噪声分布如图2-图3所示。图2-图3中方形区域表示精子，圆形区域表示杂质固点噪声。由图2-图3可见，标准序列中杂质固点噪声远少于130号序列，而表1中标准序列的rnso为1.11，130号序列的rnso为4.37，实验证明rnso指标可以反映显微环境的洁净度。
从精子提取效率与显微环境洁净度方面说明本发明的优点。
(1)精子提取效率之比：为比较本文固点噪声定位算法对静态精子的提取效率，从错检率和漏检率两方面与最大熵法和最大类间差法的分割结果进行比较，并用错检率和漏检率之和famr作为算法对静态精子提取效率的衡量指标。三种算法的famr比较结果如表2所示。famr越小，说明算法提取静态精子的效率越高。表2中“×”号表示精子不可检出的情况，说明最大类间差法无法提取精子的情况比较严重。由表2数据可见本发明提出的固点噪声提取算法具有较低的famr指标。
(2)显微环境洁净度衡量：基于噪声调节因子的固点噪声提取算法可将精子和杂质固点噪声分别提取出来，能够计算非精比指标，而目前基于传统分割算法的精子提取算法尚缺乏对显微环境洁净度量化的研究。
图1(b)中，噪声调节因子λ为1，动态精子运动轨迹及固点噪声混杂在一起；图1(c)中，λ为6，动态精子轨迹基本被过滤，剩下精子固点噪声与杂质固点噪声的混合颗粒；图1(d)中，λ为10，杂质等较小颗粒基本被过滤，剩下的主要为颗粒较大的精子固点噪声。说明通过调节λ的大小，可以将杂质固点噪声和精子固点噪声分别提取出来。
表1实验用显微序列非精比指标
表1给出7个实验用显微序列的非精比指标。rnso越大，表明显微图像含杂质固点噪声的数量越多。
图2为表1中标准实验序列的精子与杂质固点噪声分布图，图3为表1中130号序列的精子与杂质固点噪声分布图。其中，圆形区域数量越多说明杂质固点噪声越多，显微环境洁净度越差，对应的rnso指标越大。
表2三种算法的famr指标
注：表中“×”号表示精子不可检出。
表2中“×”号表示精子不可检出的情况，说明最大类间差法无法提取精子的情况比较严重。famr越小，说明算法提取精子的效率越高。
上述虽然结合附图对本发明的具体实施方式进行了描述，但并非对本发明保护范围的限制，所属领域技术人员应该明白，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。