
ORIGINAL ARTICLE 



Year : 2006  Volume
: 31
 Issue : 4  Page : 262268 

Empirical formula for the prediction of off axis ratios and isodose curves for a treatment planning system
Surajit Pal^{1}, R Ravishankar^{2}, RP Sharma^{3}, G Muthukrishnan^{4}, Dilip Kr Ray^{5}, SN Roy^{6}, DK Srivastava^{1}
^{1} Variable Energy Cyclotron Centre, 1/AF Bidhan Nagar, Kolkata, India ^{2} Health Physics Unit, VIP Road, Kolkata, India ^{3} L  404 VIP Enclave, VIP Road, Kolkata, India ^{4} B12, Adhunika Cooperative, 164/3A Lake Gardens, Kolkata, India ^{5} Chittaranjan National Cancer Institute, Kolkata, India ^{6} Department of Physics, VisvaBharati, Santiniketan, India
Correspondence Address: Surajit Pal Variable Energy Cyclotron Centre, 1/AF Bidhan Nagar, Kolkata  700 064 India
Source of Support: None, Conflict of Interest: None  Check 
DOI: 10.4103/09716203.29196
Abstract   
A mathematical model has been developed for prediction of off axis ratio (OAR), using Wood  Saxon term used to represent nuclear potential. This method has been satisfactorily applied for predicting OAR in case of ^{60}Co γ rays and high energy
Xrays. Investigations are considered upto a depth of 25 cm in the case of 4MV LINAC for which measurements were carried out in our laboratory using indigenously developed Radiation Field Analyzer. For 60Co γrays as well as 6 and 18MV LINAC beams we could get offaxis profiles only upto 20 cm. The shift δ between measured and predicted OAR is within ±2 mm except for 20 cm depth near the falling edge of the penumbra, where it is 2.80 mm. Software has been developed in Visual Basic 6 on Windows platform to plot Isodose curves, which is based on the mathematical modeling of OAR and central axis percentage depth dose.
Keywords: Isodose curves, offaxis ratio, percentage depth dose, treatment planning
How to cite this article: Pal S, Ravishankar R, Sharma R P, Muthukrishnan G, Ray DK, Roy S N, Srivastava D K. Empirical formula for the prediction of off axis ratios and isodose curves for a treatment planning system. J Med Phys 2006;31:2628 
How to cite this URL: Pal S, Ravishankar R, Sharma R P, Muthukrishnan G, Ray DK, Roy S N, Srivastava D K. Empirical formula for the prediction of off axis ratios and isodose curves for a treatment planning system. J Med Phys [serial online] 2006 [cited 2019 May 20];31:2628. Available from: http://www.jmp.org.in/text.asp?2006/31/4/262/29196 
The use of empirical formula for computing doses along the central axis and off axis ratio (OAR) along the transverse axis is advantageous in many ways. It reduces the databank requirements, which leads to a decrease in number of measurements to be made. Since one can predict the dose distributions for any field size, rectangular field sizes present no problem. In the case of highenergy xray machines, the shape of OAR curves and hence Isodose curves will depend, to a large extent, on the design of the beam flatteners.[1] It is well recognized among the clinical medical physicists that published Isodose curves can not be used and that the dose calculation performed must be specific to the teletherapy machine used, therefore, empirical formulations become quite necessary. The empirical formulations, so developed, should be able to reproduce the physical situation with a reasonable accuracy. It should also be possible to adopt the formula for any treatment machine. By reasonable accuracy we mean the distance in millimeters measured along the transverse axis between the points where the measured and calculated doses are same. According to the work of Thomas[2] and Khan,[3] a maximum shift of 2 mm is taken as sufficiently accurate, which we have adopted.
Various methods have been developed for the prediction of off axis ratios.[2],[4],[5] Usually in all these formulations, two or three sets of formulae are used to cover the penumbra and umbra regions etc. Kornelson[4] suggested the use of FermiDirac distribution function[6] to represent the OAR in the case of moving field treatment (SAD technique). But when we applied the same for fixed field cases (SSD technique) it could not be fitted to the FermiDirac Distribution. We have, therefore, developed a method in our laboratory based on WoodSaxon term,[7] which is generally used to represent nuclear potential. Here, a single formula can predict the off axis distributions covering the whole region. The method is based on the formula, developed for ^{ 60}sub Co therapy beam by Ravishankar[8] and we extended the same for the case of 4MV LINAC and other highenergy machines, thus demonstrating the versatility of the method.
In order to arrive at a suitable treatment plan, we need both central axis percentage depth dose and off axis ratio. In the earlier publication,[9] the method to calculate central axis percentage depth dose (CAPDD) using buildup concept has been explained. We present here the method to calculate OAR and also the development of treatment plan.
Materials and Methods   
The OAR is the ratio of offaxis dose to the central axis dose at the same depth. The plot of OAR vs. transverse distance resembles the formulation of Wood  Saxon term, extensively used in nuclear physics.[7] The off axis ratio R can be expressed as
where x_{ 0} is the half width of the square field at depth d below the surface of the water medium and is given by
and x is the off axis distance, f is the SSD, n is the field width at the depth of maximum dose ( dm). We have tested Eqn. (1) only for square fields. One has to study the validity for rectangular fields. The factor w is expressed as follows:
where S is the source diameter and fc is the source to collimator distance. For a given depth and field size, w should be strictly constant. But due to experimental errors or constructional details in the beam flattener, w varies with x . As in the studies with central axis depth dose distributions,[9] it is always possible to find a value of w , which will fit the OAR to the accuracy mentioned earlier. This is the basis of the formalism.
From Eqn. (1) we find the values of w for different OAR's. From these values of w , a suitable weighted average value of w is chosen so that it fulfills our requirement for the reproduction of OAR's. The w value was weighted over the OAR values, w (weighted) =εw x R /ε R . The w values at x = x_{ 0} or x close to x_{ 0} are omitted as w ®¥. Using this value of w , R was calculated for field sizes 5x5, 10x10, 15x15 and 20x20 cm^{ 2} and for depths 5, 10, 15 and 20 cm. The required data for all the field sizes and depths were taken from the measurements made at 4MV Medical LINAC Jeevan Jyoti2.[10] [Figure  1] shows the comparison between measured and calculated values for the various field sizes at a depth of 10 cm. [Table  1] gives the shift (δ mm) along the off axis between the points and position where the measured dose and the calculated dose are the same. Measurements are carried out in a water phantom (Radiation Field Analyzer) using 0.125cc, PTW make, semiflex ionization chamber. Depth of maximum dose is 1 cm for 4 MV xray beam. It is clear from [Table  1] that d lies within 2 mm.
For 4MV Medical LINAC of RRMC, f = 100 cm, fc = 38 cm and S = 0.2 cm. [Table  2] gives the w values (calculated and fitted) as a function of depth for a field size 10x10 cm^{ 2}.
The k values are again fitted against field sizes by the following equation
where n is the field size and a and b are constants. For our case a = 1.127 and b = 1.0544.
Hence using Eqs. (3) and (4) we can estimate w value for any field size and depth to predict OAR. This enables us to generate isodose curves as explained in the subsection Isodose Curves.
6 0 and 18 MV 0 x 0 rays
0The method was applied in the OAR of 6 and 18MV Xrays produced from Siemens LINAC (Mevatron, KD2), installed at the Chittaranjan National Cancer Institute, Kolkata for various field sizes. Depths of maximum dose is 1.4 and 3.4 cm for 6 and 18 MV xray beam, respectively. For Siemens LINAC (Mevatron, KD2), f = 100 cm, fc = 39.2 cm and S = 0.2 cm. [Table  3] gives the value of parameters a and b for 4, 6 and 18 MV xrays as calculated by using Eqn. (4). [Figure  2][Figure  3][Figure  4] show the comparison between measured and calculated values of OAR for various field sizes for 6 and 18 MV xrays. [Table  2] gives the w values (calculated and fitted) as a function of depth for a field size 10x10 cm[2] for these energies.
^{60}Co source
We have also applied the method for ^{ 60} Co Teletherapy machine Picker C2000 unit. The necessary input is taken from Van de Geijn[11] for a field size of 10x10 cm^{ 2} for various depths. [Figure  5] compares the typical measured and the calculated OAR values for a depth of 5 cm and the agreement between the two is quite reasonable. In this case the values are fitted by the following relationship.
Where A is 1.4736 and k is 1.0866. The values of f and fc considered for the calculation are 50 and 27 cm, respectively.
[Table  4] gives the calculated and fitted values of w for various depths.
Isodose curves in water medium
From the literature it is observed that the Decrement Line Method[12] and FanLine Method,[13] are used to generate isodose curves. However, we have used another approach where a semiempirical analytical method is applied to generate isodose curves. The basic equations used are given below.
In general, total dose D is represented by[14]
where B is the dose buildup factor and D_{ 0} is the contributions from the primary radiation. B is given by:
s represents the ratio of scattered to incident primary radiation,
We have used the same approach to represent the Central Axis Percentage Depth Dose (CAPDD)[9] in the same form as follows:
Total dose D_{ 1} at a depth of d[1] cm beyond depth of maximum dose is given by:
Total dose D_{ 2} at depth d_{ 2} cm (depth of maximum dose) is given by:
where D 0_{ 1} and D0_{ 2} are primary dose contributions and s_{ 1} and s_{ 2} are scattered components.
.^{.}. Relative dose
This can be approximated to
since s_{1} and s_{2} are relatively small.
D_{ 0} corresponds to the dose value of 0x0 cm^{ 2} field and is obtained by extrapolation of measured/published values of percentage depth doses for different depths. D_{ 0} values at different depths are found to decrease exponentially.
s represents the ratio of scattered to incident primary radiation, r is the depth inside the water phantom expressed in terms of mean free paths; a_{ 1} and k_{ 1} are constants. Equations (11) and (12) are valid for all depths greater than the depth of maximum dose. Also in the energy range considered there is no appreciable difference between the kerma and the absorbed dose.
The off axis distance, corresponding to the required OAR from Eqn. (1) is given by:
Here R´ is the inverse of OAR, which is normalised w.r.t. the percentage depth dose ratio along the central axis at the required depth using Eqn. (11). Eqn. (13) is valid for single fixed field only. We generated (x,y) coordinates for 10x10 cm^{ 2} field size and various percentages (590%). [Figure  6] shows the Isodose curves for the 4MV xrays and [Figure  7] gives the same for 6 and 18 MV xrays.
The total dose distributions for multiple fields treatment plan is obtained by calculating doses at each grid point of a matrix. The total area considered has been divided into grids with Cartesian coordinate system. The dose at the grid point, with a grid spacing of 0.2 cm, was then calculated from Eqs. (1) and (11), which give OAR and CAPDD, respectively. This procedure was repeated for all the fields for a given treatment plan and the total dose at each grid point was obtained by summing up the doses and then it was normalized against the maximum dose in percentage. [Figure  8] represents a typical treatment plan with an arbitrary patient contour and target for four orthogonal fields of 10x10 cm^{ 2}, where shaded portion represents the target volume.
Contour corrections are also included in this planning system. This is obtained by finding out air gap/extra tissue thickness due to shape of the contour and then multiplying by a suitable correction factor to each grid point.[15],[16]
Computer program
A computer program is developed in Visual Basic 6 on Windows platform to process the data faster and plotting of OAR and isodose curves. The algorithm is as follows:
i) Reading of measured OAR, off axis distance, field size, depth of maximum dose, source to surface distance, depth of measurement, source to collimator distance and source diameter, if applicable.
ii) Calculation of w for each depth and field size.
iii) Least square flitting of w with depth to calculate the constant k for each field size.
iv) Least square fitting of k with field size in the Eqn. (4) to find the constants a and b .
v) Finally, the values of a and b are used to calculate OAR for any field size and depth.
vi) For ^{ 60} Co beam least square fitting of w with depth is used to calculate not only k but also A , another constant from Eqn. (5).
vii) Knowing the above constants OAR's are calculated. Similarly, Central Axis Depth Doses are calculated through the Eqn. (11). The product of these two quantities gives the percentage depth dose for each grid point.
viii)Dose for each grid point is plotted in the 'Picture Box' to generate and display isodose curves.
ix) Treatment planning for multiple fields is executed by summing up dose values at each grid point from each field and the same is normalized against the maximum dose.
Results and discussion   
The factor w should be constant for a particular depth and field size for all off axis distances. But in actual practice it has been observed that w changes with off axis distance. To avoid any complexity, weighted average of w is taken and this takes care of variations of both w and OAR. In our experiments w is calculated for three depths for 6 and 18 MV xrays and five depths for 4 MV xrays. In both the measurements good agreement between calculated and measured OAR is obtained.
Referring to [Figure  1][Figure  2][Figure  3][Figure  4][Figure  5], the agreement between the measured and calculated OAR in the central region is within ±2%. In the penumbra region, the shift in the particular OAR is within 2 mm. We used this criterion of Khan[3] as this satisfies all field sizes and depths.
The fitted and calculated values of w are shown in the [Table  2], it is clear that the variation between fitted and calculated value is quite significant. However, this is not affecting the final result of off axis ratios, adversely. [Figure  9] is the graphical representation of actual and fitted values of w with respect to field size for 18 MV xray beam at 10 cm depth. [Table  5] gives the percentage deviation in dose along central axis between the isodose calculated by our algorithm and the same from a commercial planning system. Comparison of 4 MV xrays could not be considered since commercial values for 4 MV are not available locally.
Conclusion   
This is a simple technique for the prediction of off axis ratios, which is finally used for the development of treatment planning system along with mathematical modeling of central axis percentage depth dose technique developed earlier.[9] The system can be used by the cancer centers of our country that do not have access to sophisticated treatment planning systems due to high cost and maintenance. The system was evaluated at the Chittaranjan National Cancer Institute, Kolkata and the same is found to be accurate and user friendly. The software in the present form does not include corrections for heterogeneity of cancer patients as well as for irregular fields.
References   
1.  Webstar EW, Tsien KC. Atlas of radiation dose distributions. IAEA, Vienna 1965;1:331. 
2.  Thomas RL, Haybittle JL. A widely applicable method for computing dose distributions for external megavoltage beam. Br J Radiol 1975;48:74954. [PUBMED] 
3.  Khan FM. The physics of radiation therapy, 202. Lippincott Williams and Wilkins: Philadelphia; 2003. 
4.  Kornelsen RO, Young ME. Empirical equations for the representation of depth dose data for computerized treatment planning. Br J Radiol 1975;48:73948. [PUBMED] 
5.  Paul RB, Jayachandran CA. Development of a teletherapy treatment planning software. J Med Phys 1996;21:539. 
6.  Kittel C. Introduction to Solid State Physics. John Wiley and Sons: New York; 1971. p. 2445, 713. 
7.  Woods RD, Saxon DS. Diffuse surface optical model for nucleonnuclei scattering. Phys Rev 1954;95:5778. 
8.  Ravishankar R, Sengupta M, Bandyopadhyay T, Muthukrishnan G, Ghose AM. Semiempirical depth dose distribution formula for 60Co grays. Radiat Phys Chem 1994;44:1924. 
9.  Pal S, Ravishankar R, Sharma RP, Muthukrishnan G, Roy SN, Srivastava DK, et al . Prediction of central axis percentage depth dose distributions using buildup factor concept for a treatment planning system. J Med Phys 2003;28:704. 
10.  Pal S, Sharma RP, Muthukrishnan G. Dosimetric studies with 4MV Medical LINAC  Jeevan Jyoti2 at Cancer Centre Calcutta. Radiat Protect Environ 2001;24:3547. 
11.  van DE Geijn J. The computation of two and three dimensional dose distributions in cobalt 60 teletherapy. Br J Radiol 1965;38:36977. [PUBMED] 
12.  Orchard PG. Decrement lines: A new presentation of data in cobalt60 beam dosimetry. Br J Radiol 1964;37:75663. [PUBMED] 
13.  Milan J, Bentley RE. The storage and manipulation of radiation dose data in small digital computer. Br J Radiol 1974;47:11521. [PUBMED] 
14.  Attix FA, Roesch WC. Radiation dosimetry. Vol  I, Chap. 3. Academic Press: New York and London; 1968. p. 93155. 
15.  Johns HE, Cunningham JR. The physics of radiology. Charles C Thomas: Springfield; 1983. p. 3858. 
16.  Giessen PH vd. A method of calculating the isodose shift in correcting for oblique incidence in radiotherapy. Br J Radiol 1973;46:97882. [PUBMED] 
Figures
[Figure  1], [Figure  2], [Figure  3], [Figure  4], [Figure  5], [Figure  6], [Figure  7], [Figure  8], [Figure  9] Tables
[Table  1], [Table  2], [Table  3], [Table  4], [Table  5]
