[go: up one dir, main page]

CN110610490A - A method for locating leukocytes in images of diseased cells - Google Patents

A method for locating leukocytes in images of diseased cells Download PDF

Info

Publication number
CN110610490A
CN110610490A CN201910856026.9A CN201910856026A CN110610490A CN 110610490 A CN110610490 A CN 110610490A CN 201910856026 A CN201910856026 A CN 201910856026A CN 110610490 A CN110610490 A CN 110610490A
Authority
CN
China
Prior art keywords
image
cell
white blood
threshold
segmentation
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
CN201910856026.9A
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.)
Harbin University of Science and Technology
Original Assignee
Harbin University of Science and Technology
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 Harbin University of Science and Technology filed Critical Harbin University of Science and Technology
Priority to CN201910856026.9A priority Critical patent/CN110610490A/en
Publication of CN110610490A publication Critical patent/CN110610490A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/20Linear translation of whole images or parts thereof, e.g. panning
    • 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/136Segmentation; Edge detection involving thresholding
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/70Determining position or orientation of objects or cameras
    • 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/10024Color image
    • 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

Landscapes

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

Abstract

一种病变细胞图像中白细胞的定位方法,涉及医学诊断领域,在混有背景、红细胞、白细胞浆和白细胞核的显微图像中,通过阈值分割算法分离出细胞核,利用开运算重构去除孤立点,用膨胀处理去除无关区域,细胞核的分离实现了白细胞的初步定位。对于非黏连白细胞,直接使用阈值分割算法对白细胞浆进行分割。对于黏连白细胞,分割白细胞浆时,先用最大类间方差算法将白细胞背景和目标分成两部分,使错分概率最小,其次,利用Sobel求取梯度图像,根据获取的梯度幅值做分水岭变换,最终对获取的灰度图像做基于重建的开—闭运算,因为在不影响目标图像的前提下,去除图像的碎片效果更好,实现白细胞分割。本发明用于医学临床诊断,提高医务人员诊断效率。A method for locating leukocytes in a diseased cell image, which relates to the field of medical diagnosis. In a microscopic image mixed with background, red blood cells, leukocyte plasma and leukocyte nuclei, the nuclei are separated by a threshold segmentation algorithm, and the isolated points are removed by open operation reconstruction. , with expansion treatment to remove irrelevant areas, and isolation of nuclei to achieve preliminary localization of leukocytes. For non-adherent leukocytes, the leukocyte plasma was segmented directly using a threshold segmentation algorithm. For adherent leukocytes, when segmenting leukocyte plasma, first use the maximum inter-class variance algorithm to divide the leukocyte background and the target into two parts to minimize the probability of misclassification. Second, use Sobel to obtain the gradient image, and perform watershed transformation according to the obtained gradient amplitude. , and finally perform the reconstruction-based open-close operation on the obtained grayscale image, because without affecting the target image, the effect of removing image fragments is better, and white blood cell segmentation can be achieved. The invention is used for medical clinical diagnosis and improves the diagnosis efficiency of medical personnel.

Description

一种病变细胞图像中白细胞的定位方法A method for locating leukocytes in images of diseased cells

技术领域technical field

本发明涉及医学图像处理领域,具体涉及了一种病变细胞图像中白细胞的定位方法。The invention relates to the field of medical image processing, in particular to a method for locating leukocytes in a diseased cell image.

背景技术Background technique

医学图像处理技术对医学科研及临床实践的作用和影像日益增大,其结果使临床医生对人体内部病变部位的观察更直接、更清晰,准确率也更高。因此,医学图像处理技术一直受到国内外有关专家的高度重视,本文也在医学处理方面提出了一种白细胞的自动定位方法,通过在混有红细胞和血小板的数据图像中将白细胞细胞核分割提取出来,初步获得白细胞的位置,根据面积最大和细胞核的位置获得白细胞的具体位置,并将白细胞分割出来。The role and imaging of medical image processing technology in medical research and clinical practice is increasing day by day, and as a result, clinicians can observe the internal lesions of the human body more directly, clearly, and with higher accuracy. Therefore, medical image processing technology has always been highly valued by relevant experts at home and abroad. In this paper, an automatic leukocyte localization method is also proposed in medical processing. The location of the leukocytes was initially obtained, and the specific location of the leukocytes was obtained according to the largest area and the location of the nucleus, and the leukocytes were segmented.

