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

Méthode efficace d'assignation de tissus humains par tomodensitométrie à double énergie

Pour analyser les images en tomodensitométrie, une méthode stœchiométrique est gé- néralement utilisée. Une courbe relie les unités Hounsfield d’une image à la densité électronique du milieu. La tomodensitométrie à double énergie permet d’obtenir des informations supplémentaires sur ces images. Une méthode stœchiométrique a été dé- veloppée pour permettre de déterminer les valeurs de densité électronique et de numéro atomique effectif à partir d’une paire d’images d’un tomodensitomètre à double énergie. Le but de cette recherche est de développer une nouvelle méthode d’identification de tissus en utilisant ces paramètres extraits en tomodensitométrie à double énergie. Cette nouvelle méthode est comparée avec la méthode standard de tomodensitométrie à simple énergie. Par ailleurs, l’impact dosimétrique de bien identifier un tissu est déterminé. Des simulations Monte Carlo permettent d’utiliser des fantômes numériques dont tous les paramètres sont connus. Les différents fantômes utilisés permettent d’étalonner les méthodes stœchiométriques, de comparer la polyvalence et la robustesse des méthodes d’identification de tissus double énergie et simple énergie, ainsi que de comparer les distributions de dose dans des fantômes uniformes de mêmes densités, mais de compo- sitions différentes. La méthode utilisant la tomodensitométrie à double énergie fournit des valeurs de densi- tés électroniques plus exactes, quelles que soient les conditions étudiées. Cette méthode s’avère également plus robuste aux variations de densité des tissus. L’impact dosimé- trique d’une bonne identification de tissus devient important pour des traitements aux énergies plus faibles, donc aux énergies d’imagerie et de curiethérapie.

Determination of Very Slow Diffusion of Paramagnetic Ions by NMR Spectroscopy

In this paper, we propose a new method of measuring the very slow paramagnetic ion diffusion coefficient using a commercial high-resolution spectrometer. If there are distinct paramagnetic ions influencing the hydrogen nuclear magnetic relaxation time differently, their diffusion coefficients can be measured separately. A cylindrical phantom filled with Fricke xylenol gel solution and irradiated with gamma rays was used to validate the method. The Fricke xylenol gel solution was prepared with 270 Bloom porcine gelatin, the phantom was irradiated with gamma rays originated from a (60)Co source and a high-resolution 200 MHz nuclear magnetic resonance (NMR) spectrometer was used to obtain the phantom (1)H profile in the presence of a linear magnetic field gradient. By observing the temporal evolution of the phantom NMR profile, an apparent ferric ion diffusion coefficient of 0.50 mu m(2)/ms due to ferric ions diffusion was obtained. In any medical process where the ionizing radiation is used, the dose planning and the dose delivery are the key elements for the patient safety and success of treatment. These points become even more important in modern conformal radio therapy techniques, such as stereotactic radiosurgery, where the delivered dose in a single session of treatment can be an order of magnitude higher than the regular doses of radiotherapy. Several methods have been proposed to obtain the three-dimensional (3-D) dose distribution. Recently, we proposed an alternative method for the 3-D radiation dose mapping, where the ionizing radiation modifies the local relative concentration of Fe(2+)/Fe(3+) in a phantom containing Fricke gel and this variation is associated to the MR image intensity. The smearing of the intensity gradient is proportional to the diffusion coefficient of the Fe(3+) and Fe(2+) in the phantom. There are several methods for measurement of the ionic diffusion using NMR, however, they are applicable when the diffusion is not very slow.

Experimental evidence and model explanation for cell population characteristics modification when applying sequential photodynamic therapy

We present experimental evidence of the existence of cell variability in terms of threshold light dose for Hep G2 (liver cancer cells) cultured. Using a theoretical model to describe the effects caused by successive photodynamic therapy (PDT) sessions, and based on the consequences of a partial response we introduce the threshold dose distribution concept within a tumor. The experimental model consists in a stack of flasks, and simulates subsequent layers of a tissue exposed to PDT application. The result indicates that cells from the same culture could respond in different ways to similar PDT induced-damages. Moreover, the consequence is a partial killing of the cells submitted to PDT, and the death fraction decreased at each in vitro PDT session. To demonstrate the occurrence of cell population modification as a response to PDT, we constructed a simple theoretical model and assumed that the threshold dose distribution for a cell population of a tumor is represented by a modified Gaussian distribution.

