Retinal Image Analysis and its use in Medical Applications

Retina located in the back of the eye is not only a vital part of human sight, but also contains valuable information that can be used in biometric security applications, or for the diagnosis of certain diseases. In order to analyze this information from retinal images, its features of blood vessels, microaneurysms and the optic disc require extraction and detection respectively. We propose a method to extract vessels called MF-FDOG. MF-FDOG consists of using two filters, Matched Filter (MF) and the first-order derivative of Gaussian (FDOG). The vessel map is extracted by applying a threshold to the response of MF, which is adaptively adjusted by the mean response of FDOG. This method allows us to better distinguish vessel objects from non-vessel objects. Microaneurysm (MA) detection is accomplished with two proposed algorithms, Multi-scale Correlation Filtering (MSCF) and Dictionary Learning (DL) with Sparse Representation Classifier (SRC). MSCF is hierarchical in nature, consisting of two levels: coarse level microaneurysm candidate detection and fine level true microaneurysm detection. In the first level, all possible microaneurysm candidates are found while the second level extracts features from each candidate and compares them to a discrimination table for decision (MA or non-MA). In Dictionary Learning with Sparse Representation Classifier, MA and non-MA objects are extracted from images and used to learn two dictionaries, MA and non-MA. Sparse Representation Classifier is then applied to each MA candidate object detected beforehand, using the two dictionaries to determine class membership. The detection result is further improved by adding a class discrimination term into the Dictionary Learning model. This approach is known as Centralized Dictionary Learning (CDL) with Sparse Representation Classifier. The optic disc (OD) is an important anatomical feature in retinal images, and its detection is vital for developing automated screening programs. Currently, there is no algorithm designed to automatically detect the OD in fundus images captured from Asians, which are larger and have thicker vessels compared to Caucasians. We propose such a method to complement current algorithms using two steps: OD vessel candidate detection and OD vessel candidate matching. The proposed extraction/detection approaches are tested in medical applications, specifically the case study of detecting diabetic retinopathy (DR). DR is a complication of diabetes that damages the retina and can lead to blindness. There are four stages of DR and is a leading cause of sight loss in industrialized nations. Using MF-FDOG, blood vessels were extracted from DR images, while DR images fed into MSCF and Dictionary and Centralized Dictionary Learning with Sparse Representation Classifier produced good microaneurysm detection results. Using a new database consisting of only Asian DR patients, we successfully tested our OD detection method. As part of future work we intend to improve existing methods such as enhancing low contrast microaneurysms and better scale selection. In additional, we will extract other features from the retina, develop a generalized OD detection method, apply Dictionary Learning with Sparse Representation Classifier to vessel extraction, and use the new image database to carry out more experiments in medical applications

A Rule Based Segmentation Approaches to Extract Retinal Blood Vessels in Fundus Image

The physiological structures of the retinal blood vessel are one of the key features that visible in the retinal images and contain the information associate with the anatomical abnormalities. It is accepted all over the world to judge the cardiovascular and retinal disease. To avoid the risk of visual impairment, appropriate vessel segmentation is mandatory. Here has proposed a segmentation algorithm that efficiently extracts the blood vessels from the retinal fundus image. The proposed segmentation algorithm is performed Lab and Principle Component (PC) based gray level conversion, Contrast Limited Adaptive Histogram Equalization (CLAHE), morphological operations, Local Property-Based Pixel Correction (LPBPC). For appropriate detection proposed vessels correction algorithm LPBPC that check the feature of the vessels and remove the wrong vessel detection. To measure the appropriateness of the proposed algorithm, the experimental results are compared with the corresponding ground truth images. The experimental results have shown that the proposed blood vessel algorithm is more accurate than the existing algorithms

Analysis of Retinal Image Data to Support Glaucoma Diagnosis

