CN114897843A - Overlapping chromosome segmentation method - Google Patents

Overlapping chromosome segmentation method Download PDF

Info

Publication number
CN114897843A
CN114897843A CN202210548618.6A CN202210548618A CN114897843A CN 114897843 A CN114897843 A CN 114897843A CN 202210548618 A CN202210548618 A CN 202210548618A CN 114897843 A CN114897843 A CN 114897843A
Authority
CN
China
Prior art keywords
points
chromosome
skeleton
point
cutting
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202210548618.6A
Other languages
Chinese (zh)
Inventor
崔霄
季肖辉
张志锋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhengzhou University of Light Industry
Original Assignee
Zhengzhou University of Light Industry
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhengzhou University of Light Industry filed Critical Zhengzhou University of Light Industry
Priority to CN202210548618.6A priority Critical patent/CN114897843A/en
Publication of CN114897843A publication Critical patent/CN114897843A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/12Edge-based segmentation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/13Edge detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10056Microscopic image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20084Artificial neural networks [ANN]

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Quality & Reliability (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

The invention provides a method for segmenting overlapping chromosomes, which comprises the following steps of S100, carrying out edge extraction on an overlapping chromosome image, obtaining edge outlines of all overlapping chromosomes, and adding all edge outline points into a set of cutting points to be selected; step S200, extracting a skeleton, a skeleton end point and a skeleton intersection point of the overlapped chromosome image; step S300, determining the cutting points of the overlapped chromosomes according to the edge contour of the overlapped chromosomes obtained in the step S100 and the chromosome skeleton, the skeleton end points and the intersection points obtained in the step S200; and S400, automatically segmenting the overlapped chromosomes according to the cutting points. According to the contour and skeleton-based overlapped chromosome segmentation method, all contour points are regarded as cut points and used as candidate sets by combining the relation between the contour and the skeleton of the chromosome by utilizing the overlapped features of the overlapped chromosomes, a series of complex processes of curvature calculation are omitted, and the optimal cut points are accurately positioned by combining Euclidean distance with the skeleton of the chromosome.

Description

一种重叠染色体分割方法A method for segmenting overlapping chromosomes

技术领域technical field

本发明涉及染色体核型图像分析技术领域,尤其是涉及一种重叠染色体分割方法。The invention relates to the technical field of chromosome karyotype image analysis, in particular to a method for segmenting overlapping chromosomes.

背景技术Background technique

染色体是生命遗传物质的载体,以丝状或棒状形态存在于细胞核中,在细胞有丝分裂中期,可通过染色制片和显微成像,观察其数目和形态结构特征,染色体异常会导致遗传性疾病、新生儿缺陷等病症,染色体核型分析是寻找染色体表型数目和结构异常,如数目不符、片段异常(重复、缺失、易位、倒位等)、形态异常(染色体长度或随体大小与正常染色体不符)等异常情况,医生通过染色体核型分析,对患者做出精确的诊疗。Chromosomes are the carriers of the genetic material of life. They exist in the nucleus in the form of filaments or rods. In the middle stage of cell mitosis, the number and morphological characteristics of chromosomes can be observed by staining and microscopic imaging. Chromosomal abnormalities can lead to genetic diseases, For diseases such as neonatal defects, karyotype analysis is to look for abnormal chromosome phenotype number and structure, such as number inconsistency, fragment abnormalities (duplication, deletion, translocation, inversion, etc.), morphological abnormalities (chromosomal length or satellite size and normal Chromosome inconsistency) and other abnormal conditions, doctors make accurate diagnosis and treatment of patients through chromosomal karyotype analysis.

染色体核型图像中染色体具有不可预测性,表现为易粘连、边缘模糊、重叠嵌套多等特点,现有染色体核型图像处理算法对染色体分割、计数、消冗效果较差,精度有待提高,同时,染色体核型分析过程中需要人工全程参与,自动化程度低,染色体核型分析工作效率低下、劳动强度大、专业人才培养周期长,通过利用数字图像处理技术解决染色体计数、分割、分类等问题,可以提高工作人员的分析效率、降低劳动强度、减少专业人才的培养周期。Chromosomes in karyotype images are unpredictable, characterized by easy adhesion, blurred edges, and many overlapping and nesting. The existing karyotype image processing algorithms have poor effects on chromosome segmentation, counting, and redundancy elimination, and the accuracy needs to be improved. At the same time, the whole process of chromosome karyotype analysis requires manual participation in the whole process, the degree of automation is low, the work efficiency of chromosome karyotype analysis is low, labor-intensive, and the training cycle of professionals is long. The use of digital image processing technology solves problems such as chromosome counting, segmentation, and classification. , which can improve the analysis efficiency of staff, reduce labor intensity, and reduce the training cycle of professional talents.

发明内容SUMMARY OF THE INVENTION

为了克服背景技术中的不足,本发明公开了一种重叠染色体分割方法,本发明有效提高染色体核型检测分析的准确性和高效性。In order to overcome the deficiencies in the background technology, the present invention discloses a method for segmenting overlapping chromosomes, which effectively improves the accuracy and efficiency of chromosome karyotype detection and analysis.

为了实现所述发明目的,本发明采用如下技术方案:In order to realize the purpose of the invention, the present invention adopts the following technical solutions:

一种重叠染色体分割方法,所述的分割方法具体包括以下步骤:A method for segmentation of overlapping chromosomes, the segmentation method specifically comprises the following steps:

步骤S100,对重叠染色体图像进行边缘提取,并获得所有重叠染色体的边缘轮廓,将所述的边缘轮廓点全部加入到待选切割点集合;Step S100, perform edge extraction on the overlapping chromosome images, and obtain the edge contours of all overlapping chromosomes, and add all the edge contour points to the set of cutting points to be selected;

步骤S200,对重叠染色体图像提取其骨架、骨架端点和骨架交叉点;Step S200, extracting its skeleton, skeleton endpoints and skeleton intersections from the overlapping chromosome images;

步骤S300,根据步骤S100获取的重叠染色体边缘轮廓和步骤S200获取的染色体骨架及骨架端点、交叉点确定重叠染色体的切割点;Step S300, determining the cutting point of the overlapping chromosome according to the edge contour of the overlapping chromosome obtained in step S100 and the chromosome skeleton, skeleton endpoints, and intersection points obtained in step S200;

步骤S400,根据步骤S300所述的切割点通过自动切割方法对所述重叠染色体进行自动分割。Step S400, the overlapping chromosomes are automatically segmented by an automatic cutting method according to the cutting point described in step S300.

所述步骤S200具体包括以下步骤:The step S200 specifically includes the following steps:

步骤S201,采用细化算法将重叠染色体图像细化为仅有一条像素点的染色体骨架,并提取该重叠染色体骨架;Step S201, using a thinning algorithm to refine the overlapping chromosome image into a chromosome skeleton with only one pixel, and extract the overlapping chromosome skeleton;

步骤S202,根据邻近像素点数量定位骨架端点和交叉点;Step S202, locate skeleton endpoints and intersections according to the number of adjacent pixels;

所述步骤S300具体包括以下步骤:The step S300 specifically includes the following steps:

S301,根据步骤S200所述的染色体骨架和骨架交叉点,以骨架交叉点为中心,计算待选切割点到中心点的欧式距离,为此距离进行设定为,去除待选切割点集合中大于此距离的所有轮廓点;S301, according to the chromosome skeleton and the skeleton intersection described in step S200, taking the skeleton intersection as the center, calculate the Euclidean distance from the cutting point to be selected to the center point, and set this distance as, excluding the set of cutting points to be selected that is greater than all contour points at this distance;

S302,根据S301的计算结果,将待选集合点根据聚类算法重新分为3个或4个子集合,统计待选切割点的集合数量,并计算每个子集合中的距离步骤S301所述的中心点的欧氏距离,选取距离最小的待选切割点确定为重叠染色体的最优切割点,计入切割点集合。S302, according to the calculation result of S301, re-divide the set points to be selected into 3 or 4 subsets according to the clustering algorithm, count the number of sets of cutting points to be selected, and calculate the distance in each subset from the center described in step S301 The Euclidean distance of the point, the candidate cutting point with the smallest distance is selected as the optimal cutting point of overlapping chromosomes, and it is included in the cutting point set.

S303,根据所述S302的计算结果,判断最优切割点的数量是否为4个,如果是,进行S304,否则进行S305;S303, according to the calculation result of S302, determine whether the number of optimal cutting points is 4, if so, go to S304, otherwise go to S305;

S304,将4个连接点两两连接,共有3种连接方式,依次判断4个切割点的连接线是否交叉,如果未交叉,即依次通过此2种连接方式分别切割重叠染色体图像。S304, connect the four connection points two by two, there are three connection methods, and sequentially determine whether the connecting lines of the four cutting points cross, if not, cut the overlapping chromosome images respectively through the two connection methods in turn.

S305,分别计算所述3个切割点的曲率,取曲率较小的两个点作为候选切割点A、B,选取曲率较大的切割点C与所述S301所述的中心点O连接成为一条线段OC;S305, calculate the curvatures of the three cutting points respectively, take two points with smaller curvatures as candidate cutting points A and B, and select the cutting point C with larger curvature to connect with the center point O described in S301 to form a line line segment OC;

S306,根据步骤S403所述的线段OC,计算线段OC的斜率并作平行线L1平行于OC,L1经过候选切割点A与染色体轮廓相交于点D1,另作平行线L2平行于OC,L2经过候选切割点B与染色体轮廓相交于点D2,将D1与D2作为候选切割点。S306, according to the line segment OC described in step S403, calculate the slope of the line segment OC and draw a parallel line L1 parallel to OC, L1 passes through the candidate cutting point A and intersects the chromosome outline at point D1, and another parallel line L2 is parallel to OC, L2 passes through The candidate cutting point B intersects the chromosome outline at point D2, and D1 and D2 are used as candidate cutting points.

S307,根据步骤S306所述的方法,获取到4个候选切割点分别为A、B、D1、D2,返回步骤S303。S307, according to the method described in step S306, four candidate cutting points are acquired as A, B, D1, and D2, respectively, and the process returns to step S303.

S400,根据步骤S300所述的切割点通过自动切割方法对所述重叠染色体进行自动分割即可。S400, the overlapping chromosomes may be automatically segmented by an automatic cutting method according to the cutting point described in step S300.

由于采用了上述技术方案,本发明具有如下有益效果:Owing to adopting the above-mentioned technical scheme, the present invention has the following beneficial effects:

本发明所述的一种重叠染色体分割方法,包含步骤S100,对重叠染色体图像进行边缘提取,并获得所有重叠染色体的边缘轮廓,将所述的边缘轮廓点全部加入到待选切割点集合;步骤S200,对重叠染色体图像提取其骨架、骨架端点和骨架交叉点;步骤S300,根据步骤S100获取的重叠染色体边缘轮廓和步骤S200获取的染色体骨架及骨架端点、交叉点确定重叠染色体的切割点;步骤S400,根据所述切割点对所述重叠染色体进行自动分割。本发明提出的基于轮廓和骨架的重叠染色体分割方法,利用重叠染色体的重叠特征通过结合染色体的轮廓和骨架的关系,将所有轮廓点看作切割点并作为候选集合,免去一系列曲率计算的复杂过程,通过欧氏距离结合染色体骨架精准定位最佳切割点。A method for dividing overlapping chromosomes according to the present invention includes step S100, performing edge extraction on overlapping chromosome images, obtaining edge contours of all overlapping chromosomes, and adding all the edge contour points to the set of cutting points to be selected; step S200, extract its skeleton, skeleton endpoints and skeleton intersections from the overlapping chromosome images; Step S300, determine the cutting points of the overlapping chromosomes according to the edge contour of the overlapping chromosomes obtained in step S100 and the chromosome skeleton, skeleton endpoints, and intersections obtained in step S200; step S400. Automatically segment the overlapping chromosomes according to the cutting points. The overlapping chromosome segmentation method based on contour and skeleton proposed by the present invention utilizes the overlapping features of overlapping chromosomes to combine the relationship between the contour and the skeleton of the chromosome, and regards all contour points as cutting points and as candidate sets, eliminating the need for a series of curvature calculations. It is a complex process that accurately locates the best cutting point through Euclidean distance combined with chromosome skeleton.

附图说明Description of drawings

图1为本发明提出的基于轮廓和骨架的重叠染色体核型图像分割方法流程图;Fig. 1 is a flowchart of a method for segmenting overlapping chromosome karyotype images based on contour and skeleton proposed by the present invention;

图2为本发明使用的重叠染色体边缘检测方法边缘检测效果图;Fig. 2 is the edge detection effect diagram of the overlapping chromosome edge detection method used in the present invention;

图3为本发明实施例1 “+”字形重叠染色体分割效果图及染色体自动分割方法流程图;3 is a flowchart of an effect diagram of “+”-shaped overlapping chromosome segmentation and a method for automatic chromosome segmentation according to Embodiment 1 of the present invention;

图4为本发明实施例2 “T”字形重叠染色体分割效果图及染色体自动分割方法流程图。FIG. 4 is an effect diagram of segmentation of “T”-shaped overlapping chromosomes and a flowchart of an automatic chromosome segmentation method according to Embodiment 2 of the present invention.

具体实施方式Detailed ways

通过下面的实施例可以详细的解释本发明,公开本发明的目的旨在保护本发明范围内的一切技术改进。The present invention can be explained in detail through the following examples, and the purpose of disclosing the present invention is to protect all technical improvements within the scope of the present invention.

结合附图1~4所述的一种重叠染色体分割方法,是基于轮廓和骨架的重叠染色体核型图像分割方法,所述分割方法具体包括以下步骤:A kind of overlapping chromosome segmentation method described in conjunction with accompanying drawing 1~4, is the overlapping chromosome karyotype image segmentation method based on outline and skeleton, and described segmentation method specifically comprises the following steps:

步骤S100,对重叠染色体图像进行边缘提取,并获得所有重叠染色体的边缘轮廓,将所述的边缘轮廓点全部加入到待选切割点集合;Step S100, perform edge extraction on the overlapping chromosome images, and obtain the edge contours of all overlapping chromosomes, and add all the edge contour points to the set of cutting points to be selected;

步骤S200,对重叠染色体图像提取其骨架、骨架端点和骨架交叉点;Step S200, extracting its skeleton, skeleton endpoints and skeleton intersections from the overlapping chromosome images;

步骤S300,根据步骤S100获取的重叠染色体边缘轮廓和步骤S200获取的染色体骨架及骨架端点、交叉点确定重叠染色体的切割点;Step S300, determining the cutting point of the overlapping chromosome according to the edge contour of the overlapping chromosome obtained in step S100 and the chromosome skeleton, skeleton endpoints, and intersection points obtained in step S200;

步骤S400,根据步骤S300所述的切割点通过自动切割方法对所述重叠染色体进行自动分割。Step S400, the overlapping chromosomes are automatically segmented by an automatic cutting method according to the cutting point described in step S300.

其中关于本发明所述的重叠染色体自动分割方法的详细步骤,下面通过两个具体的实施例对所述方法进行说明。Regarding the detailed steps of the automatic segmentation method for overlapping chromosomes according to the present invention, the method will be described below through two specific embodiments.

实施例1:Example 1:

虽然重叠染色体的形态各异,但是基本属于两种类型,包括“+”字型和“T”字型或两种类型的叠加形态,因此以下步骤将通过展示基于“+”字型重叠染色体应用本发明所述的自动分割方法。Although the shapes of overlapping chromosomes are different, they basically belong to two types, including "+" shape and "T" shape or the superimposed shape of the two types, so the following steps will show the application of overlapping chromosomes based on "+" shape The automatic segmentation method of the present invention.

本发明提出的重叠染色体分割算法,是基于轮廓和骨架的重叠染色体核型图像分割方法,所述分割方法具体包括以下步骤:The overlapping chromosome segmentation algorithm proposed by the present invention is an image segmentation method of overlapping chromosome karyotypes based on contours and skeletons, and the segmentation method specifically includes the following steps:

步骤1,在进行步骤S100边缘轮廓提取之前还需要对重叠染色体图进行预处理,首先通过平滑滤波对图像进行降噪,去除噪声影响,然后通过直方图均衡化对图像进行图像增强,提高染色体图像与背景的对比度,接下来进行二值化操作,利用阈值将染色体核型图像前景与背景分离,去除大部分的无关区域,得到预处理后的重叠染色体核型图,接下来进行步骤S100的具体步骤,利用预处理后的重叠染色体核型图进行边缘轮廓提取,此步骤具体包括以下具体步骤:Step 1: Before performing the edge contour extraction in step S100, the overlapping chromosome map needs to be preprocessed. First, the image is denoised by smoothing filtering to remove the influence of noise, and then the image is enhanced by histogram equalization to improve the chromosome image. Contrast with the background, then perform a binarization operation, use a threshold to separate the foreground and background of the karyotype image, remove most of the irrelevant regions, and obtain the preprocessed overlapping karyotype map, and then perform step S100. step, using the preprocessed overlapping chromosome karyotype map to perform edge contour extraction, and this step specifically includes the following specific steps:

步骤1.1,首先采用高斯滤波算法将算子与原图进行卷积运算,对原图进行平滑滤波,得到消除高斯噪声的输出图像。Step 1.1, first use the Gaussian filtering algorithm to perform a convolution operation on the operator and the original image, and perform smooth filtering on the original image to obtain an output image with Gaussian noise removed.

步骤1.2,然后将步骤1.1的输出图像与梯度算子进行卷积运算,通过计算得到梯度幅值及梯度方向。Step 1.2, then perform a convolution operation on the output image of step 1.1 and the gradient operator, and obtain the gradient magnitude and gradient direction through calculation.

步骤1.3,接下来对进行非极大值抑制处理,以图像中像素为原点,分别按45度、135度作过原点的双向分割线,将原点像素值与沿着分割线上相邻两个像素的像素值对比,若原点像素值较大则保留;否则,归零。Step 1.3, the next step is to perform non-maximum suppression processing. Taking the pixel in the image as the origin, the bidirectional dividing line passing through the origin is at 45 degrees and 135 degrees, respectively. The pixel value comparison of the pixel, if the pixel value of the origin is larger, it is retained; otherwise, it is returned to zero.

步骤1.4,最后进行双阈值滞后化阈值处理,设定高低(Th、Tl)阈值计算因子,将高于Th的标记为强边缘点,小于Tl的点置0,介于高低阈值之间的像素点且与强边缘点不相连的点置为0,最终获取到重叠染色体得边缘轮廓,将全部轮廓点记入候选切割点集合P,附图1为染色体核型图边缘检测效果图。Step 1.4: Finally, double-threshold hysteresis threshold processing is performed, and the high and low (Th, Tl) threshold calculation factors are set, and the points higher than Th are marked as strong edge points, the points less than Tl are set to 0, and the pixels between the high and low thresholds are set to 0. The point that is not connected to the strong edge point is set to 0, and finally the edge contour of the overlapping chromosome is obtained, and all the contour points are recorded in the candidate cutting point set P. Figure 1 shows the effect of the edge detection of the chromosome karyotype map.

步骤2,根据步骤1预处理后的重叠染色体核型图像使用以下方法提取染色体骨架,此方法具体步骤包括:首先根据细化算法逐步计算,将染色体图像细化为染色体骨架,然后根据像素点的连接关系定位到染色体的骨架端点,依据端点的三侧均没有像素值的划分标准,并记录在端点集合中,然后根据像素点的判别方法,提取到骨架交叉点的所在位置,并记录到骨架的交叉点集合中。In step 2, the following method is used to extract the chromosome skeleton according to the preprocessed overlapping chromosome karyotype image in step 1. The specific steps of this method include: firstly calculate step by step according to the thinning algorithm, and refine the chromosome image into the chromosome skeleton, and then according to the pixel point The connection relationship is located at the skeleton end point of the chromosome, according to the division standard of no pixel value on the three sides of the end point, and recorded in the end point set, and then according to the pixel point discrimination method, the position of the skeleton intersection point is extracted and recorded in the skeleton in the intersection set.

步骤3,结合染色体二值图像边缘轮廓和骨架交叉点筛选分割点,通过筛选得到精准的重叠染色体分割点,具体步骤包含:Step 3, combine the edge contour of the chromosome binary image and the skeleton intersection to screen the segmentation points, and obtain accurate overlapping chromosome segmentation points through screening. The specific steps include:

步骤3.1,根据步骤1以及步骤2所述的染色体骨架和骨架交叉点,以骨架交叉点为中心,计算待选切割点到中心点的欧式距离,为此距离进行设定,去除待选切割点集合中大于此距离的所有轮廓点,得到筛选后的待选切割点,如附图2第4幅图片所示,为根据步骤3.1进行分割后筛选得到的待选切割点示意图。Step 3.1, according to the chromosome skeleton and the skeleton intersection described in Step 1 and Step 2, take the skeleton intersection as the center, calculate the Euclidean distance from the cutting point to be selected to the center point, set this distance, and remove the cutting point to be selected All contour points in the set that are larger than this distance, obtain the selected cutting points after screening, as shown in the fourth picture in Figure 2, which is a schematic diagram of the selected cutting points obtained after dividing according to step 3.1.

步骤3.2,根据步骤3.1的计算结果,将待选集合点根据聚类算法重新分为4个子集合,统计待选切割点的集合数量,并计算每个子集合中的距离步骤S301所述的中心点的欧氏距离,选取距离最小的待选切割点确定为重叠染色体的最优切割点,计入切割点集合,如附图2第5幅图片所示,为根据步骤3.2所得的最优切割点示意图。Step 3.2, according to the calculation result of step 3.1, re-divide the set points to be selected into 4 subsets according to the clustering algorithm, count the number of sets of cutting points to be selected, and calculate the distance in each subset from the center point described in step S301. The Euclidean distance of , select the candidate cutting point with the smallest distance and determine it as the optimal cutting point of overlapping chromosomes, which is included in the cutting point set, as shown in the fifth picture in Figure 2, which is the optimal cutting point obtained according to step 3.2 Schematic.

步骤3.3,判断最优切割点是否为4个,若待选最优切割点为4个,执行步骤3.4。Step 3.3, determine whether there are 4 optimal cutting points, if there are 4 optimal cutting points to be selected, go to step 3.4.

步骤3.4,将4个连接点两两连接,共有3种连接方式,依次判断4个切割点的连接线是否交叉,如果未交叉,即依次通过此2种连接方式分别切割重叠染色体图像。Step 3.4, connect the 4 connection points two by two, there are 3 connection methods, and then judge whether the connection lines of the 4 cutting points cross, if not, then cut the overlapping chromosome images respectively through the two connection methods.

步骤4,根据步骤3.3所述的最优切割点通过自动切割方法对所述重叠染色体进行自动分割。Step 4, according to the optimal cutting point described in Step 3.3, the overlapping chromosomes are automatically segmented by an automatic cutting method.

实施例2:Example 2:

下面通过展示基于“T”字型重叠染色体的实施例2,在应用本发明中所述的染色体核型图像自动分割方法过程中,具体步骤包括以下内容:Below by showing the embodiment 2 based on "T" type overlapping chromosomes, in the process of applying the automatic segmentation method for chromosome karyotype images described in the present invention, the specific steps include the following:

步骤1,在进行步骤S100边缘轮廓提取之前还需要对重叠染色体图进行预处理。首先通过平滑滤波对图像进行降噪,去除噪声影响,然后通过直方图均衡化对图像进行图像增强,提高染色体图像与背景的对比度,接下来进行二值化操作,利用阈值将染色体核型图像前景与背景分离,去除大部分的无关区域,得到预处理后的重叠染色体核型图,接下来进行步骤S100的具体步骤,利用预处理后的重叠染色体核型图进行边缘轮廓提取,此步骤具体包括以下具体步骤:In step 1, the overlapping chromosome map needs to be preprocessed before the edge contour extraction in step S100 is performed. First, the image is denoised by smoothing filtering to remove the influence of noise, and then the image is enhanced by histogram equalization to improve the contrast between the chromosome image and the background. Separating from the background, removing most of the irrelevant regions, and obtaining the preprocessed overlapping chromosome karyotype map, and then performing the specific steps of step S100, using the preprocessed overlapping chromosome karyotype map to perform edge contour extraction, and this step specifically includes The following specific steps:

步骤1.1,首先采用高斯滤波算法将算子与原图进行卷积运算,对原图进行平滑滤波,得到消除高斯噪声的输出图像。Step 1.1, first use the Gaussian filtering algorithm to perform a convolution operation on the operator and the original image, and perform smooth filtering on the original image to obtain an output image with Gaussian noise removed.

步骤1.2,然后将步骤1.1的输出图像与梯度算子进行卷积运算,通过计算得到梯度幅值及梯度方向。Step 1.2, then perform a convolution operation on the output image of step 1.1 and the gradient operator, and obtain the gradient magnitude and gradient direction through calculation.

步骤1.3,接下来对进行非极大值抑制处理,以图像中像素为原点,分别按45度、135度作过原点的双向分割线,将原点像素值与沿着分割线上相邻两个像素的像素值对比,若原点像素值较大则保留;否则,归零。Step 1.3, the next step is to perform non-maximum suppression processing. Taking the pixel in the image as the origin, the bidirectional dividing line passing through the origin is at 45 degrees and 135 degrees, respectively. The pixel value comparison of the pixel, if the pixel value of the origin is larger, it is retained; otherwise, it is returned to zero.

步骤1.4,最后进行双阈值滞后化阈值处理,设定高低(Th、Tl)阈值计算因子,将高于Th的标记为强边缘点,小于Tl的点置0,介于高低阈值之间的像素点且与强边缘点不相连的点置为0,最终获取到重叠染色体得边缘轮廓,将全部轮廓点记入候选切割点集合P,附图1为染色体核型图边缘检测效果图。Step 1.4: Finally, double-threshold hysteresis threshold processing is performed, and the high and low (Th, Tl) threshold calculation factors are set, and the points higher than Th are marked as strong edge points, the points less than Tl are set to 0, and the pixels between the high and low thresholds are set to 0. The point that is not connected to the strong edge point is set to 0, and finally the edge contour of the overlapping chromosome is obtained, and all the contour points are recorded in the candidate cutting point set P. Figure 1 shows the effect of the edge detection of the chromosome karyotype map.

步骤2,根据步骤1预处理后的重叠染色体核型图像使用以下方法提取染色体骨架,此方法具体步骤包括:首先根据细化算法逐步计算,将染色体图像细化为染色体骨架,然后根据像素点的连接关系定位到染色体的骨架端点,依据端点的三侧均没有像素值的划分标准,并记录在端点集合中,然后根据像素点的判别方法,提取到骨架交叉点的所在位置,并记录到骨架的交叉点集合中。In step 2, the following method is used to extract the chromosome skeleton according to the preprocessed overlapping chromosome karyotype image in step 1. The specific steps of this method include: firstly calculate step by step according to the thinning algorithm, and refine the chromosome image into the chromosome skeleton, and then according to the pixel point The connection relationship is located at the skeleton end point of the chromosome, according to the division standard of no pixel value on the three sides of the end point, and recorded in the end point set, and then according to the pixel point discrimination method, the position of the skeleton intersection point is extracted and recorded in the skeleton in the intersection set.

步骤3,结合染色体二值图像边缘轮廓和骨架交叉点筛选分割点,通过筛选得到精准的重叠染色体分割点,具体步骤包含:Step 3, combine the edge contour of the chromosome binary image and the skeleton intersection to screen the segmentation points, and obtain accurate overlapping chromosome segmentation points through screening. The specific steps include:

步骤3.1,根据步骤1以及步骤2所述的染色体骨架和骨架交叉点,以骨架交叉点为中心,计算待选切割点到中心点的欧式距离,为此距离进行设定,去除待选切割点集合中大于此距离的所有轮廓点,得到筛选后的待选切割点,如图3中第4幅图片所示,为根据步骤3.1进行筛选后得到的待选切割点示意图。Step 3.1, according to the chromosome skeleton and the skeleton intersection described in Step 1 and Step 2, take the skeleton intersection as the center, calculate the Euclidean distance from the cutting point to be selected to the center point, set this distance, and remove the cutting point to be selected All the contour points in the set that are greater than this distance get the selected cutting points after screening, as shown in the fourth picture in Figure 3, which is a schematic diagram of the selected cutting points obtained after screening according to step 3.1.

步骤3.2,根据步骤3.1的计算结果,将待选集合点根据聚类算法重新分为3个子集合,统计待选切割点的集合数量,并计算每个子集合中的距离步骤S301所述的中心点的欧氏距离,选取距离最小的待选切割点确定为重叠染色体的最优切割点,计入切割点集合,如附图3第5幅图片所示,为根据步骤3.2所得的最优切割点示意图。Step 3.2, according to the calculation result of step 3.1, re-divide the set points to be selected into 3 subsets according to the clustering algorithm, count the number of sets of cutting points to be selected, and calculate the distance in each subset from the center point described in step S301. The Euclidean distance, select the candidate cutting point with the smallest distance to determine the optimal cutting point of overlapping chromosomes, and be included in the cutting point set. Schematic.

步骤3.3,判断最优切割点是否为4个,若待选最优切割点为4个,执行步骤3.4,根据实施例2所述,步骤3.2只能得到3个最优切割点,因此执行步骤3.5。Step 3.3, determine whether there are 4 optimal cutting points, if the number of optimal cutting points to be selected is 4, go to step 3.4. According to Embodiment 2, only 3 optimal cutting points can be obtained in step 3.2, so perform step 3.2. 3.5.

步骤3.4,将4个连接点两两连接,共有3种连接方式,依次判断4个切割点的连接线是否交叉,如果未交叉,即依次通过此2种连接方式分别切割重叠染色体图像,并执行步骤4。Step 3.4, connect the 4 connection points in pairs, there are 3 connection methods, and then judge whether the connection lines of the 4 cutting points intersect, if not, cut the overlapping chromosome images in turn through these 2 connection methods, and execute Step 4.

步骤3.5,分别计算所述3个切割点的曲率,取曲率较小的两个点作为候选切割点A、B,选取曲率较大的切割点C与步骤2中所得到的中心点O连接成为一条线段OC;Step 3.5: Calculate the curvature of the three cutting points respectively, take the two points with smaller curvature as candidate cutting points A and B, and select the cutting point C with larger curvature to connect with the center point O obtained in step 2 as A line segment OC;

步骤3.6,根据步骤3.5所得到的线段OC,计算线段OC的斜率并作平行线L1平行于OC,L1经过候选切割点A与染色体轮廓相交于点D1,另作平行线L2平行于OC,L2经过候选切割点B与染色体轮廓相交于点D2,将D1与D2作为候选切割点,并返回步骤3.4。Step 3.6, according to the line segment OC obtained in step 3.5, calculate the slope of the line segment OC and make a parallel line L1 parallel to OC, L1 passes through the candidate cutting point A and intersects the chromosome outline at point D1, and another parallel line L2 is parallel to OC, L2 After the candidate cutting point B intersects the chromosome outline at point D2, take D1 and D2 as candidate cutting points, and return to step 3.4.

步骤4,根据步骤3.4所述的最优切割点连接方法,通过自动切割方法对所述重叠染色体进行自动分割。In step 4, according to the optimal cutting point connection method described in step 3.4, the overlapping chromosomes are automatically segmented by an automatic cutting method.

本发明未详述部分为现有技术,尽管结合优选实施方案具体展示和介绍了本发明,具体实现该技术方案方法和途径很多,以上所述仅是本发明的优选实施方式,但所属领域的技术人员应该明白,在不脱离所附权利要求书所限定的本发明的精神和范围内,在形式上和细节上可以对本发明做出各种变化,均为本发明的保护范围。The part that is not described in detail in the present invention is the prior art. Although the present invention is specifically shown and introduced in conjunction with the preferred embodiments, there are many methods and approaches to realize the technical solution. The above are only the preferred embodiments of the present invention, but the It should be understood by those skilled in the art that various changes can be made to the present invention in form and detail without departing from the spirit and scope of the present invention defined by the appended claims, which are all within the protection scope of the present invention.

Claims (7)

1. A method for segmentation of overlapping chromosome karyotypes, characterized in that automatic segmentation is performed on overlapping chromosomes, said segmentation method comprising the steps of:
step S100, carrying out edge extraction on the overlapped chromosome image, obtaining edge outlines of all overlapped chromosomes, and adding all edge outline points into a set of cutting points to be selected;
step S200, extracting a skeleton, a skeleton end point and a skeleton intersection point of the overlapped chromosome image;
step S300, determining the cutting points of the overlapped chromosomes according to the edge contour of the overlapped chromosomes obtained in the step S100 and the chromosome skeleton, the skeleton end points and the intersection points obtained in the step S200;
and S400, automatically segmenting the overlapped chromosomes through an automatic segmentation method according to the segmentation points in the step S300.
2. The method of overlapping chromosome segmentation as set forth in claim 1, wherein: the method used in the overlapping chromosome edge contour extraction method in step S100 is a contour extraction method based on the Canny operator.
3. The method of overlapping chromosome segmentation as set forth in claim 1, wherein: the above-mentioned
The method comprises the following steps of extracting a skeleton, skeleton endpoints and skeleton intersections of overlapped chromosome images, wherein the step S200 specifically comprises the following steps:
step S201, a thinning algorithm is adopted to thin the overlapped chromosome image into a chromosome skeleton with only one pixel point, and the overlapped chromosome skeleton is extracted;
step S202, positioning the skeleton endpoint, namely the number of pixels in the neighborhood is 1 according to the number of the refined adjacent pixels;
step S203, positioning the framework cross point according to the number of the neighborhood pixel points.
4. The method of overlapping chromosome segmentation as set forth in claim 1, wherein: the automatic segmentation method of the overlapped chromosomes is in a shape of a plus sign and a T or the superposition of the two types.
5. The method of overlapping chromosome segmentation according to claim 4, wherein: the method for determining the overlapped chromosome cutting point by using the overlapped chromosome edge contour method obtained in the step S100 and the chromosome skeleton, the skeleton end points and the intersection points obtained in the step S200 specifically includes the following steps:
s301, calculating the Euclidean distance from the cutting point to be selected to the central point by taking the skeleton intersection point as the center according to the chromosome skeleton and the skeleton intersection point in the step S200, setting the distance, and removing all contour points which are greater than the distance in the cutting point set to be selected;
s302, according to the calculation result of S301, subdividing the to-be-selected set points into 3 or 4 subsets according to a clustering algorithm, counting the set number of the to-be-selected cut points, calculating Euclidean distances between the to-be-selected cut points and the central point in the step S301 in each subset, selecting the to-be-selected cut point with the minimum distance, determining the to-be-selected cut point as the optimal cut point of the overlapped chromosomes, and counting the to-be-selected cut point into a cut point set;
s303, judging whether the number of the optimal cutting points is 4 or not according to the calculation result of the S302, if so, carrying out S304, and if not, carrying out S305;
s304, connecting the 4 connecting points pairwise, wherein 3 connecting modes are provided, sequentially judging whether the connecting lines of the 4 cutting points are crossed, and if not, respectively cutting the overlapped chromosome images sequentially through the 2 connecting modes;
s305, respectively calculating curvatures of the 3 cutting points, taking two points with smaller curvatures as candidate cutting points A, B, and selecting a cutting point C with larger curvature to be connected with the central point O in the S301 to form a line segment OC;
s306, according to the line segment OC in the step S403, calculating the slope of the line segment OC, drawing a parallel line L1 parallel to the OC, making L1 pass through a candidate cut point A to intersect with the chromosome contour at a point D1, drawing a parallel line L2 parallel to the OC, making L2 pass through a candidate cut point B to intersect with the chromosome contour at a point D2, and taking D1 and D2 as candidate cut points;
s307, according to the method in the step S306, acquiring that the 4 candidate cutting points are A, B, D1 and D2 respectively, and returning to the step S303;
s400, automatically segmenting the overlapped chromosomes according to the segmentation points in the step S300 by an automatic segmentation method.
6. The method of overlapping chromosome segmentation as set forth in claim 5, wherein: the automatic segmentation method based on the application of the '+' font overlapped chromosome, namely the overlapped chromosome karyotype image segmentation method based on the contour and the skeleton, specifically comprises the following steps:
step 1, before performing the edge contour extraction in step S100, preprocessing an overlapping chromosome map, first performing noise reduction on an image by smooth filtering to remove noise influence, then performing image enhancement on the image by histogram equalization to improve the contrast between a chromosome image and a background, next performing binarization operation, separating a foreground of a chromosome karyotype image from the background by using a threshold value, removing most of irrelevant areas to obtain a preprocessed overlapping chromosome karyotype map, next performing the specific steps in step S100, and performing the edge contour extraction by using the preprocessed overlapping chromosome karyotype map, wherein the specific steps specifically include the following specific steps:
step 1.1, firstly, carrying out convolution operation on an operator and an original image by adopting a Gaussian filtering algorithm, and carrying out smooth filtering on the original image to obtain an output image with Gaussian noise eliminated;
step 1.2, performing convolution operation on the output image in the step 1.1 and a gradient operator, and obtaining a gradient amplitude and a gradient direction through calculation;
step 1.3, then carrying out non-maximum value suppression treatment, taking pixels in the image as an origin, respectively making bidirectional dividing lines passing through the origin at 45 degrees and 135 degrees, comparing the pixel value of the origin with pixel values of two adjacent pixels along the dividing lines, and if the pixel value of the origin is larger, keeping the original; otherwise, returning to zero;
step 1.4, finally, carrying out double-threshold hysteresis threshold processing, setting high and low (Th and Tl) threshold calculation factors, marking the points higher than Th as strong edge points, setting the points smaller than Tl as 0, setting the points between the high and low thresholds and the points which are not connected with the strong edge points as 0, finally obtaining overlapped chromosomes to obtain edge contours, and recording all contour points into a candidate cut point set P;
step 2, extracting a chromosome skeleton according to the overlapping chromosome karyotype image preprocessed in the step 1.1-1.4 by using the following method, wherein the method specifically comprises the following steps: firstly, gradually calculating according to a thinning algorithm, thinning a chromosome image into a chromosome skeleton, then positioning skeleton endpoints of the chromosome according to the connection relation of pixel points, recording in an endpoint set according to a division standard that no pixel value exists on three sides of the endpoints, then extracting the positions of skeleton intersections according to a pixel point discrimination method, and recording in an intersection set of the skeleton;
and 3, screening the segmentation points by combining the chromosome binary image edge contour and the skeleton intersection point, and obtaining accurate overlapped chromosome segmentation points through screening, wherein the method specifically comprises the following steps:
step 3.1, calculating the Euclidean distance from the cutting point to be selected to the central point by taking the skeleton intersection point as the center according to the chromosome skeleton and the skeleton intersection point in the steps 1.1-1.4 and 2, setting the Euclidean distance for the distance, removing all contour points which are greater than the distance in the cutting point set to be selected, and obtaining the screened cutting point to be selected;
step 3.2, according to the calculation result of the step 3.1, the set points to be selected are divided into 4 subsets again according to a clustering algorithm, the set number of the cutting points to be selected is counted, the Euclidean distance from the central point in the step S301 in each subset is calculated, the cutting point to be selected with the minimum distance is selected to be determined as the optimal cutting point of the overlapped chromosomes, and the optimal cutting point is counted into a cutting point set;
step 3.3, judging whether the optimal cutting points are 4, and if the optimal cutting points to be selected are 4, executing the step 3.4;
step 3.4, connecting the 4 connecting points pairwise, wherein 3 connecting modes are provided, sequentially judging whether the connecting lines of the 4 cutting points are crossed, and if not, respectively cutting the overlapped chromosome images sequentially through the 2 connecting modes;
and 4, automatically segmenting the overlapped chromosomes by an automatic segmentation method according to the optimal segmentation points in the step 3.3.
7. The method of overlapping chromosome segmentation as set forth in claim 5, wherein: the application based on the T-shaped overlapping chromosome in the automatic segmentation method process of the chromosome karyotype image comprises the following specific steps:
step 1, before performing the edge contour extraction in step S100, preprocessing an overlapping chromosome map, first performing noise reduction on an image by smooth filtering to remove noise influence, then performing image enhancement on the image by histogram equalization to improve the contrast between a chromosome image and a background, next performing binarization operation, separating a foreground of a chromosome karyotype image from the background by using a threshold value, removing most of irrelevant areas to obtain a preprocessed overlapping chromosome karyotype map, next performing the specific steps in step S100, and performing the edge contour extraction by using the preprocessed overlapping chromosome karyotype map, wherein the specific steps specifically include the following specific steps:
step 1.1, firstly, carrying out convolution operation on an operator and an original image by adopting a Gaussian filtering algorithm, and carrying out smooth filtering on the original image to obtain an output image with Gaussian noise eliminated;
step 1.2, performing convolution operation on the output image in the step 1.1 and a gradient operator, and obtaining a gradient amplitude and a gradient direction through calculation;
step 1.3, then carrying out non-maximum value suppression treatment, taking pixels in the image as an origin, respectively making bidirectional dividing lines passing through the origin at 45 degrees and 135 degrees, comparing the pixel value of the origin with pixel values of two adjacent pixels along the dividing lines, and if the pixel value of the origin is larger, keeping the original; otherwise, returning to zero;
step 1.4, finally, carrying out double-threshold hysteresis threshold processing, setting high and low (Th and Tl) threshold calculation factors, marking the points higher than Th as strong edge points, setting the points smaller than Tl as 0, setting the points between the high and low thresholds and the points which are not connected with the strong edge points as 0, finally obtaining overlapped chromosomes to obtain edge contours, and recording all contour points into a candidate cut point set P;
step 2, extracting a chromosome skeleton according to the overlapping chromosome karyotype image preprocessed in the step 1.1-1.4 by using the following method, wherein the method specifically comprises the following steps: firstly, gradually calculating according to a thinning algorithm, thinning a chromosome image into a chromosome skeleton, then positioning skeleton endpoints of the chromosome according to the connection relation of pixel points, recording in an endpoint set according to a division standard that no pixel value exists on three sides of the endpoints, then extracting the positions of skeleton intersections according to a pixel point discrimination method, and recording in an intersection set of the skeleton;
and 3, screening the segmentation points by combining the chromosome binary image edge contour and the skeleton intersection point, and obtaining accurate overlapped chromosome segmentation points through screening, wherein the method specifically comprises the following steps:
step 3.1, calculating the Euclidean distance from the cutting point to be selected to the central point by taking the skeleton intersection point as the center according to the chromosome skeleton and the skeleton intersection point in the steps 1.1-1.4 and 2, setting the Euclidean distance for the distance, removing all contour points which are greater than the distance in the cutting point set to be selected, and obtaining the screened cutting point to be selected;
step 3.2, according to the calculation result of the step 3.1, the set points to be selected are divided into 3 subsets again according to a clustering algorithm, the set number of the cutting points to be selected is counted, the Euclidean distance between the set points to be selected and the central point in the step S301 in each subset is calculated, the cutting point to be selected with the minimum distance is selected to be determined as the optimal cutting point of the overlapped chromosomes, and the optimal cutting point is counted into a cutting point set;
step 3.3, judging whether the optimal cutting points are 4, if the optimal cutting points to be selected are 4, executing step 3.4, and according to the embodiment 2, only 3 optimal cutting points can be obtained in the step 3.2, so that the step 3.5 is executed;
step 3.4, connecting the 4 connecting points pairwise, wherein 3 connecting modes are provided, sequentially judging whether the connecting lines of the 4 cutting points are crossed, if not, respectively cutting the overlapped chromosome images sequentially through the 2 connecting modes, and executing the step 4; step 3.5, respectively calculating the curvatures of the 3 cutting points, taking two points with smaller curvatures as candidate cutting points A, B, and selecting a cutting point C with larger curvature to be connected with the central point O obtained in the step 2 to form a line segment OC;
step 3.6, according to the line segment OC obtained in step 3.5, calculating the slope of the line segment OC, drawing a parallel line L1 parallel to OC, wherein L1 crosses the chromosome contour through the candidate cut point a at a point D1, drawing a parallel line L2 parallel to OC, wherein L2 crosses the chromosome contour through the candidate cut point B at a point D2, taking D1 and D2 as candidate cut points, and returning to step 3.4;
and 4, automatically segmenting the overlapped chromosomes by an automatic segmentation method according to the optimal cutting point connection method in the step 3.4.
CN202210548618.6A 2022-05-20 2022-05-20 Overlapping chromosome segmentation method Pending CN114897843A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210548618.6A CN114897843A (en) 2022-05-20 2022-05-20 Overlapping chromosome segmentation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210548618.6A CN114897843A (en) 2022-05-20 2022-05-20 Overlapping chromosome segmentation method

Publications (1)

Publication Number Publication Date
CN114897843A true CN114897843A (en) 2022-08-12

Family

ID=82723418

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210548618.6A Pending CN114897843A (en) 2022-05-20 2022-05-20 Overlapping chromosome segmentation method

Country Status (1)

Country Link
CN (1) CN114897843A (en)

Similar Documents

Publication Publication Date Title
WO2022063198A1 (en) Lung image processing method, apparatus and device
WO2022063199A1 (en) Pulmonary nodule automatic detection method, apparatus and computer system
CN106056118B (en) A kind of identification method of counting for cell
CN104992445B (en) A kind of automatic division method of CT images pulmonary parenchyma
CN112215800B (en) Overlapping Chromosome Identification and Segmentation Method Based on Machine Learning
CN104392460B (en) A Segmentation Method of Adhesive Leukocytes Based on Nuclear Marker Watershed Transform
CN106651846B (en) Segmentation method of retinal blood vessel images
CN103984958B (en) Cervical cancer cell dividing method and system
CN108133476B (en) Method and system for automatically detecting pulmonary nodules
CN107492088B (en) Automatic identification and statistics method for white blood cells in gynecological microscopic image
CN111402267A (en) Method, device and terminal for segmentation of epithelial cell nuclei in pathological images of prostate cancer
CN107798679A (en) Breast molybdenum target image breast area is split and tufa formation method
CN110544262B (en) A method of image segmentation of cervical cells based on machine vision
CN111369530A (en) CT image pulmonary nodule rapid screening method based on deep learning
CN104778442A (en) Automatic segmentation and counting method of retina cell fluorescence microscopic image
CN108765409A (en) A kind of screening technique of the candidate nodule based on CT images
CN110400287A (en) System and method for detecting tumor invasion edge and center in IHC staining images of colorectal cancer
CN116718599B (en) A method for measuring apparent crack length based on three-dimensional point cloud data
CN111524154A (en) An image-based automatic segmentation method for tunnel segments
CN112396618B (en) Grain boundary extraction and grain size measurement method based on image processing
CN104933723B (en) Tongue Image Segmentation Method Based on Sparse Representation
CN114693672A (en) A method for removing skin glands and nipples from mammography images based on image processing
CN108257118B (en) Fracture adhesion segmentation method based on normal corrosion and random walk
CN113080843B (en) Meibomian gland image-based gland extraction method and quantitative analysis method
CN106203456A (en) Coal dust Algorithm for Overlapping Granule separation method based on improved differential evolution particle cluster algorithm

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination