Preliminary Experience with Three Alternative Motion Sensors for 0.55 Tesla MR Imaging.

Due to limitations in current motion tracking technologies and increasing interest in alternative sensors for motion tracking both inside and outside the MRI system, in this study we share our preliminary experience with three alternative sensors utilizing diverse technologies and interactions with tissue to monitor motion of the body surface, respiratory-related motion of major organs, and non-respiratory motion of deep-seated organs. These consist of (1) a Pilot-Tone RF transmitter combined with deep learning algorithms for tracking liver motion, (2) a single-channel ultrasound transducer with deep learning for monitoring bladder motion, and (3) a 3D Time-of-Flight camera for observing the motion of the anterior torso surface. Additionally, we demonstrate the capability of these sensors to simultaneously capture motion data outside the MRI environment, which is particularly relevant for procedures like radiation therapy, where motion status could be related to previously characterized cyclical anatomical data. Our findings indicate that the ultrasound sensor can track motion in deep-seated organs (bladder) as well as respiratory-related motion. The Time-of-Flight camera offers ease of interpretation and performs well in detecting surface motion (respiration). The Pilot-Tone demonstrates efficacy in tracking bulk respiratory motion and motion of major organs (liver). Simultaneous use of all three sensors could provide complementary motion information outside the MRI bore, providing potential value for motion tracking during position-sensitive treatments such as radiation therapy.

External Hardware and Sensors, for Improved MRI.

Complex engineered systems are often equipped with suites of sensors and ancillary devices that monitor their performance and maintenance needs. MRI scanners are no different in this regard. Some of the ancillary devices available to support MRI equipment, the ones of particular interest here, have the distinction of actually participating in the image acquisition process itself. Most commonly, such devices are used to monitor physiological motion or variations in the scanner's imaging fields, allowing the imaging and/or reconstruction process to adapt as imaging conditions change. "Classic" examples include electrocardiography (ECG) leads and respiratory bellows to monitor cardiac and respiratory motion, which have been standard equipment in scan rooms since the early days of MRI. Since then, many additional sensors and devices have been proposed to support MRI acquisitions. The main physical properties that they measure may be primarily "mechanical" (eg acceleration, speed, and torque), "acoustic" (sound and ultrasound), "optical" (light and infrared), or "electromagnetic" in nature. A review of these ancillary devices, as currently available in clinical and research settings, is presented here. In our opinion, these devices are not in competition with each other: as long as they provide useful and unique information, do not interfere with each other and are not prohibitively cumbersome to use, they might find their proper place in future suites of sensors. In time, MRI acquisitions will likely include a plurality of complementary signals. A little like the microbiome that provides genetic diversity to organisms, these devices can provide signal diversity to MRI acquisitions and enrich measurements. Machine-learning (ML) algorithms are well suited at combining diverse input signals toward coherent outputs, and they could make use of all such information toward improved MRI capabilities. EVIDENCE LEVEL: 2 TECHNICAL EFFICACY: Stage 1.

A relaxometry method that emphasizes practicality and availability.

PURPOSE: Although many methods have been proposed to quantitatively map the main MRI parameters (e.g., T1 , T2 , C × M0 ), these methods often involve special sequences not readily available on clinical scanners and/or may require long scan times. In contrast, the proposed method can readily run on most scanners, offer flexible tradeoffs between scan time and image quality, and map MRI parameters jointly to ensure spatial alignment. METHODS: The approach is based on the multi-shot spin-echo (SE) EPI sequence. The corresponding signal equation was derived and strategies for solving it were developed. As usual with multi-shot EPI, scan time can readily be traded-off against image quality by adjusting the echo train length. Validation was performed against reference relaxometry methods, in gel phantoms with varying concentrations of gadobutrol and gadoterate meglumine contrast agents. In vivo examples are further presented, from 3 neuroradiology patients. RESULTS: Bland-Altman analysis was performed: for T2 , as compared to 2D SE, bias was 0.29 ms and the 95% limits of agreement ranged from -1.15 to +1.73 ms. For T1 , compared to inversion-recovery SE (and MOLLI), bias was -20.2 ms (and -14.5 ms) and the limits of agreement ranged from -62.4 to +22.0 ms (and -53.8 to +24.9 ms). The mean relative T1 error between the proposed method and each of the 2 reference methods was similar to that of the reference methods among themselves. CONCLUSION: In the constellation of existing relaxometry methods, the proposed method is meant to stand out in terms of its practicality and availability.

