CN114897843A - Overlapping chromosome segmentation method - Google Patents
Overlapping chromosome segmentation method Download PDFInfo
- 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
Links
- 210000000349 chromosome Anatomy 0.000 title claims abstract description 182
- 238000000034 method Methods 0.000 title claims abstract description 76
- 230000011218 segmentation Effects 0.000 title claims abstract description 45
- 238000004364 calculation method Methods 0.000 claims abstract description 18
- 238000000605 extraction Methods 0.000 claims abstract description 15
- 238000001914 filtration Methods 0.000 claims description 12
- 238000012216 screening Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 8
- 230000002457 bidirectional effect Effects 0.000 claims description 4
- 238000012850 discrimination method Methods 0.000 claims description 4
- 230000001629 suppression Effects 0.000 claims description 4
- 238000003709 image segmentation Methods 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims 2
- 238000004458 analytical method Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000003708 edge detection Methods 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 2
- 230000002759 chromosomal effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 230000000717 retained effect Effects 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 206010008805 Chromosomal abnormalities Diseases 0.000 description 1
- 208000031404 Chromosome Aberrations Diseases 0.000 description 1
- 208000026350 Inborn Genetic disease Diseases 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 239000012634 fragment Substances 0.000 description 1
- 208000016361 genetic disease Diseases 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000011278 mitosis Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 230000008722 morphological abnormality Effects 0.000 description 1
- 210000004940 nucleus Anatomy 0.000 description 1
- 102000004169 proteins and genes Human genes 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 238000010186 staining Methods 0.000 description 1
- 230000005945 translocation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/12—Edge-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/13—Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10056—Microscopic image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20084—Artificial 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
Description
技术领域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)
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) |
-
2022
- 2022-05-20 CN CN202210548618.6A patent/CN114897843A/en active Pending
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 |