传统的白细胞分割方法采用人工设计的特征,受光照等外界因素的影响,对于不同的白细胞图像使用的方法较为局限。在部分场景中传统的方法无法得到令人满意的结果,鲁棒性较差。The traditional leukocyte segmentation method adopts artificially designed features, which are affected by external factors such as illumination, and the methods used for different leukocyte images are relatively limited. In some scenarios, traditional methods cannot obtain satisfactory results and have poor robustness.

本发明提出对医学病变细胞图像中白细胞的定位方法,该方法通过阈值分割可以将黏连细胞和非黏连细胞的细胞核从显微图像中分离出来,对于黏连的白细胞采用最大类间方差算法,将图像分成两部分,再通过分水岭算法分割白细胞浆,对于黏连的细胞简单一些,可直接通过阈值分割的算法将其分离出来。病变细胞图像中白细胞的定位方法的提出主要目的是为了使白细胞在临床诊断中定位更准确,使临床医生对人体内部病变部位的观察更直接、清晰,准确率也更高,从而使得医护人员的工作效率得以提高。The present invention proposes a method for locating leukocytes in medical diseased cell images. The method can separate the nuclei of adherent cells and non-adherent cells from the microscopic image through threshold segmentation, and uses the maximum inter-class variance algorithm for adherent leukocytes. , the image is divided into two parts, and then the leukocyte plasma is segmented by the watershed algorithm. For the adhering cells, it is simpler and can be directly separated by the threshold segmentation algorithm. The main purpose of the method for locating leukocytes in diseased cell images is to make leukocytes more accurate in clinical diagnosis, so that clinicians can observe the diseased parts in the human body more directly, clearly, and with higher accuracy, thereby making the medical staff more accurate. Work efficiency is improved.

发明内容SUMMARY OF THE INVENTION

为了解决背景技术中存在的问题,本发明的目的在于提供了一种病变细胞图像中白细胞的定位方法。In order to solve the problems existing in the background art, the purpose of the present invention is to provide a method for locating leukocytes in an image of diseased cells.

本发明通过阈值分割算法分离出细胞核,利用开运算重构去除孤立点,用膨胀处理去除无关区域,细胞核的分离实现了白细胞的初步定位。对于非黏连白细胞,直接使用阈值分割算法对白细胞浆进行分割。对于黏连白细胞,分割白细胞浆时,先用最大类间方差算法将白细胞背景和目标分成两部分,使错分概率最小,其次,利用Sobel求取梯度图像,根据获取的梯度幅值做分水岭变换,最终对获取的灰度图像做基于重建的开—闭运算,使去除图像的碎片效果更好,实现白细胞分割。The invention separates cell nuclei through a threshold segmentation algorithm, uses open operation to reconstruct and removes isolated points, and uses expansion processing to remove irrelevant areas, and the separation of cell nuclei realizes the preliminary positioning of white blood cells. For non-adherent leukocytes, the leukocyte plasma was segmented directly using a threshold segmentation algorithm. For adherent leukocytes, when segmenting leukocyte plasma, first use the maximum inter-class variance algorithm to divide the leukocyte background and the target into two parts to minimize the probability of misclassification. Second, use Sobel to obtain the gradient image, and perform watershed transformation according to the obtained gradient amplitude. , and finally perform the reconstruction-based open-close operation on the acquired grayscale image, so that the image fragment removal effect is better, and the white blood cell segmentation is realized.

本发明采用的技术方案包括如下步骤:The technical scheme adopted in the present invention comprises the following steps:

(1)获取100张白细胞图像数据集,将采集的显微图像由彩色图像转换为灰度图像;(1) Acquire 100 leukocyte image data sets, and convert the collected microscopic images from color images to grayscale images;

(2)将获取的灰度图像进行图像增强,用于突出图像的特征;(2) performing image enhancement on the acquired grayscale image to highlight the features of the image;

(3)将得到的图像增强后的图像用基于阈值分割的方法将白细胞的细胞核分离出来,实现白细胞的初步定位;(3) Separating the nuclei of leukocytes from the obtained image after image enhancement by a method based on threshold segmentation, so as to realize the preliminary positioning of leukocytes;

(4)去除获取白细胞细胞核图像的孤立点;(4) removing the isolated points where the image of leukocyte nuclei is obtained;

(5)去除细胞核图像孤立点图像无关区域,从而获取完整的白细胞细胞核,以上步骤同适应于黏连细胞和非黏连细胞提取细胞核过程;(5) Remove the irrelevant area of the isolated point image of the nucleus image, so as to obtain the complete leukocyte nucleus, the above steps are also suitable for the process of extracting the nucleus of the adherent cells and non-adherent cells;

