Optic nerve head three-dimensional shape analysis

We present a method for optic nerve head (ONH) 3-D shape analysis from retinal optical coherence tomography (OCT). The possibility to noninvasively acquire in vivo high-resolution 3-D volumes of the ONH using spectral domain OCT drives the need to develop tools that quantify the shape of this structure and extract information for clinical applications. The presented method automatically generates a 3-D ONH model and then allows the computation of several 3-D parameters describing the ONH. The method starts with a high-resolution OCT volume scan as input. From this scan, the model-defining inner limiting membrane (ILM) as inner surface and the retinal pigment epithelium as outer surface are segmented, and the Bruch's membrane opening (BMO) as the model origin is detected. Based on the generated ONH model by triangulated 3-D surface reconstruction, different parameters (areas, volumes, annular surface ring, minimum distances) of different ONH regions can then be computed. Additionally, the bending energy (roughness) in the BMO region on the ILM surface and 3-D BMO-MRW surface area are computed. We show that our method is reliable and robust across a large variety of ONH topologies (specific to this structure) and present a first clinical application

The morphogenesis of the zebrafish eye, including a fate map of the optic vesicle

We have examined the morphogenesis of the zebrafish eye, from the flat optic vesicle at 16 hours post fertilization (hpf) to the functional hemispheric eye at 72 hpf. We have produced three-dimensional reconstructions from semithin sections, measured volumes and areas, and produced a fate map by labeling clusters of cells at 14–15 hpf and finding them in the 24 hpf eye cup. Both volume and area increased sevenfold, with different schedules. Initially (16–33 hpf), area increased but volume remained constant; later (33–72 hpf) both increased. When the volume remained constant, the presumptive pigmented epithelium (PE) shrank and the presumptive neural retina (NR) enlarged. The fate map revealed that during 14–24 hpf cells changed layers, moving from the PE into the NR, probably through involution around the margin of the eye. The transformation of the flat epithelial layers of the vesicle into their cup-shaped counterparts in the eye was also accompanied by cellular rearrangements; most cells in a cluster labeled in the vesicle remained neighbors in the eye cup, but occasionally they were separated widely. This description of normal zebrafish eye development provides explanations for some mutant phenotypes and for the effects of altered retinoic acid. Dev Dyn;218:175–188. © 2000 Wiley-Liss, Inc.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/35166/1/15_ftp.pd

Surface Denoising based on The Variation of Normals and Retinal Shape Analysis

