Development of a unique phantom to assess the geometric accuracy of magnetic resonance imaging for stereotactic localization

Phantoms for Quantitative Body MRI: a review and discussion of the phantom value

In this paper, we review the value of phantoms for body MRI in the context of their uses for quantitative MRI methods research, clinical trials, and clinical imaging. Certain uses of phantoms are common throughout the body MRI community, including measuring bias, assessing reproducibility, and training. In addition to these uses, phantoms in body MRI methods research are used for novel methods development and the design of motion compensation and mitigation techniques. For clinical trials, phantoms are an essential part of quality management strategies, facilitating the conduct of ethically sound, reliable, and regulatorily compliant clinical research of both novel MRI methods and therapeutic agents. In the clinic, phantoms are used for development of protocols, mitigation of cost, quality control, and radiotherapy. We briefly review phantoms developed for quantitative body MRI, and finally, we review open questions regarding the most effective use of a phantom for body MRI.

Assessment of Geometric Distortion in Six Clinical Scanners Using a 3D-Printed Grid Phantom

A cost-effective regularly structured three-dimensional (3D) printed grid phantom was developed to enable the quantification of machine-related magnetic resonance (MR) distortion. This phantom contains reference features, “point-like” objects, or vertices, which resulted from the intersection of mesh edges in 3D space. 3D distortions maps were computed by comparing the locations of corresponding features in both MR and computer tomography (CT) data sets using normalized cross correlation. Results are reported for six MRI scanners at both 1.5 T and 3.0 T field strengths within our institution. Mean Euclidean distance error for all MR volumes in this study, was less than 2 mm. The maximum detected error for the six scanners ranged from 2.4 mm to 6.9 mm. The conclusions in this study agree well with previous studies that indicated that MRI is quite accurate near the centre of the field but is more spatially inaccurate toward the edges of the magnetic field.

Nonlinear interactions in the thalamocortical loop in essential tremor: A model-based frequency domain analysis

There is increasing evidence to suggest that essential tremor has a central origin. Different structures appear to be part of the central tremorogenic network, including the motor cortex, the thalamus and the cerebellum. Some studies using electroencephalogram (EEG) and magnetoencephalography (MEG) show linear association in the tremor frequency between the motor cortex and the contralateral tremor electromyography (EMG). Additionally, high thalamomuscular coherence is found with the use of thalamic local field potential (LFP) recordings and tremulous EMG in patients undergoing surgery for deep brain stimulation (DBS). Despite a well-established reciprocal anatomical connection between the thalamus and cortex, the functional association between the two structures during "tremor-on" periods remains elusive. Thalamic (Vim) LFPs, ipsilateral scalp EEG from the sensorimotor cortex and contralateral tremor arm EMG recordings were obtained from two patients with essential tremor who had undergone successful surgery for DBS. Coherence analysis shows a strong linear association between thalamic LFPs and contralateral tremor EMG, but the relationship between the EEG and the thalamus is much less clear. These measurements were then analyzed by constructing a novel parametric nonlinear autoregressive with exogenous input (NARX) model. This new approach uncovered two distinct and not overlapping frequency "channels" of communication between Vim thalamus and the ipsilateral motor cortex, defining robustly "tremor-on" versus "tremor-off" states. The associated estimated nonlinear time lags also showed non-overlapping values between the two states, with longer corticothalamic lags (exceeding 50ms) in the tremor active state, suggesting involvement of an indirect multisynaptic loop. The results reveal the importance of the nonlinear interactions between cortical and subcortical areas in the central motor network of essential tremor. This work is important because it demonstrates for the first time that in essential tremor the functional interrelationships between the cortex and thalamus should not be sought exclusively within individual frequencies but more importantly between cross-frequency nonlinear interactions. Should our results be successfully reproduced on a bigger cohort of patients with essential tremor, our approach could be used to create an on-demand closed-loop DBS device, able to automatically activate when the tremor is on.

Continuous table acquisition MRI for radiotherapy treatment planning: distortion assessment with a new extended 3D volumetric phantom

PURPOSE

Accurate geometry is required for radiotherapy treatment planning (RTP). When considering the use of magnetic resonance imaging (MRI) for RTP, geometric distortions observed in the acquired images should be considered. While scanner technology and vendor supplied correction algorithms provide some correction, large distortions are still present in images, even when considering considerably smaller scan lengths than those typically acquired with CT in conventional RTP. This study investigates MRI acquisition with a moving table compared with static scans for potential geometric benefits for RTP.

METHODS

A full field of view (FOV) phantom (diameter 500 mm; length 513 mm) was developed for measuring geometric distortions in MR images over volumes pertinent to RTP. The phantom consisted of layers of refined plastic within which vitamin E capsules were inserted. The phantom was scanned on CT to provide the geometric gold standard and on MRI, with differences in capsule location determining the distortion. MRI images were acquired with two techniques. For the first method, standard static table acquisitions were considered. Both 2D and 3D acquisition techniques were investigated. With the second technique, images were acquired with a moving table. The same sequence was acquired with a static table and then with table speeds of 1.1 mm/s and 2 mm/s. All of the MR images acquired were registered to the CT dataset using a deformable B-spline registration with the resulting deformation fields providing the distortion information for each acquisition.

RESULTS

MR images acquired with the moving table enabled imaging of the whole phantom length while images acquired with a static table were only able to image 50%-70% of the phantom length of 513 mm. Maximum distortion values were reduced across a larger volume when imaging with a moving table. Increased table speed resulted in a larger contribution of distortion from gradient nonlinearities in the through-plane direction and an increased blurring of capsule images, resulting in an apparent capsule volume increase by up to 170% in extreme axial FOV regions. Blurring increased with table speed and in the central regions of the phantom, geometric distortion was less for static table acquisitions compared to a table speed of 2 mm/s over the same volume. Overall, the best geometric accuracy was achieved with a table speed of 1.1 mm/s.

CONCLUSIONS

The phantom designed enables full FOV imaging for distortion assessment for the purposes of RTP. MRI acquisition with a moving table extends the imaging volume in the z direction with reduced distortions which could be useful particularly if considering MR-only planning. If utilizing MR images to provide additional soft tissue information to the planning CT, standard acquisition sequences over a smaller volume would avoid introducing additional blurring or distortions from the through-plane table movement.

