Numerical modeling of needle tip artifacts in MR gradient echo imaging

High-resolution simulation of B0 field conditions in the human heart from segmented computed tomography images

B₀ inhomogeneity leads to imaging artifacts in cardiac magnetic resonance imaging (MRI), in particular dark band artifacts with steady-state free precession pulse sequences. The limited spatial resolution of MR-derived in vivo B₀ maps and the lack of population data prevent systematic analysis of the problem at hand and the development of optimized B₀ shim strategies. We used readily available clinical computed tomography (CT) images to simulate the B₀ conditions in the human heart at high spatial resolution. Calculated B₀ fields showed consistency with MRI-based B₀ measurements. The B₀ maps for both the simulations and in vivo measurements showed local field inhomogeneities in the vicinity of lung tips with dominant Z3 spherical harmonic terms in the field distribution. The presented simulation approach allows for the derivation of B₀ field conditions at high spatial resolution from CT images and enables the development of subject- and population-specific B₀ shim strategies for the human heart.

Artifact reduction of coaxial needles in magnetic resonance imaging-guided abdominal interventions at 1.5 T: a phantom study

Needle artifacts pose a major limitation for MRI-guided interventions, as they impact the visually perceived needle size and needle-to-target-distance. The objective of this agar liver phantom study was to establish an experimental basis to understand and reduce needle artifact formation during MRI-guided abdominal interventions. Using a vendor-specific prototype fluoroscopic T1-weighted gradient echo sequence with real-time multiplanar acquisition at 1.5 T, the influence of 6 parameters (flip angle, bandwidth, matrix, slice thickness, read-out direction, intervention angle relative to B₀) on artifact formation of 4 different coaxial MR-compatible coaxial needles (Nitinol, 16G–22G) was investigated. As one parameter was modified, the others remained constant. For each individual parameter variation, 2 independent and blinded readers rated artifact diameters at 2 predefined positions (15 mm distance from the perceived needle tip and at 50% of the needle length). Differences between the experimental subgroups were assessed by Bonferroni-corrected non-parametric tests. Correlations between continuous variables were expressed by the Bravais–Pearson coefficient and interrater reliability was quantified using the intraclass classification coefficient. Needle artifact size increased gradually with increasing flip angles (p = 0.002) as well as increasing intervention angles (p < 0.001). Artifact diameters differed significantly between the chosen matrix sizes (p = 0.002) while modifying bandwidth, readout direction, and slice thickness showed no significant differences. Interrater reliability was high (intraclass correlation coefficient 0.776–0.910). To minimize needle artifacts in MRI-guided abdominal interventions while maintaining optimal visibility of the coaxial needle, we suggest medium-range flip angles and low intervention angles relative to B₀.

Minimal artifact actively shimmed metallic needles in MRI

PURPOSE

Signal voids caused by metallic needles pose visualization and monitoring challenges in many MRI applications. In this work, we explore a solution to this problem in the form of an active shim insert that fits inside a needle and corrects the field disturbance (ΔB₀ ) caused by the needle outside of it.

METHODS

The ΔB₀ induced by a 4 mm outside-diameter titanium needle at 3T is modeled and a two-coil orthogonal shim set is designed and fabricated to shim the ΔB₀ . Signal recovery around the needle is assessed in multiple orientations in a water phantom with four different pulse sequences. Phase stability around the needle is assessed in an ex-vivo porcine tissue dynamic gradient echo experiment with and without shimming. Additionally, heating of the shim insert is assessed under 8 min of continuous operation with 1A current and concurrent imaging.

RESULTS

An average recovery of ~63% of lost signal around the needle across orientations is shown with active shimming with a maximum current of 1.172 A. Signal recovery and correction of the underlying ΔB₀ is shown to be independent of imaging sequence. Needle-induced phase gradients outside the perceptible signal void are also minimized with active shimming. Temperature rise of up to 0.9° Celsius is noted over 8 min of continuous 1A active shimming operation.

CONCLUSION

A sequence independent method for minimization of metallic needle induced signal loss using an active shim insert is presented. The method has potential benefits in a range of qualitative and quantitative interventional MRI applications.

MR Imaging‒Guided Intervention: Evaluation of MR Conditional Biopsy and Ablation Needle Tip Artifacts at 3T Using a Balanced Fast Field Echo Sequence

PURPOSE

To systematically investigate artifacts produced by biopsy and ablation needles imaged at various trajectories with respect to the static magnetic field (B0).

