High-resolution, large dynamic range field map estimation

Patient specific distortion detection and mitigation in MR images used for stereotactic radiosurgery

Objective. In MRI-based radiation therapy planning, mitigating patient-specific distortion with standard high bandwidth scans can result in unnecessary sacrifices of signal to noise ratio. This study investigates a technique for distortion detection and mitigation on a patient specific basis. Approach. Fast B0 mapping was performed using a previously developed technique for high-resolution, large dynamic range field mapping without the need for phase unwrapping algorithms. A phantom study was performed to validate the method. Distortion mitigation was validated by reducing geometric distortion with increased acquisition bandwidth and confirmed by both the B0 mapping technique and manual measurements. Images and contours from 25 brain stereotactic radiosurgery patients and 95 targets were analyzed to estimate the range of geometric distortions expected in the brain and to estimate bandwidth required to keep all treatment targets within the ±0.5 mm iso-distortion contour. Main Results. The phantom study showed, at 3 T, the technique can measure distortions with a mean absolute error of 0.12 mm (0.18 ppm), and a maximum error of 0.37 mm (0.6 ppm). For image acquisition at 3 T and 1.0 mm resolution, mean absolute distortion under 0.5 mm in patients required bandwidths from 109 to 200 Hz px−1 for patients with the least and most distortion, respectively. Maximum absolute distortion under 0.5 mm required bandwidths from 120 to 390 Hz px−1. Significance. The method for B0 mapping was shown to be valid and may be applied to assess distortion clinically. Future work will adapt the readout bandwidth to prospectively mitigate distortion with the goal to improve radiosurgery treatment outcomes by reducing healthy tissue exposure.

Multi-echo gradient-recalled-echo phase unwrapping using a Nyquist sampled virtual echo train in the presence of high-field gradients

PURPOSE

To develop a spatio-temporal approach to accurately unwrap multi-echo gradient-recalled echo phase in the presence of high-field gradients.

THEORY AND METHODS

Using the virtual echo-based Nyquist sampled (VENyS) algorithm, the temporal unwrapping procedure is modified by introduction of one or more virtual echoes between the first lower and the immediate higher echo, so as to reinstate the Nyquist condition at locations with high-field gradients. An iterative extension of the VENyS algorithm maintains spatial continuity by adjusting the phase rotations to make the neighborhood phase differences less than π. The algorithm is evaluated using simulated data, Gadolinium contrast-doped phantom, and in vivo brain, abdomen, and chest data sets acquired at 3 T and 9.4 T. The unwrapping performance is compared with the standard temporal unwrapping algorithm used in the morphology-enabled dipole inversion-QSM pipeline as a benchmark for validation.

RESULTS

Quantitative evaluation using numerical phantom showed significant reduction in unwrapping errors in regions of large field gradients, and the unwrapped phase revealed an exact match with the linear concentration profile of vials in a gadolinium contrast-doped phantom data acquired at 9.4 T. Without the need for additional spatial unwrapping, the iterative VENyS algorithm was able to generate spatially continuous phase images. Application to in vivo data resulted in better unwrapping performance, especially in regions with large susceptibility changes such as the air/tissue interface.

CONCLUSION

The iterative VENyS algorithm serves as a robust unwrapping method for multi-echo gradient-recalled echo phase in the presence of high-field gradients.

Iterative Restoration of the Fringe Phase (REFRASE) for QSM

In quantitative susceptibility mapping (QSM), reconstructed results can be critically biased by misinterpreted or missing phase data near the edges of the brain support originating from the non-local relationship between field and susceptibility. These data either have to be excluded or corrected before further processing can take place. To address this, our iterative restoration of the fringe phase (REFRASE) approach simultaneously enhances the accuracy of multi-echo phase data QSM maps and the extent of the area available for evaluation. Data loss caused by strong local phase gradients near the surface of the brain support is recovered within the original phase data using harmonic and dipole-based fields extrapolated from a robust support region toward an extended brain mask. Over several iterations, phase data are rectified prior to the application of further QSM processing steps. The concept is successfully validated on numerical phantoms and brain scans from a cohort of volunteers. The increased extent of the mask and improved numerical stability within the segmented globus pallidus confirm the efficacy of the presented method in comparison to traditional evaluation.