Clinical Experience with Pedunculopontine Nucleus Stimulation in Conditions with Nigrostriatal Disconnection

BACKGROUND

The pedunculopontine nucleus (PPN) is a part of the mesencephalic locomotor region and, in recent years, it has been considered a new surgical target for deep brain stimulation (DBS) for movement disorders including atypical parkinsonian syndromes such as progressive supranuclear palsy (PSP) and multiple system atrophy. Involvement of the PPN may play an important role in gait impairment in these disorders and the development of PPN DBS could potentially provide treatment for this disabling problem. However, the role of the PPN and the specific pathways involved in gait control and other motor functions are poorly understood.

METHODS

We present a chronological account of our group's experience in the use of PPN DBS. This entails the treatment of four patients with disabling movement disorders who all exhibited either marked damage or disconnection of the nigro-striatal pathway.

RESULTS

Within our series, the results were variable in that 2 of the 4 patients benefited greatly from DBS but the other 2 did not.

CONCLUSIONS

Our findings suggest that in carefully selected patients, PPN DBS can potentially alleviate symptoms due to dopaminergic striatal inactivity; symptoms that are typically resistant to stimulation of other subcortical targets used for parkinsonian syndromes and movement disorders.

Assessment and characterization of the total geometric uncertainty in Gamma Knife radiosurgery using polymer gels

PURPOSE

This work proposes and implements an experimental methodology, based on polymer gels, for assessing the total geometric uncertainty and characterizing its contributors in Gamma Knife (GK) radiosurgery.

METHODS

A treatment plan consisting of 26, 4-mm GK single shot dose distributions, covering an extended region of the Leksell stereotactic space, was prepared and delivered to a polymer gel filled polymethyl methacrylate (PMMA) head phantom (16 cm diameter) used to accurately reproduce every link in the GK treatment chain. The center of each shot served as a "control point" in the assessment of the GK total geometric uncertainty, which depends on (a) the spatial dose delivery uncertainty of the PERFEXION GK unit used in this work, (b) the spatial distortions inherent in MR images commonly used for target delineation, and (c) the geometric uncertainty contributor associated with the image registration procedure performed by the Leksell GammaPlan (LGP) treatment planning system (TPS), in the case that registration is directly based on the apparent fiducial locations depicted in each MR image by the N-shaped rods on the Leksell localization box. The irradiated phantom was MR imaged at 1.5 T employing a T2-weighted pulse sequence. Four image series were acquired by alternating the frequency encoding axis and reversing the read gradient polarity, thus allowing the characterization of the MR-related spatial distortions.

RESULTS

MR spatial distortions stemming from main field (B0) inhomogeneity as well as from susceptibility and chemical shift phenomena (also known as sequence dependent distortions) were found to be of the order of 0.5 mm, while those owing to gradient nonlinearities (also known as sequence independent distortions) were found to increase with distance from the MR scanner isocenter extending up to 0.47 mm at an Euclidean distance of 69.6 mm. Regarding the LGP image registration procedure, the corresponding average contribution to the total geometric uncertainty ranged from 0.34 to 0.80 mm. The average total geometric uncertainty, which also includes the GK spatial dose delivery uncertainty, was found equal to (0.88 ± 0.16), (0.88 ± 0.26), (1.02 ± 0.09), and (1.15 ± 0.24) mm for the MR image series acquired with the read gradient polarity (direction) set toward right, left, posterior, and anterior, respectively.

CONCLUSIONS

The implemented methodology seems capable of assessing the total geometric uncertainty, as well as of characterizing its contributors, ascribed to the entire GK treatment delivery (i.e., from MR imaging to GK dose delivery) for an extended region of the Leksell stereotactic space. Results obtained indicate that the selection of both the frequency encoding axis and the read gradient polarity during MRI acquisition may affect the magnitude as well as the spatial components of the total geometric uncertainty.

Human periventricular grey somatosensory evoked potentials suggest rostrocaudally inverted somatotopy

BACKGROUND

Somatosensory homunculi have been demonstrated in primary somatosensory cortex and ventral posterior thalamus but not periaqueductal and periventricular grey matter (PAVG), a therapeutic target for deep brain stimulation (DBS) in chronic pain.

AIMS

The study is an investigation of somatotopic representation in PAVG and assessment for a somatosensory homunculus.

METHODS

Five human subjects were investigated using electrical somatosensory stimulation and deep brain macroelectrode recording. DBS were implanted in the contralateral PAVG. Cutaneous arm, leg and face regions were stimulated while event-related potentials were recorded from deep brain electrodes. Electrode contact positions were mapped using MRI and brain atlas information.

RESULTS

Monopolar P1 somatosensory evoked potential amplitudes were highest and onset latencies shortest in contralateral caudal PAVG with facial stimulation and rostral with leg stimulation, in agreement with reported subjective sensation during intra-operative electrode advancement.

CONCLUSIONS

A rostrocaudally inverted somatosensory homunculus exists in the human PAVG region. Objective human evidence of PAVG somatotopy increases understanding of a brainstem region important to pain and autonomic control that is a clinical target for both pharmacological and neurosurgical therapies. Such knowledge may assist DBS target localisation for neuropathic pain syndromes related to particular body regions like brachial plexopathies, anaesthesia dolorosa and phantom limb pain.

Deep brain stimulation for pain

Deep brain stimulation (DBS) is a neurosurgical intervention whose efficacy, safety, and utility have been shown in the treatment of movement disorders. For the treatment of chronic pain refractory to medical therapies, many prospective case series have been reported, but few have published findings from patients treated during the past decade using current standards of neuroimaging and stimulator technology. We summarize the history, science, selection, assessment, surgery, and personal clinical experience of DBS of the ventral posterior thalamus, periventricular/periaqueductal gray matter, and, latterly, the rostral anterior cingulate cortex (Cg24) in 100 patients treated now at two centers (John Radcliffe Hospital, Oxford, UK, and Hospital de São João, Porto, Portugal) over 12 years. Several experienced centers continue DBS for chronic pain with success in selected patients, in particular those with pain after amputation, brachial plexus injury, stroke, and cephalalgias including anesthesia dolorosa. Other successes include pain after multiple sclerosis and spine injury. Somatotopic coverage during awake surgery is important in our technique, with cingulate DBS considered for whole-body pain or after unsuccessful DBS of other targets. Findings discussed from neuroimaging modalities, invasive neurophysiological insights from local field potential recording, and autonomic assessments may translate into improved patient selection and enhanced efficacy, encouraging larger clinical trials.

