Automated detection and classification of brain injury lesions on structural magnetic resonance imaging

El daño cerebral adquirido (DCA) es un problema social y sanitario grave, de magnitud creciente y de una gran complejidad diagnóstica y terapéutica. Su elevada incidencia, junto con el aumento de la supervivencia de los pacientes, una vez superada la fase aguda, lo convierten también en un problema de alta prevalencia. En concreto, según la Organización Mundial de la Salud (OMS) el DCA estará entre las 10 causas más comunes de discapacidad en el año 2020. La neurorrehabilitación permite mejorar el déficit tanto cognitivo como funcional y aumentar la autonomía de las personas con DCA. Con la incorporación de nuevas soluciones tecnológicas al proceso de neurorrehabilitación se pretende alcanzar un nuevo paradigma donde se puedan diseñar tratamientos que sean intensivos, personalizados, monitorizados y basados en la evidencia. Ya que son estas cuatro características las que aseguran que los tratamientos son eficaces. A diferencia de la mayor parte de las disciplinas médicas, no existen asociaciones de síntomas y signos de la alteración cognitiva que faciliten la orientación terapéutica. Actualmente, los tratamientos de neurorrehabilitación se diseñan en base a los resultados obtenidos en una batería de evaluación neuropsicológica que evalúa el nivel de afectación de cada una de las funciones cognitivas (memoria, atención, funciones ejecutivas, etc.). La línea de investigación en la que se enmarca este trabajo de investigación pretende diseñar y desarrollar un perfil cognitivo basado no sólo en el resultado obtenido en esa batería de test, sino también en información teórica que engloba tanto estructuras anatómicas como relaciones funcionales e información anatómica obtenida de los estudios de imagen. De esta forma, el perfil cognitivo utilizado para diseñar los tratamientos integra información personalizada y basada en la evidencia. Las técnicas de neuroimagen representan una herramienta fundamental en la identificación de lesiones para la generación de estos perfiles cognitivos. La aproximación clásica utilizada en la identificación de lesiones consiste en delinear manualmente regiones anatómicas cerebrales. Esta aproximación presenta diversos problemas relacionados con inconsistencias de criterio entre distintos clínicos, reproducibilidad y tiempo. Por tanto, la automatización de este procedimiento es fundamental para asegurar una extracción objetiva de información. La delineación automática de regiones anatómicas se realiza mediante el registro tanto contra atlas como contra otros estudios de imagen de distintos sujetos. Sin embargo, los cambios patológicos asociados al DCA están siempre asociados a anormalidades de intensidad y/o cambios en la localización de las estructuras. Este hecho provoca que los algoritmos de registro tradicionales basados en intensidad no funcionen correctamente y requieran la intervención del clínico para seleccionar ciertos puntos (que en esta tesis hemos denominado puntos singulares). Además estos algoritmos tampoco permiten que se produzcan deformaciones grandes deslocalizadas. Hecho que también puede ocurrir ante la presencia de lesiones provocadas por un accidente cerebrovascular (ACV) o un traumatismo craneoencefálico (TCE). Esta tesis se centra en el diseño, desarrollo e implementación de una metodología para la detección automática de estructuras lesionadas que integra algoritmos cuyo objetivo principal es generar resultados que puedan ser reproducibles y objetivos. Esta metodología se divide en cuatro etapas: pre-procesado, identificación de puntos singulares, registro y detección de lesiones. Los trabajos y resultados alcanzados en esta tesis son los siguientes: Pre-procesado. En esta primera etapa el objetivo es homogeneizar todos los datos de entrada con el objetivo de poder extraer conclusiones válidas de los resultados obtenidos. Esta etapa, por tanto, tiene un gran impacto en los resultados finales. Se compone de tres operaciones: eliminación del cráneo, normalización en intensidad y normalización espacial. Identificación de puntos singulares. El objetivo de esta etapa es automatizar la identificación de puntos anatómicos (puntos singulares). Esta etapa equivale a la identificación manual de puntos anatómicos por parte del clínico, permitiendo: identificar un mayor número de puntos lo que se traduce en mayor información; eliminar el factor asociado a la variabilidad inter-sujeto, por tanto, los resultados son reproducibles y objetivos; y elimina el tiempo invertido en el marcado manual de puntos. Este trabajo de investigación propone un algoritmo de identificación de puntos singulares (descriptor) basado en una solución multi-detector y que contiene información multi-paramétrica: espacial y asociada a la intensidad. Este algoritmo ha sido contrastado con otros algoritmos similares encontrados en el estado del arte. Registro. En esta etapa se pretenden poner en concordancia espacial dos estudios de imagen de sujetos/pacientes distintos. El algoritmo propuesto en este trabajo de investigación está basado en descriptores y su principal objetivo es el cálculo de un campo vectorial que permita introducir deformaciones deslocalizadas en la imagen (en distintas regiones de la imagen) y tan grandes como indique el vector de deformación asociado. El algoritmo propuesto ha sido comparado con otros algoritmos de registro utilizados en aplicaciones de neuroimagen que se utilizan con estudios de sujetos control. Los resultados obtenidos son prometedores y representan un nuevo contexto para la identificación automática de estructuras. Identificación de lesiones. En esta última etapa se identifican aquellas estructuras cuyas características asociadas a la localización espacial y al área o volumen han sido modificadas con respecto a una situación de normalidad. Para ello se realiza un estudio estadístico del atlas que se vaya a utilizar y se establecen los parámetros estadísticos de normalidad asociados a la localización y al área. En función de las estructuras delineadas en el atlas, se podrán identificar más o menos estructuras anatómicas, siendo nuestra metodología independiente del atlas seleccionado. En general, esta tesis doctoral corrobora las hipótesis de investigación postuladas relativas a la identificación automática de lesiones utilizando estudios de imagen médica estructural, concretamente estudios de resonancia magnética. Basándose en estos cimientos, se han abrir nuevos campos de investigación que contribuyan a la mejora en la detección de lesiones. ABSTRACT Brain injury constitutes a serious social and health problem of increasing magnitude and of great diagnostic and therapeutic complexity. Its high incidence and survival rate, after the initial critical phases, makes it a prevalent problem that needs to be addressed. In particular, according to the World Health Organization (WHO), brain injury will be among the 10 most common causes of disability by 2020. Neurorehabilitation improves both cognitive and functional deficits and increases the autonomy of brain injury patients. The incorporation of new technologies to the neurorehabilitation tries to reach a new paradigm focused on designing intensive, personalized, monitored and evidence-based treatments. Since these four characteristics ensure the effectivity of treatments. Contrary to most medical disciplines, it is not possible to link symptoms and cognitive disorder syndromes, to assist the therapist. Currently, neurorehabilitation treatments are planned considering the results obtained from a neuropsychological assessment battery, which evaluates the functional impairment of each cognitive function (memory, attention, executive functions, etc.). The research line, on which this PhD falls under, aims to design and develop a cognitive profile based not only on the results obtained in the assessment battery, but also on theoretical information that includes both anatomical structures and functional relationships and anatomical information obtained from medical imaging studies, such as magnetic resonance. Therefore, the cognitive profile used to design these treatments integrates information personalized and evidence-based. Neuroimaging techniques represent an essential tool to identify lesions and generate this type of cognitive dysfunctional profiles. Manual delineation of brain anatomical regions is the classical approach to identify brain anatomical regions. Manual approaches present several problems related to inconsistencies across different clinicians, time and repeatability. Automated delineation is done by registering brains to one another or to a template. However, when imaging studies contain lesions, there are several intensity abnormalities and location alterations that reduce the performance of most of the registration algorithms based on intensity parameters. Thus, specialists may have to manually interact with imaging studies to select landmarks (called singular points in this PhD) or identify regions of interest. These two solutions have the same inconvenient than manual approaches, mentioned before. Moreover, these registration algorithms do not allow large and distributed deformations. This type of deformations may also appear when a stroke or a traumatic brain injury (TBI) occur. This PhD is focused on the design, development and implementation of a new methodology to automatically identify lesions in anatomical structures. This methodology integrates algorithms whose main objective is to generate objective and reproducible results. It is divided into four stages: pre-processing, singular points identification, registration and lesion detection. Pre-processing stage. In this first stage, the aim is to standardize all input data in order to be able to draw valid conclusions from the results. Therefore, this stage has a direct impact on the final results. It consists of three steps: skull-stripping, spatial and intensity normalization. Singular points identification. This stage aims to automatize the identification of anatomical points (singular points). It involves the manual identification of anatomical points by the clinician. This automatic identification allows to identify a greater number of points which results in more information; to remove the factor associated to inter-subject variability and thus, the results are reproducible and objective; and to eliminate the time spent on manual marking. This PhD proposed an algorithm to automatically identify singular points (descriptor) based on a multi-detector approach. This algorithm contains multi-parametric (spatial and intensity) information. This algorithm has been compared with other similar algorithms found on the state of the art. Registration. The goal of this stage is to put in spatial correspondence two imaging studies of different subjects/patients. The algorithm proposed in this PhD is based on descriptors. Its main objective is to compute a vector field to introduce distributed deformations (changes in different imaging regions), as large as the deformation vector indicates. The proposed algorithm has been compared with other registration algorithms used on different neuroimaging applications which are used with control subjects. The obtained results are promising and they represent a new context for the automatic identification of anatomical structures. Lesion identification. This final stage aims to identify those anatomical structures whose characteristics associated to spatial location and area or volume has been modified with respect to a normal state. A statistical study of the atlas to be used is performed to establish which are the statistical parameters associated to the normal state. The anatomical structures that may be identified depend on the selected anatomical structures identified on the atlas. The proposed methodology is independent from the selected atlas. Overall, this PhD corroborates the investigated research hypotheses regarding the automatic identification of lesions based on structural medical imaging studies (resonance magnetic studies). Based on these foundations, new research fields to improve the automatic identification of lesions in brain injury can be proposed.

