WO2020093563A1 - Medical image processing method, system, device, and storage medium - Google Patents

Medical image processing method, system, device, and storage medium Download PDF

Info

Publication number
WO2020093563A1
WO2020093563A1 PCT/CN2018/124660 CN2018124660W WO2020093563A1 WO 2020093563 A1 WO2020093563 A1 WO 2020093563A1 CN 2018124660 W CN2018124660 W CN 2018124660W WO 2020093563 A1 WO2020093563 A1 WO 2020093563A1
Authority
WO
WIPO (PCT)
Prior art keywords
image
lesion
medical image
abnormality
abnormal
Prior art date
Application number
PCT/CN2018/124660
Other languages
French (fr)
Chinese (zh)
Inventor
徐勇
罗笑玲
蒲祖辉
牟丽莎
胡吉英
Original Assignee
哈尔滨工业大学(深圳)
深圳市第二人民医院
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 哈尔滨工业大学(深圳), 深圳市第二人民医院 filed Critical 哈尔滨工业大学(深圳)
Publication of WO2020093563A1 publication Critical patent/WO2020093563A1/en

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/24Classification techniques
    • 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
    • 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/20Special algorithmic details
    • G06T2207/20084Artificial neural networks [ANN]
    • 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/30041Eye; Retina; Ophthalmic
    • 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
    • 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/30101Blood vessel; Artery; Vein; Vascular