Measurement of MRI scanner performance with the ADNI phantom

The objectives of this study are as follows: to describe practical implementation challenges of multisite, multivendor quantitative studies; to describe the MRI phantom and analysis software used in the Alzheimer's Disease Neuroimaging Initiative (ADNI) study, illustrate the utility of the system for measuring scanner performance, the ability to assess gradient field nonlinearity corrections: and to recover human brain images without geometric scaling errors in multisite studies. ADNI is a large multicenter study with each center having its own copy of the phantom. The design of the phantom and analysis software are presented as results from predistribution systematics studies and results from field experience with the phantom at 58 enrolling ADNI sites over a 3 year period. The estimated coefficients of variation intrinsic to measurements of geometry in a single phantom are in the range of 3-5 parts in 10(4). Phantom measurements accurately detect linear and nonlinear scaling in images. Gradient unwarping methods are readily assessed by phantom nonlinearity measurements. Phantom-based scaling correction reduces observed geometric drift in human images by one-third or more. Repair or replacement of phantoms between scans, however, is a confounding factor. The ADNI phantom can be used to assess both scanner performance and the validity of postprocessing image corrections in order to reduce systematic errors in human images. Reduced measurement errors should decrease measurement bias and increase statistical power for measurements of rates of change in the brain structure in AD treatment trials. Perhaps the greatest practical value of incorporating ADNI phantom measurements in a multisite study is to identify scanner errors through central monitoring. This approach has resulted in identification of system errors including sites misidentification of their own gradient hardware and the disabling of autoshim, and a miscalibrated laser alignment light. If undetected, these errors would have contributed to imprecision in quantitative metrics at over 25% of all enrolling ADNI sites.

Avoiding the ventricle: a simple step to improve accuracy of anatomical targeting during deep brain stimulation

OBJECT

The authors examined the accuracy of anatomical targeting during electrode implantation for deep brain stimulation in functional neurosurgical procedures. Special attention was focused on the impact that ventricular involvement of the electrode trajectory had on targeting accuracy.

METHODS

The targeting error during electrode placement was assessed in 162 electrodes implanted in 109 patients at 2 centers. The targeting error was calculated as the shortest distance from the intended stereotactic coordinates to the final electrode trajectory as defined on postoperative stereotactic imaging. The trajectory of these electrodes in relation to the lateral ventricles was also analyzed on postoperative images.

RESULTS

The trajectory of 68 electrodes involved the ventricle. The targeting error for all electrodes was calculated: the mean +/- SD and the 95% CI of the mean was 1.5 +/- 1.0 and 0.1 mm, respectively. The same calculations for targeting error for electrode trajectories that did not involve the ventricle were 1.2 +/- 0.7 and 0.1 mm. A significantly larger targeting error was seen in trajectories that involved the ventricle (1.9 +/- 1.1 and 0.3 mm; p < 0.001). Thirty electrodes (19%) required multiple passes before final electrode implantation on the basis of physiological and/or clinical observations. There was a significant association between an increased requirement for multiple brain passes and ventricular involvement in the trajectory (p < 0.01).

CONCLUSIONS

Planning an electrode trajectory that avoids the ventricles is a simple precaution that significantly improves the accuracy of anatomical targeting during electrode placement for deep brain stimulation. Avoidance of the ventricles appears to reduce the need for multiple passes through the brain to reach the desired target as defined by clinical and physiological observations.

Design, construction and evaluation of an anthropomorphic head phantom with realistic susceptibility artifacts

PURPOSE

To design and construct an anthropomorphic head phantom using materials of appropriate magnetic susceptibility and air spaces of realistic dimensions, with the aim of reproducing the susceptibility artifacts found in the human brain.

MATERIALS AND METHODS

The phantom is based on a plastic skull filled with MnCl2-doped water. Materials to mimic soft tissue (wax) and bone (plastic skull) were chosen based on mass susceptibility measurements using a superconducting quantum interference device (SQUID) magnetometer. The phantom was designed for and evaluated at 4.7T using field mapping and echo-planar imaging (EPI).

RESULTS

The main magnetic field (B0) maps of the phantom resemble those of four volunteers' brains and have similar standard deviations (SDs). Maps of the B0 field gradients in the phantom and real brains are also similar. The phantom has relaxation times close to those of brain tissue at 4.7T. Gradient-echo (GE)-EPI images of the phantom suffer from susceptibility artifacts comparable to those in real heads and at anatomically realistic locations.

CONCLUSION

The phantom is a useful tool for evaluating and comparing different susceptibility artifact reduction techniques. The phantom could also be used to test CT-MRI coregistration in the presence of susceptibility artifacts since the water-filled brain cavity is both CT- and MR-visible.

Regional cerebral perfusion differences between periventricular grey, thalamic and dual target deep brain stimulation for chronic neuropathic pain

Regional cerebral blood flow changes were evaluated in different subcortical brain targets following deep brain stimulation (DBS) for chronic pain. Three patients with intractable neuropathic pain were assessed; one had stimulating electrodes in the ventroposterolateral thalamic nucleus (VPL), one in the periventricular grey (PVG) area, and one had electrodes in both targets. Pain relief was achieved in all patients. Cerebral perfusion was measured by single-photon emission computed tomography to determine the effects of DBS. Comparison was made between individual scans using subtraction analysis. DBS consistently increased perfusion in the posterior subcortical region between VPL and PVG, regardless of the site of stimulation. Furthermore, thalamic and dual target DBS increased thalamic perfusion, yet PVG DBS decreased perfusion in the PVG-containing midbrain region and thalamus. Dual target stimulation decreased anterior cingulate and insular cortex perfusion. The study demonstrates regional differences in cerebral perfusion between three accepted and efficacious targets for analgesic DBS.