B0 shimming for in vivo magnetic resonance spectroscopy: Experts' consensus recommendations

Magnetic resonance spectroscopy (MRS) and spectroscopic imaging (MRSI) allow the chemical analysis of physiological processes in vivo and provide powerful tools in the life sciences and for clinical diagnostics. Excellent homogeneity of the static B₀ magnetic field over the object of interest is essential for achieving high-quality spectral results and quantitative metabolic measurements. The experimental minimization of B₀ variation is performed in a process called B₀ shimming. In this article, we summarize the concepts of B₀ field shimming using spherical harmonic shimming techniques, specific strategies for B₀ homogenization and crucial factors to consider for implementation and use in both brain and body. In addition, experts' recommendations are provided for minimum requirements for B₀ shim hardware and evaluation criteria for the primary outcome of adequate B₀ shimming for MRS and MRSI, such as the water spectroscopic linewidth.

Exploratory study of geometric distortion correction of prostate diffusion-weighted imaging using B0 map acquisition

BACKGROUND

Evaluation of prostate MRI relies on diffusion-weighted imaging (DWI), commonly distorted by susceptibility artifacts, thereby creating a need for approaches to correct such distortion.

PURPOSE

To compare geometric distortion on prostate MRI between standard DWI and a geometric distortion correction method for DWI described as static distortion correction DWI (SDC DWI).

STUDY TYPE

Retrospective case study.

POPULATION

Thirty patients (ages 31-81 years) undergoing prostate MRI.

SEQUENCE

Geometric distortions from echo planar imaging were corrected with the SDC DWI protocol, which first acquires a B₀ -field map to estimate geometric distortions.

ASSESSMENT

Contours of the prostate were placed on axial T₂ -weighted imaging (T₂ WI) as an anatomic standard. Pixel shifts and apparent diffusion coefficient (ADC) values were compared between prostate contours applied to the SDC DWI and standard DWI sequences. Detailed characterization of the impact of SDC DWI was performed in three representative patients.

STATISTICAL TESTS

One-way analysis of variance (ANOVA) test, Spearman correlation test, and Bland-Altman plots were calculated.

RESULTS

There was significantly greater overlap of the SDC DWI prostate region of interest (ROI) with T₂ WI than standard DWI with T₂ WI (10.56 cm² ± 3.14, P < 0.05). R² of ADC values from standard DWI vs. SDC DWI in the 30 patients ranged from 0.02-0.94 (mean 0.60). A patient without susceptibility artifact demonstrated minimal pixel shift ranging from 0.6-1.3 mm and high correlation of ADC values (R² = 0.89) between SDC DWI and standard DWI. A patient with rectal gas showed greater pixel shift (range: -2.5 to -0.5 mm) and less ADC value correlation (R² = 0.69). A patient with a pelvic phlebolith adjacent to the prostate showed an even greater pixel shift (range: 10-16 mm) and decreased ADC correlation (R² = 0.21).

DATA CONCLUSION

SDC DWI appears to correct for susceptibility-related pixel shifts in the prostate compared with standard DWI, which may have value for assessing prostate lesions obscured by geometric warping. Level of Evidence 4 Technical Efficacy Stage 1 J. Magn. Reson. Imaging 2019;50:1614-1619.

The use of Fourier-domain analyses for unwrapping phase images of low SNR

PURPOSE

We report a new postprocessing procedure that uses Fourier-domain data analyses to improve the accuracy and reliability of phase unwrapping for MRI data of low SNR.

METHODS

The developed method first identifies the Fourier-domain energy peak locations corresponding to different image-domain areas from which a robust measurement of image-domain phase gradients can be obtained even for MRI data of low SNR. The phase-gradient information measured from critical brain regions using the above-mentioned Fourier-domain analysis is then combined with the conventional temporal-domain or spatial-domain phase-unwrapping procedure to remove phase wraps. The developed method was tested with MRI data obtained from 30 healthy adult volunteers, and its performance was quantitatively evaluated.

