For real-time 3D visualization of respiratory organ motion for MRI-guided therapy, a new adaptive 4D MR imaging method based on navigator echo and multiple gating windows was developed. This method was designed to acquire a time series of volumetric 3D images of a cyclically moving organ, enabling therapy to be guided by synchronizing the 4D image with the actual organ motion in real time. The proposed method was implemented in an open-configuration 0.5T clinical MR scanner. To evaluate the feasibility and determine optimal imaging conditions, studies were conducted with a phantom, volunteers, and a patient. In the phantom study the root mean square (RMS) position error in the 4D image of the cyclically moving phantom was 1.9 mm and the imaging time was approximately 10 min when the 4D image had six frames. In the patient study, 4D images were successfully acquired under clinical conditions and a liver tumor was discriminated in the series of frames. The image quality was affected by the relations among the encoding direction, the slice orientation, and the direction of motion of the target organ. In conclusion, this study has shown that the proposed method is feasible and capable of providing a real-time dynamic 3D atlas for surgical navigation with sufficient accuracy and image quality.
We report an automated method to simultaneously monitor blood-oxygenation-level-dependent (BOLD) MR signals from multiple cortical areas in real-time. Individual brain anatomy was normalized and registered to a pre-segmented atlas in standardized anatomical space. Subsequently, using real-time fMRI (rtfMRI) data acquisition, localized BOLD signals were measured and displayed from user-selected areas labeled with anatomical and Brodmann's Area (BA) nomenclature. The method was tested on healthy volunteers during the performance of hand motor and internal speech generation tasks employing a trial-based design. Our data normalization and registration algorithm, along with image reconstruction, movement correction and a data display routine were executed with enough processing and communication bandwidth necessary for real-time operation. Task-specific BOLD signals were observed from the hand motor and language areas. One of the study participants was allowed to freely engage in hand clenching tasks, and associated brain activities were detected from the motor-related neural substrates without prior knowledge of the task onset time. The proposed method may be applied to various applications such as neurofeedback, brain-computer-interface, and functional mapping for surgical planning where real-time monitoring of region-specific brain activity is needed.
Diffusion imaging with high-b factors, high spatial resolution and cerebrospinal fluid signal suppression was performed in order to characterize the biexponential nature of the diffusion-related signal decay with b-factor in normal cortical gray and deep gray matter (GM). Integration of inversion pulses with a line scan diffusion imaging sequence resulted in 91% cerebrospinal fluid signal suppression, permitting accurate measurement of the fast diffusion coefficient in cortical GM (1.142+/-0.106 microm2/ms) and revealing a marked similarity with that found in frontal white matter (WM) (1.155+/-0.046 microm2/ms). The reversal of contrast between GM and WM at low vs high b-factors is shown to be due to a significantly faster slow diffusion coefficient in cortical GM (0.338+/-0.027 microm2/ms) than in frontal WM (0.125+/-0.014 microm2/ms). The same characteristic diffusion differences between GM and WM are observed in other brain tissue structures. The relative component size showed nonsignificant differences among all tissues investigated. Cellular architecture in GM and WM are fundamentally different and may explain the two- to threefold higher slow diffusion coefficient in GM.
This work investigated the effect of ultrasonic frequency on the threshold for blood-brain barrier (BBB) disruption induced by ultrasound pulses combined with an ultrasound contrast agent. Experiments were performed in rabbits using pulsed sonications at 2.04 MHz with peak pressure amplitudes ranging from 0.3 to 2.3 MPa. BBB disruption was evaluated using contrast-enhanced magnetic resonance imaging. The threshold for BBB disruption was estimated using probit regression. Representative samples with similar amounts of contrast enhancement were examined in light microscopy. Results from these experiments were compared with data from previous studies that used ultrasound frequencies between 0.26 and 1.63 MHz. We found that the BBB disruption threshold (value where the probability for disruption was estimated to be 50%) expressed in terms of the peak negative pressure amplitude increased as a function of the frequency. It appeared to be constant, however, when the exposures were expressed as a function of the mechanical index (peak negative pressure amplitude estimated in situ divided by square root of frequency). Regression of data from all frequencies resulted in an estimated mechanical index threshold of 0.46 (95% confidence intervals: 0.42 to 0.50). Histologic examination of representative samples with similar amounts of blood-brain barrier disruption found that the number of regions containing extravasated red blood cells per unit area was substantially lower on average for lower ultrasound frequencies. This data suggests that the mechanical index is a meaningful metric for ultrasound-induced blood-brain barrier disruption, at least for when other parameters that are not taken into account by the mechanical index are not varied. It also suggests that lower frequency sonication produces less red blood cell extravasation per unit area.
Parallel MRI (pMRI) achieves imaging acceleration by partially substituting gradient-encoding steps with spatial information contained in the component coils of the acquisition array. Variable-density subsampling in pMRI was previously shown to yield improved two-dimensional (2D) imaging in comparison to uniform subsampling, but has yet to be used routinely in clinical practice. In an effort to reduce acquisition time for 3D fast spin-echo (3D-FSE) sequences, this work explores a specific nonuniform sampling scheme for 3D imaging, subsampling along two phase-encoding (PE) directions on a rectilinear grid. We use two reconstruction methods-2D-GRAPPA-Operator and 2D-SPACE RIP-and present a comparison between them. We show that high-quality images can be reconstructed using both techniques. To evaluate the proposed sampling method and reconstruction schemes, results via simulation, phantom study, and in vivo 3D human data are shown. We find that fewer artifacts can be seen in the 2D-SPACE RIP reconstructions than in 2D-GRAPPA-Operator reconstructions, with comparable reconstruction times.
OBJECTIVE: Preoperative magnetic resonance imaging (MRI), functional MRI, diffusion tensor MRI, magnetic resonance spectroscopy, and positron-emission tomographic scans may be aligned to intraoperative MRI to enhance visualization and navigation during image-guided neurosurgery. However, several effects (both machine- and patient-induced distortions) lead to significant geometric distortion of intraoperative MRI. Therefore, a precise alignment of these image modalities requires correction of the geometric distortion. We propose and evaluate a novel method to compensate for the geometric distortion of intraoperative 0.5-T MRI in image-guided neurosurgery. METHODS: In this initial pilot study, 11 neurosurgical procedures were prospectively enrolled. The scheme used to correct the geometric distortion is based on a nonrigid registration algorithm introduced by our group. This registration scheme uses image features to establish correspondence between images. It estimates a smooth geometric distortion compensation field by regularizing the displacements estimated at the correspondences. A patient-specific linear elastic material model is used to achieve the regularization. The geometry of intraoperative images (0.5 T) is changed so that the images match the preoperative MRI scans (3 T). RESULTS: We compared the alignment between preoperative and intraoperative imaging using 1) only rigid registration without correction of the geometric distortion, and 2) rigid registration and compensation for the geometric distortion. We evaluated the success of the geometric distortion correction algorithm by measuring the Hausdorff distance between boundaries in the 3-T and 0.5-T MRIs after rigid registration alone and with the addition of geometric distortion correction of the 0.5-T MRI. Overall, the mean magnitude of the geometric distortion measured on the intraoperative images is 10.3 mm with a minimum of 2.91 mm and a maximum of 21.5 mm. The measured accuracy of the geometric distortion compensation algorithm is 1.93 mm. There is a statistically significant difference between the accuracy of the alignment of preoperative and intraoperative images, both with and without the correction of geometric distortion (P < 0.001). CONCLUSION: The major contributions of this study are 1) identification of geometric distortion of intraoperative images relative to preoperative images, 2) measurement of the geometric distortion, 3) application of nonrigid registration to compensate for geometric distortion during neurosurgery, 4) measurement of residual distortion after geometric distortion correction, and 5) phantom study to quantify geometric distortion.
We created an image-guided robot system to assist with skull base drilling by integrating a robot, a commercial navigation system, and an open source visualization platform. The objective of this procedure is to create a cavity in the skull base to allow access for neurosurgical interventions. The motivation for introducing an image-guided robot is to improve safety by preventing the surgeon from accidentally damaging critical structures during the drilling procedure. Our approach is to attach the cutting tool to the robot end-effector and operate the robot in a cooperative control mode, where robot motion is determined from the forces and torques applied by the surgeon. We employ "virtual fixtures" to constrain the motion of the cutting tool so that it remains in the safe zone that was defined on a preoperative CT scan. This paper presents the system design and the results of phantom and cadaveric experiments. Both experiments have demonstrated the feasibility of the system, with average overcut error at about 1 mm and maximum errors at 2.5 mm.
The combination of the imaging abilities of magnetic resonance imaging (MRI) with the ability to delivery energy to targets deep in the body noninvasively with focused ultrasound presents a disruptive technology with the potential to significantly affect healthcare. MRI offers precise targeting, visualization, and quantification of temperature changes and the ability to immediately evaluate the treatment. By exploiting different mechanisms, focused ultrasound offers a range of therapies, ranging from thermal ablation to targeted drug delivery. This article reviews recent preclinical and tests clinical of this technology.
Previous studies have investigated a potential method for targeted drug delivery in the central nervous system that uses focused ultrasound bursts combined with an ultrasound contrast agent to temporarily disrupt the blood-brain barrier (BBB). The purpose of this work was to investigate the integrity of the tight junctions (TJs) in rat brain microvessels after this BBB disruption. Ultrasound bursts (1.5-MHz) in combination with a gas contrast agent (Optison) was applied at two locations in the brain in 25 rats to induce BBB disruption. Using immunoelectron microscopy, the distributions of the TJ-specific transmembrane proteins occludin, claudin-1, claudin-5, and of submembranous ZO-1 were examined at 1, 2, 4, 6 and 24 h after sonication. A quantitative evaluation of the protein expression was made by counting the number of immunosignals per micrometer in the junctional clefts. BBB disruption at the sonicated locations was confirmed by the leakage of i.v. administered horseradish peroxidase (HRP, m.w. 40,000 Da) and lanthanum chloride (La(3+), m.w. approximately 139 Da). Leakage of these agents was observed at 1 and 2 h and, in a few vessels, at 4 h after ultrasound application. These changes were paralleled by the apparent disintegration of the TJ complexes, as evidenced by the redistribution and loss of the immunosignals for occludin, claudin-5 and ZO-1. Claudin-1 seemed less involved. At 6 and 24 h after sonication, no HRP or lanthanum leakage was observed and the barrier function of the TJs, as indicated by the localization and density of immunosignals, appeared to be completely restored. This study provides the first direct evidence that ultrasound bursts combined with a gas contrast agent cause disassembling of the TJ molecular structure, leading to loss of the junctional barrier functions in brain microvessels. The BBB disruption appears to last up to 4 h after sonication and permits the paracellular passage of agents with molecular weights up to at least 40 kDa. These promising features can be exploited in the future development of this method that could enable the delivery of drugs, antibodies or genes to targeted locations in the brain.
Previously, it was shown that low-intensity focused ultrasound pulses applied along with an ultrasound contrast agent results in temporary blood-brain barrier (BBB) disruption. This effect could be used for targeted drug delivery in the central nervous system. This study examined the effects of burst length, pulse repetition frequency (PRF), and ultrasound contrast agent dose on the resulting BBB disruption. One hundred nonoverlapping brain locations were sonicated through a craniotomy in experiments in 26 rabbits (ultrasound frequency: 0.69 MHz, burst: 0.1, 1, 10 ms, PRF: 0.5, 1, 2, 5 Hz, duration: 20 s, peak negative pressure amplitude: 0.1 to 1.5 MPa, Optison dosage 50, 100, 250 microl/kg). For each sonication, BBB disruption was evaluated using contrast-enhanced magnetic resonance imaging. The BBB disruption threshold (the pressure amplitude yielding a 50% probability for BBB disruption) was determined using probit regression for the three burst lengths tested. Tissue effects were examined in light microscopy for representative locations with similar amounts of contrast enhancement from each group. While changing the PRF or the Optison dosage did not result in a significant difference in the magnitude of the BBB disruption (p > 0.05), reducing the burst length resulted in significantly less contrast enhancement (p < 0.01). The BBB disruption thresholds were estimated to be 0.69, 0.47 and 0.36 MPa for 0.1, 1 and 10 ms bursts, respectively. No difference was detected in histology between any experimental groups. This data suggests that over the range of parameters tested, BBB disruption is not affected by PRF or ultrasound contrast agent dose. However, both the BBB disruption magnitude and its threshold depend on the burst length.
PURPOSE: To clinically assess a previously described method (Rieke et.al., Magn Reson Med 2004) to produce more motion-robust MRI-based temperature images using data acquired during MRI-guided focused ultrasound surgery (MRgFUS) of uterine fibroids. MATERIALS AND METHODS: The method ("referenceless thermometry") uses surface fitting in nonheated regions of individual phase images to extrapolate and then remove background phase variations that are unrelated to temperature changes. We tested this method using images from 100 sonications selected from 33 patient MRgFUS treatments. Temperature measurements and thermal dose contours estimated with the referenceless method were compared with those produced with the standard phase-difference technique. Fitting accuracy and noise level were also measured. RESULTS: In 92/100 sonications, the difference between the two measurements was less than 3 degrees C. The average difference in the measurements was 1.5 +/- 1.4 degrees C. Small motion artifacts were observed in the phase-difference imaging when the difference was greater than 3 degrees C. The method failed in two cases. The mean absolute error in the surface fit in baseline images corresponded to a temperature error of 0.8 +/- 1.4 degrees C. The noise level was approximately 40% lower than the phase-difference method. Thermal dose contours calculated from the two methods agreed well on average. CONCLUSION: Based on the small error when compared with the standard technique, this method appears to be adequate for temperature monitoring of MRgFUS in uterine fibroids and may prove useful for monitoring temperature changes in moving organs.
A method for validating the start-to-end accuracy of a 3-D ultrasound (US)-based patient positioning system for radiotherapy is described. A radiosensitive polymer gel is used to record the actual dose delivered to a rigid phantom after being positioned using 3-D US guidance. Comparison of the delivered dose with the treatment plan allows accuracy of the entire radiotherapy treatment process, from simulation to 3-D US guidance, and finally delivery of radiation, to be evaluated. The 3-D US patient positioning system has a number of features for achieving high accuracy and reducing operator dependence. These include using tracked 3-D US scans of the target anatomy acquired using a dedicated 3-D ultrasound probe during both the simulation and treatment sessions, automatic 3-D US-to-US registration and use of infrared LED (IRED) markers of the optical position-sensing system for registering simulation computed tomography to US data. The mean target localization accuracy of this system was 2.5 mm for four target locations inside the phantom, compared with 1.6 mm obtained using the conventional patient positioning method of laser alignment. Because the phantom is rigid, this represents the best possible set-up accuracy of the system. Thus, these results suggest that 3-D US-based target localization is practically feasible and potentially capable of increasing the accuracy of patient positioning for radiotherapy in sites where day-to-day organ shifts are greater than 1 mm in magnitude.
In the last decade, high intensity focused ultrasound (HIFU) has gained popularity as a minimally invasive and noninvasive therapeutic tool for treatment of cancers, arrhythmias, and other medical conditions. HIFU therapy is often guided by magnetic resonance imaging (MRI), which provides anatomical images for therapeutic device placement, temperature maps for treatment guidance, and postoperative evaluation of the region of interest. While piezoelectric transducers are dominantly used for MR-guided HIFU, capacitive micromachined ultrasonic transducers (CMUTs) show competitive advantages, such as ease of fabrication, integration with electronics, improved efficiency, and reduction of self-heating. In this paper, we will show our first results of an unfocused CMUT transducer monitored by MR-temperature maps. This 2.51 mm by 2.32 mm, unfocused CMUT heated a HIFU phantom by 14 degrees C in 2.5 min. This temperature rise was successfully monitored by MR thermometry in a 3.0 T General Electric scanner.
Ultrasound-based methods for temperature monitoring could greatly assist focused ultrasound visualization and treatment planning based on sound speed-induced change in phase as a function of temperature. A method is presented that uses reflex transmission integration, planar projection, and tomographic reconstruction techniques to visualize phase contrast by measuring the sound field before and after heat deposition. Results from experiments and numerical simulations employing a through-transmission setup are presented to demonstrate feasibility of using phase contrast methods for identifying temperature change. A 1.088-MHz focused transducer was used to interrogate a medium with a phase contrast feature, following measurement of the baseline reference field with a hydrophone. A thermal plume in water and a tissue phantom with multiple water columns was used in separate experiments to produce a phase contrast. The reference and phase contrast field scans were numerically backprojected and the phase difference correctly identified the position and orientation of the features. The peak temperature reconstructed from the phase shift was within 0.2 degrees C of the measured temperature in the plume. Simulated results were in good agreement with experimental results. Finally, employment of reflex transmission imaging techniques for adopting a pulse-echo arrangement was simulated, and its future experimental application is discussed.
We describe a method for correcting the distortions present in echo planar images (EPI) and registering the EPI to structural MRI. A fieldmap is predicted from an air / tissue segmentation of the MRI using a perturbation method and subsequently used to unwarp the EPI data. Shim and other missing parameters are estimated by registration. We obtain results that are similar to those obtained using fieldmaps, however neither fieldmaps, nor knowledge of shim coefficients is required.
Language fMRI has been used to study brain regions involved in language processing and has been applied to pre-surgical language mapping. However, in order to provide clinicians with optimal information, the sensitivity and specificity of language fMRI needs to be improved. Type II error of failing to reach statistical significance when the language activations are genuinely present may be particularly relevant to pre-surgical planning, by falsely indicating low surgical risk in areas where no activations are shown. Furthermore, since the execution of language paradigms involves cognitive processes other than language function per se, the conventional general linear model (GLM) method may identify non-language-specific activations. In this study, we assessed an exploratory approach, independent component analysis (ICA), as a potential complementary method to the inferential GLM method in language mapping applications. We specifically investigated whether this approach might reduce type II error as well as generate more language-specific maps. Fourteen right-handed healthy subjects were studied with fMRI during two word generation tasks. A similarity analysis across tasks was proposed to select components of interest. Union analysis was performed on the language-specific components to increase sensitivity, and conjunction analysis was performed to identify language areas more likely to be essential. Compared with GLM, ICA identified more activated voxels in the putative language areas, and signals from other sources were isolated into different components. Encouraging results from one brain tumor patient are also presented. ICA may be used as a complementary tool to GLM in improving pre-surgical language mapping.
The aim of this study is to evaluate the feasibility of using harmonic cancellation for a therapeutic ultrasound transducer excited by a switched-mode power converter without an additional output filter. A switching waveform without the third harmonic was created by cascading two switched-mode power inverter modules at which their output waveforms were pi/3 phase shifted from each other. A PSPICE simulation model for the power converter output stage was developed. The simulated results were in good agreement with the measurement. The waveform and harmonic contents of the acoustic pressure generated by a 1-MHz, self-focused piezoelectric transducer with and without harmonic cancellation have been evaluated. Measured results indicated that the acoustic third harmonicto- fundamental ratio at the focus was small (-48 dB) with harmonic cancellation, compared to that without harmonic cancellation (-20 dB). The measured acoustic levels of the fifth harmonic for both cases with and without harmonic cancellation also were small (-46 dB) compared to the fundamental. This study shows that it is viable to drive a piezoelectric ultrasound transducer using a switched-mode power converter without the requirement of an additional output filter in many high-intensity focused ultrasound (HIFU) applications.
The partial Fourier gradient-echo echo planar imaging (EPI) technique makes it possible to acquire high-resolution functional MRI (fMRI) data at an optimal echo time. This technique is especially important for fMRI studies at high magnetic fields, where the optimal echo time is short and may not be achieved with a full Fourier acquisition scheme. In addition, it has been shown that partial Fourier EPI provides better anatomic resolvability than full Fourier EPI. However, the partial Fourier gradient-echo EPI may be degraded by artifacts that are not usually seen in other types of imaging. Those unique artifacts in partial Fourier gradient-echo EPI, to our knowledge, have not yet been systematically evaluated. Here we use the k-space energy spectrum analysis method to understand and characterize two types of partial Fourier EPI artifacts. Our studies show that Type 1 artifact, originating from k-space energy loss, cannot be corrected with pure postprocessing, and Type 2 artifact can be eliminated with an improved reconstruction method. We propose a novel algorithm, that combines images obtained from two or more reconstruction schemes guided by k-space energy spectrum analysis, to generate partial Fourier EPI with greatly reduced Type 2 artifact. Quality control procedures for avoiding Type 1 artifact in partial Fourier EPI are also discussed.
Independent component analysis (ICA) of fMRI data generates session/individual specific brain activation maps without a priori assumptions regarding the timing or pattern of the blood-oxygenation-level-dependent (BOLD) signal responses. However, because of a random permutation among output components, ICA does not offer a straightforward solution for the inference of group-level activation. In this study, we present an independent vector analysis (IVA) method to address the permutation problem during fMRI group data analysis. In comparison to ICA, IVA offers an analysis of additional dependent components, which were assigned for use in the automated grouping of dependent activation patterns across subjects. Upon testing using simulated trial-based fMRI data, our proposed method was applied to real fMRI data employing both a single-trial task-paradigm (right hand motor clenching and internal speech generation tasks) and a three-trial task-paradigm (right hand motor imagery task). A generalized linear model (GLM) and the group ICA of the fMRI toolbox (GIFT) were also applied to the same data set for comparison to IVA. Compared to GLM, IVA successfully captured activation patterns even when the functional areas showed variable hemodynamic responses that deviated from a hypothesized response. We also showed that IVA effectively inferred group-activation patterns of unknown origins without the requirement for a pre-processing stage (such as data concatenation in ICA-based GIFT). IVA can be used as a potential alternative or an adjunct to current ICA-based fMRI group processing methods.
BACKGROUND: We developed an image-guided robot system to provide mechanical assistance for skull base drilling, which is performed to gain access for some neurosurgical interventions, such as tumour resection. The motivation for introducing this robot was to improve safety by preventing the surgeon from accidentally damaging critical neurovascular structures during the drilling procedure.
METHODS: We integrated a Stealthstation navigation system, a NeuroMate robotic arm with a six-degree-of-freedom force sensor, and the 3D Slicer visualization software to allow the robotic arm to be used in a navigated, cooperatively-controlled fashion by the surgeon. We employed virtual fixtures to constrain the motion of the robot-held cutting tool, so that it remained in the safe zone that was defined on a preoperative CT scan.
RESULTS: We performed experiments on both foam skull and cadaver heads. The results for foam blocks cut using different registrations yielded an average placement error of 0.6 mm and an average dimensional error of 0.6 mm. We drilled the posterior porus acusticus in three cadaver heads and concluded that the robot-assisted procedure is clinically feasible and provides some ergonomic benefits, such as stabilizing the drill. We obtained postoperative CT scans of the cadaver heads to assess the accuracy and found that some bone outside the virtual fixture boundary was cut. The typical overcut was 1-2 mm, with a maximum overcut of about 3 mm.
CONCLUSIONS: The image-guided cooperatively-controlled robot system can improve the safety and ergonomics of skull base drilling by stabilizing the drill and enforcing virtual fixtures to protect critical neurovascular structures. The next step is to improve the accuracy so that the overcut can be reduced to a more clinically acceptable value of about 1 mm.
OBJECTIVE: Accurate biopsy sampling of the suspected lesions is critical for the diagnosis and clinical management of prostate cancer. Transperineal in-bore MRI-guided prostate biopsy (tpMRgBx) is a targeted biopsy technique that was shown to be safe, efficient, and accurate. Our goal was to develop an open source software platform to support evaluation, refinement, and translation of this biopsy approach. METHODS: We developed SliceTracker, a 3D Slicer extension to support tpMRgBx. We followed modular design of the implementation to enable customization of the interface and interchange of image segmentation and registration components to assess their effect on the processing time, precision, and accuracy of the biopsy needle placement. The platform and supporting documentation were developed to enable the use of software by an operator with minimal technical training to facilitate translation. Retrospective evaluation studied registration accuracy, effect of the prostate segmentation approach, and re-identification time of biopsy targets. Prospective evaluation focused on the total procedure time and biopsy targeting error (BTE). RESULTS: Evaluation utilized data from 73 retrospective and ten prospective tpMRgBx cases. Mean landmark registration error for retrospective evaluation was 1.88 ± 2.63 mm, and was not sensitive to the approach used for prostate gland segmentation. Prospectively, we observed target re-identification time of 4.60 ± 2.40 min and BTE of 2.40 ± 0.98 mm. CONCLUSION: SliceTracker is modular and extensible open source platform for supporting image processing aspects of the tpMRgBx procedure. It has been successfully utilized to support clinical research procedures at our site.
Patient-mounted needle guide devices for percutaneous ablation are vulnerable to patient motion. The objective of this study is to develop and evaluate a software system for an MRI-compatible patient-mounted needle guide device that can adaptively compensate for displacement of the device due to patient motion using a novel image-based automatic device-to-image registration technique. We have developed a software system for an MRI-compatible patient-mounted needle guide device for percutaneous ablation. It features fully-automated image-based device-to-image registration to track the device position, and a device controller to adjust the needle trajectory to compensate for the displacement of the device. We performed: (a) a phantom study using a clinical MR scanner to evaluate registration performance; (b) simulations using intraoperative time-series MR data acquired in 20 clinical cases of MRI-guided renal cryoablations to assess its impact on motion compensation; and (c) a pilot clinical study in three patients to test its feasibility during the clinical procedure. FRE, TRE, and success rate of device-to-image registration were [Formula: see text] mm, [Formula: see text] mm, and 98.3% for the phantom images. The simulation study showed that the motion compensation reduced the targeting error for needle placement from 8.2 mm to 5.4 mm (p < 0.0005) in patients under general anesthesia (GA), and from 14.4 mm to 10.0 mm ([Formula: see text]) in patients under monitored anesthesia care (MAC). The pilot study showed that the software registered the device successfully in a clinical setting. Our simulation study demonstrated that the software system could significantly improve targeting accuracy in patients treated under both MAC and GA. Intraprocedural image-based device-to-image registration was feasible.
PURPOSE: To develop and evaluate an approach to estimate the respiratory-induced motion of lesions in the chest and abdomen. MATERIALS AND METHODS: The proposed approach uses the motion of an initial reference needle inserted into a moving organ to estimate the lesion (target) displacement that is caused by respiration. The needles position is measured using an inertial measurement unit (IMU) sensor externally attached to the hub of an initially placed reference needle. Data obtained from the IMU sensor and the target motion are used to train a learning-based approach to estimate the position of the moving target. An experimental platform was designed to mimic respiratory motion of the liver. Liver motion profiles of human subjects provided inputs to the experimental platform. Variables including the insertion angle, target depth, target motion velocity and target proximity to the reference needle were evaluated by measuring the error of the estimated target position and processing time. RESULTS: The mean error of estimation of the target position ranged between 0.86 and 1.29 mm. The processing maximum training and testing time was 5 ms which is suitable for real-time target motion estimation using the needle position sensor. CONCLUSION: The external motion of an initially placed reference needle inserted into a moving organ can be used as a surrogate, measurable and accessible signal to estimate in real-time the position of a moving target caused by respiration; this technique could then be used to guide the placement of subsequently inserted needles directly into the target.
Brain shift during tumor resection compromises the spatial validity of registered preoperative imaging data that is critical to image-guided procedures. One current clinical solution to mitigate the effects is to reimage using intraoperative magnetic resonance (iMR) imaging. Although iMR has demonstrated benefits in accounting for preoperative-to-intraoperative tissue changes, its cost and encumbrance have limited its widespread adoption. While iMR will likely continue to be employed for challenging cases, a cost-effective model-based brain shift compensation strategy is desirable as a complementary technology for standard resections. We performed a retrospective study of [Formula: see text] tumor resection cases, comparing iMR measurements with intraoperative brain shift compensation predicted by our model-based strategy, driven by sparse intraoperative cortical surface data. For quantitative assessment, homologous subsurface targets near the tumors were selected on preoperative MR and iMR images. Once rigidly registered, intraoperative shift measurements were determined and subsequently compared to model-predicted counterparts as estimated by the brain shift correction framework. When considering moderate and high shift ([Formula: see text], [Formula: see text] measurements per case), the alignment error due to brain shift reduced from [Formula: see text] to [Formula: see text], representing [Formula: see text] correction. These first steps toward validation are promising for model-based strategies.
OBJECTIVE: The purpose of this article is to report our intermediate to long-term outcomes with image-guided percutaneous hepatic tumor cryoablation and to evaluate its technical success, technique efficacy, local tumor progression, and adverse event rate. MATERIALS AND METHODS: Between 1998 and 2014, 299 hepatic tumors (243 metastases and 56 primary tumors; mean diameter, 2.5 cm; median diameter, 2.2 cm; range, 0.3-7.8 cm) in 186 patients (95 women; mean age, 60.9 years; range, 29-88 years) underwent cryoablation during 236 procedures using CT (n = 126), MRI (n = 100), or PET/CT (n = 10) guidance. Technical success, technique efficacy at 3 months, local tumor progression (mean follow-up, 2.5 years; range, 2 months to 14.6 years), and adverse event rates were calculated. RESULTS: The technical success rate was 94.6% (279/295). The technique efficacy rate was 89.5% (231/258) and was greater for tumors smaller than 4 cm (93.4%; 213/228) than for larger tumors (60.0%; 18/30) (p < 0.0001). Local tumor progression occurred in 23.3% (60/258) of tumors and was significantly more common after the treatment of tumors 4 cm or larger (63.3%; 19/30) compared with smaller tumors (18.0%; 41/228) (p < 0.0001). Adverse events followed 33.8% (80/236) of procedures and were grade 3-5 in 10.6% (25/236) of cases. Grade 3 or greater adverse events more commonly followed the treatment of larger tumors (19.5%; 8/41) compared with smaller tumors (8.7%; 17/195) (p = 0.04). CONCLUSION: Image-guided percutaneous cryoablation of hepatic tumors is efficacious; however, tumors smaller than 4 cm are more likely to be treated successfully and without an adverse event.