Autonomic neurosurgery: from microvascular decompression to image guided stimulation

The paper reviews mechanisms underlying autonomic disorders, with a focus on cardiovascular dysfunction. Neurosurgical approaches are described for medically refractory hypertension and orthostatic hypotension. After review of microvascular decompression of the rostral ventrolateral medulla, stereotactic CT and MRI guided deep brain stimulation of the periaqueductal grey matter (PAG) is evaluated. Results are presented from patient studies showing reductions in blood pressure with ventral PAG stimulation and increases in blood pressure with dorsal PAG stimulation. A rationale for the treatment of autonomic disorders by neurosurgical intervention is discussed.

Neuronavigation and surgery of intracerebral tumours

Approximately four decades after the successful clinical introduction of framebased stereotactic neurosurgery by Spiegel and Wycis, frameless stereotaxy emerged to enable more elaborate image guidance in open neurosurgical procedures. Frameless stereotaxy, or neuronavigation, relies on one of several different localizing techniques to determine the position of an operative instrument relative to the surgical field, without the need for a coordinate frame rigidly fixed to the patients' skull. Currently, most systems are based on the optical triangulation of infrared light sources fixed to the surgical instrument. In its essence, a navigation system is a three-dimensional digitiser that correlates its measurements to a reference data set, i.e. a preoperatively acquired CT or MRI image stack. This correlation is achieved through a patient-to-image registration procedure resulting in a mathematical transformation matrix mapping each position in 'world space' onto 'image space'. Thus, throughout the remainder of the surgical procedure, the position of the surgical instrument can be demonstrated on a computer screen, relative to the CT or MRI images. Though neuronavigation has become a routinely used addition to the neurosurgical armamentarium, its impact on surgical results has not yet been examined sufficiently. Therefore, the surgeon is left to decide on a case-by-case basis whether to perform surgery with or without neuronavigation. Future challenges lie in improvement of the interface between the surgeon and the neuronavigator and in reducing the brainshift error, i.e. inaccuracy introduced by changes in tissue positions after image acquisition.

Effect of GPi DBS on functional imaging of the brain in dystonia

Five patients with idiopathic dystonic conditions, treated successfully with deep brain stimulation (DBS) of the globus pallidus internus (GPi), were studied using single-photon emission tomography (SPET) in order to evaluate brain perfusion in the presence and absence of DBS. Comparison was made between the "on" and "off" DBS scans on an individual basis and also as part of a group analysis. Whilst the individual data suggested great regional variation in cerebral perfusion between individuals, the results of the group analysis revealed several topographically similar areas of the brain where relative hyperperfusion in the absence of DBS was common to all patients. Based on these results we postulate on possible mechanisms for this phenomenon.

Distortion of Magnetic Resonance Images Used in Gamma Knife Radiosurgery Treatment Planning: Implications for Acoustic Neuroma Outcomes

OBJECTIVE

To quantify the image distortion of our series of acoustic neuromas treated with gamma knife radiosurgery.

STUDY DESIGN

Retrospective chart and digital radiographic file review with quantitative assessment of gamma knife treatment plans.

SETTING

Tertiary referral center.

PATIENTS

Patients undergoing gamma knife radiosurgery for the treatment of acoustic neuromas.

INTERVENTION

Gamma knife radiosurgery.

MAIN OUTCOME MEASURES

Gamma knife treatment plans containing magnetic resonance images were reviewed at each axial, sagittal, and coronal slice. The length of the greatest displacement of the treatment plan was measured and the volume of the treatment plan that fell outside of the internal auditory canal calculated. Known clinical measurements of audiometric, vestibular, facial, and trigeminal nerve functions were then compared with current measurements of tumor size.

RESULTS

Twenty-two of the 23 patients had measurable image shifts on the axial images. The range of the image shift was 0 to 5.8 mm, with a mean shift of 1.92 +/- 1.29 mm (+/- standard deviation). Tumor volumes of the treatment plan that fell outside of the internal auditory canal ranged from 0 to 414 mm, with a mean of 90.5 mm. The mean percentage that fell outside of the internal auditory canal was 16.7% of total tumor volume (range, 2.4-77.6%). We could not draw any consistent correlations between degree of image shift and continued tumor growth or objective examination values.

CONCLUSION

We have demonstrated a small but potentially significant shift in the treatment plan of gamma knife radiosurgery when based on magnetic resonance images. Although the image shift does not seem to affect the growth of the acoustic neuromas or auditory or facial nerve function, longer term follow-up is required to fully appreciate the true impact of this image shift.

Geometric distortion in clinical MRI systems Part I: evaluation using a 3D phantom

Recently, a 3D phantom that can provide a comprehensive and accurate measurement of the geometric distortion in MRI has been developed. Using this phantom, a full assessment of the geometric distortion in a number of clinical MRI systems (GE and Siemens) has been carried out and detailed results are presented in this paper. As expected, the main source of geometric distortion in modern superconducting MRI systems arises from the gradient field nonlinearity. Significantly large distortions with maximum absolute geometric errors ranged between 10 and 25 mm within a volume of 240 x 240 x 240 mm(3) were observed when imaging with the new generation of gradient systems that employs shorter coils. By comparison, the geometric distortion was much less in the older-generation gradient systems. With the vendor's correction method, the geometric distortion measured was significantly reduced but only within the plane in which these 2D correction methods were applied. Distortion along the axis normal to the plane was, as expected, virtually unchanged. Two-dimensional correction methods are a convenient approach and in principle they are the only methods that can be applied to correct geometric distortion in a single slice or in multiple noncontiguous slices. However, these methods only provide an incomplete solution to the problem and their value can be significantly reduced if the distortion along the normal of the correction plane is not small.

Three-dimensional dose verification of the clinical application of gamma knife stereotactic radiosurgery using polymer gel and MRI

