Reconstruction of thermoacoustic emission sources induced by proton irradiation using numerical time reversal

Objective.Mapping of dose delivery in proton beam therapy can potentially be performed by analyzing thermoacoustic emissions measured by ultrasound arrays. Here, a method is derived and demonstrated for spatial mapping of thermoacoustic sources using numerical time reversal, simulating re-transmission of measured emissions into the medium.Approach.Spatial distributions of thermoacoustic emission sources are shown to be approximated by the analytic-signal form of the time-reversed acoustic field, evaluated at the time of the initial proton pulse. Given calibration of the array sensitivity and knowledge of tissue properties, this approach approximately reconstructs the acoustic source amplitude, equal to the product of the time derivative of the radiation dose rate, mass density, and Grüneisen parameter. This approach was implemented using two models for acoustic fields of the array elements, one modeling elements as line sources and the other as rectangular radiators. Thermoacoustic source reconstructions employed previously reported measurements of emissions from proton energy deposition in tissue-mimicking phantoms. For a phantom incorporating a bone layer, reconstructions accounted for the higher sound speed in bone. Dependence of reconstruction quality on array aperture size and signal-to-noise ratio was consistent with previous acoustic simulation studies.Main results.Thermoacoustic source distributions were successfully reconstructed from acoustic emissions measured by a linear ultrasound array. Spatial resolution of reconstructions was significantly improved in the azimuthal (array) direction by incorporation of array element diffraction. Source localization agreed well with Monte Carlo simulations of energy deposition, and was improved by incorporating effects of inhomogeneous sound speed.Significance.The presented numerical time reversal approach reconstructs thermoacoustic sources from proton beam radiation, based on straightforward processing of acoustic emissions measured by ultrasound arrays. This approach may be useful for ranging and dosimetry of clinical proton beams, if acoustic emissions of sufficient amplitude and bandwidth can be generated by therapeutic proton sources.

A Weighted Head Accelerator Mechanism (WHAM) for visualizing brain rheology using magnetic resonance imaging

BACKGROUND

A device for moving the head during MR imaging, called a Weighted Head Accelerator Mechanism (WHAM), rotates the head of a supine subject within programmable rotation limits and acceleration profiles. The WHAM can be used with custom MRI sequences to visualize the deformation and recoil of in vivo brain parenchyma with high temporal resolution, allowing element-wise calculation of strain and shear forces in the brain. Unlike previous devices, the WHAM can be configured to provide a wide range of motion and acceleration profiles.

NEW METHOD

The WHAM was calibrated using a high-speed camera on a laboratory bench and in 1.5 Tesla and 3.0 Tesla MRI scanners using gel phantoms and human subjects. The MR imaging studies employed a spatial spin-saturation tagging sub-sequence, followed by serial image acquisition. In these studies, 256 images were acquired with a temporal resolution of 2.56 ms. Deformation of the brain was quantified by following the spatial tags in the images.

RESULTS

MR imaging showed that the WHAM drove quantifiable brain motions using g forces less than those typically observed in day-to-day activities, with peak accelerations of ∼250 rad/sec².

COMPARISON WITH EXISTING METHODS

The peak pre-contact accelerations and velocities achieved by the WHAM device in this study are both higher than devices used in previous studies, while also allowing for modification of these factors.

CONCLUSIONS

MR imaging performed with the WHAM provides a direct method to visualize and quantify "brain slosh" in response to rotational acceleration. Consequently, this approach might find utility in evaluating strategies to protect the brain from mild traumatic brain injury (mTBI).

Intracardiac MR imaging (ICMRI) guiding-sheath with amplified expandable-tip imaging and MR-tracking for navigation and arrythmia ablation monitoring: Swine testing at 1.5 and 3T

PURPOSE

Develop a deflectable intracardiac MR imaging (ICMRI) guiding-sheath to accelerate imaging during MR-guided electrophysiological (EP) interventions for radiofrequency (500 kHz) ablation (RFA) of arrythmia. Requirements include imaging at three to five times surface-coil SNR in cardiac chambers, vascular insertion, steerable-active-navigation into cardiac chambers, operation with ablation catheters, and safe levels of MR-induced heating.

METHODS

ICMRI's 6 mm outer-diameter (OD) metallic-braided shaft had a 2.6 mm OD internal lumen for ablation-catheter insertion. Miniature-Baluns (MBaluns) on ICMRI's 1 m shaft reduced body-coil-induced heating. Distal section was a folded "star"-shaped imaging-coil mounted on an expandable frame, with an integrated miniature low-noise-amplifier overcoming cable losses. A handle-activated movable-shaft expanded imaging-coil to 35 mm OD for imaging within cardiac-chambers. Four MR-tracking micro-coils enabled navigation and motion-compensation, assuming a tetrahedron-shape when expanded. A second handle-lever enabled distal-tip deflection. ICMRI with a protruding deflectable EP catheter were used for MR-tracked navigation and RFA using a dedicated 3D-slicer user-interface. ICMRI was tested at 3T and 1.5T in swine to evaluate (a) heating, (b) cardiac-chamber access, (c) imaging field-of-view and SNR, and (d) intraprocedural RFA lesion monitoring.

RESULTS

The 3T and 1.5T imaging SNR demonstrated >400% SNR boost over a 4 × 4 × 4 cm3 FOV in the heart, relative to body and spine arrays. ICMRI with MBaluns met ASTM/IEC heating limits during navigation. Tip-deflection allowed navigating ICMRI and EP catheter into atria and ventricles. Acute-lesion long-inversion-time-T1-weighted 3D-imaging (TWILITE) ablation-monitoring using ICMRI required 5:30 min, half the time needed with surface arrays alone.

CONCLUSION

ICMRI assisted EP-catheter navigation to difficult targets and accelerated RFA monitoring.

Human upper-airway respiratory airflow: In vivo comparison of computational fluid dynamics simulations and hyperpolarized 129Xe phase contrast MRI velocimetry

