[发明专利]一种卵裂阶段的细胞计数定位方法

权利要求书

1.一种卵裂阶段的细胞计数定位方法，其特征在于，包括如下步骤：

S1：获取卵裂阶段的细胞显微放大图像；

S2：将显微放大图像转换为灰度图像，并进行灰度变换；

S3：利用Otsu算法寻找灰度变换后的图像的最佳阈值；

S4：利用求得的最佳阈值将灰度变换后的图像二值化；

S5：对二值化图像进行形态学变换；

S6：对形态学变换后的图像进行轮廓检测；

S7：对检测到的轮廓进行多边形逼近；

S8：检测质心点，并计算轮廓上每一点到质心的距离；

S9：比较轮廓上各点到质心距离的变化情况，将连续增加到连续减小的拐点标记为远心点，将连续减小到连续增加的拐点标记为近心点，将近心点中距离最远的点标记为针尖，远心点个数即为细胞凸包数，通过细胞凸包数即可判定卵裂细胞处于第几个发育阶段。

2.根据权利要求1所述的卵裂阶段的细胞计数定位方法，其特征在于：所述步骤S2中的灰度变换为灰度线性拉伸。

3.根据权利要求1所述的卵裂阶段的细胞计数定位方法，其特征在于：所述步骤S3中Otsu算法是利用类间方差作为判据，选取使类间方差最大的图像灰度值作为最佳阈值。

4.根据权利要求1所述的卵裂阶段的细胞计数定位方法，其特征在于：所述步骤S4中二值化变换函数为：

其中，T为最佳阈值。

5.根据权利要求1所述的卵裂阶段的细胞计数定位方法，其特征在于：所述步骤S5中形态学变换为腐蚀和膨胀，使用的结构元素为3×3。

6.根据权利要求1所述的卵裂阶段的细胞计数定位方法，其特征在于：所述步骤S7中多边形逼近的具体方法为：设原来的数字曲线为P＝{p1，p2，p3，p4，…，pn}；指定的逼近误差门限阈值为ε，逼近多边形的顶点上界为μ，逼近后多边形为P’＝{p1’，p2’，p3’，p4’，···，pm’}，则逼近后的多边形应该满足以下条件：

其中Value(P)表示多边形的边长。

7.根据权利要求1所述的卵裂阶段的细胞计数定位方法，其特征在于：所述步骤S8中质心检测的原理如下：不规则区域的矩，表示把一个归一化的灰度级图像函数理解为一个二维随机变量的概率密度，这个随机变量的属性可以用统计特征--矩来描述，通过假设非零的像素值表示区域，矩可以用于二值或灰度级的区域描述：其中x，y，i，j是区域点的坐标,具体是在数字图像中的像素坐标，f(i，j)是图像区域的灰度值，令Xc，Yc表示区域重心的坐标，则：

其中，在二值图像的情况下，M₀₀表示区域的面积。