PURPOSE: To quantify needle placement accuracy of magnetic resonance image (MRI)-guided core needle biopsy of the prostate. MATERIALS AND METHODS: A total of 10 biopsies were performed with 18-gauge (G) core biopsy needle via a percutaneous transperineal approach. Needle placement error was assessed by comparing the coordinates of preplanned targets with the needle tip measured from the intraprocedural coherent gradient echo images. The source of these errors was subsequently investigated by measuring displacement caused by needle deflection and needle susceptibility artifact shift in controlled phantom studies. Needle placement error due to misalignment of the needle template guide was also evaluated. RESULTS: The mean and standard deviation (SD) of errors in targeted biopsies was 6.5 +/- 3.5 mm. Phantom experiments showed significant placement error due to needle deflection with a needle with an asymmetrically beveled tip (3.2-8.7 mm depending on tissue type) but significantly smaller error with a symmetrical bevel (0.6-1.1 mm). Needle susceptibility artifacts observed a shift of 1.6 +/- 0.4 mm from the true needle axis. Misalignment of the needle template guide contributed an error of 1.5 +/- 0.3 mm. CONCLUSION: Needle placement error was clinically significant in MRI-guided biopsy for diagnosis of prostate cancer. Needle placement error due to needle deflection was the most significant cause of error, especially for needles with an asymmetrical bevel.
Traditional ultrasound imaging methods rely on the bandwidth and center frequency of transduction to achieve axial and radial image resolution, respectively. In this study, a new modality for spatially localizing scattering targets in a two-dimensional field is presented. In this method, the bandwidth of field excitation is high, and the center frequency is lowered such that the corresponding wavelengths are substantially larger than the target profiles. Furthermore, full two-dimensional field measurements are obtained with single send-receive sequences, demonstrating a substantial simplification of the traditional scanning techniques. Field reconstruction is based on temporal-spectral cross-correlations between measured backscatter data and a library of region of interest (ROI) backscatter data measured a priori. The transducer design is based upon a wedge-shaped geometry, which was shown to yield spatially frequency-separated bandwidths of up to 156% with center frequencies of 1.38 MHz. Initial results with these send-and-receive transducer parameters and cylindrical reflection targets in a 10-mm x 10-mm ROI demonstrate two-dimensional target localization to within 0.5 mm. Spatial localization of point scatterers is demonstrated for single and multiple scattering sites.
A method is described for detecting scattering in two-dimensions using an unfocused ultrasound field created from a continuously driven source array. The frequency of each element on the array is unique, resulting in a field that is highly variant as a function of both time and position. The scattered signal is then received by a single receiving line. The method, as currently written, is valid under the first order Born approximation. To demonstrate the approach, a series of simulations within the frequency range of 0.10-1.25 MHz are performed and compared with a simulated B-Scan in the same frequency range. The method is found to be superior in resolving closely spaced objects, discerning 1.4 mm separation in the radial and 0.5-mm separation in the axial direction. The method was also better able to determine object size, resolving scatters less than 10% of wavelength associated with the center frequency.
We have developed a method to use low-intensity focused ultrasound pulses combined with an ultrasound contrast agent to produce temporary blood-brain barrier disruption (BBBD). This method could provide a means for the targeted delivery of drugs or imaging agents into the brain. In all our previous work, we used Optison as the ultrasound contrast agent. The purpose of this study was to test the feasibility of using the contrast agent Definity for BBBD. A total of 36 non-overlapping locations were sonicated through a craniotomy in experiments in the brains of nine rabbits (four locations per rabbit; ultrasound [US] frequency: 0.69 MHz; burst: 10 ms; pulse repetition frequency (PRF): 1 Hz; duration: 20 s). The peak negative pressure amplitude ranged from 0.2 to 1.5 MPa. An additional 11 locations were sonicated using Optison at pressure amplitude of 0.5 MPa. Definity and Optison dosages were the same as those used clinically for ultrasound imaging: 10 and 50 microl/kg, respectively. The probability for BBBD (determined using MRI contrast agent enhancement) as a function of pressure amplitude was similar to that found earlier with Optison. For both agents, the probability was estimated to be 50% at 0.4 MPa using probit regression. Histologic examination revealed small, isolated areas of extravasated erythrocytes in some locations. At 0.8 MPa and higher, these areas were sometimes accompanied by tiny (dimensions of 100 microm or less) regions of damaged brain parenchyma. The magnitude of the BBBD was larger with Optison than with Definity at 0.5 MPa (signal enhancement: 13.3% +/- 4.4% vs. 8.4% +/- 4.9%; p = 0.04). In addition, more areas with extravasated erythrocytes were observed with Optison (5.0 +/- 3.5 vs. 1.4 +/- 1.9 areas with extravasation in histology section with largest effect; p = 0.03). We concluded that BBBD is possible using Definity at the dosage of contrast agent and the acoustic parameters tested in this study. The probability for BBBD as a function of pressure amplitude and the type of acute tissue effects were similar to what has been observed using Optison. However, under the experimental conditions used in this study, Optison produced a larger effect for the same acoustic pressure amplitude.
The logarithm of the odds ratio (LogOdds) is frequently used in areas such as artificial neural networks, economics, and biology, as an alternative representation of probabilities. Here, we use LogOdds to place probabilistic atlases in a linear vector space. This representation has several useful properties for medical imaging. For example, it not only encodes the shape of multiple anatomical structures but also captures some information concerning uncertainty. We demonstrate that the resulting vector space operations of addition and scalar multiplication have natural probabilistic interpretations. We discuss several examples for placing label maps into the space of LogOdds. First, we relate signed distance maps, a widely used implicit shape representation, to LogOdds and compare it to an alternative that is based on smoothing by spatial Gaussians. We find that the LogOdds approach better preserves shapes in a complex multiple object setting. In the second example, we capture the uncertainty of boundary locations by mapping multiple label maps of the same object into the LogOdds space. Third, we define a framework for non-convex interpolations among atlases that capture different time points in the aging process of a population. We evaluate the accuracy of our representation by generating a deformable shape atlas that captures the variations of anatomical shapes across a population. The deformable atlas is the result of a principal component analysis within the LogOdds space. This atlas is integrated into an existing segmentation approach for MR images. We compare the performance of the resulting implementation in segmenting 20 test cases to a similar approach that uses a more standard shape model that is based on signed distance maps. On this data set, the Bayesian classification model with our new representation outperformed the other approaches in segmenting subcortical structures.
PURPOSE: To prospectively assess patient response (after 12 months) to magnetic resonance (MR) imaging-guided focused ultrasound surgery in treatment of uterine leiomyomas by using two treatment protocols.
MATERIALS AND METHODS: This prospective clinical trial was approved by institutional review boards and was HIPAA compliant. After giving informed consent, patients with symptomatic leiomyomas were consecutively enrolled and treated at one of five U.S. centers by using an original or a modified protocol. Outcomes were assessed with the symptom severity score (SSS) obtained at baseline and 3, 6, and 12 months after treatment. Adverse events (AEs) were recorded. Statistical analysis included Student t test, Fisher exact test, analysis of covariance, Spearman correlation, and logistic regression.
RESULTS: One hundred sixty patients had a mean SSS of 62.1 +/- 16.3 (standard deviation) at baseline, which decreased to 35.5 +/- 19.5 at 3 months (P<.001) and to 32.3 +/- 19.8 at 6 months (P<.001) and was 32.7 +/- 21.0 at 12 months (P<.001). Ninety-six patients (mean age, 46.0 years +/- 4.6) were treated with an original protocol, and 64 (mean age, 45.9 years +/- 3.9) were treated with a modified protocol. Patients in the modified group had a significantly greater SSS decrease at 3 months (P=.037) than those in the original group, and 73% of those in the original group and 91% of those in the modified group reported a significant decrease in SSS (of 10 points or greater) at 12 months. No serious AEs were recorded. Fewer AEs were reported in the modified group than in the original group (25% vs 13% reporting no event). Of evaluable patients, fewer in the modified group chose alternative treatment (28%) than in the original group (37%).
CONCLUSION: MR imaging-guided focused ultrasound surgery results in symptomatic improvement, sustained to 12 months after treatment. Treatment with a modified protocol results in greater clinical effectiveness and fewer AEs.
RATIONALE AND OBJECTIVES: We introduce a validation framework for the segmentation of brain tumors from magnetic resonance (MR) images. A novel unsupervised semiautomatic brain tumor segmentation algorithm is also presented.
MATERIALS AND METHODS: The proposed framework consists of 1) T1-weighted MR images of patients with brain tumors, 2) segmentation of brain tumors performed by four independent experts, 3) segmentation of brain tumors generated by a semiautomatic algorithm, and 4) a software tool that estimates the performance of segmentation algorithms.
RESULTS: We demonstrate the validation of the novel segmentation algorithm within the proposed framework. We show its performance and compare it with existent segmentation. The image datasets and software are available at http://www.brain-tumor-repository.org/.
CONCLUSIONS: We present an Internet resource that provides access to MR brain tumor image data and segmentation that can be openly used by the research community. Its purpose is to encourage the development and evaluation of segmentation methods by providing raw test and image data, human expert segmentation results, and methods for comparing segmentation results.
RATIONALE AND OBJECTIVES: To perform a retrospective, quantitative assessment of the anatomic relationship between intra-axial, supratentorial, primary brain tumors, and adjacent white matter fiber tracts based on anatomic and diffusion tensor magnetic resonance imaging (MRI). We hypothesized that white matter infiltration may be common among different types of tumor.
MATERIAL AND METHODS: Preoperative, anatomic (T1- and T2-weighted), and LINESCAN diffusion tensor MRI were obtained in 12 patients harboring supratentorial gliomas (World Health Organization [WHO] Grades II and III). The two imaging modalities were rigidly registered. The tumors were manually segmented from the T1- and T2-weighted MRI, and their volume calculated. A three-dimensional tractography was performed in each case. A second segmentation and volume measurement was performed on the tumor regions intersecting adjacent white matter fiber tracts. Statistical methods included summary statistics to examine the fraction of tumor volume infiltrating adjacent white matter.
RESULTS: There were five patients with low-grade oligodendroglioma (WHO Grade II), one with low-grade mixed oligoastrocytoma (WHO Grade II), one with ganglioglioma, two with low-grade astrocytoma (WHO Grade II), and three with anaplastic astrocytoma (WHO Grade III). We identified white matter tracts infiltrated by tumor in all 12 cases. The median tumor volume (+/- standard deviation) in our patient population was 42.5 +/- 28.9 mL. The median tumor volume (+/- standard deviation) infiltrating white matter fiber tracts was 5.2 +/- 9.9 mL. The median percentage of tumor volume infiltrating white matter fiber tracts was 21.4% +/- 9.7%.
CONCLUSIONS: The information provided by diffusion tensor imaging combined with anatomic MRI might be useful for neurosurgical planning and intraoperative guidance. Our results confirm previous reports that extensive white matter infiltration by primary brain tumors is a common occurrence. However, prospective, large population studies are required to definitively clarify this issue, and how infiltration relates to histologic tumor type, tumor size, and location.
A classical neural tract tracer, WGA-HRP, was injected at multiple sites within the brain of a macaque monkey. Histological sections of the labeled fiber tracts were reconstructed in 3D, and the fibers were segmented and registered with the anatomical post-mortem MRI from the same animal. Fiber tracing along the same pathways was performed on the DTI data using a classical diffusion tracing technique. The fibers derived from the DTI were compared with those segmented from the histology in order to evaluate the performance of DTI fiber tracing. While there was generally good agreement between the two methods, our results reveal certain limitations of DTI tractography, particularly at regions of fiber tract crossing or bifurcation.
Some clinical applications, such as surgical planning, require volumetric models of anatomical structures represented as a set of tetrahedra. A practical method of constructing anatomical models from medical images is presented. The method starts with a set of contours segmented from the medical images by a clinician and produces a model that has high fidelity with the contours. Unlike most modeling methods, the contours are not restricted to lie on parallel planes. The main steps are a 3D Delaunay tetrahedralization, culling of non-object tetrahedra, and refinement of the tetrahedral mesh. The result is a high-quality set of tetrahedra whose surface points are guaranteed to match the original contours. The key is to use the distance map and bit volume structures that were created along with the contours. The method is demonstrated on computed tomography, MRI and 3D ultrasound data. Models of 170,000 tetrahedra are constructed on a standard workstation in approximately 10s. A comparison with related methods is also provided.
Echo-planar imaging (EPI) is widely used in functional MRI studies. It is well known that EPI quality is usually degraded by geometric distortions, when there exist susceptibility field inhomogeneities. EPI distortions may be corrected if the field maps are available. It is possible to estimate the susceptibility field gradients from the phase reconstruction of a single-TE EPI image, after a successful phase-unwrapping procedure. However, in regions affected by pronounced field gradients, the phase-unwrapping of a single-TE image may fail, and therefore the estimated field maps may be incorrect. It has been reported that the field inhomogeneity may be calculated more reliably from T2*-weighted images corresponding to multiple TEs. However, the multi-TE MRI field mapping increases the scan time. Furthermore, the measured field maps may be invalid if the subject's position changes during dynamic scans. To overcome the limitations in conventional field mapping approaches, a novel k-space energy spectrum analysis algorithm is developed, which quantifies the spatially dependent echo-shifting effect and the susceptibility field gradients directly from the k-space data of single-TE gradient-echo EPI. Using the k-space energy spectrum analysis, susceptibility field gradients can be reliably measured without phase-unwrapping, and EPI distortions can be corrected without extra field mapping scans or pulse sequence modification. The reported technique can be used to retrospectively improve the image quality of the previously acquired EPI and functional MRI data, provided that the complex-domain k-space data are still available.
Detailed measurements of water diffusion within the prostate over an extended b-factor range were performed to assess whether the standard assumption of monoexponential signal decay is appropriate in this organ. From nine men undergoing prostate MR staging examinations at 1.5 T, a single 10-mm-thick axial slice was scanned with a line scan diffusion imaging sequence in which 14 equally spaced b factors from 5 to 3,500 s/mm(2) were sampled along three orthogonal diffusion sensitization directions in 6 min. Due to the combination of long scan time and limited volume coverage associated with the multi-b-factor, multidirectional sampling, the slice was chosen online from the available T2-weighted axial images with the specific goal of enabling the sampling of presumed noncancerous regions of interest (ROIs) within the central gland (CG) and peripheral zone (PZ). Histology from prescan biopsy (n=9) and postsurgical resection (n=4) was subsequently employed to help confirm that the ROIs sampled were noncancerous. The CG ROIs were characterized from the T2-weighted images as primarily mixtures of glandular and stromal benign prostatic hyperplasia, which is prevalent in this population. The water signal decays with b factor from all ROIs were clearly non-monoexponential and better served with bi- vs. monoexponential fits, as tested using chi(2)-based F test analyses. Fits to biexponential decay functions yielded intersubject fast diffusion component fractions in the order of 0.73+/-0.08 for both CG and PZ ROIs, fast diffusion coefficients of 2.68+/-0.39 and 2.52+/-0.38 microm(2)/ms and slow diffusion coefficients of 0.44+/-0.16 and 0.23+/-0.16 um(2)/ms for CG and PZ ROIs, respectively. The difference between the slow diffusion coefficients within CG and PZ was statistically significant as assessed with a Mann-Whitney nonparametric test (P<.05). We conclude that a monoexponential model for water diffusion decay in prostate tissue is inadequate when a large range of b factors is sampled and that biexponential analyses are better suited for characterizing prostate diffusion decay curves.
Previously, activation of vesicular transport in the brain microvasculature was shown to be one of the mechanisms of focused ultrasound-induced blood-brain barrier (BBB) opening. In the present study, we aimed to estimate the rate of the transendothelial vesicular traffic after focused ultrasound sonication in the rabbit brain, using ultrastructural morphometry and horseradish peroxidase (HRP) as a tracer. In the capillaries, the mean endothelial pinocytotic densities (the number of HRP-containing vesicles per microm(2) of the cell cytoplasm) were 0.9 and 1.05 vesicles/microm(2) 1 h after sonication with ultrasound frequencies of 0.69 and 0.26 MHz, respectively. In the arterioles, these densities were 1.63 and 2.43 vesicles/microm(2), values 1.8 and 2.3 times higher. In control locations, the densities were 0.7 and 0.14 vesicles/microm(2) for capillaries and arterioles, respectively. A small number of HRP-positive vesicles were observed in the venules. Focal delivery of HRP tracer was also observed in light microscopy. The results indicate that the precapillary microvessels play an important role in macromolecular transcytoplasmic traffic through the ultrasound-induced BBB modulation, which should be considered in the future development of trans-BBB drug delivery strategies.
Two image reconstruction methods currently dominate parallel MR imaging: SENSE and GRAPPA. While both seek to reconstruct images from subsampled multi-channel MRI data, there exist fundamental differences between the two. In particular, SENSE reconstructs an image of the excited spin-density directly whereas GRAPPA reconstructs estimates of the fully sampled raw coil data and then combines them to obtain an image. In this work we show that these differences can be exploited such that each method can compliment the other. In the case of SENSE, which requires an estimate of the coil sensitivity map before reconstruction, one can use GRAPPA to improve the coil sensitivity estimates. Alternatively, using coil sensitivity estimates and the SENSE reconstruction equations, one can improve the GRAPPA reconstruction parameter estimation. Together, these approaches can provide higher image quality than either method alone.
Unaliasing by Fourier-encoding the overlaps using the temporal dimension (UNFOLD) is a method to reduce the data acquisition burden in dynamic MRI. The method works by forcing aliased signals to behave in specific ways through time, so that these unwanted signals can be detected and removed. Unexpected events in time, such as displacements caused by breathing, have the potential to disturb the temporal strategy and may affect UNFOLD's ability to suppress aliasing artifacts. This work presents an extension of the UNFOLD method to accommodate temporal encoding disruptions. While the main type of disruption considered here comes from respiratory motion, other types of disruption can be envisioned, such as departures from the usual UNFOLD k-space sampling scheme. This extended version of UNFOLD was incorporated into UNFOLD-sensitivity encoding (UNFOLD-SENSE), and should also be applicable to closely related methods such as temporal SENSE (TSENSE), k-t Broaduse Linear Acquisition Speed up Technique (k-t BLAST), and k-t SENSE. Five patients were imaged with a modified version of a myocardial-perfusion sequence, and UNFOLD was used either alone or in conjunction with SENSE to obtain an acceleration of 2.0 (in three patients) or 3.0 (in two patients). In both cases this extended version of UNFOLD was able to suppress artifacts caused by the presence of breathing motion.
MR diffusion tensor imaging (DTI) of the brain and spine provides a unique tool for both visualizing directionality and assessing intactness of white matter fiber tracts in vivo. At the spatial resolution of clinical MRI, much of primate white matter is composed of interdigitating fibers. Analyses based on an assumed single diffusion tensor per voxel yield important information about the average diffusion in the voxel but fail to reveal structure in the presence of crossing tracts. Until today, all clinical scans assume only one tensor, causing potential serious errors in tractography. Since high angular resolution imaging remains, so far, untenable for routine clinical use, a method is proposed whereby the single-tensor field is augmented with additional information gleaned from standard clinical DTI. The method effectively resolves two distinct tract directions within voxels, in which only two tracts are assumed to exist. The underlying constrained two-tensor model is fitted in two stages, utilizing the information present in the single-tensor fit. As a result, the necessary MRI time can be drastically reduced when compared with other approaches, enabling widespread clinical use. Upon evaluation in simulations and application to in vivo human brain DTI data, the method appears to be robust and practical and, if correctly applied, could elucidate tract directions at critical points of uncertainty.
The segmentation of newborn brain MRI is important for assessing and directing treatment options for premature infants at risk for developmental disorders, abnormalities, or even death. Segmentation of infant brain MRI is particularly challenging when compared with the segmentation of images acquired from older children and adults. We sought to develop a fully automated segmentation strategy and present here a Bayesian approach utilizing an atlas of priors derived from previous segmentations and a new scheme for automatically selecting and iteratively refining classifier training data using the STAPLE algorithm. Results have been validated by comparison to hand-drawn segmentations.
Simon P Dimaio, Neculai Archip, Nobuhiko Hata, Ion-Florin Talos, Simon K Warfield, Amit Majumdar, Nathan McDannold, Kullervo Hynynen, Paul R Morrison, William M Wells, Daniel F Kacher, Randy E Ellis, Alexandra J Golby, Peter M Black, Ferenc A Jolesz, and Ron Kikinis. 2006. “Image-guided Neurosurgery at Brigham and Women’s Hospital.” IEEE Eng Med Biol Mag, 25, 5, Pp. 67-73.
We report a functional magnetic resonance imaging method to deliver task-specific brain activities as biofeedback signals to guide individuals to increase cortical activity in auditory areas during sound stimulation. A total of 11 study participants underwent multiple functional magnetic resonance imaging scan sessions, while the changes in the activated cortical volume within the primary and secondary auditory areas were fed back to them between scan sessions. On the basis of the feedback information, participants attempted to increase the number of significant voxels during the subsequent trial sessions by adjusting their level of attention to the auditory stimuli. Results showed that the group of individuals who received the feedback were able to increase the activation volume and blood oxygenation level-dependent signal to a greater degree than the control group.
Histologic effects of focused ultrasound (FUS) exposures combined with an ultrasound contrast agent (Optison) were investigated to examine whether the lesions were dominated by apoptosis or necrosis. The rabbit brains (n = 17) were sonicated (1.5 MHz, peak rarefactional pressure amplitude: 1.4 to 8.8 MPa) after Optison was injected intravenously (IV). MRI and light microscopy were used to examine tissue effects. To detect apoptosis, TUNEL staining based on labeling of DNA strand breaks was used. The average number of apoptotic and necrotic cells in 300 x 220 microm microscopic fields were counted in 18 representative lesions. Lesions in the rabbit brains were created at lowered acoustic power levels when FUS was combined with Optison. In histology, the lesions exhibited red blood cell extravasations and destruction of blood vessels. At 4 h after sonication, the lesions lost many cells, and the remaining cells exhibited both necrotic and apoptotic features. Overall, apoptosis dominated; there were, on average, 32.3 +/- 13.2 apoptotic cells per microscopic field compared with only 5.1 +/- 3.4 necrotic cells per field. In conclusion, FUS combined with Optison could produce lesions that are dominated by apoptosis, presumably induced primarily via ischemia after cavitation-produced damage to the brain vasculature.
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.