Dosimetria e análise de incertezas em braquiterapia ginecológica

Mestrado em Radiações Aplicadas às tecnologias da Saúde

Monte Carlo simulations for dosimetric verification in photon and electron beam radiotherapy

Dissertação para obtenção do Grau de Doutor em Engenharia Biomédica

Model-based Segmentation and Fusion of 3D Computed Tomography and 3D Ultrasound of the Eye for Radiotherapy Planning

Computed Tomography (CT) represents the standard imaging modality for tumor volume delineation for radiotherapy treatment planning of retinoblastoma despite some inherent limitations. CT scan is very useful in providing information on physical density for dose calculation and morphological volumetric information but presents a low sensitivity in assessing the tumor viability. On the other hand, 3D ultrasound (US) allows a highly accurate definition of the tumor volume thanks to its high spatial resolution but it is not currently integrated in the treatment planning but used only for diagnosis and follow-up. Our ultimate goal is an automatic segmentation of gross tumor volume (GTV) in the 3D US, the segmentation of the organs at risk (OAR) in the CT and the registration of both modalities. In this paper, we present some preliminary results in this direction. We present 3D active contour-based segmentation of the eye ball and the lens in CT images; the presented approach incorporates the prior knowledge of the anatomy by using a 3D geometrical eye model. The automated segmentation results are validated by comparing with manual segmentations. Then, we present two approaches for the fusion of 3D CT and US images: (i) landmark-based transformation, and (ii) object-based transformation that makes use of eye ball contour information on CT and US images.

Accidental and experimental Closantel intoxication in Uruguayan sheep

Abstract An outbreak of Closantel intoxication in sheep in Uruguay is described. The outbreak occurred in a group of 1300 weaning lambs treated orally with a 10% solution of Closantel. One hundred forty eight lambs showed clinical signs of intoxication and 14 died. The clinical signs included mydriasis, nystagmus, and negative pupillary reflex, bilateral blindness, bump into objects, and lateral movement of the head. No macroscopic lesions were observed. The histological lesions of the retina were cytoplasmic vacuolization in ganglion cells and in cells of the inner and outer nuclear layers with different degrees of atrophy. Vacuolization and axonal degeneration were observed in the optic nerve, with multifocal areas of fibrosis and infiltration by lymphocytes and Gitter cells. To reproduce the intoxication, four sheep were given two, four and 10 times the therapeutic dose of Closantel (0.1g/kg of BW). Only the animals receiving 10 times the recommended dose showed clinical signs. The histological examination of the lesions in experimental sheep showed similar results to those described in the accidental outbreak, except for the absence of optic nerve fibrosis and inflammation, characterizing an acute phase. Axonal myelin sheaths loss, fibroblasts and collagen fibers were observed in the ultrastructural study of the optic nerve of accidental intoxicated animals. The optic nerve of experimentally intoxicated animals had vacuoles that separated the myelin sheaths of axons. To prevent outbreaks it is suggested to weigh the animals before Closantel administration to avoid errors in dose calculation.

É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.

Validation of the Swiss Monte Carlo Plan for a static and dynamic 6 MV photon beam

Monte Carlo (MC) based dose calculations can compute dose distributions with an accuracy surpassing that of conventional algorithms used in radiotherapy, especially in regions of tissue inhomogeneities and surface discontinuities. The Swiss Monte Carlo Plan (SMCP) is a GUI-based framework for photon MC treatment planning (MCTP) interfaced to the Eclipse treatment planning system (TPS). As for any dose calculation algorithm, also the MCTP needs to be commissioned and validated before using the algorithm for clinical cases. Aim of this study is the investigation of a 6 MV beam for clinical situations within the framework of the SMCP. In this respect, all parts i.e. open fields and all the clinically available beam modifiers have to be configured so that the calculated dose distributions match the corresponding measurements. Dose distributions for the 6 MV beam were simulated in a water phantom using a phase space source above the beam modifiers. The VMC++ code was used for the radiation transport through the beam modifiers (jaws, wedges, block and multileaf collimator (MLC)) as well as for the calculation of the dose distributions within the phantom. The voxel size of the dose distributions was 2mm in all directions. The statistical uncertainty of the calculated dose distributions was below 0.4%. Simulated depth dose curves and dose profiles in terms of [Gy/MU] for static and dynamic fields were compared with the corresponding measurements using dose difference and γ analysis. For the dose difference criterion of ±1% of D(max) and the distance to agreement criterion of ±1 mm, the γ analysis showed an excellent agreement between measurements and simulations for all static open and MLC fields. The tuning of the density and the thickness for all hard wedges lead to an agreement with the corresponding measurements within 1% or 1mm. Similar results have been achieved for the block. For the validation of the tuned hard wedges, a very good agreement between calculated and measured dose distributions was achieved using a 1%/1mm criteria for the γ analysis. The calculated dose distributions of the enhanced dynamic wedges (10°, 15°, 20°, 25°, 30°, 45° and 60°) met the criteria of 1%/1mm when compared with the measurements for all situations considered. For the IMRT fields all compared measured dose values agreed with the calculated dose values within a 2% dose difference or within 1 mm distance. The SMCP has been successfully validated for a static and dynamic 6 MV photon beam, thus resulting in accurate dose calculations suitable for applications in clinical cases.

An efficient framework for photon Monte Carlo treatment planning

