Internal three-dimensional strains in human intervertebral discs under axial compression quantified noninvasively by magnetic resonance imaging and image registration.

Study objectives were to develop, validate, and apply a method to measure three-dimensional (3D) internal strains in intact human discs under axial compression. A custom-built loading device applied compression and permitted load-relaxation outside of the magnet while also maintaining compression and hydration during imaging. Strain was measured through registration of 300 µm isotropic resolution images. Excellent registration accuracy was achieved, with 94% and 65% overlap of disc volume and lamellae compared to manual segmentation, and an average Hausdorff, a measure of distance error, of 0.03 and 0.12 mm for disc volume and lamellae boundaries, respectively. Strain maps enabled qualitative visualization and quantitative regional annulus fibrosus (AF) strain analysis. Axial and circumferential strains were highest in the lateral AF and lowest in the anterior and posterior AF. Radial strains were lowest in the lateral AF, but highly variable. Overall, this study provided new methods that will be valuable in the design and evaluation surgical procedures and therapeutic interventions.

Human L3L4 intervertebral disc mean 3D shape, modes of variation, and their relationship to degeneration.

Intervertebral disc mechanics are affected by both disc shape and disc degeneration, which in turn each affect the other; disc mechanics additionally have a role in the etiology of disc degeneration. Finite element analysis (FEA) is a favored tool to investigate these relationships, but limited data for intervertebral disc 3D shape has forced the use of simplified or single-subject geometries, with the effect of inter-individual shape variation investigated only in specialized studies. Similarly, most data on disc shape variation with degeneration is based on 2D mid-sagittal images, which incompletely define 3D shape changes. Therefore, the objective of this study was to quantify inter-individual disc shape variation in 3D, classify this variation into independently-occurring modes using a statistical shape model, and identify correlations between disc shape and degeneration. Three-dimensional disc shapes were obtained from MRI of 13 human male cadaver L3L4 discs. An average disc shape and four major modes of shape variation (representing 90% of the variance) were identified. The first mode represented disc axial area and was significantly correlated to degeneration (R(2)=0.44), indicating larger axial area in degenerate discs. Disc height variation occurred in three distinct modes, each also involving non-height variation. The statistical shape model provides an average L3L4 disc shape for FEA that is fully defined in 3D, and makes it convenient to generate a set of shapes with which to represent aggregate inter-individual variation. Degeneration grade-specific shapes can also be generated. To facilitate application, the model is included in this paper׳s supplemental content.

Evaluation of intervertebral disc cartilaginous endplate structure using magnetic resonance imaging.

PURPOSE: The cartilaginous endplate (CEP) is a thin layer of hyaline cartilage positioned between the vertebral endplate and nucleus pulposus (NP) that functions both as a mechanical barrier and as a gateway for nutrient transport into the disc. Despite its critical role in disc nutrition and degeneration, the morphology of the CEP has not been well characterized. The objective of this study was to visualize and report observations of the CEP three-dimensional morphology, and quantify CEP thickness using an MRI FLASH (fast low-angle shot) pulse sequence. METHODS: MR imaging of ex vivo human cadaveric lumbar spine segments (N = 17) was performed in a 7T MRI scanner with sequence parameters that were selected by utilizing high-resolution T1 mapping, and an analytical MRI signal model to optimize image contrast between CEP and NP. The CEP thickness at five locations along the mid-sagittal AP direction (center, 5 mm, 10 mm off-center towards anterior and posterior) was measured, and analyzed using two-way ANOVA and a post hoc Bonferonni test. For further investigation, six in vivo volunteers were imaged with a similar sequence in a 3T MRI scanner. In addition, decalcified and undecalcified histology was performed, which confirmed that the FLASH sequence successfully detected the CEP. RESULTS: CEP thickness determined by MRI in the mid-sagittal plane across all lumbar disc levels and locations was 0.77 ± 0.24 mm ex vivo. The CEP thickness was not different across disc levels, but was thinner toward the center of the disc. CONCLUSIONS: This study demonstrates the potential of MRI FLASH imaging for structural quantification of the CEP geometry, which may be developed as a technique to evaluate changes in the CEP with disc degeneration in future applications.

Superelliptical insert gradient coil with a field-modifying layer for breast imaging.

