Modelado de colimadores multilámina en un sistema de planificación de radioterapia con haces externos de fotones. / Modeling multileaf collimators in radiotherapy planning system with external photons.

Alarcón Paredes, Alberto (2016) Modelado de colimadores multilámina en un sistema de planificación de radioterapia con haces externos de fotones. / Modeling multileaf collimators in radiotherapy planning system with external photons. Maestría en Física Médica, Universidad Nacional de Cuyo, Instituto Balseiro.

[img]
Vista previa
PDF (Tesis)
Español
18Mb

Resumen en español

La determinacón precisa de la transmisión de rayos X a través de un colimador multiláminas, es importante para la correcta planificación del tratamiento, especialmente en técnicas con intensidad modulada. Este tipo de técnicas presentan tiempos de irradiación mayores, respecto a técnicas convencionales. Por lo tanto, la dosis debido a la radiación de fuga a través de colimador multiláminas, puede llegar a niveles significativos de dosis en regiones que no se desea irradiar, como ser los órganos en riesgo donde la dosis debe minimizarse. El objetivo de este trabajo fue modelar distribuciones de fluencia energética para todo tipo de forma de campo generados por el colimador multiláminas y/o colimadores secundarios. Para cumplir con este objetivo se desarrollo un modelo para el cálculo de la fluencia energética considerando la influencia del sistema de colimación total, para el acelerador lineal Varian Clinac iX de 6 MV con colimadores multiláminas Varian Millenium 120, usando un algoritmo tipo ray tracing para el calculo de la atenuación y el modelo de fuentes virtuales para el modelado de la fluencia energética. El paso crucial en el desarrollo del algoritmo de ray tracing fue considerar la exacta geometría de las tres principales láminas que componen el colimador multiláminas, con detalles geométricos tales como el agujero del tornillo de conducción, el carril de soporte y las puntas redondeadas de las láminas. Esta geometría es útil para encontrar las intersecciones de rayos con las representaciones poliédricas de cada tipo de lámina, utilizando procedimientos geométricos, tales como, la intersección de una recta con planos y cilindros. Además, se utilizaron técnicas de geometría computacional, por ejemplo, el método de winding number para analizar, con alta eficiencia computacional, si las intersecciones encontradas pertenecen a los límites establecidos de las láminas. Las intersecciones encontradas sirven para calcular la distancia recorrida del rayo dentro de una o varias láminas y así poder calcular la atenuación correspondiente. La evaluación de este procedimiento para una cantidad elevada de rayos produce un mapa de atenuación propia del colimador multilámina, ubicado en el plano isocéntrico. Este mapa de atenuación y las ventanas de atenuación producidas por los colimadores secundarios representan la atenuación total del sistema de colimación del acelerador. Para modelar la fluencia energética, se implementó el modelo de fuentes virtuales basados en trabajos previos y re-validado en este trabajo. La distribución de fluencia energética final considerando el tamaño finito de las fuentes es obtenida, convolucionando el mapa de atenuación con las distribuciones de intensidad, propias del modelo de fuentes virtuales. Todos los algoritmos fueron desarrollados en Python, para su posterior optimización. Las distribuciones de fluencia energética obtenidas con el modelo propuesto fueron verificados con el código de simulación Monte Carlo BEAMnrc, mediante la validación de perfiles de fluencia para campos cuadrados e irregulares formados por los colimadores secundarios y el colimador multiláminas. Los resultados muestran que existe una buena correlación entre el modelo y Monte Carlo. La comparación de los perfiles de radiación de fuga entre láminas y transmisión a través de las puntas redondeadas, predijeron con buena precisión los detalles esperados. También, se simularon distribuciones de intensidad modulada usando técnicas estáticas con buenos resultados cualitativos. Basados en los resultados obtenidos con el algoritmo presentado en este trabajo, se puede afirmar que es una herramienta útil para cálculo de fluencia energética para cualquier forma o tamaño de campo, con varias posibles aplicaciones clínicas, incluyendo el control de calidad del colimador multiláminas y el continuo desarrollo de modelos para el cálculo de la dosis basados en fluencia energética. Se deja como trabajos futuros el acoplamiento de este método a técnicas de Convolución/Superposición para el cálculo de dosis y la implementación de algoritmos de optimización de posicionamiento de las láminas para realizar planificación inversa asociada a técnicas moduladas.

Resumen en inglés

