CN111724391B - Lymph cancer image fine segmentation method based on dynamic threshold - Google Patents

Lymph cancer image fine segmentation method based on dynamic threshold Download PDF

Info

Publication number
CN111724391B
CN111724391B CN202010425083.4A CN202010425083A CN111724391B CN 111724391 B CN111724391 B CN 111724391B CN 202010425083 A CN202010425083 A CN 202010425083A CN 111724391 B CN111724391 B CN 111724391B
Authority
CN
China
Prior art keywords
cluster
lymph cancer
suv
clusters
threshold value
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.)
Active
Application number
CN202010425083.4A
Other languages
Chinese (zh)
Other versions
CN111724391A (en
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.)
Zhejiang University of Technology ZJUT
Original Assignee
Zhejiang University of Technology ZJUT
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 Zhejiang University of Technology ZJUT filed Critical Zhejiang University of Technology ZJUT
Priority to CN202010425083.4A priority Critical patent/CN111724391B/en
Publication of CN111724391A publication Critical patent/CN111724391A/en
Application granted granted Critical
Publication of CN111724391B publication Critical patent/CN111724391B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection
    • G06T7/11Region-based segmentation
    • 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/136Segmentation; Edge detection involving thresholding
    • 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/20081Training; Learning
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30096Tumor; Lesion
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/10Internal combustion engine [ICE] based vehicles
    • Y02T10/40Engine management systems

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 Analysis (AREA)
  • Image Processing (AREA)

Abstract

A lymph cancer image fine segmentation method based on dynamic threshold, use the method of the double prediction precision to guarantee the accuracy that ROI chooses at first, and guarantee that the selective area includes the whole lymph cancer; secondly, regarding a single lymph cancer focus newly appeared in the optimization process as a new ROI area, and continuously optimizing the edge by using an independent threshold criterion so as to adjust the edge fine segmentation of the single lymph cancer focus in the dense area. According to the method, the initial segmentation result of the single focus image in the dense region is optimized by adopting a dynamic ROI threshold method, and a more accurate lymph cancer image is obtained by utilizing double-prediction precision information.

Description

