

ORIGINAL ARTICLE 



Year : 2014  Volume
: 39
 Issue : 2  Page : 8592 

Estimation of distance error by fuzzy set theory required for strength determination of HDR ^{192} Ir brachytherapy sources
Sudhir Kumar^{1}, D Datta^{2}, SD Sharma^{1}, G Chourasiya^{1}, D. A. R. Babu^{1}, DN Sharma^{3}
^{1} Radiological Physics and Advisory Division, Bhabha Atomic Research Centre, CTCRS, Anushaktinagar, Maharashtra, India ^{2} Health Physics Division, Bhabha Atomic Research Centre, Trombay, Mumbai, Maharashtra, India ^{3} Health Safety and Environment Group, Bhabha Atomic Research Centre, Trombay, Mumbai, Maharashtra, India
Date of Submission  10Sep2013 
Date of Decision  20Jan2014 
Date of Acceptance  26Mar2014 
Date of Web Publication  23Apr2014 
Correspondence Address: Sudhir Kumar Radiological Physics and Advisory Division, Bhabha Atomic Research Centre, CT and CRS Building, Anushaktinagar, Mumbai 400 094, Maharashtra India
Source of Support: None, Conflict of Interest: None  Check 
DOI: 10.4103/09716203.131281
Abstract   
Verification of the strength of high dose rate (HDR) ^{192} Ir brachytherapy sources on receipt from the vendor is an important component of institutional quality assurance program. Either reference airkerma rate (RAKR) or airkerma strength (AKS) is the recommended quantity to specify the strength of gammaemitting brachytherapy sources. The use of Farmertype cylindrical ionization chamber of sensitive volume 0.6 cm ^{3} is one of the recommended methods for measuring RAKR of HDR ^{192} Ir brachytherapy sources. While using the cylindrical chamber method, it is required to determine the positioning error of the ionization chamber with respect to the source which is called the distance error. An attempt has been made to apply the fuzzy set theory to estimate the subjective uncertainty associated with the distance error. A simplified approach of applying this fuzzy set theory has been proposed in the quantification of uncertainty associated with the distance error. In order to express the uncertainty in the framework of fuzzy sets, the uncertainty index was estimated and was found to be within 2.5%, which further indicates that the possibility of error in measuring such distance may be of this order. It is observed that the relative distance l _{i} estimated by analytical method and fuzzy set theoretic approach are consistent with each other. The crisp values of l _{i} estimated using analytical method lie within the bounds computed using fuzzy set theory. This indicates that l _{i} values estimated using analytical methods are within 2.5% uncertainty. This value of uncertainty in distance measurement should be incorporated in the uncertainty budget, while estimating the expanded uncertainty in HDR ^{192} Ir source strength measurement.
Keywords: Brachytherapy, farmertype ionization chamber, fuzzy set theory, HDR ^{192} Ir source
How to cite this article: Kumar S, Datta D, Sharma S D, Chourasiya G, Babu D, Sharma D N. Estimation of distance error by fuzzy set theory required for strength determination of HDR ^{192} Ir brachytherapy sources. J Med Phys 2014;39:8592 
How to cite this URL: Kumar S, Datta D, Sharma S D, Chourasiya G, Babu D, Sharma D N. Estimation of distance error by fuzzy set theory required for strength determination of HDR ^{192} Ir brachytherapy sources. J Med Phys [serial online] 2014 [cited 2020 May 25];39:8592. Available from: http://www.jmp.org.in/text.asp?2014/39/2/85/131281 
Introduction   
The use of high dose rate (HDR) remote afterloading brachytherapy units are rapidly increasing in many countries around the world. Verifying the strength of HDR ^{192} Ir brachytherapy sources on receipt from the vendor is an important component of institutional quality assurance program. ^{[1],[2]} The recommended quantity to specify the strength of gammaemitting brachytherapy sources is either reference airkerma rate (RAKR) or airkerma strength (AKS). RAKR is the AKR to air, in air, at a reference distance of 1 m, corrected for attenuation and scattering; and refers to the quantity determined along the transverse bisector of the source. AKS is the AKR in air at a given distance corrected for attenuation and scattering and that is multiplied by the square of the given distance. ^{[1],[2],[3],[4],[5],[6]} Calibration of the ^{192} Ir sources used in HDR remote afterloading brachytherapy units is carried out either by using a thimble ionization chamber (inair jig method) or by using a welltype ionization chamber. A Farmertype cylindrical ionization chamber of nominal sensitive volume of 0.6 cm ^{3} is frequently used for inair calibration of HDR ^{192} Ir brachytherapy sources in addition to a suitable welltype ionization chamber. ^{[7],[8],[9]} The European Society for Therapeutic Radiology and Oncology (ESTRO) also recommends the use of thimble ionization chamber for calibration of HDR ^{192} Ir brachytherapy sources at hospitals. ^{[10]} Although welltype ionization chambers are preferred over cylindrical chambers for calibration of HDR ^{192} Ir brachytherapy sources due to ease in its use and reproducibility of source positioning, the Farmertype cylindrical ionization chamber is also used for RAKR or AKS measurement of HDR ^{192} Ir brachytherapy sources. ^{[11],[12],[13],[14]} This is due to the fact that cylindrical ionization chambers are readily available in the hospitals and in case of nonavailability of a welltype ionization chamber, the use of cylindrical ionization chamber is an obvious choice. It has also been demonstrated by Stump et al., ^{[15]} that the RAKR measured by Farmer type and welltype ionization chambers for different HDR sources are comparable within 0.5%.
A 370 GBq (10 Ci) ^{192} Ir source provides an ionization current of only about 1 × 10^{−11} A in a 1.0 cm ^{3} ionization chamber at a distance of 20 cm. ^{[16]} It is true that very near to a brachytherapy source, the radiation intensity changes very rapidly due to inversesquare law. A 0.1 cm error in a 10 cm distance causes a 2% error in calibration. ^{[17]} Small errors in positioning the chamber can translate into large errors in the estimation of source strength. Increasing the separation between centers of the chamber and the source will improve the measurement accuracy. However, this will result in proportionate reduction in the current, leading to larger percentage contributions by leakage current and gammaray scattering from the room surroundings and poor reproducibility. Getting closer of course worsens the distance error and requires a large geometric correction ^{[18],[19],[20]} for the size and shape of the ionization chamber.
The 7 distance method is recommended as a standard method to maximize the accuracy in measurement of the strength of HDR brachytherapy sources by using cylindrical ionization chamber. ^{[8],[9],[10],[15],[16]} While using the cylindrical chamber method, it is required to determine the positioning error of the ionization chamber with respect to the source, which is commonly called as the distance error. Earlier, we have developed the analytical methods to estimate the distance error required to determine the source strength using 7 distance method by cylindrical ionization chamber. ^{[21],[22]} As further research in this work, an attempt has been made to apply the fuzzy set theory to estimate the subjective uncertainty associated with the distance error, which is the subject matter of this paper.
Fuzzy set theory has been applied for risk analysis ^{[23]} and image analysis ^{[24]} in the domain of medical dosimetry. In view of these applications, we have proposed an approach of applying this fuzzy set theory in the quantification of uncertainty associated with the distance error required for the measurement of strength of HDR ^{192} Ir brachytherapy sources. While using, Farmertype cylindrical ionization chamber to measure the strength of a brachytherapy source, distance has to be measured accurately. Since the distance measured possesses some error during the measurement and the input components are imprecise, fuzzy set theory is an appropriate tool to quantify the uncertainty due to such ambiguity present in input component. ^{[25]} The fuzzy set theory is strictly applicable where there is insufficient information in the measured data.
Materials and Methods   
Multiple distance measurement technique
The microSelectronHDR unit from Nucletron was used in this work. This unit uses an old design micro Selectron ^{192} Ir HDR brachytherapy source (Nucletron B V, Veenendaal, Netherlands) with 370 GBq (10 Ci) nominal activity to treat brachytherapy patients with HDR comparable to teletherapy. The old microSelectron HDR source is cylindrical in geometry with 0.6 mm active diameter and 3.5 mm active length. PTW 30001 0.6 cm ^{3} Farmertype ionization chamber (PTW, Freiburg, Germany) was used in this work. Further details about the old microSelectron HDR source and the PTW 30001 ionization chamber are available elsewhere. ^{[26],[27]}
To determine experimentally, the RAKR of an HDR ^{192} Ir brachytherapy source using a Farmertype cylindrical ionization chamber, a multiple distance measurement technique was used. This measurement has historically been made at seven separate distances. Thus the technique has been termed the '7 distance' measurement (7DM). While using cylindrical ionization chamber (0.6 cm ^{3} ) for measurement of the strength of HDR ^{192} Ir brachytherapy sources, it is necessary to estimate three items, viz. (i) The positioning error of the ionization chamber with respect to the source which is commonly called 'distance error (±)', (ii) the contribution of scatter radiation (M _{s} ) from the floor, walls ceiling, and other material in the treatment room, and (iii) a proportionality constant. The 7DM was suggested to determine these parameters and thereafter the strength of HDR ^{192} Ir brachytherapy sources. ^{[8],[9],[15],[16]} Although Kumar et al., ^{[21],[22]} has described in detail the procedure for measuring the RAKR of HDR ^{192} Ir brachytherapy sources, a brief review of the method is given below for the sake of completeness in description.
In 7DM method, the output of the source in air is measured at seven different distances each corresponding to a meter reading M _{d}, which is the sum of primary and scattered radiation:
where d _{i} is the apparent distance between source and the chamber centers, c is the offset error in the measurement of the distance and commonly known as 'distance error'; f is a proportionality constant which is independent of distance. On solving Equation 2, one may obtain the following functional form for relative distance l _{i} between the successive measurement points (i = 1, 2,…, 6) ^{[22]}
Fuzzy set theory
Zadeh ^{[25],[29]} introduced the fuzzy set as a class of object with a continuum of grades of membership. In contrast to classical crisp sets where a set is defined by either membership or nonmembership, the fuzzy approach relates to a grade of membership between [0,1], defined in terms of the membership function of a fuzzy number. Hence, the classical notion of binary membership has been modified for the representation of uncertainty in data. The details about fuzzy set may be found elsewhere, ^{[29]} however, for the sake of completeness, a brief description of the definition of a fuzzy set and its fundamental properties pertaining to the topic of the present work is described here. This is a paradigm shift in which the crisp variable is fuzzified through a membership function or a linguistic variable depending upon the specific problem. Strictly speaking, alphacut theory of the fuzzy set ^{[29]} is adopted to compute the uncertainty associated with the distance.
Basic concept
A fuzzy set A is denoted by an ordered set of pairs (x, μ (x)), where, the element x μ X (crisp value) of a specific universe and μ (x) denotes the degree of membership, μ (x) μ [0,1]. A membership function can be of any shape depending on the type of a fuzzy set it belongs to. The only condition a membership function must satisfy is it should vary between 0 and 1. The membership function of a fuzzy set, A is defined in the form of a triangular or trapezoidal fuzzy number as shown in [Figure 1]a and b. The analytical form of the triangular membership function is depicted in Equation 13
Basically, alphacut is an interval and in practice, interval arithmetic operation ^{[30]} is carried out for obtaining the membership value of the output of a model containing the fuzzy input. The present paper, applies the alphacut value of the fuzzy set.  Figure 1: Pictorial representation of the membership function of a fuzzy set (a) triangular, (b) trapezoidal, and (c) support of a triangular fuzzy number
Click here to view 
Implementation of alphacut of a fuzzy set
In the present problem, fuzzy set theory has been applied to estimate the relative distance, l _{i} as shown in Equation 3. Each parameter in Equation 3 was treated as triangular fuzzy number because, the experimental determination provides the most likely value with the two extreme bounds scattered by the error obtained during measurement. The computation scheme of the membership value of the relative distance l _{i} is as follows:
The parameters, f, M _{di}, M _{0} , and M _{s} are taken into account as triangular fuzzy number. Alphacut representation ranged from 0 to 1 with an increment of 0.1 of these parameters was applied in Equation 3. In order to derive the membership value of the relative distance, l _{i} using the alphacut representation of fuzzy parameters as mentioned above, we used the interval arithmetic operations and the algorithm for computation. The details of this computation are given in the Appendix.
Support and uncertainty of a fuzzy set
A fuzzy set having triangular membership function is always characterized by its support and height. If height of a fuzzy set is 1 and if that fuzzy set is bounded by two extremes, then it is called a triangular fuzzy number. Support of such a triangular fuzzy number is defined as the range of the extremes at alphacut = 0 as shown in [Figure 1]c. From [Figure 1]c, we can write the support of a fuzzy set as S = (R−P), where, R and P are the two extreme bounds. In order to express the uncertainty in the framework of fuzzy sets, we define uncertainty index ^{[30]} as the ratio of the support to the most likely value (crisp value at membership equal to 1). Again, from [Figure 1]c, uncertainty index of the given fuzzy set is written as U = (S/Q), where, Q is the most likely value. The uncertainty index for each relative distance measured experimentally was estimated.
Results and Discussion   
We have estimated the uncertainty of the relative distance, l _{i} [= (d _{i} d _{0} )] for experimentally measured distances such as 5, 10, 15, 20, 25, and 30 cm and the corresponding membership functions are shown in [Figure 2]af. It can be interpreted from [Figure 2]af that the membership function μ (l _{i}) of the distance (l _{i} ) for each measurement distance is turned out to be a triangular in shape because the initial consideration of the subjectivebased uncertain parameters are taken into consideration as "around the measured value". On the contrary, had this consideration been within the phrase of "approximately lying between two different distances", we would have obtained the shape of the membership function of the output as trapezoidal. Since the measurement uncertainty is always quoted at one sigma level, fuzzy set theorybased approach of uncertainty quantification is also quoted at an equivalent level, and here this is 0.5 alphacut value of the fuzzy set. That is why, in this work, 0.5 alphacut of the output fuzzy set is considered to express the subjective uncertainty. Results of alphacut = 0.5 of the fuzzy distance (l _{i} ) along with experimental and analytical values of (l _{i}) are given in [Table 1] and it can be seen that analytical value as well as the experimentally measured relative distance lie within the bounds of the subjective uncertainty of the l _{i}. Support of each triangular membership values corresponding to each experimental distance and the associated uncertainty index are further shown in [Table 2]. It can be seen from [Table 2] that the uncertainty indices remain same for all the experimentallymeasured distances indicating that each and every triangular fuzzy membership function is normalized and convex. Maximum value of the uncertainty index is found to be within 2.5%, which further indicates that the possibility of error in measuring such distance may be of this order. It is worth mentioning here that in absence of a quantified value of uncertainty associated with distance measurement, the overall uncertainty in source strength measurement has been estimated earlier by ignoring the uncertainty associated with distance error. ^{[15]} This work reveals that the uncertainty associated with distance error should not be ignored while preparing the uncertainty budget in the determination of HDR source strength.  Figure 2: Pictorial representation of membership value μ (li) and relative distance li (cm) [= (di − d0)] for the distance (a) 5, (b) 10, (c) 15, (d) 20, (e) 25, and (f) 30 cm. Here, ◊ is lower bound and ƒ¢ is upper bound
Click here to view 
 Table 1: Comparison of experimentally recorded, analytically calculated, and fuzzy set theory computed values of li [= ( did0)]