Expanding Applications of Pulmonary MRI in the Clinical Evaluation of Lung Disorders: Fleischner Society Position Paper.

Pulmonary MRI provides structural and quantitative functional images of the lungs without ionizing radiation, but it has had limited clinical use due to low signal intensity from the lung parenchyma. The lack of radiation makes pulmonary MRI an ideal modality for pediatric examinations, pregnant women, and patients requiring serial and longitudinal follow-up. Fortunately, recent MRI techniques, including ultrashort echo time and zero echo time, are expanding clinical opportunities for pulmonary MRI. With the use of multicoil parallel acquisitions and acceleration methods, these techniques make pulmonary MRI practical for evaluating lung parenchymal and pulmonary vascular diseases. The purpose of this Fleischner Society position paper is to familiarize radiologists and other interested clinicians with these advances in pulmonary MRI and to stratify the Society recommendations for the clinical use of pulmonary MRI into three categories: (a) suggested for current clinical use, (b) promising but requiring further validation or regulatory approval, and (c) appropriate for research investigations. This position paper also provides recommendations for vendors and infrastructure, identifies methods for hypothesis-driven research, and suggests opportunities for prospective, randomized multicenter trials to investigate and validate lung MRI methods.

Multi-pathway multi-echo acquisition and neural contrast translation to generate a variety of quantitative and qualitative image contrasts.

PURPOSE: Clinical exams typically involve acquiring many different image contrasts to help discriminate healthy from diseased states. Ideally, 3D quantitative maps of all of the main MR parameters would be obtained for improved tissue characterization. Using data from a 7-min whole-brain multi-pathway multi-echo (MPME) scan, we aimed to synthesize several 3D quantitative maps (T1 and T2 ) and qualitative contrasts (MPRAGE, FLAIR, T1 -weighted, T2 -weighted, and proton density [PD]-weighted). The ability of MPME acquisitions to capture large amounts of information in a relatively short amount of time suggests it may help reduce the duration of neuro MR exams. METHODS: Eight healthy volunteers were imaged at 3.0T using a 3D isotropic (1.2 mm) MPME sequence. Spin-echo, MPRAGE, and FLAIR scans were performed for training and validation. MPME signals were interpreted through neural networks for predictions of different quantitative and qualitative contrasts. Predictions were compared to reference values at voxel and region-of-interest levels. RESULTS: Mean absolute errors (MAEs) for T1 and T2 maps were 216 ms and 11 ms, respectively. In ROIs containing white matter (WM) and thalamus tissues, the mean T1 /T2 predicted values were 899/62 ms and 1139/58 ms, consistent with reference values of 850/66 ms and 1126/58 ms, respectively. For qualitative contrasts, signals were normalized to those of WM, and MAEs for MPRAGE, FLAIR, T1 -weighted, T2 -weighted, and PD-weighted contrasts were 0.14, 0.15, 0.13, 0.16, and 0.05, respectively. CONCLUSIONS: Using an MPME sequence and neural-network contrast translation, whole-brain results were obtained with a variety of quantitative and qualitative contrast in ~6.8 min.

Improved PRF-based MR thermometry using k-space energy spectrum analysis.

PURPOSE: Proton resonance frequency (PRF) thermometry encodes information in the phase of MRI signals. A multiplicative factor converts phase changes into temperature changes, and this factor includes the TE. However, phase variations caused by B0 and/or B1 inhomogeneities can effectively change TE in ways that vary from pixel to pixel. This work presents how spatial phase variations affect temperature maps and how to correct for corresponding errors. METHODS: A method called "k-space energy spectrum analysis" was used to map regions in the object domain to regions in the k-space domain. Focused ultrasound heating experiments were performed in tissue-mimicking gel phantoms under two scenarios: with and without proper shimming. The second scenario, with deliberately de-adjusted shimming, was meant to emulate B0 inhomogeneities in a controlled manner. The TE errors were mapped and compensated for using k-space energy spectrum analysis, and corrected results were compared with reference results. Furthermore, a volunteer was recruited to help evaluate the magnitude of the errors being corrected. RESULTS: The in vivo abdominal results showed that the TE and heating errors being corrected can readily exceed 10%. In phantom results, a linear regression between reference and corrected temperatures results provided a slope of 0.971 and R2 of 0.9964. Analysis based on the Bland-Altman method provided a bias of -0.0977°C and 95% limits of agreement that were 0.75°C apart. CONCLUSION: Spatially varying TE errors, such as caused by B0 and/or B1 inhomogeneities, can be detected and corrected using the k-space energy spectrum analysis method, for increased accuracy in proton resonance frequency thermometry.