This work seeks to verify multi-shot clinical applications of stereotactic radiosurgery with a Leksell Gamma Knife model C unit employing a polymer gel-MRI based experimental procedure, which has already been shown to be capable of verifying the precision and accuracy of dose delivery in single-shot gamma knife applications. The treatment plan studied in the present work resembles a clinical treatment case of pituitary adenoma using four 8 mm and one 14 mm collimator helmet shots to deliver a prescription dose of 15 Gy to the 50% isodose line (30 Gy maximum dose). For the experimental dose verification of the treatment plan, the same criteria as those used in the clinical treatment planning evaluation were employed. These included comparison of measured and GammaPlan calculated data, in terms of percentage isodose contours on axial, coronal and sagittal planes, as well as 3D plan evaluation criteria such as dose-volume histograms for the target volume, target coverage and conformity indices. Measured percentage isodose contours compared favourably with calculated ones despite individual point fluctuations at low dose contours (e.g., 20%) mainly due to the effect of T2 measurement uncertainty on dose resolution. Dose-volume histogram data were also found in a good agreement while the experimental results for the percentage target coverage and conformity index were 94% and 1.17 relative to corresponding GammaPlan calculations of 96% and 1.12, respectively. Overall, polymer gel results verified the planned dose distribution within experimental uncertainties and uncertainty related to the digitization process of selected GammaPlan output data.

Dose verification of single shot gamma knife applications using VIPAR polymer gel and MRI

This work describes an experimental procedure with potential to assess the overall accuracy associated with gamma knife clinical applications, from patient imaging and dosimetry planning to patient positioning and dose delivery using the automated positioning system of a Leksell Gamma Knife model C. The VIPAR polymer gel-MRI dosimetry method is employed due to its inherent three-dimensional feature and linear dose response over the range of gamma knife applications. Different polymer gel vials were irradiated with single shot gamma knife treatment plans using each of the four available collimator helmets to deliver a maximum dose of 30 Gy. Percentage relative dose results are presented not only in the form of one-dimensional profiles but also planar isocontours and isosurfaces in three dimensions. Experimental results are compared with corresponding Gammaplan treatment planning system calculations as well as acceptance test radiochromic film measurements. A good agreement, within the experimental uncertainty, is observed between measured and expected dose distributions. This experimental uncertainty is of the order of one imaging pixel in the MRI gel readout session (<1 mm) and allows for the verification of single shot gamma knife applications in terms of acceptance specifications for precision in beam alignment and accuracy. Averaging net R(2) results in the dose plateau of the 4 mm and 18 mm collimator irradiated gel vials, which were MR scanned in the same session, provides a crude estimate of the 4 mm output factor which agrees within errors with the default value of 0.870.

Correspondence of microelectrode mapping with magnetic resonance imaging for subthalamic nucleus procedures

BACKGROUND

Magnetic resonance imaging (MRI) and microelectrode recording (MER) are commonly used to guide stereotactic procedures on the subthalamic nucleus (STN). Little is known about the correlation between the position of the STN as seen on MRI and that as determined by MER mapping. We compared these in 10 patients with Parkinson's disease.

METHODS

The position of the STN was determined by intraoperative MER findings and stereotactic axial T2 magnetic resonance images with 2-mm slice thickness. Images were reconstructed in a 3-dimensional workstation. The anterior, posterior, medial, lateral, dorsal, and ventral borders of the STN defined with the MRI were measured relative to the midcommissural point. The location of STN activity during MER was reconstructed relative to the midcommissural point for comparison.

RESULTS

Twenty-nine tracks recorded with microelectrodes provided clear spans of STN-like activity in 18 STN nuclei. The coordinates of MER were, in general, within the borders of the STN defined with the MRI. However, when analyzed individually, some of the tracks had STN-like activity outside the borders of the MRI-defined nucleus (mostly <1 mm). Three tracks had STN-like activity recorded between 2 and 3 mm more anterior than the anterior border of the nucleus defined with the MRI.

CONCLUSIONS

There was a good correlation between MER and the borders of the STN defined in the MRI, except for the anterior-posterior axis, in which MER indicated that the STN extended more anteriorly than as suggested by MRI. This should be taken into account in STN surgery.

A proposed scheme for comprehensive characterization of the measured geometric distortion in magnetic resonance imaging using a three-dimensional phantom

Recently, a 3-dimensional phantom that can provide a comprehensive, accurate and complete measurement of the geometric distortion in MRI has been developed. In this paper, a scheme for characterizing the measured geometric distortion using the 3-D phantom is described. In the proposed scheme, a number of quantitative measures are developed and used to characterize the geometric distortion. These measures encompass the overall and spatial aspects of the geometric distortion. Two specific types of volume of interest, rectangular parallelepipeds (including cubes) and spheres are considered in the proposed scheme. As an illustration, characterization of the geometric distortion in a Siemens 1.5T Sonata MRI system using the proposed scheme is presented. As shown, the proposed scheme provides a comprehensive assessment of the geometric distortion. The scheme can be potentially used as a standard procedure for the assessment of geometric distortion in MRI.

Analyzing 3-tesla magnetic resonance imaging units for implementation in radiosurgery

OBJECT

The limiting factor affecting accuracy during gamma knife surgery is image quality. The new generation of magnetic resonance (MR) imaging units with field strength up to 3 teslas promise superior image quality for anatomical resolution and contrast. There are, however, questions about chemical shifts or susceptibility effects, which are the subject of this paper.

METHODS

The 3-tesla MR imaging unit (Siemens Trio) was analyzed and compared with a 1-tesla unit (Siemens Magnetom Expert) and to a 1.5-tesla unit (Philips Gyroscan). Evaluation of the magnitude of error was performed within transverse slices in two orientations (axial/coronal) by using a cylindrical phantom with an embedded grid. Deviations were determined for 21 targets in a slab phantom with known geometrical positions within the stereotactic frame. Distortions caused by chemical shift and/or susceptibility effects were analyzed in a head phantom. Inhouse software was used for data analyses. The mean deviation was less than 0.3 mm in axial and less than 0.4 mm in coronal orientations. For the known targets the maximum deviation was 1.16 mm. By optimizing these parameters in the protocol these inaccuracies could be reduced to less than 1.1 mm. Due to inhomogeneities a shift in the z direction of up to 1.5 mm was observed for a dataset, which was shown to be compressed by 1.2 mm.

CONCLUSIONS

