Development of a compact MRI system for trabecular bone volume fraction measurements

Assessment of Calcaneal Spongy Bone Magnetic Resonance Characteristics in Women: A Comparison between Measures Obtained at 0.3 T, 1.5 T, and 3.0 T

BACKGROUND

There is a growing interest in bone tissue MRI and an even greater interest in using low-cost MR scanners. However, the characteristics of bone MRI remain to be fully defined, especially at low field strength. This study aimed to characterize the signal-to-noise ratio (SNR), T2, and T2* in spongy bone at 0.3 T, 1.5 T, and 3.0 T. Furthermore, relaxation times were characterized as a function of bone-marrow lipid/water ratio content and trabecular bone density.

METHODS

Thirty-two women in total underwent an MR-imaging investigation of the calcaneus at 0.3 T, 1.5 T, and 3.0 T. MR-spectroscopy was performed at 3.0 T to assess the fat/water ratio. SNR, T2, and T2* were quantified in distinct calcaneal regions (ST, TC, and CC). ANOVA and Pearson correlation statistics were used.

RESULTS

SNR increase depends on the magnetic field strength, acquisition sequence, and calcaneal location. T2* was different at 3.0 T and 1.5 T in ST, TC, and CC. Relaxation times decrease as much as the magnetic field strength increases. The significant linear correlation between relaxation times and fat/water found in healthy young is lost in osteoporotic subjects.

CONCLUSION

The results have implications for the possible use of relaxation vs. lipid/water marrow content for bone quality assessment and the development of quantitative MRI diagnostics at low field strength.

Photoacoustic Imaging and Characterization of Bone in Medicine: Overview, Applications, and Outlook

Photoacoustic techniques have shown promise in identifying molecular changes in bone tissue and visualizing tissue microstructure. This capability represents significant advantages over gold standards (i.e., dual-energy X-ray absorptiometry) for bone evaluation without requiring ionizing radiation. Instead, photoacoustic imaging uses light to penetrate through bone, followed by acoustic pressure generation, resulting in highly sensitive optical absorption contrast in deep biological tissues. This review covers multiple bone-related photoacoustic imaging contributions to clinical applications, spanning bone cancer, joint pathologies, spinal disorders, osteoporosis, bone-related surgical guidance, consolidation monitoring, and transsphenoidal and transcranial imaging. We also present a summary of photoacoustic-based techniques for characterizing biomechanical properties of bone, including temperature, guided waves, spectral parameters, and spectroscopy. We conclude with a future outlook based on the current state of technological developments, recent achievements, and possible new directions.

A Fast GPU-optimized 3D MRI Simulator for Arbitrary k-space Sampling

Purpose:

To develop a fast 3D MRI simulator for arbitrary k-space sampling using a graphical processing unit (GPU) and demonstrate its performance by comparing simulation and experimental results in a real MRI system.

Materials and Methods:

A fast 3D MRI simulator using a GeForce GTX 1080 GPU (NVIDIA Corporation, Santa Clara, CA, USA) was developed using C++ and the CUDA 8.0 platform (NVIDIA Corporation). The unique advantage of this simulator was that it could use the same pulse sequence as used in the experiment. The performance of the MRI simulator was measured using two GTX 1080 GPUs and 3D Cones sequences. The MRI simulation results for 3D non-Cartesian sampling trajectories like 3D Cones sequences using a numerical 3D phantom were compared with the experimental results obtained with a real MRI system and a real 3D phantom.

Results:

The performance of the MRI simulator was about 3800–4900 gigaflops for 128- to 4-shot 3D Cones sequences with 256³ voxels, which was about 60% of the performance of the previous MRI simulator optimized for Cartesian sampling calculated for a Cartesian sampling gradient-echo sequence with 256³ voxels. The effects of the static magnetic field inhomogeneity, radio-frequency field inhomogeneity, gradient field nonlinearity, and fast repetition times on the MR images were reproduced in the simulated images as observed in the experimental images.

Conclusion:

The 3D MRI simulator developed for arbitrary k-space sampling optimized using GPUs is a powerful tool for the development and evaluation of advanced imaging sequences including both Cartesian and non-Cartesian k-space sampling.

BlochSolver: A GPU-optimized fast 3D MRI simulator for experimentally compatible pulse sequences

A magnetic resonance imaging (MRI) simulator, which reproduces MRI experiments using computers, has been developed using two graphic-processor-unit (GPU) boards (GTX 1080). The MRI simulator was developed to run according to pulse sequences used in experiments. Experiments and simulations were performed to demonstrate the usefulness of the MRI simulator for three types of pulse sequences, namely, three-dimensional (3D) gradient-echo, 3D radio-frequency spoiled gradient-echo, and gradient-echo multislice with practical matrix sizes. The results demonstrated that the calculation speed using two GPU boards was typically about 7 TFLOPS and about 14 times faster than the calculation speed using CPUs (two 18-core Xeons). We also found that MR images acquired by experiment could be reproduced using an appropriate number of subvoxels, and that 3D isotropic and two-dimensional multislice imaging experiments for practical matrix sizes could be simulated using the MRI simulator. Therefore, we concluded that such powerful MRI simulators are expected to become an indispensable tool for MRI research and development.

Bone volume-to-total volume ratio measured in trabecular bone by single-sided NMR devices

PURPOSE

Reduced bone strength is associated with a loss of bone mass, usually evaluated by dual-energy X-ray absorptiometry, although it is known that the bone microstructure also affects the bone strength. Here, a method is proposed to measure (in laboratory) the bone volume-to-total volume ratio by single-sided NMR scanners, which is related to the microstructure of the trabecular bone.

METHODS

Three single-sided scanners were used on animal bone samples. These low-field, mobile, low-cost devices are able to detect the NMR signal, regardless of the sample sizes, without the use of ionizing radiations, with the further advantage of signal localization offered by their intrinsic magnetic field gradients.

RESULTS

The performance of the different single-sided scanners have been discussed. The results have been compared with bone volume-to-total volume ratio by micro CT and MRI, obtaining consistent values.

CONCLUSIONS

Our results demonstrate the feasibility of the method for laboratory analyses, which are useful for measurements like porosity on bone specimens. This can be considered as the first step to develop an NMR method based on the use of a mobile single-sided device, for the diagnosis of osteoporosis, through the acquisition of the signal from the appendicular skeleton, allowing for low-cost, wide screening campaigns. Magn Reson Med 79:501-510, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Non-destructive Analysis of Oil-Contaminated Soil Core Samples by X-ray Computed Tomography and Low-Field Nuclear Magnetic Resonance Relaxometry: a Case Study

Non-destructive measurements of contaminated soil core samples are desirable prior to destructive measurements because they allow obtaining gross information from the core samples without touching harmful chemical species. Medical X-ray computed tomography (CT) and time-domain low-field nuclear magnetic resonance (NMR) relaxometry were applied to non-destructive measurements of sandy soil core samples from a real site contaminated with heavy oil. The medical CT visualized the spatial distribution of the bulk density averaged over the voxel of 0.31 × 0.31 × 2 mm(3). The obtained CT images clearly showed an increase in the bulk density with increasing depth. Coupled analysis with in situ time-domain reflectometry logging suggests that this increase is derived from an increase in the water volume fraction of soils with depth (i.e., unsaturated to saturated transition). This was confirmed by supplementary analysis using high-resolution micro-focus X-ray CT at a resolution of ∼10 μm, which directly imaged the increase in pore water with depth. NMR transverse relaxation waveforms of protons were acquired non-destructively at 2.7 MHz by the Carr-Purcell-Meiboom-Gill (CPMG) pulse sequence. The nature of viscous petroleum molecules having short transverse relaxation times (T2) compared to water molecules enabled us to distinguish the water-saturated portion from the oil-contaminated portion in the core sample using an M(0)-T2 plot, where M(0) is the initial amplitude of the CPMG signal. The present study demonstrates that non-destructive core measurements by medical X-ray CT and low-field NMR provide information on the groundwater saturation level and oil-contaminated intervals, which is useful for constructing an adequate plan for subsequent destructive laboratory measurements of cores.

Development of a local electromagnetic shielding for an extremity magnetic resonance imaging system

A local radio frequency (rf) shielding consisting of a Cu plate and an LC balun circuit has been developed for a compact magnetic resonance imaging (MRI) system with a 0.3 T permanent magnet. Performance of the local rf shielding was evaluated using an artificial external noise source irradiating a human subject whose hand was inserted into the rf coil of the MRI system. Power spectra of the rf signal detected through the rf coil demonstrated that the local rf shield achieved 30.1 dB external noise suppression. With the local rf shielding, a MRI of the subject's hand was performed using a three-dimensional gradient-echo sequence. Anatomical structures of the subject's hand were clearly visualized. It was concluded that the local rf shielding could be used for the compact MRI system instead of a rf shielded room.

Quantitative MRI of parallel changes of articular cartilage and underlying trabecular bone in degeneration

OBJECTIVE

To determine the interrelations between degenerative changes in articular cartilage and underlying trabecular bone during development of osteoarthritis and to test the ability of quantitative magnetic resonance imaging (MRI) to detect those changes.

METHODS

Human cadaver patellae were investigated with quantitative MRI methods, T(2) and dGEMRIC, at 1.5T. Same measurements for isolated cartilage samples were performed at 9.4T. Bone samples, taken at sites matched with cartilage analyses, were measured with MRI and peripheral quantitative computed tomography (pQCT). Mechanical and quantitative microscopic methods were also utilized for both cartilage and bone samples.

RESULTS

Significant differences were found between the samples with different stages of degeneration in mechanical properties, T(2) at 1.5T and proteoglycan (PG) content of articular cartilage. dGEMRIC at 9.4T discerned samples with advanced degeneration from the others. Bone variables measured with pQCT discerned samples with no or minimal and advanced degeneration, and mechanical properties of trabecular bone discerned samples with no or minimal degeneration from the others. Significant linear correlations were found between the bone and cartilage parameters. Characteristically, associations between variables were stronger within the samples with no or minimal degeneration compared to all samples.

CONCLUSIONS

Quantitative MRI variables, especially T(2) relaxation time of articular cartilage, may be feasible surrogate markers for early and advanced osteoarthritic changes in joint tissues, including decreased elastic moduli, PG and collagen contents of cartilage and mineral density and volume fraction of trabecular bone. Further work is required to resolve the relaxation mechanisms at clinically applicable field strengths.

Single-chip pulse programmer for magnetic resonance imaging using a 32-bit microcontroller

A magnetic resonance imaging (MRI) pulse programmer has been developed using a single-chip microcontroller (ADmicroC7026). The microcontroller includes all the components required for the MRI pulse programmer: a 32-bit RISC CPU core, 62 kbytes of flash memory, 8 kbytes of SRAM, two 32-bit timers, four 12-bit DA converters, and 40 bits of general purpose I/O. An evaluation board for the microcontroller was connected to a host personal computer (PC), an MRI transceiver, and a gradient driver using interface circuitry. Target (embedded) and host PC programs were developed to enable MRI pulse sequence generation by the microcontroller. The pulse programmer achieved a (nominal) time resolution of approximately 100 ns and a minimum time delay between successive events of approximately 9 micros. Imaging experiments using the pulse programmer demonstrated the effectiveness of our approach.

Development of a compact MRI system for measuring the trabecular bone microstructure of the finger

A compact MRI system for measuring the trabecular bone (TB) microstructure of the finger using a high-field-strength (1.0T) permanent magnet was developed. The entire system was installed in a 0.6 mx1.2 m space. One male and 36 female subjects participated in the imaging experiments. The TB of the distal phalanx of the middle finger was imaged at a voxel resolution of (160 microm)3 using a three-dimensional (3D) driven equilibrium spin-echo (SE) imaging sequence (imaging time=approximately 14 min). The image data sets obtained yielded two distinct peaks for the bone and marrow in image intensity histograms when no motion was present. The structural parameters obtained through 3D image analysis show that this compact system is potentially useful for evaluating bone quality.

Development of a compact MRI system for trabecular bone microstructure measurements of the distal radius

A compact MRI system for trabecular bone (TB) microstructure measurements of the distal radius was developed using a 1.0 T permanent magnet and a compact MRI console. TB microstructure of the distal radius was clearly visualized using a three-dimensional (3D) driven equilibrium spin-echo (DESE) sequence in 23 min. The image obtained had a sufficient spatial resolution (150 microm x 150 microm x 500 microm) and signal-to-noise ratio (SNR) (approximately 10) for 3D bone microstructure analysis. The system demonstrated the feasibility of using a permanent magnet compact MRI system as a clinical instrument for bone microstructure measurements of the distal radius.

Prediction of mechanical properties of trabecular bone using quantitative MRI