MATERIALS AND METHODS

An acrylic phantom was scanned using a rapid balanced fast field echo sequence with 3.0-T magnetic resonance imaging. A 15-gauge microwave needle, a 17-gauge cryoneedle, and an 18-gauge coaxial biopsy needle were imaged in sagittal and axial planes, in 7 different orientations to B0 (0°, 15°, 30°, 45°, 60°, 75°, and 90°). For 4 angles (15°, 30°, 60°, and 75°), images were acquired with the slice orientation aligned to the needle angulation, resulting in the frequency encoding direction being parallel to the needle's long axis for the sagittal slice and perpendicular to the needle angulation for the axial acquisition. The artifact length at the needle tip and maximum artifact width were recorded.

RESULTS

No significant difference was noted in mean artifact length for the cryoneedle (13 mm; 95% confidence interval [CI], 7-19) and coaxial biopsy needle (8 mm; 95% CI, 5-10; P = .08). The mean artifact length was significantly smaller for the microwave ablation needle (1 mm; 95% CI, 0-2; P < .05). The mean artifact width was highest for the coaxial needle (17 mm; 95% CI, 14-19) and significantly higher than the cryoneedle (12 mm; 95% CI, 10-15; P = .024) and microwave ablation needle (8 mm; 95% CI, 6-10; P < .01). The needle tip artifact was significantly smaller when the slice orientation was aligned to the needle angulation for the coaxial and cryoablation needles (P < .01).

CONCLUSIONS

Needle tip artifact length and width increase with increasing angulation to the static field. At large angles (>15°), the needle tip position can be predicted better from images acquired when the slice orientation is aligned to the needle's angulation.

Modeling of active shimming of metallic needles for interventional MRI

PURPOSE

Artifacts caused by large magnetic susceptibility differences between metallic needles and tissue are a persistent problem in many interventional MRI applications. The signal void caused by the needle can hide procedure targets and prevent accurate image-based monitoring. In this paper, a solution to this problem is presented in the form of an active shim insert inspired from degaussing coils used in naval vessels, that is designed to correct the field disturbance (ΔB₀ ) caused by the needle.

METHODS

The ΔB₀ induced by a 10 gauge hollow single-beveled titanium needle at 3T is modeled in different orientations. A set of 63 orthogonal coil pairs with unique tip paths are evaluated for shimming performance, and an optimal coil pair is chosen. Shimming performance and current demands are evaluated over a range of needle orientations.

RESULTS

Robust correction of the titanium needle induced ΔB₀ is predicted using a flat no-loop coil combined with an orthogonal 1½ turn loop coil angled at the bevel angle for most orientations, with currents well below 1 amp per coil. Reductions in ΔB₀ standard deviations with shimming ranged from ~49% to ~10% depending on needle orientation, with performance worsening as the needle is aligned more along B₀ .

CONCLUSION

Simulations predict that it is possible to minimize metallic probe induced ΔB₀ and signal losses using externally supplied direct current shim coil inserts in arbitrary orientations for potential benefits in many interventional MRI applications.

Computation of Magnetic Field Distortions and Impact on T2*-weighted MRI with Application to Magnetic Susceptibility Parameter Estimation

A two-step numerical computation of T₂* signal weighting maps in gradient echo magnetic resonance imaging in the presence of an object with varied susceptibility property is presented. In the first step, the magnetic scalar potential is computed for an arbitrary 2D magnetic susceptibility distribution using an algebraic solver. The corresponding magnetic field disturbance is computed from the magnetic scalar potential. In the second step, nonlinear operations are used to compute T₂* from the magnetic field disturbance and then to generate a map of T₂* signal weighting. The linearity of the first step of the solution process is used to implement a superposition of basis solutions approach that increases computational efficiency. Superposition of basis solutions, computed from a system composed of a single node of differing magnetic susceptibility from the surround, herein referred to as the base system, is found to provide an accurate estimation of the scalar potential for arbitrary susceptibility distributions. Afterwards, nonlinear computation of the T₂* signal weighting maps can be performed. The properties of the algebraic magnetic scalar potential solver are discussed in this work. Finally, the linearity of the magnetic scalar potential solver is used to estimate the magnetic susceptibility of various objects from in vitro MR-imaging data acquired at 9.4T.

In vitro artifact assessment of an MR-compatible, microwave antenna device for percutaneous tumor ablation with fluoroscopic MRI-sequences