The 3-tesla imaging unit showed superior anatomical contrast and resolution in comparison with the established 1-tesla and 1.5-tesla units; however, due to the high field strength the field within the head coil is very sensitive to inhomogeneities and therefore 3-tesla imaging data will have be handled with care.

Does new magnetic resonance imaging technology provide better geometrical accuracy during stereotactic imaging?

Object.The authors sought to compare the accuracy of stereotactic target imaging using the Siemens 1T EXPERT and 1.5T SYMPHONY magnetic resonance (MR) units.

Methods.A water-filled cylindrical Perspex phantom with axial and coronal inserts containing grids of glass rods was fixed in the Leksell stereotactic frame and subjected to MR imaging in Siemens 1T EXPERT and Siemens 1.5T SYMPHONY units. Identical sequences were used for each unit. The images were transferred to the GammaPlan treatment planning system. Deviations between stereotactic coordinates based on MR images and estimated real geometrical positions given by the construction of the phantom insert were evaluated for each study. The deviations were further investigated as a function of the MR unit used, MR sequence, the image orientation, and the spatial position of measured points in the investigated volume.

Conclusions.Larger distortions were observed when using the SYMPHONY 1.5T unit than those with the EXPERT 1T unit. Typical average distortion in EXPERT 1T was not more than 0.6 mm and 0.9 mm for axial and coronal images, respectively. Typical mean distortion for SYMPHONY 1.5T was not more than 1 mm and 1.3 mm for axial and coronal images, respectively. The image sequence affected the distortions in both units. Coronal T2-weighted spin-echo images performed in subthalamic imaging produced the largest distortions of 2.6 mm and 3 mm in the EXPERT 1T and SYMPHONY 1.5T, respectively. Larger distortions were observed in coronal slices than in axial slices in both units, and this effect was more pronounced in SYMPHONY 1.5T. Noncentrally located slice positions in the investigated volume of the phantom were associated with larger distortions.

Does new magnetic resonance imaging technology provide better geometrical accuracy during stereotactic imaging?

Object. The authors sought to compare the accuracy of stereotactic target imaging using the Siemens 1T EXPERT and 1.5T SYMPHONY magnetic resonance (MR) units.

Methods. A water-filled cylindrical Perspex phantom with axial and coronal inserts containing grids of glass rods was fixed in the Leksell stereotactic frame and subjected to MR imaging in Siemens 1T EXPERT and Siemens 1.5T SYMPHONY units. Identical sequences were used for each unit. The images were transferred to the GammaPlan treatment planning system. Deviations between stereotactic coordinates based on MR images and estimated real geometrical positions given by the construction of the phantom insert were evaluated for each study. The deviations were further investigated as a function of the MR unit used, MR sequence, the image orientation, and the spatial position of measured points in the investigated volume.

Conclusions. Larger distortions were observed when using the SYMPHONY 1.5T unit than those with the EXPERT 1T unit. Typical average distortion in EXPERT 1T was not more than 0.6 mm and 0.9 mm for axial and coronal images, respectively. Typical mean distortion for SYMPHONY 1.5T was not more than 1 mm and 1.3 mm for axial and coronal images, respectively. The image sequence affected the distortions in both units. Coronal T2-weighted spin-echo images performed in subthalamic imaging produced the largest distortions of 2.6 mm and 3 mm in the EXPERT 1T and SYMPHONY 1.5T, respectively. Larger distortions were observed in coronal slices than in axial slices in both units, and this effect was more pronounced in SYMPHONY 1.5T. Noncentrally located slice positions in the investigated volume of the phantom were associated with larger distortions.

Analyzing 3-tesla magnetic resonance imaging units for implementation in radiosurgery

Object. The limiting factor affecting accuracy during gamma knife surgery is image quality. The new generation of magnetic resonance (MR) imaging units with field strength up to 3 teslas promise superior image quality for anatomical resolution and contrast. There are, however, questions about chemical shifts or susceptibility effects, which are the subject of this paper.

Methods. The 3-tesla MR imaging unit (Siemens Trio) was analyzed and compared with a 1-tesla unit (Siemens Magnetom Expert) and to a 1.5-tesla unit (Philips Gyroscan). Evaluation of the magnitude of error was performed within transverse slices in two orientations (axial/coronal) by using a cylindrical phantom with an embedded grid. Deviations were determined for 21 targets in a slab phantom with known geometrical positions within the stereotactic frame. Distortions caused by chemical shift and/or susceptibility effects were analyzed in a head phantom. Inhouse software was used for data analyses.

The mean deviation was less than 0.3 mm in axial and less than 0.4 mm in coronal orientations. For the known targets the maximum deviation was 1.16 mm. By optimizing these parameters in the protocol these inaccuracies could be reduced to less than 1.1 mm. Due to inhomogeneities a shift in the z direction of up to 1.5 mm was observed for a dataset, which was shown to be compressed by 1.2 mm.

Conclusions. The 3-tesla imaging unit showed superior anatomical contrast and resolution in comparison with the established 1-tesla and 1.5-tesla units; however, due to the high field strength the field within the head coil is very sensitive to inhomogeneities and therefore 3-tesla imaging data will have be handled with care.

Senile chorea treated by deep brain stimulation: a clinical, neurophysiological and functional imaging study

We report on a patient with senile chorea, treated with deep brain stimulation of the left globus pallidus internus and subsequently the left ventralis oralis posterior nucleus of the thalamus. Deep brain field potential recordings and functional imaging using single photon emission tomography enabled us to suggest pathophysiological mechanisms for the symptoms.

A realistic brain tissue phantom for intraparenchymal infusion studies

Object. The goal of this study was to validate a simple, inexpensive, and robust model system to be used as an in vitro surrogate for in vivo brain tissues in preclinical and exploratory studies of infusion-based intraparenchymal drug and cell delivery.

Methods. Agarose gels of varying concentrations and porcine brain were tested to determine the infusion characteristics of several different catheters at flow rates of 0.5 and 1 µl per minute by using bromophenol blue (BPB) dye (molecular weight [MW] ∼690) and gadodiamide (MW ∼573). Magnetic resonance (MR) imaging and videomicroscopy were used to measure the distribution of these infusates, with a simultaneous measurement of infusion pressures. In addition, the forces of catheter penetration and movement through gel and brain were measured.