A complete dosimetric characterization of two Sr-90-Y-90 dermatologic applicators

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Two dimensional computer simulation of plasma immersion ion implantation

The biggest advantage of plasma immersion ion implantation (PIII) is the capability of treating objects with irregular geometry without complex manipulation of the target holder. The effectiveness of this approach relies on the uniformity of the incident ion dose. Unfortunately, perfect dose uniformity is usually difficult to achieve when treating samples of complex shape. The problems arise from the non-uniform plasma density and expansion of plasma sheath. A particle-in-cell computer simulation is used to study the time-dependent evolution of the plasma sheath surrounding two-dimensional objects during process of plasma immersion ion implantation. Before starting the implantation phase, steady-state nitrogen plasma is established inside the simulation volume by using ionization of gas precursor with primary electrons. The plasma self-consistently evolves to a non-uniform density distribution, which is used as initial density distribution for the implantation phase. As a result, we can obtain a more realistic description of the plasma sheath expansion and dynamics. Ion current density on the target, average impact energy, and trajectories of the implanted ions were calculated for three geometrical shapes. Large deviations from the uniform dose distribution have been observed for targets with irregular shapes. In addition, effect of secondary electron emission has been included in our simulation and no qualitative modifications to the sheath dynamics have been noticed. However, the energetic secondary electrons change drastically the plasma net balance and also pose significant X-ray hazard. Finally, an axial magnetic field has been added to the calculations and the possibility for magnetic insulation of secondary electrons has been proven.

Avaliação das condições de uso de um equipamento de tomografia computadorizada e otimização de protocolos de crânio empregados na rotina do HCFMB-UNESP

Almost forty years computed tomography (CT) has been one of the most powerful tools in diagnostic imaging. However, this modality delivers relatively high doses to their patients. It is known that the inappropriate use and unnecessary radiation may be associated with a significant risk of cancer, especially in pediatric patients. Moreover, the quality assurance in CT, provided and required by Portaria 453/98 and the guide of the RE 1016/05, ensures that the images generated by computer tomography provide reliable diagnostic information with doses as low as reasonably achievable. This research aimed to make the quality control (QC) of CT equipment; establish a better relationship between dose and noise on the image to protocols of skull CT according to the study of optimization proposed in 2005 by Daros; and assess the dose distribution in different cranial organs for protocols of adult and pediatric use in the routine of the Department of Diagnostic Imaging of HCFMB-UNESP. The equipment used for testing QC, optimization and dosimetry was a third generation tomograph GE Sytec 3000i

Análise quantitativa dos parâmetros físicos e radiométricos de procedimetnos radioterapicos em tumores de mama

The paper presents the radiometric parameters determined by the medical physicist during routine radiotherapy planning service in cases of breast cancer . The contours of the breast volume in patients undergoing radiation breast tumors at the Hospital das Clinicas, Faculty of Medicine , UNESP, Botucatu ( HCFMB ) during the year 2012 were analyzed . In order to analyze the influence of physical and radiometric parameters for the determination of the dose distribution of irradiated breast volume , four measurements of isodose curves were prepared in four different heights breast , and compared with the isodose curves plotted computationally . In the routine of planning , the medical physicist must determine the isodose curve that gives the best dose distribution homogeneity in the irradiated volume . The choice of the treatment plan can be done by dedicated computer systems , which require significantly costly investments available services having better financial support . In the Service of Medical Physics , Department of Radiotherapy , HC FMB , we use a two-dimensional software for determination of isodose curves , however , this software is out of date and frequently becomes inoperable due to the lack of maintenance and it is a closed system without feasibility of interference from computer professionals . This fact requires manual preparation of isodose curves , which are subject to uncertainties due to the subjectivity in the clinical interpretation of medical radiation oncologist and medical physicist responsible for planning , plus dispendiar significant calculation time . The choice of the optimal isodose curve depends on the energy of the radiation beam , the geometry and dimensions of the irradiated area . The contours of the breast studied in this work evaluations showed that , for a given energy input , such as the energy of 1.25 MeV of gamma radiation Unit Telecobaltoterapia , the determination of the percentage depth dose ( PDP ) ...

Braquiterapia com sementes de iodo-125: manufatura e tratamento

