This study aims to analyze the asymmetry between both eyes of the same patient for the early diagnosis of glaucoma. Two imaging modalities, retinal fundus images and optical coherence tomographies ...(OCTs), have been considered in order to compare their different capabilities for glaucoma detection. From retinal fundus images, the difference between cup/disc ratio and the width of the optic rim has been extracted. Analogously, the thickness of the retinal nerve fiber layer has been measured in spectral-domain optical coherence tomographies. These measurements have been considered as asymmetry characteristics between eyes in the modeling of decision trees and support vector machines for the classification of healthy and glaucoma patients. The main contribution of this work is indeed the use of different classification models with both imaging modalities to jointly exploit the strengths of each of these modalities for the same diagnostic purpose based on the asymmetry characteristics between the eyes of the patient. The results show that the optimized classification models provide better performance with OCT asymmetry features between both eyes (sensitivity 80.9%, specificity 88.2%, precision 66.7%, accuracy 86.5%) than with those extracted from retinographies, although a linear relationship has been found between certain asymmetry features extracted from both imaging modalities. Therefore, the resulting performance of the models based on asymmetry features proves their ability to differentiate healthy from glaucoma patients using those metrics. Models trained from fundus characteristics are a useful option as a glaucoma screening method in the healthy population, although with lower performance than those trained from the thickness of the peripapillary retinal nerve fiber layer. In both imaging modalities, the asymmetry of morphological characteristics can be used as a glaucoma indicator, as detailed in this work.
Purpose: The aim of this study was to analyze the relevance of asymmetry features between both eyes of the same patient for glaucoma screening using optical coherence tomography. Methods: ...Spectral-domain optical coherence tomography was used to estimate the thickness of the peripapillary retinal nerve fiber layer in both eyes of the patients in the study. These measurements were collected in a dataset from healthy and glaucoma patients. Several metrics for asymmetry in the retinal nerve fiber layer thickness between the two eyes were then proposed. These metrics were evaluated using the dataset by performing a statistical analysis to assess their significance as relevant features in the diagnosis of glaucoma. Finally, the usefulness of these asymmetry features was demonstrated by designing supervised machine learning models that can be used for the early diagnosis of glaucoma. Results: Machine learning models were designed and optimized, specifically decision trees, based on the values of proposed asymmetry metrics. The use of these models on the dataset provided good classification of the patients (accuracy 88%, sensitivity 70%, specificity 93% and precision 75%). Conclusions: The obtained machine learning models based on retinal nerve fiber layer asymmetry are simple but effective methods which offer a good trade-off in classification of patients and simplicity. The fast binary classification relies on a few asymmetry values of the retinal nerve fiber layer thickness, allowing their use in the daily clinical practice for glaucoma screening.
Glaucoma is a neurodegenerative disease process that leads to progressive damage of the optic nerve to produce visual impairment and blindness. Spectral-domain OCT technology enables peripapillary ...circular scans of the retina and the measurement of the thickness of the retinal nerve fiber layer (RNFL) for the assessment of the disease status or progression in glaucoma patients. This paper describes a new approach to segment and measure the retinal nerve fiber layer in peripapillary OCT images. The proposed method consists of two stages. In the first one, morphological operators robustly detect the coarse location of the layer boundaries, despite the speckle noise and diverse artifacts in the OCT image. In the second stage, deformable models are initialized with the results of the previous stage to perform a fine segmentation of the boundaries, providing an accurate measurement of the entire RNFL. The results of the RNFL segmentation were qualitatively assessed by ophthalmologists, and the measurements of the thickness of the RNFL were quantitatively compared with those provided by the OCT inbuilt software as well as the state-of-the-art methods.
This article presents a hybrid framework for efficient and accurate orientation estimation. The proposed scheme combines the single orientation information given by a novel method and the multiple ...orientation information provided by a bank of linear orientated morphological openings. The single orientations are estimated by means of an energy-minimization Gaussian filtering which solves the drawback related to phase changes of other methods. After describing the formulation of these two approaches for estimating the existing orientations in the pixels of an image, several strategies have been analyzed to fuse and discriminate the information of both orientation vector fields in the resulting hybrid orientation vector field. The objective of the proposed hybrid method is to reduce the computational cost involved in calculating multiple orientations only in those pixels where they exist while maintaining the accuracy provided by the single orientation method in the remaining pixels. To this end, strategies ranging from a threshold in the multiple orientation vector field to a convolutional neural network trained with a set of patterns specifically designed to detect pixels with multiple orientations, passing through the Harris corner detector, have been tested to identify those pixels where multiple orientations exist. Results on natural and synthetic images show the accuracy and the computational efficiency achieved by the proposed hybrid framework to provide the vector field with single and multiple orientations.
Glaucoma is one of the leading causes of blindness worldwide and Optical Coherence Tomography (OCT) is the quintessential imaging technique for its detection. Unlike most of the state-of-the-art ...studies focused on glaucoma detection, in this paper, we propose, for the first time, a novel framework for glaucoma grading using raw circumpapillary B-scans. In particular, we set out a new OCT-based hybrid network which combines hand-driven and deep learning algorithms. An OCT-specific descriptor is proposed to extract hand-crafted features related to the retinal nerve fibre layer (RNFL). In parallel, an innovative CNN is developed using skip-connections to include tailored residual and attention modules to refine the automatic features of the latent space. The proposed architecture is used as a backbone to conduct a novel few-shot learning based on static and dynamic prototypical networks. The k-shot paradigm is redefined giving rise to a supervised end-to-end system which provides substantial improvements discriminating between healthy, early and advanced glaucoma samples. The training and evaluation processes of the dynamic prototypical network are addressed from two fused databases acquired via Heidelberg Spectralis system. Validation and testing results reach a categorical accuracy of 0.9459 and 0.8788 for glaucoma grading, respectively. Besides, the high performance reported by the proposed model for glaucoma detection deserves a special mention. The findings from the class activation maps are directly in line with the clinicians' opinion since the heatmaps pointed out the RNFL as the most relevant structure for glaucoma diagnosis.
•Circumpapillary OCT images in the spectral domain are used for the first time to grade the glaucoma severity.•The few-shot methodology is redefined in the prototypical paradigm by proposing an optimized k-shot supervised learning which allows exploiting the labelled information.•Different prototypical-based solutions are conducted and compared with conventional approaches for glaucoma grading.•The proposed fully supervised prototypical neural network outperforms the previous results achieved in the state-of-the-art for glaucoma detection and glaucoma grading.•The heat maps extracted are directly in line with the clinician’s opinion, since the highlighted regions of the B-scans correspond to the interesting areas in which the experts focus for glaucoma diagnosis.
Texture feature extraction is an important task in image processing and computer vision. Typical applications include automated inspection, image retrieval or medical analysis. In this paper we ...propose a noise robust and rotation invariant approach to texture feature analysis and classification. The proposed framework is based on a simple texture feature extraction step which decomposes the image structures by means of a family of orientated filters. The output is an orientation density vector which is called Orientation Feature Vector (OFV). The OFV can be used as the input feature vector to a texture classifier, and in addition, by using an interpolation step it is possible to extract the main orientations of the texture from the OFD, providing additional high level features for image analysis. In this work, three families of filters have been studied in the texture feature extraction step. The experimental results show the ability of the proposed framework in classification problems, improving more than 20% the results of other state-of-the-art methods when a high level of Gaussian noise is considered.
In the pioneering work, the radiation diagrams of leaky wave antenna arrays can achieve attenuation nulls and gains at specific angles by manually placing the zeros and the poles in the Z domain of ...the corresponding discrete linear time‐invariant (LTI) system. This handcrafted design procedure does not allow radiation diagrams with wide beams since the interaction between poles involved in the wide beam and their corresponding leaky modes cannot be easily handled. To overcome this limitation, this paper describes a novel method for designing radiation diagrams of leaky‐wave antenna arrays based on the theory of IIR discrete filters. The proposed method relies on the design of discrete filters with the prototypes of analog low‐pass filters defined by Butterworth and Chebyshev type I polynomials, whose roots along with the bilinear transformation provide the location of the poles and the zeros of the discrete LTI system and, therefore, the parameters of the leaky‐wave antenna array. Results with different designs and a comparison with other approaches show the utility and effectiveness of this novel method to design wide‐beam leaky‐wave antenna arrays.
Approach based on the relationship between the zeros and poles of a discrete linear time‐invariant (LTI) low‐pass filter designed using Butterworth and Chebyshev type I polynomials and the parameters of a leaky‐wave antenna (LWA) array. The method is able to produce a pass‐band without ripple, joint control on the pass‐band and rejection‐band, as well as a high attenuation level in the rejection‐band.
Glaucoma is one of the ophthalmological diseases that frequently causes loss of vision in today's society. Previous studies assess which anatomical parameters of the optic nerve can be predictive of ...glaucomatous damage, but to date there is no test that by itself has sufficient sensitivity and specificity to diagnose this disease. This work provides a public dataset with medical data and fundus images of both eyes of the same patient. Segmentations of the cup and optic disc, as well as the labeling of the patients based on the evaluation of clinical data are also provided. The dataset has been tested with a neural network to classify healthy and glaucoma patients. Specifically, the ResNet-50 has been used as the basis to classify patients using information from each eye independently as well as using the joint information from both eyes of each patient. Results provide the baseline metrics, with the aim of promoting research in the early detection of glaucoma based on the joint analysis of both eyes of the same patient.
This paper deals with the theory and applications of spatially-variant discrete mathematical morphology. We review and formalize the definition of spatially variant dilation/erosion and ...opening/closing for binary and gray-level images using exclusively the structuring function, without resorting to complement. This theoretical framework allows to build morphological operators whose structuring elements can locally adapt their shape and orientation across the dominant direction of the structures in the image. The shape and orientation of the structuring element at each pixel are extracted from the image under study: the orientation is given by means of a diffusion process of the average square gradient field, which regularizes and extends the orientation information from the edges of the objects to the homogeneous areas of the image; and the shape of the orientated structuring elements can be linear or it can be given by the distance to relevant edges of the objects. The proposed filters are used on binary and gray-level images for enhancement of anisotropic features such as coherent, flow-like structures. Results of spatially-variant erosions/dilations and openings/closings-based filters prove the validity of this theoretical sound and novel approach.
This paper addresses the formulation of adaptive morphological filters based on spatially-variant structuring elements. The adaptivity of these filters is achieved by modifying the shape and ...orientation of the structuring elements according to a multiple orientation vector field. This vector field is provided by means of a bank of directional openings which can take into account the possible multiple orientations of the contours in the image. After reviewing and formalizing the definition of the spatially-variant dilation, erosion, opening and closing, the proposed structuring elements are described. These spatially-variant structuring elements are based on ellipses which vary over the image domain adapting locally their orientation according to the multiple orientation vector field and their shape (the eccentricity of the ellipses) according to the distance to relevant contours of the objects. The proposed adaptive morphological filters are used on gray-level images and are compared with spatially-invariant filters, with spatially-variant filters based on a single orientation vector field, and with adaptive morphological bilateral filters. Results show that the morphological filters based on a multiple orientation vector field are more adept at enhancing and preserving structures which contains more than one orientation.