Currently photon Monte Carlo treatment planning (MCTP) for a patient stored in the patient database of a treatment planning system (TPS) can usually only be performed using a cumbersome multi-step procedure where many user interactions are needed. This means automation is needed for usage in clinical routine. In addition, because of the long computing time in MCTP, optimization of the MC calculations is essential. For these purposes a new graphical user interface (GUI)-based photon MC environment has been developed resulting in a very flexible framework. By this means appropriate MC transport methods are assigned to different geometric regions by still benefiting from the features included in the TPS. In order to provide a flexible MC environment, the MC particle transport has been divided into different parts: the source, beam modifiers and the patient. The source part includes the phase-space source, source models and full MC transport through the treatment head. The beam modifier part consists of one module for each beam modifier. To simulate the radiation transport through each individual beam modifier, one out of three full MC transport codes can be selected independently. Additionally, for each beam modifier a simple or an exact geometry can be chosen. Thereby, different complexity levels of radiation transport are applied during the simulation. For the patient dose calculation, two different MC codes are available. A special plug-in in Eclipse providing all necessary information by means of Dicom streams was used to start the developed MC GUI. The implementation of this framework separates the MC transport from the geometry and the modules pass the particles in memory; hence, no files are used as the interface. The implementation is realized for 6 and 15 MV beams of a Varian Clinac 2300 C/D. Several applications demonstrate the usefulness of the framework. Apart from applications dealing with the beam modifiers, two patient cases are shown. Thereby, comparisons are performed between MC calculated dose distributions and those calculated by a pencil beam or the AAA algorithm. Interfacing this flexible and efficient MC environment with Eclipse allows a widespread use for all kinds of investigations from timing and benchmarking studies to clinical patient studies. Additionally, it is possible to add modules keeping the system highly flexible and efficient.

Photon-beam subsource sensitivity to the initial electron-beam parameters

One limitation to the widespread implementation of Monte Carlo (MC) patient dose-calculation algorithms for radiotherapy is the lack of a general and accurate source model of the accelerator radiation source. Our aim in this work is to investigate the sensitivity of the photon-beam subsource distributions in a MC source model (with target, primary collimator, and flattening filter photon subsources and an electron subsource) for 6- and 18-MV photon beams when the energy and radial distributions of initial electrons striking a linac target change. For this purpose, phase-space data (PSD) was calculated for various mean electron energies striking the target, various normally distributed electron energy spread, and various normally distributed electron radial intensity distributions. All PSD was analyzed in terms of energy, fluence, and energy fluence distributions, which were compared between the different parameter sets. The energy spread was found to have a negligible influence on the subsource distributions. The mean energy and radial intensity significantly changed the target subsource distribution shapes and intensities. For the primary collimator and flattening filter subsources, the distribution shapes of the fluence and energy fluence changed little for different mean electron energies striking the target, however, their relative intensity compared with the target subsource change, which can be accounted for by a scaling factor. This study indicates that adjustments to MC source models can likely be limited to adjusting the target subsource in conjunction with scaling the relative intensity and energy spectrum of the primary collimator, flattening filter, and electron subsources when the energy and radial distributions of the initial electron-beam change.

Monte Carlo source model for photon beam radiotherapy: photon source characteristics

A major barrier to widespread clinical implementation of Monte Carlo dose calculation is the difficulty in characterizing the radiation source within a generalized source model. This work aims to develop a generalized three-component source model (target, primary collimator, flattening filter) for 6- and 18-MV photon beams that match full phase-space data (PSD). Subsource by subsource comparison of dose distributions, using either source PSD or the source model as input, allows accurate source characterization and has the potential to ease the commissioning procedure, since it is possible to obtain information about which subsource needs to be tuned. This source model is unique in that, compared to previous source models, it retains additional correlations among PS variables, which improves accuracy at nonstandard source-to-surface distances (SSDs). In our study, three-dimensional (3D) dose calculations were performed for SSDs ranging from 50 to 200 cm and for field sizes from 1 x 1 to 30 x 30 cm2 as well as a 10 x 10 cm2 field 5 cm off axis in each direction. The 3D dose distributions, using either full PSD or the source model as input, were compared in terms of dose-difference and distance-to-agreement. With this model, over 99% of the voxels agreed within +/-1% or 1 mm for the target, within 2% or 2 mm for the primary collimator, and within +/-2.5% or 2 mm for the flattening filter in all cases studied. For the dose distributions, 99% of the dose voxels agreed within 1% or 1 mm when the combined source model-including a charged particle source and the full PSD as input-was used. The accurate and general characterization of each photon source and knowledge of the subsource dose distributions should facilitate source model commissioning procedures by allowing scaling the histogram distributions representing the subsources to be tuned.

Monte Carlo simulation of a dynamic MLC: implementation and applications

PURPOSE: Study of behavior and influence of a multileaf collimator (MLC) on dose calculation, verification, and portal energy spectra in the case of intensity-modulated fields obtained with a step-and-shoot or a dynamic technique. METHODS: The 80-leaf MLC for the Varian Clinac 2300 C/D was implemented in a previously developed Monte Carlo (MC) based multiple source model (MSM) for a 6 MV photon beam. Using this model and the MC program GEANT, dose distributions, energy fluence maps and energy spectra at different portal planes were calculated for three different MLC applications. RESULTS: The comparison of MC-calculated dose distributions in the phantom and portal plane, with those measured with films showed an agreement within 3% and 1.5 mm for all cases studied. The deviations mainly occur in the extremes of the intensity modulation. The MC method allows to investigate, among other aspects, dose components, energy fluence maps, tongue-and-groove effects and energy spectra at portal planes. CONCLUSION: The MSM together with the implementation of the MLC is appropriate for a number of investigations in intensity-modulated radiation therapy (IMRT).