Agarose gel at a 0.6% concentration closely resembles in vivo brain with respect to several critical physical characteristics. The ratio of distribution volume to infusion volume of agarose was 10 compared with 7.1 for brain. The infusion pressure of the gel demonstrated profiles similar in configuration and magnitude to those of the brain (plateau pressures 10–20 mm Hg). Gadodiamide infusion in agarose closely resembled that in the brain, as documented using T1-weighted MR imaging. Gadodiamide distribution in agarose gel was virtually identical to that of BPB dye, as documented by MR imaging and videomicroscopy. The force profile for insertion of a silastic catheter into agarose gel was similar in magnitude and configuration to the force profile for insertion into the brain. Careful insertion of the cannula using a stereotactic guide is critical to minimize irregularity and backflow of infusate distribution.

Conclusions. Agarose gel (0.6%) is a useful surrogate for in vivo brain in exploratory studies of convection-enhanced delivery.

A novel phantom and method for comprehensive 3-dimensional measurement and correction of geometric distortion in magnetic resonance imaging

A phantom that can be used for mapping geometric distortion in magnetic resonance imaging (MRI) is described. This phantom provides an array of densely distributed control points in three-dimensional (3D) space. These points form the basis of a comprehensive measurement method to correct for geometric distortion in MR images arising principally from gradient field non-linearity and magnet field inhomogeneity. The phantom was designed based on the concept that a point in space can be defined using three orthogonal planes. This novel design approach allows for as many control points as desired. Employing this novel design, a highly accurate method has been developed that enables the positions of the control points to be measured to sub-voxel accuracy. The phantom described in this paper was constructed to fit into a body coil of a MRI scanner, (external dimensions of the phantom were: 310 mm x 310 mm x 310 mm), and it contained 10,830 control points. With this phantom, the mean errors in the measured coordinates of the control points were on the order of 0.1 mm or less, which were less than one tenth of the voxel's dimensions of the phantom image. The calculated three-dimensional distortion map, i.e., the differences between the image positions and true positions of the control points, can then be used to compensate for geometric distortion for a full image restoration. It is anticipated that this novel method will have an impact on the applicability of MRI in both clinical and research settings, especially in areas where geometric accuracy is highly required, such as in MR neuro-imaging.

Post-operative progress of dystonia patients following globus pallidus internus deep brain stimulation

In the current era of functional surgery for movement disorders, deep brain stimulation (DBS) of the globus pallidus internus (GPi) is emerging as the favoured intervention for patients with dystonia. Here we report our results in 20 patients with medically intractable dystonia treated with GPi stimulation. The series comprised 14 patients with generalized dystonia and six with spasmodic torticollis. Although comparisons were limited by differences in their respective neurological rating scales, chronic DBS clearly benefited both patient groups. Data conveying the rate of change in neurological function following intervention are also presented, demonstrating the gradual but progressive and sustained nature of improvement following stimulation of the GPi in dystonic patients.

Application of polymer gel dosimetry in gamma knife radiosurgery

Object. The purpose of this study was to investigate the use of a polymer gel—based dosimeter for the evaluation of geometric and dosimetric inaccuracies during gamma knife radiosurgery and during the irradiation of an experimental animal.

Methods. A polymer gel dosimeter, based on acrylic monomers, was used for experiments conducted in this study. The accuracy of the dosimeter was evaluated on a Siemens EXPERT 1-tesla scanner in the transmitter/receiver head coil with the use of a multiecho sequence with 16 echoes, TE 22.5 to 360 msec, TR 2000 msec, slice thickness 2 mm, field of view 255 mm, and a pixel size of 0.5 × 0.5 mm2. Two experiments were conducted. First, the head phantom containing the polymer gel dosimeter was irradiated using 4-, 8-, 14-, and 18-mm isocenters. Second, a specially designed rat phantom was irradiated by four 4-mm isocenters. The dose profiles in the x, y, and z axes were calculated in the treatment planning system and measured with the polymer gel dosimeter and the results were compared.

There was good agreement between the measured and calculated dose profiles. The maximum deviation in the spatial position of the center of measured and calculated dose profiles was 0.5 mm in the head phantom and 1 mm in the rat phantom. The maximum deviation in the width of the selected reference isodose of measured profiles was 1.2 mm in the head phantom and 1.1 mm in the rat phantom.

Conclusions. The use of the polymer gel—based dosimeter for the verification of stereotactic procedures has advantages compared with other dosimetric systems. The dosimeter itself is tissue equivalent. Three-dimensional dose distributions can be measured and the dosimeter allows simulation of the therapeutic procedures.

Effect of subpixel magnetic resonance imaging shifts on radiosurgical dosimetry for vestibular schwannoma

Object. The purpose of this study was to evaluate subpixel magnetic resonance (MR) imaging shifts of intracanalicular vestibular schwannomas (VSs) with respect to the internal auditory canal (IAC) as documented on computerized tomography (CT) scanning and to investigate the source of imaging-related localization errors in radiosurgery as well as the effect of such shifts on the dosimetry for small targets.

Methods. A shift of the stereotactic coordinates of intracanalicular VSs between those determined on MR imaging and those on CT scanning represents an error in localization. A shift vector places the tumor within the IAC and measures the CT scan/MR image discrepancy. The shift vectors were measured in a series of 15 largely intracanalicular VSs (all < 1.5 cm3 in volume). Using dose volume histogram measurements, the overlap between shifted and unshifted tumors and radiosurgical treatment plans were measured. Using plastic and bone phantoms and thermoluminescent dosimetry measurements, the correspondence between CT and MR imaging targets and treatments delivered using the Leksell gamma knife were measured. Combining these measurements, the correspondence between intended and actual treatments was measured.

Conclusions. The delivery of radiation to CT-imaged targets was accurate to the limits of measurement (∼ 0.1 mm). The MR imaging shifts seen in the y axis averaged 0.9 mm and in the z axis 0.8 mm. The corresponding percentage of tumor coverage with respect to apparent target shift decreased from 98 to 77%. This represents a significant potential error when targets are defined solely by MR imaging.