According to the National Institute of Cancer – INCA, 466.730 new cancer cases will occur in Brazil in 2009. Prostate and Lung cancer in man are the most incident types (in exception of the non-melanoma cancer). The brachytherapy with 125-iodine sources is an important method of prostate cancer treatment. The implant with iodine-125 seeds uses aproximately 100 seeds that are imported impossibilitating the treatment in large scale. For this reason, a multidisciplinary group was created at the Energetic and Nuclear Research Institute – Radiation Technology Center (IPEN –CTR / SP) to develop a national 125-iodine source and established a facility for local production. The seeds manufacture in Brazil will diminish the cost of treatment and prostate brachytherapy will be offered to more patients. This work aim is to observe and discuss the methods used in seeds manufacture there are being developed in Brazil and to present an prostate cancer case folloied in A.C. Camargo Hospital. The 125-iodine is adsorbed in an silver wire, then deposited at titanium coat. The weld is made by an process of plasm sealing. The seeds goes through several test to guarantee that there is no leakage. The result is an high quality and cheaper product. The implant tecnique is an fast and save procedure. The medical physicst preplan the case to stipulate the quantity of seeds there will be used. At the dat of the implant the medical physicst replan the procedure due to changes of volume in the organ. That assure the correct dose distribution in the target. Besides, the 125-iodine low energy is absorbed at the volume of interrest saving sourronding healthy tissues such as the rectum and the urethra

Análise da viabilidade de simulações com feixes de prótons de 100, 150, 200 e 250 MeV em alvos heterogêneos

In the recent years, the use of proton beams in radiotherapy has been an outstanding progress (SMITH, 2006). Up to now, computed tomography (CT) is a prerequisite for treatment planning in this kind of therapy because it provides the electron density distribution required for calculation of dose and the interval of doses. However, the use of CT images for proton treatment planning ignores fundamental differences in physical interaction processes between photons and protons and is, therefore, potentially inaccurate (SADROZINSKI, 2004). Proton CT (pCT) can in principle directly measure the density distribution needed in a patient for the dose distribution (SCHULTE, et al, 2004). One important problem that should be solved is the implementation of image reconstruction algorithms. In this sense, it is necessary to know how the presence of materials with different density and composition interfere in the energy deposition by ionization and coulomb excitation, during its trajectory. The study was conducted in two stages, was used in both the program SRIM (The Stopping and Range of Ions in Matter) to perform simulations of the interaction of proton beams with pencil beam type. In the first step we used the energies in the range of 100-250 MeV (ZIEGLER, 1999). The targets were set to 50 mm in length for the beam of 100 MeV, due to its interaction with the target, and short-range, and 70 mm for 150, 200 and 250 MeV The target was composed of liquid water and a layer of 6 mm cortical bone (ICRP). It were made 9 simulations varying the position of the heterogeneity of 5 mm. In the second step the energy of 250 MeV was taken out from the simulations, due to its greater energy and less interaction. The targets were diminished to 50 mm thick to standardize the simulations. The layer of bone was divided into two equal parts and both were put in the ends of the target... (Complete abstract click electronic access below)

Métodos de Calibração e Dosimetria de Aplicadores de Betaterapia