(6)对于无黏连细胞的显微图像,根据白细胞核的初步定位和细胞面积最大,利用阈值分割的方法将白细胞浆直接分割出来;(6) For the microscopic images of non-adherent cells, according to the initial location of the leukocyte nucleus and the largest cell area, the leukocyte plasma is directly segmented by the threshold segmentation method;

(7)对于黏连细胞的显微图像,按图像的灰度特性,将图像分成背景和目标两部分,使错分概率最小;(7) For the microscopic image of the adherent cells, according to the grayscale characteristics of the image, the image is divided into two parts: background and target, so as to minimize the probability of misclassification;

(8)求取步骤7中获得的图像的梯度图像,同时将获取的梯度图像的梯度幅值做分水岭变换;(8) Obtain the gradient image of the image obtained in step 7, and simultaneously perform watershed transformation on the gradient magnitude of the obtained gradient image;

(9)将步骤8获得的灰度图像进行基于重建的开—闭运算,目的是在不影响目标图像的前提下,去除图像中的碎片效果更好,进而提取出黏连细胞的细胞浆,实现白细胞的分割。(9) Perform the reconstruction-based open-close operation on the grayscale image obtained in step 8, the purpose is to remove the debris in the image better without affecting the target image, and then extract the cytoplasm of the adherent cells, To achieve the segmentation of leukocytes.

所述的数据集选用同时混有背景、白细胞浆、红细胞和白细胞核的图像数据集。The data set is an image data set mixed with background, leukocyte plasma, erythrocytes and leukocyte nuclei at the same time.

所述步骤(1)具体为:Described step (1) is specifically:

(1.1)将获取的白细胞图像利用色彩转换的方式将彩色图像转换为灰色图像。(1.1) Convert the acquired white blood cell image into a gray image by means of color conversion.

所述步骤(2)具体为:Described step (2) is specifically:

(2.1)将步骤(1)中获取灰度图像使用线性变换的方法采用以下公式使图像增强,用于突出图像的特征:(2.1) The method of using linear transformation to obtain the grayscale image in step (1) adopts the following formula to enhance the image to highlight the features of the image:

G=F(g)=k·g+bG=F(g)=k·g+b

其中,当k>1时,代表输出图像对比度增大;当k<1时,代表输出图像对比度降低;当 k=1,b≠0时,代表仅使输出图像的灰度值上移或下移,其效果是使整个图像更亮或更暗。一般情况下,线性变换都是将某个较小的灰度范围拉伸到较大的灰度范围,因此常成称为灰度拉伸。设图像的灰度变化范围为[gmin,gmax],通过一个函数F,使范围扩展到 [Gmin,Gmax],则有线性变换:Among them, when k>1, it means that the contrast of the output image increases; when k<1, it means that the contrast of the output image decreases; when k=1, b≠0, it means that only the gray value of the output image is moved up or down Shift, which has the effect of making the entire image lighter or darker. In general, linear transformation stretches a smaller grayscale range to a larger grayscale range, so it is often called grayscale stretching. Let the grayscale variation range of the image be [g min , g max ], and through a function F to extend the range to [G min , G max ], there is a linear transformation:

即变换函数就是通过两个点(gmin,Gmin)和(gmax,Gmax)的直线方程。That is, the transformation function is the equation of a straight line passing through two points (g min , G min ) and (g max , G max ).

所述步骤(3)具体为:Described step (3) is specifically:

(3.1)将步骤(2)中得到的图像增强后的图像进行合适阈值的选取,用基于阈值分割的方法使用如下公式将白细胞的细胞核分离出来,实现白细胞的初步定位:(3.1) Select a suitable threshold for the image after image enhancement obtained in step (2), and use the method based on threshold segmentation to separate the nuclei of leukocytes using the following formula to achieve the initial positioning of leukocytes:

其中,T代表阈值,g(x,y)代表像素点(x,y)原来的灰度值,G(x,y)代表像素(x,y)变换后(增强后) 的灰度值。阈值是在分割时作为区分目标与背景像素的门限,大于或等于阈值的像素属于目标,而其他属于背景,反之亦然。Among them, T represents the threshold, g(x, y) represents the original gray value of the pixel (x, y), and G(x, y) represents the transformed (enhanced) gray value of the pixel (x, y). The threshold is used as the threshold to distinguish the target and background pixels during segmentation. The pixels greater than or equal to the threshold belong to the target, while the others belong to the background, and vice versa.