OBJECTIVE

To evaluate artifact configuration and diameters of a magnetic resonance (MR) compatible microwave (MW) applicator using near-realtime MR-fluoroscopic sequences for percutaneous tumor ablation procedures.

MATERIAL AND METHODS

Two MW applicators (14 G and 16 G) were tested in an ex-vivo phantom at 1.5 T with two 3 D fluoroscopic sequences: T1-weighted spoiled Gradient Echo (GRE) and T1/T2-weighted Steady State Free Precession (SSFP) sequence. Applicator orientation to main magnetic field (B₀), slice orientation and phase encoding direction (PED) were systematically varied. The influence of these variables was assessed with ANOVA and post-hoc testing.

RESULTS

The artifact was homogenous along the whole length of both antennas with all tested parameters. The tip artifact diameter of the 16 G antenna measured 6.9 ± 1.0 mm, the shaft artifact diameter 8.6 ± 1.2 mm and the Tip Location Error (TLE) was 1.5 ± 1.2 mm.The tip artifact diameter of the 14 G antenna measured 7.7 ± 1.2 mm, the shaft artifact diameter 9.6 ± 1.5 mm and TLE was 1.6 ± 1.2 mm. Orientation to B₀ had no statistically significant influence on tip artifact diameters (16 G: p = .55; 14 G: p = .07) or TLE (16 G: p = .93; 14 G: p = .26). GRE sequences slightly overestimated the antenna length with TLE(16 G) = 2.6 ± 0.5 mm and TLE(14 G) = 2.7 ± 0.7 mm.

CONCLUSIONS

The MR-compatible MW applicator's artifact seems adequate with an acceptable TLE for safe applicator positioning during near-realtime fluoroscopic MR-guidance.

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.

Radial MRI with variable echo times: reducing the orientation dependency of susceptibility artifacts of an MR-safe guidewire

OBJECTIVES

Guidewires are indispensable tools for intravascular MR-guided interventions. Recently, an MR-safe guidewire made from a glass-fiber/epoxy compound material with embedded iron particles was developed. The size of the induced susceptibility artifact, and thus the guidewire's visibility, depends on its orientation against B ₀. We present a radial acquisition scheme with variable echo times that aims to reduce the artifact's orientation dependency.

MATERIALS AND METHODS

The radial acquisition scheme uses sine-squared modulated echo times depending on the physical direction of the spoke to balance the susceptibility artifact of the guidewire. The acquisition scheme was studied in simulations based on dipole fields and in phantom experiments for different orientations of the guidewire against B ₀. The simulated and measured artifact widths were quantitatively compared.

RESULTS

Compared to acquisitions with non-variable echo times, the proposed acquisition scheme shows a reduced angular variability. For the two main orientations (i.e., parallel and perpendicular to B ₀), the ratio of the artifact widths was reduced from about 2.2 (perpendicular vs. parallel) to about 1.2 with the variable echo time approach.

CONCLUSION

The reduction of the orientation dependency of the guidewire's artifact via sine-squared varying echo times could be verified in simulations and measurements. The more balanced artifact allows for a better overall visibility of the guidewire.

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.

How MRI compatible is "MRI compatible"? A systematic comparison of artifacts caused by biopsy needles at 3.0 and 1.5 T

PURPOSE

This study was designed to systematically investigate artifacts caused by interventional needles recommended for use in MRI, with focus on field strength, needle/mandrin type, orientation and sequence.

METHODS

Eight different MRI compatible needles were placed in porcine tissue and examined at 1.5 and 3.0 T with balanced-steady-state-free-precession (B-SSFP) and T1-weighted-spoiled-gradient-echo (T1-SPGR) sequences in different orientations to B0. Artifact diameters with regards to the primary, inner, and secondary, outer artifacts were assessed and statistically evaluated.

RESULTS

The types and degree of artifacts varied considerably, especially between different mandrin types even for the same needles. Orientation of the needle in the magnetic field was another main contributor to the artifact dimensions. Less important factors were the type of pulse sequence and field strength. Artifacts ranged from 0.7 mm (steel, 0°, B-SSFP, 3.0 T, inner) to 71.4 mm (nitinol, 90°, B-SSFP, 1.5 T, outer). Inner artifact diameters in B-SSFP were slightly larger (8.2 ± 5.7 mm) than those in T1-SPGR (7.6 ± 5.4 mm) and comparable between 1.5 and 3.0 T (e.g., 8.0 vs. 8.4 mm, B-SSFP).