Fundus kamera je široce dostupné zobrazovací zařízení, které umožňuje relativně rychlé a nenákladné vyšetření zadního segmentu oka – sítnice. Z těchto důvodů se mnoho výzkumných pracovišť zaměřuje právě na vývoj automatických metod diagnostiky nemocí sítnice s využitím fundus fotografií. Tato dizertační práce analyzuje současný stav vědeckého poznání v oblasti diagnostiky glaukomu s využitím fundus kamery a navrhuje novou metodiku hodnocení vrstvy nervových vláken (VNV) na sítnici pomocí texturní analýzy. Spolu s touto metodikou je navržena metoda segmentace cévního řečiště sítnice, jakožto další hodnotný příspěvek k současnému stavu řešené problematiky. Segmentace cévního řečiště rovněž slouží jako nezbytný krok předcházející analýzu VNV. Vedle toho práce publikuje novou volně dostupnou databázi snímků sítnice se zlatými standardy pro účely hodnocení automatických metod segmentace cévního řečiště.Fundus camera is widely available imaging device enabling fast and cheap examination of the human retina. Hence, many researchers focus on development of automatic methods towards assessment of various retinal diseases via fundus images. This dissertation summarizes recent state-of-the-art in the field of glaucoma diagnosis using fundus camera and proposes a novel methodology for assessment of the retinal nerve fiber layer (RNFL) via texture analysis. Along with it, a method for the retinal blood vessel segmentation is introduced as an additional valuable contribution to the recent state-of-the-art in the field of retinal image processing. Segmentation of the blood vessels also serves as a necessary step preceding evaluation of the RNFL via the proposed methodology. In addition, a new publicly available high-resolution retinal image database with gold standard data is introduced as a novel opportunity for other researches to evaluate their segmentation algorithms.

Retinal vessel segmentation using textons

Segmenting vessels from retinal images, like segmentation in many other medical image domains, is a challenging task, as there is no unified way that can be adopted to extract the vessels accurately. However, it is the most critical stage in automatic assessment of various forms of diseases (e.g. Glaucoma, Age-related macular degeneration, diabetic retinopathy and cardiovascular diseases etc.). Our research aims to investigate retinal image segmentation approaches based on textons as they provide a compact description of texture that can be learnt from a training set. This thesis presents a brief review of those diseases and also includes their current situations, future trends and techniques used for their automatic diagnosis in routine clinical applications. The importance of retinal vessel segmentation is particularly emphasized in such applications. An extensive review of previous work on retinal vessel segmentation and salient texture analysis methods is presented. Five automatic retinal vessel segmentation methods are proposed in this thesis. The first method focuses on addressing the problem of removing pathological anomalies (Drusen, exudates) for retinal vessel segmentation, which have been identified by other researchers as a problem and a common source of error. The results show that the modified method shows some improvement compared to a previously published method. The second novel supervised segmentation method employs textons. We propose a new filter bank (MR11) that includes bar detectors for vascular feature extraction and other kernels to detect edges and photometric variations in the image. The k-means clustering algorithm is adopted for texton generation based on the vessel and non-vessel elements which are identified by ground truth. The third improved supervised method is developed based on the second one, in which textons are generated by k-means clustering and texton maps representing vessels are derived by back projecting pixel clusters onto hand labelled ground truth. A further step is implemented to ensure that the best combinations of textons are represented in the map and subsequently used to identify vessels in the test set. The experimental results on two benchmark datasets show that our proposed method performs well compared to other published work and the results of human experts. A further test of our system on an independent set of optical fundus images verified its consistent performance. The statistical analysis on experimental results also reveals that it is possible to train unified textons for retinal vessel segmentation. In the fourth method a novel scheme using Gabor filter bank for vessel feature extraction is proposed. The ii method is inspired by the human visual system. Machine learning is used to optimize the Gabor filter parameters. The experimental results demonstrate that our method significantly enhances the true positive rate while maintaining a level of specificity that is comparable with other approaches. Finally, we proposed a new unsupervised texton based retinal vessel segmentation method using derivative of SIFT and multi-scale Gabor filers. The lack of sufficient quantities of hand labelled ground truth and the high level of variability in ground truth labels amongst experts provides the motivation for this approach. The evaluation results reveal that our unsupervised segmentation method is comparable with the best other supervised methods and other best state of the art methods