Many MRI applications such as dynamic contrast-enhanced MRI of the breast require high spatial and temporal resolution and can benefit from improved gradient performance, e.g., increased gradient strength and reduced gradient rise time. The improved gradient performance required to achieve high spatial and temporal resolution for this application may be achieved by using local insert gradients specifically designed for a target anatomy. Current flat gradient systems cannot create an imaging volume large enough to accommodate both breasts; further, their gradient fields are not homogeneous, dropping off rapidly with distance from the gradient coil surface. To attain an imaging volume adequate for bilateral breast MRI, a planar local gradient system design has been modified into a superellipse shape, creating homogeneous gradient volumes that are 182% (Gx), 57% (Gy), and 75% (Gz) wider (left/right direction) than those of the corresponding standard planar gradient. Adding an additional field-modifying gradient winding results in an additional improvement of the homogeneous gradient field near the gradient coil surface over the already enlarged homogeneous gradient volumes of the superelliptical gradients (67%, 89%, and 214% for Gx, Gy, and Gz respectively). A prototype y-gradient insert has been built to demonstrate imaging and implementation characteristics of the superellipse gradient in a 3 T MRI system.

Design, Fabrication and Testing of an Insertable Double-Imaging-Region Gradient Coil.

We have constructed a small-bore insertable gradient coil with two linear gradient imaging regions and interfaced it with an MRI scanner. We have also constructed an RF system capable of transmitting or receiving in both regions simultaneously.Designs for conductor placement for two-region X-, Y- and Z-gradient coils were optimized by simulated annealing. Wire patterns for each axis were chosen that gave low inductance, reasonable homogeneity over a large imaging volume and high efficiency (gradient field per-unit-current).Imaging was performed on a Siemens 3T TIM Trio scanner equipped with three additional gradient amplifier channels and a second RF/gradient array controller. Phantoms were placed in the two imaging regions as well as the central non-imaging region to test gradient homogeneity and crosstalk between regions. Images acquired simultaneously in the two regions showed very little signal crosstalk between imaging regions and even less signal from the central, non-imaging region.When combined with an overlapping single-region gradient insert, extended field-of-view (FOV) imaging will be possible without moving the table or the subject and without increasing nerve stimulation. Construction and testing of a two-region gradient coil insert is a necessary intermediate step as a proof of concept for an extended field of view, contiguous, three-region human-sized gradient system.

Magnetic Resonance Imaging with Composite (Dual) Gradients.

The tradeoff between gradient performance factors, size of the imaging region, and physiological factors such as nerve stimulation typically leads to compromises in gradient design and ultimately suboptimal imaging performance. Local gradient systems can add some performance flexibility, but are cumbersome to set up and remove. In nearly all conventional MRI systems, the use of local gradients precludes the use of the more homogeneous whole body gradients. This paper presents the concept of dynamically selectable composite gradient systems where local gradients and whole body gradients can be selected independently and simultaneously. The relative performance of whole body, insert, and composite gradients is predicted for echoplanar (EPI), turbo spin echo (TSE), and steady state free precession (SSFP). A realization of the concept is presented.

Automatic detection of three-dimensional vascular tree centerlines and bifurcations in high-resolution magnetic resonance angiography.

OBJECTIVES: We sought to develop a simple and robust algorithm capable of automatically detecting centerlines and bifurcations of a three-dimensional (3D) vascular bed. MATERIALS AND METHODS: After necessary preprocessing, an appropriate cost function is computed for all vessel voxels and Dijkstra's minimum-cost-path algorithm is implemented. By back tracing all the minimum-cost paths, centerlines and bifurcation are detected. The detected paths are then split into segments between adjacent nodes (bifurcations or vessel end-points) and smoothed by curve fitting. RESULTS: Application of the algorithm to both simulated 3D vessels and 3D magnetic resonance angiography (MRA) images of an actual intracranial arterial tree produced well-centered vessel skeletons. Quantitative assessment of the algorithm was performed. For the simulated data, the root mean square error for centerline detection is about half a voxel. For the human intracranial MRA data, the sensitivity, positive predictive value (PPV), and accuracy of bifurcation detection were calculated for different cost functions. The best case gave a sensitivity of 91.4%, a PPV of 91.4%, and an RMS error of 1.7 voxels. CONCLUSIONS: To the extent that imperfections are eliminated from the segmented image, the algorithm is effective and robust in automatic and accurate detection of centerlines and bifurcations. The cost function and algorithm used are demonstrated to be an improvement over similar algorithms in the literature.