Lymph cancer image fine segmentation method based on dynamic threshold
Technical Field
The invention relates to a lymph cancer image fine segmentation method based on a dynamic threshold value.
Technical Field
PET imaging of 18-fluorodeoxyglucose (18F-FDG) is a major means and tool for medical analysis of lymphoid cancers, and in PET images, standardized Uptake Value (SUV) is widely used to locate and segment lymphoid cancers. The PET image resolution of lymphoma is very low, and even an experienced expert, there are slight differences in manual identification of the same lymphoma, so in many cases, it is common to use 41% of the local maximum as a measure. However, using local thresholding as a segmentation criterion generally requires additional provision of regions of interest, and thresholding cannot be used in whole-body PET data. The diversity of lymphoma is clearly not satisfied by using the thresholding method directly. By being limited to a single maximum, surrounding images of lymphoid cancer masses are easily ignored, and some normal tissue images are mistaken for images of lymphoid cancer.
Disclosure of Invention
In order to overcome the defects of the prior art and solve the problem that the edge of a local focus image of some rough segmentation results is under-segmented, the invention provides a method for finely segmenting a lymph cancer image based on a dynamic threshold, and the initial segmentation result of a single focus image in a dense region is optimized by adopting a dynamic ROI threshold method.
In order to solve the above technical problems, the present invention can provide the following technical solutions:
a method for fine segmentation of a lymph cancer image based on a dynamic threshold, the method comprising the steps of:
step 1, obtaining an approximate position preparation stage, wherein probs refers to network output probability, T1 is a low threshold value, T2 is a high threshold value, and the process is as follows:
1.1 Obtaining a coordinate cluster C1 with lesion probability greater than T1 in probs;
1.2 C2 coordinates where the probability of foci in probs is greater than T2 were obtained.
1.3 For each C1, it is anded with C2, i.e., C1 ≈ C2, if not leaving C1 for the empty set, otherwise discarding.
Step 2, a local optimization stage, namely optimizing each cluster in the C1, wherein the cluster is a single focus set, and the process is as follows:
2.1 For a cluster of total pixel redundancy N1, if SUV minimum < SUV maximum 0.41: deleting the lowest value coordinate; otherwise, storing the cluster and performing the next cluster (the cluster is optimized);
2.2 For the total pixel points of the cluster less than the rest N1, storing the cluster and performing next cluster (the cluster is optimized);
2.3 Detect whether to classify into two clusters, delete the original cluster and add two clusters after splitting to C1 if split; not splitting: the next iteration continues to optimize the cluster;
2.4 Update the ROI;
2.5 ) repeat 2.1) -2.4) until C1 is
Figure BDA0002498400750000021
Further, in the step 1, the approximate position of the lymphoma is obtained, the low-precision of the lymphoma is indicated in a low threshold value T1 to ensure complete inclusion, and the high-precision of the lymphoma is indicated in a high threshold value T2 to ensure the accuracy of the focus, so that a good segmentation effect can be obtained in an undersegmented edge region to the maximum extent; excluding points with low SUV values until the lowest brightness value is greater than 41% of the highest brightness value; however, the remaining regions are monitored at regular intervals during the exclusion process, and if the region is split into two degree independent regions, there is no link between the two regions, and each region continues to optimize individually according to its own internal threshold, up to a 41% stopping criterion.
The technical conception of the invention is as follows: firstly, a double prediction precision method is used for ensuring the accuracy of ROI selection and ensuring that a selected region contains the whole lymph cancer; secondly, regarding a single lymph cancer focus newly appeared in the optimization process as a new ROI area, and continuously optimizing the edge by using an independent threshold criterion so as to adjust the edge fine segmentation of the single lymph cancer focus in the dense area.
The invention has the following beneficial effects: and acquiring a more accurate lymph cancer image by using the double prediction precision information.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
Detailed Description
The invention is further described below with reference to the accompanying drawings.
Referring to fig. 1, a method for performing a fine segmentation of a lymph cancer image based on a dynamic threshold value is a potential method for predicting a local lesion region on the basis of a rough lymph cancer region which can be obtained by a network prediction. Regions of high confidence are used as the primary criteria for candidate regions, while regions of low confidence are used primarily to include potential lesion regions as well. However, if there is only a low confidence prediction at one location, then this region will also be excluded. The focus position can be screened from the prediction result through the coincidence of high and low confidence degrees, and the omission is reduced.
Referring to fig. 1, the method for finely segmenting the lymph cancer image based on the dynamic threshold comprises the following steps:
step 1, obtaining an approximate position preparation stage, wherein probs refers to network output probability, T1 is a low threshold value, T2 is a high threshold value, and the process is as follows:
1.1 Obtaining a coordinate cluster C1 with lesion probability greater than T1 in probs;
1.2 C2 coordinates where the probability of foci in probs is greater than T2 were obtained.
1.3 For each C1, it is anded with C2, i.e., C1 ≈ C2, if not leaving C1 for the empty set, otherwise discarding.
Step 2, a local optimization stage, namely optimizing each cluster in the C1, wherein the cluster is a single focus set, and the process is as follows:
2.1 For a cluster of total pixel redundancy N1, if SUV minimum < SUV maximum 0.41: deleting the lowest value coordinate; otherwise, storing the cluster and performing next cluster (the cluster is optimized);
2.2 For the total pixel points of the cluster less than the rest N1, storing the cluster and performing next cluster (the cluster is optimized);
2.3 Detecting whether the cluster is classified into two clusters, and if the cluster is split, deleting the original cluster and adding the two clusters after splitting into C1; not splitting: the next iteration continues to optimize the cluster;
2.4 Update the ROI;
2.5 ) repeat 2.1) -2.4) until C1 is
Figure BDA0002498400750000031
Further, in the step 1, the approximate position of the lymphoma is obtained, the low-precision of the lymphoma is indicated in a low threshold value T1 to ensure complete inclusion, and the high-precision of the lymphoma is indicated in a high threshold value T2 to ensure the accuracy of the focus, so that a good segmentation effect can be obtained in an undersegmented edge region to the maximum extent; excluding points with low SUV values until the lowest brightness value is greater than 41% of the highest brightness value; however, the remaining regions are monitored at regular intervals during the exclusion process, and if the region is split into two degree independent regions, there is no link between the two regions, and each region continues to optimize individually according to its own internal threshold, up to a 41% stopping criterion.
In the scheme of this embodiment, in step 1, a candidate region is determined by two confidence levels; in the step 2, the ROI is updated through continuous iteration, so that different local thresholds are found to be more suitable for complex case characteristics of the lymph cancer. Thus, the good segmentation effect can be obtained in the marginal area of the lymph cancer undersegmented to the maximum extent.