所述步骤(4)具体为:Described step (4) is specifically:

(4.1)将步骤(3)中提取出的细胞核用开运算重构的方法去除细胞核图像的孤立点。(4.1) The cell nucleus extracted in step (3) is reconstructed by open operation to remove the isolated points of the cell nucleus image.

所述步骤(5)具体为:Described step (5) is specifically:

(5.1)将步骤(4)中获得的去除孤立点图像,进行膨胀处理用如下公式去除无关区域:(5.1) The isolated point image obtained in step (4) is subjected to expansion processing to remove irrelevant regions with the following formula:

其中,表示膨胀,表示空集,B为结构元素。in, means expansion, Represents an empty set, and B is a structuring element.

所述步骤(6)具体为:Described step (6) is specifically:

(6.1)无白细胞黏连时,可根据细胞核的初步定位以及面积最大原理,通过阈值分割算法直接将白细胞浆分割出来。(6.1) When there is no leukocyte adhesion, the leukocyte plasma can be directly segmented by the threshold segmentation algorithm according to the initial location of the nucleus and the principle of maximum area.

所述步骤(7)具体为:Described step (7) is specifically:

(7.1)对于黏连细胞的显微图像,首先用最大类间方差算法,利用以下公式按图像的灰度特性,将图像分成背景和目标两部分,使错分概率最小:(7.1) For the microscopic image of adherent cells, first use the maximum inter-class variance algorithm, and use the following formula to divide the image into two parts: background and target according to the grayscale characteristics of the image, so as to minimize the probability of misclassification:

g=w0(μ0-μ)^2+w1(μ1-μ)^2g=w0(μ0-μ)^2+w1(μ1-μ)^2

其中,对于图像I(x,y),前景(及目标)和背景的分割阈值记作T,属于前景的像素点数占整幅图像的比例记为w0,其平均灰度μ0;背景像素点数占整幅图像的比例为w1,其平均灰度为μ1,图像的总平均灰度记为μ,类间方差记为g。Among them, for the image I(x, y), the segmentation threshold of the foreground (and the target) and the background is denoted as T, the proportion of the pixels belonging to the foreground to the entire image is denoted as w0, and the average gray level μ0; the number of background pixels accounts for The scale of the whole image is w1, its average gray level is μ1, the total average gray level of the image is recorded as μ, and the inter-class variance is recorded as g.

所述步骤(8)具体为:Described step (8) is specifically:

(8.1)将步骤(7)获得的图像利用Sobel求取梯度图像,同时对获取的梯度图像的梯度幅值做分水岭变换。(8.1) Use Sobel to obtain the gradient image of the image obtained in step (7), and perform watershed transformation on the gradient magnitude of the obtained gradient image.

所述步骤(9)具体为:Described step (9) is specifically:

(9.1)将步骤(8)获得的灰度图像进行基于重建的开—闭运算,去除图像中的碎片,进而提取出黏连细胞的细胞浆,实现白细胞的分割。(9.1) Perform an open-close operation based on reconstruction on the grayscale image obtained in step (8) to remove debris in the image, and then extract the cytoplasm of adherent cells to achieve leukocyte segmentation.

本发明的有益效果:Beneficial effects of the present invention:

本发明与传统白细胞定位方法,克服了以往的运算效率慢的特点,鲁棒性更强。在血液细胞分割方面具有良好的优势可以降低全自动血液分析仪的成本,辅助临床诊断,因此具有重要的临床意义和广阔的发展前景。The present invention and the traditional leukocyte localization method overcome the previous feature of slow operation efficiency and have stronger robustness. It has good advantages in blood cell segmentation, which can reduce the cost of automatic blood analyzers and assist clinical diagnosis, so it has important clinical significance and broad development prospects.

附图说明Description of drawings

本发明的白细胞图像的分割结果示意图:图1为非黏连细胞彩色图像,图8为黏连细胞彩色图像,图2为图1灰度图像结果图,图3为图2图像增强后的图像结果图,图4为图3阈值分割后的图像结果图,图5为图4开运算重构后的图像结果图,图6为图5膨胀处理后的细胞核图像结果图,图12为图8分割后细胞核图像结果图,图7为图1阈值分割后的非黏连细胞的白细胞浆图像结果图;图9为图8灰度图像的Sobel梯度算法后的图像结果图,图10为图9分水岭算法后的图像结果图,图11为图10重建开—闭运算后的图像结果图。Schematic diagram of the segmentation result of the white blood cell image of the present invention: Figure 1 is a color image of non-adherent cells, Figure 8 is a color image of adherent cells, Figure 2 is the result of the grayscale image of Figure 1, and Figure 3 is the image of Figure 2 after image enhancement Figure 4 is the result of the image after threshold segmentation in Figure 3, Figure 5 is the result of the image reconstructed by the open operation in Figure 4, Figure 6 is the result of the nucleus image after the expansion processing of Figure 5, and Figure 12 is Figure 8 Figure 7 is the result of the leukocyte plasma image of the non-adherent cells after threshold segmentation in Figure 1; Figure 9 is the image result of the grayscale image of Figure 8 after the Sobel gradient algorithm, and Figure 10 is Figure 9 The image result after the watershed algorithm, Figure 11 is the image result after the reconstruction of the open-close operation in Figure 10.

