Experimental evaluation of MCDTK, the Monte Carlo DICOM Tool-Kit

The Monte Carlo DICOM Tool-Kit (MCDTK) is a software suite designed for treatment plan dose verification, using the BEAMnrc and DOSXYZnrc Monte Carlo codes. MCDTK converts DICOM-format treatment plan information into Monte Carlo input files and compares the results of Monte Carlo treatment simulations with conventional treatment planning dose calculations. In this study, a treatment is planned using a commercial treatment planning system, delivered to a pelvis phantom containing ten thermoluminescent dosimeters and simulated using BEAMnrc and DOSXYZnrc using inputs derived from MCDTK. The dosimetric accuracy of the Monte Carlo data is then evaluated via comparisons with the dose distribution obtained from the treatment planning system as well as the in-phantom point dose measurements. The simulated beam arrangement produced by MCDTK is found to be in geometric agreement with the planned treatment. An isodose display generated from the Monte Carlo data by MCDTK shows general agreement with the isodose display obtained from the treatment planning system, except for small regions around density heterogeneities in the phantom, where the pencil-beam dose calculation performed by the treatment planning systemis likely to be less accurate. All point dose measurements agree with the Monte Carlo data obtained using MCDTK, within confidence limits, and all except one of these point dose measurements show closer agreement with theMonte Carlo data than with the doses calculated by the treatment planning system. This study provides a simple demonstration of the geometric and dosimetric accuracy ofMonte Carlo simulations based on information from MCDTK.

Radiotherapy Monte Carlo simulation using cloud computing technology

Cloud computing allows for vast computational resources to be leveraged quickly and easily in bursts as and when required. Here we describe a technique that allows for Monte Carlo radiotherapy dose calculations to be performed using GEANT4 and executed in the cloud, with relative simulation cost and completion time evaluated as a function of machine count. As expected, simulation completion time decreases as 1=n for n parallel machines, and relative simulation cost is found to be optimal where n is a factor of the total simulation time in hours. Using the technique, we demonstrate the potential usefulness of cloud computing as a solution for rapid Monte Carlo simulation for radiotherapy dose calculation without the need for dedicated local computer hardware as a proof of principal. Funding source Cancer Australia (Department of Health and Ageing) Research Grant 614217

Faster Monte Carlo simulation of radiotherapy geometries

Using Monte Carlo simulation for radiotherapy dose calculation can provide more accurate results when compared to the analytical methods usually found in modern treatment planning systems, especially in regions with a high degree of inhomogeneity. These more accurate results acquired using Monte Carlo simulation however, often require orders of magnitude more calculation time so as to attain high precision, thereby reducing its utility within the clinical environment. This work aims to improve the utility of Monte Carlo simulation within the clinical environment by developing techniques which enable faster Monte Carlo simulation of radiotherapy geometries. This is achieved principally through the use new high performance computing environments and simpler alternative, yet equivalent representations of complex geometries. Firstly the use of cloud computing technology and it application to radiotherapy dose calculation is demonstrated. As with other super-computer like environments, the time to complete a simulation decreases as 1=n with increasing n cloud based computers performing the calculation in parallel. Unlike traditional super computer infrastructure however, there is no initial outlay of cost, only modest ongoing usage fees; the simulations described in the following are performed using this cloud computing technology. The definition of geometry within the chosen Monte Carlo simulation environment - Geometry & Tracking 4 (GEANT4) in this case - is also addressed in this work. At the simulation implementation level, a new computer aided design interface is presented for use with GEANT4 enabling direct coupling between manufactured parts and their equivalent in the simulation environment, which is of particular importance when defining linear accelerator treatment head geometry. Further, a new technique for navigating tessellated or meshed geometries is described, allowing for up to 3 orders of magnitude performance improvement with the use of tetrahedral meshes in place of complex triangular surface meshes. The technique has application in the definition of both mechanical parts in a geometry as well as patient geometry. Static patient CT datasets like those found in typical radiotherapy treatment plans are often very large and present a significant performance penalty on a Monte Carlo simulation. By extracting the regions of interest in a radiotherapy treatment plan, and representing them in a mesh based form similar to those used in computer aided design, the above mentioned optimisation techniques can be used so as to reduce the time required to navigation the patient geometry in the simulation environment. Results presented in this work show that these equivalent yet much simplified patient geometry representations enable significant performance improvements over simulations that consider raw CT datasets alone. Furthermore, this mesh based representation allows for direct manipulation of the geometry enabling motion augmentation for time dependant dose calculation for example. Finally, an experimental dosimetry technique is described which allows the validation of time dependant Monte Carlo simulation, like the ones made possible by the afore mentioned patient geometry definition. A bespoke organic plastic scintillator dose rate meter is embedded in a gel dosimeter thereby enabling simultaneous 3D dose distribution and dose rate measurement. This work demonstrates the effectiveness of applying alternative and equivalent geometry definitions to complex geometries for the purposes of Monte Carlo simulation performance improvement. Additionally, these alternative geometry definitions allow for manipulations to be performed on otherwise static and rigid geometry.