Relative Magnetic Force Measures and Their Potential Role in MRI Safety Practice.

BACKGROUND: Magnetic field markings are occasionally used at MRI sites to provide visual feedback of magnetic field strength at locations within the MRI scan room for safety purposes. In addition to magnetic field line markings, relative magnetic force, or ratio of magnetic to gravitational forces on an object, may be considered a useful complementary metric to quantify the risk associated with bringing objects containing ferromagnetic material into the magnetic field. PURPOSE: To develop and validate methods for deriving useful relative magnetic-force measures including a simple force index for application to MRI safety. STUDY TYPE: Phantom. PHANTOM: A special-purpose rig was built to experimentally measure relative magnetic forces on small ferromagnetic objects. FIELD STRENGTH: Ranging from 1.5T to 7T. ASSESSMENT: Quantitative comparisons were made between theoretical and measured relative magnetic forces on six objects containing ferromagnetic material: a piece of iron, a paper clip, a Kelly clamp, nail clippers, a cell phone, and a small permanent magnet. STATISTICAL TESTS: An analysis based on the Bland-Altman method was employed. RESULTS: After correction of the 1.5T data to account for assumed positioning errors of the test rig, limits of agreement between measured and estimated relative forces in the four MRI systems were ±0.16, where a relative force of 1.0 indicates that the magnetic force is equal to gravitation force. There was no significant bias in the data (P < = 0.05). DATA CONCLUSION: Accurate measures of relative magnetic forces on ferromagnetic objects can be derived for MRI safety purposes. LEVEL OF EVIDENCE: 1 Technical Efficacy Stage: 1 J. Magn. Reson. Imaging 2020;51:1260-1271.

Multipathway multi-echo (MPME) imaging: all main MR parameters mapped based on a single 3D scan.

PURPOSE: Quantitative parameter maps, as opposed to qualitative grayscale images, may represent the future of diagnostic MRI. A new quantitative MRI method is introduced here that requires a single 3D acquisition, allowing good spatial coverage to be achieved in relatively short scan times. METHODS: A multipathway multi-echo sequence was developed, and at least 3 pathways with 2 TEs were needed to generate T1 , T2 , T2* , B1+ , and B0 maps. The method required the central k-space region to be sampled twice, with the same sequence but with 2 very different nominal flip angle settings. Consequently, scan time was only slightly longer than that of a single scan. The multipathway multi-echo data were reconstructed into parameter maps, for phantom as well as brain acquisitions, in 5 healthy volunteers at 3 T. Spatial resolution, matrix size, and FOV were 1.2 × 1.0 × 1.2 mm3 , 160 × 192 × 160, and 19.2 × 19.2 × 19.2 cm3 (whole brain), acquired in 11.5 minutes with minimal acceleration. Validation was performed against T1 , T2 , and T2* maps calculated from gradient-echo and spin-echo data. RESULTS: In Bland-Altman plots, bias and limits of agreement for T1 and T2 results in vivo and in phantom were -2.9/±125.5 ms (T1 in vivo), -4.8/±20.8 ms (T2 in vivo), -1.5/±18.1 ms (T1 in phantom), and -5.3/±7.4 ms (T2 in phantom), for regions of interest including given brain structures or phantom compartments. Due to relatively high noise levels, the current implementation of the approach may prove more useful for region of interest-based as opposed to pixel-based interpretation. CONCLUSIONS: We proposed a novel approach to quantitatively map MR parameters based on a multipathway multi-echo acquisition.

The physics of MRI safety.