RESULTS

The developed Fourier-domain analysis could robustly quantify image-domain phase gradients even for MRI data with low SNR (e.g., SNR ≃ 2). Experimental results show that the Fourier-domain analyses could further reduce phase wrap artifact in data produced by the conventional temporal-domain or spatial-domain phase-unwrapping procedures.

CONCLUSION

Our results demonstrate that the developed phase-unwrapping method can reduce residual phase wraps resulting from conventional procedures in critical brain regions (e.g., near the air-tissue interfaces) and should prove valuable for studies that require accurate measurements of MRI phase values, such as QSM, B₀ field mapping, and temperature mapping.

Seed prioritized unwrapping (SPUN) for MR phase imaging

BACKGROUND

Region-growing-based phase unwrapping methods have the potential for lossless phase aliasing removal, but generally suffer from unwrapping error propagation associated with discontinuous phase and/or long calculation times. The tradeoff point between robustness and efficiency of phase unwrapping methods in the region-growing category requires improvement.

PURPOSE

To demonstrate an accurate, robust, and efficient region-growing phase unwrapping method for MR phase imaging applications.

STUDY TYPE

Prospective. SUBJECTS, PHANTOM: normal human subjects (10) / brain surgery patients (2) / water phantoms / computer simulation.

FIELD STRENGTH/SEQUENCE

3 T/gradient echo sequences (2D and 3D).

ASSESSMENT

A seed prioritized unwrapping (SPUN) method was developed based on single-region growing, prioritizing only a portion (eg, 100 seeds or 1% seeds) of available seed voxels based on continuity quality during each region-growing iteration. Computer simulation, phantom, and in vivo brain and pelvis scans were performed. The error rates, seed percentages, and calculation times were recorded and reported. SPUN unwrapped phase images were visually evaluated and compared with Laplacian unwrapped results.

STATISTICAL TESTS

Monte Carlo simulation was performed on a 3D dipole phase model with a signal-to-noise ratio (SNR) of 1-9 dB, to obtain the mean and standard deviation of calculation error rates and calculation times.

RESULTS

Simulation revealed a very robust unwrapping performance of SPUN, reaching an error rate of <0.4% even with SNR as low as 1 dB. For all in vivo data, SPUN was able to robustly unwrap the phase images of modest SNR and complex morphology with visually minimal errors and fast calculation speed (eg, <4 min for 368 × 312 × 128 data) when using a proper seed priority number, eg, N_sp = 1 or 10 voxels for 2D and N_sp = 1% for 3D data.

DATA CONCLUSION

SPUN offers very robust and fast region-growing-based phase unwrapping, and does not require any tissue masking or segmentation, nor poses a limitation over imaging parameters.

LEVEL OF EVIDENCE

3 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2019;50:62-70.

Robust selective weighted field mapping using multi-echo gradient echo-based MRI

In 3D gradient echo (GRE) and echo planar imaging (EPI), strong macroscopic field gradients are observed at air/tissue interfaces. The respective field gradients lead to an apparent increase in intravoxel dephasing, and, subsequently, to signal loss or image distortion. We propose an analytical approximation and a consequent method to compute low and high resolution field maps over all field map regimes (small and large echo spacing). A number of approaches which compute field maps from reconstructed phase data rely upon optimized linear least square fit and complex division approaches owing to the simplicity of their implementation. Most of these techniques, however, have historically considered only the phase signal when computing off-resonance maps while ignoring magnitude data. This latter may be of notable interest since the presence of noise is well depicted and interpreted. The presence of noise and phase aliasing that increase with increasing echo time (TE) and echo spacing (ΔTE) may seriously challenge the off-resonance map accuracy. These techniques still remain subject to the trade-off during the choice of GRE sequences, TE and ΔTE. In this work, we explore a novel model that considers any type of TE and ΔTE regime (small or large) and high phase wraps complexity. The field offset is weighted by the magnitude signal decay quality, to make the field mapping procedure as noise independent as possible. The performance of the proposed method was tested using simulated, experimental phantoms and in vivo human studies. The proposed approach markedly outperforms conventional techniques. It provides a correction equivalent to that of the conventional techniques in regions with high SNR ([Formula: see text]20), yielding a mean error of about 0.1 Hz, but appearing more robust in regions with low SNR ([Formula: see text]10), such as near the sinus cavity and at the very edge of the brain (mean error less than 1 Hz), where phase wraps and noise are highly present. The proposed technique shows promise to enhance field map generation over any acquisition regime and in regions of both high and low SNR and it can be easily implemented for rapid computation and used in a clinical setting.

