Next Article in Journal
Radiomics-Guided Precision Radiation Therapy in Head and Neck Squamous Cell Carcinoma
Previous Article in Journal
Improving the Accuracy of Bone-Scintigraphy Imaging Analysis Using the Skeletal Count Index: A Study Based on Human Trial Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A 3D Superposition Approximation for Gamma Knife Dose Calculation

1
Department of Physics and Astronomy, Louisiana State University, Baton Rouge, LA 70803, USA
2
Physics Department, Mary Bird Perkins Cancer Center, Baton Rouge, LA 70809, USA
3
Physics and Biology in Medicine, David Geffen School of Medicine, University of California, Los Angeles, CA 90095, USA
*
Author to whom correspondence should be addressed.
Radiation 2025, 5(1), 6; https://doi.org/10.3390/radiation5010006
Submission received: 21 November 2024 / Revised: 11 January 2025 / Accepted: 15 January 2025 / Published: 20 January 2025

Simple Summary

Gamma Knife stereotactic radiosurgery is a precise radiation therapy technique that delivers a high dose of radiation to treat patients with brain abnormalities. The treatment planning and dose optimization process can be complicated and time-consuming, requiring careful calculations to ensure that radiation is precisely and accurately delivered to the target while minimizing complications in healthy brain tissues. This study proposes a straightforward method to approximate the distribution of radiation dose in three dimensions outside of the standard planning software, GammaPlan. By deploying a simple sector-based superposition method, we aim to create a tool for evaluating treatment plans and improve flexibility in treatment design. This new approach could help clinicians and researchers assess plan quality with less reliance on proprietary software, potentially enhancing the overall treatment planning process.

Abstract

Effective dose calculation is essential for optimizing Gamma Knife (GK) stereotactic radiosurgery (SRS) treatment plans. Modern GK systems allow independent sector activation, enabling complex dose distributions per shot. This study presents a dose approximation method designed to account for shot flexibility and generate 3D doses external to GammaPlan. A treatment plan was created with the TMR10 calculation for individual sector activations using a Radiosurgery Head Phantom. The resulting dose arrays established a basis set of sector-specific distributions, which were then referenced by shot parameters from the plan, allowing dose accumulation through superposition. This superposition approximation (SA) was compared to the original TMR10 using the Dice Similarity Coefficient (DSC), 95% Hausdorff Distance (HD95), and GK deliverability metrics: coverage, selectivity, and gradient index, across an isodose normalization range from 10% to 90%. In a cohort of 30 patients with 71 targets, strong agreement was observed between TMR10 and SA in the clinically used 50–60% isodose range, with DSC above 85% and HD95 under 2.18 mm. The average differences for the coverage, selectivity, and gradient index were 0.014, 0.008, and 0.118, respectively. This method accurately approximates TMR10 calculations within clinically relevant ranges, offering an external tool to assess 3D dose distributions for GK treatment plans.

1. Introduction