具体实施方式Detailed ways

下面结合附图和实施例对本发明进行进一步说明。The present invention will be further described below with reference to the accompanying drawings and embodiments.

本发明实施例如下:Examples of the present invention are as follows:

步骤1:选用100张同时混有背景、白细胞浆、红细胞和白细胞核的图像数据集,将采集的数据集图像由彩色图像转换为灰度图像。Step 1: Select 100 image datasets mixed with background, leukocyte plasma, erythrocytes and leukocyte nuclei, and convert the collected dataset images from color images to grayscale images.

(1.1)将获取的白细胞图像利用色彩转换的方式将彩色图像转换为灰色图像。(1.1) Convert the acquired white blood cell image into a gray image by means of color conversion.

步骤2:将获取的灰度图像进行图像增强,用于突出图像的特征。Step 2: Perform image enhancement on the acquired grayscale image to highlight the features of the image.

(2.1)将步骤1中获取灰度图像使用线性变换的方法采用以下公式使图像增强,用于突出图像的特征:(2.1) Using the method of linear transformation to obtain the grayscale image in step 1, the following formula is used to enhance the image to highlight the features of the image:

G=F(g)=k·g+bG=F(g)=k·g+b

其中,当k>1时,代表输出图像对比度增大;当k<1时,代表输出图像对比度降低;当 k=1,b≠0时,代表仅使输出图像的灰度值上移或下移,其效果是使整个图像更亮或更暗。设图像的灰度变化范围为[gmin,gmax],通过一个函数F,使范围扩展到[Gmin,Gmax],则有线性变换:Among them, when k>1, it means that the contrast of the output image increases; when k<1, it means that the contrast of the output image decreases; when k=1, b≠0, it means that only the gray value of the output image is moved up or down Shift, which has the effect of making the entire image lighter or darker. Let the grayscale variation range of the image be [g min , g max ], and through a function F to extend the range to [G min , G max ], there is a linear transformation:

即变换函数就是通过两个点(gmin,Gmin)和(gmax,Gmax)的直线方程。That is, the transformation function is the equation of a straight line passing through two points (g min , G min ) and (g max , G max ).

步骤3:将得到的图像增强后的图像用基于阈值分割的方法将白细胞的细胞核分离出来,实现白细胞的初步定位。Step 3: The obtained image after image enhancement is used to separate the nuclei of the leukocytes by a method based on threshold segmentation, so as to realize the preliminary positioning of the leukocytes.

(3.1)将步骤2中得到的图像增强后的图像进行合适阈值的选取,用基于阈值分割的方法使用如下公式将白细胞的细胞核分离出来,实现白细胞的初步定位:(3.1) Select the appropriate threshold for the enhanced image obtained in step 2, and use the threshold-based segmentation method to separate the nuclei of the leukocytes using the following formula to achieve the initial positioning of the leukocytes:

其中,T代表阈值,g(x,y)代表像素点(x,y)原来的灰度值,G(x,y)代表像素(x,y)变换后(增强后) 的灰度值。阈值是在分割时作为区分目标与背景像素的门限,大于或等于阈值的像素属于目标,而其他属于背景,反之亦然。Among them, T represents the threshold, g(x, y) represents the original gray value of the pixel (x, y), and G(x, y) represents the transformed (enhanced) gray value of the pixel (x, y). The threshold is used as the threshold to distinguish the target and background pixels during segmentation. The pixels greater than or equal to the threshold belong to the target, while the others belong to the background, and vice versa.

步骤4:去除获取白细胞细胞核图像的孤立点。Step 4: Remove isolated spots where images of leukocyte nuclei were obtained.

(4.1)将步骤3中提取出的细胞核用开运算重构的方法去除细胞核图像的孤立点。(4.1) The cell nucleus extracted in step 3 is reconstructed by open operation to remove the isolated points of the cell nucleus image.

