Source: http://docplayer.es/1966386-Calculos-dosimetricos-mediante-codigo-monte-carlo-a-partir-de-imagenes-de-pet-ct.html
Timestamp: 2016-10-28 10:28:30+00:00

Document:
⭐CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT.
CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT.
Download "CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT."
Carmelo Venegas Martínez
1 TESIS MAESTRÍA EN FÍSICA MÉDICA CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT. Ing. Juan Camilo Bedoya Tobon Lic. Sergio L. Mosconi DIRECTOR Ing. Juan Camilo Bedoya Tobon MAESTRANDO Instituto Balseiro Comisión Nacional de Energía Atómica Universidad Nacional de Cuyo Diciembre 20112 A Dios, a mi madre, mis hermanos y mi novia3 REFERENCIAS Lista de Abreviaturas ASCII BGO CPU CT DCS FDG FEE FOV GATE Geant ICRP LMF MC MIRD MN PDF PET RT ROIs T a Código Estándar Estadounidense para el Intercambio de Información (American Standard Code for Information Interchange) Germanato de Bismuto Unidad central de procesamiento (central processing unit) Tomografía computada Sección eficaz diferencial (Differential cross section) Fluorodesoxiglucosa Front End Electronics Campo de visión (Field of views) Aplicación de Genat4 para emisión tomográfica (Geant4 Application for Emission Tomography) Geometría y seguimiento (GEometry ANd Tracking) Comisión Internacional de Protección Radiológica List Mode Format Monte Carlo Medical Internal Radiation Dose Medicina Nuclear Función de densidad de probabilidad (Probability Density Function) Tomografía por emisión positrónica Radioterapia Región de interés Tiempo de Adquisición4 ÍNDICE Summary... 7 Resumen... 9 Capítulo Introducción a la dosimetría interna Historia de la Dosimetría Interna Formalismos MIRD Tasa de dosis absorbida Energía Promedio Fracción Absorbida Factor "S" Cálculo de la Actividad Integrada Cálculo del Tiempo de Residencia Métodos actuales para realizar cálculos en Dosimetría Interna Dosimetría Interna en Diagnóstico y en Terapia con Radiofármacos Capítulo Simulaciones de Monte Carlo para el cálculo dosimétrico en PET Números Aleatorios Método de muestreo analógico o detallado Método de función de distribución (pdf) Método de rechazo Método Mixto Muestreo no analógico o condensado Métodos específicos para fotones Métodos específicos para electrones Métodos generales Transporte de fotones Transporte de electrones y positrones Colisiones elásticas e inelásticas Bremsstrahlung Aniquilación positrónica Implementación Códigos de Monte Carlo utilizados para aplicaciones de Medicina Nuclear Capítulo Plataforma de simulación GATE Características Generales5 Arquitectura del software Visualización Detalle Geometría El sistema Detectores sensibles Procesos físicos Producción de cortes Fuentes Voxelización de fuentes y fantomas Administración de tiempo y movimiento Digitalización y parámetros de lectura Experimento Cálculo de Dosis Capítulo PET/CT GE Discovery STE Especificaciones técnicas Modelo de simulación Geometría Procesos físicos y protocolos de adquisición Datos de salida y configuración de reconstrucción Experimento Capítulo Segmentación y procesamiento de imágenes Elección de órganos de estudio Órganos de mayor concentración de actividad y de mayor eliminación biológica Órganos más radiosensibles Descripción de la fuente Características físicas del flúor Proceso de simulación del fantoma Capítulo Discusión Y Resultados Tamaño de Matriz Actividad Tiempo de adquisición Algoritmo de simulación6 6.5- Segmentación Mapas Dosimétricos Incertidumbres de los niveles de dosis GATE vs. MIRD Capítulo Conclusiones Apéndice Definición de la geometría del escáner Apéndice Definición de los procesos físicos y de los sistemas de Adquisición Apéndice Librería de materiales de GATE Agradecimientos Referencias Referencias de software7 SUMMARY SUMMARY In this thesis, Carlo simulations have been configured as an alternative method for calculating organ doses as a result of irradiation from radiopharmaceuticals used for PET studies (usually FDG-18). These simulations have been made in GATE platform, free software and developed by the creators Geant4. This simulator has two great advantages over homologous platforms as PENELOPE or MNCP. The first benefit is to have an easy to understand programming language. The second advantage is to allow a realistic assessment of the patient using the information contained in CAT scans and PET studies. GATE reads the geometry of the patient, assigns the corresponding attenuation coefficients and determines the activity levels contained in the organ by voxelizing these studies and using translators range, that contain the information from these two voxel range parameters. In order to ensure a good resolving dosimetric map, tests for simulation optimization were made varying the phantom matrix size, the acquisition time, the activity and the simulation algorithms. All tests were done in three different patients. Each patient represents a different anatomical region. The regions for analysis were: pelvis, chest and head. These regions were selected because of their activity concentration. The selected body parts for the study were the organs with greater collecting activity, greater elimination and more radiation-sensivility. Dosimetric maps obtained as a result of the simulations have the same dimensions as the CT image matrix, but in a different format. In this case the images formats were float 4 bits. The GATE reported the dose voxel values, considering the voxel energy absorbed and the mass of it. Finally, we compared reported dose values by GATE and the MIRD calculated values. It was seen that MIRD method tends to underestimate in many cases the dose levels compared with the values obtained with GATE, making a false and unrealistic distribution of the activity in the specific organ. 78 Keywords: CT, PET, internal dosimetry, voxellized, maps, dosimetry, Monte Carlo 89 RESUMEN RESUMEN En esta tesis se han realizado simulaciones de Monte Carlo como método alternativo para el cálculo de dosis en órgano, producto de la irradiación de radiofarmacos usados para estudios PET, (por lo general FDG-18). Estas simulaciones han sido hechas en la plataforma GATE, un software libre, desarrollado por los creadores Geant4. Este simulador posee dos grandes ventajas frente a plataformas homólogas como PENELOPE ó MNCP, la primera de esta es poseer un lenguaje de programación fácil de entender. La segunda es permitir una evaluación realista del paciente, utilizando la información contenida en los estudios tomográficos y de PET. GATE lee la geometría del paciente, asigna los coeficientes de atenuación correspondientes y determina los niveles de actividad contenidos en órgano a través de la voxelización de dichos estudios, utilizando traductores de rango, que contienen la información de estos dos parámetros por rango de voxel. En busca de un método para la optimización de las simulaciones, garantizando una buena resolución de los mapas dosimétricos, se hicieron pruebas variando el tamaño de la matriz del fantoma, el tiempo de adquisición, los valores de actividad y los algoritmos de simulación propios del programa. Todas las pruebas fueron hechas en tres pacientes diferentes. Cada uno de los pacientes representa una región anatómica diferente. Las regiones de análisis fueron: pelvis, tórax y cabeza y cuello, regiones de gran importancia por su alta concentración de actividad. Como órganos de estudio fueron elegidos los órganos más captantes y de mayor eliminación y los mas radio-sensible. Los mapas dosimétricos obtenidos como resultado de las simulaciones poseen las mismas dimensiones que la matriz de las imágenes CT, pero formato diferente, para este caso las imágenes son tipo flotantes de 4 bits. El GATE reporta valores de dosis por voxel, evaluando la energía absorbida y la masa de cada uno de estos. Finalmente, se compararon los valores de dosis reportados por GATE, con los valores calculados con el método MIRD. De la cual se pudo concluir que el método MIRD 910 tiende a subestimar en muchos de los casos los niveles de dosis frente a los valores obtenido con GATE, ya que hace una distribución errónea y poco realista de actividad en órgano. Palabras claves: Tomografía, PET, dosimetría interna, voxelización, mapas dosimétricos, Monte Carlo 1011 CAPÍTULO 1 INTRODUCCIÓN A LA DOSIMETRÍA INTERNA Tras el avance de la tecnología y el desarrollo de nuevo equipos diagnósticos en el campo de la medicina nuclear, el uso de radiofármacos como bio-marcadores a aumentado abruptamente, y con este la necesidad de llevar un registro de la dosis recibida por el paciente en dichos exámenes. En Medicina Nuclear, la aplicación correcta de los principios de radioprotección depende siempre de un balance muy preciso entre los riesgos potenciales de exposición a las radiaciones ionizantes versus los posibles beneficios que esta exposición tenga asociados(1). En aplicaciones que involucran el uso de "emisores internos" (la administración de "fuentes abiertas" o radiofármacos que pasan a interactuar directamente con el metabolismo humano), la dosis recibida por el paciente no se puede medir de forma directa. Por lo tanto, la estimación de la dosis absorbida generalmente requiere el uso12 CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT de ecuaciones y modelos matemáticos que simulen el metabolismo humano, y más específicamente, el del paciente en estudio.(2) Por lo cual se recurre a los formalismos matemáticos desarrollados en el campo de la dosimetría interna. La dosimetría interna es "... la ciencia utilizada para medir, calcular, estimar, predecir y cuantificar la energía radiante absorbida por la ionización y excitación de los átomos en los tejidos humanos como resultado de la emisión producto de la deposición internas de radionucleido " [Raabe, 1994] Son muchos los modelos que pueden ser utilizados para el cálculo dosimétrico, pero de todos ellos el más usado, por su fácil implementación en los diferentes servicios de medicina nuclear, ha sido el sistema MIRD (Medical Internal Radiation Dose). Este método permite calcular la dosis, D, depositada en un órgano blanco desde uno o más órganos fuente en el organismo Historia de la Dosimetría Interna A continuación se puede seguir una breve reseña cronológica de los hechos que condujeron al desarrollo de la Dosimetría Interna (2)(3): 1986: Descubrimiento de la Radiactividad. 1930: Descubrimiento de los Aceleradores Lineales. 1942: Leónidas Marinelli escribe el 1er Trabajo de Dosimetría Interna titulado "Dosage Determination with Radioactive Isotopes". 1948: L. Marinelli y Edith Quimby, "Practical Considerations in therapy and protection". 1953: Loevinger propicia el 1er Simposio/clases de cálculo de dosis. 1956: "Radiation Dosimetry" 1st edition (Incluyó detalles de Dosimetría Interna). 1964: 1er Meeting de la "Nuclear Medical Society Committee". 1968: "Radiation Dosimetry" 2nd edition. 1st MIRD Pamphlet. 1969: 1st "Radiation Doses and Effects". Oak Ridge. RIDIC. 1213 CAPÍTULO 1 INTRODUCCIÓN A LA DOSIMETRÍA INTERNA Finales de los '80: Aparición de PET, el cual sumado a TAC, RMI y SPECT permite realizar una dosimetría interna paciente específica. 1998: El primer prototipo PET-TAC se desarrolló en la Universidad de Pittsburgh 1999: 16th MIRD Pamphlet. Guía principal para dosimetría interna basada en imágenes planares o tomográficas. En los últimos años con el desarrollo y la implementación de los nuevos equipos híbridos tales como el PET/CT, la dosimetría interna ha sufrido un gran crecimiento, producto de la diminución de los tiempos de estudio y la posibilidad de fusionar de una forma más precisa y sencilla estudios metabólicos y anatómicos Formalismos MIRD El sistema dosimétrico MIRD (Medical Internal Radiation Dose), fue introducido por la Sociedad de Medicina Nuclear de EEUU en 1960 como herramienta para la estimación de la dosis en órganos debido a la incorporación de material radiactivo (4)(5). A continuación se presenta la metodología MIRD para el cálculo de dosis absorbidas en distintos tejidos Tasa de dosis absorbida Para el cálculo de la tasa de dosis absorbida se considera un volumen infinito de material o tejido equivalente, que contiene una distribución uniforme de material radiactivo. Si toda la energía emitida es absorbida, podemos conocer la energía absorbida por unidad de masa y por unidad de tiempo usando el siguiente razonamiento (4)(6): (1.1) Considerando que la actividad se define como el número de transiciones por unidad de tiempo, tenemos entonces: (1.2) 1314 CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT Donde A es la actividad y E es la energía promedio por transición, suponiendo que existe más de un tipo de transición o decaimiento para el radionucleido dado, o bien la existencia de más de un radionucleido. Si toda la energía emitida es absorbida en el material, esto considerando un material de dimensiones infinitas tenemos: (1.3) Reescribiendo la ecuación (1.3), utilizando los términos propios del modelo MIRD, tenemos: (1.4) Donde = tasa de dosis absorbida, A=cantidad de actividad, m= masa del tejido, = actividad por unidad de masa o concentración, k=es la constante de proporcionalidad que dará la tasa de dosis en las unidades deseadas y = energía promedio emitida por transición nuclear. Si expresa: A (actividad)= Bq, m (masa)= kg y E (energía)= MeV, la tasa de dosis se expresa en Gy/segundo al multiplicar por k = 1.6E-13. Ahora, si son cambiadas las unidades de tal forma que A = Ci, m = g y E = MeV, La tasa de dosis quedará expresada en rad/hr multiplicando por k = Cuando un radionucleido decae, este puede emitir diferente tipos de radiación (fotones o partículas), la fracción de transformación en la cual se emite cada tipo se conoce como abundancia y es definido con el símbolo. Entonces para una dada i,. Del anterior párrafo, la tasa de dosis total para un volumen grande de tejido puede ser calculada de la siguiente manera, donde el subíndice i se refiere a cada radiación: 1415 CAPÍTULO 1 INTRODUCCIÓN A LA DOSIMETRÍA INTERNA (1.5).2.2- Energía Promedio Para cada radionucleido la energía promedio es una constante, que al ser multiplicada por la constante k, da origen a una nueva constante designada en el esquema de MIRD por la letra mayúscula griega. Así, la tasa de dosis queda expresada como: (1.6) Fracción Absorbida Como fue discutido con anterioridad el formalismo anterior solo es válido para volúmenes infinitos, en donde toda la radiación no alcanza a salir de estos e interactuar con el medio, sin embargo, el cuerpo humano y sus distintos órganos no cumplen los criterios de un volumen infinito en términos de los rangos de algunos tipos de radiación. Por esta razón se hace necesario clasificar las radiaciones en dos tipos: (a) no penetrantes, que se refiere a todas las formas de radiación que son fácilmente atenuada, es decir, la energía se deposita en las inmediaciones de la fuente y (b) penetrante, que se refiere a las radiaciones que pueden viajar largas distancias antes de interactuar y depositar sus energías. Figura 1.1 Figura 1.1. Comportamiento de los fotones, los electrones y las partículas alfa en el interior del cuerpo. (7) 1516 CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT Las partículas alfa, partículas beta, y los electrones se clasifican generalmente como radiaciones no penetrantes, debido a su corto alcance en el tejido, mientras los fotones (rayos gamma y rayos X con energías superiores a 10 kev), debido a la forma en que interactúan con los tejidos, pueden entregar parte de su energía fuera del órgano de origen y son clasificados como radiaciones penetrante. Los rayos X de energías menores a 10 kev se consideran no penetrantes. Ante esto, se ve necesario agregar un nuevo factor a la ecuación de tasa de dosis para tener en cuenta la energía emitida desde el órgano fuente, que no es absorbida en el tejido de interés u órgano blanco. Este factor representa la fracción de energía emitida por un órgano fuente que es absorbida en un órgano blanco y se llama la fracción absorbida,, y se expresa así (el blanco y la fuente pueden estar en la misma región)(6): (1.7) La tasa de dosis equivalente para un blanco que ha sido irradiado por una fuente se convierte entonces en: (1.8) Donde es la región blanco, región fuente e i representa el tipo de radiación. La tasa de dosis para un blanco considerando todas las fuentes posibles, puede ser expresada como: (1.9) Una fracción absorbida por separado se requiere para cada tipo de radiación emitida desde la fuente a un determinado blanco. Otro término introducido en el formalismo MIRD es la fracción específica absorbida, que representa la fracción absorbida por unidad de masa de los órganos 1617 CAPÍTULO 1 INTRODUCCIÓN A LA DOSIMETRÍA INTERNA blancos y se simboliza por la letra griega. La fracción específica absorbida por un blanco que es irradiado por una fuente puede ser representada por =.(6) El modelo MIRD, toma como constante las masas de los tejidos blanco, haciendo una correlación entre paciente-fantoma, por lo que la tasa de dosis de dicho órgano producto de la irradiación de un órgano fuente, queda expresada como: (1.10) La tasa de dosis total para un órgano blanco puede ser calculada como: (1.11) Factor "S" Eligiendo el modelo matemático apropiado, los parámetros, pueden ser tomados como constantes, que es denominado como el factor S y que se encuentra tabulado para los 117 radionucleido y las posibles combinaciones órgano fuente-blanco, que han sido medidas en un fantoma de un hombre de referencia. El factor S está definido como:(6) (1.12) Entonces, la ecuación de la tasa de dosis ha sido simplificada a: (1.13) En los casos particulares en los que la masa de los órganos blanco considerados se alejen significativamente de los datos del hombre o mujer de referencia, se pueden modificar añadiendo un factor de corrección, que no es más que la relación entre la masa del órgano estándar y la masa real del órgano, estimada por TAC, RMI o ecografía, o bien por palpación o por experiencia. (2) 1718 Cálculo de la Actividad Integrada CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT Si deseamos obtener la dosis absorbida por un paciente en cada órgano o tejido, debemos integrar la tasa de dosis, en general entre 0 e infinito, aunque puede ser en cualquier lapso de tiempo. Así: (1.14) En general, un radiofármaco incorporado en el metabolismo humano tiene un decaimiento diferente del físico, llamado efectivo. (1.15) El esquema MIRD usa el término actividad acumulada, para representar la integral de la actividad sobre el tiempo. es una medida del número total de desintegraciones radiactivas que ocurren durante el periodo de tiempo en el que está presente la radiación en el órgano fuente. Este intervalo de tiempo está definido, de tal forma que el tiempo inicial es el instante en el que se suministra el radionucleido, mientras que el tiempo final infinito simula cuánto tarda en desaparece el radionucleido del órgano fuente. Donde: (1.16) Con esta convención, la ecuación de la dosis absorbida en el órgano blanco k producto de la actividad en los órganos fuentes h se puede escribir como: (1.17) 1819 CAPÍTULO 1 INTRODUCCIÓN A LA DOSIMETRÍA INTERNA Cálculo del Tiempo de Residencia Con el valor de actividad inyectada y el valor de la actividad integrada, es posible calcular un valor que ha sido y es aún muy popular, llamado Tiempo de Residencia. Así: (1.18) Donde es el Tiempo de Residencia, el cual tiene unidades de horas o segundos. Los programas de cálculo de dosis, como el MIRDose, usan el valor del Tiempo de Residencia o de la Actividad Integrada para obtener las dosis en los diferentes órganos y tejidos. 1.3-Métodos actuales para realizar cálculos en Dosimetría Interna En la actualidad existen diversas formas de realizar Dosimetría Interna, distinguiéndose estas por el tipo de variables que utilizan para su estimación. Existen así: (2) La Dosimetría Interna Biológica o Genética, en la cual la estimación de dosis se realiza a través del control del número de mutaciones, aberraciones cromosómicas y otras expresiones celulares y sub-celulares, en muestras de tejidos del paciente. La Dosimetría Interna Bioquímica o Clínica, en la cual se realizan estimaciones de dosis fundamentadas en la evolución de la población de especies celulares (por ejemplo células de la sangre tales como eritrocitos, leucocitos y/o plaquetas) en el tiempo, como así también de sustancias liberadas a la sangre, en la orina, en las heces, e incluso en sudor, respiración o lágrimas, luego de la internalización de una fuente abierta. La Dosimetría Interna Estándar o Física, en la cual la estimación de dosis se realiza a través del uso de fantomas que simulan el cuerpo humano, y 1920 CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT suponiendo que las fracciones de incorporación (uptake) en cada órgano son de cierto valor específico (también estimado). La Dosimetría Interna Físico-Matemática, en la cual se suman a lo anterior, modelos matemáticos compartiméntales que tratan de simular la dinámica del metabolismo de los diferentes tejidos y su evolución en el tiempo. La Dosimetría Interna por Imágenes 2D y 3D, la cual contiene características de las anteriores, pero los datos de actividad incorporada provienen de estimaciones sobre imágenes del paciente. Estas estimaciones se vuelcan después en fantomas estándar de acuerdo a la metodología MIRD. La Dosimetría Interna por Fantomas Digitales Pseudo-Elásticos, en la cual se agrega a lo anterior la característica de que los fantomas se pueden adaptar de una forma más fiel a la anatomía del paciente, siguiendo las estructuras del mismo que se visualizan en imágenes como las de TAC o RM. Método utilizado por ejemplo en las simulaciones de Monte Carlo. La Micro-dosimetría Interna, la cual posee características del tipo anterior, pero suma además la corrección del concepto de uniformidad de distribución del radiofármaco en el tejido. Esta corrección permite discriminar a pequeña escala la distribución del radiofármaco en el tejido, y por lo tanto corregir las diferencias de actividad incorporada dentro del órgano fuente. En la práctica la dosimetría realizada en un servicio puede ser una mezcla de los tipos presentados anteriormente. Esto dependerá de los recursos con los que cuente cada centro, como así también del uso y precisión de los resultados que se deseen obtener. Esta tesis tiene como objetivo entonces realizar mediante simuladores de Monte Carlo mapas dosimétricos en base a imágenes obtenidas por el PET/CT, para poder realizar simulaciones lo más realista posible. 1.4-Dosimetría Interna en Diagnóstico y en Terapia con Radiofármacos Existen diferencias fundamentales en cuanto al tipo de información que se maneja en Dosimetría Interna aplicada a Terapia con Radiofármacos comparado con el caso del estudio de las dosis que se reciben durante una Práctica Diagnóstica (1). En este último 2021 CAPÍTULO 1 INTRODUCCIÓN A LA DOSIMETRÍA INTERNA caso se plantea la necesidad de optimizar la aplicación que involucre fuentes abiertas, sin que se altere la sensibilidad diagnóstica de la práctica, y al mismo tiempo sin que se superen los limites de dosis de algún/os tejido/s. Se deberá atender no solo el monto de actividad usada en el estudio sino también el uso de bloqueadores específicos que permitan disponer de dicha actividad solo en los órganos de interés. También se plantea la diferencia de que el valor de este tipo de dosimetría interna es solo investigativo y, una vez que ha sido determinado el protocolo a seguir, no se aplica de forma rutinaria. (2) Resumiendo, se puede decir que el concepto dosimétrico se aplica al desarrollo de un nuevo radiofármaco o de un nuevo protocolo para una Práctica Diagnóstica con el objeto de: Evaluar el impacto radiológico del nuevo radiofármaco (efectos adversos). Establecer la relación riesgo/beneficio. Establecer las estrategias comunicacionales en el equipo multidisciplinario para la aplicación en la clínica diaria. La Dosimetría Interna aplicada a Terapia con Fuentes Abiertas es de aplicación rutinaria, ya que se debe realizar un cálculo de dosis para cada paciente en cada situación clínica específica (en el caso anterior se buscan personas cercanas al hombre o mujer de referencia, o bien se toman estudios poblacionales, obteniendo finalmente resultados para la media). Además de esto, como en algunos tratamientos se requiere un rápido asesoramiento de dosis para saber la dirección terapéutica a seguir (como en el caso de los tratamientos con I131), no se pueden obtener demasiadas muestras o imágenes para llegar a una conclusión (por una cuestión de tiempo y también de poca variación de la información). De esta manera, los cálculos de dosis deben realizarse de la forma más precisa posible (aunque todavía no existen límites como en Radioterapia) y además con cierta celeridad como para que no se desfase del resto de los puntos del protocolo. (2) Existen algunos radiofármacos que, por su vida media, no permiten obtener la suficiente cantidad de muestras en el tiempo como para que los cálculos sean los 2122 CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT suficientemente precisos, por lo que para alcanzar conclusiones de valor se deben hacer alguna clase de suposiciones. En otros casos, los servicios de Medicina Nuclear se encuentran saturados de trabajo, por lo que se hace muy difícil adquirir una secuencia temporal de imágenes lo suficientemente extensa como para que los resultados que de ella se infieran sean absolutamente precisos. 2223 CAPÍTULO 2 SIMULACIONES DE MONTE CARLO PARA EL CÁLCULO DOSIMÉTRICO EN PET Las simulaciones de Monte Carlo pueden describirse como métodos estadísticos que utilizan números aleatorios como base para la solución de problemas relacionados con procesos estocásticos. El concepto general de las simulaciones de Monte Carlo es la creación de un modelo, lo más parecido posible al sistema real y simular las interacciones con el sistema, basado en las conocimiento a priori de la probabilidad de ocurrencia durante el muestreo aleatorio de las funciones de densidad de probabilidad (PDF). Debido a la naturaleza estocástica de los procesos de radiación y detección de emisiones, el método de Monte Carlo es especialmente interesante para la física médica en áreas como la radioterapia, la protección radiológica y medicina nuclear [Andreo, 1991].(8) En la actualidad, las simulaciones en Monte Carlo son una herramienta esencial en la medicina nuclear para el desarrollo de nuevo marcadores moleculares, 2324 CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT la optimización de los protocolos de adquisición existentes y la implementación de nuevos algoritmos de reconstrucción y corrección de dispersión. A la fecha, han sido desarrollados muchos simuladores, entre algunos de los más importantes se encuentran PENELOPE, MNCP y GEANT4, los cuales son muy complejos y requieren demasiado tiempo para su correcto aprendizaje Números Aleatorios Los números aleatorios son de gran importancia para realizar simulaciones de sistemas físicos a través de modelos estadísticos, ya que son los iniciadores del proceso. Algoritmos matemáticos han sido desarrollados para generar dicha secuencia de números aleatorios, estos algoritmos son activados cada vez que comienza una simulación. Estos algoritmos deben cumplir: La generación de números aleatorios debe ser rápida Las secuencias generadas no pueden estar correlacionados entre sí El período de la secuencia debe ser largo Método de muestreo analógico o detallado. El método analógico de Monte Carlo intenta simular el desarrollo estadístico completo de diferentes procesos físicos, tales como el decaimiento radioactivo, al asumir un gran número de historias de partículas, N, e incluirlas en un determinada subconjunto, para poder considerar cada sub-conjunto dentro una distribución normal. Para un cálculo dado, la incertidumbre estimada es proporcional a la inversa de la raíz cuadrada del número de historias simuladas. La eficiencia "de un cálculo de Monte Carlo por lo tanto, puede definirse como (9) (2.1) 2425 CAPÍTULO 2 SIMULACIONES DE MONTE CARLO PARA CÁLCULO DOSIMÉTRICO EN PET Donde T es el tiempo de cálculo requerido para obtener una estimación de la varianza. Para un gran N, debe ser constante siempre que el método de cálculo que sigue siendo el mismo. El sistema de imágenes y/o distribución del código fuente pueden ser descrito en términos de funciones de densidad de probabilidad (pdf). Estos pdfs, complementado con cálculos adicionales, describir la evolución del sistema en general (espacio, energía, tiempo). El objetivo del método de Monte Carlo es simular el sistema de imágenes mediante un muestreo aleatorio de los pdfs y realizar los cálculos complementarios que sean necesarios para describir la evolución del sistema. Para obtener una variable estocástica que sigua una determinada pdf, pueden ser utilizados tres métodos de muestreo: el método de función de distribución, el método de rechazo y el método mixto [Zaidi, 1999; Ljungberg, 1998; Lázaro, 2003] Método de función de distribución (pdf) Este método se puede utilizarse, solo si la inversa de la función de distribución acumulada es fácil de obtener. Una función de distribución acumulada es construida desde la integral f(x) sobre el intervalo [a, x] de acuerdo a: (2.2) Entonces se distribuye uniformemente en [0, 1], el valor de muestreo x se puede obtener mediante la sustitución de en la ecuación (pdf) por un número aleatorio uniforme R, es decir Método de rechazo Cuando es demasiado complicado obtener la inversa de la función de distribución es recomendable utilizar la técnica de rechazo, que sigue los siguientes pasos: 2526 CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT Definir una función normalizada, donde es el valor máximo de. Muestrear dos números aleatorios distribuidos uniformemente y. Calcular x utilizando la ecuación. Si es igual o inferior a entonces x es aceptada como un valor de muestra, de lo contrario un nuevo valor de x es muestreado. La principal desventaja de este método es que se requiere la utilización de dos números al azar por ensayo y realizar muchas pruebas, las cuales son función del área bajo la curva de Método Mixto Cuando los dos métodos anteriores no son prácticos, el método mixto, que combina los dos opciones, se puede utilizar (9). Supongamos que la pdf se puede factorizar como: (2.3) Donde es una función invertible y es una función plana, que contiene la mayor parte de la complejidad matemática. El método consiste en los siguientes pasos: (i) Normalizar produciendo tal que (ii) Normalizar produciendo de tal manera que para x entre. (iii) Usar el método directo para seleccionar una x con como la pdf. (iv) Utilizar x y aplicar el método de rechazo, utilizando, es decir, elegiendo un número aleatorio R, si, se acepta x, en caso contrario se vuelve al paso (iii) Muestreo no analógico o condensado Las simulaciones de MC son poco prácticas, computacionalmente hablando, ya que se requiere un gran tiempo de cálculo para obtener resultados estadísticamente 2627 CAPÍTULO 2 SIMULACIONES DE MONTE CARLO PARA CÁLCULO DOSIMÉTRICO EN PET relevantes. Para solventar este inconveniente se recurre a la denominada simulación condensada, cuyo fundamento se encuentra en las teorías de dispersión múltiple. La idea consiste, muy brevemente, en simular el efecto global de un número elevado de interacciones mediante un único suceso artificial, para esto se recurren a técnicas que manipulen el sistema de muestro tales como: registrar solo regiones que tenga una gran contribución a la simulación, ajustar los límites del sistema o realizar ponderación. Las técnicas de reducción de muestro pueden dividirse en tres categorías: aquellos que solo afectan el transporte del fotón, las que solo afectan el transporte electrones y otros métodos más generales. Es importante señalar que debido a la naturaleza Poisson de la distribución de la actividad en la medicina nuclear, estas aproximaciones de reducción de la varianza pueden dar lugar a desviaciones estadísticamente relevante por lo que los resultados del método podrán llegar a ser deficientes. (9) Métodos específicos para fotones. Interacción Forzada. En una simulación de Monte Carlo analógica, los fotones son seguidos a través del objeto hasta que bien escapen del objeto, sean absorbidos o su energía caiga por debajo de un umbral seleccionado. La función de probabilidad de una interacción del fotón está dada por: (10) (2.4) La probabilidad de que un fotón viaje una distancia d o menor, se da por: (2.5) Para determinar la longitud del paso de muestreo, el número aleatorio uniforme R se sustituye por y se resuelve el problema para d: (2.6) 2728 Entonces la distancia máxima, CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT que el fotón viaja antes de la interacción es infinita y el número de caminos libres medios a través de la geometría en cualquier situación práctica es finito, por lo que hay una gran probabilidad de que los fotones dejen la geometría de interés sin interactuar. Para aumentar la precisión estadística en el cálculo de la energía impartida, los fotones son forzados a interactuar mediante la asignación de una distancia finita simular. (10), por ejemplo, el grosor del detector que se desea Transformación exponencial, ruleta rusa y la división de las partículas. La transformación exponencial es una técnica de reducción de varianza utilizado para sesgar el procedimiento de muestreo y para dar una mayor interacción en las regiones de interés, mejorando la eficiencia del cálculo en estas regiones. Métodos como la ruleta rusa y la división se utilizan a menudo junto con la transformación exponencial, aunque siguen siendo eficaces cuando se usan de forma independiente. En la ruleta rusa, un número al azar se selecciona y se compara con un umbral. Si el número al azar resulta ser menor que el número de partículas, se le permite sobrevivir, pero su peso debe modificarse. En la división de las partículas, una partícula procedente de una región de interés se puede dividir en N partículas, cada una con una nueva ponderación Métodos específicos para electrones Rango de rechazo del electrón. Una diferencia fundamental entre el transporte de fotones y el de electrones en el código de simulación, es que los fotones viajan distancias relativamente largas antes de interactuar, mientras que el camino de los electrones es interrumpido no sólo por los límites geográficos, sino también por los "pasos" de dispersión múltiple. El rango de rechazo de los electrones, significa que los electrones con rango menor que la distancia más cercana a la frontera o a la región de interés en la simulación se rechaza para ahorrar tiempo de computación. Reducción de parámetros en los de paso de electrones. Este algoritmo permite el uso de pequeños pasos de electrones en las proximidades de las interfaces y fronteras y pasos grandes en otros lugares.(11) Sus componentes son: un algoritmo de 2829 CAPÍTULO 2 SIMULACIONES DE MONTE CARLO PARA CÁLCULO DOSIMÉTRICO EN PET corrección de longitud de paso que se basa en la teoría de dispersión múltiple de Moliere, que tenga en cuenta las diferencias entre la longitud del camino recto y la longitud total de camino curvo para cada paso de electrones, un algoritmo de correlación lateral que tenga en cuenta transportes laterales, y un algoritmo de cruce fronterizo, que asegure que los electrones son transportados con precisión en las cercanías de interfaces Métodos generales Muestro Correlacionado. La técnica de muestreo correlacionado puede ser utilizada en el transporte de los fotones y electrones. Es especialmente eficaz para el cálculo de cocientes o diferencias entre dos cantidades que son casi iguales. La idea básica es que las simulaciones de las geometrías de interés se mantengan en estrecha correlación, de modo que la mayor parte de las fluctuaciones estadísticas, se cancelará en los cocientes y las diferencias, ya que las verdaderas diferencias entre dos geometrías se reflejan mejor a través de estas operaciones. Las incertidumbres de cálculo de los cocientes y las diferencias obtenidas con muestreo de correlación son, en general, menores que los obtenidos en las simulaciones no correlacionadas.(10) Uso de la simetría geometría. El uso de algunas de la simetría inherente de la geometría puede incrementar considerable la eficiencia. Si la configuración de la fuente y blanco está contenida en geometrías que simulen cilindros planas, esféricas o cónicas, el uso de simetrías es más evidente. Otros usos de la simetría son menos obvios, pero el ahorro en tiempo de cálculo es sumamente valioso Transporte de fotones Para los problemas de transporte de radiación, el modelo computacional incluye la geometría y las especificaciones de materiales. Cada código de computadora contiene una base de datos de las cantidades obtenidas experimentalmente, que se conoce como secciones eficaces, que determinan la probabilidad de que una partícula interactué con el medio a través del cual se transporta. Cada sección se define por el tipo y la energía de la partícula incidente y el tipo de interacción que sufre. Estas secciones parciales se suman para formar una sección eficaz total. Los datos de la sección eficaz de cada tipo 2930 CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT de interacción de intereses deben ser suministrados para cada material presente. El modelo también se compone de algoritmos utilizados para calcular el resultado de las interacciones (cambios en la energía de las partículas, dirección, etc.) sobre la base de los principios físicos que describen la interacción de la radiación con la materia y los datos de secciones eficaces proporcionados. Por lo tanto, es sumamente importante la utilización de un modelo de transporte preciso. (10) Cuando un fotón (que tiene una energía inferior a 1 MeV para los isótopos de interés en la medicina nuclear) pasa a través de la materia, cualquiera de los tres procesos de interacción: efecto fotoeléctrico, dispersión coherente y dispersión Compton, puede ocurrir. La probabilidad de que un fotón de una dada energía E sometidos a la absorción o dispersión al atravesar una capa de material Z, puede expresarse cuantitativamente en términos de un coeficiente de atenuación lineal µ (cm -1 ), que depende de la densidad del material (g cm -3 ):. (10) (1.7) La sección eficaz total del efecto fotoeléctrico para una dada energía, E, es calculada utilizando las diferentes bibliotecas, por medio de interpolación. El fotón incidente es absorbido y un electrón se emite en la misma dirección del fotón incidente. La energía cinética de los electrones se define como la diferencia entre la energía del fotón original y la energía de enlace del electrón. La partícula emitida del subnivel, es nuevamente muestreado aleatoriamente, utilizando para ello los datos de la sección eficaz de todos los subniveles. Cuando se simula la dispersión Compton de un fotón desde un electrón atómico a otro lado, una sección eficaz empírica se utiliza y el muestro del estado final se realiza por una variante del método de MC mixto. La sección eficaz diferencial por átomo según la ecuación de Klein-Nishina, está definida como (8)(12): (1.8) Donde: 3031 CAPÍTULO 2 SIMULACIONES DE MONTE CARLO PARA CÁLCULO DOSIMÉTRICO EN PET = radio clásico del electrón; = masa del electrón; = energía del fotón incidente; =energía del fotón dispersado;. Asumiendo una colisión elástica, el ángulo de dispersión formula Compton es definido por la (1.9) dado por: El valor de corresponde a la mínima energía del fotón (retro-dispersión) es (1.10) Entonces Se puede afirmar que: = (1.11) Donde 3132 CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT (1.12) y son funciones de densidad de probabilidad definidas en el intervalo, y se establece como: (1.13) Siendo la función de rechazo. Dado un conjunto de tres números aleatorios r, r, r uniformemente distribuidos en un intervalo, el procedimiento de muestreo para es el siguiente: I. Decidir si se desea una muestra de o : Si selecciona en caso contrario selecciona ; II. Muestrear de la distribución de o : Para : Para : III. Se Calcula donde ; IV. Se prueba la función de rechazo: si acepto. En caso contrario regreso a (I) Transporte de electrones y positrones La simulación del transporte de electrones y positrones es mucho más difícil de realizar que del transporte de fotones. La razón principal es que la pérdida promedio de energía de un electrón en una sola interacción es muy pequeña (del orden de unas decenas de ev) y, por tanto, electrones de alta energía sufren un gran número de colisiones antes de termalizarse. Por esta razón la mayoría de los códigos de Monte Carlo recurren a la teoría de múltiples dispersiones en el caso de electrones y positrones 3233 CAPÍTULO 2 SIMULACIONES DE MONTE CARLO PARA CÁLCULO DOSIMÉTRICO EN PET de alta energía, que permiten la simulación del efecto global de un gran número de eventos en un segmento de trayectoria de una longitud dada. (13) Cabe señalar que la teoría de dispersión múltiple utilizada en las simulaciones condensadas son sólo aproximaciones [Berger y Wang 1988, Kawrakow y Bielajew 1998] y pueden dar lugar a distorsiones sistemáticas de los resultados cuando el número de interacciones experimentadas por las partículas transportadas es pequeño, como ocurre, por ejemplo, en la transmisión a través de capas delgadas. (13) Las interacciones dominantes de electrones y positrones con los átomos son la dispersión elástica, colisiones inelásticas y la emisión de radiación de frenado, Figura 2.1. Los positrones también se someten a la aniquilación, ya sea en vuelo o en reposo. Figura 2.1. Interacciones básicas de los positrones y electrones con la materia.(14) Colisiones elásticas e inelásticas Los mecanismos de pérdida de energía dominante para los electrones y positrones con energías intermedias y bajas son las colisiones inelásticas, es decir, las interacciones que producen excitaciones y ionizaciones electrónicas en el medio. (14) 3334 CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT En sentido estricto, todos los diferentes tipos de interacción de electrones deben ser considerados en detalle en una simulación de Monte Carlo, tal como se hace con las interacciones de fotones. Una de las técnicas utilizadas para este objetivo, es la llamada simulación microscópica (15), la cual usa electrones de baja energía. El modelo más simple utilizado por este método se basa en la simulación de todos los eventos individuales de dispersión, calculando la longitud de paso del electrón entre colisiones consecutivas con el camino libre medio elástico. Las pérdidas de energía se determina a partir de la teoría de Bethe del stopping power y la pérdida de energía dispersada se representa por aproximaciones. Estos modelos se han mejorado, tomando en cuenta las colisiones inelásticas. Cuando la energía del electrón aumenta por encima de unas pocas decenas de kev, el gran número de interacciones que se producen durante la desaceleración de electrones hace que el método antes mencionado, sea poco realista para simular todas las interacciones físicas. Por lo cual se recurre al uso de la teoría de las colisiones múltiples, en donde se tiene en cuenta una serie de "fotografías" de las historias de los electrones, tomadas en cierto tiempo o intervalos de energía, las interacciones individuales de diferentes tipos se consideran entre dos fotografías consecutivas. Este régimen de agrupación constituye el macroscópico (15) o técnica de condensación de historias. En los últimos años, los algoritmos de transporte de electrones más avanzados, combinan los modelos macroscópicos y microscópicos, donde se utilizan estos últimos en las situaciones límite que requieren sofisticadas el cruce de algoritmos.(16) Bremsstrahlung Como resultado de la aceleración causada por el campo electrostático de los átomos, los electrones rápidos (o positrones) emiten radiación de frenado. En cada evento de bremsstrahlung, un electrón con energía cinética E genera un fotón de energía W, que toma valores en el intervalo de 0 a E. El proceso se describe con un sección eficaz diferencial (DCS) atómica, la diferencia en la pérdida de energía W, la dirección final del proyectil y la dirección del fotón emitido (Koch y Motz, 1959; Tsai, 1974). 3435 CAPÍTULO 2 SIMULACIONES DE MONTE CARLO PARA CÁLCULO DOSIMÉTRICO EN PET La práctica habitual en la simulación de Monte Carlo es muestrear la pérdida de energía de la distribución de una variable obtenida mediante la integración de la DCS en las otras variables. Esto permite la generación de W con facilidad, pero la información sobre las distribuciones angulares se pierde por completo y tiene que ser recuperado a partir de aproximaciones adecuadas.(14) Aniquilación positrónica Los positrones pueden aniquilarse con los electrones en el medio, emitiendo dos fotones, para garantizar asi la conservación de la energía y el momento. En la mayoría de los códigos de Monte Carlo, la aniquilación de positrones es descrita por el supuesto de que los electrones son libres y en están reposo, por lo tanto sin tener en cuenta efectos de ligadura de los electrones, y teniendo en cuenta sólo la emisión de dos fotones, puede considerarse que la aniquilación de un positrón con energía cinética E a través de la emisión de dos fotones con energías E + y E -, que sumadas son equivalente E+2m e c 2, entonces completamente caracterizado por la cantidad: (13) (1.14) Donde el subíndice - denota el fotón de menor energía. (13) Implementación Las simulaciones de Monte Carlo del transporte de electrones, deben ser consideradas en situaciones cuando los rangos de las partículas sean comparables con las dimensiones o resolución espacial del problema a resolver o la contribución de la radiación de frenado juegue un papel significativo en la cantidad de energía depositada localmente o transportada. Debido a las bajas energías relacionadas con la medicina nuclear, se asume que los electrones secundarios depositan toda su energía en el punto de la interacción y la producción de bremsstrahlung es insignificante. Por lo tanto, el transporte de electrones no ha recibido particular atención en las aplicaciones de 3536 CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT medicina nuclear del método de Monte Carlo. Sin embargo, la simulación del transporte de electrones, es de vital importancia para los cálculos de dosimetría. (16) En el caso de la dosimetría, el interés se basa en el cálculo de la razón del stopping power, medio-detector, donde el uso del método de Monte Carlo para obtener espectros de electrones, junto con el análisis de las teorías de cavidades ha permitido la determinación de las relaciones de stopping power que incluyen todas las generaciones de los electrones y los fotones producidos durante la simulación del transporte de electrones o fotones incidentes Códigos de Monte Carlo utilizados para aplicaciones de Medicina Nuclear Existe gran cantidad de software que permiten realizar simulaciones en Monte Carlo. En la Tabla 2.1 son enumerados alguno de los más utilizados junto con una breve descripción de sus características clave. La mayoría de los paquetes mencionados a continuación se ejecutan en distintas plataformas y sistemas operativos, y están disponibles de forma gratuita CÓDIGO MC DESCRIPCIÓN GENERAL PENELOPE EGS4 MCNP GEANT ETRAN Cuenta con un modelo detallado para el tratamiento de las secciones transversales a bajas energías. Diseñados para el transporte de electrones y fotones en el rango de 1 kev a varios TeV. Los usuarios de bienes y servicios ambientales creó el código de viga para modelar las fuentes de radiación y DOSXYZ de cálculo de dosis Originalmente desarrollado para la física de reactores y el transporte de neutrones. Posteriormente, se incluyó el transporte de electrones (desde la versión 4 a) para evaluar dosis en geometrías complejas (superficies de hasta 4 grados, estructuras repetidas) Originalmente desarrollado para la física de partículas de altas energías, permite el transporte de partículas múltiples. (Modelos de bajo consumo, nivel y viejo Penélope). La simulación para el cálculo de la dosimetría interna no está específicamente incluida y requiere una gran cantidad de programación por el usuario en C++. Berger describe la técnica de historia condensada de transporte de electrones, que es la base de todos los códigos actuales de transporte Tabla 2.1. Las principales características de los códigos de MC de dominio público utilizados en las aplicaciones de medicina nuclear.(10) 3637 CAPÍTULO 2 SIMULACIONES DE MONTE CARLO PARA CÁLCULO DOSIMÉTRICO EN PET. 3738 CAPÍTULO 3 PLATAFORMA DE SIMULACIÓN GATE GATE, es un simulador de Monte Carlo que utiliza las librerías de GEANT4 y que fue adaptada para su fácil implementación en el campo de la medicina nuclear. Entre algunas de las ventajas que presenta en las simulaciones de PET y SPET sobresale la capacidad de modelar fenómenos dependientes del tiempo tales como detectores y fenómenos de decaimiento de la fuente, lo que permite la obtención de simulaciones en tiempo real. Combinando las ventajas que ofrecidas por GEANT4 (procesos físicos validados, geometrías bien definidas y poderosas herramientas de visualización) y el uso de nuevas herramientas desarrolladas para el análisis de fantomas y la representación de imágenes de pacientes en condiciones reales (fuentes y fantomas voxelizados), GATE permite al usuario lograr simulaciones precisas y reales de equipos, como el PET.39 CAPÍTULO 3 PLATAFORMA DE SIMULACIÓN: GATE Adicionalmente la capacidad que posee para sincronizar todos los componentes dependientes al tiemp, en función de la evolución de la actividad de la fuente, subdividiendo la adquisición en pasos de tiempo, en donde se almacena la cinemática de decaimiento del radioisótopo, permite modelar adecuadamente procesos tales como tasa de conteo, tiempo muerto y coincidencias al azar. GATE se ha desarrollado desde 2001 por un consorcio de institutos internacionales/laboratorios sin ningún tipo de apoyo financiero específico, y ha sido publicado en Un lanzamiento público del GATE licenciado bajo la GNU Licencia Pública General Menor (LGPL) se puede descargar en la dirección electrónica: El GATE, es la estructura de simulación que constituye la base para los estudios llevados a cabo a través de esta tesis. La investigación presentada en este trabajo se ha desarrollado en la versión de GATE 6.1.0, publicada el primero de marzo del 2011, que incluye algunas herramientas específicas para la radioterapia (RT) de aplicaciones [GATE, 2010b]. La documentación en la que se basa este capítulo en la documentación de la puerta, los comentarios del código fuente y el código de GATE en sí mismo [Jan et al.,2008; GATE, 2010a; Jan et al., 2004b] Características Generales Arquitectura del software Una simulación es general dividida en 7 pasos, estos son: Visualización y nivel de detalle. Geometría del escáner y el fantoma. Selección de los procesos físico. Determinación del tipo fuente. Elección de parámetros de digitalización. Obtención de datos. Experimento. 3940 CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT Los cuatro primeros pasos son los encargados de definir al sistema, es decir, son los comandos de inicialización de la simulación, los tres últimos pasos terminan de dar la forma a la simulación, pero no pueden modificar los parámetros de inicialización de la simulación Visualización Hay varias herramientas disponibles para la visualización online (OpenGL, VRML, DAWN) que pueden ser activadas en función de las opciones de visualización seleccionada por el usuario. La visualización en línea es una herramienta útil en el desarrollo de nuevas geometrías, permitiendo comprobar visualmente la geometría del escáner (posiciones, las superposiciones de volúmenes físicos, etc.). Una vez que la geometría se controla y se quiere ejecutar una simulación completa, se recomienda desactivar la visualización en línea con el fin de no sobrecargar la CPU Detalle Para cada módulo de simulación, se puede establecer un nivel de detalle entre 0 y 2. Cuanto más alto sea el nivel de detalle, mayor será el nivel de información devuelta por GATE. Por defecto, el nivel de detalle se establece en 0, pero si se quiere seguir con detalle cada paso de la simulación, se puede establecer un valor más altos Geometría Todos los volúmenes están unidos entre sí siguiendo una estructura de árbol, donde cada rama representa un volumen. La base del árbol está representada por el volumen madre o mundo (The world) figura 3.1, que establece el marco experimental de la simulación. (17) El mundo es una caja centrada en el origen y de tamaño variable, que define los límites de la simulación, ya que el seguimiento a una partícula termina, cuando está escapa de su volumen. Esté, encierra varios sub-volúmenes que son llamados como volumen hija (daughter). Cada volumen se caracteriza por su forma, tamaño, posición y material. 4041 CAPÍTULO 3 PLATAFORMA DE SIMULACIÓN: GATE Figura 3.1. Volumen de madre (the world).(17) El sistema El sistema, proporciona un modelo de geometría predefinido para simular diferentes tipos de equipos. GATE dispone por defecto de plantillas para cabezales de SPECT, cilindros de PET, ECAT, CPET y OPET, que permiten simular la mayoría de los dispositivos de imágenes tomográficos y de PET existentes. Cada sistema tiene un número determinado de niveles con una organización jerárquica específica (tipo árbol) y un formato de salida de datos específicos. A Cada nivel se le debe asignar una geometría de volumen físico, que estará dentro de un volumen superior; GATE ha fijado el número de niveles propios de un sistema, las diferentes especificaciones son mostradas en la Tabla 3.1. cumplir:(18) Para realizar una conexión adecuada entre el sistema y su geometría, se debe Primero es necesario definir la estructura geométrica, con las restricciones pertinentes del sistema. Luego, el sistema geométrico debe ser introducido o atado, en el proceso de simulación con los comandos propios del GATE. Tabla 3.2. Por último, debe ser definido el formato de salida, para el análisis de datos. 4142 CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT SISTEMA COMPONENTES FORMA SALIDAS DISPONIBLE level1 level2 Escáner level3 Geometría no definida Salida básica: Ascii o ROOT, conteo de coincidencia level4 solo para escáneres PET level5 Escáner CT block caja Imagen binaria CT pixel caja Salida básica: Ascii o ROOT rsector caja Cilindro del module caja Salida básica: Ascii, ROOT and Raw PET submode caja Especifica: LMF crystal caja layer caja CPET crystal cilindro Salida básica: Ascii, ROOT and Raw Cabezal de crystal Geometría no definida Salida básica: Ascii, ROOT and Raw SPECT pixel Especifica: Set de proyecciones o interfile ecat block caja Salida básica: Ascii, ROOT and Raw crystal caja Especifica: Sinograma o Ecat7 ecataccel block caja Salida básica: Ascii, ROOT and Raw crystal caja Especifica: Sinograma o Ecat7 rsector caja module caja Salida básica: Ascii, ROOT and Raw OPET submode caja Especifica: LMF crystal layer caja cuña Tabla 3.1. Descripción de los diferentes sistemas disponibles en GATE. La segunda columna se muestra los nombres con los cuales los volúmenes hijas son llamadas en la simulación.(17) 4243 CAPÍTULO 3 PLATAFORMA DE SIMULACIÓN: GATE Sistema Comandos para atar Profundidad level1 1 level2 2 Escáner level3 3 level4 4 level5 5 Escáner CT block 1 pixel 2 rsector 1 module 2 Cilindro del PET submode 3 crystal 4 layer 5 CPET crystal 1 Cabezal de SPECT crystal 1 pixel 2 ecat Block 1 crystal 2 ecataccel Block 1 crystal 2 rsector 1 module 2 OPET submode 3 crystal 4 layer 5 Tabla 3.2. Palabras claves que corresponden a los componentes del sistema que pueden ser utilizados para atar la geometría a la simulación.(17) 4344 CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT Detectores sensibles Los detectores sensibles, se utilizan para almacenar información sobre las interacciones de una partícula con la materia, hits, usando información de los pasos a lo largo de la trayectoria de una partícula. Un hit es una foto de una interacción física de una trayectoria de una partícula en una región sensible de un detector, dicho concepto es ilustrado en la figura3.2. GATE, solo registra y almacena la información relacionada a los hits de los volúmenes que están conectados a un detector sensible. Toda la información sobre las interacciones que ocurren en los volúmenes no sensibles se pierde. Dos detectores sensibles se definen en GATE: El crystalsd permite registrar información sobre las interacciones dentro de los volúmenes pertenecientes a un escáner, por ejemplo, los cristales o colimadores. El phantomsd, puede ser utilizado para registrar información sobre las interacciones Compton y Rayleigh que tienen lugar en el campo de visión del escáner. Esta información se utiliza para estimar si un fotón que llega a un detector es directo o es producto de dispersión Compton. Así, en PET, el phantomsd es actualmente la única manera de discriminar los fotones dispersos de los directos. Figura 3.2. Interacción de una partícula con un detector sensible (17) 4445 CAPÍTULO 3 PLATAFORMA DE SIMULACIÓN: GATE Procesos físicos Todas las interacciones electromagnéticas utilizadas en GATE se derivan de GEANT4. [Agostinelli et al., 2003; Allison et al., 2006]. Existen tres modelos disponibles para simular estos procesos: Procesos de energía estándar-std (Rango efectivo 1 kev-100 TeV). Procesos de baja energía - LowE (Rango efectivo 250 ev-100 GeV). Procesos basados en Penélope (Rango efectivo 250 ev-1 GeV) Todos estos modelos pueden manipular electrones, positrones, rayos gammas, rayos X, fotones ópticos, muones, hadrones e iones (17). Los diferentes procesos físicos disponibles para los diferentes modelos son: Efecto Fotoeléctrico Efecto Compton. Dispersión Rayleight. Bremsstrahlung. Producción de Pares. Ionización (electrón - positrón, hadrones-iones, iones). Dispersiones simples y múltiples. Neutrones. Aniquilación positrónica Producción de cortes GATE posee funciones de corte que limitan la producción de partículas secundarias. Los cortes pueden ser de dos tipos, por rango o por energía. Los cortes por rango, definen una distancia mínima, en la cual se permite producción de partículas. Este tipo de corte no es aplicable a los límites de volúmenes cerrados, ni para algunas partículas. Los cortes por energía, definen una energía mínima, en la cual se permite la producción de partículas. 4546 CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT Fuentes Para definir una fuente en GATE es necesario definir el tipo de partícula a utilizar (radionúclido, positrones, gamma) y sus características (la posición, la distribución angular de irradiación, el espectro energético y la distribución de actividad). En cada nuevo evento la fuente directora decide al azar el decaimiento a seguir y genera para esto, una o más partículas primarias. Los datos de la voxelización del paciente o fantoma pueden ser utilizados como fuentes para reproducir adquisiciones realistas: los datos del PET se convierten en los niveles de actividad, y los datos tomográficos se convierten en los mapas de atenuación dentro de la definición de los materiales mediante un traductor análogo, del cual se hablara más adelante.(8) Voxelización de fuentes y fantomas GATE ofrece dos métodos para describir múltiples copias de un volumen dentro del volumen madre: creación de replicas y parametrización de volúmenes. La parametrización de volúmenes ofrece grandes ventajas sobre el método de creación de replicas, al permitir la variación del tamaño, forma y material del voxel, por lo que es el método recomendado por GATE. Para la parametrización de los volúmenes es necesario introducir las imágenes tomografías y de PET en el programa, la primera definirá la anatomía del paciente y la segunda la distribución de actividad. La lectura de la imagen, requiere sus datos sean enteros sin signo de 16 bits y que se encuentre en formato ASCII o interfile. En el primer modo es necesario declarar (mediante scripts) que el fantoma voxelizado se introducirá a partir de datos tipo ASCII, antes de insertar el archivo ASCII que contiene todos los números necesarios para todos los píxeles (usando la opción imagereader). Además es necesario especificar los materiales de atenuación usando un traductor tabulado, que interpretar los materiales de la tabla de datos. El fantoma se convierte entonces en definiciones de materiales. 4647 CAPÍTULO 3 PLATAFORMA DE SIMULACIÓN: GATE La entrada Interfile es más fácil de usar. En este caso es necesario usar el comando InterfileReader para leer el archivo de datos Interfile. De la misma manera que el formato ASCII, con el lector Interfile los usuarios tienen que especificar los materiales del fantoma, en este caso mediante el uso de un traductor de rango. En el traductor de rango son específicos los materiales a partir de las bibliotecas de Geant4 y el rango de voxel, en el cual se encuentran las estructura anatómicas a definir. La base de datos de material se puede actualizar con los materiales del fantoma o con nuevos materiales relacionados con la geometría del escáner en uso (GateMaterials.db, es el archivo que contiene todos los parámetros de los materiales requeridos por Geant4 para calcular la interacción de la sección eficaz). (8) Administración de tiempo y movimiento Una de las características más distintivas de GATE es su capacidad para manejar fenómenos dependientes del tiempo [Strul et al., 2003; Santin et al., 2003]. La sincronización de la cinética de la fuente con los movimientos de la geometría, permite simular de una forma realista las condiciones de adquisición, incluyendo los movimientos del paciente (respiración y movimientos cardiacos), la rotación del escáner y los cambios de la distribución de actividad a lo largo del tiempo.(8) Tratar con el tiempo en GATE incluye: Describir los movimientos asociados con los volúmenes físicos que son descritos por el detector y el fantoma. Descripción de la fuente radioactiva. Especificación de los tiempos de comienzo y finalización de adquisición, que son equivalentes a los tiempos de inicio y finalización del experimento real. Todos los movimientos están sincronizados con la evolución de las actividades de la fuente mediante la subdivisión del tiempo de adquisición en pasos de tiempo más pequeños (típicamente del orden de minutos u horas). Dentro de cada uno de estos pasos la fuente comienza a decaer exponencialmente desde un paso de tiempo al próximo 4748 Digitalización y parámetros de lectura CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT El sistema de digitalización, es uno de los mecanismos que utiliza GATE para asegurar al usuario que la simulación es lo más realista posible. Esté se encarga de simular el comportamiento de los detectores del escáner y la cadena de procesamiento de señales, que es utilizada en la clasificación de los diferentes tipos de coincidencias (reales, dispersadas y múltiples). GATE hace uso de Geant4 para evaluar el comportamiento de las partículas generadas que son transportadas a través de los diferentes materiales del fantoma. Para esto, imita las interacciones físicas entre las partículas y la materia. La información generada durante este proceso es utilizado por GATE para simular los pulsos del detector, que corresponden a los datos observados. El digitalizador representa la serie de pasos y los filtros que componen este proceso. Un típico flujo de datos para un evento se compone de la siguiente manera: a. Una partícula es generada, con sus respectivos parámetros, tales como tipo inicial, tiempo, movimiento y energía. b. Una trayectoria elemental es elegida aleatoriamente, la cual correspondiente a la trayectoria de una partícula entre interacciones discretas (es decir, fotoeléctrico, compton, etc.). Durante este paso los cambios de energía y momento de la partícula son calculados. La longitud del paso depende de la naturaleza de la interacción, el tipo de la partícula y el material, etc. c. Si un paso ocurre dentro de un volumen correspondiente a un detector sensible, la información de la interacción partícula materia es almacenada. Este conjunto de datos se conoce como un hit. d. Los pasos a y b se repiten hasta que la energía de la partícula sea tan baja como el valor predeterminado o hasta que la posición de la partícula se encuentra fuera de los límites fijados. e. La cantidad de energía depositada en un cristal se filtra por el módulo digitalizador. La salida del digitalizador corresponde a la señal después de haber sido procesada por el Front End Electronics (FEE). El FEE transforma la energía de un hit en el valor digital final, a través de un proceso llamado 4849 CAPÍTULO 3 PLATAFORMA DE SIMULACIÓN: GATE digitalización, que se realiza por la parte digitalizadora de la arquitectura de GATE. El valor final obtenido tras filtrar por un conjunto de estos módulos se llama una single. Cada valor transitorio, entre dos módulos, se llama pulso. Este proceso se repite para cada evento en la simulación con el fin de producir una o más series de singles. Una vez que esta lista se crea, una segunda etapa de procesamiento se pueden insertar para ordenar la lista de singles en coincidencias. Para ello, el algoritmo busca en esta lista, los singles que son detectan dentro de un intervalo de tiempo dado (llamados eventos de coincidentes). Finalmente, los datos de las coincidencias pueden ser filtradas en la salida para imitar cualquier posible pérdida de datos que pudiera producirse en el circuito de coincidencia lógica o durante el transporte de datos. (17) Para simular los valores digitales (pulsos) que resultan de la salida de la FEE, los métodos de muestreo de la señal deben ser especificados. Con el fin de hacer esto, una serie de módulos digitalizadores están disponibles y se describen a continuación. El papel del digitalizador es construir, a partir de la información de los hits, la física observable que comprende la energía, la posición y el tiempo de detección para cada partícula. Además, el digitalizador debe implementar la lógica necesaria para simularlas coincidencias en las simulaciones del PET. El uso típico del módulo digitalizador incluye las siguientes acciones: Simular la respuesta del detector, Simular el esquema de lectura y Simular la lógica de disparo. La digitalización consiste en una serie de procesadores de señal, figura 3.3. La salida en cada paso a lo largo de la serie se define como un pulso. En la figura 3.4 se ilustran las acciones tanto de modulo sumador como el de lectura. Después de sumar los hits y leer los pulsos, los hits se reagrupan en pulsos y entonces se suman los pulsos, los módulos restantes de la cadena de digitalizadores, transforman estos pulsos en fenómenos físicos observables por el escáner: singles. Estos singles 4950 CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT simulan de una forma realista los fenómenos físicos observables por el detector como respuesta de la interacción de una partícula con él. (17) Figura 3.3. El digitalizador está organizado como una cadena de módulos que comienzan con un hit y termina con un single, que representa los fenómenos físicos que pueden ser observables por el detector.(17) Figura 3.4. Acciones del modulo sumador y de lectura. El módulo de sumador transforma los hits en un pulso en cada volumen individual, a continuación, el módulo de lectura suma estos pulsos en uno solo pulso en el nivel de profundidad definido por el usuario. (17) Experimento Durante este último paso de la simulación, es fijado el tiempo de duración de la adquisición de la misma, definiendo el comienzo y el final del proceso. Este tiempo, en general puede ser subdividido en fracciones de tiempo o proyecciones, por medio de la ecuación: (3.1) 5051 CAPÍTULO 3 PLATAFORMA DE SIMULACIÓN: GATE Donde: número de proyecciones. Hora de inicio de la simulación. Hora de finalización del estudio. Duración del slice. Esta característica es muy útil en GATE, ya que la geometría se actualiza sólo entre dos intervalos de tiempo. Esto ofrece la posibilidad de tener en cuenta los movimientos de las fuentes o los detectores mediante la subdivisión de una proyección en intervalos de tiempo, que permanecen constantes para una misma geometría Cálculo de Dosis GATE calcula la dosis depositada por voxel a través de la expresión: (3.2) Donde el sub-índice i, corresponde a cada uno de los voxeles que conforman el fantoma o órgano de estudio, es la energía promedio por voxel y m es la masa de cada voxel. Para determinar la dosis colectada depositada en el fantoma, se debe relacionar el detector sensible con el fantoma a través del comando attachvoxelphantomsd y se añade el módulo de salida de la dosis:(17) /gate/anyname/attachvoxelphantomsd /gate/anyname/addoutput outputmodulename El archivo de salida es un archivo binario (formato de número es de 4 bytes float) que contiene la dosis en cgy y posee las mismas dimensiones que el fantoma. 5152 CAPÍTULO 4 PET/CT GE DISCOVERY STE Figura 4.1. Equipo PET/CT GE Discovery STE53 CAPÍTULO 4 PET/CT GE DISCOVERY STE 4.1- Especificaciones técnicas El PEC/CT Discovery STE es un escáner comercial que consta de 24 anillos detectores: cada uno hecho de 35 bloques detectores. Cada bloque detector está compuesto por 13,440 cristales de BGO con dimensiones de 4.7 x 6.3 x 30 mm 3. Esté tiene un campo de visión axial del PET 70 cm y un diámetro de 88.6 cm Modelo de simulación Las instrucciones del modelo de simulación (desde la geometría al protocolo de adquisición), para la plataforma GATE serán mostradas a continuación. En la tabla 4.1 se resume la geometría del sistema y las condiciones de adquisición. Especificaciones Número de imágenes planas 47 FOV axial (cm) 15,7 Sensibilidad 2D (cps/kbq) 2 Sensibilidad 3D (cps/kbq) 8,5 Fracción dispersada 2D (%) 19 Fracción dispersada 3D (%) 35 Ventana de coincidencias 2D (nseg) 11,05 Ventana de coincidencias 3D (nseg) 9,75 Máxima tasa de cuentas 2D 49 kbq/cc) 95 Máxima tasa de cuentas 3D 12 kbq/cc) 75 Máxima tasa de cuentas 3D clínica 6 kbq/cc) 57 Resolución transaxial, fuente lineal a 1 cm para 2D y 3D (mm) 5 Resolución transaxial, fuente lineal a 10 cm para 2D (mm) 5,7 Resolución transaxial, fuente lineal a 10 cm para 3D (mm) 5,6 Resolución transaxial, fuente lineal a 20 cm para 2D y 3D (mm) 6,3 Resolución axial, fuente lineal a 1 cm para 2D (mm) 4,8 Resolución axial, fuente lineal a 1 cm para 3D (mm) 5,6 Resolución axial, fuente lineal a 10 cm para 2D y 3D (mm) 5,9 Tabla 4.1. Especificaciones del PET/CT GE Discovery STE 5354 CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT Geometría Como plantilla para la simulación se utilizo el sistema del PET cilíndrico. La relación entre los componentes físicos del escáner y los niveles jerárquicos del sistema es descrito en la tabla 4.2. La geometría del escáner es mostrada en la figura 4.2. En el apéndice 1 se muestra el código utilizado para la definición de la geometría del escáner utilizado en este trabajo Componentes Físicos Sistema de PET Cilíndrico Cabezal rsector Modulo module Cristal crystal BGO layer0 Tabla 4.2. Relación entre los componentes físicos del escáner y los niveles jerárquicos del sistema Figura 4.2. Ilustraciones de la geometría del PET Discovery STE, dada por GATE Procesos físicos y protocolos de adquisición Para la descripción de la dispersión Compton y el efecto fotoeléctrico se uso el paquete de procesos electromagnético de energía estándar y el paquete de baja energía fue aplicado para la dispersión Rayleigh. Para acelerar la simulación, se establecieron rangos de corte equivalentes a 0.01 mm, tanto en el fantoma como en el cristal, por debajo de los cuales electrones, fotones o positrones no pueden generar partículas secundarias. 5455 CAPÍTULO 4 PET/CT GE DISCOVERY STE Toda la simulación fue hecha con todos los parámetros que normalmente se utilizan para una adquisición real. Por defecto GATE, utiliza una distribución Gaussina y una distribución exponencial para definir el ruido de fondo: la primera distribución, llamado energy_distrib, es una Gaussiana centrada en 450 kev con una desviación estándar de 1 kev, mientras que la segunda, dt_distrib, es una distribución exponencial con una potencia de 7,57 ms. En apéndice 2 se muestra el código utilizado. A la distribución de Gauss se le atribuye los niveles de energía del fondo, y a la exponencial todos los fenómenos que puedan ocurrir en el intervalo de tiempo entre dos eventos de fondo consecutivos. Un blurring de 25% del pico de máxima energía (511 kev) es aplicado a la energía depositada en los cristales BGO. Un fotón es detectado y procesado sólo si su energía se sitúa entre 350 y 650 kev. Igualmente fue aplicado un tiempo muerto de 250ns, tipo no paralysabel (solo guarda los datos que se están procesando, los demás son descartados) en los eventos individuales. La ventana de coincidencia fue definida en un tiempo de 9.75 ns y la ventana de retraso se definió con un offset de 500 ns Datos de salida y configuración de reconstrucción En GATE, hay varios tipos de formato de salida, que pueden ser activado de manera predeterminada. Para todos los sistemas están disponibles los archivos tipo ASCII y ROOT. El formato ROOT tiene los mismos datos que el formato ASCII, pero estos son comprimidos automáticamente, por lo que es recomendado en caso de alta estadística. Un ejemplo de la salida ROOT, es descrita a continuación. Durante la adquisición de la simulación todos los datos se almacenan en el archivo Fantoma.root: # D A T A O U T P U T # /gate/output/root/enable /gate/output/root/setfilename Fantoma /gate/output/root/setroothitflag 0 /gate/output/root/setrootsinglesflag 1 /gate/output/root/setrootcutsinglesflag 1 5556 CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT /gate/output/root/setrootcoincidencesflag 1 /gate/output/root/setrootdelayflag 0 /gate/output/root/setrootfinalcoincflag 1 /gate/output/root/setrootntupleflag 0 /gate/output/root/setsaverndmflag 0 El formato Root es un macro para el procesamiento de datos, desarrollado en el CERN, pensado en investigación sobre física de alta energía [ROOT, 2010; Brun & Rademakers, 1997]. El sistema ROOT proporciona un conjunto de marcos orientado a objetos, construida en C + +, con toda la funcionalidad necesaria para gestionar y analizar grandes cantidades de datos de una manera muy eficiente. (8) Experimento Como se explico en la sección , en el experimento se define el tiempo que dura la simulación, para nuestro caso trabajaremos con un tiempo de simulación de 60 seg, subdividido en intervalos de 10 seg. Un ejemplo de la línea de comando utilizada se muestra a continuación: /gate/application/settimeslice /gate/application/settimestart /gate/application/settimestop /gate/application/startdaq 10 s 0 s 60 s 5657 CAPÍTULO 5 SEGMENTACIÓN Y PROCESAMIENTO DE IMÁGENES La ventaja más importante que brindan las simulaciones es la posibilidad de tener un modelo realista de la anatomía del sujeto y las funciones fisiológicas a través de datos obtenidos de imágenes generadas en resonancia, tomografía o PET. Conceptualmente, el propósito de un fantoma físico o informático es representar a un órgano o región de interés, para poder modelar la biodistribución de un radio trazador en particular y la composición química del medio de dispersión, que absorbe y dispersa la radiación emitida de manera similar a los tejidos biológicos. En otras palabras, un fantoma es un modelo matemático diseñado para representar, con la mayor precisión posible, un órgano o tejido del cuerpo, un sistema de órganos, o el cuerpo entero. (8)58 CAPÍTULO 5 SEGMENTACIÓN Y PROCESAMIENTO DE IMÁGENES 5.1- Elección de órganos de estudio. La correcta elección de los órganos de estudio es fundamental cuando se habla de dosimetría interna, para obtener resultados relevantes, hablando desde el ámbito de la protección radiológica. Los órganos de estudio seleccionado fueron: Los órganos de mayor concentración de actividad y de mayor eliminación biológica. Los órganos mas radiosensibles Órganos de mayor concentración de actividad y de mayor eliminación biológica Los órganos de mayor captación y eliminación biológica fueron tomados de las imágenes de PET de varios pacientes, figura 5.1, bajo la supervisión de los diferentes médicos del servicio. Los órganos elegidos son nombrados en la Tabla 5.1. Figura 5.1. Selección de órganos fuentes a partir del estudio PET. 5859 CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT TIPO ÓRGANO Captación Cerebro (Órgano más captante de FDG) Hígado Intestino (colón derecho, zona del ciego) Corazón Timo (principalmente en los jóvenes) Eliminación Vejiga Riñones Tabla 5.1. Clasificación de órganos fuentes, según su tipo: captantes o de eliminación Órganos más radiosensibles El termino radio-sensibilidad se refiere a la sensibilidad relativa de células, tejidos, órganos, organismos o cualquier otra sustancia viviente ante los efectos de la radiación. La radiosensibilidad de los tejidos difiere mucho entre ellos. Existen tejidos con alta radiosensibilidad, tales como: Médula ósea y timo. Epidermis, mucosa orofaríngea. Epitelio intestinal, epitelio del cristalino. Gónadas. Tiroides. Existen también tejidos con radiosensibilidad intermedia, tales como: Riñón. Pulmón. Finalmente, existen tejidos de baja radiosensibilidad, tales como: 5960 CAPÍTULO 5 SEGMENTACIÓN Y PROCESAMIENTO DE IMÁGENES Músculo esquelético y cardíaco. Hueso, tejido conectivo blando. Sistema nervioso. En el ICRP 103, dicha clasificación se hace a través del factor W T. El factor W T, es el factor por el que se pondera la dosis equivalente en un órgano o tejido T para representar la contribución relativa de ese órgano o tejido, al detrimento total en la salud, como respuesta a una exposición total del cuerpo (ICRP 1991b). La ponderación se efectúa de forma que:(19) (5.1) En la Tabla 5.2 se proveen los órganos y tejidos para los cuales los valores del W T están especificados Tejido W T ΣWT Medula ósea, colon, pulmón, estómago, mama, resto de los tejidos * 0,12 0,72 Gónadas 0,08 0,08 Vejiga, esófago, hígado, tiroides 0,04 0,16 Superficie del hueso, cerebro, glándulas salivales, piel 0,01 0,04 Total 1 Tabla 5.2. Factores de ponderación de los tejidos recomendados. * Resto de los Tejidos: Adrenales, región extra torácica (ET), vesícula, corazón, riñones, nódulos linfáticos, músculo, mucosa oral, páncreas, próstata, intestino delgado, bazo, timo, útero/cérvix. (19) 5.2- Descripción de la fuente. La mayoría de los radionucleidos emisores de positrones presentan vidas medias extremadamente cortas, por lo cual, es necesario que su producción sea en un sitio cercano al sitio de utilización. Todos son producidos por ciclotrón en vez de reactor, dado que los emisores de positrones son deficientes de neutrones. Los más usados son el Flúor-18, el Carbono-11, el Oxígeno-15 y el Nitrógeno-13, todos los cuales son capaces de marcar moléculas orgánicas fisiológicamente utilizadas por el metabolismo celular. 6061 CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT Para el desarrollo del presente trabajo fue utilizada la 2-[ 18 F]-flúor-2-desoxi-Dglucosa o [ 18 F]-FDG. La 18F-fludesoxiglucosa (FDG) es un radiofármaco constituido por un análogo de la glucosa 2-[18F]-fluoro-2-desoxi-D-glucosa unido al isótopo radiactivo flúor-18. Su administración se realiza de forma intravenosa. Su fórmula molecular es C 6 H FO 5 y su estructura química es esquematizada en la figura 5.2: Figura 5.2. Estructura química del [ 18 F]-FDG. La 18-FDG es sin duda el radiofármaco PET más importante. Esto se debe no sólo a su aplicación al estudio de patologías muy diversas, sino también a sus características metabólicas, ya que, la FDG al igual que la glucosa, atraviesa la barrera hemato-encefálica y entra fácilmente en las células y a su rapidez de su síntesis. (20) Características físicas del flúor-18 El flúor (F-18) decae a través de emisión positrónica (β+) y tiene una vida media de 109,7 minutos, su esquema de decaimiento se muestra en la figura 5.3. Los fotones usados para el diagnóstico por imagen son los fotones gammas de 511 kev resultantes de la aniquilación de un positrón con un electrón del medio. Tabla 5.3. [Kocher, D.C. "Radioactive Decay Tables" DOE/TIC-11026, 89 (1981)] Radiación/Emisión % por desintegración Energía Media (kev) Positrón (β + ) 96,73 249,8 Gammas producto de la aniquilación de un positrón 193, Tabla 5.3. Datos principales acerca de la emisión de F-18 [Kocher, D.C. "Radioactive Decay Tables" DOE/TIC-11026, 89 (1981)]. 6162 CAPÍTULO 5 SEGMENTACIÓN Y PROCESAMIENTO DE IMÁGENES Figura 5.3. Esquema de decaimiento del F-18. [NIST, Standard Reference Data] El flúor-18 es sin duda el radionucleido que cuantitativamente más se utiliza en tomografía de emisión de positrones. No obstante, este hecho se debe casi exclusivamente a la utilización de la 18-FDG. Entre algunas de sus características nucleares y químicas más importantes, se encuentran: (20) Puede producirse en cantidades de hasta varios curios sin dificultad Su energía de emisión del positrón es de tan sólo 0,64 MeV (la menor de todos los emisores de positrones con aplicación en PET), lo que conlleva diversas ventajas: i. La dosis recibida por el paciente será menor que en los casos en que la energía de emisión del positrón sea mayor; dado que la energía depositada hasta la aniquilación será menor. ii. La distancia recorrida por el positrón hasta su aniquilación será asimismo reducida, lo que permitirá obtener imágenes con mayor resolución porque la incertidumbre entre el punto de desintegración del radionucleido y de aniquilación del positrón será menor. En la desintegración del 18-F no se emiten rayos gamma que puedan interferir la detección de los fotones de aniquilación, ni partículas (β, α) que puedan suponer un incremento en la dosis recibida por el paciente. 6263 CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT Su período de semi-desintegración (109,8 minutos) permite: i. Llevar a cabo síntesis complejas. ii. Aplicar protocolos de estudio mediante PET que se prolonguen varias horas, lo que facilita la realización de estudios de farmacocinética (estudios dinámicos) y análisis de metabólicos. iii. El transporte de radiofármacos marcados con flúor-18 a centros satélite u hospitales que dispongan de un tomógrafo PET o una cámara gamma de coincidencia Proceso de simulación del fantoma. Tal como fue mencionado en el capítulo 3, GATE permite la voxelización de los datos de los pacientes a través de la información presente en los estudios tomográficos y de PET. La lectura de las imágenes es solo posible, si éstas son tipo interfile (la información del fantoma está asociado a dos archivos: un archivo cabecera donde se guardan todos los datos del estudio y del paciente, y un archivo imagen, donde se almacenan las imágenes) y sus datos son enteros sin signo de 16 bits. En la mayoría de los servicios estos estudios son guardados como imágenes DICOM de 16 bits enteros con signo, por lo que debe realizarse un cambio de formato; conversión que puede ser hecha con ayuda de los softwars (X)MedCon, PMOD e ImageJ ; los dos primeros son los encargados de convertir las imágenes DICOM a interfile e ImagenJ, modifica el tipo de datos del estudio. Antes de comenzar el proceso de voxelización, es necesario elegir el algoritmo de simulación. GATE ofrece tres algoritmos de simulación: parameterizedboxmatrix, compressedmatrix y regularmatrix. El primero de estos, busca en todos los voxeles del fantoma, el voxel de mayor probabilidad para interactuar con la partícula, analizando cada uno de los voxeles de fantoma. El compressedmatrix genera un fantoma comprimido de tamaño de voxel variable. En este algoritmo de compresión, todos los voxeles adyacentes del mismo material se fusionan para formar el voxel rectangular más grande posible. 6364 CAPÍTULO 5 SEGMENTACIÓN Y PROCESAMIENTO DE IMÁGENES Y el regularmatrix, el cual únicamente tiene en cuenta los voxeles vecinos para determinar las posibles interacciones de los fotones, reduciendo drásticamente el tiempo de CPU necesario. Una vez seleccionado el algoritmo de simulación comienza el proceso de voxelización. Inicialmente, los datos de los estudios de PET, que se encuentran en escala de gris, son convertidos a niveles de actividad utilizando un traductor de rango, el cual contiene la actividad por rango de voxel, obtenida de medidas hechas en las imágenes y que dependen de la actividad inicial suministrada al paciente. Posteriormente estos niveles de actividad se utilizan para determinar el número de partículas primarias de cada voxel. Un ejemplo de un macro para la lectura de imágenes Interfile como fuente de distribuciones, se detalla a continuación, donde la fuente de emisión se definió para flúor-18, considerando prioritarios los fotones cuya aniquilación se generan a 180 grados. La fuente no se encuentra confinada a ningún órgano, la distribución espacial de los puntos de emisión es especificada a partir de la forma de la fuente, y la distribución angular de los ángulos de emisión es isótropa.(8) /gate/source/addsource sourcef18 voxel /gate/source/sourcef18/reader/insert interfile /gate/source/sourcef18/interfilereader/translator/insert range /gate/source/sourcef18/interfilereader/rangetranslator/readtable activity_range_brain1.dat /gate/source/sourcef18/interfilereader/rangetranslator/describe 1 /gate/source/sourcef18/interfilereader/readfile RiverosPET.h33 /gate/source/sourcef18/setposition mm /gate/source/sourcef18/settype backtoback /gate/source/sourcef18/gps/particle e+ /gate/source/sourcef18/setforcedunstableflag true /gate/source/sourcef18/gps/energytype Fluor18 6465 CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT /gate/source/sourcef18/setforcedhalflife 6586 s /gate/source/sourcef18/gps/monoenergy MeV /gate/source/sourcef18/gps/confine NULL /gate/source/sourcef18/gps/angtype iso /gate/source/sourcef18/dump 1 /gate/source/list Un ejemplo del traductor de rango de actividad se muestra a continuación. La primera línea define el número de subdivisiones utilizadas y las líneas siguientes definen los intervalos de voxel y sus correspondientes actividades, en becquereles. Como resultado se obtiene que la fuente voxelizada esta formada por una serie de valores discretos de actividad En el caso del fantoma utilizado en este trabajo, cada órgano fue segmentado y posterizado (disminución de la escala de gris), con ayuda de los software Amira e ImageJ, respectivamente, para poder lograr una fácil definición de la concentración de actividad y del coeficiente de atenuación tejido a usar. Una vez definida la distribución de la actividad en el fantoma a través de los datos del PET, es necesario realizar de nuevo el proceso anteriormente mencionado para asignar los materiales y coeficientes de atenuación a los diferentes órganos que conforman el fantoma, a través de la voxelización de las imágenes tomográficas. GATE ofrece una amplia biblioteca de materiales, apéndice 3, que contiene toda la información necesaria para poder asignar propiedades nucleares a los diferentes set de datos, dichas bibliotecas son las mismas que utiliza Geant4. Cada material, está definido por una combinación de elementos, los cuales son definidos como en la tabla periódica, es decir, se almacena el nombre, símbolo, número atómico y la masa molar de cada elemento. La buena definición de estos materiales es crítica, ya que esta información es 6566 CAPÍTULO 5 SEGMENTACIÓN Y PROCESAMIENTO DE IMÁGENES utilizada para determinar todas las interacciones de partículas que tienen lugar durante la simulación. Estas combinaciones de elementos requieren la definición de cuatro parámetros adicionales, estos son: el nombre del material, la densidad, elemento constitutivo y su abundancia individual. (8) Un ejemplo del mapa de atenuación utilizado para la simulación en GATE se muestra a continuación: /gate/world/daughters/name phantom /gate/world/daughters/insert regularmatrix /gate/phantom/geometry/insertreader interfile /gate/phantom/interfilereader/inserttranslator range /gate/phantom/interfilereader/rangetranslator/readtable range_atten_brain1.dat /gate/phantom/interfilereader/rangetranslator/describe 1 /gate/phantom/interfilereader/readfile RiverosCT.h33 /gate/phantom/placement/settranslation mm /gate/phantom/attachvoxelphantomsd El traductor utilizado para el mapa de atenuación, el cual, ahora da información de los diferentes materiales por voxel, se muestra a continuación. La primera línea define el número de subdivisiones utilizadas y las líneas siguientes definen los intervalos de voxel y el nombre del material correspondiente Air Body Water Body PMMA Body PMMA Body Water Body 6667 CÁLCULOS DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT Water Body Water Body Intestine Body Brain Body 6768 CAPÍTULO 6 DISCUSIÓN Y RESULTADOS Las simulaciones de Monte Carlo son métodos eficaces para la elaboración de mapas dosimétricos útiles en el área de medicina nuclear, ya que permite recrear de una forma realista los datos del paciente o fantoma y realizar una evaluación voxel a voxel, eliminando las limitantes que presenta el método MIRD, método válido únicamente cuando se tiene distribuciones homogéneas de actividad en órgano y cuando las dimensiones de dicho órganos son infinitas, algo que no se cumple cuando hablamos de un paciente, por lo que se hacen demasiadas aproximaciones. Son muchas las ventajas que posee este método estadístico, pero estas son en muchos casos opacados por dos grandes desventajas: su complejidad y alto consumo de tiempo de cálculo requerido para obtener resultados con bajas incertidumbres. La primera de estas es solucionada por GATE, ya que posee un lenguaje simple y fácil de implementar.69 CAPÍTULO 6 DISCUSIÓN Y RESULTADOS CÁLCULOS En este capítulo se tratara de analizar posibles soluciones para resolver el tiempo de cálculo, y poder obtener simulaciones con una resolución aceptable. El tiempo de cálculo depende de cinco factores: El tamaño de la matriz de adquisición de los estudios de PET y tomografía. La cantidad de actividad que se ingresa en el traductor. El tiempo de adquisición. El algoritmo de reconstrucción. Segmentación. De la buena combinación de estos factores depende en gran parte cuanto tiempo tardará la simulación Tamaño de Matriz. El tamaño de la matriz del mapa dosimétrico es dependiente del tamaño de la matriz del fantoma. Cuando se habla del tamaño de la matriz, se hace referencia a la resolución con la que puede distinguirse la dosis en cada órgano, ya que este tamaño, determina las dimensiones del voxel. Si la matriz es pequeña, el tamaño de voxel es grande y la exactitudes de la medida es baja, al contarse con poca estadística. Por esta razón, este parámetro debe ser elegido con precaución. Igualmente, el tamaño de voxel interviene en la definición de los traductores de rango tanto para los coeficientes atenuación como para la distribución de actividad. Mientras más grande sea el voxel, es decir menor sea el tamaño de matriz, las aproximaciones hechas por el programa para la evaluación de las trayectorias de los fotones entre voxeles es mayor. Para evaluar la dependencia entre la resolución de la imagen y el tamaño de la matriz, se hicieron pruebas con diferentes tamaños de matriz, 64, 128 y 512, conservando los mismos datos de entrada (actividad, tiempo de adquisición y paciente) y así poder tener criterios en la elección de dicho parámetro. En la figura 6.1 se ven los resultados obtenidos. 6970 DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT a. b. c. Figura 6.1.Mapas dosimétricos de tórax construido con la misma actividad y con igual tiempo de adquisición (1 seg.), variando el tamaño de matriz de entrada. a. Matriz de entrada de 64 x 64, b. Matriz de entrada de 128 x 128, c. Matriz de entrada de 512 x 512. También, fue calculado el tiempo que tarda cada simulación al variar el tamaño de matriz. Los resultados se muestran en la Tabla 6.1 Matriz Tiempo (min) 64x64 0,42 128x128 0,85 512x512 3,42 Tabla 6.1.Medidas de tiempo de duración para cada simulación variando el tamaño de matriz y conservando los demás parámetros como constantes. 7071 CAPÍTULO 6 DISCUSIÓN Y RESULTADOS CÁLCULOS De los resultados obtenidos en la figura 6.1 y en la tabla 6.1, se observa que aun cuando el tiempo de cálculo es mucho mayor para el caso de la matriz de 512 x 512, su resolución también es mayor, razón por la cual se decidió sacrificar tiempo de cálculo para garantizar una buena resolución en los mapas dosimétricos, por lo que fue elegida como matriz de entrada para el fantoma, la matriz de Actividad. Otro factor determinante cuando hablamos de tiempo de cálculo, es el número de partículas a considerar. Un gran número de partículas indica un gran consumo de tiempo en la simulación, pero una baja incertidumbre. La determinación del número de eventos a analizar viene implícita en el archivo traductor de rango utilizado para la definición de los rangos de actividad o cantidad de partículas por órgano. Por ejemplo, para un examen típico de PET/CT, al paciente se le administra entre 3mCi-20mCi de [ 18 F]-FDG, según sea su peso, un 10% del radio-marcador es absorbido por el cerebro y el resto se distribuye en los demás órganos. Realizando cuentas al aire, se puede suponer que cada órgano tiene una actividad de [ 18 F]-FDG promedio del orden de los µci, lo cual equivaldría a tener aproximadamente en 40 kbq o eventos por región de interés. Desafortunadamente, al no contar con la infraestructura necesaria para realiza un cálculo con actividades de este orden, se recurrió a un análisis cualitativo de dosis, utilizando niveles de actividad inferiores a los µci. Valores que serán especificados posteriormente Tiempo de adquisición. El tiempo de adquisición, como se menciono en el capítulo 3, se refiere a la duración virtual de la simulación, la cual debe ser equivalente a los tiempos de adquisición reales. En Fuesmen, una camilla o corte del PET tiene una duración entre 60 a 420 seg., dicha variación es determinada por el tipo de lección, la región de interés y la experiencia del físico y el médico. La relación que existe entre el tiempo de adquisición y tiempo de cálculo, fue analizada y se determino que existe una proporcionalidad lineal, entre ambos 7172 DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT parámetros, Tabla 6.2, es decir, si una simulación con un tiempo de adquisición de 1 seg. tarda alrededor de 1 hora, al duplicar el tiempo de adquisición el tiempo de cálculo también lo hace. Tiempo de Adquisición (s) Tiempo de cálculo (h) 0,1 0, Tabla 6.2 Tiempos de cálculo promedio de las diferentes simulaciones Algoritmo de simulación. Como fue mencionado en el capítulo 5, GATE ofrece tres algoritmos de simulación: parameterizedboxmatrix, compressedmatrix y regularmatrix, cada uno de los cuales tiene sus pros y sus contras. El primero de estos, al buscar en todos los voxeles del fantoma, el voxel de mayor probabilidad para interactuar con la partícula, requiere un alto consumo de CPU y demanda de mucho tiempo de cálculo. El compressedmatrix al generar fantomas comprimidos reduciendo el tiempo de cálculo y la memoria necesaria, a costa de una baja resolución. Por último, el regularmatrix al evaluar solamente los vecinos del voxel que contiene el fotón, reduce el tiempo de cálculo, manteniendo una buena resolución. Por esta razón, el algoritmo que fue seleccionado para realizar la simulación es el regularmatrix, ya que garantiza buena resolución con tiempos de cálculo relativamente bajos. 7273 CAPÍTULO 6 DISCUSIÓN Y RESULTADOS CÁLCULOS 6.5- Segmentación. Para garantizar la correcta voxelización de los datos del paciente es necesario realizar una segmentación previa. Esta operación, es fundamental para identificar los diferentes órganos de interés que serán evaluados dentro de los mapas dosimétricos. En la figura 6.2 se muestra el proceso de segmentación que se realizo, con ayuda del software Amira. a. b. c. Figura 6.2.Proceso de segmentación de un estudio de PET de tórax, a. Imagen original, b. selección de órganos de interés, a través del Amira, c. Segmentación completa del estudio PET. Los datos de entrada para la elaboración de los mapas dosimétricos, fueron tomados de tres pacientes diferentes. Cada uno de los pacientes representa una región anatómica diferente, dicho corte se realizo para poder obtener datos dosimétrico 7374 DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT estadísticamente relevantes en el menor tiempo posible. Las regiones de análisis fueron: pelvis, tórax y cabeza y cuello, regiones de gran importancia por su alta concentración de actividad. En las Tablas 6.3 son enumerados los órganos de estudio con su respectivo voxel, para cada uno de los casos. Región Anatómica Órgano Voxel Cabeza Y Cuello Tórax Pelvis Cerebro 339 Cristalino 269 Garganta 256 Pulmones 310 Tiroides 349 Medula 246 Pulmones 32 Medula 216 Riñones 305 Bazo 203 Estomago 254 Hígado 246 Intestino 373 Corazón 187 Intestino 381 Recto 210 Próstata 309 Vejiga 174 Orina 183 Medula 398 Tabla 6.3 Órganos de estudio para cada una de las regiones de interés. Una vez definido cuales son los órganos de estudio, es necesario definir cuáles de ellos poseen actividad (órganos fuente). Estos, son fácilmente reconocidos en las imágenes de PET. En tabla 6.4, se muestran los órganos seleccionados como fuentes con su respectiva actividad. En la figura 6.3, se muestra cada una de las regiones con su respectiva segmentación. Es importante resalta el hecho, que ante la complejidad en la segmentación, el intestino no fue segmentado en su totalidad y los genitales no fueron tomados. 7475 CAPÍTULO 6 DISCUSIÓN Y RESULTADOS CÁLCULOS Región Anatómica Órgano Fuente Voxel Actividad (Bq) Cabeza Y Cuello Cerebro Tórax Pelvis Corazón Riñones Hígado Intestino Orina Intestino Tabla 6.4 Órganos fuentes de cada una de las regiones de interés a estudiar, con su respectiva actividad. a. b. c. Figura 6.3. Segmentación de las diferentes regiones de estudio, a. Pelvis, b. Cabeza y Cuello, c. Tórax 7576 DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT 6.6- Mapas Dosimétricos. Todos los parámetros antes mencionados son claves para la formación de los mapas dosimétricos, una mala combinación de estos, se verá reflejado en mapas erróneos y tiempo de cálculo extremadamente largos. En la Tabla 6.5, se muestran los parámetros seleccionados, para la elaboración del presente trabajo. Parámetro Descripción Tamaño de matriz Fue elegido siguiendo el protocolo utilizado en FUESMEN, donde la matriz del CT es de 512x512 y la matriz de PET es de 128 x 128. Actividad Tiempo de adquisición Algoritmo de simulación Fueron utilizadas actividades del orden de los kbq. (Posteriormente, serán detallados estos valores) Para poder realizar un análisis temporal de la dosis fueron tomados cuatro diferentes tiempos 100 mseg, 1 seg. 10 seg y 60seg En aras de optimizar la simulación, el algoritmo seleccionado para la simulación es "RegularMatrix" Tabla 6.5.Parametros utilizados para la optimización de las simulaciones de Monte Carlo. Una vez seleccionado todos los parámetros de simulación y segmentadas las diferentes regiones anatómicas, se dispuso a iniciar el proceso de simulación. Los resultados obtenidos para las diferentes regiones anatómicas y distintos tiempos de adquisición son mostrados en las figuras 6.4 a 6.6. Igualmente en las Tablas 6.6 a 6.8 son presentados los niveles de dosis por órgano, valores, obtenidos con ayuda del software Amide. Los mapas dosimétricos construidos por GATE, son archivos tipo binario flotantes de 4 bytes, las unidades de la dosis se encuentran cgy. La matriz de imagen tiene las mismas dimensiones que el fantoma. Para el caso de estudio de la presente tesis únicamente se tuvo en cuenta el decaimiento físico de la fuente, en el cálculo dosimétrico. 7677 CAPÍTULO 6 DISCUSIÓN Y RESULTADOS CÁLCULOS Como se menciono en el capítulo 3, el software GATE, calcula niveles dosimétricos voxel a voxel a través de la ecuación 3.2, lo cual permite obtener un valor de dosis mucho más preciso y realista. Partiendo de estos valores es posible calcular el valor de dosis por órgano, lo cual es el fin último de estas simulaciones. La dosis en un determinado órgano se calcula como: (6.1) Donde es la energía absorbida por órgano y m es la masa del órgano de estudio. Expandiendo esta expresión, se tiene: (6.2) Donde el subíndice, i, equivale a los voxel que conforman el órgano, entonces y, son la energía absorbida y la masa, respectivamente, de cada voxel. Como se está evaluando la dosis por órgano, se puede asumir que la masa de cada voxel es igual, por lo tanto, donde N es el número de voxel dentro del órgano. Rescribiendo la ecuación 6.2, se tiene que: (6.3) Al tomar la masa de cada voxel como una constante, ésta puede entrar en la sumatoria, entonces la sumatoria de energía se convierte en una sumatoria de dosis por voxel, por lo que la ecuación 6.3, es equivalente a = (6.4) Así la dosis por órgano, no es más que el promedio de la dosis por voxel, de cada órgano. La expresión 6.4 solo es válida para casos donde la masa sea constante. 7778 DOSIMÉTRICOS MEDIANTE CÓDIGO MONTE CARLO A PARTIR DE IMÁGENES DE PET/CT a. b. c. d. Figura 6.4. Mapas dosimétricos de cabeza y cuello para diferentes tiempos de adquisición (T a ), a. T a = 100 ms, b. T a = 1 s, c. T a = 10 s, d. T a = 60 s Cabeza y cuello Dosis (pgy) 100 ms 1 s 10 s 60s Cerebro 7,90 77, Cristalino 8,17 77, Medula 3,24 62, Garganta 6,56 65, Tiroides 4,21 41, Pulmones 2,82 27, Tabla 6.6. Cálculo dosimétrico por órgano, para la región de la cabeza y cuello a diferentes tiempos de adquisición. 7879 CAPÍTULO 6 DISCUSIÓN Y RESULTADOS CÁLCULOS a. b. c. d. Figura 6.5. Mapas dosimétricos de tórax para diferentes tiempos de adquisición (T a ), a. T a = 100 ms, b. T a = 1 s, c. T a = 10 s, d. T a = 60 s Tórax Dosis (pgy) 100 ms 1 s 10 s 60s Corazón 1, Medula 1, Hígado 0, Bazo 1, Estomago 0, Riñones 0, Intestino 0, Pulmones 1, Tabla 6.7. Cálculo dosimétrico por órgano, para la región torácica a diferentes tiempos de adquisición. 79 Mostrar más
SIMULACION DE EVENTOS CLAVE : LII 210 PROFESOR: MTRO. ALEJANDRO SALAZAR GUERRERO 1 1. LA SIMULACIÒN DE EVENTOS DISCRETOS. 1.1 Definición de Eventos Discretos. 1.2 Estructura de la simulación de eventos Más detalles SIMULACIÓN. Texto Base de Teoría. Alfonso Urquía Moraleda
SIMULACIÓN Texto Base de Teoría Alfonso Urquía Moraleda Departamento de Informática y Automática Escuela Técnica Superior de Ingeniería Informática, UNED Juan del Rosal 16, 28040 Madrid, España E-mail: Más detalles OPTIMIZACIÓN DE PROTOCOLO CT, EN PET-CT, DE CUERPO ENTERO
OPTIMIZACIÓN DE PROTOCOLO CT, EN PET-CT, DE CUERPO ENTERO Fredys Santos Gutiérrez 1, Mauro Namías 2 1 Área Control de Calidad y Protección Radiológica (ACCPR) Caja Costarricense de Seguro Social (CCSS) Más detalles UNIVERSIDAD NACIONAL DE TRUJILLO
UNIVERSIDAD NACIONAL DE TRUJILLO FACULTAD DE CIENCIAS FISICAS Y MATEMATICAS ESCUELA ACADEMICO PROFESIONAL DE INFORMÁTICA DESARROLLO DE UN SISTEMA INTELIGENTE BASADO EN ALGORITMOS GENETICOS PARA INCREMENTAR Más detalles Guía de Seguridad del CSN GS-11.03
Guía de Seguridad del CSN GS-11.03 Sobre Metodología para la evaluación del impacto radiológico de las industrias NORM 12 de diciembre de 2012 ÍNDICE PREÁMBULO... 3 1. OBJETO... 4 2. ÁMBITO DE APLICACIÓN... Más detalles INTRODUCCIÓN MARCO PROBLÉMICO JUSTIFICACIÓN 2. OBJETIVOS 2.1 OBJETIVO GENERAL 2.2 OBJETIVOS ESPECÍFICOS 3. MARCO REFERENCIAL 3.
Universidad de Valladolid E. U. DE INFORMÁTICA (SEGOVIA) Ingeniería Técnica en Informática de Gestión Modelos de Fiabilidad del Software Alumno: Jesús Plaza Rubio Tutores: Pablo Fernández Blanco Pilar Más detalles Proyecto Fin de Carrera
Proyecto Fin de Carrera CÁLCULO COMPUTACIONAL DE PERFILES AERODINÁMICOS OS EN AEROTURBINA Víctor M. Martínez Fernández Director: Antonio Pascau Benito Departamento de Ciencia y Tecnología de Materiales Más detalles Universidad Tecnológica de Panamá. Facultad de Ingeniería Industrial. Departamento de Producción
http://www.eugenio-soler-galan.es/aleatoriedad Trabajo de investigación I Profesor: Benjamín Ramos Álvarez junio de 2008 Página 2 de 54 INDICE 1. MOTIVACIÓN... pág. 4 2. DÓNDE ENCONTRAMOS ALEATORIEDAD?... Más detalles AUTORES OTROS INTEGRANTES DEL PROYECTO
AUTORES Dr. Gabriel González Sprinberg. Profesor titular del Instituto de Física de la Facultad de Ciencias (Universidad de la República, Uruguay). Cursó sus estudios de licenciatura y doctorado en física Más detalles SIMULACIÓN DEL FLUJO TURBULENTO EN UN CODO DE 90 DE TUBERIA DE SECCIÓN CIRCULAR MANUEL ALEJANDRO TABORDA CEBALLOS
SIMULACIÓN DEL FLUJO TURBULENTO EN UN CODO DE 90 DE TUBERIA DE SECCIÓN CIRCULAR MANUEL ALEJANDRO TABORDA CEBALLOS UNIVERSIDAD AUTONOMA DE OCCIDENTE FACULTAD DE INGENIERIA DEPARTAMENTO DE ENERGETICA Y MECANICA Más detalles FUNDAMENTACIÓN EN RIESGO CON RISK SIMULATOR
FUNDAMENTACIÓN EN RIESGO CON RISK SIMULATOR TABLA DE CONTENIDO 1. RIESGO... 3 2. TIPOS DE RIESGO... 4 2.1 RIESGO DE MERCADO... 4 2.1.1 Valor en Riesgo... 6 2.1.2 Cópulas... 6 2.1.3 Teoría del Valor Extremo... Más detalles Cartagena de Indias, 16 de julio de 2013
Cartagena de Indias, 16 de julio de 2013 Señores: COMITÉ EVALUADOR DE PROYECTOS DE GRADO Programa de Ingeniería Civil Facultad de Ingeniería Universidad de Cartagena E. S. D. Cordial saludo, Por medio Más detalles Diseño de una estrategia para la gestión energética de un vehículo híbrido enchufable
UNIVERSIDAD DE ALMERÍA ESCUELA POLITÉCNICA SUPERIOR Y FA. CC. EXPERIMENTALES Titulación de Ingeniería Técnica Industrial Mecánica Diseño de una estrategia para la gestión energética de un vehículo híbrido Más detalles UNIVERSIDAD SIMÓN BOLÍVAR DECANATO DE ESTUDIOS PROFESIONALES COORDINACIÓN DE INGENIERÍA QUÍMICA
UNIVERSIDAD SIMÓN BOLÍVAR DECANATO DE ESTUDIOS PROFESIONALES COORDINACIÓN DE INGENIERÍA QUÍMICA ESTUDIO FLUIDODINÁMICO DE MEZCLA BIFÁSICA EN MEDIDOR TIPO VENTURI VERTICAL Por: Carlos Pedron PROYECTO DE Más detalles Evaluación de datos de medición Guía para la Expresión de la Incertidumbre de Medida
Evaluación de datos de medición Guía para la Expresión de la Incertidumbre de Medida k A + B EDICIÓN DIGITAL en español (traducción ª Ed. Sept. 008) Primera edición Septiembre 008 (original en inglés) Más detalles Informática Aplicada a la Ingeniería de Procesos I (Orientación I) Data Mining
Universidad Tecnológica Nacional, Facultad Regional Rosario Departamento de Ingeniería Química Grupo de Investigación Aplicada a la ingeniería Química (GIAIQ) Informática Aplicada a la Ingeniería de Procesos Más detalles UNIVERSIDAD NACIONAL AUTÓNOMA DE MÉXICO
UNIVERSIDAD NACIONAL AUTÓNOMA DE MÉXICO FACTULTAD DE ESTUDIOS SUPERIORES ACATLÁN MINERIA DE DATOS CON APLICACIONES TESIS QUE PARA OBTENER EL TITULO DE LICENCIADA EN MATEMÁTICAS APLICADAS Y COMPUTACIÓN Más detalles "The solutions are all simple, after you have arrived at them. But they're simple only when you know already what they are. -Robert M.
INSTITUTO POLITÉCNICO NACIONAL CENTRO DE INVESTIGACIÓN Y DESARROLLO DE TECNOLOGÍA DIGITAL MAESTRÍA EN CIENCIAS EN SISTEMAS DIGITALES PLANEACIÓN DE TRAYECTORIAS PARA ROBÓTICA MÓVIL MEDIANTE OPTIMIZACIÓN Más detalles Este documento espera reflejar todo el trabajo realizado en este tiempo y servir de guía y explicación de la aplicación realizada.
. INTRODUCCIÓN Este documento espera reflejar todo el trabajo realizado en este tiempo y servir de guía y explicación de la aplicación realizada. Inicialmente este proyecto forma parte del proyecto global Más detalles Modelado y simulación de un robot Robucar TT
Proyecto Fin de Carrera Ingeniería Industrial Automatización Industrial y Robótica Modelado y simulación de un robot Robucar TT Carlos Gracia Sola Director: Luis Montano Gella Departamento de Informática Más detalles Bioestadística: Unidad de Demografía UNIVERSIDAD DE LA FRONTERA FACULTAD DE MEDICINA DEPARTAMENTO SALUD PUBLICA.
Bioestadística: Unidad de Demografía UNIVERSIDAD DE LA FRONTERA FACULTAD DE MEDICINA DEPARTAMENTO SALUD PUBLICA. DINÁMICA DE POBLACIONES VARIABLES DEMOGRAFICAS Los hechos más importantes que determinan Más detalles UNIVERSIDAD AUTÓNOMA DE BAJA CALIFORNIA SUR ÁREA DE CONOCIMIENTO DE CIENCIAS DEL MAR DEPARTAMENTO ACADÉMICO DE SISTEMASCOMPUTACIONALES TESIS
UNIVERSIDAD AUSTRAL DE CHILE FACULTAD DE CIENCIAS DE LA INGENIERÍA ESCUELA DE INGENIERÍA CIVIL EN INFORMÁTICA Aplicación de Data Mining sobre un GIS de caracterización forestal. TESIS DE GRADO PARA OPTAR Más detalles Algoritmos Evolutivos
Algoritmos Evolutivos 2003-11-20 Cristóbal Romero, Pedro González y Sebastián Ventura 1.1 Introducción En general, cualquier tarea abstracta a realizar puede considerarse como la resolución de un problema Más detalles CAPÍTULO 4 ARITMÉTICA Y REPRESENTACIÓN DE LA INFORMACIÓN EN EL COMPUTADOR
CAPÍTULO 4 ARITMÉTICA Y REPRESENTACIÓN DE LA INFORMACIÓN EN EL COMPUTADOR Dos de los aspectos básicos que se presentan en el tratamiento de la información son cómo representarla (de lo cual dependerá sus Más detalles Francisco J. Ferrer Jesús Aguilar Joaquín Peña. Departamento de Lenguajes y Sistemas Informáticos Facultad de Informática Universidad de Sevilla
Francisco J. Ferrer Jesús Aguilar Joaquín Peña Departamento de Lenguajes y Sistemas Informáticos Facultad de Informática Universidad de Sevilla INFORME TÉCNICO: LSI-2000-08 Julio, 2000 1. Introducción Más detalles TREBALL FI DE CARRERA TRABAJO FINAL DE CARRERA
TREBALL FI DE CARRERA TRABAJO FINAL DE CARRERA Evaluación y optimización de los servicios de un taller de mantenimiento mediante el software de simulación ARENA. Nombre estudiante: David Llerena Santana Más detalles IMPLANTACIÓN SISTEMA DE GESTIÓN TOTAL DE LA CALIDAD EN UNA FÁBRICA DE MOTORES

References: resolución 
 resolución 
 Resolución 
 Resolución 
 Resolución 
 Resolución 
 Resolución 
 Resolución 
 Resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución