Development and operation of a pixel segmented liquid-filled linear array for radiotherapy quality assurance

A liquid isooctane (C$_{8}$H$_{18}$) filled ionization linear array for radiotherapy quality assurance has been designed, built and tested. The detector consists of 128 pixels, each of them with an area of 1.7 mm $\times$ 1.7 mm and a gap of 0.5 mm. The small pixel size makes the detector ideal for high gradient beam profiles like those present in Intensity Modulated Radiation Therapy (IMRT) and radiosurgery. As read-out electronics we use the X-Ray Data Acquisition System (XDAS) with the Xchip developed by the CCLRC. Studies concerning the collection efficiency dependence on the polarization voltage and on the dose rate have been made in order to optimize the device operation. In the first tests we have studied dose rate and energy dependences, and signal reproducibility. Dose rate dependence was found lower than 2.5 % up to 5 Gy min$^{-1}$, and energy dependence lower than 2.1 % up to 20 cm depth in solid water. Output factors and penumbras for several rectangular fields have been measured with the linear array and were compared with the results obtained with a 0.125 cm$^{3}$ air ionization chamber and radiographic film, respectively. Finally, we have acquired profiles for an IMRT field and for a virtual wedge. These profiles have also been compared with radiographic film measurements. All the comparisons show a good correspondence. Signal reproducibility was within a 2% during the test period (around three months). The device has proved its capability to verify on-line therapy beams with good spatial resolution and signal to noise ratio.Comment: 16 pages, 12 figures Submitted to Phys. Med. Bio

PTV margin for dose escalated radiation therapy of prostate cancer with daily on-line realignment using internal fiducial markers: Monte Carlo approach and dose population histogram (DPH) analysis.

Using internal fiducial markers and electronic portal imaging (EPID) to realign patients has been shown to significantly reduce positioning uncertainties in prostate radiation treatment. This creates the possibility of decreasing the planning target volume (PTV) margin added on the clinical target volume (CTV), which in turn may allow for dose escalation. We compared the outcome of two plans: 70Gy/35fx, 10 mm PTV margin without patient realignment (Reference Plan) vs. 78Gy/39fx, 5 mm PTV margin with patient realignment (Escalated Plan). Four-field-oblique (gantry angles 35 degrees, 90 degrees, 270 degrees, 176 degrees, 325 degrees) beam arrangement was used. Monte Carlo code was used to simulate the daily organ motion. Dose to each organ was calculated. Tumour control probability (TCP) and the effective dose to critical organ (Deff) were calculated using the biologically normalized dose-volume histograms. By comparing the biological factors, we found that the prescription dose can be escalated to 78Gy/39fx with a 5 mm PTV margin when using internal fiducial markers and EPID. Based on the available dose-response data for intermediate risk prostate patients, this will result in a 20% increase of local control and significantly reduced rectal complications provided that less serial dose-volume behaviour of rectum is proven

Radiation therapy calculations using an on-demand virtual cluster via cloud computing

Computer hardware costs are the limiting factor in producing highly accurate radiation dose calculations on convenient time scales. Because of this, large-scale, full Monte Carlo simulations and other resource intensive algorithms are often considered infeasible for clinical settings. The emerging cloud computing paradigm promises to fundamentally alter the economics of such calculations by providing relatively cheap, on-demand, pay-as-you-go computing resources over the Internet. We believe that cloud computing will usher in a new era, in which very large scale calculations will be routinely performed by clinics and researchers using cloud-based resources. In this research, several proof-of-concept radiation therapy calculations were successfully performed on a cloud-based virtual Monte Carlo cluster. Performance evaluations were made of a distributed processing framework developed specifically for this project. The expected 1/n performance was observed with some caveats. The economics of cloud-based virtual computing clusters versus traditional in-house hardware is also discussed. For most situations, cloud computing can provide a substantial cost savings for distributed calculations.Comment: 12 pages, 4 figure

Four-dimensional Cone Beam CT Reconstruction and Enhancement using a Temporal Non-Local Means Method

