[发明专利]一种血液肿瘤染色体核型自动分析方法及系统

权利要求书

1.一种血液肿瘤染色体核型自动分析方法，其特征在于，包括以下步骤：

S1：采集染色体图像；

S2：对采集的染色体图像进行预处理；

S3：根据预处理得到的染色体图像对染色体进行分割；

S4：对分割完成的染色体进行特征参数的提取；

S5：根据特征参数的提取结果，对染色体进行分类，完成染色体核型自动分析。

2.根据权利要求1所述的一种血液肿瘤染色体核型自动分析方法，其特征在于，在所述步骤S2中，所述预处理包括灰度变换、图像滤波、图像分割、边缘检测与边界跟踪；其中：

所述灰度变换采用线性单值函数公式进行变换，具体表示为：

由于少部分像素分布在小于a和大于b的区间内,此时可用下式作变换：

其中，原染色体图像F(j,k)的灰度范围为[a,b]，线性变换后图像F'(j,k)的范围为[a',b']；

图像滤波采用中值滤波法：对窗口内的所有像素的灰度值进行排序，取排序结果的中间值作为原窗口中心点处的像素的灰度值；

图像分割采用灰度阈值分割法，并应用迭代法选取阈值；

边缘检测采用Sobel算子算法进行处理；

所述边界跟踪采用Freeman链码跟踪算法处理。

3.根据权利要求2所述的一种血液肿瘤染色体核型自动分析方法，其特征在于，在所述步骤S3中，所述对染色体进行分割具体包括以下步骤：

S31：采用三点链码计算方法寻找染色体边界凹点以作为候选分割点，具体为：

第一步，求相对链码RelativeCode[i]，即出该点的链码值减去进入该点的链码值；

第二步，求绝对链码AbsoluteCode[i]，即将前一点的绝对链码加上当前点的相对链码；

第三步，求相对链码和SumCode[i]，即将当前点与相邻前两点的绝对链码和相加；

第四步，求链码差DiffCode[i]，即将当前点的后第三个点的链码和减去当前点的链码和；

当链码差小于0时，表示该点为凹点，且负值越小越凹陷，当链码差大于0时，表示该点为凸点，且正值越大越凸出；

S32：通过链码差的阈值判断寻找合适的分割点，具体计算过程为：

D_k＝{i|i＝0，1，2，…，N}

C_k＝{j|DiffCode[j]T_C，j∈D_k}

其中，D_k表示该条染色体轮廓上每个点的序号集合，C_k为该条染色体所有满足阈值条件的点集，T_c为负值；

S33：连接分割点完成分割，需满足下列阈值条件：

S＝{(x,y)|L(x,y)＜T_d∧Arc(x，y)T_a}

其中，S为最终确立的分割点集，L(·)计算两点间的距离，T_d为分割路径阈值，Arc(·)计算两点间的弧长，T_a为染色体边界弧长阈值；连接S内每个二元元素所表示的点，完成粘连分割。