The main risks associated with magnetic resonance imaging (MRI) have been extensively reported and studied; for example, everyday objects may turn into projectiles, energy deposition can cause burns, varying fields can induce nerve stimulation, and loud noises can lead to auditory loss. The present review article is geared toward providing intuition about the physical mechanisms that give rise to these risks. On the one hand, excellent literature already exists on the practical aspect of risk management, with clinical workflow and recommendations. On the other hand, excellent technical articles also exist that explain these risks from basic principles of electromagnetism. We felt that an underserved niche might be found between the two, ie, somewhere between basic science and practical advice, to help develop intuition about electromagnetism that might prove of practical value when working around MR scanners. Following a wide-ranging introduction, risks originating from the main magnetic field, the excitation RF electromagnetic field, and switching of the imaging gradients will be presented in turn. LEVEL OF EVIDENCE: 5 Technical Efficacy: 1 J. Magn. Reson. Imaging 2018;47:28-43.

Application of a forward model of axisymmetric shear wave propagation in viscoelastic media to shear wave elastography.

A simple but general solution of Navier's equation for axisymmetric shear wave propagation in a homogeneous isotropic viscoelastic medium is presented. It is well-suited for use as a forward model for some acoustic radiation force impulse based shear wave elastography applications because it does not require precise knowledge of the strength of the source, nor its spatial or temporal distribution. Instead, it depends on two assumptions: (1) the source distribution is axisymmetric and confined to a small region near the axis of symmetry, and (2) the propagation medium is isotropic and homogeneous. The model accounts for the vector polarization of shear waves and exactly represents geometric spreading of the shear wavefield, whether spherical, cylindrical, or neither. It makes no assumption about the frequency dependence of material parameters, i.e., it is material-model independent. Validation using measured shear wavefields excited by acoustic radiation force in a homogeneous gelatin sample show that the model accounts for well over 90% of the measured wavefield "energy." An optimal fit of the model to simulated shear wavefields with noise in a homogeneous viscoelastic medium enables estimation of both the shear storage modulus and shear wave attenuation to within 1%.

Fast diffusion imaging with high angular resolution.

PURPOSE: High angular resolution diffusion imaging (HARDI) is a well-established method to help reveal the architecture of nerve bundles, but long scan times and geometric distortions inherent to echo planar imaging (EPI) have limited its integration into clinical protocols. METHODS: A fast imaging method is proposed here that combines accelerated multishot diffusion imaging (AMDI), multiplexed sensitivity encoding (MUSE), and crossing fiber angular resolution of intravoxel structure (CFARI) to reduce spatial distortions and reduce total scan time. A multishot EPI sequence was used to improve geometrical fidelity as compared to a single-shot EPI acquisition, and acceleration in both k-space and diffusion sampling enabled reductions in scan time. The method is regularized and self-navigated for motion correction. Seven volunteers were scanned in this study, including four with volumetric whole brain acquisitions. RESULTS: The average similarity of microstructural orientations between undersampled datasets and their fully sampled counterparts was above 85%, with scan times below 5 min for whole-brain acquisitions. Up to 2.7-fold scan time acceleration along with four-fold distortion reduction was achieved. CONCLUSION: The proposed imaging strategy can generate HARDI results with relatively good geometrical fidelity and low scan duration, which may help facilitate the transition of HARDI from a successful research tool to a practical clinical one. Magn Reson Med 77:696-706, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Dual-Pathway sequences for MR thermometry: When and where to use them.

PURPOSE: Dual-pathway sequences have been proposed to help improve the temperature-to-noise ratio (TNR) in MR thermometry. The present work establishes how much of an improvement these so-called "PSIF-FISP" sequences may bring in various organs and tissues. METHODS: Simulations and TNR calculations were validated against analytical equations, phantom, abdomen, and brain scans. Relative TNRs for PSIF-FISP, as compared to a dual-FISP reference standard, were calculated for flip angle (FA) = 1 to 85 º and repetition time (TR) = 6 to 60 ms, for gray matter, white matter, cervix, endometrium, myometrium, prostate, kidney medulla and cortex, bone marrow, pancreas, spleen, muscle, and liver tissues. RESULTS: PSIF-FISP was TNR superior in the kidney, pelvis, spleen, or gray matter at most tested TR and FA settings, and benefits increased at shorter TRs. PSIF-FISP was TNR superior in other tissues, e.g., liver, muscle, pancreas, for only short TR settings (20 ms or less). The TNR benefits of PSIF-FISP increased slightly with FA, and strongly with decreasing TR. Up to two- to three-fold reductions in TR with 20% TNR gains were achievable. In any given tissue, TNR performance is expected to further improve with heating, due to changes in relaxation rates. CONCLUSION: Dual-pathway PSIF-FISP can improve TNR and acquisition speed over standard gradient-recalled echo sequences, but optimal acquisition parameters are tissue dependent. Magn Reson Med 77:1193-1200, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Hybrid MRI-Ultrasound acquisitions, and scannerless real-time imaging.