Computational fluid dynamics (CFD) simulations of respiratory airflow have the potential to change the clinical assessment of regional airway function in health and disease, in pulmonary medicine and otolaryngology. For example, in diseases where multiple sites of airway obstruction occur, such as obstructive sleep apnea (OSA), CFD simulations can identify which sites of obstruction contribute most to airway resistance and may therefore be candidate sites for airway surgery. The main barrier to clinical uptake of respiratory CFD to date has been the difficulty in validating CFD results against a clinical gold standard. Invasive instrumentation of the upper airway to measure respiratory airflow velocity or pressure can disrupt the airflow and alter the subject's natural breathing patterns. Therefore, in this study, we instead propose phase contrast (PC) velocimetry magnetic resonance imaging (MRI) of inhaled hyperpolarized 129Xe gas as a non-invasive reference to which airflow velocities calculated via CFD can be compared. To that end, we performed subject-specific CFD simulations in airway models derived from 1H MRI, and using respiratory flowrate measurements acquired synchronously with MRI. Airflow velocity vectors calculated by CFD simulations were then qualitatively and quantitatively compared to velocity maps derived from PC velocimetry MRI of inhaled hyperpolarized 129Xe gas. The results show both techniques produce similar spatial distributions of high velocity regions in the anterior-posterior and foot-head directions, indicating good qualitative agreement. Statistically significant correlations and low Bland-Altman bias between the local velocity values produced by the two techniques indicates quantitative agreement. This preliminary in vivo comparison of respiratory airway CFD and PC MRI of hyperpolarized 129Xe gas demonstrates the feasibility of PC MRI as a technique to validate respiratory CFD and forms the basis for further comprehensive validation studies. This study is therefore a first step in the pathway towards clinical adoption of respiratory CFD.

MR Angiography Series: Fundamentals of Contrast-enhanced MR Angiography

The Society for Magnetic Resonance Angiography (SMRA) is a group of researchers and clinicians who are passionate about the benefits of MR angiography (MRA) but understand its challenges. Their mission is to study MRA, continually improve and innovate for the benefit of patients, and most important, educate the medical community so they can take full advantage of the benefits of MRA and overcome its challenges. In support of that mission, the authors have created a series of self-learning modules on behalf of the SMRA to demystify MRA protocols and help the reader perform patient-friendly high-quality MRA on a routine basis in clinical practice. The full digital presentation is available online. ^©RSNA, 2021.

Fusing acceleration and saturation techniques with wave amplitude labeling of time-shifted zeniths MR elastography

PURPOSE

To design a new 2D gradient recalled echo MR elastography (MRE) pulse sequence with inflow saturation for measuring liver stiffness in half the breath-hold time compared to standard of care (SC) 2D GRE MRE sequences.

METHODS

FASTWALTZ (fusing acceleration and saturation techniques with wave amplitude labeling of time-shifted zeniths) MRE employs an interleaved dual TR strategy with wave amplitude labeling and compressed SENSE undersampling to reduce breath-hold time while incorporating inflow saturation to suppress flow artifacts. The sequence was implemented and compared with SC MRE both in phantoms and in vivo in 5 asymptomatic volunteers. Stiffness values, region of interest size, and breath-hold times were compared between sequences.

RESULTS

Stiffness values were comparable between FASTWALTZ and SC MRE for both phantoms and in-vivo data. In volunteers, the group mean stiffness values at 60 Hz and region of interest size were 1.96 ± 0.30 kilopascals and 2279 ± 516 mm² for SC MRE, and 1.95 ± 0.29 kilopascals and 2061 ± 464 mm² for FASTWALTZ. Breath-hold duration for FASTWALTZ was 6.3 s compared to 13.3 s for SC MRE.

CONCLUSION

FASTWALTZ provides comparable stiffness values in half the breath-hold time compared to SC MRE and may have clinical benefits in patients with limited breath-holding capacity.

MRI Conditional Actively Tracked Metallic Electrophysiology Catheters and Guidewires With Miniature Tethered Radio-Frequency Traps: Theory, Design, and Validation

OBJECTIVE

Cardiovascular interventional devices typically have long metallic braids or backbones to aid in steerability and pushability. However, electromagnetic coupling of metallic-based cardiovascular interventional devices with the radiofrequency (RF) fields present during Magnetic Resonance Imaging (MRI) can make a device unsafe for use in an MRI scanner. We aimed to develop MRI conditional actively-tracked cardiovascular interventional devices by sufficiently attenuating induced currents on the metallic braid/tube and internal-cabling using miniaturized resonant floating RF traps (MBaluns).

METHOD

MBaluns were designed for placement at multiple locations along a conducting cardiovascular device to prevent the establishment of standing waves and to dissipate RF-induced energy. The MBaluns were constructed with loosely-wound solenoids to be sensitive to transverse magnetic fields created by both surface currents on the device's metallic backbone and common-mode currents on internal cables. Electromagnetic simulations were used to optimize MBalun parameters. Following optimization, two different MBalun designs were applied to MR-actively-tracked metallic guidewires and metallic-braided electrophysiology ablation catheters. Control-devices were constructed without MBaluns. MBalun performance was validated using network-analyzer quantification of current attenuation, electromagnetic Specific-Absorption-Rate (SAR) analysis, thermal tests during high SAR pulse sequences, and MRI-guided cardiovascular navigation in swine.

RESULTS

Electromagnetic SAR simulations resulted in ≈20 dB attenuation at the tip of the wire using six successive MBaluns. Network-analyzer tests confirmed ∼17 dB/MBalun surface-current attenuation. Thermal tests indicated temperature decreases of 5.9 °C in the MBalun-equipped guidewire tip. Both devices allowed rapid vascular navigation resulting from good torquability and MR-Tracking visibility.

CONCLUSION

MBaluns increased device diameter by 20%, relative to conventional devices, providing a spatially-efficient means to prevent heating during MRI.

SIGNIFICANCE

MBaluns allow use of long metallic components, which improves mechanical performance in active MR-guided interventional devices.

Assessing the relationship between movement and airflow in the upper airway using computational fluid dynamics with motion determined from magnetic resonance imaging

BACKGROUND

Computational fluid dynamics simulations of respiratory airflow in the upper airway reveal clinically relevant information, including sites of local resistance, inhaled particle deposition, and the effect of pathological constrictions. Unlike previous simulations, which have been performed on rigid anatomical models from static medical imaging, this work utilises ciné imaging during respiration to create dynamic models and more closely represent airway physiology.

METHODS

Airway movement maps were obtained from non-rigid image registration of fast-cine MRI and applied to high-spatial-resolution airway surface models. Breathing flowrates were recorded simultaneously with imaging. These data formed the boundary conditions for large eddy simulation computations of the airflow from exterior mask to bronchi. Simulations with rigid geometries were performed to demonstrate the resulting airflow differences between airflow simulations in rigid and dynamic airways.