Signal changes in the spinal cord of the rat after injection of formalin into the hindpaw: Characterization using functional magnetic resonance imaging

Changes in metabolism and local circulation occur in the spinal cord during peripheral noxious stimulation. Evidence is presented that this stimulation also causes signal intensity alterations in functional magnetic resonance images of the spinal cord during formalin-induced pain. These results indicate the potential of functional magnetic resonance imaging in assessing noninvasively the extent and intensity of spinal cord excitation in this well characterized pain model. Therefore, the aim of this study was to establish functional magnetic resonance imaging as a noninvasive method to characterize temporal changes in the spinal cord after a single injection of 50 μl of formalin subcutaneously into the hindpaw of the anesthetized rat. This challenge produced a biphasic licking activity in the freely moving conscious animal. Images of the spinal cord were acquired within 2 min, enabling monitoring of the site and the temporal evolution of the signal changes during the development of formalin-induced hyperalgesia without the need of any surgical procedure. The time course of changes in the spinal cord functional image in the isoflurane-anesthetized animal was similar to that obtained from behavioral experiments. Also, comparable physiological data, control experiments, and the inhibition of a response through application of the local anesthetic agent lidocaine indicate that the signal changes observed after formalin injection were specifically related to excitability changes in the relevant segments of the lumbar spinal cord. This approach could be useful to characterize different models of pain and hyperalgesia and, more importantly, to evaluate effects of analgesic drugs.

