Desarrollo de un programa de simulación basado en el método de Montecarlo para el cálculo de dosis con maniquíes divididos en voxels. Aplicaciones en tomografía computarizada

Abstract

Among the different radiodiagnostic techniques, computed tomography (CT) represents the major contribution to the collective dose. In order to have available a tool to assess and ascertain the imparted dose upon CT examinations, we have developed and validated a method to calculate radiation absorbed dose in CT examinations from images of phantoms and standard patients by using a voxel-based Monte Carlo simulation method.<br/>Simulations and measurements of radiation dose were conducted in air, within cylindrical phantoms of both circular and elliptic cross sections, and in an anthropomorphic phantom, with the simulations being made on volumes based on their own CT images. The radiation beam geometry (focus-to-axis distance, field of view, collimation, and primary and beam shaper filtration) and the X-ray spectral distribution for each scanner under study were implemented for the simulation. Simulations were performed for both single axial scan and a series of contiguous scans covering the complete phantom-object volume. All the axial CT images were reduced to a 256 x 256 matrix and stacked in a volume. Patient images were segmented prior to simulation of radiation transport by using three categories (lung, muscle and bone). To test the method, values of CT specific dose quantities for CT exposure simulations were calculated and compared to those derived from measurements performed free-in-air and within cylindrical phantoms. To complete the scope of the work, series of CT scans of the head, neck, and trunk of an anthropomorphic phantom and standard patients were simulated. <br/>Using the EGS4 code, written in an extended FORTRAN called Mortran3, we have estimated first the radiation dose in the individual voxels; radiation transport was simulated by means of the EGS4 code system. The simulation of photon transport included incoherent and coherent scattering, photoelectric absorption with creation of fluorescent photons or Auger electrons. The photon interaction cross sections and the electronic density for each material or compound were prepared by the PEGS4 pre-processor. The electron interaction processes considered were inelastic collisions with atomic electrons, scattering effects by the nucleus, and bremsstrahlung losses that are negligible for the energies considered. The low energy cutoff for photon transport was 5 keV, whereas for electrons was fixed at 30 keV, with a ESTEPE (fraction of the energy lost in a step) value of 1%. The conversion of energy fluence to dose rate was made using the mass energy-absorption coefficients from Hubbell and Seltzer. The number of histories in each simulation was 3×106 per each 1 cm slice.<br/>Secondly, the dose CT indices for peripheral and central positions and the dose to bladder were estimated as well, the later being used to test organ doses in the anthropomorphic phantom. The dose measurements were performed by using a "pencil" ionisation chamber and 30 thermoluminiscent dosimeters. The comparison between simulated and measured dose data shows differences below 5% in all cases.<br/>The designed program allows contrasted results to be achieved, and might be an adequate tool for the analysis of doses imparted to individual patients with different models of CT scanner.

De las diferentes técnicas usadas con fines radiodiagnósticos, las exploraciones con equipos de tomografía computarizada (TC) suponen la mayor contribución a la dosis colectiva, ya que las dosis impartidas en cada exploración de TC son significativamente más elevadas que las de los estudios radiológicos convencionales. Con el fin de disponer de una herramienta para conocer las dosis impartidas a los pacientes en estudios de TC, se ha desarrollado y validado un programa de simulación para estimar la dosis absorbida en volúmenes divididos en voxels usando técnicas de Montecarlo. Para obtener los "objetos" tridimensionales divididos en voxels que representen a cada maniquí o paciente, se han transformado las propias imágenes de TC mediante la segmentación en materiales u órganos, y posteriormente se han apilado a lo largo del eje longitudinal. El haz se ha simulado teniendo en cuenta los espectros de emisión de rayos X, la filtración total del equipo (filtración primaria y filtro de forma) y reproduciendo la geometría concreta de irradiación para los tres equipos utilizados (dos helicoidales y un multicorte). Mediante el código EGS4 se han reproducido los procesos relevantes (interacción y transporte de las partículas) en el rango de energías de los equipos de TC, y se ha estimado la energía depositada en cada voxel. <br/>Durante el proceso de validación se han realizado estudios de dosis absorbida en aire, en maniquíes cilíndricos de sección circular y elíptica, y en un maniquí antropomórfico, construidos a partir de sus propias imágenes de TC. Se han comparado los resultados de la simulación con los de medidas realizadas en maniquíes físicos utilizando una cámara de ionización o dosimetría por termoluminiscencia. El análisis de las diferencias y de las incertidumbres atribuidas a los resultados ha mostrado que ambas son compatibles. <br/>Una vez validado el programa de simulación se ha aplicado a tres estudios concretos, permitiendo estimar la dosis promedio por rodaja, y la dosis promedio equivalente en distintos órganos de dos maniquíes antropomórficos y de seis pacientes seleccionados para representar las distintas morfologías de ambos sexos. También se han asignado equivalencias entre maniquíes y grupos de pacientes pediátricos de edades comprendidas entre los 0 a los 15 años. Finalmente se ha aplicado el método a un caso real de protección radiológica mediante la utilización de protectores con bismuto destinados a reducir la dosis en ojos, tiroides y mamas.<br/>El programa diseñado, que se puede adaptar a cada estudio particular, permite obtener resultados contrastados y puede ser una herramienta adecuada para analizar las dosis impartidas a pacientes individuales en cada protocolo y equipo de TC concreto. Además, el programa se muestra como herramienta eficaz para analizar los parámetros que influyen en la distribución de dosis, ya sean geométricos, del haz, o relacionados con los modos de adquisición del equipo.