Click here to view 
Conclusion   
Uncertainty of the positioning error, the so called "distance error", of the ionization chamber with respect to the source was evaluated. Fuzzy set theory was applied for this evaluation due to the subjectivity involved in the experimental facility. Uncertainty in the possible input parameters was addressed as triangular fuzzy number. Propagation of uncertainty of the input parameters is carried out on the basis of the model described in this work (see subsection: Algorithm to compute alphacut distance representation of distance) via the alphacut of a fuzzy set. The crisp values of l _{i} estimated using analytical method lie within the bounds computed using fuzzy set theory. This indicates that l _{i} values estimated using analytical methods are within 2.5% uncertainty. This value of uncertainty in distance measurement should be incorporated in the uncertainty budget, while estimating the expanded uncertainty in HDR ^{192} Ir source strength measurement.
Acknowledgment   
The authors wish to express their gratitude to Shri HS Kushwaha, Distinguished Scientist and EX. Director of Health, Safety and Environment Group, BARC for his valuable suggestions and technical discussions during preparation of this manuscript.
Appendix   
Algorithm to compute alphacut representation of distance
References   
1.  Nath R, Anderson LL, Meli JA, Olch AJ, Stitt JA, Williamson JF. Code of practice for brachytherapy physics: Report of the AAPM Radiation Therapy Committee Task Group No. 56. American Association of Physicists in Medicine. Med Phys 1997;24:155798. 
2.  Kubo HD, Glasgow GP, Pethel TD, Thomadsen BR, Williamson JF. High doserate brachytherapy treatment delivery: Report of the AAPM Radiation Therapy Committee Task Group No. 59. Med Phys 1998;25:375403. 
3.  International Commission on Radiation Unit and Measurements. Dose and volume specification for reporting interstitial therapy, Report No. 58. ICRU, Bethesda, USA; 1997. 
4.  International Commission on Radiation Unit and Measurements. Dose and volume specification for reporting intracavitary therapy in gynecology, Report No. 38. Washington, DC; 1985. 
5.  Specification of brachytherapy sources. Memorandum from the British Committee on Radiation Units and Measurements. Br J Radiol 1984;57:9412. [PUBMED] 
6.  International Commission on Radiation Unit and Measurements. Radiation quantities and units, Report No. 33. ICRU, Washington DC; 1980. 
7.  Marechal MH, de Almeida CE, Ferreira IH, Sibata CH. Experimental derivation of wall correction factors for ionization chambers used in high dose rate 192Ir source calibration. Med Phys 2002;29:15. 
8.  International Atomic Energy Agency. Calibration of photon and beta ray sources used in Brachytherapy, TECDOC 1274. IAEA, Vienna; 2002. 
9.  International Atomic Energy Agency. Calibration of brachytherapy sources: Guidelines of standardized procedures for the calibration of brachytherapy sources at secondary standard dosimetry laboratories and hospitals, TECDOC 1079. IAEA, Vienna; 1999. 
10.  European Society for Therapeutic Radiology and Oncology. A practical guide to quality control of brachytherapy equipment, ESTRO booklet no. 8; 2004. 
11.  Douysset G, Gouriou J, Delaunay F, DeWerd L, Stump K, Micka J, et al. Laboratoire National Henri Becquerel, University of Wisconsin Accredited Dosimetry Calibration Laboratory. Comparison of dosimetric standards of USA and France for HDR brachytherapy. Phys Med Biol 2005;50:196178. 
12.  Patel NP, Majumdar B, Vijayan V. Study of scattered radiation for inair calibration by a multipledistance method using ionization chambers and an HDR 192Ir brachytherapy source. Br J Radiol 2006;79:34752. 
13.  Douysset G, Sander T, Gouriou J, Nutbrown R. Comparison of airkerma standards of LNELNHB and NPL for 192Ir HDR brachytherapy sources: EUROMET project no 814. Phys Med Biol 2008;53:N8597. 
14.  Soares CG, Douysset G, Mitch MG. Primary standards and dosimetry protocols for brachytherapy sources. Metrologia 2009;46:S8098. 
15.  Stump KE, Dewerd LA, Micka JA, Anderson DR. Calibration of new high dose rate 192Ir sources. Med Phys 2002;29:14838. 
16.  Goetsch SJ, Attix FH, Pearson DW, Thomadsen BR. Calibration of 192Ir highdoserate afterloading system. Med Phys 1991;18:4627. 
17.  Goetsch SJ, Attix FH, Dewerd LA, Thomadsen BR. A new reentrant ionization chamber for the calibration of iridium192 high dose rate sources. Int J Radiat Oncol Bio Phys 1992;24:16770. 
18.  Bielajew AF. Correction factors for thickwalled ionization chambers in pointsource photon beams. Phys Med Biol 1990;35:50116. 
19.  Bielajew AF. An analytical theory of the pointsource nonuniformity correction factor for thickwalled ionization chambers in photon beams. Phys Med Biol 1990;35:51738. 
20.  Kondo VS, Randolph ML. Effect of finite size of ionization chambers on measurement of small photon sources. Radiat Res 1960;13:3760. 
21.  Kumar S, Srinivasan P, Sharma SD, Mayya YS. A simplified analytical approach to estimate the parameters required for strength determination of HDR 192Ir brachytherapy sources using Farmertype ionization chamber. Appl Radiat Isot 2012;70:2829. 
22.  Kumar S, Srinivasan P, Sharma SD, Subbaiah KV, Mayya YS. Evaluation of scatter contribution and distance error by Iterative methods for strength determination of HDR 192Ir brachytherapy source. Med Dosim 2010;35:2307. 
23.  Castiglia F, Giardina M, Tomarchio E. Risk analysis using fuzzy set theory of the accidental exposure of medical staff during brachytherapy procedures. J Radiol Prot 2010;30:4962. 
24.  Park SB, Monroe JI, Yao M, Machtay M, Sohn JW. Composite radiation dose representation using fuzzy set theory. Inf Sci 2012;187:2042015. 
25.  Ross TJ. Fuzzy logic with engineering applications. McGrawHill Inc.; 1997. 
26.  Williamson JF, Li Z. Monte Carlo aided dosimetry of the microselectron pulsed and high doserate 192Ir sources. Med Phys 1995;22:80919. 
27.  Pena J, SanchezDoblado JF, Capote R, Terron JA, Gomez F. Monte Carlo correction factors for a Farmer 0.6 cm ^{3} ion chamber dose measurement in the buildup region of the 6 MV clinical beam. Phys Med Biol 2006;51:152332. 
28.  Williamson JF, Nath R. Clinical implementation of AAPM Task Group 32 recommendations on brachytherapy source strength specification. Med Phys 1991;18:43948. 
29.  Dubois D, Prade H. editors. Fuzzy sets and systems: Theory and applications. New York: Academic Press; 1980. 
30.  Klir GJ, Yuan B. Fuzzy sets and fuzzy logic: Theory and applications, prenticehall, upper saddle river. New Jersey; 1994. 
[Figure 1], [Figure 2]
[Table 1], [Table 2]