FINDINGS

In the analysed rapid breathing manoeuvre, incorporating airway movement significantly changed the findings of the CFD simulations. Peak resistance increased by 19.8% and occurred earlier in the breath. Overall pressure loss decreased by 19.2%, and the proportion of flow in the mouth increased by 13.0%. Airway wall motion was out-of-phase with the air pressure force, demonstrating the presence of neuromuscular motion. In total, the anatomy did 25.2% more work on the air than vice versa.

INTERPRETATIONS

Realistic movement of the airway is incorporated into CFD simulations of airflow in the upper airway for the first time. This motion is vital to producing clinically relevant computational models of respiratory airflow and will allow novel analysis of dynamic conditions, such as sleep apnoea.

Mapping and correcting hyperpolarized magnetization decay with radial keyhole imaging

PURPOSE

Hyperpolarized (HP) media enable biomedical imaging applications that cannot be achieved with conventional MRI contrast agents. Unfortunately, quantifying HP images is challenging, because relaxation and radio-frequency pulsing generate spatially varying signal decay during acquisition. We demonstrate that, by combining center-out k-space sampling with postacquisition keyhole reconstruction, voxel-by-voxel maps of regional HP magnetization decay can be generated with no additional data collection.

THEORY AND METHODS

Digital phantom, HP ¹²⁹ Xe phantom, and in vivo ¹²⁹ Xe human (N = 4 healthy; N = 2 with cystic fibrosis) imaging was performed using radial sampling. Datasets were reconstructed using a postacquisition keyhole approach in which 2 temporally resolved images were created and used to generate maps of regional magnetization decay following a simple analytical model.

RESULTS

Mean, keyhole-derived decay terms showed excellent agreement with the decay used in simulations (R² = 0.996) and with global attenuation terms in HP ¹²⁹ Xe phantom imaging (R² > 0.97). Mean regional decay from in vivo imaging agreed well with global decay values and displayed spatial heterogeneity that matched expected variations in flip angle and oxygen partial pressure. Moreover, these maps could be used to correct variable signal decay across the image volume.

CONCLUSIONS

We have demonstrated that center-out trajectories combined with keyhole reconstruction can be used to map regional HP signal decay and to quantitatively correct images. This approach may be used to improve the accuracy of quantitative measures obtained from hyperpolarized media. Although validated with gaseous HP ¹²⁹ Xe in this work, this technique can be generalized to any hyperpolarized agent.

A novel method to generate dynamic boundary conditions for airway CFD by mapping upper airway movement with non-rigid registration of dynamic and static MRI

Computational fluid dynamics (CFD) simulations of airflow in the human airways have the potential to provide a great deal of information that can aid clinicians in case management and surgical decision making, such as airway resistance, energy expenditure, airflow distribution, heat and moisture transfer, and particle deposition, as well as the change in each of these due to surgical interventions. However, the clinical relevance of CFD simulations has been limited to date, as previous models either did not incorporate neuromuscular motion or any motion at all. Many common airway pathologies, such as obstructive sleep apnea (OSA) and tracheomalacia, involve large movements of the structures surrounding the airway, such as the tongue and soft palate. Airway wall motion may be due to many factors including neuromuscular motion, internal aerodynamic forces, and external forces such as gravity. Therefore, to realistically model these airway diseases, a method is required to derive the airway wall motion, whatever the cause, and apply it as a boundary condition to CFD simulations. This paper presents and validates a novel method of capturing in vivo motion of airway walls from magnetic resonance images with high spatiotemporal resolution, through a novel combination of non-rigid image, surface, and surface-normal-vector registration. Coupled with image-synchronous pneumotachography, this technique provides the necessary boundary conditions for dynamic CFD simulations of breathing, allowing the effect of the airway's complex motion to be calculated for the first time, in both normal subjects and those with conditions such as OSA.

Prospective Clinical Implementation of a Novel Magnetic Resonance Tracking Device for Real-Time Brachytherapy Catheter Positioning

PURPOSE

We designed and built dedicated active magnetic resonance (MR)-tracked (MRTR) stylets. We explored the role of MRTR in a prospective clinical trial.

METHODS AND MATERIALS

Eleven gynecologic cancer patients underwent MRTR to rapidly optimize interstitial catheter placement. MRTR catheter tip location and orientation were computed and overlaid on images displayed on in-room monitors at rates of 6 to 16 frames per second. Three modes of actively tracked navigation were analyzed: coarse navigation to the approximate region around the tumor; fine-tuning, bringing the stylets to the desired location; and pullback, with MRTR stylets rapidly withdrawn from within the catheters, providing catheter trajectories for radiation treatment planning (RTP). Catheters with conventional stylets were inserted, forming baseline locations. MRTR stylets were substituted, and catheter navigation was performed by a clinician working inside the MRI bore, using monitor feedback.

RESULTS

Coarse navigation allowed repositioning of the MRTR catheters tips by 16 mm (mean), relative to baseline, in 14 ± 5 s/catheter (mean ± standard deviation [SD]). The fine-tuning mode repositioned the catheter tips by a further 12 mm, in 24 ± 17 s/catheter. Pullback mode provided catheter trajectories with RTP point resolution of ∼1.5 mm, in 1 to 9 s/catheter.

CONCLUSIONS

MRTR-based navigation resulted in rapid and optimal placement of interstitial brachytherapy catheters. Catheters were repositioned compared with the initial insertion without tracking. In pullback mode, catheter trajectories matched computed tomographic precision, enabling their use for RTP.

Neonatal imaging using an on-site small footprint MR scanner

With its soft-tissue definition, multiplanar capabilities and advanced imaging techniques, magnetic resonance imaging (MRI) for neonatal care can provide better understanding of pathology, allowing for improved care and counseling to families. However, MR imaging in neonates is often difficult due to patient instability and the complex support necessary for survival. In our institution, we have installed a small footprint magnet in the neonatal intensive care unit (NICU) to minimize patient risks and provide the ability to perform MR imaging safely in this population. With this system, we have been able to provide more information with regard to central nervous system disorders, abdominal pathology, and pulmonary and airway abnormalities, and have performed postmortem imaging as an alternative or supplement to pathological autopsy. In our experience, an MR scanner situated within the NICU has allowed for safer and more expedited imaging of this vulnerable population.

Evaluation of an active magnetic resonance tracking system for interstitial brachytherapy