CONCLUSIONS

Although all were sold as "MR compatible," the artifacts differed greatly between needle types, and even more so for different mandrins. The results suggest an empirical approach to the needle choice based on lesion type and approach angle.

Single shot MR tagging to quantify local tissue deformation during MRI-guided needle interventions: a feasibility study

PURPOSE

In MRI-guided needle interventions, such as biopsies and brachytherapy, tissue deformation caused by needle movement may result in localization errors and thus hamper the outcome of the procedure. Monitoring the local tissue deformation provides the ability to compensate for it, e.g., by increasing the needle insertion depth. Fast MR scans are useful to track the needle, but cannot be used to quantify local tissue deformation, in case anatomical landmarks are missing. Artificial landmarks can be created by MR tagging. This method provides a spatial saturation pattern (tag) in the tissue. Deformation of this pattern reflects the tissue motion between tag creation and tag imaging. As the needle movement is nonperiodic, k-space cannot be acquired with a multishot approach, like is usually done for cardiac imaging. Hence, a single shot MR tagging sequence is needed, which entails tag creation, needle movement and tag acquisition. In this study, the feasibility of single shot MR tagging for MRI-guided needle interventions in phantom and volunteer experiments is shown.

METHODS

Four different experiments were performed on a 1.5 T MR scanner: the first to quantify translations, the second to quantify rotations, the third to mimic a needle intervention, and the fourth to investigate the tag persistence in a volunteer. The tag pattern is created by a 1331 composite pulse. A balanced steady state free precession sequence is used for imaging. To minimize undesired changes in contrast or sharpness of the tag pattern, we chose a relatively small flip angle and a short imaging time in all experiments. In the volunteer experiments, we modified the sequence to also be able to inspect the influence of the used k-space sampling profile and the flip angle on the temporal persistence of image contrast and tag pattern. In all scans, head or surface coils were used for signal reception.

RESULTS

In all experiments, the tag pattern was clearly visible and could be used to quantify the local tissue deformation caused by (needle) movement. Strong correlations between the actual and measured (angular) phantom motions were obtained. In the needle intervention experiment, the tag lines were perfectly horizontal when there was no needle movement. With needle movement, local tissue displacements up to 5 mm were observed. Volunteer's anatomy could be discriminated, despite the tag pattern. The tag pattern in the prostate, for example, could still be read in all tagging images acquired 2 s after creating the tag pattern. With optimized scan parameters the tag persistence was even longer. The best image tag contrast was obtained using a large flip angle and the profile order low-high, although the image was slightly blurred.

CONCLUSIONS

This study demonstrates that single shot MR tagging can be used to quantify tissue deformation caused by needle movement. The in-vivo tag persistence is sufficient to enable the application of the tagging sequence during MRI-guided needle interventions in patients.

Magnetic field distribution and signal decay in functional mri in very high fields (up to 9.4 T) using monte carlo diffusion modeling

Extravascular signal decay rate R2 or R2 * as a function of blood oxygenation, geometry, and field strength was calculated using a Monte Carlo (MC) algorithm for a wider parameter range than hitherto by others. The relaxation rates of gradient-recalled-echo (GRE) and Hahn-spin-echo (HSE) imaging in the presence of blood vessels (ranging from capillaries to veins) have been computed for a wide range of field strengths up to 9.4T and 50% blood deoxygenation. The maximum HSE decay was found to be shifted to lower radii in higher compared to lower field strengths. For GRE, however, the relaxation rate was greatest for large vessels at any field strength. In addition, assessments of computational reliability have been carried out by investigating the influence of the time step, the Monte Carlo step procedure, boundary conditions, the number of angles between the vessel and the exterior field B0, the influence of neighboring vessels having the same orientation as the central vessel, and the number of proton spins. The results were compared with those obtained from a field distribution of the vessel computed by an analytic formula describing the field distribution of an ideal object (an infinitely long cylinder). It was found that the time step is not critical for values equal to or lower than 200 microseconds. The choice of the MC step procedure (three-dimensional Gaussian diffusion, constant one- or three-dimensional diffusion step) also failed to influence the results significantly; in contrast, the free boundary conditions, as well as taking too few angles into account, did introduce errors. Next neighbor vessels with the same orientation as the main vessel did not contribute significantly to signal decay. The total number of particles simulated was also found to play a minor role in computing R2/ R2 *.

