

TECHNICAL NOTE 



Year : 2015  Volume
: 40
 Issue : 4  Page : 240245 

Derivative based sensitivity analysis of gamma index
Biplab Sarkar^{1}, Anirudh Pradhan^{2}, T Ganesh^{3}
^{1} Department of Radiation Oncology, Fortis Memorial Research Institute, Gurgaon, Haryana; Department of Physics, GLA University, Mathura, Uttar Pradesh, India ^{2} Department of Mathematics, GLA University, Mathura, Uttar Pradesh, India ^{3} Department of Radiation Oncology, Fortis Memorial Research Institute, Gurgaon, Haryana, India
Date of Submission  20Jun2015 
Date of Decision  27Sep2015 
Date of Acceptance  28Sep2015 
Date of Web Publication  1Dec2015 
Correspondence Address: Biplab Sarkar Department of Radiation Oncology, Fortis Memorial Research Institute, Gurgaon, Haryana India
Source of Support: None, Conflict of Interest: None  Check 
DOI: 10.4103/09716203.170789
Abstract   
Originally developed as a tool for patientspecific quality assurance in advanced treatment delivery methods to compare between measured and calculated dose distributions, the gamma index (γ) concept was later extended to compare between any two dose distributions. It takes into effect both the dose difference (DD) and distancetoagreement (DTA) measurements in the comparison. Its strength lies in its capability to give a quantitative value for the analysis, unlike other methods. For every point on the reference curve, if there is at least one point in the evaluated curve that satisfies the pass criteria (e.g., δDD = 1%, δDTA = 1 mm), the point is included in the quantitative score as "pass." Gamma analysis does not account for the gradient of the evaluated curve  it looks at only the minimum gamma value, and if it is <1, then the point passes, no matter what the gradient of evaluated curve is. In this work, an attempt has been made to present a derivativebased method for the identification of dose gradient. A mathematically derived reference profile (RP) representing the penumbral region of 6 MV 10 cm × 10 cm field was generated from an error function. A general test profile (GTP) was created from this RP by introducing 1 mm distance error and 1% dose error at each point. This was considered as the first of the two evaluated curves. By its nature, this curve is a smooth curve and would satisfy the pass criteria for all points in it. The second evaluated profile was generated as a sawtooth test profile (STTP) which again would satisfy the pass criteria for every point on the RP. However, being a sawtooth curve, it is not a smooth one and would be obviously poor when compared with the smooth profile. Considering the smooth GTP as an acceptable profile when it passed the gamma pass criteria (1% DD and 1 mm DTA) against the RP, the first and second order derivatives of the DDs (δD', δD") between these two curves were derived and used as the boundary values for evaluating the STTP against the RP. Even though the STTP passed the simple gamma pass criteria, it was found failing at many locations when the derivatives were used as the boundary values. The proposed derivativebased method can identify a noisy curve and can prove to be a useful tool for improving the sensitivity of the gamma index.
Keywords: Derivative, dose distribution comparison, gamma index, radiotherapy, sensitivity, sensitivity analysis
How to cite this article: Sarkar B, Pradhan A, Ganesh T. Derivative based sensitivity analysis of gamma index. J Med Phys 2015;40:2405 
Introduction   
Modern day radiation treatment delivery often involves highly complex dose distributions coupled with a large number of monitoring units from the linear accelerator. This necessitates patient specific quality assurance (QA) tests for every treatment plan before its execution. Planar dose verification in the form of comparison between measured dose distribution using a planar detector and calculated dose distribution from the treatment planning system (TPS) is the highly preferred method for patientspecific QA. Although there are several ways of comparing the two dose distributions such as simple visual inspection, superpositioning of isodose lines, the gamma index method has become the most popular comparison method since it incorporates the quantitative element in the analysis besides the qualitative aspect. The method, initially introduced to compare measured and calculated dose distributions, has been extended to make a comparison between any two dose distributions  measured versus measured, measured versus calculated, and calculated versus calculated. Thus, its application has extended beyond simple patient specific QA in intensity modulated radiation therapy, volumetric modulated arc therapy (VMAT) to the commission of TPSs and beam matching between linear accelerators.^{[1]}
Earlier evaluation of the dose distribution for TPS was calculated by superimposing the isodose distributions, either manually with the isodose distributions on a light box or by software tools. This method was purely qualitative rather than quantitative. van Dyk et al. in 1993, proposed a system wherein the isodose distribution was subdivided into regions of high and lowdose gradients, each with a different acceptance criterion.^{[2]} Vensalaar et al. described another way of qualitative evaluation of a dose distribution and proposed acceptance criteria along the lines proposed by van Dyk et al.^{[3]} In both the abovementioned methods, the doses are compared directly in the lowgradient regions, with a specified limit of the difference between the measured and calculated doses.^{[2],[3]} In the highgradient region such as penumbra, Venselaar et al. proposed that the acceptable deviation in the dose could be as high as 30–50% depending upon the complexity of the geometry.^{[3]} Comparison of dose difference (DD) is avoided in the highgradient region since a small spatial error between test and reference distribution will lead to a large difference in dose. The concept of distancetoagreement (DTA) was introduced to determine the acceptability of the dose distribution in these regions.^{[4]} DTA is the distance between a reference data point and the nearest point in the test dose distribution which exhibits dose specified by %DD. It is meaningful to specify DTA for the highgradient region and percentage %DD for the lowgradient region. Therefore, when used concurrently, DTA and %DD are complementary to each other in the high and lowgradient region, respectively. The composite analysis using DTA and %DD was developed by Harms et al.^{[5]} which had its roots on the methodology proposed by Shiu et al.^{[6]} Cheng et al. applied this concept of using a passfail criterion incorporating both the DTA and DD.^{[7]} Each test point was investigated to establish if both the DTA and DD exceeded the specified tolerances. The inadequacy of this methodology was that it gave a binary distribution only. It did not give any unique numerical index which represented the degree of accuracy between the reference and evaluated distributions. In 1998, Low et al. proposed the gamma index (γ) method.^{[1]} This method still uses the concepts of %DD and DTA, but additionally it gives a quantitative value for each point in the evaluated distribution that is computed from the %DD, DTA, and the selected acceptance criteria. To date, it has remained as the most often used methodology for comparing two dose distributions. The techniques of onedimensional (1998) and twodimensional gamma analysis (2003) are wellstudied areas and have been studied in depth by several investigators.^{[1],[8],[9]} However, one known disadvantage of the gamma index method is that it is sensitive to normalization and does not recognize dose gradient, if at least a single point from test curve falls within the ellipsoid created by %DD and DTA around the reference point. Li et al. proposed a surfacebased method to recognize the dose gradient which the gamma index could not detect.^{[10]}
In radiotherapy use of derivatives is commonly employed and can be found in different applications: For example, a technique such as gradient search method is frequently used in inverse plan optimization during the treatment planning. Similarly, dose gradient concept is used in beam matching between two or more linear accelerators.^{[11]} In general, any gradient can be recognized and quantified using different order derivatives. Therefore, different order derivatives are useful tools to handle the short comings of the gamma index in terms of dose gradient.
The choice of the order of derivative (first order, second order, etc.) for evaluation of dose distribution depends on the dose gradient in the specific case. In a flat beam geometry, derivative up to second order is sufficient for the evaluation of dose distribution.^{[11]}
In this article, we present a theoretical formulation based on the first and second order dose derivatives that can help to distinguish between smooth curves and highly irregular noisy curves. While both the curves satisfy the gamma tolerance criteria, the derivatives show up unacceptable differences in the noisy curve. The method is limited to the onedimensional gamma analysis only.
Materials and Methods   
The method of obtaining the functional form of the penumbral region from the error function has been described by Low et al.^{[1]} In this work, we applied the same method using the fitting coefficient for our machine profile and the matched data.
Where, D(x_{i}) is the dose at any arbitrary position x_{i}; b_{1}, b_{2}, and a are the fitting coefficients, t is the collimator transmission, erf(x) is error function having value ranging between −1 and +1 at −∞ and +∞, respectively. Using the value a = 0.22, t = 0.02, b_{1} = 0.456 cm ^{–1}, b_{2} = 5 cm ^{–1}, the calculated profile from equation (i) and measured profile matched within 0.3 mm DTA and 0.3% DD. All the dose profiles used in this work were such calculated profiles or their derivatives. Throughout this work, a gamma pass criteria of δDTA = 1 mm and δDD = 1% was followed.
The method presented here uses a comparison between reference and evaluated (test) profile through their DDs and the first and second derivatives of dose in space. The reference profile (RP) was generated using equation (i), representing a typical profile of a 10 cm × 10 cm photon beam field from a linear accelerator. The profile has two distinct regions, a lowgradient region (the flat portion) and a highgradient region (the penumbral region). In the lowgradient region, DD between the reference and the evaluated (test) profile is important. Since dose derivatives barely have any role to play in this region, this region was not considered in this study. After removing the lowgradient region, midpoint on the highgradient part of the curve was chosen as the origin of a new coordinate system, and the entire curve was normalized at this point in such a way to have all other points in the curve within ±1.^{[1]} The chosen point divided the curve into a positive half and a negative half while perfectly preserving the functional nature of the curve.
The first derivative (D') in space with respect to x_{i} (called the dose slope) represents the spatial rate of variation of dose, and the second derivative (D"), called the curvature, gives the spatial rate of change of dose slope. [Figure 1] represents a dose profile, its first and second derivatives in space. For simplicity, only half of the profile is shown. The behavior of the derivative on the other half of the profile will be reciprocal. As shown in [Figure 1], the profile is subdivided into the first derivative significant area (slope) and second derivative significant area (curvature). First derivative will be significant around the region where x = 0 and second derivative in the bending parts of the profile where the steep slope region of the profile bends to merge with the plateau region. It is demarcated in [Figure 1] as the second derivative significant area. There is no clear demarcation between these first and second derivative significant areas and they overlap with each other. Hence, in the overlap region, both derivatives will be tested. More detailed characteristics of first and second derivatives of dose profile are described elsewhere.^{[11]}  Figure 1: Profile and its first and second derivative in space. In general, mathematics first derivative termed as slope, and the second derivative termed as curvature. The first derivative significant area and second derivative significant area is demarcated in the figure
Click here to view 
To investigate, whether the conventional gamma analysis between the RP and a sawtooth test profile (STTP) is an adequate test for quantitative evaluation of the dose matching between the two curves, we first compared the RP with a general test profile (GTP) using the first and second order spatial derivatives. Since the GTP was an acceptable smooth curve that satisfied the 1 mm DTA and 1% DD criteria against the RP, these spatial derivatives were used as boundary conditions for subsequent evaluation of STTP against RP.
The STTP was generated in such a way that it passed the conventional 1% DD and 1 mm DTA gamma criteria for all x_{i} as shown in [Figure 2]a when compared against RP. [Figure 2]b shows the minimum %DD for 0 mm and ±1 mm DTA search yielding γ ≤ ±1 for all points. Thus, the STTP is an acceptable curve when one goes only by the conventional gamma method.^{[1],[12]} However, it should ideally be not acceptable because of its sawtooth nature. The fact that such an unacceptable profile passes a stringent gamma criteria goes on to prove that the gamma index does not detect the slope of dose or positional changes if these are within certain limits.  Figure 2: (a) Reference profile and sawtooth test profile and their 0 mm, 1 mm, and −1 mm distancetoagreement search (b) Minimum dose difference for 0 mm, +1 mm, and −1 mm distancetoagreement search yielding γ ≤ ±1 for all points
Click here to view 
Results   
[Figure 3]a shows the RP and GTP, and the maximum and minimum variation of %DD (δD) between the two profiles. Similarly, [Figure 3]b and [Figure 3]c show first and second derivative of RP and its maximum and minimum variation of D' (δD') and D" (δD") differences, respectively. The profiles for maximum and minimum difference for D, D', and D" obtained by 0 and ±1 mm DTA search are shown in [Figure 3]a,[Figure 3]b,[Figure 3]c, respectively. Therefore, the test profile used here (GTP) is acceptable for the 1% DD and 1 mm DTA (γ ≤ 1) search criteria when tested against RP. Consequently, these derivatives can serve as limits for testing derivatives of any other profile satisfying γ ≤ 1 for limits of 1% DD and 1 mm DTA, the difference in dose profile, its slope and curvature with those of the RP should lie within these limits profiles.  Figure 3: (a) Maximum and minimum D difference between reference and general test profile for 0 mm, +1 mm, and −1 mm distancetoagreement (b) Maximum and minimum D' difference between reference and general test profile for 0 mm, +1 mm, and −1 mm distancetoagreement (c) Maximum and minimum D" difference between reference and general test profile for 0 mm, 1 mm, and −1 mm distancetoagreement
Click here to view 
The difference in dose (δD) and derivatives (δD' and δD") between these two curves (GTP and RP) give the upper and lower limits of the differences for dose and its derivatives as shown in [Figure 3]a,[Figure 3]b,[Figure 3]c. If the differences δD', δD" between the reference and the evaluated profiles are outside the limits specified in [Figure 3]b and [Figure 3]c, then the evaluated profile is not acceptable even though it may be passing the gamma test as in [Figure 3]a. To investigate the acceptability of STTP, the curve was examined with the derivatives as the upper and lower boundary limits. For the curve to pass, the first and second order derivatives of the STTP should fall within these limits. The results are shown in [Figure 4]a and [Figure 4]b. These figures clearly indicate the derivatives of RPSTTP combination fall beyond the limits derived from the RPGTP combination. From [Figure 4]a it is evident that the limits are exceeded at the following points: 1 mm, −2 mm, −5 mm, and −3 mm, −4 mm, 2 mm, 3 mm, respectively. Similarly from [Figure 4]b, it is evident that the limits are exceeded at 0 mm, −5 mm, −6 mm, 5 mm, and −3 mm, 4 mm, respectively. The results demonstrate that it is possible to identify large fluctuations in the evaluated profile by means of derivatives which was not otherwise possible with the simple gamma test based on DTA and %DD.  Figure 4: (a) D' differences for evaluated profile (reference profile + sawtooth test profile) in terms of limiting profile (reference profile + general test profile) for 0 mm, 1 mm, and −1 mm distancetoagreement (b) D" differences for evaluated profile (reference profile + sawtooth test profile) in terms of limiting profile (reference profile + general test profile) for 0 mm, 1 mm, and −1 mm distancetoagreement search
Click here to view 
Discussion   
Gamma analysis is a mathematically rigorous technique to compare two or more dose distributions in one or more dimensions. Matching of gamma index within a specified value is a necessary boundary condition; nevertheless, it is not the sufficient boundary condition for all situations, especially in evaluating a noisy curve. In our study, gamma analysis of an acceptable profile (the GTP) and an unacceptable profile (STTP) with respect to a mathematically generated RP resulted in a value of γ = 1 for both cases indicating both curves passed the test. However, the slope of the STTP suggested that it was not an acceptable dose profile. For steepest part of the test profile, difference in the first derivative indicated whether the dose profile was within the acceptable error limits or not. Similarly, in the bending portions of the profile where the slope changes rapidly, evaluation of D" at different points helped us to correctly identify whether the evaluated curve is an acceptable one or not.
The other drawback of gamma index is that it underestimates the results against no noise condition.^{[1]} Noise in dose distribution has a profound impact on gamma calculation. The impact depends on whether the noise is in the reference or evaluated distribution. If the evaluated distribution has noise, and if the reference distribution has no noise then it will be added to γ distribution in the same ratio to the normalized dose noise.^{[1]} It is a known fact that differentiation increases noise, the effect being predominant in the lowgradient region. The mathematical property of the derivative quantifies the gradient, which was exploited in this study. As shown in [Figure 2]a, the dose profile (signal) gets saturated at y = ±1, whereas the first derivative peaks at x = 0 where the rate of change of dose (dose slope) is the highest. The second derivative peaks at two positions, going up to ±6 where the dose curve bends (dose curvature), clearly indicating where the rate of change of slope is at its highest. The true characteristics of the dose curve at the said positions can thus only be determined uniquely by derivatives and not the simple gamma test.
Tremendous advancement in the radiotherapy equipment and delivery technique in recent times such as linear accelerator based stereotactic radiosurgery and radiotherapy using miniature MLC, VMAT, Tomotherapy, and CyberKnife make the dose distributions fall more sharply with a very highdegree of dose and fluence modulation, leading to dose escalation attempts for the tumor volume and better sparing of the normal tissues. It may not always be possible to quantify such highgradient distribution using the simple gamma passfail criteria. Physically, a dose distribution is continuous in character, but because of the limitation of our measuring devices only a discreet dose distribution can be obtained. Gamma analysis is also discreet in nature, and that is attributed to the inherent character of the dose measurement. It is possible to qualitatively identify a distribution within a specified limit by gamma analysis applied to its dose, first derivative and second derivative differences.
The two components of gamma index DTA and %DD are complementary to each other in high and lowgradient region, respectively. Derivative gamma index is effective in the highgradient region where DTA is prominent. The evaluation using a derivative gamma index is more stringent than using DTA. In the recent times, several modifications have been tested on the gamma index methodology.^{[10],[13]} Li et al. described a surfacebased distance method for gamma index. They established the relation between DTA and %DD using dose gradient factor, or the reciprocal of the mean dose gradient.^{[10]} In our study, we exploited similar characteristic of the dose distribution and gamma index. Li et al. also established a minimum detectable error of the gamma index with a given set of parameters using a step function.^{[10]} In our study, we used a hypothetical profile (STTP) to establish the insensitivity of the gamma index.
Further studies will be carried out to examine the usefulness of derivative gamma method for twodimensional cases. A study of the method's reliability and sensitivity to noise will also be investigated and presented in future articles.
Conclusion   
It is known that derivatives are the most suited mathematical functions to detect the slope and curvature. In derivativebased gamma evaluation method presented here, these characteristics were successfully used to identify the unacceptable nature of a noisy curve such as STTP which otherwise would have passed a routine gamma analysis. Use of derivative increases the minimum detectable error of gamma index.
Financial support and sponsorship
Nil
Conflicts of interest
There are no conflicts of interest.
References   
1.  Low DA, Harms WB, Mutic S, Purdy JA. A technique for the quantitative evaluation of dose distributions. Med Phys 1998;25:65661. 
2.  van Dyk J, Barnett RB, Cygler JE, Shragge PC. Commissioning and quality assurance of treatment planning computers. Int J Radiat Oncol Biol Phys 1993;26:26173. 
3.  Venselaar J, Welleweerd H, Mijnheer B. Tolerances for the accuracy of photon beam dose calculations of treatment planning systems. Radiother Oncol 2001;60:191201. 
4.  Hogstrom KR, Mills MD, Meyer JA, Palta JR, Mellenberg DE, Meoz RT, et al. Dosimetric evaluation of a pencilbeam algorithm for electrons employing a twodimensional heterogeneity correction. Int J Radiat Oncol Biol Phys 1984;10:5619. [ PUBMED] 
5.  Harms WB Sr, Low DA, Wong JW, Purdy JA. A software tool for the quantitative evaluation of 3D dose calculation algorithms. Med Phys 1998;25:18306. 
6.  Shiu AS, Tung S, Hogstrom KR, Wong JW, Gerber RL, Harms WB, et al. Verification data for electron beam dose algorithms. Med Phys 1992;19:62336. 
7.  Cheng A, Harms WB Sr, Gerber RL, Wong JW, Purdy JA. Systematic verification of a threedimensional electron beam dose calculation algorithm. Med Phys 1996;23:68593. 
8.  Low DA, Dempsey JF. Evaluation of the gamma dose distribution comparison method. Med Phys 2003;30:245564. 
9.  Low DA, Mutic S, Dempsey JF, Gerber RL, Bosch WR, Perez CA, et al. Quantitative dosimetric verification of an IMRT planning and delivery system. Radiother Oncol 1998;49:30516. 
10.  Li H, Dong L, Zhang L, Yang JN, Gillin MT, Zhu XR. Toward a better understanding of the gamma index: Investigation of parameters with a surfacebased distance method. Med Phys 2011;38:673041. 
11.  Sarkar B, Manikandan A, Nandy M, Gossman MS, Sureka CS, Ray A, et al. A mathematical approach to beam matching. Br J Radiol 2013;86:20130238. 
12.  Prescribing, recording, and reporting photonbeam intensitymodulated radiation therapy (IMRT): ICRU Report 83. J ICRU 2010 Apr;10(1):NP. doi: 10.1093/jicru/ndq002. 
13.  Clasie BM, Sharp GC, Seco J, Flanz JB, Kooy HM. Numerical solutions of the γindex in two and three dimensions. Phys Med Biol 2012;57:698197. 
[Figure 1], [Figure 2], [Figure 3], [Figure 4]
