Dosimetric accuracy of the Convolution algorithm for Leksell Gamma Plan radiosurgery treatment planning: Evaluation in the presence of clinically relevant inhomogeneities.

PURPOSE: The Leksell Gamma Plan Convolution algorithm (LGP-Convolution) has not been widely adopted. This mainly stems from the higher calculated beam-on times relative to the standard ray tracing-based LGP-TMR10 dose calculation algorithm. This study aims to evaluate the accuracy of the LGP-Convolution in scenarios where the treated lesions are in the vicinity of or encompassed by bone and/or air inhomogeneities. METHODS: The solid water dosimetry phantom provided by the vendor was modified to include bone and air inhomogeneities. Two treatment planning scenarios were investigated involving a single shot and multiple shots, respectively. Treatment planning and dose prescription were performed using the LGP-Convolution algorithm. Triple channel film dosimetry was performed using GafChromic EBT3 films calibrated in terms of absorbed dose to water in a 60 Co beam. Monte Carlo (MC) simulation dosimetry was also performed in the inhomogeneous experimental geometry using the EGSnrc MC platform and a previously validated sector-based phase-space source model. MC simulations were also employed to determine correction factors required for converting EBT3 measurements at points within the bone and air inhomogeneities from dose-to-water values to the corresponding dose to medium values. RESULTS AND CONCLUSIONS: EBT3 dose to medium correction factors ranged with field size (4, 8, or 16 mm) within 0.941-0.946 for bone and 0.745-0.749 for air inhomogeneities. An excellent agreement was found between the LGP-Convolution calculations with corresponding EBT3 and MC dose to medium results at all measurement points, except those located inside the air inhomogeneity. The latter is of no clinical importance and excluding them yielded gamma index passing rates of nearly 100% for 3% local dose difference and 1 mm distance-to-agreement criteria. The excellent agreement observed between LGP-Convolution calculations and film as well as MC results of dose to medium indicates that the latter is the quantity reported by the LGP-Convolution.

Cone-Beam CT image contrast and attenuation-map linearity improvement (CALI) for brain stereotactic radiosurgery procedures.