Claims (1)

1. A lymph cancer image fine segmentation method based on dynamic threshold is characterized by comprising the following steps:
step 1, obtaining an approximate position preparation stage, wherein probs refers to network output probability, T1 is a low threshold value, T2 is a high threshold value, and the process is as follows:
1.1 Obtaining a coordinate cluster C1 with lesion probability greater than T1 in probs;
1.2 Obtaining a coordinate C2 with a lesion probability in probs greater than T2;
1.3 For each C1, make it and C2 phase, namely C1 ^ C2, if not for the empty set to reserve C1, otherwise abandon;
step 2, a local optimization stage, namely optimizing each cluster in the C1, wherein the cluster is a single focus set, and the process is as follows:
2.1 For clusters with more than N1 total pixels, if SUV minimum < SUV maximum 0.41: deleting the lowest value coordinate; otherwise, storing the cluster and carrying out the next cluster;
2.2 For the total pixel points of the cluster are less than N1, storing the cluster and carrying out the next cluster;
2.3 Detecting whether the cluster is classified into two clusters, and if the cluster is split, deleting the original cluster and adding the two clusters after splitting into C1; not splitting: the next iteration continues to optimize the cluster;
2.4 Update the ROI;
2.5 ) repeat 2.1) -2.4) until C1 is
Figure FDA0004058155650000011
In the step 1, the approximate position of the lymph cancer is obtained, the low precision of the lymph cancer is indicated in a low threshold value T1 to ensure the completeness of the lymph cancer, and the high precision is indicated in a high threshold value T2 to ensure the accuracy of the focus, so that a good segmentation effect can be obtained in an undersegmented edge region to the maximum extent; in the step 2, points with low SUV values are excluded until the SUV lowest value is greater than 41% of the SUV maximum value; however, the remaining clusters are monitored at regular intervals during the exclusion process, if the cluster is split into two independent clusters, the two clusters are not linked to each other, and each cluster is individually optimized according to its own internal stopping criteria until the SUV minimum in the cluster is greater than 41% of the SUV maximum.
CN202010425083.4A 2020-05-19 2020-05-19 Lymph cancer image fine segmentation method based on dynamic threshold Active CN111724391B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010425083.4A CN111724391B (en) 2020-05-19 2020-05-19 Lymph cancer image fine segmentation method based on dynamic threshold

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010425083.4A CN111724391B (en) 2020-05-19 2020-05-19 Lymph cancer image fine segmentation method based on dynamic threshold

Publications (2)

Publication Number Publication Date
CN111724391A CN111724391A (en) 2020-09-29
CN111724391B true CN111724391B (en) 2023-04-07

Family

ID=72564622

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010425083.4A Active CN111724391B (en) 2020-05-19 2020-05-19 Lymph cancer image fine segmentation method based on dynamic threshold

Country Status (1)