PURPOSE

In gynecologic cancers, magnetic resonance (MR) imaging is the modality of choice for visualizing tumors and their surroundings because of superior soft-tissue contrast. Real-time MR guidance of catheter placement in interstitial brachytherapy facilitates target coverage, and would be further improved by providing intraprocedural estimates of dosimetric coverage. A major obstacle to intraprocedural dosimetry is the time needed for catheter trajectory reconstruction. Herein the authors evaluate an active MR tracking (MRTR) system which provides rapid catheter tip localization and trajectory reconstruction. The authors assess the reliability and spatial accuracy of the MRTR system in comparison to standard catheter digitization using magnetic resonance imaging (MRI) and CT.

METHODS

The MRTR system includes a stylet with microcoils mounted on its shaft, which can be inserted into brachytherapy catheters and tracked by a dedicated MRTR sequence. Catheter tip localization errors of the MRTR system and their dependence on catheter locations and orientation inside the MR scanner were quantified with a water phantom. The distances between the tracked tip positions of the MRTR stylet and the predefined ground-truth tip positions were calculated for measurements performed at seven locations and with nine orientations. To evaluate catheter trajectory reconstruction, fifteen brachytherapy catheters were placed into a gel phantom with an embedded catheter fixation framework, with parallel or crossed paths. The MRTR stylet was then inserted sequentially into each catheter. During the removal of the MRTR stylet from within each catheter, a MRTR measurement was performed at 40 Hz to acquire the instantaneous stylet tip position, resulting in a series of three-dimensional (3D) positions along the catheter's trajectory. A 3D polynomial curve was fit to the tracked positions for each catheter, and equally spaced dwell points were then generated along the curve. High-resolution 3D MRI of the phantom was performed followed by catheter digitization based on the catheter's imaging artifacts. The catheter trajectory error was characterized in terms of the mean distance between corresponding dwell points in MRTR-generated catheter trajectory and MRI-based catheter digitization. The MRTR-based catheter trajectory reconstruction process was also performed on three gynecologic cancer patients, and then compared with catheter digitization based on MRI and CT.

RESULTS

The catheter tip localization error increased as the MRTR stylet moved further off-center and as the stylet's orientation deviated from the main magnetic field direction. Fifteen catheters' trajectories were reconstructed by MRTR. Compared with MRI-based digitization, the mean 3D error of MRTR-generated trajectories was 1.5 ± 0.5 mm with an in-plane error of 0.7 ± 0.2 mm and a tip error of 1.7 ± 0.5 mm. MRTR resolved ambiguity in catheter assignment due to crossed catheter paths, which is a common problem in image-based catheter digitization. In the patient studies, the MRTR-generated catheter trajectory was consistent with digitization based on both MRI and CT.

CONCLUSIONS

The MRTR system provides accurate catheter tip localization and trajectory reconstruction in the MR environment. Relative to the image-based methods, it improves the speed, safety, and reliability of the catheter trajectory reconstruction in interstitial brachytherapy. MRTR may enable in-procedural dosimetric evaluation of implant target coverage.

A 1.5T MRI-conditional 12-lead electrocardiogram for MRI and intra-MR intervention

PURPOSE

High-fidelity 12-lead electrocardiogram (ECG) is important for physiological monitoring of patients during MR-guided intervention and cardiac MRI. Issues in obtaining noncorrupted ECGs inside MRI include a superimposed magneto-hydro-dynamic voltage, gradient switching-induced voltages, and radiofrequency heating. These problems increase with magnetic field. The aim of this study is to develop and clinically validate a 1.5T MRI-conditional 12-lead ECG system.

METHODS

The system was constructed with transmission lines to reduce radiofrequency induction and switching circuits to remove induced voltages. Adaptive filters, trained by 12-lead measurements outside MRI and in two orientations inside MRI, were used to remove the magneto-hydro-dynamic voltage. The system was tested on 10 (one exercising) volunteers and four arrhythmia patients.

RESULTS

Switching circuits removed most imaging-induced voltages (residual noise <3% of the R-wave). Magneto-hydro-dynamic voltage removal provided intra-MRI ECGs that varied by <3.8% from those outside the MRI, preserving the true S-wave to T-wave segment. In premature ventricular contraction (PVC) patients, clean ECGs separated premature ventricular contraction and sinus rhythm beats. Measured heating was <1.5°C. The system reliably acquired multiphase (steady-state free precession) wall-motion-cine and phase-contrast-cine scans, including subjects in whom 4-lead gating failed. The system required a minimum repetition time of 4 ms to allow robust ECG processing.

CONCLUSION

High-fidelity intra-MRI 12-lead ECG is possible.

Voltage-based device tracking in a 1.5 Tesla MRI during imaging: initial validation in swine models

PURPOSE

Voltage-based device-tracking (VDT) systems are commonly used for tracking invasive devices in electrophysiological cardiac-arrhythmia therapy. During electrophysiological procedures, electro-anatomic mapping workstations provide guidance by integrating VDT location and intracardiac electrocardiogram information with X-ray, computerized tomography, ultrasound, and MR images. MR assists navigation, mapping, and radiofrequency ablation. Multimodality interventions require multiple patient transfers between an MRI and the X-ray/ultrasound electrophysiological suite, increasing the likelihood of patient-motion and image misregistration. An MRI-compatible VDT system may increase efficiency, as there is currently no single method to track devices both inside and outside the MRI scanner.

METHODS

An MRI-compatible VDT system was constructed by modifying a commercial system. Hardware was added to reduce MRI gradient-ramp and radiofrequency unblanking pulse interference. VDT patches and cables were modified to reduce heating. Five swine cardiac VDT electro-anatomic mapping interventions were performed, navigating inside and thereafter outside the MRI.

RESULTS

Three-catheter VDT interventions were performed at >12 frames per second both inside and outside the MRI scanner with <3 mm error. Catheters were followed on VDT- and MRI-derived maps. Simultaneous VDT and imaging was possible in repetition time >32 ms sequences with <0.5 mm errors, and <5% MRI signal-to-noise ratio (SNR) loss. At shorter repetition times, only intracardiac electrocardiogram was reliable. Radiofrequency heating was <1.5°C.

CONCLUSION

An MRI-compatible VDT system is feasible.

A novel acoustically quiet coil for neonatal MRI system