A Contrast and Attenuation-map Linearity Improvement (CALI) framework is proposed for cone-beam CT (CBCT) images used for brain stereotactic radiosurgery (SRS). The proposed framework is tailored to improve soft tissue contrast of a new point-of-care image-guided SRS system that employs a challenging half cone beam geometry, but can be readily reproduced on any CBCT platform. CALI includes a pre- and post-processing step. In pre-processing we apply a shading and beam hardening artifact correction to the projections, and in post-processing step we correct the dome/capping artifact on reconstructed images caused by the spatial variations in X-ray energy generated by the bowtie-filter. The shading reduction together with the beam hardening and dome artifact correction algorithms aim to improve the linearity and accuracy of the CT-numbers (CT#). The CALI framework was evaluated using CatPhan to quantify linearity, contrast-to-noise (CNR), and CT# accuracy, as well as subjectively on patient images acquired on a clinical system. Linearity of the reconstructed attenuation-map was improved from 0.80 to 0.95. The CT# mean absolute measurement error was reduced from 76.1 to 26.9 HU. The CNR of the acrylic insert in the sensitometry module was improved from 1.8 to 7.8. The resulting clinical brain images showed substantial improvements in soft tissue contrast visibility, revealing structures such as ventricles which were otherwise undetectable in the original clinical images obtained from the system. The proposed reconstruction framework also improved CT# accuracy compared to the original images acquired on the system. For frameless image-guided SRS, improving soft tissue visibility can facilitate evaluation of MR to CBCT co-registration. Moreover, more accurate CT# may enable the use of CBCT for daily dose delivery measurements.

Variability in target delineation for cavernous sinus meningioma and anaplastic astrocytoma in stereotactic radiosurgery with Leksell Gamma Knife Perfexion.

BACKGROUND: Radiosurgery clinical practice relays on empirical observations and the experience of the practitioners involved in determining and delineating the target and therefore variability in target delineation might be expected for all the radiosurgery approaches, independent of the technique and the equipment used for delivering the treatment. The main aim of this study was to quantify the variability of target delineation for two radiosurgery targets expected to be difficult to delineate. The secondary aim was to investigate the dosimetric implications with respect to the plan conformity. The primary aim of the study has therefore a very general character, not being bound to one specific radiosurgery technique. MATERIALS AND METHODS: Twenty radiosurgery centers were asked to delineate one cavernous sinus meningioma and one astrocytoma and to plan the treatments for Leksell Gamma Knife Perfexion. The analysis of the delineated targets was based on the calculated 50 % agreement volume, AV50. The AV50 was compared to each delineated target by the concordance index and discordance index. The differences in location, size, and shape of the delineated targets were also analyzed using the encompassing volume compared to the common volume, i.e., the AV100, of all delineated structures. RESULTS: Target delineation led to major differences between the participating centers and therefore the AV50 was small in comparison to each delineated target volume. For meningioma, the AV50 was 5.90 cm(3), the AV100 was 2.60 cm(3), and the encompassing volume was 13.14 cm(3). For astrocytoma, the AV50 was 2.06 cm(3) while the AV100 was extremely small, only 0.05 cm(3), and the encompassing volume was 43.27 cm(3). These variations translate into corresponding discrepancies in plan conformity. CONCLUSIONS: Significant differences in shape, size, and location between the targets included in this study were identified and therefore the clinical implications of these differences should be further investigated.

Physical dose validation of dynamic treatment for Gamma Knife radiosurgery.

BACKGROUND: Dynamic treatment in Gamma Knife (GK) radiosurgery systems delivers radiation continuously with couch movement, as opposed to stationary step-and-shoot treatment where radiation is paused when moving between isocenters. Previous studies have shown the potential for dynamic GK treatment to give faster treatment times and improved dose conformity and homogeneity. However, these studies focused only on computational simulations and lack physical validation. PURPOSE: This study aims conduct dynamic treatment dosimetric validation with physical experimental measurements. The experiments aim to (1) address assumptions made with computational studies, such as the validity of treating a continuous path as discretised points, (2) investigate uncertainties in translating computed plans to actual treatment, and (3) determine ideal treatment planning parameters, such as interval distance for the path discretization, collimator change limitations, and minimum isocenter treatment times. METHODS: This study uses a GK ICON treatment delivery machine, and a motion phantom custom-made to attach to the machine's mask adapter and move in 1D superior-inferior motion. Phantom positioning is first verified through comparisons against couch motion and computed doses. For dynamic treatment experiments, the phantom is moved through a program that first reads the desired treatment plan isocenters' position, time, and collimator sizes, then carries out the motion continuously while the treatment machine delivers radiation. Measurements are done with increasing levels of complexity: varying speed, varying collimator sizes, varying both speed and collimator sizes, then extends the same measurements to simulated 2D motion by combining phantom and couch motion. Dose comparisons between phantom motion radiation measurements and either couch motion measurements or dose calculations are analyzed with 2 mm/2% and 1 mm/2% gamma indices, using both local and global gamma index calculations. RESULTS: Phantom positional experiments show a high accuracy, with global gamma indices for all dose comparisons ≥ $\ge $ 99%. Discretization level to approximate continuous path as discrete points show the good dose matches with dose calculations when using 1 and 2-mm gaps. Complex 1D motion, including varying speed, collimator sizes, or both, as well as 2D motion with the same complexities, all show good dose matches with dose calculations: the scores are ≥ $\ge $ 92.0% for the strictest 1 mm/2% local gamma index calculation, ≥ $\ge $ 99.8% for 2 mm/2% local gamma index, and ≥ $\ge $ 97.0% for all global gamma indices. Five simulated 2D treatments with optimized plans scored highly as well, with all gamma index scores ≥ $\ge $ 95.3% when compared to stationary treatment, and scores ≥ $\ge $ 97.9% when compared to plan calculated dose. CONCLUSIONS: Dynamic treatment computational studies are validated, with dynamic treatment shown to be physically feasible and deliverable with high accuracy. A 2-mm discretization level in treatment planning is proposed as the best option for shorter dose calculation times while maintaining dose accuracy. Our experimental method enables dynamic treatment measurements using the existing clinical workflow, which may be replicated in other centers, and future studies may include 2D or 3D motion experiments, or planning studies to further quantify potential indication-specific benefits.

A dose calculation algorithm for radiosurgery treatment room shielding optimization.

BACKGROUND: Methods described in the literature for designing shielding for treatment rooms for radiotherapy systems often involve assumptions that lead to overestimations of the wall thicknesses required to meet dose rate constraints outside the room. The Leksell Gamma Knife (LGK) has built-in shielding that results in primarily scattered photons leaking into the room. The field of leakage radiation, therefore, has a wide spectrum of energies, up to the primary energies of cobalt-60, and is highly anisotropic, making standard site planning methods difficult to adapt to the LGK. Monte Carlo (MC) simulations are an alternative for dose calculations but are computationally expensive. PURPOSE: The aim was to develop a dose calculation algorithm for fast estimations of dose rates outside the barriers of a treatment room. The algorithm could then be employed in iterative methods to optimize treatment room parameters such as wall thicknesses and position of the LKG unit. METHODS: The algorithm uses pre-calculated MC simulation data in two steps. First, it uses phase spaces that describe single photons in the radiation field around the LGK. Using each photon's position and direction, they are raytraced to the outside of the room. Based on their energy, angle of incidence to the barrier, and the barrier's thickness, a depth-dose profile is sampled from a pre-calculated library of profiles. The contribution from each photon is added to the space outside the barriers, giving the total dose distribution. The barrier thicknesses are optimized by iteratively running the algorithm and finding the minimum thickness required to keep the resulting maximum dose rate outside below a given upper limit. RESULTS: Comparisons between dose distributions from the dose algorithm and full MC simulations of two typical rooms show good agreement, with the maximum dose rate outside each barrier being estimated within what represents a 2 cm error in concrete barrier thickness. The algorithm is successfully used in optimizing the barrier thicknesses and the results show potential decreases in most barriers' thickness, compared to a conventional method. CONCLUSIONS: The algorithm is a promising alternative to conventional barrier design methods, eliminating several shortcomings of the latter whilst being significantly faster than a full MC simulation. The flexibility of the algorithm would further enable solving more complex cases, such as optimizing the LGK position or barrier material to further reduce material use and cost.

Plan Assessment Metrics for Dose Painting in Stereotactic Radiosurgery.

Purpose: As radiation therapy treatment precision increases with advancements in imaging and radiation delivery, dose painting treatment becomes increasingly feasible, where targets receive a nonuniform radiation dose. The high precision of stereotactic radiosurgery (SRS) makes it a good candidate for dose painting treatments, but no suitable metrics to assess dose painting SRS plans exist. Existing dose painting assessment metrics weigh target overdose and underdose equally but are unsuited for SRS plans, which typically avoid target underdose more. Current SRS metrics also prioritize reducing healthy tissue dose through selectivity and dose fall-off, and these metrics assume single prescriptions. We propose a set of metrics for dose painting SRS that would meet clinical needs and are calculated with nonuniform dose painting prescriptions. Methods and Materials: Sample dose painting SRS prescriptions are first created from Gamma Knife SRS cases, apparent diffusion coefficient magnetic resonance images, and various image-to-prescription functions. Treatment plans are found through semi-infinite linear programming optimization and using clinically determined isocenters, then assessed with existing and proposed metrics. Modified versions of SRS metrics are proposed, including coverage, selectivity, conformity, efficiency, and gradient indices. Quality factor, a current dose painting metric, is applied both without changes and with modifications. A new metric, integral dose ratio, is proposed as a measure of target overdose. Results: The merits of existing and modified metrics are demonstrated and discussed. A modified conformity index using mean or minimum prescription dose would be suitable for dose painting SRS with integral or maximum boost methods, respectively. Either modified efficiency index is a suitable replacement for the existing gradient index. Conclusions: The proposed modified SRS metrics are appropriate measures of plan quality for dose painting SRS plans and have the advantage of giving equal values as the original SRS metrics when applied to single-prescription plans.

Characterization of the radiation field surrounding the Leksell Gamma Knife® and shielding applications.

The aim of this study is to improve the characterization and modeling of the radiation field surrounding the Leksell Gamma Knife®-PerfexionTM. The improved characterization of the radiation field enables more accurate shielding calculations to be performed for the areas adjacent to the treatment room. With the aid of a high-purity germanium detector and a satellite dose rate meter, γ-ray spectra and ambient dose equivalent H*(10) data were acquired at various locations in the field of a Leksell Gamma Knife unit in a treatment room at Karolinska University Hospital, Sweden. These measurements were used to validate the results of the PEGASOS Monte Carlo simulation system with a PENELOPE kernel. The levels of the radiation that passes through the shielding of the machine (leakage radiation) are shown to be much lower than what is suggested by various bodies, e.g. the National Council on Radiation Protection and Measurements, to be used when calculating radiation shielding barriers. The results clearly indicate that Monte Carlo simulations may be used in structural shielding design calculations for γ rays from the Leksell Gamma Knife.

Radiobiological Evaluation of Combined Gamma Knife Radiosurgery and Hyperthermia for Pediatric Neuro-Oncology.

Combining radiotherapy (RT) with hyperthermia (HT) has been proven effective in the treatment of a wide range of tumours, but the combination of externally delivered, focused heat and stereotactic radiosurgery has never been investigated. We explore the potential of such treatment enhancement via radiobiological modelling, specifically via the linear-quadratic (LQ) model adapted to thermoradiotherapy through modulating the radiosensitivity of temperature-dependent parameters. We extend this well-established model by incorporating oxygenation effects. To illustrate the methodology, we present a clinically relevant application in pediatric oncology, which is novel in two ways. First, it deals with medulloblastoma, the most common malignant brain tumour in children, a type of brain tumour not previously reported in the literature of thermoradiotherapy studies. Second, it makes use of the Gamma Knife for the radiotherapy part, thereby being the first of its kind in this context. Quantitative metrics like the biologically effective dose (BED) and the tumour control probability (TCP) are used to assess the efficacy of the combined plan.

8+ Year Performance of the Gamma Knife Perfexion/Icon Patient Positioning System and Possibilities for Preemptive Fault Detection Using Statistical Process Control.

BACKGROUND: The large fractional doses, steep dose gradients, and small targets found in intracranial radiosurgery require extremely low beam delivery uncertainty. In the case of Gamma Knife radiosurgery (GKRS), this includes minimizing patient positioning system (PPS) positioning uncertainty. Existing QA techniques are recipe based, and feature point in time pass/fail tolerances. However, modern treatment machines, including the Gamma Knife Perfexion/Icon systems, record extensive internal data in treatment logs. These data can be analyzed through statistical process control (SPC) methods which are designed to detect changes in process behavior. The purpose of this study was to characterize the long-term (8+ year) performance of a Perfexion/Icon unit and use SPC methods to determine if performance changes could be detected at levels lower than existing QA and internal manufacturer performance tolerances. METHODS: In-house software was developed to parse Perfexion/Icon log-files and store relevant information on shot delivery in a relational database. A last-in, first-out (LIFO) queuing algorithm was created to heuristically match messages associated with a given delivered shot. Filtering criteria were developed to filter QA and uncompleted shots. The resulting matched shots were extracted. Achieved versus planned PPS position was determined for each PPS motor as well as for the vector magnitude difference in PPS position. Exponentially weighted moving average (EWMA) control charts were plotted to determine when process behavior changed over time. RESULTS: 53833 shots were delivered over an 8+ year span in the study. The mean vector magnitude PPS difference was 32.7 µm, with 97.5% of all shots within 70.1 µm. Several changes in PPS positioning behavior were observed over time, corresponding with control system faults on several occasions requiring PPS recalibration. EWMA control charts clearly demonstrate that these faults could be identified and possibly predicted as many as 3 years before there were faults beyond control system tolerance. CONCLUSION: The PPS of Gamma Knife Perfexion/Icon systems has extremely low positioning uncertainties. EWMA control chart method can be utilized to track PPS performance over time and can potentially detect changes in performance that may indicate a component requiring maintenance. This would allow planned service visits to mitigate problems and prevent unplanned downtime.

Using Monte-Carlo-simulated radiation transport to calculate dose distribution in rats before irradiation with Leksell Gamma Knife 4C: technical note.

BACKGROUND: Gamma knife surgery (GKS) is used at subnecrotic doses for temporal lobe epilepsy (TLE) treatment. Rat models of TLE have been used to probe the mechanisms underlying GKS. Previous GKS studies on rats have used the Leksell GammaPlan (LGP) treatment planning system to determine the irradiation time to achieve the dose to deliver. Since LGP is not designed for such small structures, it is important to calibrate the system for the rat brain. METHODS: We have used a Monte Carlo simulation (MCS) radiation transport scheme, with CT data as anatomical and tissue-specific information, to simulate the dose distribution in a rat brain when using a Leksell Gamma Knife. RESULTS: We show how dose distributions obtained by MCS quantitatively compare to those predicted by LGP, and discuss whether LGP should be used for studies involving rats. The energy deposited when using the 4-mm collimators was calculated for targets on both sides of the rat brain in the dorsal hippocampus, which allowed us to determine the exact time to irradiate rats with a given dose. CONCLUSION: The MCS method used in this study can easily be used for future GKS studies on small animals when accurate dose distributions are required.

Technical Principles of Dual-Energy Cone Beam Computed Tomography and Clinical Applications for Radiation Therapy.

PURPOSE: Medical imaging is an indispensable tool in radiotherapy for dose planning, image guidance and treatment monitoring. Cone beam CT (CBCT) is a low dose imaging technique with high spatial resolution capability as a direct by-product of using flat-panel detectors. However, certain issues such as x-ray scatter, beam hardening and other artifacts limit its utility to the verification of patient positioning using image-guided radiotherapy. METHODS AND MATERIALS: Dual-energy (DE)-CBCT has recently demonstrated promise as an improved tool for tumor visualization in benchtop applications. It has the potential to improve soft-tissue contrast and reduce artifacts caused by beam hardening and metal. In this review, the practical aspects of developing a DE-CBCT based clinical and technical workflow are presented based on existing DE-CBCT literature and concepts adapted from the well-established library of work in DE-CT. Furthermore, the potential applications of DE-CBCT on its future role in radiotherapy are discussed. RESULTS AND CONCLUSIONS: Based on current literature and an investigation of future applications, there is a clear potential for DE-CBCT technologies to be incorporated into radiotherapy. The applications of DE-CBCT include (but are not limited to): adaptive radiotherapy, brachytherapy, proton therapy, radiomics and theranostics.

Investigation of intracranial peripheral dose arising from the treatment of large lesions with Leksell GammaKnife Perfexion.

This investigation involves quantifying the extent of intracranial peripheral dose arising from simulated targets situated in the skull-base or upper-spine region using the Leksell GammaKnife Perfexion treatment unit. For each of three spherical target volumes--denoted as Vs (4 cm3), VM (18 cm3), and VL (60 cm3)--three treatment plans were manually generated, one for each of the three collimator sizes--4, 8, and 16 mm. Each of the plans was delivered to a spherical dosimetry phantom with an insert containing EBT Gafchromic film. The total dose at 70 mm from the targets' edges, %D(70 mm), was measured as a function of elevation angle and expressed as a percentage of the prescription dose. The film insert was placed centered in the median sagittal plane (Leksell X = 100) and %D(70 mm) was measured for the angular range from 0 degree (superior/along Z axis) to 90 degrees (anterior/along Y axis). For a given collimator i, the irradiation time ti to treat a spherical target of volume V using the 50% isodose line was observed to follow a power-law relationship of the form ti = Ai(V/ Vi)n where Ai was the maximum dose divided by collimator dose rate and Vi was the volume encompassed by the 50% isodose line for a single shot. The mean value of n was 0.61 (range: 0.61-0.62). Along the superior (Z) direction (angle=0 degree) and up to angles of around 30 degrees, the %D(70 mm) was always highest for the 4 mm plans, followed by the 8 mm, followed by the 16 mm. In this angular range, the maximum measured %D(70 mm) was 1.7% of the prescription dose. The intracranial peripheral dose along the superior direction (combined scatter and leakage dose) resulting from irradiation of upper-spine or base-of-skull lesions is measured to be less than 2% of the prescription dose, even for very large (60 cm3) targets. The results of this study indicate that, for a given target volume, treatment plans consisting of only 4 mm shots yield larger peripheral dose in the superior direction than 8 mm shot only plans, which in turn yield larger peripheral dose than 16 mm shot only plans.

A simple and effective method for validation and measurement of collimator output factors for Leksell Gamma Knife Perfexion.

Accurate determination of collimator output factors is important for Leksell Gamma Knife radiosurgery. The new Leksell Gamma Knife Perfexion system has a completely redesigned collimator system and the collimator output factors are different from the other Leksell Gamma Knife models. In this study, a simple method was developed to validate the collimator output factors specifically for Leksell Gamma Knife Perfexion. The method uses double-shot exposures on a single film to eliminate repeated setups and the necessity to acquire dose calibration curves required for the traditional film exposure method. Using the method, the collimator output factors with respect to the 16 mm collimator were measured to be 0.929 +/- 0.009 and 0.817 +/- 0.012 for the 8 mm and the 4 mm collimator, respectively. These values are in agreement (within 2%) with the default values of 0.924 and 0.805 in the Leksell Gamma Plan treatment planning system. These values also agree with recently published results of 0.917 (8 mm collimator) and 0.818 (4 mm collimator) obtained from the traditional methods. Given the efficiency of the method, measurement and validation of the collimator output factors can be readily adopted in commissioning and quality assurance of a Leksell Gamma Knife Perfexion system.

Gamma Knife radiosurgery: Scenarios and support for re-irradiation.

Stereotactic radiosurgery (SRS) involves the focal delivery of large, cytotoxic doses of radiation to small targets within the brain, often located in close proximity to radiosensitive normal tissue structures and requiring very low procedural uncertainties to perform safely. Historically, neurosurgeons considered SRS as a one-time, single session procedure. However therapeutic advances and a better understanding of the clinical response to SRS have caused a renewal of interest in a variety of re-irradiation scenarios; including re-irradiation of the same target after prior SRS, SRS treatments after prior broad-field radiation, hypofractionated treatments, and volume-staged treatments. Re-irradiation may in some cases require even greater effort towards minimizing treatment uncertainties as compared to one-time-only treatments. Gamma Knife radiosurgery (GKRS) has evolved over time in ways that directly supports many re-irradiation scenarios while helping to minimize overall procedural uncertainty.

Clinical Image Coregistration Variability on a Dedicated Radiosurgery Unit.

BACKGROUND: On a new dedicated radiosurgery unit enabling frameless treatments, a cone-beam computed tomography (CBCT) can be used for stereotactic definition. Since magnetic resonance imaging (MRI) is used to delineate target, reproducible MRI-to-CBCT coregistration is vital for accurate target localization. OBJECTIVE: To evaluate reproducibility of image coregistration in patient images. METHODS: Three types of coregistration (source-to-target) were analyzed: (1) MRI-to-CT; (2) MRI-to-CBCT; and (3) CT-to-CBCT. For each patient (n = 15), each coregistration type was independently performed 5 to 30 times (total: 465 coregistrations). Each coregistration yielded a transformation matrix, which was subsequently applied to transform every point in the source image to stereotactic coordinates. Two metrics were measured: (1) target registration error (TRE): mean distance between the registered position of each target point and the average registration position of that point; (2) compound registration error (CRE): mean spatial difference between stereotactic coordinates using (A) MRI-to-CT-to-CBCT and (B) MRI-to-CBCT. RESULTS: The median (range) of TRE was 0.11 mm (0.06-0.22 mm), 0.17 mm (0.10-0.36 mm), and 0.12 mm (0.08-0.21 mm) for MRI-to-CT, MRI-to-CBCT, and CT-to-CBCT, respectively. The TRE for MRI-to-CBCT was statistically higher than the other 2 methods (P < .01). The median (range) of CRE was 0.44 mm (0.22-0.59 mm). The maximum point CRE between patients ranged from 0.37-1.15 mm when considering all MRI points, but reduced to 0.31-0.90 mm within the central 16 cm. The CRE varied across the image volume, and typically was minimized near the center. CONCLUSION: The variation in image coregistration is within 0.2 mm, indicating a high degree of reproducibility. The CRE varies throughout the head but is submillimeter in the central 16 cm region.

Simultaneous deblurring and iterative reconstruction of CBCT for image guided brain radiosurgery.

One of the limiting factors in cone-beam CT (CBCT) image quality is system blur, caused by detector response, x-ray source focal spot size, azimuthal blurring, and reconstruction algorithm. In this work, we develop a novel iterative reconstruction algorithm that improves spatial resolution by explicitly accounting for image unsharpness caused by different factors in the reconstruction formulation. While the model-based iterative reconstruction techniques use prior information about the detector response and x-ray source, our proposed technique uses a simple measurable blurring model. In our reconstruction algorithm, denoted as simultaneous deblurring and iterative reconstruction (SDIR), the blur kernel can be estimated using the modulation transfer function (MTF) slice of the CatPhan phantom or any other MTF phantom, such as wire phantoms. The proposed image reconstruction formulation includes two regularization terms: (1) total variation (TV) and (2) nonlocal regularization, solved with a split Bregman augmented Lagrangian iterative method. The SDIR formulation preserves edges, eases the parameter adjustments to achieve both high spatial resolution and low noise variances, and reduces the staircase effect caused by regular TV-penalized iterative algorithms. The proposed algorithm is optimized for a point-of-care head CBCT unit for image-guided radiosurgery and is tested with CatPhan phantom, an anthropomorphic head phantom, and 6 clinical brain stereotactic radiosurgery cases. Our experiments indicate that SDIR outperforms the conventional filtered back projection and TV penalized simultaneous algebraic reconstruction technique methods (represented by adaptive steepest-descent POCS algorithm, ASD-POCS) in terms of MTF and line pair resolution, and retains the favorable properties of the standard TV-based iterative reconstruction algorithms in improving the contrast and reducing the reconstruction artifacts. It improves the visibility of the high contrast details in bony areas and the brain soft-tissue. For example, the results show the ventricles and some brain folds become visible in SDIR reconstructed images and the contrast of the visible lesions is effectively improved. The line-pair resolution was improved from 12 line-pair/cm in FBP to 14 line-pair/cm in SDIR. Adjusting the parameters of the ASD-POCS to achieve 14 line-pair/cm caused the noise variance to be higher than the SDIR. Using these parameters for ASD-POCS, the MTF of FBP and ASD-POCS were very close and equal to 0.7 mm-1 which was increased to 1.2 mm-1 by SDIR, at half maximum.

Dosimetric effects of Onyx embolization on Gamma Knife arteriovenous malformation dose distributions.

OBJECTIVE Patients with arteriovenous malformations (AVMs) treated with Gamma Knife radiosurgery (GKRS) subsequent to embolization suffer from elevated local failure rates and differences in adverse radiation effects. Onyx is a common embolic material for AVMs. Onyx is formulated with tantalum, a high atomic number (Z = 73) element that has been investigated as a source of dosimetric uncertainty contributing to the less favorable clinical results. However, prior studies have not modeled the complicated anatomical and beam geometries characteristic of GKRS. This study investigated the magnitude of dose perturbation that can occur due to Onyx embolization using clinically realistic anatomical and Gamma Knife beam models. METHODS Leksell GammaPlan (LGP) was used to segment the AVM nidus and areas of Onyx from postcontrast stereotactic MRI for 7 patients treated with GKRS postembolization. The resulting contours, skull surface, and clinically selected dose distributions were exported from LGP in DICOM-RT (Digital Imaging and Communications in Medicine-radiotherapy) format. Isocenter locations and dwell times were recorded from the LGP database. Contours were converted into 3D mesh representations using commercial and in-house mesh-editing software. The resulting data were imported into a Monte Carlo (MC) dose calculation engine (Pegasos, Elekta Instruments AB) with a beam geometry for the Gamma Knife Perfexion. The MC-predicted dose distributions were calculated with Onyx assigned manufacturer-reported physical constants (MC-Onyx), and then compared with corresponding distributions in which Onyx was reassigned constants for water (MC-water). Differences in dose metrics were determined, including minimum, maximum, and mean dose to the AVM nidus; selectivity index; and target coverage. Combined differences in dose magnitude and distance to agreement were calculated as 3D Gamma analysis passing rates using tolerance criteria of 0.5%/0.5 mm, 1.0%/1.0 mm, and 3.0%/3.0 mm. RESULTS Overall, the mean percentage differences in dose metrics for MC-Onyx relative to MC-water were as follows; all data are reported as mean (SD): minimum dose to AVM = -0.7% (1.4%), mean dose to AVM = 0.1% (0.2%), maximum dose to AVM = 2.9% (5.0%), selectivity = 0.1% (0.2%), and coverage = -0.0% (0.2%). The mean percentage of voxels passing at each Gamma tolerance were as follows: 99.7% (0.1%) for 3.0%/3.0 mm, 98.2% (0.7%) for 1.0%/1.0 mm, and 52.1% (4.4%) for 0.5%/0.5 mm. CONCLUSIONS Onyx embolization appears to have a detectable effect on the delivered dose distribution. However, the small changes in dose metrics and high Gamma passing rates at 1.0%/1.0 mm tolerance suggest that these changes are unlikely to be clinically significant. Additional sources of delivery and biological uncertainty should be investigated to determine the root cause of the observed less favorable postembolization GKRS outcomes.

Application of the concept of biologically effective dose (BED) to patients with Vestibular Schwannomas treated by radiosurgery.

In the application of stereotactic radiosurgery, using the Gamma Knife, there are large variations in the overall treatment time for the same prescription dose, given in a single treatment session, for different patients. This is due to not only changes in the activity of the Cobolt-60 sources, but also to variations in the number of iso-centers used, the collimator size for a particular iso-center, and the time gap between the different iso-centers. Although frequently viewed as a single dose treatment the concept of biologically effective dose (BED), incorporating concurrent fast and a slow components of repair of sublethal damage, would imply potential variations in BED because of the influence of these different variables associated with treatment. This was investigated in 26 patients, treated for Vestibular Schwannomas, using the Series B Gamma-Knife, between 1999 and 2005. The iso-center number varied between 2 and 13, and the overall treatment time from 25.4-129.58 min. The prescription doses varied from 10-14 Gy. To obtain physical dose and dose-rates from each iso-center, in a number of locations in the region of interest, a prototype version of the Leksell GammaPlan® was used. For an individual patient, BED values varied by up to 15% for a given physical iso-dose. This was due to variation in the dose prescription at different locations on that iso-dose. Between patients there was a decline in the range of BED values as the overall treatment time increased. This increased treatment time was partly a function of the slow decline in the activity of the sources with time but predominantly due to changes in the number of iso-centers used. Thus, variations in BED values did not correlate with prescription dose but was modified by the overall treatment time.