A Robust Method for Estimating B0 Inhomogeneity Field in the Liver by Mitigating Fat Signals and Phase-Wrapping

We developed an optimized and robust method to estimate liver B0 field inhomogeneity for monitoring and correcting susceptibility-induced geometric distortion in magnetic resonance images for precision therapy. A triple-gradient-echo acquisition was optimized for the whole liver B0 field estimation within a single-exhale breath-hold scan on a 3 T scanner. To eliminate chemical-shift artifacts, fat signals were chosen in-phase between 2 echoes with an echo time difference (ΔTE) of 2.3 milliseconds. To avoid phase-wrapping, other 2 echoes provided a large field dynamic range (1/ΔTE) to cover the B0 field inhomogeneity. In addition, using high parallel imaging factor of 4 and a readout-bandwidth of 1955 Hz/pixel, an ∼18-second acquisition time for breath-held scans was achieved. A 2-step, 1-dimensional regularized method for the ΔB0 field map estimation was developed, tested and validated in phantom and patient studies. Our method was validated on a water phantom with fat components and air pockets; it yielded ΔB0-field maps that had no chemical-shift and phase-wrapping artifacts, and it had a <0.5 mm of geometric distortion near the air pockets. The ΔB0-field maps of the patients' abdominal regions were also free from phase-wrapping and chemical-shift artifacts. The maximum field inhomogeneity was found near the lung–liver interface, up to ∼300 Hz, resulting in ∼2 mm of distortions in anatomical images with a readout-bandwidth of 440 Hz/pixel. The field mapping method in the abdominal region is robust; it can be easily integrated in clinical workflow for patient-based quality control of magnetic resonance imaging geometric integrity.

MR phase imaging with bipolar acquisition

We have previously proposed a novel magnetic resonance (MR) phase imaging framework (MAGPI) based on a three-echo sequence that demonstrated substantial gains in phase signal-to-noise ratio (SNR). We improve upon the performance of MAGPI by extending the formulation to handle (i) an alternating gradient polarity (bipolar) readout scheme and (ii) an arbitrary number of echoes. We formulate the phase-imaging problem using maximum-likelihood (ML) estimation. The acquisition uses an optimized multi-echo gradient echo (MEGE) sequence. The tissue-phase estimation algorithm is a voxel-per-voxel approach, which requires no reference scans, no phase unwrapping and no spatial denoising. Unlike other methods, our bipolar readout model is general and does not make simplifying assumptions about the even-odd echo phase errors. The results show that (a) our proposed bipolar MAGPI approach improves on the phase SNR gains achieved with monopolar MAGPI and (b) the phase SNR converges with the number of echoes more rapidly with bipolar MAGPI. Importantly, bipolar MAGPI enables phase imaging in severely SNR-constrained scenarios, where monopolar MAGPI is unable to find solutions. The substantial phase SNR gains achieved with our framework are used here to (a) accelerate acquisitions (full brain 0.89 mm in-plane resolution in 2 min 30 sec) and (b) enable high-contrast high-resolution phase imaging (310 µm in-plane resolution) at clinical field strengths. Copyright © 2016 John Wiley & Sons, Ltd.

An illustrated comparison of processing methods for MR phase imaging and QSM: combining array coil signals and phase unwrapping