Evaluation of magnetic resonance imaging-compatible needles and interactive sequences for musculoskeletal interventions using an open high-field magnetic resonance imaging scanner

In this article, we study in vitro evaluation of needle artefacts and image quality for musculoskeletal laser-interventions in an open high-field magnetic resonance imaging (MRI) scanner at 1.0T with vertical field orientation. Five commercially available MRI-compatible puncture needles were assessed based on artefact characteristics in a CuSO4 phantom (0.1%) and in human cadaveric lumbar spines. First, six different interventional sequences were evaluated with varying needle orientation to the main magnetic field B0 (0 degrees to 90 degrees ) in a sequence test. Artefact width, needle-tip error, and contrast-to-noise ratio (CNR) were calculated. Second, a gradient-echo sequence used for thermometric monitoring was assessed and in varying echo times, artefact width, tip error, and signal-to-noise ratio (SNR) were measured. Artefact width and needle-tip error correlated with needle material, instrument orientation to B0, and sequence type. Fast spin-echo sequences produced the smallest needle artefacts for all needles, except for the carbon fibre needle (width <3.5 mm, tip error <2 mm) at 45 degrees to B0. Overall, the proton density-weighted spin-echo sequences had the best CNR (CNR(Muscle/Needle) >16.8). Concerning the thermometric gradient echo sequence, artefacts remained <5 mm, and the SNR reached its maximum at an echo time of 15 ms. If needle materials and sequences are accordingly combined, guidance and monitoring of musculoskeletal laser interventions may be feasible in a vertical magnetic field at 1.0T.

In vivo quantification of contrast agent concentration using the induced magnetic field for time-resolved arterial input function measurement with MRI

For pharmacokinetic modeling of tissue physiology, there is great interest in measuring the arterial input function (AIF) from dynamic contrast-enhanced (DCE) magnetic resonance imaging (MRI) using paramagnetic contrast agents. Due to relaxation effects, the measured signal is a nonlinear function of the injected contrast agent concentration and depends on sequence parameters, system calibration, and time-of-flight effects, making it difficult to accurately measure the AIF during the first pass. Paramagnetic contrast agents also affect susceptibility and modify the magnetic field in proportion to their concentration. This information is contained in the MR signal phase which is discarded in a typical image reconstruction. However, quantifying AIF through contrast agent susceptibility induced phase changes is made difficult by the fact that the induced magnetic field is nonlocal and depends upon the contrast agent spatial distribution and thus on organ and vessel shapes. In this article, the contrast agent susceptibility was quantified through inversion of magnetic field shifts using a piece-wise constant model. Its feasibility is demonstrated by a determination of the AIF from the susceptibility-induced field changes of an intravenous bolus. After in vitro validation, a time-resolved two-dimensional (2D) gradient echo scan, triggered to diastole, was performed in vivo on the aortic arch during a bolus injection of 0.1 mmol/kg Gd-DTPA. An approximate geometrical model of the aortic arch constructed from the magnitude images was used to calculate the spatial variation of the field associated with the bolus. In 14 subjects, Gd concentration curves were measured dynamically (one measurement per heart beat) and indirectly validated by independent 2D cine phase contrast flow rate measurements. Flow rate measurements using indicator conservation with this novel quantitative susceptibility imaging technique were found to be in good agreement with those obtained from the cine phase contrast measurements in all subjects. Contrary to techniques that rely on intensity, the accuracy of this signal phase based method is insensitive to factors influencing signal intensity such as flip angle, coil sensitivity, relaxation changes, and time-of-flight effects extending the range of pulse sequences and contrast doses for which quantitative DCE-MRI can be applied.

Biopsy needle tips with markers--MR compatible needles for high-precision needle tip positioning

Needle tip visualization is of high importance in magnetic resonance imaging (MRI) guided interventional procedures, for example for taking biopsies from suspicious lesions in the liver or kidney. The exact position of the needle tip is often obscured by image artifacts arising from the magnetic properties of the needle. The authors investigated two special biopsy needle tip designs using diamagnetic coatings. For common interventional MR sequences, the needle tip can be identified in the MR image by several equidistant dark spots arranged along a straight line. A dotted instead of a solid line allows for an improved control of the movement of the needle, not only if the needle is tilted toward the imaging plane, but also if the needle leaves an empty canal with signal extinction, which cannot be distinguished from the needle material itself. With the proposed design the position of the needle tip can be estimated with a precision of approximately 1 mm using conventional FLASH, FISP, and TSE sequences, as used for interventional MR. Furthermore, the size of the biopsy probe can be estimated from the artifact. In using needles with a properly designed tip coating, taking biopsies under MR control is beginning to be greatly simplified. The approach to design artifacts using diamagnetic material in combination with paramagnetic material paves the way toward new instruments and implants, suitably tailored to the needs of the interventional radiologist.