Modification of reference array attachment for image guided neurosurgery: direct cranial fixation

We have modified the cranial fixation technique of the reference array used for the Stealth (Medtronics Inc., Minneapolis, MN) image guided neurosurgical workstation to avoid rigid immobilization and to accommodate patients undergoing awake procedures. The modification allows attachment of a reference array directly to the skull prior to registration, avoiding the requirement for rigid cranial fixation. The accuracy of fiducial registration for the modified reference array was compared to the conventional reference array using a phantom system yielding similar registration results and target accuracy. The successful application of the modified system to two operative cases is described.

Quality control of the stereotactic radiosurgery procedure with the polymer-gel dosimetry

PURPOSE

To assess the entire geometric and dosimetric (relative) uncertainties of the radiosurgery procedure with the Leksell gamma knife.

MATERIALS AND METHODS

The entire Leksell gamma knife stereotactic radiosurgery treatment procedure was simulated with the use of a special water filled head phantom and polymer-gel dosimeter evaluated by nuclear magnetic resonance (NMR). A test vessel filled with the polymer-gel dosimeter was fixed in the head phantom. The phantom underwent stereotactic NMR imaging, treatment planning and then irradiation according to the treatment plan prepared exactly the same way as in the ordinary treatment procedure for a patient. The treatment plan was represented by one isocenter positioned approximately centrally in the head phantom. This procedure was subsequently repeated for all four collimators (4, 8, 14, 18mm) used on the Leksell gamma knife. Evaluation of dosimeters was performed on a Siemens EXPERT 1T NMR scanner. Dose profiles in X, Y and Z axes through the ellipsoidal shaped dose distribution were obtained to compare experimental results from the irradiated phantom with the treatment planning system calculations.

RESULTS

Reasonable agreement was observed between the treatment planning system calculations of relative dose distribution and the measured data. The maximum observed deviation in the spatial position between the center of the measured and calculated dose profiles was 0.6mm. The maximum observed difference in full width of half maximum between calculated and measured profiles was 1.2mm.

CONCLUSIONS

The use of polymer-gel dosimetry for a verification of stereotactic procedures has some unique advantages that can be summarized as follows: the dosimeter itself is tissue equivalent, three-dimensional dose distributions can be measured and the dosimeter allows simulation of the patient's procedures without any limitations.

Bilateral internal globus pallidus stimulation for the treatment of spasmodic torticollis

Three patients with spasmodic torticollis (ST) obtained substantial benefit from bilateral globus pallidus internus (GPi) stimulation. Progressive improvement in ST occurred over several months but residual cervical dystonia remained. These results corroborate those obtained by Krauss et al. on three patients with ST.

A phantom study of the geometric accuracy of computed tomographic and magnetic resonance imaging stereotactic localization with the Leksell stereotactic system

OBJECTIVE

To assess the spatial accuracy of magnetic resonance imaging (MRI) and computed tomographic stereotactic localization with the Leksell stereotactic system.

METHODS

The phantom was constructed in the shape of a box, 164 mm in each dimension, with three perpendicular arrays of solid acrylic rod, 5 mm in diameter and spaced 30 mm apart within the phantom. In this study, images from two different MRI scanners and a computed tomographic scanner were obtained using the same Leksell (Elekta Instruments, Stockholm, Sweden) head frame placement. The coordinates of the rod images in the three principal planes were measured by using a tool provided with Leksell GammaPlan software (Elekta Instruments, Norcross, GA) and were compared with the physical phantom measurements.

RESULTS

The greatest distortion was found around the periphery, and the least distortion (<1.5 mm) was present in the middle and most other areas of the phantom. In the phantom study using computed tomography, the mean values of the maximum errors for the x, y, and z axes were 1.0 mm (range, 0.2-1.3 mm), 0.4 mm (range, 0.1-0.8 mm), and 3.8 mm (range, 1.9-5.1 mm), respectively. The mean values of the maximum errors when using the Philips MRI scanner (Philips Medical Systems, Shelton, CT) were 0.9 mm (range, 0.4-1.7 mm), 0.2 mm (range, 0.0-0.7 mm), and 1.9 mm (range, 1.3-2.3 mm), respectively. Using the Siemens MRI scanner (Siemens Medical Systems, New York, NY), these values were 0.4 mm (range, 0.0-0.7 mm), 0.6 mm (range, 0.0-1.0 mm), and 1.6 mm (range, 0.8-2.0 mm), respectively. The geometric accuracy of the MRI scans when using the Siemens scanner was greatly improved after the implementation of a new software patch provided by the manufacturer. The accuracy also varied with the direction of phase encoding.

CONCLUSION

The accuracy of target localization for most intracranial lesions during stereotactic radiosurgery can be achieved within the size of a voxel, especially by using the Siemens MRI scanner at current specifications and with a new software patch. However, caution is warranted when imaging peripheral lesions, where the distortion is greatest.

An image fusion study of the geometric accuracy of magnetic resonance imaging with the Leksell stereotactic localization system

A special acrylic phantom designed for both magnetic resonance imaging (MRI) and computed tomography (CT) was used to assess the geometric accuracy of MRI-based stereotactic localization with the Leksell stereotactic head frame and localizer system. The acrylic phantom was constructed in the shape of a cube, 164 mm in each dimension, with three perpendicular arrays of solid acrylic rods, 5 mm in diameter and spaced 30 mm apart within the phantom. Images from two MR scanners and a CT scanner were obtained with the same Leksell head frame placement. Using image fusion provided by the Leksell GammaPlan (LGP) software, the coordinates of the intraphantom rod positions from two MRI scanners were compared to that of CT imaging. The geometric accuracy of MR images from the Siemens scanner was greatly improved after the implementation of a special software patch provided by the manufacturer. In general, much better accuracy was achieved in the transverse plane where images were acquired. Most distortion was found around the periphery while least distortion was present in the middle and most other parts of the phantom. For most intracranial lesions undergoing stereotactic radiosurgery, accuracy of target localization can be achieved within size of a voxel, especially with the Siemens scanner. However, extra caution should be taken for imaging of peripheral lesions where the distortion is the greatest.