The sources of betatherapy for clinical use in Brazil are, the vast majority of strontium-90, radioactive element that is not produced in the country, and therefore requires importation of international laboratories accredited by the International Atomic Energy Agency (IAEA).The use of these resources is always limited the crediting of characteristic values supplied by the manufacturer tables that provide the nominal value of activity and dose distribution to determine the irradiation time of the injury. The Institute of Nuclear Energy Research (IPEN / CNEN-SP) has recently researching the emission profile of these types of radiation sources, and some jobs are being developed with ionization chambers extrapolation for the purpose of standardizing a systematic calibration sources betatherapy. Other studies using parallel measures dosimeters (TLD's) and simulations with the Monte Carlo method. Radiological films have also been used in studies of applicators dosimetric analysis of strontium-90. This paper seeks to analyze the different methods for calibration of applicators betatherapy, already consolidated in studies by examining the advantages and disadvantages of each procedure

Advanced aspects of radiation protection in the use of particle accelerators in the medical field

In this work, the well-known MC code FLUKA was used to simulate the GE PETrace cyclotron (16.5 MeV) installed at “S. Orsola-Malpighi” University Hospital (Bologna, IT) and routinely used in the production of positron emitting radionuclides. Simulations yielded estimates of various quantities of interest, including: the effective dose distribution around the equipment; the effective number of neutron produced per incident proton and their spectral distribution; the activation of the structure of the cyclotron and the vault walls; the activation of the ambient air, in particular the production of 41Ar, the assessment of the saturation yield of radionuclides used in nuclear medicine. The simulations were validated against experimental measurements in terms of physical and transport parameters to be used at the energy range of interest in the medical field. The validated model was also extensively used in several practical applications uncluding the direct cyclotron production of non-standard radionuclides such as 99mTc, the production of medical radionuclides at TRIUMF (Vancouver, CA) TR13 cyclotron (13 MeV), the complete design of the new PET facility of “Sacro Cuore – Don Calabria” Hospital (Negrar, IT), including the ACSI TR19 (19 MeV) cyclotron, the dose field around the energy selection system (degrader) of a proton therapy cyclotron, the design of plug-doors for a new cyclotron facility, in which a 70 MeV cyclotron will be installed, and the partial decommissioning of a PET facility, including the replacement of a Scanditronix MC17 cyclotron with a new TR19 cyclotron.

Sviluppo di un metodo di radioterapia adattativa tramite immagini portali e Cone-Beam CT

Il presente lavoro, svolto presso il servizio di Fisica Sanitaria dell’Azienda Ospedaliera Universitaria di Parma, consiste nello sviluppo di un metodo innovativo di radioterapia adattativa. Il metodo è stato applicato a pazienti affetti da varie patologie, trattati con tecnica VMAT, (Volumetric Modulated Arc Therapy), altamente conformata al target. Il metodo sviluppato si compone di due fasi: nella prima fase vengono effettuate due analisi su immagini portali, di ricostruzione della dose all'isocentro e l'analisi gamma 2D. Se almeno una di queste fallisce, si interviene con la seconda fase, che vede l'acquisizione della CBCT del paziente e la taratura in densità elettronica della stessa. Si calcola dunque il piano su CBCT, previa operazione di contouring da parte del medico e, infine, si esegue l'analisi gamma 3D sulle matrici di dose calcolate sulla CT e sulla CBCT del paziente, quantificando gli indici gamma sulle strutture PTV, CTV e OAR di interesse clinico. In base ai risultati, se necessario, si può intervenire sul piano di trattamento. Le analisi gamma 2D e 3D sono state svolte avvalendosi di un software toolkit chiamato GADD-23 (Gamma Analysis on 2D and 3D Dose Distributions) implementato e sviluppato appositamente in ambiente Matlab per questo lavoro di tesi; in particolare, la realizzazione di GADD-23 è stata resa possibile grazie all'interazione con due software di tipo open-source, Elastix e CERR, specifici per l’elaborazione e la registrazione di immagini mediche. I risultati ottenuti mostrano come il metodo sviluppato sia in grado di mettere in luce cambiamenti anatomici che alcuni pazienti hanno subìto, di tipo sistematico, in cui è possibile prendere in considerazione una ripianificazione del trattamento per correggerli, o di tipo casuale, sui quali può essere utile condurre l'attenzione del medico radioterapista, sebbene non sia necessario un replanning.

Leaf transmission reduction using moving jaws for dynamic MLC IMRT

The aim of this work is to investigate to what extent it is possible to use the secondary collimator jaws to reduce the transmitted radiation through the multileaf collimator (MLC) during an intensity modulated radiation therapy (IMRT). A method is developed and introduced where the jaws follow the open window of the MLC dynamically (dJAW method). With the aid of three academic cases (Closed MLC, Sliding-gap, and Chair) and two clinical cases (prostate and head and neck) the feasibility of the dJAW method and the influence of this method on the applied dose distributions are investigated. For this purpose the treatment planning system Eclipse and the Research-Toolbox were used as well as measurements within a solid water phantom were performed. The transmitted radiation through the closed MLC leads to an inhomogeneous dose distribution. In this case, the measured dose within a plane perpendicular to the central axis differs up to 40% (referring to the maximum dose within this plane) for 6 and 15 MV. The calculated dose with Eclipse is clearly more homogeneous. For the Sliding-gap case this difference is still up to 9%. Among other things, these differences depend on the depth of the measurement within the solid water phantom and on the application method. In the Chair case, the dose in regions where no dose is desired is locally reduced by up to 50% using the dJAW method instead of the conventional method. The dose inside the chair-shaped region decreased up to 4% if the same number of monitor units (MU) as for the conventional method was applied. The undesired dose in the volume body minus the planning target volume in the clinical cases prostate and head and neck decreased up to 1.8% and 1.5%, while the number of the applied MU increased up to 3.1% and 2.8%, respectively. The new dJAW method has the potential to enhance the optimization of the conventional IMRT to a further step.

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.