Phase imaging benefits from strong susceptibility effects at very high field and the high signal-to-noise ratio (SNR) afforded by multi-channel coils. Combining the information from coils is not trivial, however, as the phase that originates in local field effects (the source of interesting contrast) is modified by the inhomogeneous sensitivity of each coil. This has historically been addressed by referencing individual coil sensitivities to that of a volume coil, but alternative approaches are required for ultra-high field systems in which no such coil is available. An additional challenge in phase imaging is that the phase that develops up to the echo time is "wrapped" into a range of 2π radians. Phase wraps need to be removed in order to reveal the underlying phase distribution of interest. Beginning with a coil combination using a homogeneous reference volume coil - the Roemer approach - which can be applied at 3 T and lower field strengths, we review alternative methods for combining single-echo and multi-echo phase images where no such reference coil is available. These are applied to high-resolution data acquired at 7 T and their effectiveness assessed via an index of agreement between phase values over channels and the contrast-to-noise ratio in combined images. The virtual receiver coil and COMPOSER approaches were both found to be computationally efficient and effective. The main features of spatial and temporal phase unwrapping methods are reviewed, placing particular emphasis on recent developments in temporal phase unwrapping and Laplacian approaches. The features and performance of these are illustrated in application to simulated and high-resolution in vivo data. Temporal unwrapping was the fastest of the methods tested and the Laplacian the most robust in images with low SNR. © 2016 The Authors. NMR in Biomedicine published by John Wiley & Sons Ltd.

Foundations of MRI phase imaging and processing for Quantitative Susceptibility Mapping (QSM)

Quantitative Susceptibility Mapping (QSM) is a novel MRI based technique that relies on estimates of the magnetic field distribution in the tissue under examination. Several sophisticated data processing steps are required to extract the magnetic field distribution from raw MRI phase measurements. The objective of this review article is to provide a general overview and to discuss several underlying assumptions and limitations of the pre-processing steps that need to be applied to MRI phase data before the final field-to-source inversion can be performed. Beginning with the fundamental relation between MRI signal and tissue magnetic susceptibility this review covers the reconstruction of magnetic field maps from multi-channel phase images, background field correction, and provides an overview of state of the art QSM solution strategies.

Automated slice-specific simultaneous z-shim method for reducing B1 inhomogeneity and susceptibility-induced signal loss with parallel transmission at 3T

PURPOSE

Through-plane susceptibility-induced signal loss in gradient recalled echo (GRE)-based sequences can considerably impair both the clinical diagnosis and functional analysis of certain brain areas. In this work, a fully automated simultaneous z-shim approach is proposed on the basis of parallel transmit (pTX) to reduce those signal dropouts at 3T.

THEORY AND METHODS

The approach uses coil-specific time-delayed excitations to impose a z-shim phase. It was extended toward B1 inhomogeneity mitigation and adequate slice-specific signal-dephasing cancellation on the basis of the prevailing B0 and B1 spatial information. Local signal recovery level and image quality preservation were analyzed using multi-slice FLASH experiments in humans and compared to the standard excitation. Spatial blood-oxygen-level-dependent (BOLD) activation coverage was further compared in breath-hold functional MRI.

RESULTS

The pTX z-shim approach recovered approximately 47% of brain areas affected by signal loss in standard excitation images across all subjects. At the same time, B1 shading effects could be substantially reduced. In these areas, BOLD activation coverage could be also increased by approximately 57%.

CONCLUSION

The proposed fully automated pTX z-shim method enables time-efficient and robust signal recovery in GRE-based sequences on a clinical scanner using two standard whole-body transmit coils.

Current progress in nanotechnology applications for diagnosis and treatment of kidney diseases

Significant progress has been made in nanomedicine, primarily in the form of nanoparticles, for theranostic applications to various diseases. A variety of materials, both organic and inorganic, have been used to develop nanoparticles with promise to achieve improved efficacy in medical applications as well as reduced systemic side effects compared to current standard of care medical practices. In particular, this article highlights the recent development and application of nanoparticles for diagnosing and treating nephropathologies.

MAGPI: A framework for maximum likelihood MR phase imaging using multiple receive coils

PURPOSE

Combining MR phase images from multiple receive coils is a challenging problem, complicated by ambiguities introduced by phase wrapping, noise, and the unknown phase-offset between the coils. Various techniques have been proposed to mitigate the effect of these ambiguities but most of the existing methods require additional reference scans and/or use ad hoc post-processing techniques that do not guarantee any optimality.