Human treelike tubular structure segmentation: A comprehensive review and future perspectives

Various structures in human physiology follow a treelike morphology, which often expresses complexity at very fine scales. Examples of such structures are intrathoracic airways, retinal blood vessels, and hepatic blood vessels. Large collections of 2D and 3D images have been made available by medical imaging modalities such as magnetic resonance imaging (MRI), computed tomography (CT), Optical coherence tomography (OCT) and ultrasound in which the spatial arrangement can be observed. Segmentation of these structures in medical imaging is of great importance since the analysis of the structure provides insights into disease diagnosis, treatment planning, and prognosis. Manually labelling extensive data by radiologists is often time-consuming and error-prone. As a result, automated or semi-automated computational models have become a popular research field of medical imaging in the past two decades, and many have been developed to date. In this survey, we aim to provide a comprehensive review of currently publicly available datasets, segmentation algorithms, and evaluation metrics. In addition, current challenges and future research directions are discussed

A Multi-Anatomical Retinal Structure Segmentation System For Automatic Eye Screening Using Morphological Adaptive Fuzzy Thresholding

Eye exam can be as efficacious as physical one in determining health concerns. Retina screening can be the very first clue to detecting a variety of hidden health issues including pre-diabetes and diabetes. Through the process of clinical diagnosis and prognosis; ophthalmologists rely heavily on the binary segmented version of retina fundus image; where the accuracy of segmented vessels, optic disc and abnormal lesions extremely affects the diagnosis accuracy which in turn affect the subsequent clinical treatment steps. This thesis proposes an automated retinal fundus image segmentation system composed of three segmentation subsystems follow same core segmentation algorithm. Despite of broad difference in features and characteristics; retinal vessels, optic disc and exudate lesions are extracted by each subsystem without the need for texture analysis or synthesis. For sake of compact diagnosis and complete clinical insight, our proposed system can detect these anatomical structures in one session with high accuracy even in pathological retina images. The proposed system uses a robust hybrid segmentation algorithm combines adaptive fuzzy thresholding and mathematical morphology. The proposed system is validated using four benchmark datasets: DRIVE and STARE (vessels), DRISHTI-GS (optic disc), and DIARETDB1 (exudates lesions). Competitive segmentation performance is achieved, outperforming a variety of up-to-date systems and demonstrating the capacity to deal with other heterogenous anatomical structures

Computational Analysis of Fundus Images: Rule-Based and Scale-Space Models