Optimizing treatment plans in radiation therapy is a challenging task across many modalities. Depending on the optimizer, this process often requires multiple dose calculations at various stages. Iterative optimizers, commonly employed in clinical practice, perform simplified dose calculations between iterations to evaluate the plan’s progress toward achieving treatment objectives [1]. For this reason, dose calculation approximations, especially those capable of calculating 3D dose distributions, are necessary to continue the development of plan optimization systems. In this study, a simplified dose calculation method was developed to approximate 3D dose distributions from plans generated for the Leksell Gamma Knife Icon (Elekta AB, Stockholm, Sweden).
Gamma Knife stereotactic radiosurgery (GK-SRS) is commonly used to treat small targets, such as brain metastases, with doses ranging from 18–24 Gy prescribed to the 50–60% isodose line [2,3,4]. Larger targets, such as meningiomas, can be treated with lower doses of 12 Gy or more for patients that are not surgical candidates [5]. Very high doses of 80 Gy or more may be used to provide relief to patients from function disorders like trigeminal neuralgia [6].
GK-SRS systems are designed with submillimeter precision. They utilize 192 cobalt-60 (Co-60) sources, each individually collimated in a non-coplanar configuration around the patient’s head. These sources are grouped into eight sectors, with each sector capable of collimating its sources independently of the other sectors. Each sector can be blocked or opened to one collimation of 4, 8, or 16 mm in size [7]. At the isocenter, the convergence of the 192 beamlets results in a highly focused dose, delivered with an accuracy of 0.15 mm. During treatment, the combination of machine parameters for patient position relative to the isocenter, sector activation of the sources, and beam-on time at that position form what is known as a “shot”. A GK treatment plan contains the collection of this information for all the shots required to treat that patient. The optimal arrangement of these plan parameters, namely the number of shots, shot position, sector activation, and shot weighting allows GK to achieve the core concept of SRS, precisely delivering high doses of radiation to each target lesion, often in a single fraction (SRS), or fractionated stereotactic radiotherapy (SRT), while simultaneously sparing surrounding healthy tissues [8,9].
Previous versions of the GK system lacked the flexibility of the modern Perfexion, Icon, and Esprit models. Rather than controlling each sector individually, large single-size collimator helmets were required to be manually plugged prior to shot delivery, resulting in treatment plans composed almost entirely of what modern systems now refer to as “classic shots” [10,11,12]. To enhance optimization routines in treatment planning, algebraic approximations leveraged this simpler delivery method of older GK systems by first approximating the dose distribution as spherically symmetric [13]. These approximations were then improved by shifting from spherical to ellipsoidal distributions [14]. These endeavors into intermediate dose calculations have allowed researchers to bypass the closed nature of the GK system while developing the original inverse optimization routine adopted by Elekta for GK [15]. While these geometric methods remain useful for approximating dose distributions in classic shots, the independent sector activations introduced by the modern models complicate the process of approximating their resultant dose distributions.
To address this enhanced complexity, dose calculation efforts have also explored the option of modeling all sources and collimation options individually, by either summing their contributions or via Monte Carlo simulations [16,17,18,19]. In this work, we developed an alternative method capable of producing clinically relevant dose approximations for modern GK systems utilizing a dose distribution superposition approach that relies on fewer individual models without the computational cost of Monte Carlo. This tool provides a simple approach to calculating doses based on key machine parameters—patient position, sector activations, and beam-on time—enabling an intermediate dose calculation external to the closed GK system. With this tool, we can batch process any arbitrary plan parameter combinations, enhancing our understanding of the dose distributions produced by planning optimizations and automation approaches in GK treatment planning.

2. Materials and Methods

2.1. Data Acquisition and Preparation

Dose arrays and their associated machine parameters were exported from the Leksell GammaPlan (Leksell GammaPlan, Version 11, Elekta AB, Stockholm, Sweden) treatment planning system (TPS) for the Gamma Knife Icon unit. Each dose array consists of a 448 × 448 × 448 array with a pixel spacing of 0.5 mm per voxel. The GammaPlan TMR10 algorithm uses a homogenous dose calculation, where the patients’ anatomy is considered as water. TMR10 accounts for the divergence of each Co-60 beamlet with the inverse square law, exponential attenuation of each beamlet as they traverse the assumed water environment, differences in the dose deposition rate between collimation sizes with output factors, and depth variations with a scaling distance. Any heterogeneities within the volume bounded by the external boundary are not accounted for with the standard TMR10 dose calculation [20]. This information is also available in documents typically associated with the Written Directive for the treatment unit, such as the GK Details Document and the Operator’s Report.
A dataset of 30 anonymized patients, previously treated with the GK system at the Mary Bird Perkins Cancer Center/Our Lady of the Lake, was curated for this study. This cohort consisted of patients with a variety of lesions treatable by GK, including 15 patients receiving treatment for multiple metastases, 5 for meningiomas, 3 for single metastases, 3 for pituitary adenomas, 3 for glioblastomas, and 1 for an oligodendroglioma [9]. The total target number among this patient cohort was 71. The physician-defined contours for each of these targets were also exported for the calculation of deliverability metrics.
A standard dose distribution base set for each collimation setting of each sector was calculated within GammaPlan using the Radiosurgery Head Phantom CIRS Model 605. The planning position was set as the (100, 100, 100) point in the GammaPlan coordinate system, corresponding to an axial plane 80 mm inferior to the apex of the phantom’s external contour. The position was roughly centered in the axial slice along the GammaPlan-measured lengths of 202 mm in the anterior–posterior direction and 156 mm in the medial–lateral direction. This base set is composed of three sizes (4, 8, and 16 mm) for each of the eight sectors, resulting in 24 single-sector dose arrays. After exporting from GammaPlan, the resulting arrays were cropped to 148 × 148 × 148 voxels for the 4 mm and 8 mm activations, and 224 × 224 × 224 voxels for the 16 mm activations. This was performed to optimize storage and processing demands, and the size of each cropped array was chosen to maintain the sub-10% low-dose region.

2.2. Sector Activation Superposition