Through the development of this thesis, starting from the curvature tensor, we have been able to understand the variation of tangent vectors to define a shape analysis operator and also a relationship between the classical shape operator and the curvature tensor on a triangular surface. In continuation, the first part of the thesis analyzed the variation of surface normals and introduced a shape analysis operator, which is further used for mesh and point set denoising. In the second part of the thesis, mathematical modeling and shape quantification algorithms are introduced for retinal shape analysis. In the first half, this thesis followed the concept of the variation of surface normals, which is termed as the normal voting tensor and derived a relation between the shape operator and the normal voting tensor. The concept of the directional and the mean curvatures is extended on the dual representation of a triangulated surface. A normal voting tensor is defined on each triangle of a geometry and termed as the element-based normal voting tensor (ENVT). Later, a deformation tensor is extracted from the ENVT and it consists of the anisotropy of a surface and the mean curvature vector is defined based on the ENVT deformation tensor. The ENVT-based mesh denoising algorithm is introduced, where the ENVT is used as a shape operator. A binary optimization technique is applied on the spectral components of the ENVT that helps the algorithm to retain sharp features in the concerned geometry and improves the convergence rate of the algorithm. Later, a stochastic analysis of the effect of noise on the triangular mesh based on the minimum edge length of the elements in the geometry is explained. It gives an upper bound to the noise standard deviation to have minimum probability for flipped element normals. The ENVT-based mesh denoising concept is extended for a point set denoising, where noisy vertex normals are filtered using the vertex-based NVT and the binary optimization. For vertex update stage in point set denoising, we added different constraints to the quadratic error metric based on features (edges and corners) or non-feature (planar) points. This thesis also investigated a robust statistics framework for face normal bilateral filtering and proposed a robust and high fidelity two-stage mesh denoising method using Tukey’s bi-weight function as a robust estimator, which stops the diffusion at sharp features and produces smooth umbilical regions. This algorithm introduced a novel vertex update scheme, which uses a differential coordinate-based Laplace operator along with an edge-face normal orthogonality constraint to produce a high-quality mesh without face normal flips and it also makes the algorithm more robust against high-intensity noise. The second half of thesis focused on the application of the proposed geometric processing algorithms on the OCT (optical coherence tomography) scan data for quantification of the human retinal shape. The retina is a part of the central nervous system and comprises a similar cellular composition as the brain. Therefore, many neurological disorders affect the retinal shape and these neuroinflammatory conditions are known to cause modifications to two important regions of the retina: the fovea and the optical nerve head (ONH). This thesis consists of an accurate and robust shape modeling of these regions to diagnose several neurological disorders by detecting the shape changes. For the fovea, a parametric modeling algorithm is introduced using Cubic Bezier and this algorithm derives several 3D shape parameters, which quantify the foveal shape with high accuracy. For the ONH, a 3D shape analysis algorithm is introduced to measure the shape variation regarding different neurological disorders. The proposed algorithm uses triangulated manifold surfaces of two different layers of the retina to derive several 3D shape parameters. The experimental results of the fovea and the ONH morphometry confirmed that these algorithms can provide an aid to diagnose several neurological disorders

Deep learning network to correct axial and coronal eye motion in 3D OCT retinal imaging

Optical Coherence Tomography (OCT) is one of the most important retinal imaging technique. However, involuntary motion artifacts still pose a major challenge in OCT imaging that compromises the quality of downstream analysis, such as retinal layer segmentation and OCT Angiography. We propose deep learning based neural networks to correct axial and coronal motion artifacts in OCT based on a single volumetric scan. The proposed method consists of two fully-convolutional neural networks that predict Z and X dimensional displacement maps sequentially in two stages. The experimental result shows that the proposed method can effectively correct motion artifacts and achieve smaller error than other methods. Specifically, the method can recover the overall curvature of the retina, and can be generalized well to various diseases and resolutions

Computerised stereoscopic measurement of the human retina

The research described herein is an investigation into the problems of obtaining useful clinical measurements from stereo photographs of the human retina through automation of the stereometric procedure by digital stereo matching and image analysis techniques. Clinical research has indicated a correlation between physical changes to the optic disc topography (the region on the retina where the optic nerve enters the eye) and the advance of eye disease such as hypertension and glaucoma. Stereoscopic photography of the human retina (or fundus, as it is called) and the subsequent measurement of the topography of the optic disc is of great potential clinical value as an aid in observing the pathogenesis of such disease, and to this end, accurate measurements of the various parameters that characterise the changing shape of the optic disc topography must be provided. Following a survey of current clinical methods for stereoscopic measurement of the optic disc, fundus image data acquisition, stereo geometry, limitations of resolution and accuracy, and other relevant physical constraints related to fundus imaging are investigated. A survey of digital stereo matching algorithms is presented and their strengths and weaknesses are explored, specifically as they relate to the suitability of the algorithm for the fundus image data. The selection of an appropriate stereo matching algorithm is discussed, and its application to four test data sets is presented in detail. A mathematical model of two-dimensional image formation is developed together with its corresponding auto-correlation function. In the presense of additive noise, the model is used as a tool for exploring key problems with respect to the stereo matching of fundus images. Specifically, measures for predicting correlation matching error are developed and applied. Such measures are shown to be of use in applications where the results of image correlation cannot be independently verified, and meaningful quantitative error measures are required. The application of these theoretical tools to the fundus image data indicate a systematic way to measure, assess and control cross-correlation error. Conclusions drawn from this research point the way forward for stereo analysis of the optic disc and highlight a number of areas which will require further research. The development of a fully automated system for diagnostic evaluation of the optic disc topography is discussed in the light of the results obtained during this research