Fundus images are one of the most important imaging examinations in modern ophthalmology because they are simple, inexpensive and, above all, noninvasive. Nowadays, the acquisition and storage of highresolution fundus images is relatively easy and fast. Therefore, fundus imaging has become a fundamental investigation in retinal lesion detection, ocular health monitoring and screening programmes. Given the large volume and clinical complexity associated with these images, their analysis and interpretation by trained clinicians becomes a timeconsuming task and is prone to human error. Therefore, there is a growing interest in developing automated approaches that are affordable and have high sensitivity and specificity. These automated approaches need to be robust if they are to be used in the general population to diagnose and track retinal diseases. To be effective, the automated systems must be able to recognize normal structures and distinguish them from pathological clinical manifestations. The main objective of the research leading to this thesis was to develop automated systems capable of recognizing and segmenting retinal anatomical structures and retinal pathological clinical manifestations associated with the most common retinal diseases. In particular, these automated algorithms were developed on the premise of robustness and efficiency to deal with the difficulties and complexity inherent in these images. Four objectives were considered in the analysis of fundus images. Segmentation of exudates, localization of the optic disc, detection of the midline of blood vessels, segmentation of the vascular network and detection of microaneurysms. In addition, we also evaluated the detection of diabetic retinopathy on fundus images using the microaneurysm detection method. An overview of the state of the art is presented to compare the performance of the developed approaches with the main methods described in the literature for each of the previously described objectives. To facilitate the comparison of methods, the state of the art has been divided into rulebased methods and machine learningbased methods. In the research reported in this paper, rulebased methods based on image processing methods were preferred over machine learningbased methods. In particular, scalespace methods proved to be effective in achieving the set goals. Two different approaches to exudate segmentation were developed. The first approach is based on scalespace curvature in combination with the local maximum of a scalespace blob detector and dynamic thresholds. The second approach is based on the analysis of the distribution function of the maximum values of the noise map in combination with morphological operators and adaptive thresholds. Both approaches perform a correct segmentation of the exudates and cope well with the uneven illumination and contrast variations in the fundus images. Optic disc localization was achieved using a new technique called cumulative sum fields, which was combined with a vascular enhancement method. The algorithm proved to be reliable and efficient, especially for pathological images. The robustness of the method was tested on 8 datasets. The detection of the midline of the blood vessels was achieved using a modified corner detector in combination with binary philtres and dynamic thresholding. Segmentation of the vascular network was achieved using a new scalespace blood vessels enhancement method. The developed methods have proven effective in detecting the midline of blood vessels and segmenting vascular networks. The microaneurysm detection method relies on a scalespace microaneurysm detection and labelling system. A new approach based on the neighbourhood of the microaneurysms was used for labelling. Microaneurysm detection enabled the assessment of diabetic retinopathy detection. The microaneurysm detection method proved to be competitive with other methods, especially with highresolution images. Diabetic retinopathy detection with the developed microaneurysm detection method showed similar performance to other methods and human experts. The results of this work show that it is possible to develop reliable and robust scalespace methods that can detect various anatomical structures and pathological features of the retina. Furthermore, the results obtained in this work show that although recent research has focused on machine learning methods, scalespace methods can achieve very competitive results and typically have greater independence from image acquisition. The methods developed in this work may also be relevant for the future definition of new descriptors and features that can significantly improve the results of automated methods.As imagens do fundo do olho são hoje um dos principais exames imagiológicos da oftalmologia moderna, pela sua simplicidade, baixo custo e acima de tudo pelo seu carácter nãoinvasivo. A aquisição e armazenamento de imagens do fundo do olho com alta resolução é também relativamente simples e rápida. Desta forma, as imagens do fundo do olho são um exame fundamental na identificação de alterações retinianas, monitorização da saúde ocular, e em programas de rastreio. Considerando o elevado volume e complexidade clínica associada a estas imagens, a análise e interpretação das mesmas por clínicos treinados tornase uma tarefa morosa e propensa a erros humanos. Assim, há um interesse crescente no desenvolvimento de abordagens automatizadas, acessíveis em custo, e com uma alta sensibilidade e especificidade. Estas devem ser robustas para serem aplicadas à população em geral no diagnóstico e seguimento de doenças retinianas. Para serem eficazes, os sistemas de análise têm que conseguir detetar e distinguir estruturas normais de sinais patológicos. O objetivo principal da investigação que levou a esta tese de doutoramento é o desenvolvimento de sistemas automáticos capazes de detetar e segmentar as estruturas anatómicas da retina, e os sinais patológicos retinianos associados às doenças retinianas mais comuns. Em particular, estes algoritmos automatizados foram desenvolvidos segundo as premissas de robustez e eficácia para lidar com as dificuldades e complexidades inerentes a estas imagens. Foram considerados quatro objetivos de análise de imagens do fundo do olho. São estes, a segmentação de exsudados, a localização do disco ótico, a deteção da linha central venosa dos vasos sanguíneos e segmentação da rede vascular, e a deteção de microaneurismas. De acrescentar que usando o método de deteção de microaneurismas, avaliouse também a capacidade de deteção da retinopatia diabética em imagens do fundo do olho. Para comparar o desempenho das metodologias desenvolvidas neste trabalho, foi realizado um levantamento do estado da arte, onde foram considerados os métodos mais relevantes descritos na literatura para cada um dos objetivos descritos anteriormente. Para facilitar a comparação entre métodos, o estado da arte foi dividido em metodologias de processamento de imagem e baseadas em aprendizagem máquina. Optouse no trabalho de investigação desenvolvido pela utilização de metodologias de análise espacial de imagem em detrimento de metodologias baseadas em aprendizagem máquina. Em particular, as metodologias baseadas no espaço de escalas mostraram ser efetivas na obtenção dos objetivos estabelecidos. Para a segmentação de exsudados foram usadas duas abordagens distintas. A primeira abordagem baseiase na curvatura em espaço de escalas em conjunto com a resposta máxima local de um detetor de manchas em espaço de escalas e limiares dinâmicos. A segunda abordagem baseiase na análise do mapa de distribuição de ruído em conjunto com operadores morfológicos e limiares adaptativos. Ambas as abordagens fazem uma segmentação dos exsudados de elevada precisão, além de lidarem eficazmente com a iluminação nãouniforme e a variação de contraste presente nas imagens do fundo do olho. A localização do disco ótico foi conseguida com uma nova técnica designada por campos de soma acumulativos, combinada com métodos de melhoramento da rede vascular. O algoritmo revela ser fiável e eficiente, particularmente em imagens patológicas. A robustez do método foi verificada pela sua avaliação em oito bases de dados. A deteção da linha central dos vasos sanguíneos foi obtida através de um detetor de cantos modificado em conjunto com filtros binários e limiares dinâmicos. A segmentação da rede vascular foi conseguida com um novo método de melhoramento de vasos sanguíneos em espaço de escalas. Os métodos desenvolvidos mostraram ser eficazes na deteção da linha central dos vasos sanguíneos e na segmentação da rede vascular. Finalmente, o método para a deteção de microaneurismas assenta num formalismo de espaço de escalas na deteção e na rotulagem dos microaneurismas. Para a rotulagem foi utilizada uma nova abordagem da vizinhança dos candidatos a microaneurismas. A deteção de microaneurismas permitiu avaliar também a deteção da retinopatia diabética. O método para a deteção de microaneurismas mostrou ser competitivo quando comparado com outros métodos, em particular em imagens de alta resolução. A deteção da retinopatia diabética exibiu um desempenho semelhante a outros métodos e a especialistas humanos. Os trabalhos descritos nesta tese mostram ser possível desenvolver uma abordagem fiável e robusta em espaço de escalas capaz de detetar diferentes estruturas anatómicas e sinais patológicos da retina. Além disso, os resultados obtidos mostram que apesar de a pesquisa mais recente concentrarse em metodologias de aprendizagem máquina, as metodologias de análise espacial apresentam resultados muito competitivos e tipicamente independentes do equipamento de aquisição das imagens. As metodologias desenvolvidas nesta tese podem ser importantes na definição de novos descritores e características, que podem melhorar significativamente o resultado de métodos automatizados