MRI acoustic exposure has the potential to elicit physiological distress and impact development in preterm and term infants. To mitigate this risk, a novel acoustically quiet coil was developed to reduce the sound pressure level experienced by neonates during MR procedures. The new coil has a conventional high-pass birdcage RF design, but is built on a framework of sound abating material. We evaluated the acoustic and MR imaging performance of the quiet coil and a conventional body coil on two small footprint NICU MRI systems. Sound pressure level and frequency response measurements were made for six standard clinical MR imaging protocols. The average sound pressure level, reported for all six imaging pulse sequences, was 82.2 dBA for the acoustically quiet coil, and 91.1 dBA for the conventional body coil. The sound pressure level values measured for the acoustically quiet coil were consistently lower, 9 dBA (range 6-10 dBA) quieter on average. The acoustic frequency response of the two coils showed a similar harmonic profile for all imaging sequences. However, the amplitude was lower for the quiet coil, by as much as 20 dBA.

Gradient-induced voltages on 12-lead ECGs during high duty-cycle MRI sequences and a method for their removal considering linear and concomitant gradient terms

PURPOSE

To restore 12-lead electrocardiographic (ECG) signal fidelity inside MRI by removing magnetic field gradient-induced voltages during high gradient duty cycle sequences.

THEORY AND METHODS

A theoretical equation was derived to provide first- and second-order electrical fields induced at individual ECG electrodes as a function of gradient fields. Experiments were performed at 3T on healthy volunteers using a customized acquisition system that captured the full amplitude and frequency response of ECGs, or a commercial recording system. The 19 equation coefficients were derived via linear regression of data from accelerated sequences and were used to compute induced voltages in real-time during full resolution sequences to remove ECG artifacts. Restored traces were evaluated relative to ones acquired without imaging.

RESULTS

Measured induced voltages were 0.7 V peak-to-peak during balanced steady state free precession (bSSFP) with the heart at the isocenter. Applying the equation during gradient echo sequencing, three-dimensional fast spin echo, and multislice bSSFP imaging restored nonsaturated traces and second-order concomitant terms showed larger contributions in electrodes further from the magnet isocenter. Equation coefficients are evaluated with high repeatability (ρ = 0.996) and are dependent on subject, sequence, and slice orientation.

CONCLUSION

Close agreement between theoretical and measured gradient-induced voltages allowed for real-time removal. Prospective estimation of sequence periods in which large induced voltages occur may allow hardware removal of these signals.

Characterization of acoustic noise in a neonatal intensive care unit MRI system

BACKGROUND

To eliminate the medical risks and logistical challenges of transporting infants from the neonatal intensive care unit (NICU) to the radiology department for magnetic resonance imaging, a small-footprint 1.5-T MRI scanner has been developed for neonatal imaging within the NICU. MRI is known to be noisy, and exposure to excessive acoustic noise has the potential to elicit physiological distress and impact development in the term and preterm infant.

OBJECTIVE

To measure and compare the acoustic noise properties of the NICU MRI system against those of a conventional 1.5-T MRI system.

MATERIALS AND METHODS

We performed sound pressure level measurements in the NICU MRI scanner and in a conventional adult-size whole-body 1.5-T MRI system. Sound pressure level measurements were made for six standard clinical MR imaging protocols.

RESULTS

The average sound pressure level value, reported in unweighted (dB) and A-weighted (dBA) decibels for all six imaging pulse sequences, was 73.8 dB and 88 dBA for the NICU scanner, and 87 dB and 98.4 dBA for the conventional MRI scanner. The sound pressure level values measured on the NICU scanner for each of the six MR imaging pulse sequences were consistently and significantly (P = 0.03) lower, with an average difference of 14.2 dB (range 10-21 dB) and 11 dBA (range 5-18 dBA). The sound pressure level frequency response of the two MR systems showed a similar harmonic structure above 200 Hz for all imaging sequences. The amplitude, however, was appreciably lower for the NICU scanner, by as much as 30 dB, for frequencies below 200 Hz.

CONCLUSION

The NICU MRI system is quieter than conventional MRI scanners, improving safety for the neonate and facilitating siting of the unit within the NICU.

Real-time active MR-tracking of metallic stylets in MR-guided radiation therapy

PURPOSE

To develop an active MR-tracking system to guide placement of metallic devices for radiation therapy.

METHODS

An actively tracked metallic stylet for brachytherapy was constructed by adding printed-circuit micro-coils to a commercial stylet. The coil design was optimized by electromagnetic simulation, and has a radio-frequency lobe pattern extending ∼5 mm beyond the strong B0 inhomogeneity region near the metal surface. An MR-tracking sequence with phase-field dithering was used to overcome residual effects of B0 and B1 inhomogeneities caused by the metal, as well as from inductive coupling to surrounding metallic stylets. The tracking system was integrated with a graphical workstation for real-time visualization. The 3 Tesla MRI catheter-insertion procedures were tested in phantoms and ex vivo animal tissue, and then performed in three patients during interstitial brachytherapy.

RESULTS

The tracking system provided high-resolution (0.6 × 0.6 × 0.6 mm(3) ) and rapid (16 to 40 frames per second, with three to one phase-field dithering directions) catheter localization in phantoms, animals, and three gynecologic cancer patients.

CONCLUSION

This is the first demonstration of active tracking of the shaft of metallic stylet in MR-guided brachytherapy. It holds the promise of assisting physicians to achieve better targeting and improving outcomes in interstitial brachytherapy.

Prospective motion correction using tracking coils

Intracavity imaging coils provide higher signal-to-noise than surface coils and have the potential to provide higher spatial resolution in shorter acquisition times. However, images from these coils suffer from physiologically induced motion artifacts, as both the anatomy and the coils move during image acquisition. We developed prospective motion-correction techniques for intracavity imaging using an array of tracking coils. The system had <50 ms latency between tracking and imaging, so that the images from the intracavity coil were acquired in a frame of reference defined by the tracking array rather than by the system's gradient coils. Two-dimensional gradient-recalled and three-dimensional electrocardiogram-gated inversion-recovery-fast-gradient-echo sequences were tested with prospective motion correction using ex vivo hearts placed on a moving platform simulating both respiratory and cardiac motion. Human abdominal tests were subsequently conducted. The tracking array provided a positional accuracy of 0.7 ± 0.5 mm, 0.6 ± 0.4 mm, and 0.1 ± 0.1 mm along the X, Y, and Z directions at a rate of 20 frames-per-second. The ex vivo and human experiments showed significant image quality improvements for both in-plane and through-plane motion correction, which although not performed in intracavity imaging, demonstrates the feasibility of implementing such a motion-correction system in a future design of combined tracking and intracavity coil.