Neue Methoden der Nachbearbeitung und Analyse retinaler optischer Kohärenztomografieaufnahmen bei neurologischen Erkrankungen

Viele neurologische Krankheiten verursachen Veränderungen in der Netzhaut, die mit Hilfe der optischen Kohärenztomography (optical coherence tomography, OCT) dargestellt werden können. Dabei entstehen viele Bilddaten, deren Auswertung zeitintensiv ist und geschultes Personal erfordert. Ziel dieser Arbeit war die Entwicklung neuer Methoden zur Vorverarbeitung und Analyse retinaler OCT-Bilddaten, um Outcome-Parameter für Studien und diagnostische Marker für neurologische Erkrankungen zu verbessern. Dazu wurden Methoden für zwei wichtige Aufnahmebereiche der Netzhaut, den Sehnervenkopf (optic nerve head, ONH) und die Makula, entwickelt. Für den ONH-Bereich wurde eine automatische Segmentierung auf Basis aktiver Konturen entwickelt, die eine akkurate Segmentierung der inneren Grenzmembran auch bei komplexer Topografie ermöglicht. Für den Bereich um die Makula entstand eine intraretinale Schichtensegmentierungspipeline, die von der Auswahl der Bilddaten über die automatische Segmentierung sowie die manuelle Nachkorrektur bis zur Ausgabe verschiedener Schichtdicken in Tabellenform reicht. Für beide Aufnahmebereiche wurden mehrere Programme entwickelt, die auf einer gemeinsamen Basis zur Verarbeitung von OCT-Daten fußen. Eines dieser Programme bietet eine grafische Oberfläche zur manuellen Verarbeitung der Bilddaten. Mit dieser Software wurden Teile der Referenzdaten manuell erstellt, die innere Grenzmembran des ONH automatisch segmentiert sowie eine komfortable Nachbearbeitung von intraretinalen Segmentierungen vorgenommen. Dies ermöglichte die automatische Auswertung morphologischer Parameter des ONH, wovon einige signifikante Unterschiede zwischen Patienten mit neurologischen Krankheiten und gesunden Kontrollen zeigten. Weiter kam die Schichtensegmentierungspipeline beim Aufbau einer normativen Datenbank sowie in einer Studie zum Zusammenhang des retinalen Schadens mit der kritischen Flimmerfrequenz zum Einsatz. Ein Teil der Software wurde als freie und quelloffene Software (free and open-source software, FOSS) und der normative Datensatz für die Verwendung in anderen Studien freigegeben. Beides wird bereits in weiteren Studien eingesetzt und wird auch die Durchführung zukünftiger Studien vereinfachen sowie die Entwicklung neuer Methoden unterstützen.Many neurological diseases cause changes in the retina, which can be visualized using optical coherence tomography (OCT). This process produces large amounts of image data. Its evaluation is time-consuming and requires medically trained personnel. This dissertation aims to develop new methods for preprocessing and analyzing retinal OCT data in order to improve outcome parameters for clinical studies and diagnostic markers for neurological diseases. For this purpose, methods concerning the regions of two landmarks of the retina, the optic nerve head (ONH) and the macula, were developed. For the ONH, an automatic segmentation method based on active contours was developed, which allows accurate segmentation of the inner limiting membrane even in complex topography. For the macular region, an intraretinal layer segmentation pipeline from image data via automatic segmentation to manual post-correction and the output of different layer thicknesses in tabular form was developed. For both, ONH and macular region, several programs were developed, which share a common basis for processing OCT data. One of these programs offers a graphical user interface for the manual processing of image data. Parts of the reference data were created manually using this software. Moreover, the inner limiting membrane of the ONH was segmented automatically and post-processing of intraretinal segmentations was performed. This allowed for automatic evaluation of morphological parameters of the ONH, some of which showed significant differences between patients with neurological diseases and the healthy control group. Furthermore, the layer segmentation pipeline was utilized to create a normative database as well as to investigate the correlation of retinal damage and critical flicker frequency. Part of the software was released as free and open-source software (FOSS) and the normative data set was released for use in other studies. Both are already being used in further studies and will also aid in future studies, as well as support the development of new methods