Magnetic resonance guidance for radiofrequency ablation of liver tumors

Image-guided thermal ablation of liver tumors is a minimally invasive treatment option. Techniques used for thermal ablation are radiofrequency (RF) ablation, laser interstitial thermotherapy (LITT), microwave (MW) ablation, high-intensity focused ultrasound (HIFU), and cryoablation. Among these techniques RF ablation attained widespread consideration. Image guidance should ensure a precise ablation therapy leading to a complete coagulation of tumor tissue without injury to critical structures. Therefore, the modality of image guidance has an important impact on the safety and efficacy of percutaneous RF ablation. The current literature regarding percutaneous RF ablation mainly describes the use of computed tomography (CT) and ultrasonography (US) guidance. In addition, interventional MR systems offer the possibility to utilize the advantages of MR imaging such as excellent soft-tissue contrast, multiplanar and interactive capabilities, and sensitivity to thermal effects during the entire RF ablation procedure. Monitoring of thermally induced coagulation by MR imaging is supportive to control the ablation procedure. MR imaging can be advantageously used to guide overlapping ablation if necessary as well as to define the endpoint of RF ablation after complete coverage of the target tissue is verified. Furthermore, monitoring of thermal effects is essential in order to prevent unintended thermal damage from critical structures surrounding the target region. Therefore, MR-guided RF ablation offers the possibility for a safe and effective therapy option in the treatment of primary and secondary hepatic malignancies. The article summarizes the role of MR guidance for RF ablation of liver tumors.

COMBINED X-RAY AND MAGNETIC RESONANCE IMAGING FACILITY

OBJECTIVE

To assess the feasibility of a hybrid imaging setup combining x-ray and magnetic resonance imaging (MRI) in the setting of both stereotactic and functional neurosurgery.

METHODS

A combined x-ray and MRI scanning facility with a trolley system for a fast patient transfer between both modalities was installed in a neurosurgical setting. A registration algorithm for fusion of MRI scans and x-ray images was derived for augmentation of fluoroscopic x-ray projection images with MRI scan data, such as anatomic structures and planned probe trajectories. Phantom measurements were obtained between both modalities for estimation of registration accuracy. Practical application of our system in stereotactic and functional neurosurgery was tested in brachytherapy, deep brain stimulation, and motor cortex stimulation.

RESULTS

Phantom measurements yielded a mean spatial deviation of 0.7 +/- 0.3 mm with a maximum deviation of 1.1 mm for MRI scans versus x-rays. Augmentation of x-ray images with MRI scan data allowed intraoperative verification of the planned trajectory and target in three types of neurosurgical procedures: positioning iodine seeds in brachytherapy in one case with cerebellar metastasis, placement of electrodes for deep brain stimulation in two cases of advanced Parkinson's disease, and placement of an epidural grid for motor cortex stimulation in two cases of intractable pain.

CONCLUSION

Combined x-ray and MRI-guided stereotactic and functional neurosurgery is feasible. Augmentation of x-ray projection images with MRI scan data, such as planned probe trajectories and MRI scan segmented anatomic structures may be beneficial for probe guidance in stereotactic and functional neurosurgery.

Magnetic resonance imaging of phase transitions in nitinol

Magnetic resonance images are prone to artifacts caused by metallic objects. Apart from being a source of image degradation, such artifacts can also provide information about the magnetic properties of the foreign object. In this work, we aim to explore the potential of magnetic resonance imaging to detect and characterize changes in magnetic properties of nitinol undergoing temperature- or strain-induced phase changes. A spin echo and a gradient echo method were used to measure the magnetization changes related to the phase transformations. Results of both methods were in agreement and in accordance with the independent measurements using a vibrating sample magnetometer. Magnetic resonance imaging turned out to be a suitable method to visualize and quantify magnetization and phase changes in situ. It is not restricted to a single imaging strategy and does not require any modification of the test object. The results indicate the potential of magnetic resonance imaging to provide direct feedback of the thermomechanical state of the alloy.

Combined passive and active shimming for in vivo MR spectroscopy at high magnetic fields

