Automatic measurement of ISNT and CDR on retinal images by means of a fast and efficient method based on mathematical morphology and active contours Rafael Verdú-Monedero 1 , Juan Morales-Sánchez 1 , Rafael Berenguer-Vidal 2 , Inmaculada Sellés-Navarro 3 , Ana Palazón-Cabanes 3 1 Universidad Politécnica de Cartagena, Cartagena 30202, España 2 Universidad Católica San Antonio, Murcia 30107, Spain 3 Hospital Universitario Reina Sofía, Murcia 30003, Spain 8th International Work-Conference on the Interplay between Natural and Artificial Computation - IWINAC 2019
Index Introduction 1 Method 2 Image acquisition and enhancement Automatic localization by morphological operations Accurate delineation by Parametric Active Contours Results 3 Conclusions 4 Verdú et al. (UPCT-UCAM-HGURS) IWINAC 2019 Thursday, 06.06.2019 2 / 25
Index Introduction 1 Method 2 Image acquisition and enhancement Automatic localization by morphological operations Accurate delineation by Parametric Active Contours Results 3 Conclusions 4 Verdú et al. (UPCT-UCAM-HGURS) IWINAC 2019 Thursday, 06.06.2019 3 / 25
Introduction Spanish National projects - Instituto de Salud Carlos III Subproject AES2017-PI17/00771 Extracción automática de características del nervio óptico de ambos ojos mediante procesado de imagen y reconocimiento de patrones para su integración en aplicación CAD de telemedicina Universidad Politécnica de Cartagena. Subproject AES2017-PI17-00821 Desarrollo y evaluación de una aplicación de telemedicina para el análisis anatómico conjunto del nervio óptico de ambos ojos Hospital Universitario Reina Sofía. Verdú et al. (UPCT-UCAM-HGURS) IWINAC 2019 Thursday, 06.06.2019 4 / 25
Introduction About Glaucoma It is undoubtedly the most prevalent disease of optic neuropathies. This clinically silent ocular disease can lead to a progressive deterioration of the visual field causing eventually total blindness. Ocular anatomical alterations appear several years before the characteristic and irreversible lesions in the visual field, as a consequence of the death of the retina ganglion cells. Initial clinical explorations that may lead to suspected glaucoma intraocular pressure levels, and the appearance of the optic nerve. In these cases, other complementary test as the assessment of the retinal fibers layer and the study of the visual field should confirm the diagnosis. Verdú et al. (UPCT-UCAM-HGURS) IWINAC 2019 Thursday, 06.06.2019 5 / 25
Introduction About Glaucoma Verdú et al. (UPCT-UCAM-HGURS) IWINAC 2019 Thursday, 06.06.2019 6 / 25
Introduction About Glaucoma Some anatomical parameters of the papilla can be used as reference indices in the early detection of glaucoma: size and shape of the optic nerve (ISNT rule), the relationship between the size of the papillary excavation and the neuroretinal ring (cup-to-disc ratio, CDR), the configuration of the excavation depth, the exit position of the papillary vessels, the presence of peripapillary hemorrhages, the nerve fiber defect or the chorioretinal atrophy, ... One of the drawbacks is the variability in size and appearance of the papilla in the normal population. Verdú et al. (UPCT-UCAM-HGURS) IWINAC 2019 Thursday, 06.06.2019 7 / 25
Introduction Reference indices based on anatomical parameters: ISNT rule To differentiate normal from glaucomatous eyes in the clinical evaluation of the optic nerve head, normal eyes follow a characteristic configuration for the thickness of the disc rim: I ≥ S ≥ N ≥ T Eyes that deviate from the ISNT rule may need close monitoring for glaucoma. Verdú et al. (UPCT-UCAM-HGURS) IWINAC 2019 Thursday, 06.06.2019 8 / 25
Introduction Reference indices based on anatomical parameters: ISNT rule The utility of this rule as a standalone criterion is not demonstrated for the early diagnosis of glaucoma. The predictability of this rule improves when it is complemented with other measurements and exploration data, such as, i.e., the cup-to-disc ratio. The presence of asymmetries in the papillary excavation between both eyes may be a suggestive factor of suffering from glaucoma. Verdú et al. (UPCT-UCAM-HGURS) IWINAC 2019 Thursday, 06.06.2019 8 / 25
Introduction Reference indices based on anatomical parameters: CDR The Cup-to-disc ratio, CDR, compares the diameter of the cup portion with the total diameter of the optic disc. Normal values of CDR are 0.3-0.4. The CDR is commonly used in ophthalmology as a measurement of the evolution of glaucoma. As glaucoma advances, the cup enlarges until it occupies most of the disc area. Verdú et al. (UPCT-UCAM-HGURS) IWINAC 2019 Thursday, 06.06.2019 9 / 25
Index Introduction 1 Method 2 Image acquisition and enhancement Automatic localization by morphological operations Accurate delineation by Parametric Active Contours Results 3 Conclusions 4 Verdú et al. (UPCT-UCAM-HGURS) IWINAC 2019 Thursday, 06.06.2019 10 / 25
Method Steps of the method Image acquisition. 1 Preprocessing: Enhancement of the dynamic range and slight 2 correction of the illumination. Automatic localization of the optic disk and excavation by means 3 of morphological operations. Accurate delineation of the optic disk and the excavation by two 4 active contours. Measurements: ISNT and CDR. 5 Verdú et al. (UPCT-UCAM-HGURS) IWINAC 2019 Thursday, 06.06.2019 11 / 25
Method RGB channels of the color image Original color image Red channel Green channel Blue channel Verdú et al. (UPCT-UCAM-HGURS) IWINAC 2019 Thursday, 06.06.2019 12 / 25
Method 1. Image acquisition The fundus images were acquired by the Ophthalmology Service of the Hospital General Universitario Reina Sofía (Murcia, Spain) by means of a Topcon TRC-NW400 non-mydriatic retinal camera. These color images (RGB and 8 bit/channel) are in DICOM format, with a picture angle of 30 ◦ and an initial resolution of 1934 × 2576 pixels (downsampled to 600 × 800 pixels before processing). 2. Image enhancement Contrast stretching for a full widening of the dynamic range. Smooth local and adaptive histogram equalization for homogenizing the illumination of the images. Verdú et al. (UPCT-UCAM-HGURS) IWINAC 2019 Thursday, 06.06.2019 13 / 25
Method 3. a) Automatic location of optic disc by morphological operators Red channel of the image, I R 0 . Closing and opening with a disk-shaped se R = 11 : I R 1 = (( I R 0 • D 11 ) ◦ D 11 ) Top-hat using a R = 99 disk-shaped se: I R 2 = ( I R 1 − ( I R 1 ◦ D 99 )) Threshold: M OD = thres ( I R 2 ) Verdú et al. (UPCT-UCAM-HGURS) IWINAC 2019 Thursday, 06.06.2019 14 / 25
Method 3. a) Automatic location of optic disc by morphological operators Red channel of the image, I R 0 . Closing and opening with a disk-shaped se R = 11 : I R 1 = (( I R 0 • D 11 ) ◦ D 11 ) Top-hat using a R = 99 disk-shaped se: I R 2 = ( I R 1 − ( I R 1 ◦ D 99 )) Threshold: M OD = thres ( I R 2 ) Verdú et al. (UPCT-UCAM-HGURS) IWINAC 2019 Thursday, 06.06.2019 14 / 25
Method 3. a) Automatic location of optic disc by morphological operators Red channel of the image, I R 0 . Closing and opening with a disk-shaped se R = 11 : I R 1 = (( I R 0 • D 11 ) ◦ D 11 ) Top-hat using a R = 99 disk-shaped se: I R 2 = ( I R 1 − ( I R 1 ◦ D 99 )) Threshold: M OD = thres ( I R 2 ) Verdú et al. (UPCT-UCAM-HGURS) IWINAC 2019 Thursday, 06.06.2019 14 / 25
Method 3. a) Automatic location of optic disc by morphological operators Red channel of the image, I R 0 . Closing and opening with a disk-shaped se R = 11 : I R 1 = (( I R 0 • D 11 ) ◦ D 11 ) Top-hat using a R = 99 disk-shaped se: I R 2 = ( I R 1 − ( I R 1 ◦ D 99 )) Threshold: M OD = thres ( I R 2 ) Verdú et al. (UPCT-UCAM-HGURS) IWINAC 2019 Thursday, 06.06.2019 14 / 25
Method 3. a) Automatic location of optic disc by morphological operators Red channel of the image, I R 0 . Closing and opening with a disk-shaped se R = 11 : I R 1 = (( I R 0 • D 11 ) ◦ D 11 ) Top-hat using a R = 99 disk-shaped se: I R 2 = ( I R 1 − ( I R 1 ◦ D 99 )) Threshold: M OD = thres ( I R 2 ) Verdú et al. (UPCT-UCAM-HGURS) IWINAC 2019 Thursday, 06.06.2019 14 / 25
Method 3. b) Automatic location of excavation by morphological operators Blue channel of the image, I B 0 . ROI from the optic disc, I B 1 . Closing and opening I B 2 = (( I B 1 • D 11 ) ◦ D 11 ) . Top-hat: I B 3 = ( I B 2 − ( I B 2 ◦ D 51 )) Threshold: M E = thres ( I B 3 ) Verdú et al. (UPCT-UCAM-HGURS) IWINAC 2019 Thursday, 06.06.2019 15 / 25
Method 3. b) Automatic location of excavation by morphological operators Blue channel of the image, I B 0 . ROI from the optic disc, I B 1 . Closing and opening I B 2 = (( I B 1 • D 11 ) ◦ D 11 ) . Top-hat: I B 3 = ( I B 2 − ( I B 2 ◦ D 51 )) Threshold: M E = thres ( I B 3 ) Verdú et al. (UPCT-UCAM-HGURS) IWINAC 2019 Thursday, 06.06.2019 15 / 25
Method 3. b) Automatic location of excavation by morphological operators Blue channel of the image, I B 0 . ROI from the optic disc, I B 1 . Closing and opening I B 2 = (( I B 1 • D 11 ) ◦ D 11 ) . Top-hat: I B 3 = ( I B 2 − ( I B 2 ◦ D 51 )) Threshold: M E = thres ( I B 3 ) Verdú et al. (UPCT-UCAM-HGURS) IWINAC 2019 Thursday, 06.06.2019 15 / 25
Recommend
More recommend