The development of Monte Carlo techniques for the verification of radiotherapy treatments

Introduction: Recent advances in the planning and delivery of radiotherapy treatments have resulted in improvements in the accuracy and precision with which therapeutic radiation can be administered. As the complexity of the treatments increases it becomes more difficult to predict the dose distribution in the patient accurately. Monte Carlo (MC) methods have the potential to improve the accuracy of the dose calculations and are increasingly being recognised as the ‘gold standard’ for predicting dose deposition in the patient [1]. This project has three main aims: 1. To develop tools that enable the transfer of treatment plan information from the treatment planning system (TPS) to a MC dose calculation engine. 2. To develop tools for comparing the 3D dose distributions calculated by the TPS and the MC dose engine. 3. To investigate the radiobiological significance of any errors between the TPS patient dose distribution and the MC dose distribution in terms of Tumour Control Probability (TCP) and Normal Tissue Complication Probabilities (NTCP). The work presented here addresses the first two aims. Methods: (1a) Plan Importing: A database of commissioned accelerator models (Elekta Precise and Varian 2100CD) has been developed for treatment simulations in the MC system (EGSnrc/BEAMnrc). Beam descriptions can be exported from the TPS using the widespread DICOM framework, and the resultant files are parsed with the assistance of a software library (PixelMed Java DICOM Toolkit). The information in these files (such as the monitor units, the jaw positions and gantry orientation) is used to construct a plan-specific accelerator model which allows an accurate simulation of the patient treatment field. (1b) Dose Simulation: The calculation of a dose distribution requires patient CT images which are prepared for the MC simulation using a tool (CTCREATE) packaged with the system. Beam simulation results are converted to absolute dose per- MU using calibration factors recorded during the commissioning process and treatment simulation. These distributions are combined according to the MU meter settings stored in the exported plan to produce an accurate description of the prescribed dose to the patient. (2) Dose Comparison: TPS dose calculations can be obtained using either a DICOM export or by direct retrieval of binary dose files from the file system. Dose difference, gamma evaluation and normalised dose difference algorithms [2] were employed for the comparison of the TPS dose distribution and the MC dose distribution. These implementations are spatial resolution independent and able to interpolate for comparisons. Results and Discussion: The tools successfully produced Monte Carlo input files for a variety of plans exported from the Eclipse (Varian Medical Systems) and Pinnacle (Philips Medical Systems) planning systems: ranging in complexity from a single uniform square field to a five-field step and shoot IMRT treatment. The simulation of collimated beams has been verified geometrically, and validation of dose distributions in a simple body phantom (QUASAR) will follow. The developed dose comparison algorithms have also been tested with controlled dose distribution changes. Conclusion: The capability of the developed code to independently process treatment plans has been demonstrated. A number of limitations exist: only static fields are currently supported (dynamic wedges and dynamic IMRT will require further development), and the process has not been tested for planning systems other than Eclipse and Pinnacle. The tools will be used to independently assess the accuracy of the current treatment planning system dose calculation algorithms for complex treatment deliveries such as IMRT in treatment sites where patient inhomogeneities are expected to be significant. Acknowledgements: Computational resources and services used in this work were provided by the HPC and Research Support Group, Queensland University of Technology, Brisbane, Australia. Pinnacle dose parsing made possible with the help of Paul Reich, North Coast Cancer Institute, North Coast, New South Wales.

CTCombine : rotating and combining CT data with an accurate detector model to simulate radiotherapy portal imaging at non-zero beam angles

Established Monte Carlo user codes BEAMnrc and DOSXYZnrc permit the accurate and straightforward simulation of radiotherapy experiments and treatments delivered from multiple beam angles. However, when an electronic portal imaging detector (EPID) is included in these simulations, treatment delivery from non-zero beam angles becomes problematic. This study introduces CTCombine, a purpose-built code for rotating selected CT data volumes, converting CT numbers to mass densities, combining the results with model EPIDs and writing output in a form which can easily be read and used by the dose calculation code DOSXYZnrc...