THEORY AND METHODS

Here, the phase estimation problem is formulated rigorously using a maximum-likelihood (ML) approach. The proposed framework jointly designs the acquisition-processing chain: the optimized pulse sequence is a single multiecho gradient echo scan and the corresponding postprocessing algorithm is a voxel-per-voxel ML estimator of the underlying tissue phase.

RESULTS

Our proposed framework (Maximum AmbiGuity distance for Phase Imaging, MAGPI) achieves substantial improvements in the phase estimate, resulting in phase signal-to-noise ratio (SNR) gains by up to an order of magnitude compared to existing methods.

CONCLUSION

The advantages of MAGPI are: (1) ML-optimal combination of phase data from multiple receive coils, without a reference scan; (2) voxel-per-voxel ML-optimal estimation of the underlying tissue phase, without the need for phase unwrapping or image smoothing; and (3) robust dynamic estimation of channel-dependent phase-offsets.

A Joint Acquisition-Estimation Framework for MR Phase Imaging

Measuring the phase of the MR signal is faced with fundamental challenges such as phase aliasing, noise and unknown offsets of the coil array. There is a paucity of acquisition, reconstruction and estimation methods that rigorously address these challenges. This reduces the reliability of information processing in phase domain. We propose a joint acquisition-processing framework that addresses the challenges of MR phase imaging using a rigorous theoretical treatment. Our proposed solution acquires the multi-coil complex data without any increase in acquisition time. Our corresponding estimation algorithm is applied optimally voxel-per-voxel. Results show that our framework achieves performance gains up to an order of magnitude compared to existing methods.

Shaped saturation with inherent radiofrequency-power-efficient trajectory design in parallel transmission

PURPOSE

A target-pattern-driven (TD) trajectory design is introduced in combination with parallel transmit (pTX) radiofrequency (RF) pulses to provide localized suppression of unwanted signals. The design incorporates target-pattern and B1+ information to adjust denser sampling and coverage in k-space regions where the main pattern information lies. Based on this approach, two-dimensional RF spiral saturation pulses sensitive to RF power limits were applied in vivo for the first time.

THEORY AND METHODS

The TD method was compared with two state-of-the-art spiral design methods. Simulations at different spatial fidelities, acceleration factors and anatomical regions were carried out for an eight-channel pTX 3 Tesla (T) coil. Human in vivo experiments were performed on a two-channel pTX 3T scanner saturating shaped patterns in the brain, heart, and thoracic spine.

RESULTS

Using the TD trajectory, RF pulse power can be substantially reduced by up to 34% compared with other trajectory designs with the same spatial accuracy. Local and global specific absorption rates are decreased in most cases.

CONCLUSION

The TD trajectory design uses available a priori information to enhance RF power efficiency and spatial response of the RF pulses. Shaped saturation pulses show improved spatial accuracy and saturation performance. Thus, RF pulses can be designed more efficiently and can be further accelerated.

A method for unwrapping highly wrapped multi-echo phase images at very high field: UMPIRE

PURPOSE

To develop a method of unwrapping phase images from multi-echo scans that works even where there are several wraps between echoes, and which generates unwrapped phase images in addition to phase difference (PD) images.

THEORY

The difference between the echo spacings in an acquisition with three unevenly spaced echoes (an imposed delay) can be selected such that the phase evolution in that time is in the range -π to +π in all voxels of interest. Under this condition, an image of the difference between the phase evolutions in the two inter-echo periods, an estimate of ΔB0 , is free of wraps. This ΔB0 estimate can be used to identify and remove receiver phase offsets and wraps in phase images.

METHODS

The approach was tested on simulated data and high-resolution in vivo brain data acquired from six subjects at 7 Tesla.

RESULTS

The method generated wrap-free phase images. It was able to remove more wraps than is possible with PD imaging and was faster and more reliable than spatial unwrapping.

CONCLUSION

Unwrapping Multi-echo Phase Images with iRregular Echo spacings (UMPIRE) is conceptually simple, fast, reliable, and requires no fitting, thresholds, or operator intervention.