The use of high magnetic fields increases the sensitivity and spectral dispersion in magnetic resonance spectroscopy (MRS) of brain metabolites. Practical limitations arise, however, from susceptibility-induced field distortions, which are increased at higher magnetic field strengths. Solutions to this problem include optimized shimming, provided that active, i.e., electronic, shimming can operate over a sufficient range. To meet our shim requirements, which were an order of magnitude greater than the active shim capacity of our 7T MR system, we developed a combined passive and active shim approach. Simple geometries of ferromagnetic shim elements were derived and numerically optimized to generate a complete set of second-order spherical harmonic shim functions in a modular manner. The major goals of the shim design were maximization of shim field accuracy and ease of practical implementation. The theoretically optimized ferro-shim geometries were mounted on a cylindrical surface and placed inside the magnet bore, surrounding the subject's head and the RF coil. Passive shimming generated very strong shim fields and eliminated the worst of the field distortions, after which the field was further optimized by flexible and highly accurate active shimming. Here, the passive-shimming procedure was first evaluated theoretically, then applied in phantom studies and subsequently validated for in vivo 1H MRS in the macaque visual cortex. No artifacts due to the passive shim setup were observed; adjustments were reproducible between sessions. The modularity and the reduction to two pieces per shim term in this study is an important simplification that makes the method applicable also for passive shimming within single sessions. The feasibility of very strong, flexible and high-quality shimming via a combined approach of passive and active shimming is of great practical relevance for MR imaging and spectroscopy at high field strengths where shim power is limited or where shimming of specific anatomical regions inherently requires strong shim fields.

Geometry and extension of signal voids in MR images induced by aggregations of magnetically labelled cells

The ability of magnetic resonance imaging (MRI) to visualize magnetically labelled cells has attracted much attention for revealing cellular events. The present study addressed the geometry and the extension of signal voids in static signal dephasing MRI induced by aggregations of magnetically labelled cells by means of a three-dimensional numerical model. The magnetic field distortions around spherical cell aggregations were treated as equivalent to those of a magnetic dipole. Intravoxel signal dephasing and respective signal voids attributed to these field inhomogeneities were computed. Effects of cell concentration on the signal void in the plane of view were evaluated in terms of dipole magnetization. Signal void characteristics were scrutinized systematically for fundamental sequence parameters including echo time, voxel size and plane-of-view orientation. For all variables examined, significant changes in geometry as well as extension of signal voids were demonstrated. The results are of crucial importance to optimize and interpret MR images with regard to spatial accuracy as well as sensitivity to detect aggregations of labelled cells in vitro or even in vivo. It is anticipated that the dependence of the extension of signal voids on the local magnetization may be valuable for quantifying labelled cells.

Numerical simulations of intra-voxel dephasing effects and signal voids in gradient echo MR imaging using different sub-grid sizes

Signal void artifacts in gradient echo imaging are caused by the intra-voxel dephasing of the spins. Intra-voxel dephasing can be estimated by computing the field distribution on a sub-grid inside each picture element, followed by integration of all magnetization components. The strategy of computing the artifacts based on the integration of the sub-voxel signal components is presented here for different sub-grids. The coarseness of the sub-grid is directly related to computational effort. The possibility to save memory space and computing time for the dipole model by computing the field only on a sub-grid is addressed in the presented article. It is investigated as to how far computational time and memory space can be reduced by using an appropriate sub-grid. Numerical results for a model of a partially diamagnetically coated needle shaft are compared to experimental findings. In the case of a pure titanium needle, it is shown as being sufficient to compute the field distribution on a sub-grid that is at least four times coarser in each direction than the grid used to discretize the object in the related MR image. Due to three nested loops over the 3D grid, the need for memory space and time is saved by a factor 64. Deviations between measurements and simulations for the broad side of the artifact (uncompensated) and for the small side of the artifact (compensated) were 15.5%, respectively, 19.1% for orientation parallel to the exterior field, and 22.7%, respectively, 23.1% for orientation perpendicular to the exterior field.

Effect of Spatial Distribution of Magnetic Dipoles on Lamor Frequency Distribution and MR Signal Decay – a Numerical Approach Under Static Dephasing Conditions