步骤5:去除细胞核图像孤立点图像无关区域,从而获取完整的白细胞细胞核,以上步骤同适应于黏连细胞和非黏连细胞提取细胞核过程。Step 5: Remove the irrelevant area of the isolated point image of the nucleus image, so as to obtain the complete leukocyte nucleus. The above steps are also applicable to the process of extracting the nucleus of the adherent cells and non-adherent cells.

将步骤(4)中获得的去除孤立点图像,进行膨胀处理去除无关区域,公式为:The isolated point image obtained in step (4) is expanded to remove irrelevant regions, and the formula is:

其中,表示膨胀,表示空集,B为结构元素。in, means expansion, Represents an empty set, and B is a structuring element.

步骤6:对于无黏连细胞的显微图像,根据白细胞核的初步定位和细胞面积最大,利用阈值分割的方法将白细胞浆直接分割出来。Step 6: For the microscopic images of non-adherent cells, according to the initial location of the leukocyte nucleus and the largest cell area, the leukocyte plasma is directly segmented by the threshold segmentation method.

(6.1)无白细胞黏连时,可根据细胞核的初步定位以及面积最大原理,通过阈值分割算法直接将白细胞浆分割出来。(6.1) When there is no leukocyte adhesion, the leukocyte plasma can be directly segmented by the threshold segmentation algorithm according to the initial location of the nucleus and the principle of maximum area.

步骤7:对于黏连细胞的显微图像,按图像的灰度特性,将图像分成背景和目标两部分,使错分概率最小。Step 7: For the microscopic image of the adherent cells, according to the grayscale characteristics of the image, the image is divided into two parts, the background and the target, so as to minimize the probability of misclassification.

(7.1)对于黏连细胞的显微图像,首先用最大类间方差算法,利用以下公式按图像的灰度特性,将图像分成背景和目标两部分,使错分概率最小:(7.1) For the microscopic image of adherent cells, first use the maximum inter-class variance algorithm, and use the following formula to divide the image into two parts: background and target according to the grayscale characteristics of the image, so as to minimize the probability of misclassification:

g=w0(μ0-μ)^2+w1(μ1-μ)^2g=w0(μ0-μ)^2+w1(μ1-μ)^2

其中,对于图像I(x,y),前景(及目标)和背景的分割阈值记作T,属于前景的像素点数占整幅图像的比例记为w0,其平均灰度μ0;背景像素点数占整幅图像的比例为w1,其平均灰度为μ1,图像的总平均灰度记为μ,类间方差记为g。Among them, for the image I(x, y), the segmentation threshold of the foreground (and the target) and the background is denoted as T, the proportion of the pixels belonging to the foreground to the entire image is denoted as w0, and the average gray level μ0; the number of background pixels accounts for The scale of the whole image is w1, its average gray level is μ1, the total average gray level of the image is recorded as μ, and the inter-class variance is recorded as g.

步骤8:求取步骤7中获得的图像的梯度图像,同时将获取的梯度图像的梯度幅值做分水岭变换。Step 8: Obtain the gradient image of the image obtained in step 7, and perform watershed transformation on the gradient magnitude of the obtained gradient image.

步骤9:将步骤8获得的灰度图像进行基于重建的开—闭运算,去除图像中的碎片,进而提取出黏连细胞的细胞浆,实现白细胞的分割。Step 9: Perform the reconstruction-based open-close operation on the grayscale image obtained in Step 8 to remove debris in the image, and then extract the cytoplasm of the adherent cells to achieve the segmentation of white blood cells.

Claims (10)