An MRI system for imaging neonates in the NICU: initial feasibility study

BACKGROUND

Transporting premature infants from a neonatal intensive care unit (NICU) to a radiology department for MRI has medical risks and logistical challenges.

OBJECTIVE

To develop a small 1.5-T MRI system for neonatal imaging that can be easily installed in the NICU and to evaluate its performance using a sheep model of human prematurity.

MATERIALS AND METHODS

A 1.5-T MRI system designed for orthopedic use was adapted for neonatal imaging. The system was used for MRI examinations of the brain, chest and abdomen in 12 premature lambs during the first hours of life. Spin-echo, fast spin-echo and gradient-echo MR images were evaluated by two pediatric radiologists.

RESULTS

All animals remained physiologically stable throughout the imaging sessions. Animals were imaged at two or three time points. Seven brain MRI examinations were performed in seven different animals, 23 chest examinations in 12 animals and 19 abdominal examinations in 11 animals. At each anatomical location, high-quality images demonstrating good spatial resolution, signal-to-noise ratio and tissue contrast were routinely obtained within 30 min using standard clinical protocols.

CONCLUSION

Our preliminary experience demonstrates the feasibility and potential of the neonatal MRI system to provide state-of-the-art MRI capabilities within the NICU. Advantages include overall reduced cost and site demands, lower acoustic noise, improved ease of access and reduced medical risk to the neonate.

Phase-field dithering for active catheter tracking

A strategy to increase the robustness of active MR tracking of microcoils in low signal-to-noise ratio conditions was developed and tested. The method employs dephasing magnetic field gradient pulses that are applied orthogonal to the frequency-encoding gradient pulse used in conventional point-source MR tracking. In subsequent acquisitions, the orthogonal dephasing gradient pulse is rotated while maintaining a perpendicular orientation with respect to the frequency-encoding gradient. The effect of the dephasing gradient is to apply a spatially dependent phase shift in directions perpendicular to the frequency-encoding gradient. Since the desired MR signal for robust MR tracking comes from the small volume of nuclear spins near the small detection coil, the desired signal is not dramatically altered by the dephasing gradient. Undesired MR signals arising from larger volumes (e.g., due to coupling with the body coil or surface coils), on the other hand, are dephased and reduced in signal intensity. Since the approach requires no a priori knowledge of the microcoil orientation with respect to the main magnetic field, data from several orthogonal dephasing gradients are acquired and analyzed in real time. One of several selection algorithms is employed to identify the "best" data for use in the coil localization algorithm. This approach was tested in flow phantoms and animal models, with several multiplexing schemes, including the Hadamard and zero-phase reference approaches. It was found to provide improved MR tracking of untuned microcoils. It also dramatically improved MR tracking robustness in low signal-to-noise-ratio conditions and permitted tracking of microcoils that were inductively coupled to the body coil.

Electroanatomic mapping and radiofrequency ablation of porcine left atria and atrioventricular nodes using magnetic resonance catheter tracking

BACKGROUND

The MRI-compatible electrophysiology system previously used for MR-guided left ventricular electroanatomic mapping was enhanced with improved MR tracking, an MR-compatible radiofrequency ablation system and higher-resolution imaging sequences to enable mapping, ablation, and ablation monitoring in smaller cardiac structures. MR-tracked navigation was performed to the left atrium (LA) and atrioventricular (AV) node, followed by LA electroanatomic mapping and radiofrequency ablation of the pulmonary veins (PVs) and AV node.

METHODS AND RESULTS

One ventricular ablation, 7 PV ablations, 3 LA mappings, and 3 AV node ablations were conducted. Three MRI-compatible devices (ablation/mapping catheter, torqueable sheath, stimulation/pacing catheter) were used, each with 4 to 5 tracking microcoils. Transseptal puncture was performed under x-ray, with all other procedural steps performed in the MRI. Preacquired MRI roadmaps served for real-time catheter navigation. Simultaneous tracking of 3 devices was performed at 13 frames per second. LA mapping and PV radiofrequency ablation were performed using tracked ablation catheters and sheaths. Ablation points were registered and verified after ablation using 3D myocardial delayed enhancement and postmortem gross tissue examination. Complete LA electroanatomic mapping was achieved in 3 of 3 pigs, Right inferior PV circumferential ablation was achieved in 3 of 7 pigs, with incomplete isolation caused by limited catheter deflection. During AV node ablation, ventricular pacing was performed, 3 devices were simultaneously tracked, and intracardiac ECGs were displayed. 3D myocardial delayed enhancement visualized node injury 2 minutes after ablation. AV node block succeeded in 2 of 3 pigs, with 1 temporary block.

CONCLUSIONS

LA mapping, PV radiofrequency ablation, and AV node ablation were demonstrated under MRI guidance. Intraprocedural 3D myocardial delayed enhancement assessed lesion positional accuracy and dimensions.

Electroanatomic mapping of the left ventricle in a porcine model of chronic myocardial infarction with magnetic resonance-based catheter tracking

BACKGROUND

X-ray fluoroscopy constitutes the fundamental imaging modality for catheter visualization during interventional electrophysiology procedures. The minimal tissue discriminative capability of fluoroscopy is mitigated in part by the use of electroanatomic mapping systems and enhanced by the integration of preacquired 3-dimensional imaging of the heart with computed tomographic or magnetic resonance (MR) imaging. A more ideal paradigm might be to use intraprocedural MR imaging to directly image and guide catheter mapping procedures.

METHODS AND RESULTS

An MR imaging-based electroanatomic mapping system was designed to assess the feasibility of navigating catheters to the left ventricle in vivo using MR tracking of microcoils incorporated into the catheters, measuring intracardiac ventricular electrograms, and integrating this information with 3-dimensional MR angiography and myocardial delayed enhancement images to allow ventricular substrate mapping. In all animals (4 normal, and 10 chronically infarcted swine), after transseptal puncture under fluoroscopic guidance, catheters were successfully navigated to the left ventricle with MR tracking (13 to 15 frames per second) by both transseptal and retrograde aortic approaches. Electrogram artifacts related to the MR imaging gradient pulses were successfully removed with analog and digital signal processing. In all animals, it was possible to map the entire left ventricle and to project electrogram voltage amplitude maps to identify the scarred myocardium.

