

TECHNICAL NOTE 



Year : 2012  Volume
: 37
 Issue : 2  Page : 107111 

A simple calculation method for determination of equivalent square field
Seyed Ali Shafiei^{1}, Hadi Hasanzadeh^{2}, Seyed Ahmad Shafiei^{3}
^{1} Department of BiophysicsBiochemistry, School of Medicine, Rafsanjan University of Medical Sciences, Rafsanjan, Iran ^{2} Department of Medical Physics, Semnan University of Medical Sciences, Semnan, Iran ^{3} Department of Biology, University of Isfahan, Isfahan, Iran
Date of Submission  14Nov2011 
Date of Decision  16Feb2012 
Date of Acceptance  16Feb2012 
Date of Web Publication  7Apr2012 
Correspondence Address: Seyed Ali Shafiei Department of BiophysicsBiochemistry, School of Medicine, Rafsanjan University of Medical Sciences, Rafsanjan Iran
Source of Support: None, Conflict of Interest: None  Check 
DOI: 10.4103/09716203.94746
Abstract   
Determination of the equivalent square fields for rectangular and shielded fields is of great importance in radiotherapy centers and treatment planning software. This is accomplished using standard tables and empirical formulas. The goal of this paper is to present a formula based on analysis of scatter reduction due to inverse square law to obtain equivalent field. Tables are published by different agencies such as ICRU (International Commission on Radiation Units and measurements), which are based on experimental data; but there exist mathematical formulas that yield the equivalent square field of an irregular rectangular field which are used extensively in computation techniques for dose determination. These processes lead to some complicated and timeconsuming formulas for which the current study was designed. In this work, considering the portion of scattered radiation in absorbed dose at a point of measurement, a numerical formula was obtained based on which a simple formula was developed to calculate equivalent square field. Using polar coordinate and inverse square law will lead to a simple formula for calculation of equivalent field. The presented method is an analytical approach based on which one can estimate the equivalent square field of a rectangular field and may be used for a shielded field or an offaxis point. Besides, one can calculate equivalent field of rectangular field with the concept of decreased scatter radiation with inverse square law with a good approximation. This method may be useful in computing Percentage Depth Dose and TissuePhantom Ratio which are extensively used in treatment planning.
Keywords: Equivalent square field, irregular field, Percentage Depth Dose, scattered beam, TissuePhantom Ratio
How to cite this article: Shafiei SA, Hasanzadeh H, Shafiei SA. A simple calculation method for determination of equivalent square field. J Med Phys 2012;37:10711 
How to cite this URL: Shafiei SA, Hasanzadeh H, Shafiei SA. A simple calculation method for determination of equivalent square field. J Med Phys [serial online] 2012 [cited 2020 Apr 8];37:10711. Available from: http://www.jmp.org.in/text.asp?2012/37/2/107/94746 
Introduction   
Determination of equivalent square field of a rectangular or shielded field (in order to calculate PDD, TAR) is one of the prevalent methods of dose determination on the central axis of radiation beam. There exist different approaches to determine equivalent field, such as empirical table published by ICRU, ^{[1]} in which nonpresented values are obtained by interpolation or computational methods which perform via a simple empirical relation as and some other methods which are not clearly based on the concept of equivalent field. ^{[2]} The concept of equivalent square field is usually utilized to calculate percentage depth dose and similar quantities in radiotherapy, which in most dose calculation softwares and algorithms is used for irregular fields too. The abovementioned formula is used extensively in optimizing these tables by considering scattered radiation and beam quality. ^{[3],[4],[5]} Many attempts have been made to calculate depth dose independent of equivalent field, and some algorithms (e.g. ALFARD) have been developed ^{[4]} in which equivalent field is one of its Kernel's major factors ^{[6],[7],[8]} and is used to calculate absorbed dose at an offaxis point. ^{[9].} Various methods of dose calculating have been investigated by Tudor et al. in ellipsoidal fields in which it was concluded that dose calculation based on the concept of equivalent field results in a lower percentage error. ^{[10]}
The aim of the present paper is to prove a formula based on the reduction in scattered beam due to inverse square law at a point on the central axis of beam, based on which one can extract equivalent square field of a rectangular field.
Materials and Methods   
Considering a rectangular field with dimensions of a × b, absorbed dose at point A is due to primary and scattered beams. So, one can consider that absorbed dose at point A is caused by square field with dimensions of a × a and two strips at the corners with dimensions of exist at both side of center [Figure 1].
On the other hand, the a × b rectangular field is equivalent to x × x square field. So, if a × a field is omitted from x × x field (necessarily x × x > a × a) [Figure 2], scattered portion from the rest of the field is equivalent to scattered radiation resulted from the two abovementioned fields. In other words,
 Figure 2: Geometry of x × x field in which the central part is a square field with dimensions of a × a