Detection of cortical activation during averaged single trials of a cognitive task using functional magnetic resonance imaging

Functional neuroimaging studies in human subjects using positron emission tomography or functional magnetic resonance imaging (fMRI) are typically conducted by collecting data over extended time periods that contain many similar trials of a task. Here methods for acquiring fMRI data from single trials of a cognitive task are reported. In experiment one, whole brain fMRI was used to reliably detect single-trial responses in a prefrontal region within single subjects. In experiment two, higher temporal sampling of a more limited spatial field was used to measure temporal offsets between regions. Activation maps produced solely from the single-trial data were comparable to those produced from blocked runs. These findings suggest that single-trial paradigms will be able to exploit the high temporal resolution of fMRI. Such paradigms will provide experimental flexibility and time-resolved data for individual brain regions on a trial-by-trial basis.

31P magnetic resonance spectroscopy of the Sherpa heart: a phosphocreatine/adenosine triphosphate signature of metabolic defense against hypobaric hypoxia.

Of all humans thus far studied, Sherpas are considered by many high-altitude biomedical scientists as most exquisitely adapted for life under continuous hypobaric hypoxia. However, little is known about how the heart is protected in hypoxia. Hypoxia defense mechanisms in the Sherpa heart were explored by in vivo, noninvasive 31P magnetic resonance spectroscopy. Six Sherpas were examined under two experimental conditions [normoxic (21% FiO2) and hypoxic (11% FiO2) and in two adaptational states--the acclimated state (on arrival at low-altitude study sites) and the deacclimating state (4 weeks of ongoing exposure to low altitude). Four lowland subjects were used for comparison. We found that the concentration ratios of phosphocreatine (PCr)/adenosine triphosphate (ATP) were maintained at steady-state normoxic values (0.96, SEM = 0.22) that were about half those found in normoxic lowlanders (1.76, SEM = 0.03) monitored the same way at the same time. These differences in heart energetic status between Sherpas and lowlanders compared under normoxic conditions remained highly significant (P < 0.02) even after 4 weeks of deacclimation at low altitudes. In Sherpas under acute hypoxia, the heart rate increased by 20 beats per min from resting values of about 70 beats per min, and the percent saturation of hemoglobin decreased to about 75%. However, these perturbations did not alter the PCr/ATP concentration ratios, which remained at about 50% of the values expected in healthy lowlanders. Because the creatine phosphokinase reaction functions close to equilibrium, these steady-state PCr/ATP ratios presumably coincided with about 3-fold higher free adenosine diphosphate (ADP) concentrations. Higher ADP concentrations (i.e., lower [PCr]/[ATP] ratios) were interpreted to correlate with the Km values for ADP-requiring kinases of glycolysis and to reflect elevated carbohydrate contributions to heart energy needs. This metabolic organization is postulated as advantageous in hypobaria because the ATP yield per O2 molecule is 25-60% higher with glucose than with free fatty acids (the usual fuels utilized in the human heart in postfasting conditions).

Activation of single whisker barrel in rat brain localized by functional magnetic resonance imaging.

The previously established cortical representation of rat whiskers in layer IV of the cortex contains distinct cylindrical columns of cellular aggregates, which are termed barrels and correlate in a one-to-one relation to whiskers on the contralateral rat face. In the present study, functional magnetic resonance imaging (fMRI) of the rat brain was used to map whisker barrel activation during mechanical up-down movement (+/- 2.5 mm amplitude at 8 Hz) of single/multiple whisker(s). Multislice gradient echo fMRI experiments were performed at 7 T with in-plane image resolution of 220 x 220 microns, slice thickness of 1 mm, and echo time of 16 ms. Highly significant (P < 0.001) and localized contralateral regions of activation were observed upon stimulation of single/multiple whisker(s). In all experiments (n = 10), the locations of activation relative to bregma and midline were highly correlated with the neuroanatomical position of the corresponding whisker barrels, and the results were reproducible intra- and interanimal. Our results indicate that fMRI based on blood oxygenation level-dependent image contrast has the sensitivity to depict activation of a single whisker barrel in the rat brain. This noninvasive technique will supplement existing methods in the study of rat barrel cortex and should be particularly useful for the long-term investigations of central nervous system in the same animal.

Identification of human brain regions underlying responses to resistive inspiratory loading with functional magnetic resonance imaging.

Compensatory ventilatory responses to increased inspiratory loading are essential for adequate breathing regulation in a number of pulmonary diseases; however, the human brain sites mediating such responses are unknown. Midsagittal and axial images were acquired in 11 healthy volunteers during unloaded and loaded (30 cmH2O; 1 cmH2O = 98 Pa) inspiratory breathing, by using functional magnetic resonance imaging (fMRI) strategies (1.5-tesla MR; repetition time, 72 msec; echo time, 45 msec; flip angle, 30 degrees; field of view, 26 cm; slice thickness, 5 mm; number of excitations, 1; matrix, 128 x 256). Digital image subtractions and region of interest analyses revealed significantly increased fMRI signal intensity in discrete areas of the ventral and dorsal pons, interpeduncular nucleus, basal forebrain, putamen, and cerebellar regions. Upon load withdrawal, certain regions displayed a rapid fMRI signal off-transient, while in others, a slower fMRI signal decay emerged. Sustained loading elicited slow decreases in fMRI signal across activated regions, while second application of an identical load resulted in smaller signal increases compared to initial signal responses (P < 0.001). A moderate inspiratory load is associated with consistent regional activation of discrete brain locations; certain of these regions have been implicated in mediation of loaded breathing in animal models. We speculate that temporal changes in fMRI signal may indicate respiratory after-discharge and/or habituation phenomena.

Applications of spin resonance techniques to the study of radiation effects,

Includes bibliography.

Speeding up the EM algorithm for mixture model-based segmentation of magnetic resonance images

Mixture models implemented via the expectation-maximization (EM) algorithm are being increasingly used in a wide range of problems in pattern recognition such as image segmentation. However, the EM algorithm requires considerable computational time in its application to huge data sets such as a three-dimensional magnetic resonance (MR) image of over 10 million voxels. Recently, it was shown that a sparse, incremental version of the EM algorithm could improve its rate of convergence. In this paper, we show how this modified EM algorithm can be speeded up further by adopting a multiresolution kd-tree structure in performing the E-step. The proposed algorithm outperforms some other variants of the EM algorithm for segmenting MR images of the human brain. (C) 2004 Pattern Recognition Society. Published by Elsevier Ltd. All rights reserved.

A novel phantom and method for comprehensive 3-dimensional measurement and correction of geometric distortion in magnetic resonance imaging

A phantom that can be used for mapping geometric distortion in magnetic resonance imaging (MRI) is described. This phantom provides an array of densely distributed control points in three-dimensional (3D) space. These points form the basis of a comprehensive measurement method to correct for geometric distortion in MR images arising principally from gradient field non-linearity and magnet field inhomogeneity. The phantom was designed based on the concept that a point in space can be defined using three orthogonal planes. This novel design approach allows for as many control points as desired. Employing this novel design, a highly accurate method has been developed that enables the positions of the control points to be measured to sub-voxel accuracy. The phantom described in this paper was constructed to fit into a body coil of a MRI scanner, (external dimensions of the phantom were: 310 mm x 310 mm x 310 mm), and it contained 10,830 control points. With this phantom, the mean errors in the measured coordinates of the control points were on the order of 0.1 mm or less, which were less than one tenth of the voxel's dimensions of the phantom image. The calculated three-dimensional distortion map, i.e., the differences between the image positions and true positions of the control points, can then be used to compensate for geometric distortion for a full image restoration. It is anticipated that this novel method will have an impact on the applicability of MRI in both clinical and research settings. especially in areas where geometric accuracy is highly required, such as in MR neuro-imaging. (C) 2004 Elsevier Inc. All rights reserved.

A proposed scheme for comprehensive characterization of the measured geometric distortion in magnetic resonance imaging using a three-dimensional phantom

Recently, a 3-dimensional phantom that can provide a comprehensive, accurate and complete measurement of the geometric distortion in MRI has been developed. In this paper, a scheme for characterizing the measured geometric distortion using the 3-D phantom is described. In the proposed scheme, a number of quantitative measures are developed and used to characterize the geometric distortion. These measures encompass the overall and spatial aspects of the geometric distortion. Two specific types of volume of interest, rectangular parallelepipeds (including cubes) and spheres are considered in the proposed scheme. As an illustration, characterization of the geometric distortion in a Siemens 1.5T Sonata MRI system using the proposed scheme is presented. As shown, the proposed scheme provides a comprehensive assessment of the geometric distortion. The scheme can be potentially used as a standard procedure for the assessment of geometric distortion in MRI. (C) 2004 American Association of Physicists in Medicine.

Magnetic resonance imaging analysis of the upper cervical spine extensor musculature in an asymptomatic cohort: an index of fat within muscle

AIM: To establish a simple method to quantify muscle/fat constituents in cervical muscles of asymptomatic women using magnetic resonance imaging (MRI), and to determine whether there is an age effect within a defined age range. MATERIALS AND METHODS: MRI of the upper cervical spine was performed for 42 asymptomatic women aged 18-45 years. The muscle and fat signal intensities on axial spin echo T1-weighted images were quantitatively classified by taking a ratio of the pixel intensity profiles of muscle against those of intermuscular fat for the rectus capitis posterior major and minor and inferior obliquus capitis muscles bilaterally. Inter- and intra-examiner agreement was scrutinized. RESULTS: The average relative values of fat within the upper cervical musculature compared with intermuscular fat indicated that there were only slight variations in indices between the three sets of muscles. There was no significant correlation between age and fat indices. There were significant differences for the relative fat within the muscle compared with intermuscular fat and body mass index for the right rectus capitis posterior major and right and left inferior obliquus capitis muscles (p = 0.032). Intraclass correlation coefficients for intraobserver agreement ranged from 0.94 to 0.98. Inter-rater agreement of the measurements ranged from 0.75 to 0.97. CONCLUSION: A quantitative measure of muscle/fat constituents has been developed, and results of this study indicate that relative fatty infiltration is not a feature of age in the upper cervical extensor muscles of women aged 18-45 years. (C) 2005 The Royal College of Radiologists. Published by Elsevier Ltd. All rights reserved.

Method for a detailed measurement of image intensity nonuniformity in magnetic resonance imaging

In magnetic resonance imaging (MRI), the MR signal intensity can vary spatially and this spatial variation is usually referred to as MR intensity nonuniformity. Although the main source of intensity nonuniformity arises from B, inhomogeneity of the coil acting as a receiver and/or transmitter, geometric distortion also alters the MR signal intensity. It is useful on some occasions to have these two different sources be separately measured and analyzed. In this paper, we present a practical method for a detailed measurement of the MR intensity nonuniformity. This method is based on the same three-dimensional geometric phantom that was recently developed for a complete measurement of the geometric distortion in MR systems. In this paper, the contribution to the intensity nonuniformity from the geometric distortion can be estimated and thus, it provides a mechanism for estimation of the intensity nonuniformity that reflects solely the spatial characteristics arising from B-1. Additionally, a comprehensive scheme for characterization of the intensity nonuniformity based on the new measurement method is proposed. To demonstrate the method, the intensity nonuniformity in a 1.5 T Sonata MR system was measured and is used to illustrate the main features of the method. (c) 2005 American Association of Physicists in Medicine.

Magnetic properties and Spin dynamics in magnetic Molecular Rings

Molecular nanomagnets are spin clusters whose topology and magnetic interactions can be modulated at the level of the chemical synthesis. They are formed by a small number of transition metal ions coupled by the Heisenberg's exchange interactions. Each cluster is magnetically isolated from its neighbors by organic ligands, making each unit not interacting with the others. Therefore, we can investigate the magnetic properties of an isolated molecular nanomagnet by bulk measurements. The present thesis has been mostly devoted to the experimental investigation of the magnetic properties and spin dynamics of different classes of antiferromagnetic (AF) molecular rings. This study has been exploiting various techniques of investigations, such as Nuclear Magnetic Resonance (NMR), muon spin relaxation (muSR) and SQUiD magnetometry. We investigate the magnetic properties and the phonon-induced relaxation dynamics of the first regular Cr9 antiferromagnetic (AF) ring, which represents a prototype frustrated AF ring. The magnetically-open AF rings like Cr8Cd are model systems for the study of the microscopic magnetic behaviour of finite AF Heisenberg chains. In this type of system the different magnetic behaviour depends length and on the parity of the chain (odd or even). In order to study the local spin densities on the Cr sites, the Cr-NMR spectra was collected at low temperature. The experimental result confirm the theoretical predictions for the spin configuration. Finally, the study of Dy6, the first rare-earth based ring that has been ever synthesized, has been performed by AC-SQuID and muSR measurements. We found that the dynamics is characterized by more than one characteristic correlation time, whose values depend strongly on the applied field.

Three-dimensional modeling of the human eye based on magnetic resonance imaging

PURPOSE. A methodology for noninvasively characterizing the three-dimensional (3-D) shape of the complete human eye is not currently available for research into ocular diseases that have a structural substrate, such as myopia. A novel application of a magnetic resonance imaging (MRI) acquisition and analysis technique is presented that, for the first time, allows the 3-D shape of the eye to be investigated fully. METHODS. The technique involves the acquisition of a T2-weighted MRI, which is optimized to reveal the fluid-filled chambers of the eye. Automatic segmentation and meshing algorithms generate a 3-D surface model, which can be shaded with morphologic parameters such as distance from the posterior corneal pole and deviation from sphericity. Full details of the method are illustrated with data from 14 eyes of seven individuals. The spatial accuracy of the calculated models is demonstrated by comparing the MRI-derived axial lengths with values measured in the same eyes using interferometry. RESULTS. The color-coded eye models showed substantial variation in the absolute size of the 14 eyes. Variations in the sphericity of the eyes were also evident, with some appearing approximately spherical whereas others were clearly oblate and one was slightly prolate. Nasal-temporal asymmetries were noted in some subjects. CONCLUSIONS. The MRI acquisition and analysis technique allows a novel way of examining 3-D ocular shape. The ability to stratify and analyze eye shape, ocular volume, and sphericity will further extend the understanding of which specific biometric parameters predispose emmetropic children subsequently to develop myopia. Copyright © Association for Research in Vision and Ophthalmology.