The precise determination of X-ray transmission through a multileaf collimator is important for a correct treatment planning, especially in intensity modulated techniques, given the longer irradiation time compared to conventional techniques. Therefore, this type of technique can produce significant dose levels in protected regions, such as the organs at risk where the dose must be minimized, due to leakage radiation through the multileaf collimators. This work’s aim was to model the energy fluence distributions for every type of radiation field generated by the multileaf collimator and/or secondary collimators. In order to reach this goal, an energy fluence model for the Varian Clinac iX 6 MV with Millenium 120 multileaf collimator system was developed. The implemented model employs a ray tracing algorithm and includes geometric details, such as the screw hole for leaf conduction, support rail, and rounded end tips. The crucial step in the development of the ray tracing algorithm was to consider the exact geometry of the three main leaf type that make up the multileaf collimator, with geometric details such as the hole of the driving screw, the support rail and the rounded tips of the leaf. This geometry is useful for finding intersections of the rays with the polyhedral representations of each type of leaf, using geometric procedures, such as the intersection of a straight line with planes and cylinders. In addition, computational geometry techniques were used, for example, the winding number method, to analyze, with high computational efficiency, if the intersections found belonged to the geometric limits of the leaves. The intersections found are used to calculate the distance traveled by the ray within one or several leaves and thus calculate the corresponding attenuation. The evaluation of this procedure for a high amount of rays produces an attenuation map of the multileaf collimator. In order to model the energy fluence, a virtual source model was implemented based on previous works, which parameterize the various components that are present in the accelerator’s head, and allow to model the energy fluence. In turn, using the method of ray tracing, an algorithm was developed which calculates the path length travelled by the rays inside the multileaf collimator using the exact geometry of the leaves. This results allows obtaining a total attenuation map generated by the jaws and the multileaf collimator in the isocentric plane considering a point source of radiation. The final energy fluence distribution considering the sources’finite size is obtained, convolving the attenuation map with the intensity distributions, characteristic of the virtual source model. All the algorithms were developed in Python, for later optimization. As an initial code validation, the leavesındividual geometries were visualized through the rays´ıntersection points with each leaf. This test’s result shows a good visual correlation between the implemented geometry and the leavesactual shape. The energy fluence distributions obtained with the proposed model were verified with the Monte Carlo BEAMnrc simulation code, by validating fluence profiles for square and irregular fields formed by the jaws and the multileaf collimator. The results show a good correlation between the model and Monte Carlo simulations. The comparison between the interleaf, intraleaf leakage and transmission through rounded end tips accurately predicted the expected details. Also, modulated intensity distributions were simulated using static techniques with good qualitative results. Based on the results obtained with the algorithm presented in this work, it can be affirmed that it is a useful tool for calculating energy fluence for any shape or field size, with several possible clinical applications, including the multileaf collimator’s quality assurance and the continuous development of models for dose calculation based on energy fluence. Prospective work includes the coupling of this method to Convolution/Superposition techniques for dose calculation, and the implementation of algorithms for the optimization of leaves positioning to perform inverse planning associated with modulated techniques