An empty array is initialized to the same dimensions and spacing as the CBCT used to establish the patient’s stereotactic treatment space. This array is 448 × 448 × 448 voxels with 0.5 mm/voxel spacing in all three dimensions. The origin coordinate for the stereotactic treatment space is the ImagePositionPatient metadata entry for the first slice of the patient’s CBCT. Each patient’s target is accessed sequentially. For each shot in a target, the shot’s position relative to the origin coordinate is calculated. The routine then loops through each of the shot’s sector activations, retrieving the corresponding dose contributions from the base set. These contributions are scaled by their shot-weighting parameter and superimposed into the global patient array at the calculated position relative to the origin coordinate of the stereotactic treatment space. This process is repeated for all shots in the current target before moving to the next target, if the patient was treated for multiple targets.
A mathematical representation of the accumulation is provided below:
d [ x , y , z , w , S ] = j = 1 m w j i = 1 8 K j S i ,   x j ,   y j ,   z j
where   K j S i = k i ,   s m a l l k i ,   m e d i u m k i ,   l a r g e 0 ,     S i = 4 ,     S i = 8 ,     S i = 16 ,     S i = 0
D D m a x = t = 1 n d t × D m a x , t m a x ( d t )
where the process is presented in two steps. The first step, d([(x, y, z), w, S]), is an accumulation of a single target’s dose distribution as a function of all shots, indexed by j, in the target with individual center positions (xj, yj, zj), shot weighting (wj), and the eight sector activations Kj(Si). A dose is accumulated for each shot one sector at a time, from i = 1 to i = 8, by summing the appropriate dose contribution kernel from the basis set at the (xj, yj, zj) position of the accumulation array. This is repeated for each individual shot, from j = 1 to j = m, where m is the total number of shots in the target’s plan. The second step then accumulates all target contributions, indexed by t, after normalizing the dose distribution for the target and scaling it to units of Gray. The scaling value, Dmax, is equal to the prescription dose [Gy] divided by the fractional planning isodose level.
Using Figure 1 as an example, the inputs required to achieve Figure 1B would be d([(70,70,70), 1, (8,8,8,8,8,8,8,8)]) and D(40), where the (70,70,70) position is the center of the array in voxel coordinates. Similarly, Figure 1C would correspond to inputs of d([(70,70,70), 1, (8,4,0,16,0,4,8,4)]) and D(40). Finally, Figure 1D would require the information of both shots separated by a distance of 1 cm (20 voxels) along the x-axis, corresponding to inputs of d([(60,70,70), 1, (8,8,8,8,8,8,8,8)], [(80,70,70), 1, (8,4,0,16,0,4,8,4)]) and D(40).
After producing the superposition approximation (SA) dose arrays for each of the 30 patients, the original and recreated full-size arrays were cropped into additional target-centered arrays. This resulted in 71 arrays, each cropped to match the dimensions of the target box used during the original GammaPlan treatment planning. These arrays were compared as absolute doses or normalized to their maximum for relative dose comparisons. Both approaches are illustrated in Figure 2, which provides an example of using the SA method to recreate an entire target’s 3D dose. Figure 2A,B share a color bar that represents absolute dose in Gy. Figure 2A shows the dose distribution for a meningioma as calculated with TMR10. Figure 2B shows the result of the SA method for the same target. The black isodose lines in both figures represent the 52% prescription isodose volume. Pairwise comparisons and deliverability metrics were then calculated using these binary volumes. Figure 2C overlaps those volumes along the target’s central axial slice, with the TMR10 calculation shown in green, the SA method in magenta, and their intersection in white.

2.3. Approximation Evaluation

Dose arrays were approximated for all 71 targets from the 30-patient cohort. The performance of the SA was evaluated by comparing these approximations to the original TMR10 dose arrays obtained from the GammaPlan system. Due to the known limitation of the SA not correcting for depth differences, there is an inherent inability to replicate interactions between shots that would produce hotspots. In GK-SRS, a measure of dose heterogeneity within the target is not typically used to evaluate plan quality. Clinical evidence has shown that the presence of these hotspots within the target volume provides a radiobiological advantage, and penalizing these hotspots in planning optimization negatively impacts dosimetric fall-off around the target [3,4]. To conform with this, evaluation methods like gamma analysis, which would penalize differences in heterogeneous regions, were not used. Instead, our focus was to evaluate the clinical performance of the SA using two methods, especially in the 50–60% isodose level commonly used as the prescription line in planning. First, the pairwise comparisons of the Dice Similarity Coefficient (DSC) and Hausdorff Distance (HD) were each calculated at every 10% threshold from 10% to 90%, where higher thresholds provide insight into differences between hotspot regions. The second method is to compare the deliverability metrics that are commonly used in a clinical setting to evaluate plan quality. Our goal with both approaches is to demonstrate the similarity between the SA and the original TMR10 calculation.

2.3.1. Pairwise Comparisons