Cells loaded with superparamagnetic iron oxide (SPIO) cause relatively strong magnetic field distortions, implying that field position effects of neighboring SPIO loaded cells have to be accounted for. We treated SPIO loaded cells as magnetic dipoles in a homogeneous magnetic field and computed the 3D frequency distribution and the related signal decay using a numerical approach under static dephasing conditions. The volume fraction of dipoles was kept constant for all simulations. For larger randomly distributed magnetic dipoles we found a non-Lorentzian frequency distribution and a non-monoexponential signal decay whereas, for smaller dipoles, the frequency distribution was more Lorentzian and the signal decay was well fitted monoexponentially. Moreover, based on our numerical and experimental findings, we found the gradient echo signal decay due to a single SPIO labeled cell to be non-monoexponential. The numerical approach provides deeper understanding of how the spatial distribution of SPIO loaded cells affects the MR signal decay. This fact has to be considered for the in vivo quantification of SPIO loaded cells, implying that in tissues with different spatial distributions of identical SPIO concentrations, different signal decays might be observed.

Magnetic susceptibility effects on the accuracy of MR temperature monitoring by the proton resonance frequency method

PURPOSE

To evaluate the error of MR temperature assessment based on the temperature-dependent Larmor frequency shift of water protons, which can result from susceptibility effects caused by the radiofrequency (RF) applicator.

MATERIALS AND METHODS

Local frequency shifts due to RF applicator displacements were simulated numerically by means of a three-dimensional elementary dipole model. Experimental examinations using a water tank phantom equipped with a high-precision screw thread were applied to examine temperature and movement effects for five commercially available, MR-compatible RF applicators. Measurements were performed at 1.5 Tesla.

RESULTS

For single-needle electrodes perpendicular to the external field, a distortion of 0.1 ppm and 0.2 ppm was recorded at a distance of 17.5 mm and 12.5 mm, respectively, to the needle shaft. Cluster applicators and umbrella-shaped applicators caused distortions of 0.1 ppm up to distances of 36 mm. Sinusoidal dependence on applicator orientation was found with the highest values for perpendicular orientation and the lowest values for orientation parallel to the magnetic field. With a single electrode oriented perpendicular to the field at a distance of 1.5 cm and 2.0 cm, a needle displacement of 5 mm led to an error in temperature measurement of 16.3 degrees C and 7.5 degrees C, respectively.

CONCLUSION

In MR temperature measurement, displacement of the RF applicator by patient movement or breathing leads to significant errors that have to be taken into account when PRF temperature maps are used to monitor tumor ablation in the presence of paramagnetic applicators.

Compensation of magnetic field distortions from paramagnetic instruments by added diamagnetic material: Measurements and numerical simulations

In minimally invasive procedures guided by magnetic resonance (MR) imaging instruments usually are made of titanium or titanium alloys (e.g., nitinol), because other more MR-compatible materials often cannot provide sufficient mechanical properties. Artifacts depending on susceptibility arise in MR images due to incorrect spatial encoding and intravoxel dephasing and thereby hamper the surgeon's view onto the region of interest. To overcome the artifact problem, compensation of the paramagnetic properties by diamagnetic coating or filling of the instruments has been proposed in the literature. We used a numerical modeling procedure to estimate the effect of compensation. Modeling of the perturbation of the static magnetic field close to the instruments reflects the underlying problem and is much faster and cost efficient than manufacturing prototypes and measuring artifact behavior of these prototypes in the MR scanner. A numerical model based on the decomposition of the susceptibility distribution in elementary dipoles was developed by us. The program code was written object oriented to allow for both maximum computational speed and minimum random access memory. We used System International units throughout the modeling for the magnetic field, allowing absolute quantification of the magnetic field disturbance. The field outside a simulated needlelike instrument, modeled by a paramagnetic cylinder (out of titan, chi =181.1) of length 8.0 mm and of diameter 1.0 mm, coated with a diamagnetic layer (out of bismuth, chi=-165.0) of thickness 0, 0.1, 0.2, 0.3, and 0.4 mm, was found to be best compensated if the cross-sectional area of the cylinder, multiplied by the absolute susceptibility value of the cylinder material, is equal to the cross-sectional area of the coating, multiplied by the absolute susceptibility value of the coating material. At the extremity of the coated cylinder an uncompensated field distortion was found to remain. We studied various tip shapes and geometries using our computational model: Suitable diamagnetic coating or filling of paramagnetic instruments clearly reduced tip artifacts and diminished the dependency of artifact size on orientation of the instrument with respect to B0 in the numerical studies. We verified the results of the simulations by measuring coated and uncoated titanium wires in a 1.5 T MR scanner.