Techniques for quantitative magnetic resonance imaging (MRI) have been developed for non-invasive estimation of the mineral density and structure of trabecular bone. The R*(2) relaxation rate (i.e. 1/T*(2)) is sensitive to bone mineral density (BMD) via susceptibility differences between trabeculae and bone marrow, and by binarizing MRI images, structural variables, such as apparent bone volume fraction, can be assessed. In the present study, trabecular bone samples of human patellae were investigated in vitro at 1.5 T to determine the ability of MRI-derived variables (R*(2) and bone volume fraction) to predict the mechanical properties (Young's modulus, yield stress and ultimate strength). Further, the MRI variables were correlated with reference measurements of volumetric BMD and bone area fraction as determined with a clinical pQCT system. The MRI variables correlated significantly (p < 0.01) with the mechanical variables (r = 0.32-0.46), BMD (r = 0.56) and bone structure (r = 0.51). A combination of R*(2) and MRI-derived bone volume fraction further improved the prediction of yield stress and ultimate strength. Although pQCT showed a trend towards better prediction of the mechanical properties, current results demonstrate the feasibility of combined MR imaging of marrow susceptibility and bone volume fraction in predicting the mechanical strength of trabecular bone and bone mineral density.

Quantitative MRI for the assessment of bone structure and function

Osteoporosis is the most common degenerative disease in the elderly. It is characterized by low bone mass and structural deterioration of bone tissue, leading to morbidity and increased fracture risk in the hip, spine and wrist-all sites of predominantly trabecular bone. Bone densitometry, currently the standard methodology for diagnosis and treatment monitoring, has significant limitations in that it cannot provide information on the structural manifestations of the disease. Recent advances in imaging, in particular MRI, can now provide detailed insight into the architectural consequences of disease progression and regression in response to treatment. The focus of this review is on the emerging methodology of quantitative MRI for the assessment of structure and function of trabecular bone. During the past 10 years, various approaches have been explored for obtaining image-based quantitative information on trabecular architecture. Indirect methods that do not require resolution on the scale of individual trabeculae and therefore can be practiced at any skeletal location, make use of the induced magnetic fields in the intertrabecular space. These fields, which have their origin in the greater diamagnetism of bone relative to surrounding marrow, can be measured in various ways, most typically in the form of R2', the recoverable component of the total transverse relaxation rate. Alternatively, the trabecular network can be quantified by high-resolution MRI (micro-MRI), which requires resolution adequate to at least partially resolve individual trabeculae. Micro-MRI-based structure analysis is therefore technically demanding in terms of image acquisition and algorithms needed to extract the structural information under conditions of limited signal-to-noise ratio and resolution. Other requirements that must be met include motion correction and image registration, both critical for achieving the reproducibility needed in repeat studies. Key clinical applications targeted involve fracture risk prediction and evaluation of the effect of therapeutic intervention.

Low-field compact magnetic resonance imaging system for the hand and wrist in rheumatoid arthritis

PURPOSE

To investigate the feasibility of an originally developed compact MRI system for evaluating rheumatoid arthritis (RA), and determine its advantages and disadvantages as an imaging modality for evaluating RA.

MATERIALS AND METHODS

We prospectively studied 13 healthy controls with no clinical symptoms of arthritis, and 13 patients with hand and wrist pains (including pain from RA) with a 0.2 T permanent-magnet compact MR imager. All MR images were obtained while the subjects were in a sitting position. Coronal three-dimensional spin-echo T1-weighted images and coronal two-dimensional short tau inversion recovery (STIR) images were obtained with image matrix = 256 x 128 and field of view (FOV) = 20.48 cm. Plain radiograph findings and MRI findings of patients were compared.

RESULTS

In three of the patients with suspected early RA (N = 7), early RA was evaluated based on STIR images. All RA patients showed morphologic or signal intensity changes that allowed an evaluation of RA from MR findings. Four of five RA patients showed high signal intensity on STIR images in the wrist, proximal interphalangeal (PIP) joint, or metacarpophalangeal (MCP) joint that suggested synovitis. Multiple erosions in the hand and wrist were seen in four RA patients, with low signal intensity on T1-weighted images.

CONCLUSION

RA was correctly evaluated, and early RA could be identified with the compact MRI system. However, the current system has limitations, such as the nonselective STIR sequence used and magnetic field inhomogeneity.

Trabecular bone volume fraction measurements of a large number of subjects using a compact MRI

Trabecular bone volume fraction (TBVF) and speed of sound (SOS) were measured for the right calcanei of 416 female volunteers. The TBVF was measured with a compact MRI developed in our laboratory. The SOS was measured with a commercial quantitative ultrasound system. It was observed that the correlation coefficient between TBVF and SOS and that between TBVF and age varied depending on the location of region of interest (ROI) in the calcaneus. As a result, an optimum circular ROI with a diameter of 20 mm was determined so that the correlation coefficients were maximized. In the optimum ROI, transverse relaxation time (T2) of the bone marrow protons of the calcaneus was found to be concentrated in a narrow range over the subjects. This result suggested that a 50% scan time reduction in the TBVF measurements could be made by skipping the T2 correction procedure.