1. A method for locating white blood cells in a lesion cell image is characterized by comprising the following steps:
(1) acquiring 100 white blood cell image data sets, selecting the image data sets mixed with background, white blood cell plasma, red blood cells and white cell nuclei at the same time, and converting the acquired microscopic images into gray images from color images;
(2) carrying out image enhancement on the acquired gray level image for highlighting the characteristics of the image;
(3) separating the cell nucleus of the white blood cell from the obtained image after the image enhancement by using a threshold segmentation method to realize the primary positioning of the white blood cell;
(4) removing isolated points for obtaining the image of the cell nucleus of the white blood cell;
(5) removing the independent area of the isolated point image of the cell nucleus image so as to obtain the complete cell nucleus of the white cell, wherein the steps are simultaneously suitable for the process of extracting the cell nucleus from the adherent cell and the non-adherent cell;
(6) for the microscopic image without the adherent cells, directly segmenting the white cell pulp by using a threshold segmentation method according to the primary positioning of white cell nuclei and the maximum cell area;
(7) for the microscopic image of the adhesion cells, dividing the image into a background part and a target part according to the gray characteristic of the image, so that the probability of wrong division is minimum;
(8) solving a gradient image of the image obtained in the step 7, and performing watershed transformation on the gradient amplitude of the obtained gradient image;
(9) and (3) carrying out reconstruction-based opening-closing operation on the gray level image obtained in the step (8), so that the fragment removing effect in the image is better on the premise of not influencing the target image, and then the cytoplasm with adhered cells is extracted, and the segmentation of white blood cells is realized.
2. The method of claim 1, wherein the method comprises: the step (1) is specifically as follows:
and (1.1) converting the acquired white blood cell image into a gray image by using a color conversion mode.
3. The method of claim 1, wherein the method comprises: the step (2) is specifically as follows:
(2.1) the method of obtaining the gray-scale image in the step (1) uses a linear transformation to enhance the image by adopting the following formula for highlighting the characteristics of the image:
wherein whenWhen, the output image contrast is increased; when in useWhen, it represents a decrease in output image contrast; when in useIn general, linear transformation is to stretch a small gray scale range to a large gray scale range, and is often called gray scale stretching, where the gray scale variation range of an image is set to beBy a functionExtend the range toThen, linear transformation:
i.e. the transformation function is passed through two pointsAndthe equation of the straight line of (c).
4. The method of claim 1, wherein the method comprises: the step (3) is specifically as follows:
(3.1) selecting a proper threshold value from the image obtained in the step (2) after image enhancement, and separating the cell nucleus of the white blood cell by using a threshold value segmentation-based method and using the following formula to realize the primary positioning of the white blood cell:
wherein, representing a threshold, G (x, y) represents the original gray value of the pixel (x, y), G (x, y) represents the gray value of the pixel (x, y) after transformation (after enhancement), the threshold is used as a threshold for distinguishing the target from the background pixel during segmentation, the pixel larger than or equal to the threshold belongs to the target, and the other pixels belong to the background, or vice versa.
5. The method of claim 1, wherein the method comprises: the step (4) is specifically as follows:
and (4.1) removing isolated points of the cell nucleus image by using a method of open operation reconstruction of the cell nucleus extracted in the step (3).
6. The method of claim 1, wherein the method comprises: the step (5) is specifically as follows:
(5.1) carrying out expansion processing on the isolated point removed image obtained in the step (4) to remove the irrelevant area by using the following formula:
in which, in the representation of swelling,representing an empty set, B being a structural element.
7. The method of claim 1, wherein the method comprises: the step (6) is specifically as follows:
(6.1) when no white blood cells are adhered, the white blood cell pulp can be directly segmented by a threshold segmentation algorithm according to the preliminary location and the area maximization principle of cell nuclei.
8. The method of claim 1, wherein the method comprises: the step (7) is specifically as follows:
(7.1) for the microscopic image of the adherent cells, firstly, dividing the image into a background part and a target part according to the gray level characteristics of the image by using a maximum between-class variance algorithm and the following formula so as to minimize the probability of wrong division:
wherein, for image I (x, y), the segmentation threshold of foreground (and target) and background is marked as T, the proportion of the number of pixels belonging to foreground in the whole image is marked as w0, and the average gray level is mu 0; the proportion of the number of background pixels in the whole image is w1, the average gray scale is μ 1, the total average gray scale of the image is μ, and the inter-class variance is g.
9. The method of claim 1, wherein the method comprises: the step (8) is specifically as follows:
and (8.1) solving a gradient image from the image obtained in the step (7) by using Sobel, and performing watershed transformation on the gradient amplitude of the obtained gradient image.
10. The method of claim 1, wherein the method comprises: the step (9) is specifically as follows:
and (9.1) performing on-off operation based on reconstruction on the gray image obtained in the step (8), removing fragments in the image, further extracting cytoplasm with adhered cells, and realizing the segmentation of the white blood cells.
CN201910856026.9A 2019-09-11 2019-09-11 A method for locating leukocytes in images of diseased cells Pending CN110610490A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910856026.9A CN110610490A (en) 2019-09-11 2019-09-11 A method for locating leukocytes in images of diseased cells

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910856026.9A CN110610490A (en) 2019-09-11 2019-09-11 A method for locating leukocytes in images of diseased cells

Publications (1)