Four-dimensional Cone Beam Computed Tomography (4D-CBCT) has been developed to provide respiratory phase resolved volumetric imaging in image guided radiation therapy (IGRT). Inadequate number of projections in each phase bin results in low quality 4D-CBCT images with obvious streaking artifacts. In this work, we propose two novel 4D-CBCT algorithms: an iterative reconstruction algorithm and an enhancement algorithm, utilizing a temporal nonlocal means (TNLM) method. We define a TNLM energy term for a given set of 4D-CBCT images. Minimization of this term favors those 4D-CBCT images such that any anatomical features at one spatial point at one phase can be found in a nearby spatial point at neighboring phases. 4D-CBCT reconstruction is achieved by minimizing a total energy containing a data fidelity term and the TNLM energy term. As for the image enhancement, 4D-CBCT images generated by the FDK algorithm are enhanced by minimizing the TNLM function while keeping the enhanced images close to the FDK results. A forward-backward splitting algorithm and a Gauss-Jacobi iteration method are employed to solve the problems. The algorithms are implemented on GPU to achieve a high computational efficiency. The reconstruction algorithm and the enhancement algorithm generate visually similar 4D-CBCT images, both better than the FDK results. Quantitative evaluations indicate that, compared with the FDK results, our reconstruction method improves contrast-to-noise-ratio (CNR) by a factor of 2.56~3.13 and our enhancement method increases the CNR by 2.75~3.33 times. The enhancement method also removes over 80% of the streak artifacts from the FDK results. The total computation time is ~460 sec for the reconstruction algorithm and ~610 sec for the enhancement algorithm on an NVIDIA Tesla C1060 GPU card.Comment: 20 pages, 3 figures, 2 table

Improving plan quality and consistency by standardization of dose constraints in prostate cancer patients treated with CyberKnife.

Treatment plans for prostate cancer patients undergoing stereotactic body radiation therapy (SBRT) are often challenging due to the proximity of organs at risk. Today, there are no objective criteria to determine whether an optimal treatment plan has been achieved, and physicians rely on their personal experience to evaluate the plan's quality. In this study, we propose a method for determining rectal and bladder dose constraints achievable for a given patient's anatomy. We expect that this method will improve the overall plan quality and consistency, and facilitate comparison of clinical outcomes across different institutions. The 3D proximity of the organs at risk to the target is quantified by means of the expansion-intersection volume (EIV), which is defined as the intersection volume between the target and the organ at risk expanded by 5 mm. We determine a relationship between EIV and relevant dosimetric parameters, such as the volume of bladder and rectum receiving 75% of the prescription dose (V75%). This relationship can be used to establish institution-specific criteria to guide the treatment planning and evaluation process. A database of 25 prostate patients treated with CyberKnife SBRT is used to validate this approach. There is a linear correlation between EIV and V75% of bladder and rectum, confirming that the dose delivered to rectum and bladder increases with increasing extension and proximity of these organs to the target. This information can be used during the planning stage to facilitate the plan optimization process, and to standardize plan quality and consistency. We have developed a method for determining customized dose constraints for prostate patients treated with robotic SBRT. Although the results are technology specific and based on the experience of a single institution, we expect that the application of this method by other institutions will result in improved standardization of clinical practice

Automation of the Monte Carlo simulation of medical linear accelerators