CONCLUSIONS

It is possible to use MR tracking to navigate catheters to the left ventricle, to measure electrogram activity, and to render accurate 3-dimensional voltage maps in a porcine model of chronic myocardial infarction, completely in the MR imaging environment. Myocardial delayed enhancement guidance provided dense sampling of the proximity of the infarct and accurate localization of complex infarcts.

Large field-of-view real-time MRI with a 32-channel system

The emergence of parallel MRI techniques and new applications for real-time interactive MRI underscores the need to evaluate performance gained by increasing the capability of MRI phased-array systems beyond the standard four to eight high-bandwidth channels. Therefore, to explore the advantages of highly parallel MRI a 32-channel 1.5 T MRI system and 32-element torso phased arrays were designed and constructed for real-time interactive MRI. The system was assembled from multiple synchronized scanner-receiver subsystems. Software was developed to coordinate across subsystems the real-time acquisition, reconstruction, and display of 32-channel images. Real-time, large field-of-view (FOV) body-survey imaging was performed using interleaved echo-planar and single-shot fast-spin-echo pulse sequences. A new method is demonstrated for augmenting parallel image acquisition by independently offsetting the frequency of different array elements (FASSET) to variably shift their FOV. When combined with conventional parallel imaging techniques, image acceleration factors of up to 4 were investigated. The use of a large number of coils allowed the FOV to be doubled in two dimensions during rapid imaging, with no degradation of imaging time or spatial resolution. The system provides a platform for evaluating the applications of many-channel real-time MRI, and for understanding the factors that optimize the choice of array size.

Feasibility of Stent Placement in Carotid Arteries with Real-time MR Imaging Guidance in Pigs

All examinations were performed with approval from the institutional animal care and use committee of Columbia University. To assess the feasibility of real-time magnetic resonance (MR) imaging-guided neurovascular intervention in a swine model, the authors placed stents in the carotid arteries of five domestic pigs. Seven-French vascular sheaths were placed in the target carotid arteries via femoral access by using active MR tracking. Ten nitinol stents (8-10 x 20-40 mm) were successfully deployed in the target segments of carotid arteries bilaterally. MR imaging and necropsy findings confirmed stent position. Necropsy revealed no gross vascular injury. Study results demonstrated the feasibility of performing real-time MR imaging-guided neurovascular intervention by using an active-tracking technique in an animal model.

Transfemoral catheterization of carotid arteries with real-time MR imaging guidance in pigs

All procedures and protocols were approved by the institutional animal care and use committee of Columbia University. To determine whether transfemoral catheterization of the carotid arteries can be performed entirely with real-time magnetic resonance (MR) imaging guidance, the authors catheterized the carotid arteries in six domestic pigs by using active-tracking catheters and guidewires and MR tracking software created for neurovascular procedures. The carotid arteries were successfully catheterized 24 times, on average within 5 minutes after insertion of the catheter into the femoral artery. Results demonstrated the feasibility of performing transfemoral catheterization of the carotid arteries with active MR tracking devices in a conventional MR imaging unit.

Highly parallel volumetric imaging with a 32-element RF coil array

The improvement of MRI speed with parallel acquisition is ultimately an SNR-limited process. To offset acquisition- and reconstruction-related SNR losses, practical parallel imaging at high accelerations should include the use of a many-element array with a high intrinsic signal-to-noise ratio (SNR) and spatial-encoding capability, and an advantageous imaging paradigm. We present a 32-element receive-coil array and a volumetric paradigm that address the SNR challenge at high accelerations by maximally exploiting multidimensional acceleration in conjunction with noise averaging. Geometric details beyond an initial design concept for the array were determined with the guidance of simulations. Imaging with the support of 32-channel data acquisition systems produced in vivo results with up to 16-fold acceleration, including images from rapid abdominal and MRA studies.

Real-time magnetic resonance imaging for interventional neuroradiological procedures

Over the last two decades, interventional neuroradiologists have developed powerful techniques for the treatment of cerebrovascular disorders and brain tumors. Current interventional neuroradiological procedures are performed under X-ray fluoroscopy, which has allowed for high temporal and spatial resolution. However, these imaging techniques do not provide the treating physician with vital anatomic and functional information regarding vessel walls and the surrounding brain tissue. Better visualization of vessel structures and real-time information about the state of perfusion and metabolism of the surrounding brain tissue (real-time magnetic resonance arteriography, diffusion and perfusion-weighted imaging, apparent diffusion coefficient maps) would enhance safety and efficacy of neuroendovascular procedures available currently. Recent advances in magnetic resonance hardware and software have permitted significant enhancements in temporal and spatial resolution, which have resulted in the capability of visualizing anatomic structures with real-time fluoroscopy and angiography. This review outlines how real-time magnetic resonance procedures may replace conventional X-ray fluoroscopy in diagnostic and interventional neuroradiology during the next decade.

Extended field-of-view imaging with table translation and frequency sweeping

A new method for MRI of an extended field of view (FOV) has been developed and validated. The method employs concurrent MR data acquisition and patient table motion. Table motion-induced image artifacts are minimized by sweeping the frequency of the receiver at a rate matching the table's speed. Multiple regional images are collected and combined to reconstruct the full FOV. The imaging parameters and table speed are chosen to ensure that each regional image of the subject is collected while the corresponding anatomy is in the useable imaging volume of the scanner. Additional strategies are applied to further reduce field inhomogeneity-induced artifacts, especially distortions due to gradient field nonlinearity. The method is robust and can be easily incorporated into most multislice 2D and volumetric 3D imaging pulse sequences. It is anticipated that this technique will be useful for a variety of applications, including angiographic runoffs, whole-body screening, and short-magnet imaging.

The influence of out-of-plane geometry on the flow within a distal end-to-side anastomosis

This paper describes a computational and experimental investigation of flow in a proto-type model geometry of a fully occluded 45 deg distal end-to-side anastomosis. Previous investigations have considered a similar configuration where the centerlines of the bypass and host vessels lie within a plane, thereby producing a plane of symmetry within the flow. We have extended these investigations by deforming the bypass vessel out of the plane of symmetry, thereby breaking the symmetry of the flow and producing a nonplanar geometry. Experimental data were obtained using magnetic resonance imaging of flow within perspex models and computational data were obtained from simulations using a high-order spectral/hp element method. We found that the nonplanar three-dimensional flow notably alters the distribution of wall shear stress at the bed of the anastomosis, reducing the peak wall shear stress peak by approximately 10 percent when compared with the planar model. Furthermore, an increase in the absolute flux of velocity into the occluded region, proximal to the anastomosis, of 80 percent was observed in the nonplanar geometry when compared with the planar geometry.