Publication Number Publication Date
CN110610490A true CN110610490A (en) 2019-12-24

Family

ID=68892553

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910856026.9A Pending CN110610490A (en) 2019-09-11 2019-09-11 A method for locating leukocytes in images of diseased cells

Country Status (1)

Country Link
CN (1) CN110610490A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112070722A (en) * 2020-08-14 2020-12-11 厦门骁科码生物科技有限公司 A kind of fluorescence in situ hybridization cell nucleus segmentation method and system
CN119151977A (en) * 2024-11-18 2024-12-17 南京市浦口人民医院(江苏省人民医院浦口分院) Blood cell morphology contour obtaining method and system based on microscopic image

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104392460A (en) * 2014-12-12 2015-03-04 山东大学 Adherent white blood cell segmentation method based on nucleus-marked watershed transformation
CN105160668A (en) * 2015-08-26 2015-12-16 爱威科技股份有限公司 Image segmentation method and system, and cell image segmentation method and system

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104392460A (en) * 2014-12-12 2015-03-04 山东大学 Adherent white blood cell segmentation method based on nucleus-marked watershed transformation
WO2016091016A1 (en) * 2014-12-12 2016-06-16 山东大学 Nucleus marker watershed transformation-based method for splitting adhered white blood cells
CN105160668A (en) * 2015-08-26 2015-12-16 爱威科技股份有限公司 Image segmentation method and system, and cell image segmentation method and system

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
侯慧等: "基于距离变换的改进分水岭算法在白细胞图像分割中的应用", 《计算技术与自动化》 *
李舒等: "白细胞图像分割", 《辽宁科技大学学报》 *
缪慧司等: "基于HSI修正空间信息融合的彩色白细胞图像分割方法", 《生物医学工程学杂志》 *
胡彩亮: "基于指纹图的白细胞图像分割算法研究", 《中国优秀硕士学位论文全文数据库 信息科技辑(月刊)》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112070722A (en) * 2020-08-14 2020-12-11 厦门骁科码生物科技有限公司 A kind of fluorescence in situ hybridization cell nucleus segmentation method and system
CN119151977A (en) * 2024-11-18 2024-12-17 南京市浦口人民医院(江苏省人民医院浦口分院) Blood cell morphology contour obtaining method and system based on microscopic image

Similar Documents

Publication Publication Date Title
CN104484877B (en) A kind of AML cell segmentation methods clustered based on Meanshift with morphological operation
CN110675411B (en) Cervical squamous intraepithelial lesion recognition algorithm based on deep learning
CN110852350B (en) Pulmonary nodule benign and malignant classification method and system based on multi-scale migration learning
CN108537793A (en) A kind of pulmonary nodule detection method based on improved u-net networks
WO2016091016A1 (en) Nucleus marker watershed transformation-based method for splitting adhered white blood cells
Fan et al. Cancer survival prediction from whole slide images with self-supervised learning and slide consistency
Win et al. Automated segmentation of cell nuclei in cytology pleural fluid images using OTSU thresholding
WO2023045231A1 (en) Method and apparatus for facial nerve segmentation by decoupling and divide-and-conquer
Bai et al. Automatic segmentation of cervical region in colposcopic images using K-means
CN104933701B (en) The mammary glandular cell dividing method of adhesion model is removed with double strategies based on multiple dimensioned growth
CN103984958A (en) Method and system for segmenting cervical caner cells
CN111640128A (en) Cell image segmentation method based on U-Net network
CN103985119A (en) Method for partitioning cytoplasm and cell nucleuses of white blood cells in color blood cell image
CN110390678B (en) A tissue type segmentation method for colorectal cancer IHC stained images
CN104992435A (en) Cervix uteri single cell image segmentation algorithm
CN105761260A (en) Skin image affected part segmentation method
CN108320289A (en) A kind of bone marrow cell dividing method based on rarefaction representation and morphological operation
CN109033936A (en) A kind of cervical exfoliated cell core image-recognizing method
CN110544262B (en) A method of image segmentation of cervical cells based on machine vision
CN110610490A (en) A method for locating leukocytes in images of diseased cells
CN105574882B (en) Lung segmentation extracting method based on chest cross section CT images and system
CN117036288A (en) Tumor subtype diagnosis method for full-slice pathological image
CN110189327A (en) Fundus and retinal blood vessel segmentation method based on structured random forest encoder
CN112614093A (en) Breast pathology image classification method based on multi-scale space attention network
CN104933723A (en) Tongue image segmentation method based on sparse representation

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
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20191224