The Dice Similarity Coefficient (DSC), first proposed by Dice [21], is a measure of overlap between two sets. It has commonly been used to evaluate image segmentations [22,23] and is defined for a 3D implementation in this work as
D S C = 2 | V T M R 10 V S A | V T M R 10 + V S A
where VTMR10 and VSA refer to the binary volumes produced by the GammaPlan TMR10 calculation and SA, respectively. The DSC ranges from 0 to 1, with high values indicating higher similarity between the sets. The metric is sensitive to differences between the sets in both size and position [24]. Figure 2C shows two sets that achieved a high 3D DSC of 0.939 at the prescription line used in planning, indicating very strong agreement between the two dose distributions for this target.
The Hausdorff Distance (HD) provides a measure of the maximum distance between two sets [25]. It is defined as
H D = m a x ( d x , S A   ,   d y , T M R 10   )
where x and y represent a point along the surface of the volume created by thresholding the TMR10 and SA dose arrays, respectively. Following this, d(x, SA) and d(y, TMR10) represent the minimum distance from a point x or y to a point on the surface of the thresholded SA or TMR10 volume, respectively. Essentially, HD represents the largest minimum distance between the surfaces of the SA and TMR10 distributions when each is thresholded to the same percentage.
In this study, we used a 95% Hausdorff Distance (HD95), which reports the 95th percentile of the distances between points in two sets, rather than the maximum distances. With the exclusion of the largest 5% of distances, HD95 is more robust and outlier-resistant than the traditional HD metric [26].

2.3.2. Deliverability Metrics

The GammaPlan TPS reports three metrics commonly used to evaluate the quality of a GK treatment plan. These are the coverage, selectivity, and gradient index [27], each of which is defined as follows:
C o v e r a g e = T T V G T V
S e l e c t i v i t y = T T V P I V
G r a d i e n t   I n d e x = P I V 50 %   P I V
where PIV refers to the prescription isodose volume, GTV refers to the physician-defined gross target volume, TTV refers to the treated target volume (or the fraction of GTV receiving PIV), and PIV50% refers to the volume bounded by the isodose line at 50% of the prescription dose. Coverage describes how well the prescribed dose covers the target and ranges from 0 to 1. It has an ideal value of 1, with lower values representing an underdosing of the TV. Selectivity describes how well the prescribed dose is limited to the TV. It shares a range and ideal value with coverage, with lower values representing an overdosing of healthy tissue surrounding the TV. The gradient index quantifies the steepness of the dose gradient, specifically from the prescription dose to 50% of that value. Lower values indicate a steeper dose gradient and better sparing of healthy tissues.

3. Results

3.1. Results of Pairwise Comparisons

Thresholds were applied to the normalized dose arrays in 10% increments of the isodose volumes for each of the 71 targets, starting from 10% and continuing to 90% of the maximum dose. Each target was also thresholded to its specific prescription isodose value. Using these thresholded binary volumes, the DSC and HD95 were calculated between the TMR10-calculated isodose volume and the corresponding SA isodose volume. Figure 3 displays the average and standard deviation in the DSC and HD95 for all 71 targets at each threshold level. The left axis corresponds to DSC values, indicated by the black circles and bars, while the right axis corresponds to HD95 values, indicated by the cyan diamonds and bars. Threshold levels ranging from 20 to 60% demonstrate strong agreement in the DSC, with values above 0.85 (indicated by a black dashed line). In other applications, such as in evaluating the performance of auto-segmentation routines, a DSC > 0.70 is generally considered to indicate a good overlap between two volumes [22,24]. This range includes the prescription isodose used for treatment planning, denoted as “Rx”, where a typical Gamma Knife target is prescribed 18–24 Gy to an isodose line between 50 and 60% [5]. The HD95 values also show the strongest agreement in the 50–60% range, with a maximum discrepancy of 2.18 mm. At the prescription (Rx) isodose level, the maximum HD95 value reaches 2.24 mm, marked by a cyan dashed line in Figure 3.

3.2. Results of Deliverability Metrics

The same prescription isodose binary volumes used to generate the DSC and HD95 values in the ’Rx’ column of Figure 3 were also used to calculate the deliverability metrics of coverage, selectivity, and gradient index for both the original and approximated dose arrays. For coverage and selectivity, the cropped and binarized dose arrays were compared to their respective physician-defined target contours. To calculate the gradient index, additional binaries were generated at 50% of the prescribed dose level (25% of the maximum dose). Figure 4 (top row) illustrates a 2D histogram comparison of these deliverability metrics between the GammaPlan calculation and the SA. Coverage shows a close agreement between the datasets, with a TMR10 mean value of 0.973 (±0.021) and an SA mean of 0.959 (±0.044). The selectivity (0.674 ± 0.193 vs. 0.666 ± 0.202) and gradient index (3.218 ± 0.667 vs. 3.101 ± 0.642) also demonstrate strong similarity. The distributions of differences between the two calculations are included as 1D histograms in Figure 4 (bottom row). The mean and standard deviation for the coverage difference (0.014 ± 0.035), selectivity difference (0.008 ± 0.037), and gradient index difference (0.118 ± 0.192) were also calculated.