PURPOSE: To combine MRI, ultrasound, and computer science methodologies toward generating MRI contrast at the high frame rates of ultrasound, inside and even outside the MRI bore. METHODS: A small transducer, held onto the abdomen with an adhesive bandage, collected ultrasound signals during MRI. Based on these ultrasound signals and their correlations with MRI, a machine-learning algorithm created synthetic MR images at frame rates up to 100 per second. In one particular implementation, volunteers were taken out of the MRI bore with the ultrasound sensor still in place, and MR images were generated on the basis of ultrasound signal and learned correlations alone in a "scannerless" manner. RESULTS: Hybrid ultrasound-MRI data were acquired in eight separate imaging sessions. Locations of liver features, in synthetic images, were compared with those from acquired images: The mean error was 1.0 pixel (2.1 mm), with best case 0.4 and worst case 4.1 pixels (in the presence of heavy coughing). For results from outside the bore, qualitative validation involved optically tracked ultrasound imaging with/without coughing. CONCLUSION: The proposed setup can generate an accurate stream of high-speed MR images, up to 100 frames per second, inside or even outside the MR bore. Magn Reson Med 78:897-908, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Accurate field mapping in the presence of B0 inhomogeneities, applied to MR thermometry.

PURPOSE: To describe how B0 inhomogeneities can cause errors in proton resonance frequency (PRF) shift thermometry, and to correct for these errors. METHODS: With PRF thermometry, measured phase shifts are converted into temperature measurements through the use of a scaling factor proportional to the echo time, TE. However, B0 inhomogeneities can deform, spread, and translate MR echoes, potentially making the "true" echo time vary spatially within the imaged object and take on values that differ from the prescribed TE value. Acquisition and reconstruction methods able to avoid or correct for such errors are presented. RESULTS: Tests were performed in a gel phantom during sonication, and temperature measurements were made with proper shimming as well as with intentionally introduced B0 inhomogeneities. Errors caused by B0 inhomogeneities were observed, described, and corrected by the proposed methods. No statistical difference was found between the corrected results and the reference results obtained with proper shimming, while errors by more than 10% in temperature elevation were corrected for. The approach was also applied to an abdominal in vivo dataset. CONCLUSION: Field variations induce errors in measured field values, which can be detected and corrected. The approach was validated for a PRF thermometry application.

Accelerated multi-shot diffusion imaging.

PURPOSE: To reduce image distortion in MR diffusion imaging using an accelerated multi-shot method. METHODS: The proposed method exploits the fact that diffusion-encoded data tend to be sparse when represented in the kb-kd space, where kb and kd are the Fourier transform duals of b and d, the b-factor and the diffusion direction, respectively. Aliasing artifacts are displaced toward under-used regions of the kb-kd plane, allowing nonaliased signals to be recovered. A main characteristic of the proposed approach is how thoroughly the navigator information gets used during reconstruction: The phase of navigator images is used for motion correction, while the magnitude of the navigator signal in kb-kd space is used for regularization purposes. As opposed to most acceleration methods based on compressed sensing, the proposed method reduces the number of ky lines needed for each diffusion-encoded image, but not the total number of images required. Consequently, it tends to be most effective at reducing image distortion rather than reducing total scan time. RESULTS: Results are presented for three volunteers with acceleration factors ranging from 4 to 8, with and without the inclusion of parallel imaging. CONCLUSION: An accelerated motion-corrected diffusion imaging method was introduced that achieves good image quality at relatively high acceleration factors.

Incidence and severity of pulmonary embolism in COVID-19 infection: Ancestral, Alpha, Delta, and Omicron variants.