Tipo de objeto:Tesis (Maestría en Física Médica)
Información Adicional:Este trabajo fue realizado en: Laboratorio de Física Médica computacional del Centro Atómico Bariloche. Área Temática: Física médica computacional.
Palabras Clave:Collimators; Colimadores; Radiotherapy; Radioterapia; Photons; Fotones; [Multileaf collimator; Colimador multilámina]
Referencias:[1] World Health Organization et al. National cancer control programmes: policies and managerial guidelines. Medical Heathcare, 2002. [2] R Roots and S Okada. Protection of dna molecules of cultured mammalian cells from radiation-induced single-strand scissions by various alcohols and sh compounds. International Journal of Radiation Biology and Related Studies in Physics, Chemistry and Medicine, 21(4):329–342, 1972. [3] Benedick A Fraass. The development of conformal radiation therapy. Medical physics, 22(11):1911–1921, 1995. [4] David I Thwaites and John B Tuohy. Back to the future: the history and development of the clinical linear accelerator. Physics in medicine and biology, 51(13):R343, 2006. [5] Godfrey N Hounsfield. Computerized transverse axial scanning (tomography): Part 1. description of system. The British journal of radiology, 46(552):1016–1022, 1973. [6] Arthur Boyer, Peter Biggs, James Galvin, Eric Klein, Thomas LoSasso, Daniel Low, Katherine Mah, and Cedric Yu. Basic applications of multileaf collimators. report of the aapm radiation therapy committee task group no. 50. Med Phys, 2001. [7] Anders Brahme. Optimization of radiation therapy and the development of multileaf collimation. International Journal of Radiation Oncology* Biology* Physics, 25(2):373–375, 1993. [8] M Saiful Huq, Indra J Das, Todd Steinberg, and James M Galvin. A dosimetric comparison of various multileaf collimators. Physics in medicine and biology, 47(12):N159, 2002. [9] Emily Heath and Jan Seuntjens. Development and validation of a beamnrc component module for accurate monte carlo modelling of the varian dynamic millennium multileaf collimator. Physics in Medicine and Biology, 48(24):4045, 2003. [10] Yan Chen, Arthur L Boyer, and C-M Ma. Calculation of x-ray transmission through a multileaf collimator. Medical physics, 27(8):1717–1726, 2000. [11] Jun Deng, Todd Pawlicki, Yan Chen, Jinsheng Li, Steve B Jiang, and CM Ma. The mlc tongue-and-groove effect on imrt dose distributions. Physics in Medicine and Biology, 46(4):1039, 2001. [12] JM Galvin. The multileaf collimator: past, present and future. Med. Phys, 26, 1999. [13] CM Ma, JS Li, T Pawlicki, SB Jiang, J Deng, MC Lee, T Koumrian, M Luxton, and S Brain. A monte carlo dose calculation tool for radiotherapy treatment planning. Physics in Medicine and Biology, 47(10):1671, 2002. [14] M Saiful Huq, Yan Yu, Zong-Ping Chen, and N Suntharalingam. Dosimetric characteristics of a commercial multileaf collimator. Medical physics, 22(2):241–247, 1995. [15] Thomas J Jordan and Peter C Williams. The design and performance characteristics of a multileaf collimator. Physics in medicine and biology, 39(2):231, 1994. [16] Arthur L Boyer and Shidong Li. Geometric analysis of light-field position of a multileaf collimator with curved ends. Medical physics, 24(5):757–762, 1997. [17] Mark R Arnfield, Jeffrey V Siebers, Jong O Kim, QiuwenWu, Paul J Keall, and Radhe Mohan. A method for determining multileaf collimator transmission and scatter for dynamic intensity modulated radiotherapy. Medical physics, 27(10):2231–2241, 2000. [18] Thomas LoSasso, Chen-Shou Chui, and C Clifton Ling. Physical and dosimetric aspects of a multileaf collimation system used in the dynamic mode for implementing intensity modulated radiotherapy. Medical physics, 25(10):1919–1927, 1998. [19] Xiaohong Wang, Spiridon Spirou, Thomas LoSasso, J¨org Stein, Chen-Shou Chui, and Radhe Mohan. Dosimetric verification of intensity-modulated fields. Medical physics, 23(3):317–327, 1996 [20] Mgter. Diego Carrasco. Desarrollo de un m´etodo de c´alculo del tipo Convoluci´on/ Superposici´on y modelado de estructuras para haces externos de fotones. Instituto Balseiro, 2015. [21] Mgter. Mar´ıa Julieta Irazoqui. Modelado, an´alisis y evaluaci´on de magnitudes involucradas en el c´alculo de dosis mediante m´etodos de convoluci´on/superposici´on para haces externos de fotones de alta energ´ıa. Instituto Balseiro, 2014. [22] Frank Herbert Attix. Introduction to radiological physics and radiation dosimetry. John Wiley & Sons, 2008. [23] Martin J Berger, JH Hubbell, SM Seltzer, J Chang, JS Coursey, Ry Sukumar, DS Zucker, and K Olsen. Xcom: photon cross sections database. NIST Standard reference database, pages 3587–3597, 1998. [24] Anders Ahnesj¨o, Lars Weber, and Per Nilsson. Modeling transmission and scatter for photon beam attenuators. Medical physics, 22(11):1711–1720, 1995. [25] Andrew J Hanson. Geometry for n-dimensional graphics. Graphics Gems IV, pages 149–170, 1994. [26] Philip Schneider and David H Eberly. Geometric tools for computer graphics. Morgan Kaufmann, 2002. [27] Paul Heckbert. Graphics Gems IV (IBM Version). Elsevier, 1994. [28] Joseph o’Rourke. Computational geometry in C. Cambridge university press, 1998. [29] John F Hughes, Andries Van Dam, James D Foley, and Steven K Feiner. Computer graphics: principles and practice. Pearson Education, 2014. [30] John M Snyder and Alan H Barr. Ray tracing complex models containing surface tessellations, volume 21. ACM, 1987. [31] Alex Crist´ov˜ao Holanda de Oliveira, Jos´e Wilson Vieira, and Fernando Roberto Andrade Lima. Monte carlo modeling of multileaf collimators using the code geant4. Brazilian Journal of Radiation Sciences, 3(1A), 2015. [32] Guido Van Rossum et al. Python programming language. 41, 2007. [33] Mario Mustra, Kresimir Delac, and Mislav Grgic. Overview of the dicom standard. In ELMAR, 2008. 50th International Symposium, volume 1, pages 39–44. IEEE, 2008. [34] DWO Rogers, BA Faddegon, GX Ding, C-M Ma, JWe, and TR Mackie. Beam: a monte carlo code to simulate radiotherapy treatment units. Medical physics, 22(5):503– 524, 1995. [35] I Kawrakow and D Rogers. The egsnrc code system: Monte carlo simulation of electron and photon transport. Medical physics, 2000. [36] DWO Rogers, BWalters, Iwan Kawrakow, et al. Beamnrc users manual. NRC Report PIRS, 509:12, 2009. [37] Robert L Cook, Thomas Porter, and Loren Carpenter. Distributed ray tracing. 18(3):137–145, 1984. [38] Robert L Siddon. Prism representation: a 3d ray-tracing algorithm for radiotherapy applications. Physics in medicine and biology, 30(8):817, 1985. [39] A Boyer, P Biggs, J Galvin, E Klein, T LoSasso, D Low, K Mah, and C Yu. Aapm report no 72: Basic applications of multileaf collimators (madison, wi, 2001. [40] IAEA. Phase-space database for external beam radiotherapy. https://wwwnds.iaea.org/phsp/photon1/, 2015. [41] Matthias Fippel, Freddy Haryanto, Oliver Dohm, Fridtjof N¨usslin, and Stephan Kriesen. A virtual photon energy fluence model for monte carlo dose calculation. Medical Physics, 30(3):301–311, 2003. [42] H Helen Liu, T Rock Mackie, and Edwin C McCullough. A dual source photon beam model used in convolution/superposition dose calculations for clinical megavoltage x-ray beams. Medical physics, 24(12):1960–1974, 1997. [43] Jong Oh Kim, Jeffrey V Siebers, Paul J Keall, Mark R Arnfield, and Radhe Mohan. A monte carlo study of radiation transport through multileaf collimators. Medical physics, 28(12):2497–2506, 2001. [44] Si Young Jang, Oleg N Vassiliev, H Helen Liu, Radhe Mohan, and Jeffrey V Siebers. Development and commissioning of a multileaf collimator model in monte carlo dose calculations for intensity-modulated radiation therapy. Medical physics, 33(3):770– 781, 2006. [45] Michael K Fix, Werner Volken, Daniel Frei, Daniel Frauchiger, Ernst J Born, and Peter Manser. Monte carlo implementation, validation, and characterization of a 120 leaf mlc. Medical physics, 38(10):5311–5320, 2011. [46] Jeffrey V Siebers, Paul J Keall, Jong Oh Kim, and Radhe Mohan. A method for photon beam monte carlo multileaf collimator particle transport. Physics in medicine and biology, 47(17):3225, 2002. [47] JPC Van Santvoort and BJM Heijmen. Dynamic multileaf collimation withouttongueand- groove´underdosage effects. Physics in medicine and biology, 41(10):2091, 1996. [48] James M Galvin, Kerry Han, and Richard Cohen. A comparison of multileafcollimator and alloy-block field shaping. International Journal of Radiation Oncology* Biology* Physics, 40(3):721–731, 1998. [49] Ann Van Esch, J¨org Bohsung, Pekka Sorvari, Mikko Tenhunen, Marta Paiusco, Mauro Iori, Per Engstr¨om, H°akan Nystr¨om, and Dominique Pierre Huyskens. Acceptance tests and quality control (qc) procedures for the clinical implementation of intensity modulated radiotherapy (imrt) using inverse planning and the sliding window technique: experience from five radiotherapy departments. Radiotherapy and Oncology, 65(1):53–70, 2002. [50] Daniel A Low, William B Harms, Sasa Mutic, and James A Purdy. A technique for the quantitative evaluation of dose distributions. Medical physics, 25(5):656–661, 1998.
Materias:Medicina > Radioterapia
Código ID:580
Depositado Por:Tamara Cárcamo
Depositado En:19 Abr 2017 15:16
Última Modificación:19 Abr 2017 15:30

Personal del repositorio solamente: página de control del documento