|Year : 2018 | Volume
| Issue : 3 | Page : 107-112
Dosimetric effect of the gantry rotations of a novel trunk phantom using an area integration algorithm
Michael Onoriode Akpochafor1, Akintayo Daniel Omojola2, Muhammad Yaqub Habeebu1, Samuel Olaolu Adeneye1, CB Madu3, ME Ekpo4, Moses Adebayo Aweda1, Temitope Aminat Orotoye1
1 Department of Radiation Biology, Radiotherapy, Radiodiagnosis and Radiography, College of Medicine, Lagos University Teaching Hospital, Lagos, Nigeria
2 Department of Radiology, Federal Medical Center, Medical Physics Unit, Asaba, Nigeria
3 Department of Radiology, Medical Physics Unit, University College Hospital, Ibadan, Nigeria
4 Department of Physics, Faculty of Science, University of Ibadan, Ibadan, Nigeria
|Date of Submission||15-Sep-2017|
|Date of Decision||11-Nov-2017|
|Date of Acceptance||12-Jan-2018|
|Date of Web Publication||1-Jun-2018|
Mr. Akintayo Daniel Omojola
Department of Radiology, Federal Medical Center, Medical Physics Unit, Asaba, Delta State
Source of Support: None, Conflict of Interest: None
Background: Treatment planning systems (TPSs) have proved to be a useful tool in predetermining how a treatment outcome will be in radiotherapy. The accuracy of any TPS to calculate dose to any arbitrary point within a material is largely dependent on the mathematical algorithm used. Aims: The purpose of this study was to design a local trunk phantom and use the phantom to check the percentage dose accuracy of the Area Integration Algorithm of a Precise PLAN 2.16 TPS if it is in agreement with results obtained from manufacturer's verification by varying the gantry angle and whether it is within ± 5% International Commission on Radiation Units and Measurements (ICRU) minimal limit. Materials and Methods: The study was executed with a locally designed phantom made of Plexiglas with six insert and an ionization chamber port. The phantom was simulated using a HiSpeed NX/i computed tomography scanner and Precise PLAN 2.16 TPS for application of beam setup parameters. The mimicked organs for the inserts were: 25%–75% Glycerol-Water for liver, pure carboxyl methyl cellulose was used for lungs, 30%–70% Glycerol-Water for muscle, 40%–60% Glycerol-Water was used for adipose, pure Sodium hypochlorite was used for bone and pure sodium laureth sulfate (Texapon) for kidney. Results: The maximum percentage (%) deviation with a large field for six inhomogeneous inserts and with bone only homogeneous inserts were 3.4% and 2.9%, respectively. The maximum % deviation with a small field for six inhomogeneous inserts was 3.2%. The % deviation between the solid water phantom and the locally designed phantom was 3.5%. Conclusion: The Area Integration Algorithm has shown an overall accuracy of 4% below 5% ICRU minimal limit. There was no statistically significant difference in field sizes and in inhomogeneity/homogeneity, respectively. Variation exists in % deviation for small field size with parallel opposed field between our verification and the manufacturers.
Keywords: Phantom, treatment planning system, area integration algorithm, irregular field program, computed tomography
|How to cite this article:|
Akpochafor MO, Omojola AD, Habeebu MY, Adeneye SO, Madu C B, Ekpo M E, Aweda MA, Orotoye TA. Dosimetric effect of the gantry rotations of a novel trunk phantom using an area integration algorithm. J Med Sci 2018;38:107-12
|How to cite this URL:|
Akpochafor MO, Omojola AD, Habeebu MY, Adeneye SO, Madu C B, Ekpo M E, Aweda MA, Orotoye TA. Dosimetric effect of the gantry rotations of a novel trunk phantom using an area integration algorithm. J Med Sci [serial online] 2018 [cited 2018 Aug 18];38:107-12. Available from: http://www.jmedscindmc.com/text.asp?2018/38/3/107/229257
| Introduction|| |
The International Commission on Radiation Units and Measurements (ICRU) has recommended that radiation dose must be delivered to within ± 5% of the prescribed dose Although ICRU report 24 also states that recommended uncertainty in the delivered dose to a phantom (mimicked human subject) at “optimal model” should be ± 2.5% and at “minimal” or “lowest acceptable” model be ± 5% and International Atomic Energy Agency (IAEA) dosimetry limit be ± 3.0%.,, When commissioning treatment planning dose calculation algorithms in radiotherapy, the aim is often to achieve good agreement between calculated and measured doses within 1%–2% for open and wedge (block or compensators) fields in water. While this is possible to achieve using both measurement- and model-based algorithms in water phantoms, such an agreement is usually not possible for measurement-based algorithms in phantoms with heterogeneities. This can be explained based on the fact that measurement-based models are able to account for the effect of tissue inhomogeneities on the primary radiation. However, correcting for the scatter radiation is difficult since it depends on field size, beam energy and shape, location, and density of the inhomogeneities. In contrast, model-based algorithms can account for the effect of tissue inhomogeneities on the scatter radiation using the density scaling method or other approaches.,,,,
Several technique of carrying out the quality assurance of treatment planning system (TPS) has been proposed by the various authors.,,,, Likewise, the reduction of errors and uncertainties in the dose calculation plays an important role in the success of a treatment procedure. The performance and quality of any TPS is dependent on the type of algorithm used.,,,,,,,
Treatment planning requires the ability to calculate dose to any arbitrary point, within the patient, for any beam orientation. The Area Integration Algorithm (Irregular Field program) is well suited for this type of general problem. Patient tissue inhomogeneities, beam blocks, and beam compensators are included in the calculation model. The Area Integration Algorithm (Irregular Field program) requires the separation of the dose into primary and scatter components. The primary component is usually computed, and it includes the transmission through any blocks and blocking trays, beam compensators, and patient inhomogeneities. The scatter component is usually computed, and it includes the presence of blocks, beam compensators, and curvature of the patient but not patient inhomogeneities.,,
The concept of this dosimetry of Area Integration Algorithm (Irregular Field program) involves the use of tissue-maximum ratio and scatter-maximum ratio which is analogous to tissue-air ratio (TAR) and scatter-air ratio concepts. The underlining program equation of the Area Integration (Irregular Field program) which is similar to external beam program is as follows:
TRAY and TRAY2 = are the tray factors.
OUTPUT = the output factor normalized to a (0 × 0) field size at a distance source to axis distance + dose at maximum (SAD + Dmax).
FSC = the airfield size correction dependency factor, computed for equivalent square of the collimator opening
SSD = source to surface distance
Dmax = Dose at maximum
SPD = is the source to point distance of the point of calculation
X and Y = coordinates at a depth of the point of calculation
c = is the correction for the virtual location of the source. c is found from a plot of the inverse of the square root of measured dose versus distance from the source.
QF = is the off-axis beam quality factor
OCR = is the in air off central axis ratio value
TAR0 = is the zero (0 × 0) field TAR for the slant depth
SC = is defined as scatter contribution computed from the field size and block contours at the level of the point of calculation.
P = is the value of the penumbra, calculated by Wilkinson Source Model
Note: The above program equation does not support the use of wedge since wedge field was not used in this study.
This study will focus on verifying the percentage dose accuracy of the Area Integration Algorithm using both homogeneous and inhomogeneous insert by varying gantry angles for given field sizes of 5 cm × 5 cm and 22 cm × 25 cm, respectively.
| Materials and Methods|| |
The Plexiglas used was purchased from a local plastic shop and was machined together with a collaborative effect from the Department of Radiation Biology and Radiotherapy and the Department of Biomedical Engineering from College of Medicine of the University of Lagos. The glass dimension was 4 by 8 feet, a part of which was cut using a plastic cutter into six slabs each of dimension 30 cm × 30 cm. A plastic-based hardener (allplast) was used for holding one slab to another to form a cube. Seven holes were drilled on one face. Each drilled hole had a diameter of 2.5 cm gummed together using plastic-based hardener called “allplast.” Before the holes were drilled, the distance from the surface of the designed in-house phantom to the ionization chamber was 15 cm, while the distance between two diagonal insert were approximately 22 cm. The distances from one insert to the other (horizontally) was 7 cm and vertically were 18 cm. Furthermore, additional drilled hole was made at the top of the locally designed phantom to allow for easy filling of water and evacuation of water from the phantom. After these holes and distance have been marked out, another cylindrical rod made of Plexiglas material of thickness 0.2 mm, length 14.3 cm, and diameter 2.5 cm were fitted into the seven drilled holes and were held together at the tip by the “allplast” gum to avoid leakage [Figure 1]. Furthermore, the locally designed phantom was loaded with tissue-equivalent material putting into consideration the attenuation coefficients, electron densities, and the effective atomic numbers of each chemical composition. The reason for using chemical component with our local phantom was because of lack of availability of commercial readymade phantom [Table 1]. The locally designed phantom was scanned under a Hi-Speed NX/i computed tomography (CT)-scanner. Slices of images were acquired for six different tissue-equivalent materials. The scanned images were transferred to the precise PLAN Release 2.16 TPS for beam setup. The CT number (is also expressed in a standardized and convenient form as the Hounsfield unit [HU]) for the six different tissue equivalent materials was determined from the CT monitor [Table 2]. Images were transported through the Digital Imaging and Communications in Medicine (DICOM) to the Precise PLAN Release 2.16 TPS where twelve-field technique denoted as Beam (BM) 1–BM 12 were used with large field sizes covering the six inserts with no wedge used. The gantry angles in degree were: 0°, 22.5°, 45°, 90°, 135°, 157.5°, 180°, 197.5°, 215°, 270°, 315°, and 337.5°. The total dose for the twelve fields was 100 cGy, and the total monitor unit (MU) was 100 MU. The type of beam used was “simple.” The photon energy used was 6 MV, SAD was 100 cm and SSD was approximately 84 cm. Collimator angle was 0°, the diaphragm upper SAD was approximately 25 cm and lower SAD was 22 cm, giving at total area diaphragm size of 22 cm × 25 cm. Under modifiers, tray factor was 1, and no multileaf collimator (MLC) was present.
|Table 1: Equivalent tissue materials, elemental composition with their densities in grams per cubic centimetres|
Click here to view
|Table 2: Mimicked tissue equivalent material and their measured mean computed tomography number in Hounsfield unit|
Click here to view
A second scan was carried using the same protocol using large field sizes, but the insert was only bone equivalent material (assumed homogeneous medium). The acquired images from the CT simulator were also transferred to the Precise PLAN 2.16 TPS through a DICOM. Six-field technique was used denoted as BM 1–BM 6 covering the six inserts which were uniformly homogeneous with no wedge used. The gantry angles were: 0°, 45°, 90°, 180°, 225°, and 270°, respectively The total dose for the six fields was 100 cGy, and total MU was 100 MU. The type of beam used was “simple.” The photon energy used was 6 MV; SAD was 100 cm and SSD was approximately 84 cm collimator angle was 0°, the diaphragm upper SAD was approximately 25 cm, and lower SAD was 22 cm, giving at total area diaphragm size of 22 cm × 25 cm. Under modifiers, tray factor was 1, and no MLC was present.
A third scan was done following the same protocol with small fields. Six-field technique denoted as BM 1–BM 6 was used covering the six inserts with no wedge used. The gantry angles were: 0°, 45°, 90°, 180°, 225°, and 270°, respectively. The total dose for the six fields was 100 cGy, and total MU making up 100 MU was prescribed. The type of beam used was “simple.” Photon energy of 6 MV was used. SAD was 100 cm, and SSD was approximately 84 cm collimator angle was 0°, the diaphragm upper SAD was approximately 5 cm and lower SAD was 5 cm, giving at total area diaphragm size of 5 cm × 5 cm. Under modifiers, tray factor was 1, and no MLC was present. All planned images from the Precise PLAN 2.16 TPS were transferred to the Elekta Precise Linear Accelerator for treatment.
Furthermore, a simple experimental protocol was done to validate the accuracy of the designed trunk phantom. This was done by comparing dose values (in Gy) between the locally designed trunk phantom and solid water phantom (SWP) which is also a standard phantom in place of the nonavailable Rando Anderson Phantom which was used by the manufacturer. Photon energy of 6MV with SSD of 100 cm was used for the setup of the designed trunk phantom and the SWP. Six readings were also observed with the gantry angle at 0°.
The Elekta Precise Clinical Linear Accelerator was calibrated using a large water phantom before this study was done using a 6 MV photon beam to give 100 cGy (1 Gy) at 100 MU with a pre-calibrated NE 2570/1 Farmer-Type Ionization Chamber to determine the absorbed dose. Necessary corrections for temperature, pressure, polarization, recombination, etc., were effected on the Ionization Chamber response. Absorbed dose at reference depth was calculated as follows:
Where MQ is the electrometer reading (charge) corrected for temperature and pressure, ND, W is the chamber calibration factor and KQ, Q0 is the factor which corrects for difference in the response of the dosimeter at the calibration quality Q and at quality Q0 of the clinical X-ray beam according to the TRS 398 protocol of the IAEA.,
The deviation between the calculated and measured dose was obtained using the equation as:
Dc= Calculated dose
Dm= Measured dose
In addition to the equation 2, the introduction of the absolute value sign is to turn all negative values from measured absorbed dose into positive values.
Data analysis value was done using data analysis was done using SPSS 16.0 (SPSS Inc, Chicago, IL, USA). Descriptive statistics and unpaired t-test were implored at 95% level of significance. P <0.05 was considered statistically significant.
| Results|| |
A summary of the CT numbers also known as the HU for the mimicked chemical composition for liver, lung, muscle, adipose, bone, and kidney was 42 HU,–840 HU, 47 HU,–104 HU, 1480 HU, and 38 HU, respectively [Table 2].
The gantry angle, measured absorbed dose (Gy) and % deviation for twelve beams with six inhomogeneous inserts with field size of 22 cm × 25 cm at 0° was 1.0004, 22.5° was 0.9981, 45° was 0.9974, 90° was 1.0034, 135° was 0.976, 157.5° was 1.0082, 180° was 0.9693, 197.5° was 0.9722, 215° was 0.9820, 270° was 0.9661, 315° was 0.9708 and 337.5° was 0.9799 Gy respectively and corresponding % deviation was 0.04, 0.19, 0.26, 0.34, 2.40, 0.82, 3.07, 2.78, 1.80, 3.39, 2.92, and 2.01, respectively [Table 3].
|Table 3: Measured absorbed dose (Gy) and percentage deviation for twelve beams with six in homogeneous inserts with field size of 22 cm × 25 cm|
Click here to view
The gantry angle, measured absorbed dose (Gy), and % deviation for six beams with bone only homogeneous inserts with field size of 22 cm × 25 cm at 0° was 1.0293, 45° was 1.0174, 90° was 1.0234, 180° was 1.0113, 225° was 1.0127, 270° was 0.9915 Gy respectively and corresponding % deviation was 2.93, 1.74, 2.34, 1.13, 1.27, and 0.85, respectively [Table 4].
|Table 4: Measured absorbed dose (Gy) and percentage deviation for six beams with bone only in homogeneous inserts with field size of 22 cm × 25 cm|
Click here to view
The gantry angle, measured absorbed dose (Gy), and % deviation for six beams with six inhomogeneous inserts with field size of 5 cm × 5 cm at 0° was 1.0056, 45° was 0.9996, 90° was 0.9844, 180° was 0.9684, 225° was 0.9723, 270° was 0.9755 Gy respectively and corresponding % deviation was 0.56, 0.04, 1.56, 3.16, 2.77, and 2.45, respectively [Table 5].
|Table 5: Measured absorbed dose (Gy) and percentage deviation for six beams with six in homogeneous inserts with field size of 5 cm × 5 cm|
Click here to view
The comparison was made between the mean dose (Gy) of the locally designed pelvic phantom and SWP directly from the linear accelerator at a gantry angle of zero degree (0°) to further validate our result. The mean doses were 0.6065 and 0.6274 Gy, respectively, and % deviation that exists between them was 3.50% [Table 6].
|Table 6: Percentage deviation between solid water phantom and locally designed trunk phantom at gantry angle of 0°|
Click here to view
| Discussion|| |
A study has been carried out to verify the performance and accuracy of a TPS which uses an Area Integration Algorithm (Irregular Field Program) based on Clarkson Integration.,,,,, The results were within ± 5% as recommended by ICRU  and were consistent with Van Dyk whose variation was within ± 4%, Mijnheer and Brahme results were within 3%−3.5% which is consistent with our study., This study was also in line with Akpochafor et al. whose maximum % deviation was 4%, against 3.5% which was determined in this study.
Accuracy and verification of % deviation by the manufacturer using Rando Phantom (Alderson Research Laboratories, 390 Ludlow St., Stamford Conn. 06904) showed an overall accuracy 3.3% using parallel opposed fields in regions of uniform dose distribution without beam compensators, this result was closely in line with this study's parallel opposed field (at 90° and 270°) % deviation, which showed an overall accuracy 3.4%, with inhomogeneous inserts. Better accuracy was seen with bone only homogeneous inserts with uniform dose distribution. Our result was lower (1.7%) compared to the manufacturer verification (3.3%) in terms of accuracy. This variation might be attributed to the comparison made using homogeneity against inhomogeneity of the manufacturer material for verification. In the same vein, using small field sizes, verification of % deviation by the manufacturer showed an overall accuracy 1.8% using parallel opposed field in regions of uniform dose distribution with inhomogeneous insert, this result was better than that obtained in this study whose overall accuracy was 2.5%.
Maximum % deviation of 3.39% was noticed at a gantry angle of 270° (opposed field) and the least deviation 0.04 was noticed at a gantry angle of 0° (anterior field) for six inhomogeneous inserts at field size of 22 cm × 25 cm. In addition, maximum % deviation of 2.93% was noticed for bone only homogeneous inserts at gantry angle of 0° (anterior field) and least deviation 0.85% was noticed at a gantry angle of 270° (opposed field) with a large field size of 22 cm × 25 cm. Furthermore, the maximum % deviation 3.16% was noticed at a gantry angle of 180° and least the maximum % deviation 0.04% was noticed at a gantry angle of 45° with a small field size of 5 cm × 5 cm. This result indicated that with small field at 180°, the accelerator couch will contribute and affect percentage deviation. There was no statistically significant difference in % deviation between inhomogeneous and homogeneous inserts (P = 0.592). Similarly, there was no statistically significant difference in % deviation between large and small field sizes (P = 0.891) [Table 3], [Table 4], [Table 5].
In conclusion, the highest deviation was noticed when the plastic water phantom was compared to the local trunk phantom. The % deviation between both materials was 3.5% at a gantry angle of 0°. It was observed that one side of the parallel opposed field (270°) for six inhomogeneous insert with large field size, anterior field (0°) for bone only homogeneous insert with large field size and posterior field (180°) for six inhomogeneous insert with small field size showed the highest deviations. An overview of this study showed that the locally designed phantom was within approximately ± 4% [Table 6].
| Conclusion|| |
A low cost locally designed trunk phantom has been made for use in radiotherapy centers in Nigeria with about six centers using the same Precise PLAN 2.16 TPS. The validation of our result using a standard SWP in place of Rando Alderson Phantom which was not available was within acceptable range; it was seen to be below ICRU minimal limit. The designed trunk phantom showed an overall accuracy of ± 4% with the Area Integration algorithm of the TPS which also falls within the acceptable range of ± 5% set by ICRU. The designed phantom was made of inexpensive materials that are readily available with almost equal density with some organs in the body. The phantom proves highly useful for quality assurance and control test in radiotherapy.
Financial support and sponsorship
Conflicts of interest
There are no conflicts of interest.
| References|| |
Absorbed Dose Determination in External Beam Radiotherapy. IAEA Technical Reports Series. No. 398. Vienna, Austria: IAEA Publications; 2000.
International Commission on Radiation Units and Measurements. Determination of Absorbed Dose in Patient Irradiated by Means of X or Gamma Rays in Radiotherapy Procedures. ICRU Report 24. Bethesda: ICRU Publications; 1977.
Alam R, Ibbott GS, Pourang R, Nath R. Application of AAPM radiation therapy committee task group 23 test package for comparison of two treatment planning systems for photon external beam radiotherapy. Med Phys 1997;24:2043-54.
Rice RK, Mijnheer BJ, Chin LM. Benchmark measurements for lung dose corrections for X-ray beams. Int J Radiat Oncol Biol Phys 1988;15:399-409.
O'connor JE. The variation of scattered x-rays with density in an irradiated body. Phys Med Biol 1957;1:352-69.
Woo MK, Cunningham JR. The validity of the density scaling method in primary electron transport for photon and electron beams. Med Phys 1990;17:187-94.
Keall P, Hoban P. Accounting for primary electron scatter in x-ray beam convolution calculations. Med Phys 1995;22:1413-8.
Miften M, Wiesmeyer M, Krippner K. A fermi-eyges based 3d kernel scaling model to account for primary electron transport in superposition photon dose calculations. Med Phys 1998;25:A188.
Mackie TR, el-Khatib E, Battista J, Scrimger J, Van Dyk J, Cunningham JR, et al.
Lung dose corrections for 6-and 15-MV x rays. Med Phys 1985;12:327-32.
Fraass BA. Quality assurance for 3-D treatment planning. In: Palta J, Mackie TR, editos. Teletherapy: Present and Future. Madison: Advanced Medical Publishing; 1996. p. 253-318.
Fraass B, Doppke K, Hunt M, Kutcher G, Starkschall G, Stern R, et al.
American association of physicists in medicine radiation therapy committee task group 53: Quality assurance for clinical radiotherapy treatment planning. Med Phys 1998;25:1773-829.
Institute of Physics and Engineering in Medicine (IPEM). Physics Aspects of Quality Control in Radiotherapy. IPEM Report 81. New York: IPEM; 1999.
International Commission on Radiation Units and Measurements (ICRU). Dose Specifications for Reporting External Beam Therapy with Photons and Electrons. ICRU Report 29. Baltimore: Bethesda; 1978.
Jacky J, White CP. Testing a 3-D radiation therapy planning program. Int J Radiat Oncol Biol Phys 1990;18:253-61.
Cygler J, Ross J. Electron dose distributions in an anthropomorphic phantom – Verification of theraplan treatment planning algorithm. Med Dosim 1988;13:155-8.
Johns HE, Cunningham JR. The Physics of Radiology. 3rd
ed. Charles C. Thomas Publisher; 1971. p. 362-3.
Khan FM. The Physics of Radiation Therapy. Baltimore MD: Williams and Wilkins; 1984. p. 321, 787-94.
Khan FM. The Physics of Radiation Therapy. Baltimore: Lippincott Williams and Wilkins; 2003.
Mijnheer BJ, Battermann JJ, Wambersie A. What degree of accuracy is required and can be achieved in photon and neutron therapy? Radiother Oncol 1987;8:237-52.
Podgorsak EB. Radiation Oncology Physics: A Handbook for Teachers and Students. Vienna: IAEA Publication; 2005.
Shaw JE, editor. A Guide to Commissioning and Quality Control of Treatment Planning Systems. The Institution of Physics and Engineering in Medicine and Biology; 1994.
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:261-73.
Cunningham JR, Shrivastava PN, Wilkinson JM. Computer calculation of dose within irregularly shaped beam. In: Dosimetry Workshop, Hodgkin's Disease. Chicago ill: Radiological Physics Center; 1970.
Cunningham JR. Scatter-air ratios. Phys Med Biol 1972;17:42-51.
Clarkson JR. Note on depth doses in fields of irregular shapes. Br J Radiol 1941;14:265-8.
International Atomic Energy Agency (IAEA). Absorbed Dose Determination In External Beam Radiotherapy: An International Code of Practice for Dosimetry Based on Standards of Absorbed dose to Water. IAEA Technical Report Series No. 398. Vienna (Austria); 2000.
Brahme A, Chavaudra J, Landberg T, Mc Cullough EC, Nusslin F, Rawlinson JA, et al
. Accuracy requirements and quality assurance of external beam therapy with photons and electrons. Acta Oncol (Stockholm) 1988;27 Suppl 1:1-76.
Akpochafor MO, Omojola AD, Adeneye SO, Aweda MA, Oniyangi MS, Iloputaife CC. Verification of an irregular field algorithm of a treatment planning system using a locally designed pelvic phantom: A simple design low-cost phantom suitable for quality assurance and control test in radiotherapy. Int J Health Allied Sci 2017;6:39-44. [Full text]
[Table 1], [Table 2], [Table 3], [Table 4], [Table 5], [Table 6]