La consulta íntegra de la tesi, inclosos els articles no comunicats públicament per drets d'autor, es pot realitzar prèvia petició a l'Arxiu de la UPCThe main result of this thesis is a software system, called PRIMO, which simulates clinical linear accelerators and the subsequent dose distributions using the Monte Carlo method. PRIMO has the following features: (i) it is self- contained, that is, it does not require additional software libraries or coding; (ii) it includes a geometry library with most Varian and Elekta linacs; (iii) it is based on the general-purpose Monte Carlo code PENELOPE; (iv) it provides a suite of variance-reduction techniques and distributed parallel computing to enhance the simulation efficiency; (v) it is graphical user interfaced; and (vi) it is freely distributed through the website http://www.primoproject.net In order to endow PRIMO with these features the following tasks were conducted: - PRIMO was conceived with a layered structure. The topmost layer, named the GLASS, was developed in this thesis. The GLASS implements the GUI, drives all the functions of the system and performs the analysis of results. Lower layers generate geometry files, provide input data and execute the Monte Carlo simulation. - The geometry of Elekta linacs from series SU and MLCi were coded in the PRIMO system. - A geometrical model of the Varian True Beam linear accelerator was developed and validated. This model was created to surmount the limitations of the Varian distributed phase-space files and the absence of released information about the actual geometry of that machine. This geometry model was incorporated into PRIMO. - Two new variance-reduction techniques, named splitting roulette and selective splitting, were developed and validated. In a test made with an Elekta linac it was found that when both techniques are used in conjunction the simulation efficiency improves by a factor of up to 45. - A method to automatically distribute the simulation among the available CPU cores of a computer was implemented. The following investigations were done using PRIMO as a research tool : - The configu ration of the condensed history transport algorithm for charged particles in PENELOPE was optimized for linac simulation. Dose distributions in the patient were found to be particularly sensitive to the values of the transport parameters in the linac target. Use of inadequate values of these parameters may lead to an incorrect determination of the initial beam configuration or to biased dose distributions. - PRIMO was used to simulate phase-space files distributed by Varian for the True Beam linac. The results were compared with experimental data provided by five European radiotherapycenters. It was concluded thatthe latent variance and the accuracy of the phase-space files were adequate for the routine clinical practice. However, for research purposes where low statistical uncertainties are required the phase-space files are not large enough. To the best of our knowledge PRIMO is the only fully Monte Carlo-based linac and dose simulation system , addressed to research and dose verification, that does not require coding tasks from end users and is publicly available.El principal resultado de esta tesis es un sistema informático llamado PRIMO el cual simula aceleradores lineales médicos y las subsecuentes distribuciones de dosis empleando el método de Monte Carlo. PRIMO tiene las siguiente características: (i) es auto contenido, o sea no requiere de librerías de código ni de programación adicional ; (ii) incluye las geometrías de los principales modelos de aceleradores Varían y Elekta; (iii) está basado en el código Monte Carlo de propósitos generales PENELOPE; (iv) contiene un conjunto de técnicas de reducción de varianza y computación paralela distribuida para mejorar la eficiencia de simulación; (v) tiene una interfaz gráfica de usuario; y (vi) se distribuye gratis en el sitio web http://vvww.primoproject.net. Para dotar a PRIMO de esas características, se realizaron las tareas siguientes: - PRIMO se concibió con una estructura de capas. La capa superior, nombrada GLASS, fue desarrollada en esta tesis. GLASS implementa la interfazgráfica de usuario, controla todas las funciones del sistema y realiza el análisis de resultados. Las capas inferiores generan los archivos de geometría y otros datos de entrada y ejecutan la simulación Monte Carlo. - Se codificó en el sistema PRIMO la geometría de los aceleradores Elekta de las series SLi y MLC. - Se desarrolló y validó un modelo geométrico del acelerador TrueBeam de Varian. Este modelo fue creado para superar las limitaciones de los archivos de espacio de fase distribuidos por Varian, así como la ausencia de información sobre la geometría real de esta máquina. Este modelo geométrico fue incorporado en PRIMO. - Fueron desarrolladas y validadas dos nuevas técnicas de reducción de varianza nombradas splitting roulette y selective splitting. En pruebas hechas en un acelerador Elekta se encontró que cuando ambas técnicas se usan en combinación, la eficiencia de simulación mejora 45 veces. - Se implementó un método para distribuir la simulación entre los procesadores disponibles en un ordenador. Las siguientes investigaciones fueron realizadas usando PRIMO como herramienta: - Fue optimizada la configuración del algoritmo de PENELOPE para el transporte de partículas cargadas con historia condensada en la simulación del linac. Se encontró que las distribuciones de dosis en el paciente son particularmente sensibles a los valores de los parámetros de transporte usados para el target del linac. El uso de va lores inadecuados para esos parámetros puede conducir a una incorrecta determinación de la configuración del haz inicial o producir sesgos en las distribuciones de dosis. - Se utilizó PRIMO para simular archivos de espacios de fase distribuidos por Varian para el linac TrueBeam. Los resultados se compararon con datos experimentales aportados por cinco centros de radioterapia europeos. Se concluyó que la varianza latente y la exactitud de estos espacios de fase son adecuadas para la práctica clínica de rutina. Sin embargo estos espacios de fase no son suficientemente grandes para emplearse en investigaciones que requieren alcanzar una baja incertidumbre estadística. Hasta donde conocemos, PRIMO es el único sistema Monte Carlo que simula completamente el acelerador lineal y calcula la dosis absorbida, dirigido a la investigación y la verificación de dosis que no requiere del usuario tareas de codificación y está disponible públicamentePostprint (published version