IMP 重离子治癌中的剂量计算方法

Basic algorithms of biological effective dose optimization and dose distribution on CT image for the heavy ion therapy project at the Institute of Modern Physics(IMP),Chinese Academy of Sciences(CAS) are reported in this paper.Firstly,biological effective dose optimization is conducted in water.According to the relationship between CT number and water equivalent path length,an integral algorithm is used to calculate the average dose within a pixel and then the dose distribution in tissue is derived.Secondly...中文文摘：针对深部肿瘤重离子治疗临床试验的需求,首先在水介质中进行生物有效剂量的优化计算,然后根据CT图像中像素CT值与水等效长度转换系数之间的关系,结合水中的深度剂量分布曲线对每个像素进行积分得到CT图像上的生物有效剂量分布。同时介绍了基于被动式束流配送系统适形照射时的剂量确定方式,并提出二维适形放疗也应使用分层照射方式以适应治疗时的不同要求。这些方法适合目前及今后在IMP进行的重离子治癌临床试验研究中治疗计划系统的需要。

X射线在介质中的剂量分布计算

研究了放射治疗中X射线在介质中的输运过程,编程实现了基于蒙特卡罗方法的剂量计算．并在便于图形处理的软件Matlab中对光子输运结果进行了可视化处理．对X射线在均匀介质和菲均匀介质中的蒙特卡罗模拟结果与实测结果、其他蒙特卡罗软件模拟结果进行了比较,结果符合较好．实验结果表明该方法既可以获得很快的仿真速度,又能得到精确直观的剂量计算结果,为提高放射治疗水平具有重要的指导意义和应用价值．

Clinical implementation of dynamic multileaf collimation for compensated bladder treatments.

BACKGROUND AND PURPOSE: To describe the clinical implementation of dynamic multileaf collimation (DMLC). Custom compensated four-field treatments of carcinoma of the bladder have been used as a simple test site for the introduction of intensity modulated radiotherapy.MATERIALS AND METHODS: Compensating intensity modulations are calculated from computed tomography (CT) data, accounting for scattered, as well as primary radiation. Modulations are converted to multileaf collimator (MLC) leaf and jaw settings for dynamic delivery on a linear accelerator. A full dose calculation is carried out, accounting for dynamic leaf and jaw motion and transmission through these components. Before treatment, a test run of the delivery is performed and an absolute dose measurement made in a water or solid water phantom. Treatments are verified by in vivo diode measurements and real-time electronic portal imaging. RESULTS: Seven patients have been treated using DMLC. The technique improves dose homogeneity within the target volume, reducing high dose areas and compensating for loss of scatter at the beam edge. A typical total treatment time is 20 min. CONCLUSIONS: Compensated bladder treatments have proven an effective test site for DMLC in an extremely busy clinic.

IMRT for image-guided single vocal cord irradiation

PURPOSE: We have been developing an image-guided single vocal cord irradiation technique to treat patients with stage T1a glottic carcinoma. In the present study, we compared the dose coverage to the affected vocal cord and the dose delivered to the organs at risk using conventional, intensity-modulated radiotherapy (IMRT) coplanar, and IMRT non-coplanar techniques.

METHODS AND MATERIALS: For 10 patients, conventional treatment plans using two laterally opposed wedged 6-MV photon beams were calculated in XiO (Elekta-CMS treatment planning system). An in-house IMRT/beam angle optimization algorithm was used to obtain the coplanar and non-coplanar optimized beam angles. Using these angles, the IMRT plans were generated in Monaco (IMRT treatment planning system, Elekta-CMS) with the implemented Monte Carlo dose calculation algorithm. The organs at risk included the contralateral vocal cord, arytenoids, swallowing muscles, carotid arteries, and spinal cord. The prescription dose was 66 Gy in 33 fractions.

RESULTS: For the conventional plans and coplanar and non-coplanar IMRT plans, the population-averaged mean dose ± standard deviation to the planning target volume was 67 ± 1 Gy. The contralateral vocal cord dose was reduced from 66 ± 1 Gy in the conventional plans to 39 ± 8 Gy and 36 ± 6 Gy in the coplanar and non-coplanar IMRT plans, respectively. IMRT consistently reduced the doses to the other organs at risk.