Definitions

  • the invention relates to the field of image processing, in particular to a medical image processing method, system, equipment and storage medium.
  • diabetes is considered to be the basis of many health problems and late-stage disorders, that is, diabetes can cause a series of diseases and complications, including: leading to serious heart disease, diabetic retinopathy (DR) and kidney problems.
  • DR diabetic retinopathy
  • DR detection technology is particularly important, and medical image analysis is one of the research fields that have attracted great interest from scientists and physicians.
  • DR can be divided into two major categories, namely non-proliferative diabetic retinopathy (NPDR) and proliferative retinopathy (PDR). There are three subcategories of NPDR: mild NPDR, moderate NPDR, and severe NPDR.
  • the non-proliferative type is the early stage of the lesion, and the lesion is limited to the retina, manifested as microhemangioma, hemorrhage, hard and soft exudate, retinal artery and venous lesions.
  • Proliferative lesions at least partially extend beyond the inner limiting membrane, and the appearance of neovascularization is a sign of proliferative.
  • DR lesions are described in detail as follows:
  • Microaneurysms represent the most primitive perceptible signs of retinal damage.
  • the abnormal permeability of retinal blood vessels leads to the formation of microaneurysms.
  • the microaneurysm can be regarded as a red dot with sharp edges and a size between 20 ⁇ m and 200 ⁇ m, which is approximately 8.25% of the size of the disc.
  • Rigid exudates are formed by the leakage of lipoproteins and other proteins from the blood vessels of the retina. Visually, it looks like small white or yellowish-white deposits with distinct edges. Rigid exudates are usually organized in the form of a ring, usually appearing in the outer layer of the retina. Rigid exudates are usually irregular and shiny, and are found close to the edges of microaneurysms or retinal edema.
  • RNFL retinal nerve fiber layer
  • Bleeding occurs due to leakage of weak blood vessels.
  • the focus of bleeding is in the form of red dots with different densities and uneven edges, and it is found in the range of 125 ⁇ m.
  • hemorrhage is divided into two categories: flame and dot blotting.
  • the first type originates from the anterior capillary artery and appears on nerve fibers.
  • the second type of dot blot bleeding is round and smaller than a microaneurysm. Dot blot bleeding can appear on the retina at different levels, however, in most cases it will appear at the end of the capillaries' meridians.
  • Neovascularization usually indicates the atypical appearance of new blood vessels that appear on the inner surface of the retina.
  • the new blood vessels are small and repeatedly penetrate into the vitreous cavity, which reduces visual ability and makes it significantly blurry, which eventually leads to blindness.
  • Macular edema is identified as a swollen part of the retina, which usually occurs due to the permeability of abnormal retinal capillaries. Macular edema causes leakage of fluid or other solutes around the macula and severely affects vision.
  • an object of the present invention is to provide a medical image processing method, system, device, and storage medium for improving the efficiency of processing and analyzing medical images.
  • the present invention provides a medical image processing method, including the following steps:
  • the lesion image acquisition step the lesion image is extracted based on the medical image
  • the lesion score is a first abnormal probability that the lesion image belongs to different image abnormalities
  • a second abnormal probability acquisition step acquiring multiple second abnormal probabilities of the medical image according to the medical image and the lesion image, the second abnormal probability being the probability that the medical image belongs to different image abnormalities;
  • the classification step according to a plurality of the first abnormal probability, a plurality of the second abnormal probability and different weighting coefficients, the final probability that the medical image belongs to different image abnormalities is obtained according to different image abnormalities, and the largest The degree of image abnormality of the final probability is taken as the degree of image abnormality of the medical image.
  • the medical image includes a fundus photograph.
  • the step of acquiring a lesion image includes:
  • step of generating scores includes:
  • the step of acquiring the second abnormal probability includes:
  • the degree of image abnormality includes normal, mild first abnormality, moderate first abnormality, severe first abnormality, and second abnormality.
  • the weight coefficient of the first abnormal probability is 0.2
  • the weight coefficient of the second abnormal probability is 0.8.
  • the present invention provides a medical image processing system, including:
  • the lesion image acquisition unit is used to extract the lesion image according to the medical image
  • a score generating unit configured to obtain a plurality of lesion scores based on the lesion image, where the lesion score is a first abnormal probability that the lesion image belongs to different image abnormalities;
  • a second abnormal probability obtaining unit configured to obtain a plurality of second abnormal probabilities of the medical image according to the medical image and the lesion image, the second abnormal probability being the probability that the medical image belongs to different image abnormalities ;
  • the classification unit is used to obtain the final probability that the medical image belongs to different image abnormalities according to different image abnormalities according to multiple first abnormal probabilities, multiple second abnormal probabilities and different weighting coefficients
  • the image abnormality of the final probability is taken as the image abnormality of the medical image.
  • the present invention provides a medical image processing device, including:
  • At least one processor and,
  • a memory communicatively connected to the at least one processor; wherein,
  • the memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor to enable the at least one processor to execute the medical image processing method.
  • the present invention provides a computer-readable storage medium that stores computer-executable instructions that are used to cause a computer to execute the medical image processing method.
  • the present invention obtains multiple first abnormal probabilities through the focus image, on the other hand obtains multiple second abnormal probabilities through the medical image and focus image, and then according to multiple first abnormal probabilities, multiple second abnormal probabilities and different
  • the weight coefficients are used to obtain the final probability of medical images belonging to different image abnormalities according to different image abnormalities, and the image abnormality of the maximum final probability is taken as the image abnormality of medical images to realize the analysis of medical image abnormality and overcome the existing
  • the present invention also obtains multiple lesion mask images of medical images by segmenting a neural network to further obtain multiple lesion images, ensuring accurate extraction of lesion images.
  • FIG. 1 is a flowchart of a specific embodiment of a medical image processing method in the present invention.
  • FIG. 2 is a schematic diagram of a specific embodiment of a segmented neural network of a medical image processing method in the present invention.
  • a medical image processing method includes the following steps:
  • the lesion image acquisition step the lesion image is extracted based on the medical image
  • the lesion score is the first abnormal probability that the lesion image belongs to different image abnormalities. For example, if the image abnormality includes both normal and abnormal, the score generation step is to obtain the focus image. Normal or abnormal probability, get two first abnormal probabilities, and so on;
  • the second abnormal probability obtaining step multiple second abnormal probabilities of the medical image are obtained according to the medical image and the lesion image, and the second abnormal probability is the probability that the medical image belongs to different image abnormalities;
  • the classification step is to obtain the final probabilities of medical images belonging to different image abnormalities according to different image abnormalities according to multiple first abnormal probabilities, multiple second abnormal probabilities and different weighting coefficients, and take the maximum final probability image abnormality as the The degree of abnormality of medical images. That is, the first abnormal probability and the second abnormal probability of the same image abnormal degree are added according to different weight coefficients to obtain the final probability of the image abnormal degree. Since there are multiple different image abnormal degrees, multiple final probabilities of the medical image can be obtained .
  • a medical image processing method On the one hand, multiple first abnormal probabilities are obtained from focus images, on the other hand, multiple second abnormal probabilities are obtained from medical images and focus images, and then multiple first abnormal probabilities and multiple second Abnormal probability and different weight coefficients are used to obtain the final probability that medical images belong to different image abnormalities according to different image abnormalities, and the image abnormality of the maximum final probability is taken as the image abnormality of medical images to realize the abnormality of medical images.
  • the analysis improves the processing and analysis efficiency of medical images, and overcomes the existing technical problems of low efficiency in relying on the naked eye for processing and analyzing pathological images.
  • medical images include, but are not limited to, fundus photos.
  • the fundus photos are used as examples to describe the medical image processing method.
  • the degree of image abnormality includes but is not limited to normal, mild first abnormality, moderate first abnormality, severe first abnormality, and second abnormality. Taking the fundus photo as an example, the first abnormality is a non-proliferative lesion. The second abnormality is a proliferative lesion.
  • FIG. 1 is a flowchart of a specific embodiment of a medical image processing method according to the present invention; the medical image processing method will be specifically described below:
  • Step 1 Pre-process the data.
  • the data used in the present invention is a two-dimensional color fundus photograph.
  • the original picture size is 3000 ⁇ 3000.
  • the picture needs to be pre-processed, and the background and other unnecessary parts of the picture are cropped to adjust the picture size to 128 ⁇ 128.
  • Step 2 The step of acquiring lesion images, including:
  • the main task of segmenting the neural network is to realize the segmentation and extraction of multiple lesions.
  • the types of lesions include microaneurysms, rigid exudates, and neonatal retinal blood vessels. Lesions, new retinal blood vessels are proliferative lesions.
  • the three masks are the micro-aneurysm mask, the hard exudate mask, and the retinal vascular mask.
  • the mask is specifically composed of 0 and 1, the background is 0, and the target is 1.
  • micro-aneurysm mask when generating a micro-aneurysm mask, predict whether each pixel belongs to the micro-aneurysm area by segmenting the neural network. If the pixel position is marked as 1, Otherwise, the pixel position is marked as 0, and finally a micro-aneurysm mask image with the same size as the original image of 128 ⁇ 128 is obtained.
  • the mask can accurately describe the location of the target. Corresponding the mask to the original image, you can extract the image of each lesion, which is the lesion image.
  • FIG. 2 is a schematic diagram of a specific embodiment of a segmentation neural network of a medical image processing method in the present invention
  • the implementation of the segmentation neural network is mainly based on a multi-task learning framework, and there are three segmentation tasks of the segmentation neural network , Respectively, segmentation of microaneurysm, segmentation of hard exudate, and segmentation of retinal vessels.
  • the first part of the network is called the "shared layer", and the weights of these layers are common to all three tasks.
  • the parameters between these layers are calculated independently and do not participate in cross-layer sharing. In these task-specific layers, the network learns task-specific information.
  • Each of these independent task layers will generate a separate output, specifically a segmentation mask that outputs microaneurysms, hard exudates, and retinal blood vessels.
  • the above segmentation neural network can segment multiple lesion images simultaneously, and is not limited to the three lesion images in this embodiment.
  • the network Before using the segmented neural network, the network needs to be trained. Specifically, the model is first initialized, and then the sample with the marked mask data is input for training. The entire picture is used for training without patch sampling. The experiment proves that it is more effective to use the whole picture.
  • N is the number of training samples
  • L MA is the pixel-level loss of micro-aneurysm segmentation
  • L HES is the segmentation loss of hard exudate
  • L NV is the segmentation loss of retinal vessels
  • f (x, ⁇ ) is the segmentation The segmentation result predicted by the neural network, where x is a pixel of the sample and ⁇ is the learning rate
  • S 1 , S 2 , and S 3 are the labeled positions of the samples of microaneurysm, hard exudate, and retinal blood vessel
  • ⁇ ( ⁇ ) Is a regular item.
  • the trained segmentation neural network uses the trained segmentation neural network to obtain a variety of lesion mask images of the fundus photos, including micro-aneurysm mask images, hard exudate mask images, and retinal vascular mask images; according to the lesion mask images and medicine
  • the image can be segmented to obtain a variety of lesion images, including micro-aneurysm lesion images, rigid exudate lesion images, and retinal vascular lesion images.
  • Step 3 The score generation step, including:
  • the features of multiple lesions are spliced and input into the first machine learning classification algorithm to obtain multiple lesion scores.
  • the features such as the color and shape of each lesion are extracted according to the lesion image as shown in Table 2, 86-dimensional features can be extracted from each lesion image, and the lesion features of the 3 lesion images are The sequence of micro-aneurysm lesion images, rigid exudate lesion images, and retinal vascular lesion images are stitched together into a 258-dimensional feature vector, and the feature vector is input into the score generator (that is, the first machine learning classification algorithm) to obtain multiple
  • the score of each lesion is i is the category of the abnormality of the image, X is the number of the medical image, that is, the acquired lesion image belongs to different image abnormality (including normal, mild non-proliferative lesions, moderate non-proliferative lesions, severe non-proliferative lesions and proliferative Disease), that is, the score of 5 lesions can be obtained.
  • the score generator uses a random forest algorithm.
  • the score generator needs to be trained before it is used. Specifically, the 258-dimensional feature vector of the lesion feature of the training sample and the degree of lesion marked by the sample are input into the score generator for training.
  • the trained score generator can calculate the probability that the input lesion belongs to 5 image abnormalities (that is, normal, mild non-proliferative lesions, moderate non-proliferative lesions, severe non-proliferative lesions, and proliferative lesions) , Take this probability value as the lesion score of the sample.
  • Step 4 The second abnormal probability acquisition step, including:
  • the image abnormality includes normal, mild non-proliferative lesions, moderate non-proliferative lesions, severe non-proliferative lesions, and proliferative lesions;
  • the second machine learning classification algorithm is designed based on 3D convolutional neural network
  • the 3D convolutional neural network can simultaneously input the original medical images (ie, fundus photos), microaneurysm lesion images, rigid exudate lesion images, and retinal vascular lesion images, which can extract features from multiple dimensions and then perform 3D convolution To capture feature information from multiple images.
  • the 3D convolutional neural network includes a solid-wired layer (ie, hard-connected layer), 4 convolutional layers, 3 down-sampled layers, and a fully connected layer.
  • the convolved cubes of each convolution kernel are 4 pictures including the original image, the image of the microhemangioma lesion, the image of the retinal vascular lesion, and the image of the hard exudate lesion.
  • the size of each picture is 128 ⁇ 128.
  • a fixed solid-line kernel is applied to process the input graph to generate multiple channels of information, then multiple channels are processed separately, and finally the information of all channels are combined to obtain the final Feature description.
  • This solid layer actually encodes a priori knowledge of the lesion, which is better than random initialization performance.
  • L 3 layer which is down-sampled with a 3 ⁇ 3 window in the feature map of the L 2 layer, so that the same number of lesion maps with reduced spatial resolution will be obtained.
  • the L 5 layer uses a 2 ⁇ 2 downsampling window, followed by two 2D convolutions and downsampling to obtain a 128 ⁇ 1 ⁇ 1 feature vector, and the 128 dimensions are determined based on past experience.
  • the input multiple images After inputting medical images, after multiple layers of convolution and downsampling, the input multiple images are converted into a 128-dimensional feature vector, which captures the feature information of the retinal fundus photo. Then the 128-dimensional feature vector obtained is input to the softmax layer.
  • the number of nodes in the softmax layer is consistent with the number of image abnormality categories, and each node is fully connected to the 128 nodes in L 9 .
  • the 3D convolutional neural network Before using the 3D convolutional neural network, you need to train it, first initialize the model, and then input the sample marked with abnormal degree for training. Input the original image of the sample, the image of the microhemangioma lesion, and the image of the hard exudate lesion during training. 4 pictures of retinal vascular lesion images for training.
  • the 3D convolutional neural network model is continuously optimized during the training process, and the parameters are adjusted until convergence. The goal of optimization is to continuously reduce the difference between the classification results predicted by the network and the abnormal degree of sample labeling, that is, to minimize the loss function L2:
  • N is the number of training samples
  • L softmax represents the network classification loss
  • f (x, ⁇ ) is the classification result predicted by the network, where x is a sample, ⁇ is the learning rate
  • C is the label category
  • ⁇ ( ⁇ ) is the regular term.
  • a medical image processing method a score fusion mechanism is proposed, and the lesion scores obtained from each lesion image are integrated into a classification network, which improves the sensitivity and classification accuracy of each lesion.
  • the score fusion classification framework of the present invention is not limited to the three lesions detected in the present invention, and can be expanded when more lesions are detected, and the system has strong adaptability.
  • a medical image processing system including:
  • the lesion image acquisition unit is used to extract the lesion image according to the medical image
  • the score generating unit is used to obtain multiple lesion scores according to the lesion image, and the lesion score is the first abnormal probability that the lesion image belongs to different image abnormalities;
  • a second abnormal probability acquisition unit configured to acquire multiple second abnormal probabilities of the medical image according to the medical image and the lesion image, the second abnormal probability is the probability that the medical image belongs to different image abnormalities;
  • the classification unit is used to obtain the final probabilities that medical images belong to different image abnormalities according to different image abnormalities according to multiple first abnormal probabilities, multiple second abnormal probabilities and different weight coefficients, and to abnormalize the images with the largest final probability
  • the degree is regarded as the abnormal degree of the medical image.
  • a medical image processing device including:
  • At least one processor and,
  • a memory communicatively connected to the at least one processor; wherein,
  • the memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor to enable the at least one processor to execute the medical image processing method.
  • the medical image processing method refer to the description in Embodiment 1, and no further description is required.
  • a computer-readable storage medium stores computer-executable instructions for causing a computer to execute the medical image processing method.
  • the medical image processing method refer to the description in Embodiment 1, and no further description is required.

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Quality & Reliability (AREA)
  • Artificial Intelligence (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medical Informatics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Image Analysis (AREA)
  • Eye Examination Apparatus (AREA)
  • Image Processing (AREA)

Abstract

Disclosed in the present invention are a medical image processing method, a system, a device, and a storage medium. In one aspect, a plurality of first abnormality probabilities are acquired by means of a lesion image, and in another aspect, a plurality of second abnormality probabilities are acquired by means of a medical image and the lesion image. Then, on the basis of the plurality of first abnormality probabilities, the plurality of second abnormality probabilities, and different weight coefficients, and according to different image abnormality degrees, a final probability that the medical image belongs to the different image abnormality degrees is acquired, and the image abnormality degree having the greatest final probability is set as the image abnormality degree of the medical image, thus achieving analysis of the abnormality degree of the medical image, and overcoming the problems in the prior art of reliance on the naked eye and low efficiency in pathological image processing and analysis.

Description

一种医学图像处理方法、系统、设备、存储介质Medical image processing method, system, equipment and storage medium 技术领域Technical field
本发明涉及图像处理领域,尤其是一种医学图像处理方法、系统、设备、存储介质。The invention relates to the field of image processing, in particular to a medical image processing method, system, equipment and storage medium.
背景技术Background technique
过去的几十年期间糖尿病的高速增长引起了各界的关注,糖尿病引起的疾病的指数增长已成为当前医疗保健行业面临的巨大挑战。而且不幸的是,患有糖尿病引起的疾病的患者数量仍然以惊人的速度持续增长。还有更让人担忧的是,仅有大约70%的患者意识到他们患有这种疾病。从医学角度来看,糖尿病被认为是许多健康问题和后期障碍的基础,即糖尿病会引起一系列的病变和并发症,包括:导致严重心脏病,糖尿病性视网膜病变(DR)和肾脏问题等。DR是最常见的糖尿病并发症之一,它被认为是失明的最主要原因之一。研究表明,无论人口数量或社会经济背景如何,世界各个地区的糖尿病问题都在以惊人的速度增长。此外,一项研究显示,由于发展中国家的生活条件和治疗设施存在缺陷,所以近75%的DR患者属于发展中国家,而且DR患者的失明可能性比没有患有DR的人几乎高25倍。但是当前缺少糖尿病视网膜病变的医学专家,世界上大部分的DR患者不能及时得到有效的病变检测和治疗,大部分患者只有当视网膜病变已经发展到治疗变得高度复杂且有时几乎不可行的程度时才能意识到寻求治疗。然而,DR在初始阶段治愈率可以达到90%,因此及早发现DR并且得到有效治疗能都大大降低DR导致失明的风险。所以,DR检测技术显得尤为重要,医学图像分析是目前引起科学家和医师极大兴趣的研究领域之一。The rapid growth of diabetes during the past few decades has attracted attention from all walks of life. The exponential growth of diseases caused by diabetes has become a huge challenge facing the current healthcare industry. And unfortunately, the number of patients with diseases caused by diabetes continues to grow at an alarming rate. What is more worrying is that only about 70% of patients realize that they have this disease. From a medical point of view, diabetes is considered to be the basis of many health problems and late-stage disorders, that is, diabetes can cause a series of diseases and complications, including: leading to serious heart disease, diabetic retinopathy (DR) and kidney problems. DR is one of the most common complications of diabetes, and it is considered to be one of the most important causes of blindness. Studies have shown that regardless of population size or socioeconomic background, diabetes problems in all regions of the world are increasing at an alarming rate. In addition, a study showed that due to defects in living conditions and treatment facilities in developing countries, nearly 75% of DR patients belong to developing countries, and the probability of blindness of DR patients is almost 25 times higher than those without DR . However, there is currently a lack of medical experts for diabetic retinopathy. Most DR patients in the world cannot get effective lesion detection and treatment in a timely manner. Most patients only develop when retinopathy has developed to a point where the treatment becomes highly complex and sometimes almost impossible. To be aware of seeking treatment. However, the cure rate of DR in the initial stage can reach 90%, so early detection of DR and effective treatment can greatly reduce the risk of DR causing blindness. Therefore, DR detection technology is particularly important, and medical image analysis is one of the research fields that have attracted great interest from scientists and physicians.
DR可以分为两大类即非增殖性糖尿病视网膜病变(NPDR)以及增殖性视网膜病变(PDR)。NPDR有三个亚类:轻度NPDR,中度NPDR和重度NPDR。非增殖型为病变早期,病变局限于视网膜内,表现为微血管瘤、出血、硬性和软性渗出物、视网膜动脉和静脉病变。增殖型病变至少有部分伸延超过内界膜,出现新生血管是增殖型的标志。DR的病灶详细描述如下:DR can be divided into two major categories, namely non-proliferative diabetic retinopathy (NPDR) and proliferative retinopathy (PDR). There are three subcategories of NPDR: mild NPDR, moderate NPDR, and severe NPDR. The non-proliferative type is the early stage of the lesion, and the lesion is limited to the retina, manifested as microhemangioma, hemorrhage, hard and soft exudate, retinal artery and venous lesions. Proliferative lesions at least partially extend beyond the inner limiting membrane, and the appearance of neovascularization is a sign of proliferative. DR lesions are described in detail as follows:
(1)微动脉瘤(1) Microaneurysm
微动脉瘤代表了视网膜损伤最原始的可感知迹象,视网膜血管的异常通透性导致微动脉 瘤的形成。在医学图像上,它是小的,圆形的,并且具有通常在黄斑周期的暗红色斑点。微动脉瘤可以看作是一个红点,边缘锐利,尺寸大小在20μm到200μm之间,近似于光盘大小的8.25%。Microaneurysms represent the most primitive perceptible signs of retinal damage. The abnormal permeability of retinal blood vessels leads to the formation of microaneurysms. On medical images, it is small, round, and has dark red spots that are usually in the macula cycle. The microaneurysm can be regarded as a red dot with sharp edges and a size between 20 μm and 200 μm, which is approximately 8.25% of the size of the disc.
(2)硬性渗出物(2) Hard exudates
与微动脉瘤不同,硬渗出物是脂蛋白和其他一些蛋白质从视网膜血管中漏出所形成的。在视觉上,它看起来像小的白色或黄白色沉积物,具有明显的边缘。硬性渗出物通常以环状形式的形式组织,通常在视网膜外层中出现。硬性渗出物通常是不规则的并且有光泽的,并且发现的位置接近微动脉瘤或视网膜水肿的边缘。Unlike microaneurysms, hard exudates are formed by the leakage of lipoproteins and other proteins from the blood vessels of the retina. Visually, it looks like small white or yellowish-white deposits with distinct edges. Rigid exudates are usually organized in the form of a ring, usually appearing in the outer layer of the retina. Rigid exudates are usually irregular and shiny, and are found close to the edges of microaneurysms or retinal edema.
(3)软性渗出物(3) Soft exudate
一般而言,软性渗出物是由于小动脉闭塞而形成的。流向视网膜的血流量减少导致视网膜神经纤维层(RNFL)缺血,最终影响轴浆流,从而在视网膜神经节细胞轴突上积累轴浆碎片。这种积聚可以像RNFL中的蓬松白色病变一样可视化,这通常被称为软性渗出物。In general, soft exudates are formed due to occlusion of small arteries. Reduced blood flow to the retina leads to ischemia of the retinal nerve fiber layer (RNFL), which ultimately affects axoplasmic flow, thereby accumulating axoplasmic debris on the axons of retinal ganglion cells. This accumulation can be visualized like a fluffy white lesion in RNFL, which is often referred to as soft exudate.
(4)出血(4) Bleeding
出血由于弱血管的渗漏而发生,出血的病灶在于具有不同密度和不均匀边缘的红点形式,且它是在125μm的范围内被发现的。从广义上讲,出血分为两类:火焰和斑点印迹出血,其中第一种类型起源于毛细血管前小动脉并出现在神经纤维上。第二种类型斑点印迹出血是圆形的并且小于微动脉瘤。斑点印迹出血可以出现在不同水平的视网膜上,但是,它大多数情况下会出现在毛细血管的经脉末端。Bleeding occurs due to leakage of weak blood vessels. The focus of bleeding is in the form of red dots with different densities and uneven edges, and it is found in the range of 125 μm. In a broad sense, hemorrhage is divided into two categories: flame and dot blotting. The first type originates from the anterior capillary artery and appears on nerve fibers. The second type of dot blot bleeding is round and smaller than a microaneurysm. Dot blot bleeding can appear on the retina at different levels, however, in most cases it will appear at the end of the capillaries' meridians.
(5)新血管形成(5) New blood vessel formation
新血管形成通常表示出现在视网膜内表面上的新血管的非典型出现。新血管很细小并且反复渗入玻璃体腔,这会降低视觉能力并使其显着模糊,最终导致失明。Neovascularization usually indicates the atypical appearance of new blood vessels that appear on the inner surface of the retina. The new blood vessels are small and repeatedly penetrate into the vitreous cavity, which reduces visual ability and makes it significantly blurry, which eventually leads to blindness.
(6)黄斑水肿(6) Macular edema
黄斑水肿被鉴定为视网膜的肿胀部分,其通常由于异常视网膜毛细血管的渗透性而发生。黄斑水肿导致黄斑周围的液体或其他溶质泄漏且严重影响视力。Macular edema is identified as a swollen part of the retina, which usually occurs due to the permeability of abnormal retinal capillaries. Macular edema causes leakage of fluid or other solutes around the macula and severely affects vision.
因此,现有技术中,针对包括糖尿病性视网膜病变(DR)在内的病变患者,仅依靠医生人眼通过观察病理图像判断是否发生病变以及病变程度,依靠肉眼观察、分析图像非常占用时间和精力,而且效率低下,不同医生的判断还存在差异性,这一问题亟待解决。Therefore, in the prior art, for patients with lesions including diabetic retinopathy (DR), only relying on the doctor ’s eye to determine whether the lesion occurred and the degree of the lesion by observing the pathological image, relying on the naked eye to observe and analyze the image takes a lot of time and effort And, the efficiency is low, and the judgments of different doctors are still different. This problem needs to be solved urgently.
发明内容Summary of the invention
本发明旨在至少在一定程度上解决相关技术中的技术问题之一。为此,本发明的一个目的是提供一种医学图像处理方法、系统、设备、存储介质,用于提高对医学图像的处理、分析效率。The present invention aims to solve one of the technical problems in the related art at least to a certain extent. To this end, an object of the present invention is to provide a medical image processing method, system, device, and storage medium for improving the efficiency of processing and analyzing medical images.
本发明所采用的技术方案是:The technical scheme adopted by the present invention is:
第一方面,本发明提供一种医学图像处理方法,包括以下步骤:In a first aspect, the present invention provides a medical image processing method, including the following steps:
病灶图像获取步骤,根据医学图像提取病灶图像;In the lesion image acquisition step, the lesion image is extracted based on the medical image;
得分生成步骤,根据所述病灶图像获取多个病灶得分,所述病灶得分为所述病灶图像属于不同图像异常程度的第一异常概率;In the score generation step, multiple lesion scores are obtained according to the lesion image, and the lesion score is a first abnormal probability that the lesion image belongs to different image abnormalities;
第二异常概率获取步骤,根据所述医学图像和所述病灶图像获取所述医学图像的多个第二异常概率,所述第二异常概率为所述医学图像属于不同图像异常程度的概率;A second abnormal probability acquisition step, acquiring multiple second abnormal probabilities of the medical image according to the medical image and the lesion image, the second abnormal probability being the probability that the medical image belongs to different image abnormalities;
分类步骤,根据多个所述第一异常概率、多个所述第二异常概率和不同的权重系数按照不同图像异常程度获取所述医学图像属于不同图像异常程度的最终概率,并将最大的所述最终概率的图像异常程度作为所述医学图像的图像异常程度。In the classification step, according to a plurality of the first abnormal probability, a plurality of the second abnormal probability and different weighting coefficients, the final probability that the medical image belongs to different image abnormalities is obtained according to different image abnormalities, and the largest The degree of image abnormality of the final probability is taken as the degree of image abnormality of the medical image.
进一步地,所述医学图像包括眼底照片。Further, the medical image includes a fundus photograph.
进一步地,所述病灶图像获取步骤包括:Further, the step of acquiring a lesion image includes:
根据所述医学图像和分割神经网络获取多种病灶掩模图像;Acquire a variety of focus mask images based on the medical image and the segmentation neural network;
根据所述病灶掩模图像和所述医学图像获取多种所述病灶图像。Acquire multiple lesion images based on the lesion mask image and the medical image.
进一步地,所述得分生成步骤包括:Further, the step of generating scores includes:
分别提取多种所述病灶图像的病灶特征,所述病灶特征包括颜色或形状;Extracting lesion features of multiple lesion images separately, the lesion features including color or shape;
拼接多个所述病灶特征并将其输入第一机器学习分类算法获取多个所述病灶得分。Splicing multiple lesion features and inputting them into the first machine learning classification algorithm to obtain multiple lesion scores.
进一步地,所述第二异常概率获取步骤包括:Further, the step of acquiring the second abnormal probability includes:
根据所述医学图像、多种所述病灶图像和第二机器学习分类算法获取所述医学图像的多个第二异常概率。Acquire a plurality of second abnormal probabilities of the medical image according to the medical image, a plurality of the lesion images, and a second machine learning classification algorithm.
进一步地,所述图像异常程度包括正常、轻度第一异常、中度第一异常、重度第一异常和第二异常。Further, the degree of image abnormality includes normal, mild first abnormality, moderate first abnormality, severe first abnormality, and second abnormality.
进一步地,所述第一异常概率的权重系数为0.2,所述第二异常概率的权重系数为0.8。Further, the weight coefficient of the first abnormal probability is 0.2, and the weight coefficient of the second abnormal probability is 0.8.
第二方面,本发明提供一种医学图像处理系统,包括:In a second aspect, the present invention provides a medical image processing system, including:
病灶图像获取单元,用于根据医学图像提取病灶图像;The lesion image acquisition unit is used to extract the lesion image according to the medical image;
得分生成单元,用于根据所述病灶图像获取多个病灶得分,所述病灶得分为所述病灶图 像属于不同图像异常程度的第一异常概率;A score generating unit, configured to obtain a plurality of lesion scores based on the lesion image, where the lesion score is a first abnormal probability that the lesion image belongs to different image abnormalities;
第二异常概率获取单元,用于根据所述医学图像和所述病灶图像获取所述医学图像的多个第二异常概率,所述第二异常概率为所述医学图像属于不同图像异常程度的概率;A second abnormal probability obtaining unit, configured to obtain a plurality of second abnormal probabilities of the medical image according to the medical image and the lesion image, the second abnormal probability being the probability that the medical image belongs to different image abnormalities ;
分类单元,用于根据多个所述第一异常概率、多个所述第二异常概率和不同的权重系数按照不同图像异常程度获取所述医学图像属于不同图像异常程度的最终概率,并将最大的所述最终概率的图像异常程度作为所述医学图像的图像异常程度。The classification unit is used to obtain the final probability that the medical image belongs to different image abnormalities according to different image abnormalities according to multiple first abnormal probabilities, multiple second abnormal probabilities and different weighting coefficients The image abnormality of the final probability is taken as the image abnormality of the medical image.
第三方面,本发明提供一种医学图像处理设备,包括:In a third aspect, the present invention provides a medical image processing device, including:
至少一个处理器;以及,At least one processor; and,
与所述至少一个处理器通信连接的存储器;其中,A memory communicatively connected to the at least one processor; wherein,
所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器能够执行所述的医学图像处理方法。The memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor to enable the at least one processor to execute the medical image processing method.
第四方面,本发明提供一种计算机可读存储介质,所述计算机可读存储介质存储有计算机可执行指令,所述计算机可执行指令用于使计算机执行所述的医学图像处理方法。In a fourth aspect, the present invention provides a computer-readable storage medium that stores computer-executable instructions that are used to cause a computer to execute the medical image processing method.
本发明的有益效果是:The beneficial effects of the present invention are:
本发明一方面通过病灶图像获取多个第一异常概率,另一方面通过医学图像和病灶图像获取多个第二异常概率,再根据多个第一异常概率、多个第二异常概率和不同的权重系数按照不同图像异常程度获取医学图像属于不同图像异常程度的最终概率,并将最大的最终概率的图像异常程度作为医学图像的图像异常程度,实现对医学图像的异常程度进行分析,克服现有技术中存在对病理图像处理、分析依靠肉眼,效率低下的技术问题。On the one hand, the present invention obtains multiple first abnormal probabilities through the focus image, on the other hand obtains multiple second abnormal probabilities through the medical image and focus image, and then according to multiple first abnormal probabilities, multiple second abnormal probabilities and different The weight coefficients are used to obtain the final probability of medical images belonging to different image abnormalities according to different image abnormalities, and the image abnormality of the maximum final probability is taken as the image abnormality of medical images to realize the analysis of medical image abnormality and overcome the existing There is a technical problem that the processing and analysis of pathological images rely on the naked eye and the efficiency is low.
另外,本发明还通过分割神经网络获取医学图像的多种病灶掩模图像,以进一步获取多种病灶图像,确保准确提取病灶图像。In addition, the present invention also obtains multiple lesion mask images of medical images by segmenting a neural network to further obtain multiple lesion images, ensuring accurate extraction of lesion images.
附图说明BRIEF DESCRIPTION
图1是本发明中一种医学图像处理方法的一具体实施例流程图;1 is a flowchart of a specific embodiment of a medical image processing method in the present invention;
图2是本发明中一种医学图像处理方法的分割神经网络的一具体实施例示意图。2 is a schematic diagram of a specific embodiment of a segmented neural network of a medical image processing method in the present invention.
具体实施方式detailed description
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。It should be noted that the embodiments in the present application and the features in the embodiments can be combined with each other if there is no conflict.
实施例1Example 1
一种医学图像处理方法,包括以下步骤:A medical image processing method includes the following steps:
病灶图像获取步骤,根据医学图像提取病灶图像;In the lesion image acquisition step, the lesion image is extracted based on the medical image;
得分生成步骤,根据病灶图像获取多个病灶得分,病灶得分为病灶图像属于不同图像异常程度的第一异常概率,例如,图像异常程度包括正常和异常两种,则得分生成步骤是获取病灶图像属于正常或异常的概率,获得两个第一异常概率,依此类推;In the score generation step, multiple lesion scores are obtained according to the lesion image. The lesion score is the first abnormal probability that the lesion image belongs to different image abnormalities. For example, if the image abnormality includes both normal and abnormal, the score generation step is to obtain the focus image. Normal or abnormal probability, get two first abnormal probabilities, and so on;
第二异常概率获取步骤,根据医学图像和病灶图像获取医学图像的多个第二异常概率,第二异常概率为医学图像属于不同图像异常程度的概率;In the second abnormal probability obtaining step, multiple second abnormal probabilities of the medical image are obtained according to the medical image and the lesion image, and the second abnormal probability is the probability that the medical image belongs to different image abnormalities;
分类步骤,根据多个第一异常概率、多个第二异常概率和不同的权重系数按照不同图像异常程度获取医学图像属于不同图像异常程度的最终概率,并将最大的最终概率的图像异常程度作为医学图像的图像异常程度。即将同一图像异常程度的第一异常概率和第二异常概率按照不同权重系数相加获得该图像异常程度的最终概率,由于有多个不同的图像异常程度,则可以获得医学图像的多个最终概率。The classification step is to obtain the final probabilities of medical images belonging to different image abnormalities according to different image abnormalities according to multiple first abnormal probabilities, multiple second abnormal probabilities and different weighting coefficients, and take the maximum final probability image abnormality as the The degree of abnormality of medical images. That is, the first abnormal probability and the second abnormal probability of the same image abnormal degree are added according to different weight coefficients to obtain the final probability of the image abnormal degree. Since there are multiple different image abnormal degrees, multiple final probabilities of the medical image can be obtained .
一种医学图像处理方法,一方面通过病灶图像获取多个第一异常概率,另一方面通过医学图像和病灶图像获取多个第二异常概率,再根据多个第一异常概率、多个第二异常概率和不同的权重系数按照不同图像异常程度获取医学图像属于不同图像异常程度的最终概率,并将最大的最终概率的图像异常程度作为医学图像的图像异常程度,实现对医学图像的异常程度进行分析,提高医学图像的处理和分析效率,克服现有技术中存在对病理图像处理、分析依靠肉眼,效率低下的技术问题。A medical image processing method. On the one hand, multiple first abnormal probabilities are obtained from focus images, on the other hand, multiple second abnormal probabilities are obtained from medical images and focus images, and then multiple first abnormal probabilities and multiple second Abnormal probability and different weight coefficients are used to obtain the final probability that medical images belong to different image abnormalities according to different image abnormalities, and the image abnormality of the maximum final probability is taken as the image abnormality of medical images to realize the abnormality of medical images. The analysis improves the processing and analysis efficiency of medical images, and overcomes the existing technical problems of low efficiency in relying on the naked eye for processing and analyzing pathological images.
进一步地,医学图像包括但不限于眼底照片,本实施例中,以眼底照片为例对医学图像处理方法进行说明,其余种类的医学图像的处理方法,可以参照以下对眼底照片的处理方法的具体描述,对图像的处理方法类似,不再赘述。而图像异常程度包括但不限于正常、轻度第一异常、中度第一异常、重度第一异常和第二异常,其中,以眼底照片为例,则第一异常为非增殖性病变,第二异常为增殖性病变,则眼底照片的图像异常程度包括但不限于正常、轻度非增殖性病变、中度非增殖性病变、重度非增殖性病变和增殖性病变。参考图1,图1是本发明中一种医学图像处理方法的一具体实施例流程图;下面对医学图像处理方法做具体说明:Further, medical images include, but are not limited to, fundus photos. In this embodiment, the fundus photos are used as examples to describe the medical image processing method. For other types of medical image processing methods, refer to the following specific methods for processing fundus photos The description is similar to the image processing method and will not be repeated here. The degree of image abnormality includes but is not limited to normal, mild first abnormality, moderate first abnormality, severe first abnormality, and second abnormality. Taking the fundus photo as an example, the first abnormality is a non-proliferative lesion. The second abnormality is a proliferative lesion. The abnormality of the fundus image includes, but is not limited to, normal, mild non-proliferative lesion, moderate non-proliferative lesion, severe non-proliferative lesion, and proliferative lesion. Referring to FIG. 1, FIG. 1 is a flowchart of a specific embodiment of a medical image processing method according to the present invention; the medical image processing method will be specifically described below:
第一步:对数据进行预处理。本发明使用的数据为二维彩色眼底照片,原图片大小为3000×3000,需对图片进行预处理,先将图片的背景等多余部分剪裁掉,以将图片的大小调整为128×128。Step 1: Pre-process the data. The data used in the present invention is a two-dimensional color fundus photograph. The original picture size is 3000 × 3000. The picture needs to be pre-processed, and the background and other unnecessary parts of the picture are cropped to adjust the picture size to 128 × 128.
第二步:病灶图像获取步骤,包括:Step 2: The step of acquiring lesion images, including:
根据医学图像和分割神经网络获取多种病灶掩模图像;Acquire multiple lesion mask images based on medical images and segmented neural networks;
根据病灶掩模图像和医学图像获取多种病灶图像。Acquire multiple lesion images based on the lesion mask image and medical image.
其中,分割神经网络的主要任务是实现多种病灶分割和提取,本实施例中,病灶种类包括微动脉瘤、硬性渗出物以及新生视网膜血管,微动脉瘤、硬性渗出物属于非增殖性病变,新生视网膜血管属于增殖性病变。将待测图片输入到训练好的分割神经网络中,通过分割神经网络预测图片上的每个像素是属于背景、微动脉瘤、硬性渗出物还是视网膜血管,并且将预测到的结果输出成三个掩模,分别是微动脉瘤掩模、硬性渗出物掩模和视网膜血管掩模。掩模具体是由0和1组成,背景为0,目标为1,例如生成微动脉瘤掩模时,通过分割神经网络预测每一个像素是否属于微动脉瘤区域,若是该像素位置标记为1,否则该像素位置标记为0,最终得到与原图大小一样的大小为128×128的微动脉瘤掩模图。掩模可以准确描述目标的所在位置,将掩模对应到原图,即可提取到各病灶分割后的图像,也即病灶图像。Among them, the main task of segmenting the neural network is to realize the segmentation and extraction of multiple lesions. In this embodiment, the types of lesions include microaneurysms, rigid exudates, and neonatal retinal blood vessels. Lesions, new retinal blood vessels are proliferative lesions. Input the image to be tested into the trained segmentation neural network, and predict whether each pixel on the image belongs to the background, microaneurysm, hard exudate or retinal blood vessel through the segmentation neural network, and output the predicted result into three The three masks are the micro-aneurysm mask, the hard exudate mask, and the retinal vascular mask. The mask is specifically composed of 0 and 1, the background is 0, and the target is 1. For example, when generating a micro-aneurysm mask, predict whether each pixel belongs to the micro-aneurysm area by segmenting the neural network. If the pixel position is marked as 1, Otherwise, the pixel position is marked as 0, and finally a micro-aneurysm mask image with the same size as the original image of 128 × 128 is obtained. The mask can accurately describe the location of the target. Corresponding the mask to the original image, you can extract the image of each lesion, which is the lesion image.
参考表1和图2,图2是本发明中一种医学图像处理方法的分割神经网络的一具体实施例示意图;分割神经网络的实现主要基于多任务学习框架,分割神经网络的分割任务有三个,分别是分割微动脉瘤、分割硬性渗出物以及分割视网膜血管。网络的第一部分称为“共享层”,这些层的权重对于这三个任务都是通用的。接着,是每个任务的特定层,称之为“任务特定层”,这些层之间的参数都是独立计算的,不参与跨层共享。在这些特定任务层中,网络学习特定于任务的信息。每个这些独立的任务层都会生成一个单独的输出结果,具体是输出微动脉瘤、硬性渗出物、视网膜血管的分割掩模。Referring to Table 1 and FIG. 2, FIG. 2 is a schematic diagram of a specific embodiment of a segmentation neural network of a medical image processing method in the present invention; the implementation of the segmentation neural network is mainly based on a multi-task learning framework, and there are three segmentation tasks of the segmentation neural network , Respectively, segmentation of microaneurysm, segmentation of hard exudate, and segmentation of retinal vessels. The first part of the network is called the "shared layer", and the weights of these layers are common to all three tasks. Next, there is a specific layer for each task, called "task-specific layer". The parameters between these layers are calculated independently and do not participate in cross-layer sharing. In these task-specific layers, the network learns task-specific information. Each of these independent task layers will generate a separate output, specifically a segmentation mask that outputs microaneurysms, hard exudates, and retinal blood vessels.
表1分割神经网络层次结构Table 1 Segmented neural network hierarchy
Figure PCTCN2018124660-appb-000001
Figure PCTCN2018124660-appb-000001
Figure PCTCN2018124660-appb-000002
Figure PCTCN2018124660-appb-000002
上述分割神经网络能同时分割出多种病灶图像,不限于本实施例中的三种病灶图像。而在使用分割神经网络之前,需要对网络进行训练,具体地,首先初始化模型,然后输入已标注好掩模标记数据的样本进行训练,训练时采用整张图片进行训练,不进行补丁抽样,而实验证明采用全图更加有效。输入图片后,分别进行三次卷积和下采样得到特征小图(5×5×128),然后再进行一次卷积得到特征(1×1×1024),接着进行三组全连接达到三个任务的特征,三个任务分别进行反卷积和上采样,得到与原图一样大小的大图,这些图就是分割出来的掩模。分割神经网络模型在训练过程中不断优化,调整参数直至收敛,优化的目标是不断缩小网络预测出来的分割图(病灶图像)与样本标注数据之间的差异,即最小化损失函数L1:The above segmentation neural network can segment multiple lesion images simultaneously, and is not limited to the three lesion images in this embodiment. Before using the segmented neural network, the network needs to be trained. Specifically, the model is first initialized, and then the sample with the marked mask data is input for training. The entire picture is used for training without patch sampling. The experiment proves that it is more effective to use the whole picture. After inputting the picture, perform three times convolution and downsampling to get the feature small image (5 × 5 × 128), and then perform another convolution to get the feature (1 × 1 × 1024), and then perform three sets of full connection to achieve three tasks The features of the three tasks are deconvolution and upsampling, respectively, to get a large picture of the same size as the original picture, these pictures are the masks that are segmented out. The segmentation neural network model is continuously optimized during the training process, and the parameters are adjusted until convergence. The goal of optimization is to continuously reduce the difference between the segmentation map (lesion image) predicted by the network and the sample annotation data, that is, to minimize the loss function L1:
Figure PCTCN2018124660-appb-000003
Figure PCTCN2018124660-appb-000003
上式中,N为训练样本数,L MA表示微动脉瘤分割的像素级损失,L HES表示硬性渗出物的分割损失,L NV表示视网膜血管的分割损失;f(x,θ)为分割神经网络预测的分割结果,其中x为样本的一个像素,θ为学习率;S 1,S 2,S 3分别是微动脉瘤、硬性渗出物、视网膜血管的样本标注位置;Φ(θ)是正则项。 In the above formula, N is the number of training samples, L MA is the pixel-level loss of micro-aneurysm segmentation, L HES is the segmentation loss of hard exudate, and L NV is the segmentation loss of retinal vessels; f (x, θ) is the segmentation The segmentation result predicted by the neural network, where x is a pixel of the sample and θ is the learning rate; S 1 , S 2 , and S 3 are the labeled positions of the samples of microaneurysm, hard exudate, and retinal blood vessel; Φ (θ) Is a regular item.
则利用训练好的分割神经网络可以获得眼底照片的多种病灶掩模图像,包括微动脉瘤掩模图像、硬性渗出物掩模图像、以及视网膜血管掩模图像;根据病灶掩模图像和医学图像即可以分割得到多种病灶图像,分别是微动脉瘤病灶图像、硬性渗出物病灶图像、以及视网膜血管病灶图像。Then use the trained segmentation neural network to obtain a variety of lesion mask images of the fundus photos, including micro-aneurysm mask images, hard exudate mask images, and retinal vascular mask images; according to the lesion mask images and medicine The image can be segmented to obtain a variety of lesion images, including micro-aneurysm lesion images, rigid exudate lesion images, and retinal vascular lesion images.
第三步:得分生成步骤,包括:Step 3: The score generation step, including:
分别提取多种病灶图像的病灶特征,病灶特征包括颜色或形状;Extract the lesion features of multiple lesion images separately, and the lesion features include color or shape;
拼接多个病灶特征并将其输入第一机器学习分类算法以获取多个病灶得分。The features of multiple lesions are spliced and input into the first machine learning classification algorithm to obtain multiple lesion scores.
具体地,在获得各个病灶图像之后,根据病灶图像提取各病灶的颜色、形状等如表2所 示的特征,每个病灶图像可以提取出来86维特征,并将3个病灶图像的病灶特征按微动脉瘤病灶图像、硬性渗出物病灶图像、视网膜血管病灶图像的顺序拼接起来成为258维特征向量,并将特征向量输入到得分生成器(即第一机器学习分类算法)中即可获得多个病灶得分为
Figure PCTCN2018124660-appb-000004
i为图像异常程度的类别,X为医学图像的编号,即获得病灶图像属于不同图像异常程度(包括正常、轻度非增殖性病变、中度非增殖性病变、重度非增殖性病变和增殖性病变)的概率,即可以获得5个病灶得分。
Specifically, after obtaining the images of each lesion, the features such as the color and shape of each lesion are extracted according to the lesion image as shown in Table 2, 86-dimensional features can be extracted from each lesion image, and the lesion features of the 3 lesion images are The sequence of micro-aneurysm lesion images, rigid exudate lesion images, and retinal vascular lesion images are stitched together into a 258-dimensional feature vector, and the feature vector is input into the score generator (that is, the first machine learning classification algorithm) to obtain multiple The score of each lesion is
Figure PCTCN2018124660-appb-000004
i is the category of the abnormality of the image, X is the number of the medical image, that is, the acquired lesion image belongs to different image abnormality (including normal, mild non-proliferative lesions, moderate non-proliferative lesions, severe non-proliferative lesions and proliferative Disease), that is, the score of 5 lesions can be obtained.
表2病灶特征提取列表Table 2 List of lesion feature extraction
特征feature 特征描述Feature description
f 1 f 1 RGB图像上的均值Mean on RGB image
f 2 f 2 RGB图像上的方差Variance on RGB images
f 3-f 12 f 3 -f 12 RGB图像上的颜色矩Color moments on RGB images
f 13 f 13 目标区域的周长和The perimeter of the target area and
f 14 f 14 目标区域的面积和The area of the target area and
f 15-f 74 f 15 -f 74 目标区域的LBP特征LBP features of the target area
f 75-f 86 f 75 -f 86 HSV特征HSV characteristics
本实施例中,得分生成器采用随机森林算法。得分生成器在使用之前需要进行训练,具体地,将训练样本的病灶特征258维特征向量以及该样本标注的病变程度输入到得分生成器中进行训练。训练好的得分生成器通过计算可以得出输入的病灶属于5种图像异常程度(即正常、轻度非增殖性病变、中度非增殖性病变、重度非增殖性病变以及增殖性病变)的概率,取此概率值作为该样本的病灶得分。In this embodiment, the score generator uses a random forest algorithm. The score generator needs to be trained before it is used. Specifically, the 258-dimensional feature vector of the lesion feature of the training sample and the degree of lesion marked by the sample are input into the score generator for training. The trained score generator can calculate the probability that the input lesion belongs to 5 image abnormalities (that is, normal, mild non-proliferative lesions, moderate non-proliferative lesions, severe non-proliferative lesions, and proliferative lesions) , Take this probability value as the lesion score of the sample.
第四步:第二异常概率获取步骤,包括:Step 4: The second abnormal probability acquisition step, including:
根据医学图像、多种病灶图像和第二机器学习分类算法获取医学图像的多个第二异常概率。Acquire multiple second abnormal probabilities of the medical image according to the medical image, multiple lesion images, and the second machine learning classification algorithm.
本实施例中,图像异常程度包括正常、轻度非增殖性病变、中度非增殖性病变、重度非增殖性病变和增殖性病变;第二机器学习分类算法是基于3D卷积神经网络进行设计的,3D卷积神经网络可同时输入原来的医学图像(即眼底照片)以及微动脉瘤病灶图像、硬性渗出物病灶图像和视网膜血管病灶图像,能从多维度提取特征,然后进行3D卷积,以捕捉从多个图像中得到的特征信息。具体地,参考表3,3D卷积神经网络包含一个实连线层(即硬连接层)、4个卷积层、3个下采样层和一个全连接层。每个卷积核卷积的立方体是包括原图、微血管瘤病灶图像、视网膜血管病灶图像、硬性渗出物病灶图像的4张图片,每张图片的大 小是128×128。在第一层,应用了一个固定的实连线的核去对输入的图进行处理,产生多个通道的信息,然后对多个通道分别处理,最后再将所有通道的信息组合起来得到最终的特征描述。这个实线层实际上是编码了对病灶的先验知识,这比随机初始化性能要好。每张图提取五个通道的信息,分别是:R、G、B颜色通道,以及x和y方向的梯度。所以有4×5=20个病灶图。然后再用一个9×9×2的3D卷积核在五个通道的每一个通道分别进行卷积。为了增加特征图的个数即提取不同的特征,在每一个位置都采用两个不同的卷积核,这样L 1层中包含了两个特征图组,每组都包含(4-2+1)×5=15个特征图。紧接着的是下采样层L 3层,在L 2层的特征图中用3×3窗口进行下采样,这样就会得到相同数目但是空间分辨率降低的病灶图。L 4是在5个通道中分别采用9×9×2的3D卷积核计算得到。为了增加特征图个数,在每个位置都采用三个不同的卷积核,这样就可以得到6组不同的特征图,每组有(4-2+1-2+1)×5=10个特征图。L 5层用的是2×2的下采样窗口,接着再进行两次2D卷积和下采样得到128×1×1的特征向量,而128维是根据以往的经验确定的。输入医学图像后经过多层的卷积和下采样,输入的多个图像都被转化为一个128维的特征向量,这个特征向量捕捉了视网膜眼底照片的特征信息。再将得到的128维的特征向量输入到softmax层,softmax层的节点数与图像异常程度的类别数目一致,而且每个节点与L 9中的128个节点是全连接的。最终softmax层会得到一系列预测值,即医学图像属于5种图像异常程度的概率值,即第二异常概率,假设此值为
Figure PCTCN2018124660-appb-000005
其中,i(i=0,1,2,3,4)为图像异常程度类别(即正常、轻度非增殖性病变、中度非增殖性病变、重度非增殖性病变和增殖性病变),X为输入的医学图像编号。
In this embodiment, the image abnormality includes normal, mild non-proliferative lesions, moderate non-proliferative lesions, severe non-proliferative lesions, and proliferative lesions; the second machine learning classification algorithm is designed based on 3D convolutional neural network The 3D convolutional neural network can simultaneously input the original medical images (ie, fundus photos), microaneurysm lesion images, rigid exudate lesion images, and retinal vascular lesion images, which can extract features from multiple dimensions and then perform 3D convolution To capture feature information from multiple images. Specifically, referring to Table 3, the 3D convolutional neural network includes a solid-wired layer (ie, hard-connected layer), 4 convolutional layers, 3 down-sampled layers, and a fully connected layer. The convolved cubes of each convolution kernel are 4 pictures including the original image, the image of the microhemangioma lesion, the image of the retinal vascular lesion, and the image of the hard exudate lesion. The size of each picture is 128 × 128. At the first layer, a fixed solid-line kernel is applied to process the input graph to generate multiple channels of information, then multiple channels are processed separately, and finally the information of all channels are combined to obtain the final Feature description. This solid layer actually encodes a priori knowledge of the lesion, which is better than random initialization performance. Each picture extracts the information of five channels, namely: R, G, B color channels, and gradients in the x and y directions. So there are 4 × 5 = 20 lesion maps. Then, a 9 × 9 × 2 3D convolution kernel is used to perform convolution on each of the five channels. In order to increase the number of feature maps to extract different features, two different convolution kernels are used at each location, so that the L 1 layer contains two feature map groups, each group contains (4-2 + 1 ) × 5 = 15 feature maps. Immediately following is the down-sampling layer L 3 layer, which is down-sampled with a 3 × 3 window in the feature map of the L 2 layer, so that the same number of lesion maps with reduced spatial resolution will be obtained. L 4 is calculated by using 9 × 9 × 2 3D convolution kernels in 5 channels. In order to increase the number of feature maps, three different convolution kernels are used at each position, so that 6 different feature maps can be obtained, each group has (4-2 + 1-2 + 1) × 5 = 10 Feature maps. The L 5 layer uses a 2 × 2 downsampling window, followed by two 2D convolutions and downsampling to obtain a 128 × 1 × 1 feature vector, and the 128 dimensions are determined based on past experience. After inputting medical images, after multiple layers of convolution and downsampling, the input multiple images are converted into a 128-dimensional feature vector, which captures the feature information of the retinal fundus photo. Then the 128-dimensional feature vector obtained is input to the softmax layer. The number of nodes in the softmax layer is consistent with the number of image abnormality categories, and each node is fully connected to the 128 nodes in L 9 . Finally, the softmax layer will obtain a series of prediction values, that is, the probability value of the medical image belonging to the five image abnormalities, that is, the second abnormal probability, assuming this value
Figure PCTCN2018124660-appb-000005
Among them, i (i = 0,1,2,3,4) is the category of image abnormality (ie normal, mild non-proliferative lesions, moderate non-proliferative lesions, severe non-proliferative lesions and proliferative lesions), X is the input medical image number.
表3 3D卷积神经网络层次结构Table 3 3D Convolutional Neural Network Hierarchy
层数Number of layers 类型Types of 神经元/特征尺寸Neuron / feature size 卷积核/窗口尺寸Convolution kernel / window size
00 输入层Input layer 4×128×1284 × 128 × 128 ---
11 硬连接层Hard link layer 20×128×12820 × 128 × 128 ---
22 卷积层Convolutional layer 15×2×120×12015 × 2 × 120 × 120 9×9×29 × 9 × 2
33 下采样层Downsampling layer 15×2×40×4015 × 2 × 40 × 40 3×33 × 3
44 卷积层Convolutional layer 10×6×32×3210 × 6 × 32 × 32 9×9×29 × 9 × 2
55 下采样层Downsampling layer 10×6×16×1610 × 6 × 16 × 16 2×22 × 2
66 卷积层Convolutional layer 10×6×10×1010 × 6 × 10 × 10 7×77 × 7
77 下采样层Downsampling layer 10×6×5×510 × 6 × 5 × 5 2×22 × 2
88 卷积层Convolutional layer 10×6×1×110 × 6 × 1 × 1 5×55 × 5
99 全连接层Fully connected layer 128×1×1128 × 1 × 1 ---
1010 softmax层softmax layer  A  A
在使用3D卷积神经网络之前,需要对其进行训练,首先初始化模型,然后输入已标注异常程度的样本进行训练,训练时输入样本的原图、微血管瘤病灶图像、硬性渗出物病灶图像、视网膜血管病灶图像的4张图片进行训练。3D卷积神经网络模型在训练过程中不断优化,调整参数直至收敛,优化的目标是不断缩小网络预测出来的分类结果与样本标注异常程度之间的差异,即最小化损失函数L2:Before using the 3D convolutional neural network, you need to train it, first initialize the model, and then input the sample marked with abnormal degree for training. Input the original image of the sample, the image of the microhemangioma lesion, and the image of the hard exudate lesion during training. 4 pictures of retinal vascular lesion images for training. The 3D convolutional neural network model is continuously optimized during the training process, and the parameters are adjusted until convergence. The goal of optimization is to continuously reduce the difference between the classification results predicted by the network and the abnormal degree of sample labeling, that is, to minimize the loss function L2:
Figure PCTCN2018124660-appb-000006
Figure PCTCN2018124660-appb-000006
上式中,N为训练样本数,L softmax表示网络分类损失,标注数据;f(x,θ)为网络预测的分类结果,其中x为一个样本,θ为学习率;C为标注类别;Φ(θ)是正则项。 In the above formula, N is the number of training samples, L softmax represents the network classification loss, and annotated data; f (x, θ) is the classification result predicted by the network, where x is a sample, θ is the learning rate; C is the label category; Φ (θ) is the regular term.
第五步:得分融合分类,具体地,在获得第二异常概率
Figure PCTCN2018124660-appb-000007
(本实施例中有5个第二异常概率)和第一异常概率
Figure PCTCN2018124660-appb-000008
(本实施例中有5个第一异常概率,即病灶得分),根据设置不同的权重系数以获取最终的异常概率,本实施例中,第一异常概率的权重系数为0.2,第二异常概率的权重系数为0.8,则最终的异常概率,即得分融合值为
Figure PCTCN2018124660-appb-000009
Figure PCTCN2018124660-appb-000010
(本实施例中有5个得分融合值),在算出的结果中,哪个图像异常程度类别的得分融合值最高,则该类别即为医学图像的所属图像异常程度类别,即求C=argmax iF i(X)。
Step 5: Score fusion classification, specifically, in obtaining the second anomaly probability
Figure PCTCN2018124660-appb-000007
(There are 5 second abnormal probabilities in this embodiment) and the first abnormal probability
Figure PCTCN2018124660-appb-000008
(There are 5 first abnormal probabilities in this embodiment, that is, the lesion score), according to the setting of different weight coefficients to obtain the final abnormal probability, in this embodiment, the weight coefficient of the first abnormal probability is 0.2, the second abnormal probability The weight coefficient of is 0.8, then the final abnormal probability, that is, the score fusion value is
Figure PCTCN2018124660-appb-000009
Figure PCTCN2018124660-appb-000010
(There are 5 score fusion values in this embodiment). In the calculation result, which image abnormality category has the highest score fusion value, the category is the image abnormality category of the medical image, that is, C = argmax i F i (X).
一种医学图像处理方法,提出得分融合机制,将由各病灶图像获得的病灶得分融入分类网络,提高了对各病灶的敏感度和分类准确性。且本发明的得分融合分类框架不仅局限于本发明中检测出的三种病灶,待检测出更多病灶时可进行扩充,系统适应性强。A medical image processing method, a score fusion mechanism is proposed, and the lesion scores obtained from each lesion image are integrated into a classification network, which improves the sensitivity and classification accuracy of each lesion. Moreover, the score fusion classification framework of the present invention is not limited to the three lesions detected in the present invention, and can be expanded when more lesions are detected, and the system has strong adaptability.
实施例2Example 2
一种医学图像处理系统,包括:A medical image processing system, including:
病灶图像获取单元,用于根据医学图像提取病灶图像;The lesion image acquisition unit is used to extract the lesion image according to the medical image;
得分生成单元,用于根据病灶图像获取多个病灶得分,病灶得分为病灶图像属于不同图像异常程度的第一异常概率;The score generating unit is used to obtain multiple lesion scores according to the lesion image, and the lesion score is the first abnormal probability that the lesion image belongs to different image abnormalities;
第二异常概率获取单元,用于根据医学图像和病灶图像获取医学图像的多个第二异常概率,第二异常概率为医学图像属于不同图像异常程度的概率;A second abnormal probability acquisition unit, configured to acquire multiple second abnormal probabilities of the medical image according to the medical image and the lesion image, the second abnormal probability is the probability that the medical image belongs to different image abnormalities;
分类单元,用于根据多个第一异常概率、多个第二异常概率和不同的权重系数按照不同 图像异常程度获取医学图像属于不同图像异常程度的最终概率,并将最大的最终概率的图像异常程度作为医学图像的图像异常程度。The classification unit is used to obtain the final probabilities that medical images belong to different image abnormalities according to different image abnormalities according to multiple first abnormal probabilities, multiple second abnormal probabilities and different weight coefficients, and to abnormalize the images with the largest final probability The degree is regarded as the abnormal degree of the medical image.
医学图像处理系统的具体工作过程和有益效果请参照实施例1中医学图像处理方法的具体描述,不再赘述。For the specific working process and beneficial effects of the medical image processing system, please refer to the specific description of the medical image processing method in Embodiment 1, which will not be repeated here.
实施例3Example 3
一种医学图像处理设备,包括:A medical image processing device, including:
至少一个处理器;以及,At least one processor; and,
与所述至少一个处理器通信连接的存储器;其中,A memory communicatively connected to the at least one processor; wherein,
所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器能够执行所述的医学图像处理方法。医学图像处理方法的具体描述参照实施例1中的描述,不再赘述。The memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor to enable the at least one processor to execute the medical image processing method. For a detailed description of the medical image processing method, refer to the description in Embodiment 1, and no further description is required.
实施例4Example 4
一种计算机可读存储介质,所述计算机可读存储介质存储有计算机可执行指令,所述计算机可执行指令用于使计算机执行所述的医学图像处理方法。医学图像处理方法的具体描述参照实施例1中的描述,不再赘述。A computer-readable storage medium stores computer-executable instructions for causing a computer to execute the medical image processing method. For a detailed description of the medical image processing method, refer to the description in Embodiment 1, and no further description is required.
以上是对本发明的较佳实施进行了具体说明,但本发明创造并不限于所述实施例,熟悉本领域的技术人员在不违背本发明精神的前提下还可做出种种的等同变形或替换,这些等同的变形或替换均包含在本申请权利要求所限定的范围内。The above is a specific description of the preferred implementation of the present invention, but the invention is not limited to the embodiments, and those skilled in the art can make various equivalent modifications or replacements without violating the spirit of the present invention These equivalent variations or replacements are included in the scope defined by the claims of this application.

Claims (10)

  1. 一种医学图像处理方法,其特征在于,包括以下步骤:A medical image processing method, characterized in that it includes the following steps:
    病灶图像获取步骤,根据医学图像提取病灶图像;In the lesion image acquisition step, the lesion image is extracted based on the medical image;
    得分生成步骤,根据所述病灶图像获取多个病灶得分,所述病灶得分为所述病灶图像属于不同图像异常程度的第一异常概率;In the score generation step, multiple lesion scores are obtained according to the lesion image, and the lesion score is a first abnormal probability that the lesion image belongs to different image abnormalities;
    第二异常概率获取步骤,根据所述医学图像和所述病灶图像获取所述医学图像的多个第二异常概率,所述第二异常概率为所述医学图像属于不同图像异常程度的概率;A second abnormal probability acquisition step, acquiring multiple second abnormal probabilities of the medical image according to the medical image and the lesion image, the second abnormal probability being the probability that the medical image belongs to different image abnormalities;
    分类步骤,根据多个所述第一异常概率、多个所述第二异常概率和不同的权重系数按照不同图像异常程度获取所述医学图像属于不同图像异常程度的最终概率,并将最大的所述最终概率的图像异常程度作为所述医学图像的图像异常程度。In the classification step, according to a plurality of the first abnormal probability, a plurality of the second abnormal probability and different weighting coefficients, the final probability that the medical image belongs to different image abnormalities is obtained according to different image abnormalities, and the largest The degree of image abnormality of the final probability is taken as the degree of image abnormality of the medical image.
  2. 根据权利要求1所述的医学图像处理方法,其特征在于,所述医学图像包括眼底照片。The medical image processing method according to claim 1, wherein the medical image includes a fundus photograph.
  3. 根据权利要求1所述的医学图像处理方法,其特征在于,所述病灶图像获取步骤包括:The medical image processing method according to claim 1, wherein the lesion image acquisition step includes:
    根据所述医学图像和分割神经网络获取多种病灶掩模图像;Acquire a variety of focus mask images based on the medical image and the segmentation neural network;
    根据所述病灶掩模图像和所述医学图像获取多种所述病灶图像。Acquire multiple lesion images based on the lesion mask image and the medical image.
  4. 根据权利要求3所述的医学图像处理方法,其特征在于,所述得分生成步骤包括:The medical image processing method according to claim 3, wherein the score generation step includes:
    分别提取多种所述病灶图像的病灶特征,所述病灶特征包括颜色或形状;Extracting lesion features of multiple lesion images separately, the lesion features including color or shape;
    拼接多个所述病灶特征并将其输入第一机器学习分类算法获取多个所述病灶得分。Splicing multiple lesion features and inputting them into the first machine learning classification algorithm to obtain multiple lesion scores.
  5. 根据权利要求3所述的医学图像处理方法,其特征在于,所述第二异常概率获取步骤包括:The medical image processing method according to claim 3, wherein the step of acquiring the second abnormal probability includes:
    根据所述医学图像、多种所述病灶图像和第二机器学习分类算法获取所述医学图像的多个第二异常概率。Acquire a plurality of second abnormal probabilities of the medical image according to the medical image, a plurality of the lesion images, and a second machine learning classification algorithm.
  6. 根据权利要求1至5任一项所述的医学图像处理方法,其特征在于,所述图像异常程度包括正常、轻度第一异常、中度第一异常、重度第一异常和第二异常。The medical image processing method according to any one of claims 1 to 5, wherein the degree of image abnormality includes normal, mild first abnormality, moderate first abnormality, severe first abnormality, and second abnormality.
  7. 根据权利要求1至5任一项所述的医学图像处理方法,其特征在于,所述第一异常概率的权重系数为0.2,所述第二异常概率的权重系数为0.8。The medical image processing method according to any one of claims 1 to 5, wherein the weight coefficient of the first abnormal probability is 0.2, and the weight coefficient of the second abnormal probability is 0.8.
  8. 一种医学图像处理系统,其特征在于,包括:A medical image processing system, characterized in that it includes:
    病灶图像获取单元,用于根据医学图像提取病灶图像;The lesion image acquisition unit is used to extract the lesion image according to the medical image;
    得分生成单元,用于根据所述病灶图像获取多个病灶得分,所述病灶得分为所述病灶图像属于不同图像异常程度的第一异常概率;A score generating unit, configured to obtain multiple lesion scores based on the lesion image, where the lesion score is a first abnormal probability that the lesion image belongs to different image abnormalities;
    第二异常概率获取单元,用于根据所述医学图像和所述病灶图像获取所述医学图像的多个第二异常概率,所述第二异常概率为所述医学图像属于不同图像异常程度的概率;A second abnormal probability obtaining unit, configured to obtain a plurality of second abnormal probabilities of the medical image according to the medical image and the lesion image, the second abnormal probability being the probability that the medical image belongs to different image abnormalities ;
    分类单元,用于根据多个所述第一异常概率、多个所述第二异常概率和不同的权重系数按照不同图像异常程度获取所述医学图像属于不同图像异常程度的最终概率,并将最大的所述最终概率的图像异常程度作为所述医学图像的图像异常程度。The classification unit is used to obtain the final probability that the medical image belongs to different image abnormalities according to different image abnormalities according to multiple first abnormal probabilities, multiple second abnormal probabilities and different weighting coefficients The image abnormality of the final probability is taken as the image abnormality of the medical image.
  9. 一种医学图像处理设备,其特征在于,包括:A medical image processing device, characterized in that it includes:
    至少一个处理器;以及,At least one processor; and,
    与所述至少一个处理器通信连接的存储器;其中,A memory communicatively connected to the at least one processor; wherein,
    所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,以使所述至少一个处理器能够执行如权利要求1至7任一项所述的医学图像处理方法。The memory stores instructions executable by the at least one processor, and the instructions are executed by the at least one processor to enable the at least one processor to execute any one of claims 1 to 7. Medical image processing method.
  10. 一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有计算机可执行指令,所述计算机可执行指令用于使计算机执行如权利要求1至7任一项所述的医学图像处理方法。A computer-readable storage medium, characterized in that the computer-readable storage medium stores computer-executable instructions for causing a computer to perform the medicine according to any one of claims 1 to 7. Image processing method.
PCT/CN2018/124660 2018-11-07 2018-12-28 Medical image processing method, system, device, and storage medium WO2020093563A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201811320629.9A CN109671049B (en) 2018-11-07 2018-11-07 Medical image processing method, system, equipment and storage medium
CN201811320629.9 2018-11-07

Publications (1)

Publication Number Publication Date
WO2020093563A1 true WO2020093563A1 (en) 2020-05-14

Family

ID=66142041

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2018/124660 WO2020093563A1 (en) 2018-11-07 2018-12-28 Medical image processing method, system, device, and storage medium

Country Status (2)

Country Link
CN (1) CN109671049B (en)
WO (1) WO2020093563A1 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115311188A (en) * 2021-05-08 2022-11-08 数坤(北京)网络科技股份有限公司 Image identification method and device, electronic equipment and storage medium

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111144296B (en) * 2019-12-26 2023-04-18 湖南大学 Retina fundus picture classification method based on improved CNN model
CN111932562B (en) * 2020-09-22 2021-01-19 平安科技(深圳)有限公司 Image identification method and device based on CT sequence, electronic equipment and medium
CN112652394A (en) * 2021-01-14 2021-04-13 浙江工商大学 Multi-focus target detection-based retinopathy of prematurity diagnosis system
CN115816834B (en) * 2023-02-20 2023-04-25 常熟理工学院 Method and system for real-time monitoring of printing quality of printer

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013031536A1 (en) * 2011-08-30 2013-03-07 大日本印刷株式会社 Image analysis device for ophthalmic disease, image analysis method for ophthalmic disease, and image analysis program for ophthalmic disease
CN105310645A (en) * 2014-06-18 2016-02-10 佳能株式会社 Image processing apparatus and image processing method
CN106530295A (en) * 2016-11-07 2017-03-22 首都医科大学 Fundus image classification method and device of retinopathy
CN106530290A (en) * 2016-10-27 2017-03-22 朱育盼 Medical image analysis method and device
CN107045720A (en) * 2017-05-04 2017-08-15 深圳硅基智能科技有限公司 Artificial neural network and system for recognizing eye fundus image lesion
CN107145756A (en) * 2017-05-17 2017-09-08 上海辉明软件有限公司 A kind of stroke types Forecasting Methodology and device
CN107679525A (en) * 2017-11-01 2018-02-09 腾讯科技(深圳)有限公司 Image classification method, device and computer-readable recording medium
CN108734209A (en) * 2018-05-16 2018-11-02 上海鹰瞳医疗科技有限公司 Feature recognition based on more images and equipment

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9972092B2 (en) * 2016-03-31 2018-05-15 Adobe Systems Incorporated Utilizing deep learning for boundary-aware image segmentation
DE202017104953U1 (en) * 2016-08-18 2017-12-04 Google Inc. Processing fundus images using machine learning models
CN106548192B (en) * 2016-09-23 2019-08-09 北京市商汤科技开发有限公司 Image processing method, device and electronic equipment neural network based
CN108230294B (en) * 2017-06-14 2020-09-29 北京市商汤科技开发有限公司 Image detection method, image detection device, electronic equipment and storage medium
CN108010021B (en) * 2017-11-30 2021-12-10 上海联影医疗科技股份有限公司 Medical image processing system and method
CN108634934B (en) * 2018-05-07 2021-01-29 北京长木谷医疗科技有限公司 Method and apparatus for processing spinal sagittal image

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013031536A1 (en) * 2011-08-30 2013-03-07 大日本印刷株式会社 Image analysis device for ophthalmic disease, image analysis method for ophthalmic disease, and image analysis program for ophthalmic disease
CN105310645A (en) * 2014-06-18 2016-02-10 佳能株式会社 Image processing apparatus and image processing method
CN106530290A (en) * 2016-10-27 2017-03-22 朱育盼 Medical image analysis method and device
CN106530295A (en) * 2016-11-07 2017-03-22 首都医科大学 Fundus image classification method and device of retinopathy
CN107045720A (en) * 2017-05-04 2017-08-15 深圳硅基智能科技有限公司 Artificial neural network and system for recognizing eye fundus image lesion
CN107145756A (en) * 2017-05-17 2017-09-08 上海辉明软件有限公司 A kind of stroke types Forecasting Methodology and device
CN107679525A (en) * 2017-11-01 2018-02-09 腾讯科技(深圳)有限公司 Image classification method, device and computer-readable recording medium
CN108734209A (en) * 2018-05-16 2018-11-02 上海鹰瞳医疗科技有限公司 Feature recognition based on more images and equipment

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115311188A (en) * 2021-05-08 2022-11-08 数坤(北京)网络科技股份有限公司 Image identification method and device, electronic equipment and storage medium
CN115311188B (en) * 2021-05-08 2023-12-22 数坤科技股份有限公司 Image recognition method and device, electronic equipment and storage medium

Also Published As

Publication number Publication date
CN109671049A (en) 2019-04-23
CN109671049B (en) 2024-03-01

Similar Documents

Publication Publication Date Title
Li et al. Iternet: Retinal image segmentation utilizing structural redundancy in vessel networks
Li et al. A large-scale database and a CNN model for attention-based glaucoma detection
Asiri et al. Deep learning based computer-aided diagnosis systems for diabetic retinopathy: A survey
WO2020093563A1 (en) Medical image processing method, system, device, and storage medium
Tabassum et al. CDED-Net: Joint segmentation of optic disc and optic cup for glaucoma screening
Veena et al. A novel optic disc and optic cup segmentation technique to diagnose glaucoma using deep learning convolutional neural network over retinal fundus images
Sheng et al. Retinal vessel segmentation using minimum spanning superpixel tree detector
Akram et al. Automated detection of dark and bright lesions in retinal images for early detection of diabetic retinopathy
Xiang et al. Automatic segmentation of retinal layer in OCT images with choroidal neovascularization
Abràmoff et al. Retinal imaging and image analysis
Bilal et al. Survey on recent developments in automatic detection of diabetic retinopathy
Li et al. DeepRetina: layer segmentation of retina in OCT images using deep learning
Akil et al. Detection of retinal abnormalities in fundus image using CNN deep learning networks
Zhu et al. Digital image processing for ophthalmology: Detection of the optic nerve head
Bhati et al. Discriminative kernel convolution network for multi-label ophthalmic disease detection on imbalanced fundus image dataset
Suriyasekeran et al. Algorithms for diagnosis of diabetic retinopathy and diabetic macula edema-a review
Qin et al. Optic disc and cup segmentation based on deep learning
Wang et al. Accurate disease detection quantification of iris based retinal images using random implication image classifier technique
Karkuzhali et al. Robust intensity variation and inverse surface adaptive thresholding techniques for detection of optic disc and exudates in retinal fundus images
Sangeethaa Presumptive discerning of the severity level of glaucoma through clinical fundus images using hybrid PolyNet
Rahhal et al. Detection and classification of diabetic retinopathy using artificial intelligence algorithms
Samawi et al. Optic disc segmentation in retinal fundus images using morphological techniques and intensity thresholding
Engberg et al. Automated quantification of retinal microvasculature from OCT angiography using dictionary-based vessel segmentation
Khan et al. Screening fundus images to extract multiple ocular features: A unified modeling approach
Akram et al. Retinal image blood vessel segmentation

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 18939157

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 18939157

Country of ref document: EP

Kind code of ref document: A1

32PN Ep: public notification in the ep bulletin as address of the adressee cannot be established

Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 112(1) EPC (EPO FORM 1205A DATED 26.01.2022)

122 Ep: pct application non-entry in european phase

Ref document number: 18939157

Country of ref document: EP

Kind code of ref document: A1