Applications of three-dimensional printing in ophthalmology

Three-dimensional (3D) printing is increasingly used to produce customised objects and is a promising alternative to traditional manufacturing methods in diverse fields, such as dentistry and orthopaedics. Already in use in other medical specialities, adoption in ophthalmology has been limited to date. This review aims to provide an overview of 3D printing technology with respect to current and potential applications in ophthalmic practice. Medline, Embase and internet search were performed with "3D printing", "ophthalmology", "dentistry", "orthopaedics" and their synonyms used as main search terms. In addition, search terms related to clinical applications such as "surgery" and "implant" were employed. 3D printing has multiple applications in ophthalmology, including in diagnosis, surgery, prosthetics, medications and medical education. Within the past decade, researchers have produced 3D printed models of objects such as implants, prostheses, anatomical models and surgical simulators. Further development is necessary to generate optimal biomaterials for various applications, and the quality and long-term performance of 3D models needs to be validated

Review on retrospective procedures to correct retinal motion artefacts in OCT imaging

Motion artefacts from involuntary changes in eye fixation remain a major imaging issue in optical coherence tomography (OCT). This paper reviews the state-of-the-art of retrospective procedures to correct retinal motion and axial eye motion artefacts in OCT imaging. Following an overview of motion induced artefacts and correction strategies, a chronological survey of retrospective approaches since the introduction of OCT until the current days is presented. Pre-processing, registration, and validation techniques are described. The review finishes by discussing the limitations of the current techniques and the challenges to be tackled in future developments

Improving the forward model for electrical impedance tomography of brain function through rapid generation of subject specific finite element models

Electrical Impedance Tomography (EIT) is a non-invasive imaging method which allows internal electrical impedance of any conductive object to be imaged by means of current injection and surface voltage measurements through an array of externally applied electrodes. The successful generation of the image requires the simulation of the current injection patterns on either an analytical or a numerical model of the domain under examination, known as the forward model, and using the resulting voltage data in the inverse solution from which images of conductivity changes can be constructed. Recent research strongly indicates that geometric and anatomical conformance of the forward model to the subject under investigation significantly affects the quality of the images. This thesis focuses mainly on EIT of brain function and describes a novel approach for the rapid generation of patient or subject specific finite element models for use as the forward model. After introduction of the topic, methods of generating accurate finite element (FE) models using commercially available Computer-Aided Design (CAD) tools are described and show that such methods, though effective and successful, are inappropriate for time critical clinical use. The feasibility of warping or morphing a finite element mesh as a means of reducing the lead time for model generation is then presented and demonstrated. This leads on to the description of methods of acquiring and utilising known system geometry, namely the positions of electrodes and registration landmarks, to construct an accurate surface of the subject, the results of which are successfully validated. The outcome of this procedure is then used to specify boundary conditions to a mesh warping algorithm based on elastic deformation using well-established continuum mechanics procedures. The algorithm is applied to a range of source models to empirically establish optimum values for the parameters defining the problem which can successfully generate meshes of acceptable quality in terms of discretization errors and which more accurately define the geometry of the target subject. Further validation of the algorithm is performed by comparison of boundary voltages and image reconstructions from simulated and laboratory data to demonstrate that benefits in terms of image artefact reduction and localisation of conductivity changes can be gained. The processes described in the thesis are evaluated and discussed and topics of further work and application are described