CONCLUSIONS: Single vocal cord irradiation with IMRT resulted in good target coverage and provided significant sparing of the critical structures. This has the potential to improve the quality-of-life outcomes after RT and maintain the same local control rates.

Étude des artefacts en tomodensitométrie par simulation Monte Carlo

En radiothérapie, la tomodensitométrie (CT) fournit l’information anatomique du patient utile au calcul de dose durant la planification de traitement. Afin de considérer la composition hétérogène des tissus, des techniques de calcul telles que la méthode Monte Carlo sont nécessaires pour calculer la dose de manière exacte. L’importation des images CT dans un tel calcul exige que chaque voxel exprimé en unité Hounsfield (HU) soit converti en une valeur physique telle que la densité électronique (ED). Cette conversion est habituellement effectuée à l’aide d’une courbe d’étalonnage HU-ED. Une anomalie ou artefact qui apparaît dans une image CT avant l’étalonnage est susceptible d’assigner un mauvais tissu à un voxel. Ces erreurs peuvent causer une perte cruciale de fiabilité du calcul de dose. Ce travail vise à attribuer une valeur exacte aux voxels d’images CT afin d’assurer la fiabilité des calculs de dose durant la planification de traitement en radiothérapie. Pour y parvenir, une étude est réalisée sur les artefacts qui sont reproduits par simulation Monte Carlo. Pour réduire le temps de calcul, les simulations sont parallélisées et transposées sur un superordinateur. Une étude de sensibilité des nombres HU en présence d’artefacts est ensuite réalisée par une analyse statistique des histogrammes. À l’origine de nombreux artefacts, le durcissement de faisceau est étudié davantage. Une revue sur l’état de l’art en matière de correction du durcissement de faisceau est présentée suivi d’une démonstration explicite d’une correction empirique.

Cálculo de doses em exames de TC abdómen-pélvico em dois hospitais da região de Lisboa

Este trabalho tem como objectivo avaliar se os valores de dose recebida pelos doentes que realizam exames de Tomografia Computorizada (TC) Abdómen-pélvicos em dois hospitais da grande Lisboa estão de acordo com os Níveis de Referência de Diagnóstico Europeus (NRD). Foram comparadas também as unidades dosimétricas de exames realizados com corrente modelada e de exames realizados com corrente contínua em ambos os hospitais. O estudo consistiu na recolha de dados de 200 exames de TC Abdómen-pélvicos, 100 em cada hospital. No Hospital A, a média dos valores de DLPfoi 562,34 mGy.cm e a média dos valores de CTDIvol foi 12,06 mGy. No Hospital B, a média dos valores DLPe CTDIvol foi 767,14 mGy.cm e 15,02 mGy, respectivamente. No total da nossa amostra, concluiu-se que em nenhum exame os valores de unidades dosimétricas ultrapassam os valores dos NRD. Verificou-se também que os exames realizados no Hospital B implicam em média uma maior dose de radiação ionizante para o doente.

Controle de qualidade em radioterapia: dosimetria de fótons em acelerador linear – Mevatron MXT

The treatment of a tumor with ionizing radiation is an ongoing process with well differentiated stages. These ones include the tumor diagnosis and location, the decision on the treatment strategy, the absorbed dose planning and calculation, the treatment administration, the absorbed dose verification and the evaluation of results in short and long terms. The quality of a radiotherapy procedure is closely linked to factors that may be classified as clinical, such as the diagnosis, the tumor location, the treatment strategy chosen and the continuous treatment reassessment; dosimetric or physical, such as the uncertainty in the dose calculation, its optimization and verification, the suitability of the equipment to provide a radiation beam consistent with the treatment planning; finally, others which are related to the practical application of radiotherapy treatment and the handling of the patient. In order to analyze the radiotherapy quality, one should realize that the three aspects (medical, physical or dosimetric and practical application) should be considered in a combined way. This means that numerous actions of the radiotherapists, medical physicists and technicians in radiotherapy should be held jointly and their knowledge level will significantly affect the treatment quality. In this study, the main physical parameters used in dosimetry are defined as well as determined experimentally for a linear accelerator Mevatron - MXT. With this, it is intended to provide recommendations for the physical aspects of Quality Assurance (QA) in the radiotherapy treatments, and these will usually be applied by professionals in Medical Physics. In addition to these instructions, it is recommended that additional texts are prepared to address in detail the clinical aspects of the treatments QA