4. Discussion

We propose a method for recreating Gamma Knife Icon dose distributions using treatment plan parameters. By superimposing and summing sector-specific dose arrays from the base set into the positions defined within those plan parameters, dose contributions for any shot can be accumulated into a full-sized 448 × 448 × 448 array. Dose arrays created with this method were compared with the original TMR10-calculated dose arrays using pairwise comparisons for similarity between binarized arrays and the deliverability metrics reported using GammaPlan during the planning process.
The superposition dose approximation method produced dose arrays closely resembling the TMR10 distributions, especially within the isodose normalization ranges commonly used in treatment planning (50–60%). In this isodose range, the DSCs exceeded 85% and the HD95s were below 2.18 mm. A similar pattern was seen at the prescription (Rx) isodose level, with a maximum HD95 value of 2.24 mm (Figure 3). This strong agreement is further supported by the similarity in deliverability metrics between the two dose array sets, with minimal average differences observed (Figure 4, bottom).
Differences in low- and very high-dose regions arise primarily from two factors: the absence of depth correction in the SA and the exclusion of very low-dose regions from the base set arrays. GammaPlan uses the TMR10 dose calculation algorithm to account for attenuation, scatter, and beamlet divergence across the volume within the external boundary [18]. In contrast, our method superimposes sector activations based on the same reference depth for each sector in the base set. This introduces a limitation to the superposition approximation, since regions with dose levels significantly above or below the Rx dose may not accumulate accurately.

5. Conclusions

In conclusion, we introduce a new method for approximating Gamma Knife dose distributions external to the GammaPlan treatment planning software. This method works by superimposing sector-specific dose distributions into the correct positions, using the machine parameters associated with the original treatment plan. Our evaluation demonstrates that this method yields dose distributions comparable to those produced by GammaPlan within clinically relevant dose levels, though some differences were noted in the positioning and size of hotspots and low-dose regions. This tool offers a simple calculation approach that enables the external calculation of dose distributions given the plan parameters of a GK treatment and lays the groundwork for future automation approaches to GK treatment planning. In future work, this dose approximation tool can be deployed to auto-generate simulated GK plan data and associated dose distributions for use in artificial intelligence development for prediction models of plan parameters, where large datasets would require batch processing.

Author Contributions

Conceptualization, P.H.S., J.D., C.C., D.N., and D.S.; methodology, P.H.S., R.R., C.C., D.N., and D.S.; software, P.H.S., L.M.L., R.R., and D.S.; validation, P.H.S., L.M.L., and D.S.; formal analysis, P.H.S., L.M.L., and D.S.; investigation, P.H.S., L.M.L., R.R., C.C., D.N., J.D., and D.S; data curation, P.H.S., R.R., and D.S.; writing—original draft preparation, P.H.S., L.M.L., and D.S.; writing—review and editing, P.H.S., L.M.L., J.D., C.C., D.N., and D.S.; visualization, P.H.S., L.M.L., and D.S.; project administration, D.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported in part by the Dr. Charles M. Smith Medical Physics Research Endowment. This work was also partially funded by the National Science Foundation through the REU Site in Physics & Astronomy, NSF grant #1852356, at Louisiana State University.

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Institutional Review Board of Louisiana State University (protocol code E11732, approved on 30 May 2019 as Exempt, category 4b).

Informed Consent Statement

Not applicable because of the retrospective nature and the absence of patient risk in this research. All potentially identifiable information has been anonymized, in line with the above IRB exemption.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