Country Link
CN (1) CN111724391B (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104463862A (en) * 2014-11-28 2015-03-25 哈尔滨工业大学 Method for fast segmenting kidney CT sequential image
CN106997596A (en) * 2017-04-01 2017-08-01 太原理工大学 A kind of Lung neoplasm dividing method of the LBF movable contour models based on comentropy and joint vector
CN107610145A (en) * 2017-07-26 2018-01-19 同济大学 A kind of automatic pancreas dividing method based on adaptive threshold and template matches
CN110197474A (en) * 2018-03-27 2019-09-03 腾讯科技(深圳)有限公司 The training method of image processing method and device and neural network model

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010027476A1 (en) * 2008-09-03 2010-03-11 Rutgers, The State University Of New Jersey System and method for accurate and rapid identification of diseased regions on biological images with applications to disease diagnosis and prognosis

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104463862A (en) * 2014-11-28 2015-03-25 哈尔滨工业大学 Method for fast segmenting kidney CT sequential image
CN106997596A (en) * 2017-04-01 2017-08-01 太原理工大学 A kind of Lung neoplasm dividing method of the LBF movable contour models based on comentropy and joint vector
CN107610145A (en) * 2017-07-26 2018-01-19 同济大学 A kind of automatic pancreas dividing method based on adaptive threshold and template matches
CN110197474A (en) * 2018-03-27 2019-09-03 腾讯科技(深圳)有限公司 The training method of image processing method and device and neural network model

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
《A Background-based Data Enhancement Method for Lymphoma》;Haigen Hu;《2019 IEEE International Conference on Bioinformatics and Biomedicine (BIBM)》;1194-1196 *
《A Prior Knowledge Intergrated Scheme for Detection and Segmentation of Lymphomas in 3D PET Images based on DBSCAN and GAs》;Haigen Hu;《2019 IEEE International Conference on Bioinformatics and Biomedicine (BIBM)》;2413-2420 *
《基于图像分割的淋巴癌细胞提取方法研究》;张佳佳;《中国优秀硕士学位论文全文数据库信息科技辑》;第2014年卷(第06期);I138-898 *
《基于自适应滤波的淋巴细胞图像分割算法研究》;杨建菊;《计算机仿真》;第29卷(第8期);265-268 *
PET-CT用于评价食管鳞癌放疗中~(18)F-FDG高摄取区域的空间动态变化的前瞻性研究;刘琪等;《中国癌症杂志》(第02期);161-167 *

Also Published As

Publication number Publication date
CN111724391A (en) 2020-09-29

Similar Documents

Publication Publication Date Title
US9311698B2 (en) Detecting defects on a wafer using template image matching
CN111445478B (en) Automatic intracranial aneurysm region detection system and detection method for CTA image
JP5981838B2 (en) System, method and computer program product for detecting defects in inspection images
CN109461148A (en) Steel rail defect based on two-dimentional Otsu divides adaptive fast algorithm
US20220122252A1 (en) Techniques for determining tissue characteristics using multiplexed immunofluorescence imaging
CN110853005A (en) Immunohistochemical membrane staining section diagnosis method and device
CN112561849B (en) Wafer defect detection method
US8019156B2 (en) Method and device for segmenting a digital representation of cells
CN112215790A (en) KI67 index analysis method based on deep learning
CN109416749B (en) Image gray scale classification method and device and readable storage medium
CN104657980A (en) Improved multi-channel image partitioning algorithm based on Meanshift
CN112037221B (en) Multi-domain co-adaptation training method for cervical cancer TCT slice positive cell detection model
JP6733983B2 (en) Image analysis device
CN113439276A (en) Defect classification in semiconductor samples based on machine learning
CN116740072B (en) Road surface defect detection method and system based on machine vision
CN115170518A (en) Cell detection method and system based on deep learning and machine vision
US11270439B2 (en) Histogram-based method for auto segmentation of integrated circuit structures from SEM images
CN113253269B (en) SAR self-focusing method based on image classification
CN111724391B (en) Lymph cancer image fine segmentation method based on dynamic threshold
CN117911408A (en) Road pavement construction quality detection method and system
Triyani et al. Performance analysis of image segmentation for breast ultrasound images
CN117853722A (en) Steel metallographic structure segmentation method integrating superpixel information
CN115511815A (en) Cervical fluid-based cell segmentation method and system based on watershed
CN117995317B (en) Method, device and medium for estimating heavy atom position based on electron density map
CN111833361A (en) Tool and method for segmenting target size

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
GR01 Patent grant
GR01 Patent grant