Click here to view 
where S is scattered beam, x ^{2}  a ^{2} is stripped peripheral surface in [Figure 2], and a(b  a) is total stripped surface in [Figure 1].
Using Cartesian system to compute scattered radiation at point A from stripped surface [Figure 3], one can obtain an expression as follows:
 Figure 3: (a) An arbitrary rectangular field with Cartesian coordinate system parameters, (b) modified field to obtain equivalent square field and polar coordinate system parameters
Click here to view 
Solving the above integral using Mathematica math software (Wolfram Research Inc., IL, USA) redounds to the following solution:
Although the above answer seems to be complicated, since all variables in Equation (3) are numerical, the answer is easily obtained by computer for each desired value of a and b.
Omitting midmost area in [Figure 1] prevents the singularity might happen due to above integration, when r reaches to zero. The main problem arises when the effect of scattered beams in [Figure 2] is to be calculated; if the above procedure follows in Cartesian coordinate system, a nonlinear equation similar to Equation (3) obtains with respect to x in which it is not possible to eliminate x from it. This problem is resolved using polar coordinates system instead of Cartesian coordinate system. To do so, the central hollow square is divided into four trapeziums, as shown in [Figure 3], and the effect of scattered beams at point A is calculated as the following equation:
Comparing the above simple logarithmic equation with Equation (3), an expression for the scattered radiation obtains as follows:
In the above equation, x is the dimension of the equivalent square field and can be easily calculated as:
Results   
A comparison of published tables by ICRU [Table 1] with tables obtained from Equation (7) [Table 2] and from [Table 3] ^{[11]} shows that in large fields, which are usually used in radiotherapy, calculated values in [Table 1] are very close to the values in [Table 2]. For instance, if we consider a 20 × 10 field, a value of 13 is obtained from [Table 1] and values of 13.1 and 13.3 are obtained from [Table 2] and [Table 3], respectively.  Table 1: Equivalent values for fields with dimensions from 2 cm to 30 cm published by ICRU
Click here to view 
Discussion   
In large fields, values obtained from suggested method [Table 2] are much closer to the values suggested by ICRU [Table 1] than the values computed by conventional methods [Table 3], the accuracy of which has been accepted by some treatment planning software and texts. ^{[11].} Besides, many dose correction methods in their analytical part use conventional method mentioned above to calculate equivalent square field through the correction of scattered beams. However, it is clear that in small enough fields, the error will be considerably high in suggested and conventional methods [Table 2] and [Table 3]; this error has been removed by applying correction coefficients to ^{[12]}
Recommendations for researchers based on the capabilities of the above formula
The power of the proposed method with respect to the previous methods is its possibility to calculate the equivalent field for shielded fields. Most computational methods and formulas used in treatment planning for determination of equivalent square field and dose calculation provide a fixed equivalent field for a shielded field without considering the shield position in the field; this is while the scattered beams reaching point A depend on shield position and consequently the equivalent fields must be different [Figure 4]. As it is obvious from this figure, although the shield position and so the scatter is different, the equivalent square field for both conditions is the same.  Figure 4: Two fields with the same dimension and similar shields but different shield positions in the fields
Click here to view 
This will be important when working with fields and shields of large dimensions. However, this will be easily achieved by the proposed method for obtaining Equation (7) by separating a central square field to calculate equivalent field and scattered beams considering inverse square law. Besides, in multileaf collimators [Figure 5], if the dimensions of each leaf are known, this technique may be easily used to calculate the dose and equivalent square field (Equation (8)).  Figure 5: (a) The field acquired by multileaf collimator and (b) the equivalent field for an offaxis point
Click here to view 
Where a _{i} and b _{i} are dimensions of area outside the a × a square, which totally results in the desired field shape, this is a potential capability which needs more detailed work to be practically useful. The importance of the proposed method is that it enables one to calculate dose at an arbitrary offaxis point provided this point does not place on the sides of the rectangular field.
It is possible to consider the effects of beam quality by calculation of attenuation in phantom and inserting the e ^{μr} factor in Equations (2 and 4); considering the dependence of attenuation coefficient on energy and atomic number, the proposed formula might become more appropriate with results much closer to real conditions, which needs a brief simulation and phantom measurement.
Conclusion   
The presented method is an analytical approach based on which one can estimate the equivalent square field of a rectangular field and may be used for a shielded field or an offaxis point. Besides, one can calculate equivalent field of rectangular field with the concept of decreased scatter radiation with inverse square law with a good approximation. This method may be useful in computing PDD and TPR which are extensively used in treatment planning.
References   
1.  Depth dose tables for use in radiotherapy. A survey, prepared by the Scientific SubCommittee of the Hospital Physicist's Association, of central axis depthdose data measured in water or equivalent media. Br J Radiol 1978;10 Suppl 11:196. 
2.  McCurdy BM, Pistorius S. Determination of equivalent photon fields through integrated 1D convolution kernels. Phys Med Biol 1999;44:297185. [PUBMED] [FULLTEXT] 
3.  Pirotta M, Aquilina D, Bhikha T, Georg D. Comparison of the ESTRO Formalism for monitor unit calculation with a Clarkson based algorithm of a treatment planning system and a traditional "fullscatter" methodology. Zeitschrift fur Medizinische Physik 2005;15:23546. [PUBMED] 
4.  Venselaar JL, Heukelom S, Jager HN, Mijnheer BJ, van Gasteren JJ, van Kleffens HJ,et al. Is there a need for a revised table of equivalent square fields for the determination of phantom scatter correction factors? Phys Med Biol 1997;42:236981. [PUBMED] [FULLTEXT] 
5.  Xiao Y, Bjärngard BE, Reiff J. Equivalent fields and scatter integration for photon fields. Phys Med Biol 1999;44:105365. 
6.  Lombardi P, Fiorino C, Cattaneo GM, Calandrino R. Monitor unit calculation in 6 MV irregularly shaped beams  Accuracy in clinical practice. Br J Radiol 1997;70:63844. [PUBMED] [FULLTEXT] 
7.  Storchi P, Woudstra E. Calculation of the absorbed dose distribution due to irregularly shaped photon beams using pencil beam kernels derived from basic beam data. Phys Med Biol 1996;41:63756. [PUBMED] [FULLTEXT] 
8.  Stern RL, Heaton R, Fraser MW, Goddu SM, Kirby TH, Lam KL . Verification of monitor unit calculations for nonIMRT clinical radiotherapy: Report of AAPM Task Group 114. Med Phys 2011;38:50430. 
9.  Storchi P, Woudstra E. Calculation models for determining the absorbed dose in water phantoms in offaxis planes of rectangular fields of open and wedged photon beams. Phys Med Biol 1995;40:51127. [PUBMED] [FULLTEXT] 
10.  Tudor GS, Thomas SJ. Equivalent diameters of elliptical fields. Br J Radiol 2004;77:9413. [PUBMED] [FULLTEXT] 
11.  Khan FM. The physics of radiation therapy. Philadelphia, PA: Lippincott Williams and Wilkins; 2010. 
12.  Thomas SJ, Eaton DJ, Tudor GS, Twyman NI. Equivalent squares for small field dosimetry. Br J Radiol 2008;81:897901. [PUBMED] [FULLTEXT] 
[Figure 1], [Figure 2], [Figure 3], [Figure 4], [Figure 5]
[Table 1], [Table 2], [Table 3]
