In algorithms for processing diffusion tensor images, two common ingredients are interpolating tensors, and measuring the distance between them. We propose a new class of interpolation paths for tensors, termed geodesic-loxodromes, which explicitly preserve clinically important tensor attributes, such as mean diffusivity or fractional anisotropy, while using basic differential geometry to interpolate tensor orientation. This contrasts with previous Riemannian and Log-Euclidean methods that preserve the determinant. Path integrals of tangents of geodesic-loxodromes generate novel measures of over-all difference between two tensors, and of difference in shape and in orientation.
OBJECTIVE: The objective of our study was to report the occurrence of benign inflammatory nodules that develop in or near applicator tracks after percutaneous radiofrequency ablation and cryoablation of renal tumors.
CONCLUSION: Benign inflammatory nodules occur rarely after percutaneous ablation of renal tumors and may mimic tumor seeding of the applicator track.
BACKGROUND: The authors prospectively evaluated the late gastrointestinal (GI) and genitourinary (GU) toxicity and prostate-specific antigen (PSA) control of magnetic resonance imaging (MRI)-guided brachytherapy used as salvage for radiation therapy (RT) failure. METHODS: From October 2000 to October 2005, 25 men with a rising PSA level and biopsy-proven, intraprostatic cancer at least 2 years after initial RT (external beam in 13 men and brachytherapy in 12 men) who had favorable clinical features (Gleason score < or =7, PSA < 10 ng/mL, negative pelvic and bone imaging studies), received MRI-guided salvage brachytherapy to a minimum peripheral dose of 137 gray on a phase 1/2 protocol. Estimates of toxicity and cancer control were calculated using the Kaplan-Meier method. RESULTS: The median follow-up was 47 months. The 4-year estimate of grade 3 or 4 GI or GU toxicity was 30%, and 13% of patients required a colostomy and/or urostomy to repair a fistula. An interval < 4.5 years between RT courses was associated with both outcomes with a hazard ratio of 12 (95% confidence interval [95% CI], 1.4-100; P = .02) for grade 3 or 4 toxicity and 25 (95% CI, 1.1-529; P = .04) for colostomy and/or urostomy. PSA control (nadir +2 definition) was 70% at 4 years. CONCLUSIONS: The current results indicated that MRI-guided salvage brachytherapy in men who are selected based on presenting characteristics and post-failure PSA kinetics can achieve high PSA control rates, although complications requiring surgical intervention may occur in 10% to 15% of patients. Prospective randomized studies are needed to characterize the relative cancer control and toxicity after all forms of salvage local therapy.
We formalize the pair-wise registration problem in a maximum a posteriori (MAP) framework that employs a multinomial model of joint intensities with parameters for which we only have a prior distribution. To obtain an MAP estimate of the aligning transformation alone, we treat the multinomial parameters as nuisance parameters, and marginalize them out. If the prior on those is uninformative, the marginalization leads to registration by minimization of joint entropy. With an informative prior, the marginalization leads to minimization of the entropy of the data pooled with pseudo observations from the prior. In addition, we show that the marginalized objective function can be optimized by the Expectation-Maximization (EM) algorithm, which yields a simple and effective iteration for solving entropy-based registration problems. Experimentally, we demonstrate the effectiveness of the resulting EM iteration for rapidly solving a challenging intra-operative registration problem.
We present an algorithm to generate samples from probability distributions on the space of curves. We view a traditional curve evolution energy functional as a negative log probability distribution and sample from it using a Markov chain Monte Carlo (MCMC) algorithm. We define a proposal distribution by generating smooth perturbations to the normal of the curve and show how to compute the transition probabilities to ensure that the samples come from the posterior distribution. We demonstrate some advantages of sampling methods such as robustness to local minima, better characterization of multi-modal distributions, access to some measures of estimation error, and ability to easily incorporate constraints on the curve.
We sought to determine the safety and feasibility of percutaneous MRI-guided cryotherapy in the care of patients with refractory or painful metastatic lesions of soft tissue and bone adjacent to critical structures.
MATERIALS AND METHODS:
Twenty-seven biopsy-proven metastatic lesions of soft tissue (n = 17) and bone (n = 10) in 22 patients (15 men, seven women; age range, 24-85 years) were managed with MRI-guided percutaneous cryotherapy. The mean lesion diameter was 5.2 cm. Each lesion was adjacent to or encasing one or more critical structures, including bowel, bladder, and major blood vessels. A 0.5-T open interventional MRI system was used for cryoprobe placement and ice-ball monitoring. Complications were assessed for all treatments. CT or MRI was used to determine local control of 21 tumors. Pain palliation was assessed clinically in 19 cases. The mean follow-up period was 19.5 weeks.
Twenty-two (81%) of 27 tumors were managed without injury to adjacent critical structures. Two patients had transient lower extremity numbness, and two had both urinary retention and transient lower extremity paresthesia. One patient had chronic serous vaginal discharge, and one sustained a femoral neck fracture at the ablation site 6 weeks after treatment. Thirteen (62%) of the 21 tumors for which follow-up information was available either remained the same size as before treatment or regressed. Eight tumors progressed (mean local progression-free interval, 5.6 months; range, 3-18 months). Pain was palliated in 17 of 19 patients; six of the 17 experienced complete relief, and 11 had partial relief.
MRI-guided percutaneous cryotherapy for metastatic lesions of soft tissue and bone adjacent to critical structures is safe and can provide local tumor control and pain relief in most patients.
The purpose of this study is to examine the feasibility of simultaneous Magnetic Resonance Imaging (MRI) and Ultrasound (US) imaging in visualizing anatomical structures and functions in human carotid arteries. US has high frame rate in visualizing dynamic changes while high resolution MRI is capable of capturing volumetric structures with the best tissue contrast. Concurrent multi-modal image acquisition allows fusion of US Doppler flow measurement with volumetric MRI. We present a method for acquiring MR images in a known orientation with respect to US image by passive fiducial tracking and demonstrate concurrent real-time imaging in the right Common Carotid Artery (CCA) in both modalities. Preliminary results suggest that US and MRI can operate concurrently with proper shielding. Dispensability measurements are feasible on both modalities at the co-incident plane.
This paper presents a novel multiscale shape representation and segmentation algorithm based on the spherical wavelet transform. This work is motivated by the need to compactly and accurately encode variations at multiple scales in the shape representation in order to drive the segmentation and shape analysis of deep brain structures, such as the caudate nucleus or the hippocampus. Our proposed shape representation can be optimized to compactly encode shape variations in a population at the needed scale and spatial locations, enabling the construction of more descriptive, nonglobal, nonuniform shape probability priors to be included in the segmentation and shape analysis framework. In particular, this representation addresses the shortcomings of techniques that learn a global shape prior at a single scale of analysis and cannot represent fine, local variations in a population of shapes in the presence of a limited dataset. Specifically, our technique defines a multiscale parametric model of surfaces belonging to the same population using a compact set of spherical wavelets targeted to that population. We further refine the shape representation by separating into groups wavelet coefficients that describe independent global and/or local biological variations in the population, using spectral graph partitioning. We then learn a prior probability distribution induced over each group to explicitly encode these variations at different scales and spatial locations. Based on this representation, we derive a parametric active surface evolution using the multiscale prior coefficients as parameters for our optimization procedure to naturally include the prior for segmentation. Additionally, the optimization method can be applied in a coarse-to-fine manner. We apply our algorithm to two different brain structures, the caudate nucleus and the hippocampus, of interest in the study of schizophrenia. We show: 1) a reconstruction task of a test set to validate the expressiveness of our multiscale prior and 2) a segmentation task. In the reconstruction task, our results show that for a given training set size, our algorithm significantly improves the approximation of shapes in a testing set over the Point Distribution Model, which tends to oversmooth data. In the segmentation task, our validation shows our algorithm is computationally efficient and outperforms the Active Shape Model algorithm, by capturing finer shape details.
A minimalist numerical model of white matter is presented, the objective of which is to help provide a biological basis for improved diffusion tensor imaging (DTI) analysis. Water diffuses, relaxes, and exchanges in three compartments-intracellular, extracellular, and myelin sheath. Exchange between compartments is defined so as to depend on the diffusion coefficients and the compartment sizes. Based on the model, it is proposed that an additive "baseline tensor" that correlates with intraaxonal water volume be included in the computation. Anisotropy and tortuosity calculated from such analysis may correspond better to tract ultrastructure than if calculated without the baseline. According to the model, reduced extracellular volume causes increased baseline and reduced apparent diffusion. Depending on the pulse sequence, reduced permeability can cause an increase in both the baseline and apparent diffusion.
A system was developed for real-time electrocardiogram (ECG) analysis and artifact correction during magnetic resonance (MR) scanning, to improve patient monitoring and triggering of MR data acquisitions. Based on the assumption that artifact production by magnetic field gradient switching represents a linear time invariant process, a noise cancellation (NC) method is applied to ECG artifact linear prediction. This linear prediction is performed using a digital finite impulse response (FIR) matrix, that is computed employing ECG and gradient waveforms recorded during a training scan. The FIR filters are used during further scanning to predict artifacts by convolution of the gradient waveforms. Subtracting the artifacts from the raw ECG signal produces the correction with minimal delay. Validation of the system was performed both off-line, using prerecorded signals, and under actual examination conditions. The method is implemented using a specially designed Signal Analyzer and Event Controller (SAEC) computer and electronics. Real-time operation was demonstrated at 1 kHz with a delay of only 1 ms introduced by the processing. The system opens the possibility of automatic monitoring algorithms for electrophysiological signals in the MR environment.
OBJECTIVE: The usefulness of neurosurgical navigation with current visualizations is seriously compromised by brain shift, which inevitably occurs during the course of the operation, significantly degrading the precise alignment between the pre-operative MR data and the intra-operative shape of the brain. Our objectives were (i) to evaluate the feasibility of non-rigid registration that compensates for the brain deformations within the time constraints imposed by neurosurgery, and (ii) to create augmented reality visualizations of critical structural and functional brain regions during neurosurgery using pre-operatively acquired fMRI and DT-MRI. MATERIALS AND METHODS: Eleven consecutive patients with supratentorial gliomas were included in our study. All underwent surgery at our intra-operative MR imaging-guided therapy facility and have tumors in eloquent brain areas (e.g. precentral gyrus and cortico-spinal tract). Functional MRI and DT-MRI, together with MPRAGE and T2w structural MRI were acquired at 3 T prior to surgery. SPGR and T2w images were acquired with a 0.5 T magnet during each procedure. Quantitative assessment of the alignment accuracy was carried out and compared with current state-of-the-art systems based only on rigid registration. RESULTS: Alignment between pre-operative and intra-operative datasets was successfully carried out during surgery for all patients. Overall, the mean residual displacement remaining after non-rigid registration was 1.82 mm. There is a statistically significant improvement in alignment accuracy utilizing our non-rigid registration in comparison to the currently used technology (p<0.001). CONCLUSIONS: We were able to achieve intra-operative rigid and non-rigid registration of (1) pre-operative structural MRI with intra-operative T1w MRI; (2) pre-operative fMRI with intra-operative T1w MRI, and (3) pre-operative DT-MRI with intra-operative T1w MRI. The registration algorithms as implemented were sufficiently robust and rapid to meet the hard real-time constraints of intra-operative surgical decision making. The validation experiments demonstrate that we can accurately compensate for the deformation of the brain and thus can construct an augmented reality visualization to aid the surgeon.
In the United States, unenhanced CT is currently the most common imaging modality used to guide percutaneous biopsy and tumor ablation. The majority of liver tumors such as hepatocellular carcinomas are visible on contrast-enhanced CT or MRI obtained prior to the procedure. Yet, these tumors may not be visible or may have poor margin conspicuity on unenhanced CT images acquired during the procedure. Non-rigid registration has been used to align images accurately, even in the presence of organ motion. However, to date, it has not been used clinically for radiofrequency ablation (RFA), since it requires significant computational infrastructure and often these methods are not sufficient robust. We have already introduced a novel finite element based method (FEM) that is demonstrated to achieve good accuracy and robustness for the problem of brain shift in neurosurgery. In this current study, we adapt it to fuse pre-procedural MRI with intra-procedural CT of liver. We also compare its performance with conventional rigid registration and two non-rigid registration methods: b-spline and demons on 13 retrospective datasets from patients that underwent RFA at our institution. FEM non-rigid registration technique was significantly better than rigid (p < 10-5), non-rigid b-spline (p < 10-4) and demons (p < 10-4) registration techniques. The results of our study indicate that this novel technology may be used to optimize placement of RF applicator during CT-guided ablations.
Using direct cortical stimulation to map language function during awake craniotomy is a well-described and useful technique. However, the optimum neuropsychological tasks to use have not been detailed. We used both functional MRI (fMRI) and direct cortical stimulation to compare the sensitivity of two behavioral paradigms, number counting and object naming, in the demonstration of eloquent cortical language areas. Fifteen patients with left hemisphere lesions and seven healthy control subjects participated. Patients had both preoperative fMRI at 3 T and direct cortical stimulation. Patients and controls performed object naming and number counting during fMRI at 3 T. Laterality indices were calculated from the fMRI maps for the Number-counting>Object-naming and Object-naming>Number-counting contrasts. The same number-counting and object-naming paradigms were tested during awake craniotomy and assessed for sensitivity to speech disruption. In all patients during intraoperative cortical stimulation, speech disruption occurred at more sites during object naming than during number counting. Subtle speech errors were only elicited with the object-naming paradigm, whereas only speech arrest and/or hypophonia were measured using the number counting paradigm. In both patients and controls, fMRI activation maps demonstrated greater left lateralization for object naming as compared to number counting in both frontal and temporal language areas. Number counting resulted in a more bihemispheric distribution of activations than object naming. Both cortical stimulation testing and fMRI suggest that automated speech tasks such as number counting may not fully engage putative language networks and therefore are not optimal for language localization for surgical planning.
In this paper a comprehensive framework for pre-operative planning, procedural skill training, and intraoperative navigation is presented. The goal of this system is to integrate surgical simulation with surgical planning in order to improve the individual treatment of patients. Various surgical approaches and new, more complex procedures can be assessed using a safe and objective platform that will allow the physicians to explore and discuss possible risks and benefits prior to the intervention. A simulation environment extends the pre-operative planning in a natural way, as it allows for direct evaluation of the surgical approach envisioned for each case. In addition, by providing intraoperative navigation based on this simulation, surgeons can carry out the previously optimized plan with higher precision and greater confidence.
This contribution presents finite element computation of the deformation field within the brain during craniotomy-induced brain shift. The results were used to illustrate the capabilities of non-linear (i.e. accounting for both geometric and material non-linearities) finite element analysis in non-rigid registration of pre- and intra-operative magnetic resonance images of the brain. We used patient-specific hexahedron-dominant finite element mesh, together with realistic material properties for the brain tissue and appropriate contact conditions at boundaries. The model was loaded by the enforced motion of nodes (i.e. through prescribed motion of a boundary) at the brain surface in the craniotomy area. We suggest using explicit time-integration scheme for discretised equations of motion, as the computational times are much shorter and accuracy, for practical purposes, the same as in the case of implicit integration schemes. Application of the computed deformation field to register (i.e. align) the pre-operative images with the intra-operative ones indicated that the model very accurately predicts the displacements of the tumour and the lateral ventricles even for limited information about the brain surface deformation. The prediction accuracy improves when information about deformation of not only exposed (during craniotomy) but also unexposed parts of the brain surface is used when prescribing loading. However, it appears that the accuracy achieved using information only about the deformation of the exposed surface, that can be determined without intra-operative imaging, is acceptable. The presented results show that non-linear biomechanical models can complement medical image processing techniques when conducting non-rigid registration. Important advantage of such models over the previously used linear ones is that they do not require unrealistic assumptions that brain deformations are infinitesimally small and brain stress-strain relationship is linear.
A novel framework for joint clustering and point-by-point mapping of white matter fiber pathways is presented. Accurate clustering of the trajectories into fiber bundles requires point correspondence determined along the fiber pathways. This knowledge is also crucial for any tract-oriented quantitative analysis. We employ an expectation-maximization (EM) algorithm to cluster the trajectories in a Gamma mixture model context. The result of clustering is the probabilistic assignment of the fiber trajectories to each cluster, an estimate of the cluster parameters, and point correspondences along the trajectories. Point-by-point correspondence of the trajectories within a bundle is obtained by constructing a distance map and a label map from each cluster center at every iteration of the EM algorithm. This offers a time-efficient alternative to pairwise curve matching of all trajectories with respect to each cluster center. Probabilistic assignment of the trajectories to clusters is controlled by imposing a minimum threshold on the membership probabilities, to remove outliers in a principled way. The presented results confirm the efficiency and effectiveness of the proposed framework for quantitative analysis of diffusion tensor MRI.
BACKGROUND: Mapping functional areas of the brain is important for planning tumour resections. With the increased use of functional magnetic resonance imaging (fMRI) for presurgical planning, there is a need to validate that fMRI activation mapping is consistent with the mapping obtained during surgery using direct electrocortical stimulation (DECS). METHODS: A quantitative comparison of DECS and fMRI mapping techniques was performed, using a patient-specific conductivity model to find the current distribution resulting from each stimulation site. The resulting DECS stimulation map was compared to the fMRI activation map, using the maximal Dice similarity coefficient (MDSC). RESULTS: Our results show some agreement between these two mapping techniques--the stimulation site with the largest MOSC was the only site that demonstrated intra-operative effect. CONCLUSIONS: There is a substantial effort to improve the techniques used to map functional areas, particularly using fMRI. It seems likely that fMRI will eventually provide a valid non-invasive means for functional mapping.
OBJECTIVES: The determination of eloquent cortex is essential when planning neurosurgical approaches to brain lesions. This study examined the abilities of medical personnel of various backgrounds to predict the location of functional cortex using anatomical information provided by MR imaging.
PATIENTS AND METHODS: Neurosurgeons, neuroscientists, neuroradiologists, medical students and MR technologists viewed anatomical MR images acquired from patients with brain tumors and healthy controls. These five groups of raters were then asked to locate the primary motor hand, supplementary motor and primary auditory areas and their predictions were compared to fMRI data acquired from the same subjects.
RESULTS: The overall mean distance from the center of the fMRI activation was 2.38 cm. The neuroscientists performed the best and MR technologists performed the worst (mean distance from center of 1.83 and 3.04 cm, respectively, p<0.05). The difference between patients and controls was not significant. The mean distance by ROI was primary motor hand 2.03 cm, auditory area 2.06 cm and supplementary motor area 3.18 cm (p<0.05). Raters also performed best in the medial-lateral direction, compared to superior-inferior and anterior-posterior directions (mean distances from center 0.42, 1.04 and 1.81 cm, respectively). Finally, the approximate minimum fields of view necessary to capture the entire fMRI activations using the raters' predictions ranged from 5 to 15 cm, or 3 to 12 cm larger than the fMRI activations.
CONCLUSION: Medical personnel of various training perform poorly when using only anatomical information to predict the location of functional areas of cortex.
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.