We extend our gratitude to the Mary Bird Perkins Cancer Center for their ongoing support and to Jonas Fontenot for facilitating the IRB exemption.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Li, Y.; Rodrigues, A.; Li, T.; Yuan, L.; Yin, F.F.; Wu, Q.J. Impact of dose calculation accuracy during optimization on lung IMRT plan quality. J. Appl. Clin. Med. Phys. 2015, 16, 5137. [Google Scholar] [CrossRef]
  2. Reynolds, T.A.; Jensen, A.R.; Bellairs, E.E. Ozer Dose gradient index for stereotactic radiosurgery/radiation therapy. Int. J. Radiat. Oncol. Biol. Phys. 2020, 106, 604–611. [Google Scholar] [CrossRef] [PubMed]
  3. Brown, M.H.; Marcrom, S.R.; Patel, M.P.; Popple, R.A.; Travis, R.L.; McDonald, A.M.; Riley, K.O.; Markert, J.M.; Willey, C.D.; Bredel, M.; et al. Understanding the Effect of Prescription Isodose in Single-Fraction Stereotactic Radiosurgery on Plan Quality and Clinical Outcomes for Solid Brain Metastases. Neurosurgery 2023, 93, 1313–1318. [Google Scholar] [CrossRef] [PubMed]
  4. Kondziolka, D.; Shin, S.M.; Brunswick, A.; Kim, I.; Silverman, J.S. The biology of radiosurgery and its clinical applications for brain tumors. Neuro-Oncology 2015, 17, 29–44. [Google Scholar] [CrossRef]
  5. Ashley Wackym, P.; Runge-Samuelson, C.L.; Friedland, D.R. Stereotactic Radiosurgery of Skull Base Tumors. In Otologic Surgery, 3rd ed.; Brackmann, D.E., Shelton, C., Arriaga, M.A., Eds.; W.B. Saunders: Philadelphia, PA, USA, 2010; pp. 785–798. [Google Scholar]
  6. Boling, W.; Song, M.; Shih, W.; Karlsson, B. Gamma Knife Radiosurgery for Trigeminal Neuralgia: A Comparison of Dose Protocols. Brain Sci. 2019, 9, 134. [Google Scholar] [CrossRef] [PubMed]
  7. Liu, H.; Andrews, D.W.; Evans, J.J.; Werner-Wasik, M.; Yu, Y.; Dicker, A.P.; Shi, W. Plan Quality and Treatment Efficiency for Radiosurgery to Multiple Brain Metastases: Non-Coplanar RapidArc vs. Gamma Knife. Front. Oncol. 2016, 6, 26. [Google Scholar] [CrossRef] [PubMed]
  8. Schell, M.C.; Bova, F.J.; Larson, D.A.; Leavitt, D.D.; Lutz, W.R.; Podgorsak, E.B.; Wu, A. Stereotactic Radiosurgery: Report of Task Group 42, Radiation Therapy Committee; American Institute of Physics: Woodbury, New York, USA, 1995. [Google Scholar]
  9. Desai, R.; Rich, K.M. Therapeutic Role of Gamma Knife Stereotactic Radiosurgery in Neuro-Oncology. Mo. Med. 2020, 117, 33–38. [Google Scholar] [PubMed]
  10. Lindquist, C.; Paddick, I. The Leksell Gamma Knife Perfexion and comparisons with its predecessors. Neurosurgery 2007, 61 (Suppl. 3), 130–140; discussion 140-1. [Google Scholar] [PubMed]
  11. Novotny, J.; Bhatnagar, J.P.; Niranjan, A.; Quader, M.A.; Huq, M.S.; Bednarz, G.; Flickinger, J.C. Dosimetric comparison of the Leksell Gamma Knife Perfexion and 4C. J. Neurosurg. 2008, 109, 8–14. [Google Scholar] [CrossRef]
  12. Asgari, S.; Banaee, N.; Nedaie, H.A. Comparison of full width at half maximum and penumbra of different Gamma Knife models. J. Cancer Res. Ther. 2018, 14, 260–266. [Google Scholar] [CrossRef] [PubMed]
  13. Shepard, D.M.; Ferris, M.C.; Ove, R.; Ma, L. Inverse treatment planning for Gamma Knife radiosurgery. Med. Phys. 2000, 27, 2748–2756. [Google Scholar] [CrossRef] [PubMed]
  14. Ferris, M.; Lim, J.; Shepard, D. Radiosurgery Treatment Planning via Nonlinear Programming. Ann. Oper. Res. 2003, 119, 247–260. [Google Scholar] [CrossRef]
  15. Shepard, D.M.; Chin, L.S.; DiBiase, S.J.; Naqvi, S.A. Clinical implementation of an automated planning system for gamma knife radiosurgery. Int. J. Radiat. Oncol. Biol. Phys. 2003, 56, 1488–1494. [Google Scholar] [CrossRef] [PubMed]
  16. Li, K.; Ma, L. Selective source blocking for Gamma Knife radiosurgery of trigeminal neuralgia based on analytical dose modelling. Phys. Med. Biol. 2004, 49, 3455–3463. [Google Scholar] [CrossRef] [PubMed]
  17. Moskvin, V.; Timmerman, R.; DesRosiers, C.; Randall, M.; DesRosiers, P.; Dittmer, P.; Papiez, L. Monte Carlo simulation of the Leksell Gamma Knife®: II. Effects of heterogeneous versus homogeneous media for stereotactic radiosurgery. Phys. Med. Biol. 2004, 49, 4879–4895. [Google Scholar]
  18. Xiong, W.; Huang, D.; Lee, L.; Feng, J.; Morris, K.; Calugaru, E.; Burman, C.; Li, J.; Ma, C.M. Implementation of Monte Carlo Simulations for the Gamma Knife System. J. Phys. Conf. Ser. 2007, 74, 021023. [Google Scholar] [CrossRef]
  19. Yuan, J.; Lo, S.S.; Zheng, Y.; Sohn, J.W.; Sloan, A.E.; Ellis, R.; Machtay, M.; Wessels, B. Development of a Monte Carlo model for treatment planning dose verification of the Leksell Gamma Knife Perfexion radiosurgery system. J. Appl. Clin. Med. Phys. 2016, 17, 190–201. [Google Scholar] [CrossRef] [PubMed]
  20. Xu, A.; Bhatnagar, J.; Bednarz, G.; Niranjan, A.; Flickinger, J.; Lunsford, L.D.; Huq, M.S. Dose differences between the three dose calculation algorithms in Leksell GammaPlan. J. Appl. Clin. Med. Phys. 2014, 15, 4844. [Google Scholar] [CrossRef]
  21. Dice, L.R. Measures of the amount of ecologic association between species. Ecology 1945, 26, 297–302. [Google Scholar] [CrossRef]
  22. Zou, K.H.; Warfield, S.K.; Bharatha, A.; Tempany, C.M.C.; Kaus, M.R.; Haker, S.J.; Wells, W.M., 3rd; Ferenc, A.J.; Kikinis, R. Statistical validation of image segmentation quality based on a spatial overlap index. Acad. Radiol. 2004, 11, 178–189. [Google Scholar] [CrossRef] [PubMed]
  23. Taha, A.A.; Hanbury, A. Metrics for evaluating 3D medical image segmentation: Analysis, selection, and tool. BMC Med. Imaging 2015, 15, 29. [Google Scholar] [CrossRef] [PubMed]
  24. Zijdenbos, A.P.; Dawant, B.M.; Margolin, R.A.; Palmer, A.C. Morphometric analysis of white matter lesions in MR images: Method and validation. IEEE Trans. Med. Imaging 1994, 13, 716–724. [Google Scholar] [CrossRef] [PubMed]
  25. Porter, E.; Solis, D.; Bruckmeier, P.; Siddiqui, Z.A.; Zamdborg, L.; Guerrero, T. Effect of Loss Functions in Deep Learning-Based Segmentations. In Auto-Segmentation for Radiation Oncology State of the Art; Yang, J., Sharp, G.C., Gooding, M.J., Eds.; CRC Press: Boca Raton, FL, USA, 2021; pp. 133–149. [Google Scholar]
  26. Chollet, F. Deeplearning with Python; Simonand Schuster: New York, NY, USA, 2021. [Google Scholar]
  27. Torrens, M.; Chung, C.; Chung, H.T.; Hanssens, P.; Jaffray, D.; Kemeny, A.; Yu, C.P. Standardization of terminology in stereotactic radiosurgery: Report from the Standardization Committee of the International Leksell Gamma Knife Society: Special topic. J. Neurosurg. 2014, 121 (Suppl. 2), 2–15. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Example sector activations and shot interactions created using the superposition method. (A) is an example of only sector 1 being active to 8 mm, as indicated by the key in the top right of the image. (B) shows all 8 sectors activated to the same collimation setting of 8 mm, forming a classic 8 mm shot. In (C), an example shot with mixed activations is displayed. (D) shows the interaction of the shots in (B,C) when separated by a horizontal distance of 1 cm, with the 50% isodose line shown in black. “x” and “+” indicate the shot positions for the shot shown in (B,C), respectively.