Glenohumeral relationships during physiologic shoulder motion and stress testing: initial experience with open MR imaging and active imaging-plane registration

PURPOSE

To test the hypotheses that open dynamic magnetic resonance (MR) imaging can (a) be used to evaluate and define normal shoulder motion in active joint motion and muscle contraction and (b) be used in conjunction with physical examination.

MATERIALS AND METHODS

With an open-configuration, 0.5-T MR imaging system and active image-plane tracking, 10 shoulders were studied in five asymptomatic subjects to establish normal patterns of glenohumeral motion during abduction and adduction and internal and external rotation. Preliminary studies of physical examination during MR imaging, in which a physician examiner applied mechanical force to the humeral head, were also performed.

RESULTS

During abduction and adduction and internal and external rotation maneuvers with active subjects muscle contraction, the humeral head remained precisely centered on the glenoid fossa in all asymptomatic subjects, which is in agreement with findings of previous radiographic studies. Application of force to the humeral head by an examiner was associated with as much as 6 mm of anterior translation and 13 mm of posterior translation.

CONCLUSION

Dynamic MR imaging of the glenohumeral joint is possible over a wide range of physiologic motion in vertically open systems. Use of an MR tracking coil enabled accurate tracking of the anatomy of interest. These preliminary measurements of normal glenohumeral motion patterns begin to establish normal ranges of motion and constitute a necessary first step in characterizing pathologic motion in patients with common clinical problems such as instability and impingement.

Joint motion in an open MR unit using MR tracking

A system for active scan plane guidance during kinematic magnetic resonance (MR) examination of joint motion was developed utilizing an external tracking coil and MR tracking software. In a phantom study and during upright, weight-bearing, physiologic knee flexion, the external tracking coil maintained the scan plane through desired structures. Thus, MR tracking provides a robust method to guide the scan plane during MR imaging of active joint motion.

Active MR tracking on a 0.2 Tesla MR imager

An active MR tracking system was implemented on a 0.2 Tesla open MRI system. Interventional devices with receive-only microcoils at their tips were developed and investigated on the scanner. Microcoils having a diameter of about 1 mm and 20 turns were found to provide sufficient signal-to-noise ratios for stable tracking. Positional accuracy and precision were found to be acceptable under practical conditions. Simulation of MR-guided biopsy using biplane images with tracking was performed in a gelatin phantom and dog livers. Successful tracking of catheters with integrated microcoils was also demonstrated in the aorta and IVC of live dogs.

MR-guided percutaneous angioplasty: assessment of tracking safety, catheter handling and functionality

PURPOSE

Magnetic resonance (MR)-guided percutaneous vascular interventions have evolved to a practical possibility with the advent of open-configuration MR systems and real-time tracking techniques. The purpose of this study was to assess an MR-tracking percutaneous transluminal angioplasty (PTA) catheter with regard to its safety profile and functionality.

METHODS

Real-time, biplanar tracking of the PTA catheter was made possible by incorporating a small radiofrequency (RF) coil in the catheter tip and connecting it to a coaxial cable embedded in the catheter wall. To evaluate potentially hazardous thermal effects due to the incorporation of the coil, temperature measurements were performed within and around the coil under various scanning and tracking conditions at 1.5 Tesla (T). Catheter force transmission and balloon-burst pressure of the MR-tracking PTA catheter were compared with those of a standard PTA catheter. The dilatative capability of the angioplasty balloon was assessed in vitro as well as in vivo, in an isolated femoral artery segment in a swine.

RESULTS

The degree of heating at the RF coil was directly proportional to the power of the RF pulses. Heating was negligible with MR tracking, conventional spin-echo and low-flip gradient-echo sequences. Sequences with higher duty cycles, such as fast spin echo, produced harmful heating effects. Force transmission of the MR-tracking PTA catheter was slightly inferior to that of the standard PTA catheter, while balloon-burst pressures were similar to those of conventional catheters. The MR-tracking PTA catheter functioned well both in vitro and in vivo.

CONCLUSION

The in vivo use of an MR-tracking PTA catheter is safe under most scanning conditions.

Interactive coronary MRI

The acquisition of complete three-dimensional (3D), segmented gradient-echo data sets to visualize the coronary arteries can be both time consuming and sensitive to motion, even with use of multiple breath-holding or respiratory gating. An alternate hybrid approach is demonstrated here, in which real-time interactive imaging is first used to locate an optimal oblique coronary scan plane. Then, a limited number of contiguous slices are acquired around that plane within a breath-hold with use of two-dimensional (2D) segmented gradient-echo imaging. Dual inversion nulling is used to suppress fat and myocardium. Finally, if needed, a limited reformat of the data is performed to produce images from relatively long sections of the coronaries. This approach yields relatively rapid visualization of portions of the coronary tree. Several different methods are compared for interactively moving the scan plane.

Tissue-independent MR tracking of invasive devices with an internal signal source

An improved MR tracking coil design is described that provides more robust tracking performance. It is shown theoretically and experimentally that a coil equipped with an internal spin source increases the signal-to-noise ratio in comparison to a coil system without internal source. The tracking behavior is stable even in air, which is important for laparoscopic work, in which the device has to be placed inside the inflated abdomen. Two types of coil windings were tested: one with the windings perpendicular to the coil axis and the other one with the windings tilted relative to this axis. The experiments showed no significant effect on the signal-to-noise ratio between the two types. The improved MR-tracking coil design with internal source was successfully used in cholecystostomies and in laparoscopic interventions; both procedures were performed on swine.

Visualization of vascular guidewires using MR tracking

MR tracking of a vascular guidewire sized for a .035-inch (.89-mm) catheter lumen was performed. The guidewire was actively tracked by incorporation of a miniature radiofrequency (RF) receive coil built into its tip. After in vitro validation, simultaneous tracking of the guidewire and a catheter was performed in the aortic and abdominal vessels of a swine at 1.5 T. The ability to track such a small device and the ability to simultaneously track multiple devices are significant steps towards vascular interventions under MR guidance.