CN111640090A - Method for evaluating fundus image quality - Google Patents
Method for evaluating fundus image quality Download PDFInfo
- Publication number
- CN111640090A CN111640090A CN202010399562.3A CN202010399562A CN111640090A CN 111640090 A CN111640090 A CN 111640090A CN 202010399562 A CN202010399562 A CN 202010399562A CN 111640090 A CN111640090 A CN 111640090A
- Authority
- CN
- China
- Prior art keywords
- image
- retina
- fundus
- evaluating
- contrast
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/194—Segmentation; Edge detection involving foreground-background segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10024—Color image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20056—Discrete and fast Fourier transform, [DFT, FFT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20112—Image segmentation details
- G06T2207/20132—Image cropping
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30041—Eye; Retina; Ophthalmic
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30168—Image quality inspection
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Quality & Reliability (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Eye Examination Apparatus (AREA)
Abstract
The invention belongs to the field of medical image processing, and discloses a method for evaluating the quality of fundus images. Preprocessing an acquired fundus image, cutting redundant background around a retina image, and obtaining an area only containing a retina; then, color, focus, contrast, and illumination characteristics are respectively extracted based on the preprocessed fundus images and evaluated; finally, based on the confidence of the above-described feature weighted evaluation, the evaluation result of the fundus image is determined, and the cause of low imaging quality is analyzed. Thus, the method can be used to guide the retinal image acquisition process, helping to reduce the workload of ophthalmologists in screening and improving cost effectiveness.
Description
Technical Field
The invention relates to the field of medical image processing, in particular to a method for evaluating the quality of fundus images.
Background
With the development of socioeconomic development and the change of diet and living habits of Chinese people, the incidence of diabetes is on the trend of increasing year by year, and at present, China becomes the world with the largest number of diabetics. Among them, diabetic retinopathy is one of the common diabetic complications, the incidence rate of which is more than 40% in diabetic patients, and diabetic retinopathy has become one of the main causes of visual deterioration and even blindness. However, by early diagnosis and treatment of diabetic retinopathy patients, visual impairment and blindness can be prevented in more than 50% of patients. Fundus screening is therefore critical for the early diagnosis of diabetic retinopathy and the prevention of blindness.
However, the fundus camera sometimes has poor imaging quality due to defocus or insufficient illumination, and further affects subsequent image analysis and clinical diagnosis. Therefore, before further analyzing the fundus image, the fundus image quality needs to be evaluated.
Disclosure of Invention
In view of the above problems, the present invention aims to provide a method for evaluating fundus image quality. The method carries out weighted evaluation on the fundus image quality aiming at the reasons which are easy to cause fundus image quality low, such as color distortion, blurring, low contrast ratio, insufficient or too strong illumination and the like.
In order to achieve the above object, the technical solution of the present invention is as follows:
a method for evaluating the quality of eyeground image, said method comprises carrying on the preconditioning to the examinee's eyeground image that the eyeground camera gathers, cut out the redundant background around the retinal image, obtain the area comprising retina only; then, based on the preprocessed fundus images, color, focus, contrast and illumination features are respectively extracted and evaluated; and finally, performing weighted evaluation on the characteristics, calculating the confidence coefficient of the characteristics, giving an evaluation result of the fundus image and reasons for low imaging quality, and guiding the next fundus image acquisition process.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
Detailed Description
The following detailed description of specific embodiments of the present invention is provided in connection with the accompanying drawings.
A method for evaluating the quality of an image of a fundus of the eye, comprising in particular the steps of:
step 1, preprocessing a fundus image of a detected person acquired by a fundus camera, cutting redundant background around the fundus image, and obtaining a region only containing retina;
in the present embodiment, step 1 converts the retina image into a grayscale image. And obtaining a mask image by using a threshold segmentation algorithm, searching the maximum radius as a circular area in the foreground area by taking the central point of the image as the center of a circle, wherein the circular area comprises all foreground areas, and finally cutting the circular area to obtain a retina image.
Step 2, respectively extracting color, focus, contrast and illumination characteristics based on the preprocessed fundus image and evaluating the color, focus, contrast and illumination characteristics;
in the present embodiment, in the color feature evaluation process in step 2, the cut RGB color image of the retina is converted into HSV color space, and a single threshold T is set on the H channelcAnd evaluating whether the color of the retina image is distorted. When H channel is larger than threshold value TcWhen the color of the retina image is not distorted, when the H channel is less than the threshold value TcThe retinal image color is distorted.
In the focusing characteristic evaluation process in the step 2, firstly, blood vessel and optic disc characteristics of the preprocessed retina image are extracted, the blood vessel and the optic disc are divided, the density of the blood vessel and the optic disc is obtained through calculation, and whether the retina image is focused or not is determined through the density of the blood vessel and the optic disc. Wherein the blood vessel and optic disc density is obtained by the ratio of the segmented blood vessel and optic disc area to the area of the cropped retina image.
In the contrast characteristic evaluation process in the step 2, firstly, the retina image of the standard database is selected as a reference retina image, and the difference value between the reference retina image and the retina image to be tested is calculated to obtain an error retina image. The spatial domain of the error retinal image is then transformed to the frequency domain using a fast fourier transform. Then, contrast sensitivity function weighted modulation is applied to the frequency domain coefficients, and the frequency domain is transformed to a space domain by using fast Fourier transform. Finally, the root mean square error of the spatial error image is calculated, the lower the root mean square error value, the higher the contrast of the retinal image. And setting a threshold value for the root mean square value of the calculated spatial error image, wherein if the threshold value is larger than the threshold value, the retinal image contrast is low, and if the threshold value is smaller than the threshold value, the retinal image contrast is high.
In the illumination characteristic evaluation process in the step 2, the retina RGB color image after pretreatment is converted into HSV color space, and two thresholds T are set on the V channelH1And TH2To assess whether the retinal image illumination is uniform. When V channel is less than threshold TH1When the retina image is not sufficiently exposed, the retina image is not sufficiently exposed; when V channel is larger than threshold value TH2When the retina image is over-exposed; when V channel is at threshold TH1And TH2In between, the retinal image is uniformly illuminated.
Step 3, evaluating the characteristics of color, focus, contrast, illumination and the like according to w respectively1,w2,w3,w4The weights of (a) and (b) are weighted, and a quality evaluation result of the fundus image and a cause of low quality of the fundus image are given based on the confidence of the feature weighted evaluation.
Claims (7)
1. A method for evaluating the quality of a fundus image, characterized in that features of the fundus image, such as color, focus, contrast and illumination, are weighted and evaluated, and based on the confidence of the weighted evaluation of the features, an evaluation result of the fundus image is given, comprising the steps of:
preprocessing a fundus image of a detected person acquired by a fundus camera, cutting redundant background around the fundus image, and obtaining an area only containing retinas;
secondly, respectively extracting color, focus, contrast and illumination characteristics based on the preprocessed fundus image and evaluating the color, the focus, the contrast and the illumination characteristics;
and step three, performing weighted evaluation on the characteristics, calculating the confidence coefficient of the characteristics and giving an evaluation result of the fundus image.
2. The preprocessing procedure in a method for evaluating quality of a fundus image according to claim 1, wherein said fundus image is converted into a gray scale image, and a mask image is obtained by using a threshold segmentation algorithm, a circular area is found as a maximum radius from a center point of the image in a foreground area, the circular area includes all foreground areas, and finally the retinal image is obtained by cutting out the circular area.
3. The color feature evaluation process in a method for evaluating fundus image quality according to claim 1 or 2, wherein the pre-processed retina RGB color image is converted into HSV color space, and a single threshold is set on an H-channel to evaluate whether the retina image color is distorted; further, the retinal image color is not distorted when the H channel is greater than the threshold, and is distorted when the H channel is less than the threshold.
4. The focus characteristic evaluation process in a method for evaluating fundus image quality according to claim 1 or 2, characterized in that blood vessel and optic disc characteristics of the pre-processed retinal image are first extracted, the blood vessel and optic disc are divided, blood vessel and optic disc densities are calculated, and whether the retinal image is focused or not is judged by the blood vessel and optic disc densities; further, the blood vessel and optic disc densities are the ratio of the segmented blood vessel and optic disc areas to the area of the cropped retina image, respectively.
5. The contrast characteristic evaluation process in a method for evaluating the quality of a fundus image according to claim 1 or 2, wherein a retina image of a standard database is first selected as a reference retina image, a difference between the reference retina image and a retina image to be tested is calculated, and an error retina image is obtained; then transforming the error retina image from a space domain to a frequency domain by using fast Fourier transform; then applying contrast sensitivity function weighted modulation to the frequency domain coefficient, and transforming the frequency domain coefficient to a space domain by using inverse fast Fourier transform to obtain a space error image; finally, calculating the root mean square error of the spatial error image, wherein the lower the root mean square error value is, the higher the contrast of the representative retina image is; furthermore, a threshold value is set for the root mean square value of the spatial error image, the retinal image contrast is low when the threshold value is larger than the threshold value, and the retinal image contrast is high when the threshold value is smaller than the threshold value.
6. The illumination characteristic evaluation process in a method for evaluating fundus image quality according to claim 1 or 2, wherein the pre-processed retina RGB color image is converted into HSV color space, and two thresholds T are set on the V-channelH1And TH2To assess whether the retinal image illumination is uniform; further, when the V channel is less than the threshold TH1When the V channel is larger than the threshold value T, the retina image is underexposedH2When the retinal image is overexposed, when the V channel is at the threshold value TH1And TH2In between, the retinal image is uniformly illuminated.
7. A method for evaluating the quality of a fundus image according to claim 2 or 3, characterized in that the results of the evaluation of the characteristics of its colour, focus, contrast and illumination are according to w respectively1,w2,w3,w4The weighting is carried out, and finally a global retinal image quality evaluation result is obtained.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010399562.3A CN111640090B (en) | 2020-05-12 | 2020-05-12 | Method for evaluating quality of fundus images |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010399562.3A CN111640090B (en) | 2020-05-12 | 2020-05-12 | Method for evaluating quality of fundus images |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111640090A true CN111640090A (en) | 2020-09-08 |
CN111640090B CN111640090B (en) | 2022-05-31 |
Family
ID=72330045
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010399562.3A Active CN111640090B (en) | 2020-05-12 | 2020-05-12 | Method for evaluating quality of fundus images |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111640090B (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113037998A (en) * | 2021-02-04 | 2021-06-25 | 浙江大华技术股份有限公司 | Focusing method, camera device and storage medium |
CN114299037A (en) * | 2021-12-30 | 2022-04-08 | 广州极飞科技股份有限公司 | Method and device for evaluating quality of object detection result, electronic equipment and computer readable storage medium |
CN115187519A (en) * | 2022-06-21 | 2022-10-14 | 上海市计量测试技术研究院 | Image quality evaluation method, system and computer readable medium |
CN117041531A (en) * | 2023-09-04 | 2023-11-10 | 无锡维凯科技有限公司 | Mobile phone camera focusing detection method and system based on image quality evaluation |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018116321A2 (en) * | 2016-12-21 | 2018-06-28 | Braviithi Technologies Private Limited | Retinal fundus image processing method |
CN108921191A (en) * | 2018-05-25 | 2018-11-30 | 北方工业大学 | Multi-biological-feature fusion recognition method based on image quality evaluation |
CN109544540A (en) * | 2018-11-28 | 2019-03-29 | 东北大学 | A kind of diabetic retina picture quality detection method based on image analysis technology |
CN109829877A (en) * | 2018-09-20 | 2019-05-31 | 中南大学 | A kind of retinal fundus images cup disc ratio automatic evaluation method |
CN110009631A (en) * | 2019-04-15 | 2019-07-12 | 唐晓颖 | Vascular quality appraisal procedure, device, equipment and the medium of eye fundus image |
WO2019246098A1 (en) * | 2018-06-18 | 2019-12-26 | New Jersey Institute Of Technology | Method, system and apparatus for diagnostic assessment and screening of binocular dysfunctions |
CN111127298A (en) * | 2019-06-12 | 2020-05-08 | 上海大学 | Panoramic image blind quality assessment method |
CN111127425A (en) * | 2019-12-23 | 2020-05-08 | 北京至真互联网技术有限公司 | Target detection positioning method and device based on retina fundus image |
-
2020
- 2020-05-12 CN CN202010399562.3A patent/CN111640090B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018116321A2 (en) * | 2016-12-21 | 2018-06-28 | Braviithi Technologies Private Limited | Retinal fundus image processing method |
CN108921191A (en) * | 2018-05-25 | 2018-11-30 | 北方工业大学 | Multi-biological-feature fusion recognition method based on image quality evaluation |
WO2019246098A1 (en) * | 2018-06-18 | 2019-12-26 | New Jersey Institute Of Technology | Method, system and apparatus for diagnostic assessment and screening of binocular dysfunctions |
CN109829877A (en) * | 2018-09-20 | 2019-05-31 | 中南大学 | A kind of retinal fundus images cup disc ratio automatic evaluation method |
CN109544540A (en) * | 2018-11-28 | 2019-03-29 | 东北大学 | A kind of diabetic retina picture quality detection method based on image analysis technology |
CN110009631A (en) * | 2019-04-15 | 2019-07-12 | 唐晓颖 | Vascular quality appraisal procedure, device, equipment and the medium of eye fundus image |
CN111127298A (en) * | 2019-06-12 | 2020-05-08 | 上海大学 | Panoramic image blind quality assessment method |
CN111127425A (en) * | 2019-12-23 | 2020-05-08 | 北京至真互联网技术有限公司 | Target detection positioning method and device based on retina fundus image |
Non-Patent Citations (1)
Title |
---|
张悦等: "基于深度学习辅助诊断青光眼病灶检测算法及应用", 《眼科》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113037998A (en) * | 2021-02-04 | 2021-06-25 | 浙江大华技术股份有限公司 | Focusing method, camera device and storage medium |
CN113037998B (en) * | 2021-02-04 | 2023-03-24 | 浙江大华技术股份有限公司 | Focusing method, camera device and storage medium |
CN114299037A (en) * | 2021-12-30 | 2022-04-08 | 广州极飞科技股份有限公司 | Method and device for evaluating quality of object detection result, electronic equipment and computer readable storage medium |
CN114299037B (en) * | 2021-12-30 | 2023-09-01 | 广州极飞科技股份有限公司 | Quality evaluation method and device for object detection result, electronic equipment and computer readable storage medium |
CN115187519A (en) * | 2022-06-21 | 2022-10-14 | 上海市计量测试技术研究院 | Image quality evaluation method, system and computer readable medium |
CN117041531A (en) * | 2023-09-04 | 2023-11-10 | 无锡维凯科技有限公司 | Mobile phone camera focusing detection method and system based on image quality evaluation |
CN117041531B (en) * | 2023-09-04 | 2024-03-15 | 无锡维凯科技有限公司 | Mobile phone camera focusing detection method and system based on image quality evaluation |
Also Published As
Publication number | Publication date |
---|---|
CN111640090B (en) | 2022-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111640090B (en) | Method for evaluating quality of fundus images | |
CN110400289B (en) | Fundus image recognition method, fundus image recognition device, fundus image recognition apparatus, and fundus image recognition storage medium | |
CN107451998B (en) | Fundus image quality control method | |
Alwazzan et al. | A hybrid algorithm to enhance colour retinal fundus images using a Wiener filter and CLAHE | |
CN109472781B (en) | Diabetic retinopathy detection system based on serial structure segmentation | |
CN108272434B (en) | Method and device for processing fundus images | |
CN109948719B (en) | Automatic fundus image quality classification method based on residual dense module network structure | |
Esmaeili et al. | A new curvelet transform based method for extraction of red lesions in digital color retinal images | |
CN111986211A (en) | Deep learning-based ophthalmic ultrasonic automatic screening method and system | |
KR102155381B1 (en) | Method, apparatus and software program for cervical cancer decision using image analysis of artificial intelligence based technology | |
Yadav et al. | Comparative analysis of fundus image enhancement in detection of diabetic retinopathy | |
WO2010030159A2 (en) | A non invasive method for analysing the retina for ocular manifested diseases | |
Li et al. | Vessel recognition of retinal fundus images based on fully convolutional network | |
Subramanian et al. | Diabetic Retinopathy–Feature Extraction and Classification using Adaptive Super Pixel Algorithm | |
CN111899267A (en) | Retina blood vessel segmentation algorithm based on level set | |
KR20210124024A (en) | Method, apparatus and computer program for calculating quality score of fundus image data using artificial intelligence | |
Bhardwaj et al. | Appraisal of pre-processing techniques for automated detection of diabetic retinopathy | |
Zhao et al. | Attention residual convolution neural network based on U-net (AttentionResU-Net) for retina vessel segmentation | |
Dias et al. | Evaluation of retinal image gradability by image features classification | |
Atheesan et al. | Automatic glaucoma detection by using funduscopic images | |
Chakour et al. | Blood vessel segmentation of retinal fundus images using dynamic preprocessing and mathematical morphology | |
Giancardo et al. | Quality assessment of retinal fundus images using elliptical local vessel density | |
Umamageswari et al. | A novel approach for classification of diabetics from retinal image using deep learning technique | |
CN111292285B (en) | Automatic screening method for diabetes mellitus based on naive Bayes and support vector machine | |
CN114972148A (en) | Fundus image quality evaluation method, system and device |
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 |