Figure 1. Example sector activations and shot interactions created using the superposition method. (A) is an example of only sector 1 being active to 8 mm, as indicated by the key in the top right of the image. (B) shows all 8 sectors activated to the same collimation setting of 8 mm, forming a classic 8 mm shot. In (C), an example shot with mixed activations is displayed. (D) shows the interaction of the shots in (B,C) when separated by a horizontal distance of 1 cm, with the 50% isodose line shown in black. “x” and “+” indicate the shot positions for the shot shown in (B,C), respectively.
Radiation 05 00006 g001
Figure 2. An example of the full target reconstruction of the dose array using the superposition method. (A) is the TMR10-produced dose array for the treatment of a meningioma. (B) shows the superposition approximation. The black isodose line shown in both (A,B) highlights the prescription isodose of 52% for each dose distribution. (C) highlights differences between the prescription iso-volumes of (A,B) (shown in black) as an RGB formatted fusion, where green is the TMR10-calculated iso-volume, magenta is the superposition approximation iso-volume, and white is their intersection. The shown Dice Similarity Coefficient is the result of a 3D comparison for this target with these prescription volumes.
Figure 2. An example of the full target reconstruction of the dose array using the superposition method. (A) is the TMR10-produced dose array for the treatment of a meningioma. (B) shows the superposition approximation. The black isodose line shown in both (A,B) highlights the prescription isodose of 52% for each dose distribution. (C) highlights differences between the prescription iso-volumes of (A,B) (shown in black) as an RGB formatted fusion, where green is the TMR10-calculated iso-volume, magenta is the superposition approximation iso-volume, and white is their intersection. The shown Dice Similarity Coefficient is the result of a 3D comparison for this target with these prescription volumes.
Radiation 05 00006 g002
Figure 3. Average Dice Similarity Coefficient (DSC) and 95th percentile Hausdorff Distance (HD95) comparisons between the original (TMR10-produced) dose arrays and the superposition approximation (SA) dose arrays for 71 targets at 10 threshold levels. Both metrics are reported at each 10% isodose threshold between 10% and 90%, as well as at each target’s respective prescription isodose level, denoted as “Rx”. The DSC, aligned with the left vertical axis, is shown using black circles and error bars, with a black dashed line marking the 0.85 DSC level. HD95, aligned with the right vertical axis, is illustrated with cyan diamonds and error bars, with a cyan dashed line marking the maximum HD95 value at the prescription (Rx) level. Higher DSC values and lower HD95 scores indicate a closer match between the SA-produced dose and the TMR10-produced dose.
Figure 3. Average Dice Similarity Coefficient (DSC) and 95th percentile Hausdorff Distance (HD95) comparisons between the original (TMR10-produced) dose arrays and the superposition approximation (SA) dose arrays for 71 targets at 10 threshold levels. Both metrics are reported at each 10% isodose threshold between 10% and 90%, as well as at each target’s respective prescription isodose level, denoted as “Rx”. The DSC, aligned with the left vertical axis, is shown using black circles and error bars, with a black dashed line marking the 0.85 DSC level. HD95, aligned with the right vertical axis, is illustrated with cyan diamonds and error bars, with a cyan dashed line marking the maximum HD95 value at the prescription (Rx) level. Higher DSC values and lower HD95 scores indicate a closer match between the SA-produced dose and the TMR10-produced dose.
Radiation 05 00006 g003
Figure 4. (Top row) Two-dimensional histograms for the deliverability metric scores achieved by each dose distribution, with the red line representing 1-to-1 correlation. The x-axis corresponds to the original TMR10-produced dose array, while the y-axis corresponds to the SA-produced dose array. (Bottom row) One-dimensional histograms for the target-to-target difference for each of the three metrics. For coverage (A) and selectivity (B), binning was selected as every 1%. For the gradient index (C), bins were assigned every increment of 0.05.
Figure 4. (Top row) Two-dimensional histograms for the deliverability metric scores achieved by each dose distribution, with the red line representing 1-to-1 correlation. The x-axis corresponds to the original TMR10-produced dose array, while the y-axis corresponds to the SA-produced dose array. (Bottom row) One-dimensional histograms for the target-to-target difference for each of the three metrics. For coverage (A) and selectivity (B), binning was selected as every 1%. For the gradient index (C), bins were assigned every increment of 0.05.
Radiation 05 00006 g004
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Stone, P.H.; Lay, L.M.; Ramirez, R.; Neck, D.; Chu, C.; Dey, J.; Solis, D. A 3D Superposition Approximation for Gamma Knife Dose Calculation. Radiation 2025, 5, 6. https://doi.org/10.3390/radiation5010006

AMA Style

Stone PH, Lay LM, Ramirez R, Neck D, Chu C, Dey J, Solis D. A 3D Superposition Approximation for Gamma Knife Dose Calculation. Radiation. 2025; 5(1):6. https://doi.org/10.3390/radiation5010006

Chicago/Turabian Style

Stone, Payton H., Lam M. Lay, Raymi Ramirez, Daniel Neck, Connel Chu, Joyoni Dey, and David Solis. 2025. "A 3D Superposition Approximation for Gamma Knife Dose Calculation" Radiation 5, no. 1: 6. https://doi.org/10.3390/radiation5010006

APA Style

Stone, P. H., Lay, L. M., Ramirez, R., Neck, D., Chu, C., Dey, J., & Solis, D. (2025). A 3D Superposition Approximation for Gamma Knife Dose Calculation. Radiation, 5(1), 6. https://doi.org/10.3390/radiation5010006

Article Metrics

Back to TopTop