Little information is available regarding incidence and severity of pulmonary embolism (PE) across the periods of ancestral strain, Alpha, Delta, and Omicron variants. The aim of this study is to investigate the incidence and severity of PE over the dominant periods of ancestral strain and Alpha, Delta, and Omicron variants. We hypothesized that the incidence and the severity by proximity of PE in patients with the newer variants and vaccination would be decreased compared with those in ancestral and earlier variants. Patients with COVID-19 diagnosis between March 2020 and February 2022 and computed tomography pulmonary angiogram performed within a 6-week window around the diagnosis (-2 to +4 weeks) were studied retrospectively. The primary endpoints were the associations of the incidence and location of PE with the ancestral strain and each variant. Of the 720 coronavirus disease 2019 patients with computed tomography pulmonary angiogram (58.6 ± 17.2 years; 374 females), PE was diagnosed among 42/358 (12%) during the ancestral strain period, 5/60 (8%) during the Alpha variant period, 16/152 (11%) during the Delta variant period, and 13/150 (9%) during the Omicron variant period. The most proximal PE (ancestral strain vs variants) was located in the main/lobar arteries (31% vs 6%-40%), in the segmental arteries (52% vs 60%-75%), and in the subsegmental arteries (17% vs 0%-19%). There was no significant difference in both the incidence and location of PE across the periods, confirmed by multivariable logistic regression models. In summary, the incidence and severity of PE did not significantly differ across the periods of ancestral strain and Alpha, Delta, and Omicron variants.

Severe COVID-19 pneumonia leads to post-COVID-19 lung abnormalities on follow-up CT scans.

Purpose: To investigate the association of the maximal severity of pneumonia on CT scans obtained within 6-week of diagnosis with the subsequent development of post-COVID-19 lung abnormalities (Co-LA). Methods: COVID-19 patients diagnosed at our hospital between March 2020 and September 2021 were studied retrospectively. The patients were included if they had (1) at least one chest CT scan available within 6-week of diagnosis; and (2) at least one follow-up chest CT scan available ≥ 6 months after diagnosis, which were evaluated by two independent radiologists. Pneumonia Severity Categories were assigned on CT at diagnosis according to the CT patterns of pneumonia and extent as: 1) no pneumonia (Estimated Extent, 0%); 2) non-extensive pneumonia (GGO and OP, <40%); and 3) extensive pneumonia (extensive OP and DAD, >40%). Co-LA on follow-up CT scans, categorized using a 3-point Co-LA Score (0, No Co-LA; 1, Indeterminate Co-LA; and 2, Co-LA). Results: Out of 132 patients, 42 patients (32%) developed Co-LA on their follow-up CT scans 6-24 months post diagnosis. The severity of COVID-19 pneumonia was associated with Co-LA: In 47 patients with extensive pneumonia, 33 patients (70%) developed Co-LA, of whom 18 (55%) developed fibrotic Co-LA. In 52 with non-extensive pneumonia, 9 (17%) developed Co-LA: In 33 with no pneumonia, none (0%) developed Co-LA. Conclusions: Higher severity of pneumonia at diagnosis was associated with the increased risk of development of Co-LA after 6-24 months of SARS-CoV-2 infection.

Investigation of the PSF-choice method for reduced lipid contamination in prostate MR spectroscopic imaging.

The purpose of this work was to evaluate a previously proposed approach that aims to improve the point spread function (PSF) of MR spectroscopic imaging (MRSI) to avoid corruption by lipid signal arising from neighboring voxels. Retrospective spatial filtering can be used to alter the PSF; however, this either reduces spatial resolution or requires extending the acquisition in k-space at the cost of increased imaging time. Alternatively, the method evaluated here, PSF-choice, can modify the PSF localization to reduce the contamination from adjacent lipids by conforming the signal response more closely to the desired MRSI voxel grid. This is done without increasing scan time or degrading SNR of important metabolites. PSF-choice achieves improvements in spatial localization through modifications to the radiofrequency excitation pulses. An implementation of this method is reported for MRSI of the prostate, where it is demonstrated that, in 13 of 16 pilot prostate MRSI scans, intravoxel spectral contamination from lipid was significantly reduced when using PSF-choice. Phantom studies were also performed that demonstrate, compared with MRSI with standard Fourier phase encoding, out-of-voxel signal contamination of spectra was significantly reduced in MRSI with PSF-choice.