In-vivo-Tumorbiologie mittels 68 Ga-PET: Beurteilung der p-Glykoprotein-Aktivität in Weichgewebe-Tumoren sowie Diagnose und Therapie von Knochenmetastasen

Krebs ist eine der häufigsten Krankheiten und stellt eine der wichtigsten medizinischen Herausforderungen des 21. Jahrhunderts dar. Eine frühzeitige Diagnose ist dabei essentiell für eine individuell angepasste Therapie zur Verbesserung der Lebensqualität und -erwartung der Patienten. Hierbei kommen der 68Ge/68Ga-Generator und das daraus resultierende PET-Nuklid 68Ga immer stärker in den Fokus von Wissenschaft und Medizin. rnrnFür eine erfolgreiche Therapie stellt die Chemoresistenz (Multi-Drug-Resistance) zahlreicher Tumore eine schwerwiegende Komplikation dar. Für das Therapieversagen ist die Aktivierung des Transportproteins p-Glykoprotein (pGP) maßgeblich mit verantwortlich. Mit Hilfe der Schiff’schen Base [68Ga]MFL6.MZ konnte die Aktivitätsänderung von pGP unter verschiedener Beeinflussung erstmals in vivo beobachtet werden. So zeigte sich, dass sich unter azidotischen Bedingungen in Tumoren die Aktivität des pGP erhöht und somit vermehrt auch Zytostatika, die pGP-Substrate sind, aus den Tumoren transportiert werden. Durch Aufklärung der Abhängigkeit der pGP-Aktivität von dessen Signalkaskade konnte gezeigt werden, dass durch eine Blockade der MAP-Kinase p38 eine Erniedrigung der pGP-Aktivität zu verzeichnen ist. Die ebenfalls in der Signalkaskade eingebundene MAP-Kinase ERK1/2 hingegen spielt hier nur eine untergeordnete Rolle.rnrnNeben dem Versagen der Chemotherapie stellt auch die Metastasierung eines Malignoms massive Einschnitte in die Lebensqualität von Erkrankten dar. Befallen die Metastasen das Skelett eines Menschen, wird dies zumeist erst spät registriert. 68Ga-markierte Bisphosphonate bieten nun die Möglichkeit, Patienten quantitativ auf Knochenmetastasen hin untersuchen zu können. So konnten zu Beginn einfache Phosphonate wie EDTMP und DOTP nicht die nötige in vivo Stabilität bzw. hohe radiochemische Ausbeuten liefern und sind damit für die Anwendung am Menschen uninteressant. Jedoch die DOTA-basierten Bisphosphonate allen voran der Ligand BPAMD zeigen ein großes Potential. In vivo-Versuche an Ratten mit Knochenmetastasen zeigten, dass sich [68Ga]BPAMD an den Metastasen anreichert und einen sehr guten Kontrast zum gesunden Knochen darstellt. Der Tracer konnte erstmals am Menschen angewendet werden und zeigte in ausgewählten Regionen eine höhere Anreicherung als eine zuvor durchgeführte PET-Aufnahme mit [18F]Fluorid. Der Ligand BPAMD bietet außerdem den Vorteil, neben 68Ga auch andere dreiwertige Radionuklide wie das therapeutische 177Lu komplexieren zu können. Durch Studien zur Komplexbildung und Stabilität konnte auch [177Lu]BPAMD in der klinischen Anwendung erprobt werden und zeigte eine Anreicherung an den Knochenmetastasen. So ist es nun möglich, Knochenmetastasen mittels 68Ga-PET zu diagnostizieren, eine entsprechende Dosisberechnung anzustellen und anschließend mit dem gleichen Liganden eine Therapie mit [177Lu]BPAMD durchzuführen.

VMC++ validation for photon beams in the energy range of 20-1000 keV

In high energy teletherapy, VMC++ is known to be a very accurate and efficient Monte Carlo (MC) code. In principle, the MC method is also a powerful dose calculation tool in other areas in radiation oncology, e.g., brachytherapy or orthovoltage radiotherapy. However, VMC++ is not validated for the low-energy range of such applications. This work aims in the validation of the VMC++ MC code for photon beams in the energy range between 20 and 1000 keV.

The Swiss IMRT dosimetry intercomparison using a thorax phantom

In 2008, a national intensity modulated radiation therapy (IMRT) dosimetry intercomparison was carried out for all 23 radiation oncology institutions in Switzerland. It was the aim to check the treatment chain focused on the planning, dose calculation, and irradiation process.