Detection and characterisation of vessels in retinal images.

Doctor of Philosophy in Mathematics, Statistics & Computer Science. University of KwaZulu-Natal, Durban 2015.As retinopathies such as diabetic retinopathy (DR) and retinopathy of prematurity (ROP) continue to be the major causes of blindness globally, regular retinal examinations of patients can assist in the early detection of the retinopathies. The manual detection of retinal vessels is a very tedious and time consuming task as it requires about two hours to manually detect vessels in each retinal image. Automatic vessel segmentation has been helpful in achieving speed, improved diagnosis and progress monitoring of these diseases but has been challenging due to complexities such as the varying width of the retinal vessels from very large to very small, low contrast of thin vessels with respect to background and noise due to nonhomogeneous illumination in the retinal images. Although several supervised and unsupervised segmentation methods have been proposed in the literature, the segmentation of thinner vessels, connectivity loss of the vessels and time complexity remain the major challenges. In order to address these problems, this research work investigated di erent unsupervised segmentation approaches to be used in the robust detection of large and thin retinal vessels in a timely e cient manner. Firstly, this thesis conducted a study on the use of di erent global thresholding techniques combined with di erent pre-processing and post-processing techniques. Two histogram-based global thresholding techniques namely, Otsu and Isodata were able to detect large retinal vessels but fail to segment the thin vessels because these thin vessels have very low contrast and are di cult to distinguish from the background tissues using the histogram of the retinal images. Two new multi-scale approaches of computing global threshold based on inverse di erence moment and sum-entropy combined with phase congruence are investigated to improve the detection of vessels. One of the findings of this study is that the multi-scale approaches of computing global threshold combined with phase congruence based techniques improved on the detection of large vessels and some of the thin vessels. They, however, failed to maintain the width of the detected vessels. The reduction in the width of the detected large and thin vessels results in low sensitivity rates while relatively good accuracy rates were maintained. Another study on the use of fuzzy c-means and GLCM sum entropy combined on phase congruence for vessel segmentation showed that fuzzy c-means combined with phase congruence achieved a higher average accuracy rates of 0.9431 and 0.9346 but a longer running time of 27.1 seconds when compared with the multi-scale based sum entropy thresholding combined with phase congruence with the average accuracy rates of 0.9416 and 0.9318 with a running time of 10.3 seconds. The longer running time of the fuzzy c-means over the sum entropy thresholding is, however, attributed to the iterative nature of fuzzy c-means. When compared with the literature, both methods achieved considerable faster running time. This thesis investigated two novel local adaptive thresholding techniques for the segmentation of large and thin retinal vessels. The two novel local adaptive thresholding techniques applied two di erent Haralick texture features namely, local homogeneity and energy. Although these two texture features have been applied for supervised image segmentation in the literature, their novelty in this thesis lies in that they are applied using an unsupervised image segmentation approach. Each of these local adaptive thresholding techniques locally applies a multi-scale approach on each of the texture information considering the pixel of interest in relationship with its spacial neighbourhood to compute the local adaptive threshold. The localised multi-scale approach of computing the thresholds handled the challenge of the vessels' width variation. Experiments showed significant improvements in the average accuracy and average sensitivity rates of these techniques when compared with the previously discussed global thresholding methods and state of the art. The two novel local adaptive thresholding techniques achieved a higher reduction of false vessels around the border of the optic disc when compared with some of the previous techniques in the literature. These techniques also achieved a highly improved computational time of 1.9 to 3.9 seconds to segment the vessels in each retinal image when compared with the state of the art. Hence, these two novel local adaptive thresholding techniques are proposed for the segmentation of the vessels in the retinal images. This thesis further investigated the combination of di erence image and kmeans clustering technique for the segmentation of large and thin vessels in retinal images. The pre-processing phase computed a di erence image and k-means clustering technique was used for the vessel detection. While investigating this vessel segmentation method, this thesis established the need for a difference image that preserves the vessel details of the retinal image. Investigating the di erent low pass filters, median filter yielded the best di erence image required by k-means clustering for the segmentation of the retinal vessels. Experiments showed that the median filter based di erence images combined with k-means clustering technique achieved higher average accuracy and average sensitivity rates when compared with the previously discussed global thresholding methods and the state of the art. The median filter based di erence images combined with k-means clustering technique (that is, DIMDF) also achieved a higher reduction of false vessels around the border of the optic disc when compared with some previous techniques in the literature. These methods also achieved a highly improved computational time of 3.4 to 4 seconds when compared with the literature. Hence, the median filter based di erence images combined with k-means clustering technique are proposed for the segmentation of the vessels in retinal images. The characterisation of the detected vessels using tortuosity measure was also investigated in this research. Although several vessel tortuosity methods have been discussed in the literature, there is still need for an improved method that e ciently detects vessel tortuosity. The experimental study conducted in this research showed that the detection of the stationary points helps in detecting the change of direction and twists in the vessels. The combination of the vessel twist frequency obtained using the stationary points and distance metric for the computation of normalised and nonnormalised tortuosity index (TI) measure was investigated. Experimental results showed that the non-normalised TI measure had a stronger correlation with the expert's ground truth when compared with the distance metric and normalised TI measures. Hence, a non-normalised TI measure that combines the vessel twist frequency based on the stationary points and distance metric is proposed for the measurement of vessel tortuosity