RESUMO
PURPOSE: Functional lung avoidance (FLA) radiation therapy (RT) aims to minimize post-RT pulmonary toxicity by preferentially avoiding dose to high-functioning lung (HFL) regions. A common limitation is that FLA approaches do not consider the conducting architecture for gas exchange. We previously proposed the functionally weighted airway sparing (FWAS) method to spare airways connected to HFL regions, showing that it is possible to substantially reduce risk of radiation-induced airway injury. Here, we compare the performance of FLA and FWAS and propose a novel method combining both approaches. METHODS: We used breath-hold computed tomography (BHCT) and simulation 4-dimensional computed tomography (4DCT) from 12 lung stereotactic ablative radiation therapy patients. Four planning strategies were examined: (1) Conventional: no sparing other than clinical dose-volume constraints; (2) FLA: using a 4DCT-based ventilation map to delineate the HFL, plans were optimized to reduce mean dose and V13.50 in HFL; (3) FWAS: we autosegemented 11 to 13 generations of individual airways from each patient's BHCT and assigned priorities based on the relative contribution of each airway to total ventilation. We used these priorities in the optimization along with airway dose constraints, estimated as a function of airway diameter and 5% probability of collapse; and (4) FLA + FWAS: we combined information from the 2 strategies. We prioritized clinical dose constraints for organs at risk and planning target volume in all plans. We performed the evaluation in terms of ventilation preservation accounting for radiation-induced damage to both lung parenchyma and airways. RESULTS: We observed average ventilation preservation for FLA, FWAS, and FLA + FWAS as 3%, 8.5%, and 14.5% higher, respectively, than for Conventional plans for patients with ventilation preservation in Conventional plans <90%. Generalized estimated equations showed that all improvements were statistically significant (P ≤ .036). We observed no clinically relevant improvements in outcomes of the sparing techniques in patients with ventilation preservation in Conventional plans ≥90%. CONCLUSIONS: These initial results suggest that it is crucial to consider the parallel and the serial nature of the lung to improve post-radiation therapy lung function and, consequently, quality of life for patients.
Assuntos
Neoplasias Pulmonares , Lesões por Radiação , Radiocirurgia , Tomografia Computadorizada Quadridimensional/métodos , Humanos , Pulmão/diagnóstico por imagem , Neoplasias Pulmonares/diagnóstico por imagem , Neoplasias Pulmonares/radioterapia , Qualidade de Vida , Lesões por Radiação/prevenção & controle , Planejamento da Radioterapia Assistida por Computador/métodosRESUMO
We report on the design, implementation and characterization of a multi-graphic processing unit (GPU) computational platform for higher-order optimization in radiotherapy treatment planning. In collaboration with a commercial vendor (Varian Medical Systems, Palo Alto, CA), a research prototype GPU-enabled Eclipse (V13.6) workstation was configured. The hardware consisted of dual 8-core Xeon processors, 256 GB RAM and four NVIDIA Tesla K80 general purpose GPUs. We demonstrate the utility of this platform for large radiotherapy optimization problems through the development and characterization of a parallelized particle swarm optimization (PSO) four dimensional (4D) intensity modulated radiation therapy (IMRT) technique. The PSO engine was coupled to the Eclipse treatment planning system via a vendor-provided scripting interface. Specific challenges addressed in this implementation were (i) data management and (ii) non-uniform memory access (NUMA). For the former, we alternated between parameters over which the computation process was parallelized. For the latter, we reduced the amount of data required to be transferred over the NUMA bridge. The datasets examined in this study were approximately 300 GB in size, including 4D computed tomography images, anatomical structure contours and dose deposition matrices. For evaluation, we created a 4D-IMRT treatment plan for one lung cancer patient and analyzed computation speed while varying several parameters (number of respiratory phases, GPUs, PSO particles, and data matrix sizes). The optimized 4D-IMRT plan enhanced sparing of organs at risk by an average reduction of [Formula: see text] in maximum dose, compared to the clinical optimized IMRT plan, where the internal target volume was used. We validated our computation time analyses in two additional cases. The computation speed in our implementation did not monotonically increase with the number of GPUs. The optimal number of GPUs (five, in our study) is directly related to the hardware specifications. The optimization process took 35 min using 50 PSO particles, 25 iterations and 5 GPUs.
Assuntos
Tomografia Computadorizada Quadridimensional/instrumentação , Tomografia Computadorizada Quadridimensional/métodos , Neoplasias Pulmonares/radioterapia , Órgãos em Risco/efeitos da radiação , Planejamento da Radioterapia Assistida por Computador/métodos , Radioterapia de Intensidade Modulada/métodos , Humanos , Neoplasias Pulmonares/diagnóstico por imagem , Dosagem Radioterapêutica , Estudos RetrospectivosRESUMO
PURPOSE: We present a particle swarm optimization (PSO)-based technique to create deliverable four-dimensional (4D = 3D + time) intensity-modulated radiation therapy (IMRT) plans for lung stereotactic body radiotherapy (SBRT). The 4D planning concept uses respiratory motion as an additional degree of freedom to achieve further sparing of organs at risk (OARs). The 4D-IMRT plan involves the delivery of an order of magnitude more IMRT apertures (~15,000-20,000), with potentially large interaperture variations in the delivered fluence, compared to conventional (i.e., 3D) IMRT. In order to deliver the 4D plan in an efficient manner, we present an optimization-based aperture sequencing technique. METHOD: A graphic processing unit (GPU)-enabled PSO-based inverse planning engine, developed and integrated with a research version of the Eclipse (Varian, Palo Alto, CA) treatment planning system (TPS), was employed to create 4D-IMRT plans as follows. Four-dimensional computed tomography scans (4DCTs) and beam configurations from clinical treatment plans of seven lung cancer patients were retrospectively collected, and in each case, the PSO engine iteratively adjusted aperture monitor unit (MU) weights for all beam apertures across all respiratory phases to optimize OAR dose sparing while maintaining planning target volume (PTV) coverage. We calculated the transition times from each aperture to all other apertures for each beam, taking into account the maximum leaf velocity of the multileaf collimator (MLC), and developed a mixed integer optimization technique for aperture sequencing. The goal of sequencing was to maximize delivery efficiency (i.e., minimize the time required to deliver the dose map) by accounting for leaf velocity, aperture MUs, and duration of each respiratory phase. The efficiency of the proposed delivery method was compared with that of a greedy algorithm which chose only from neighboring apertures for the subsequent steps in the sequence. RESULTS: 4D-IMRT-optimized plans achieved PTV coverage comparable to clinical plans while improving OAR sparing by an average of 39.7% for D max heart, 20.5% for D max esophagus, 25.6% for D max spinal cord, and 2.1% for V 13 lung (with D max standing for maximum dose and V 13 standing for volume receiving ≥ 13 Gy). Our mixed integer optimization-based aperture sequencing enabled the delivery to be performed in fewer cycles compared to the greedy method. This reduction was 89 ± 79 cycles corresponding to an improvement of 15.94 ± 8.01%, when considering respiratory cycle duration of 4 s, and 55 ± 33 cycles corresponding to an improvement of 15.14 ± 4.45%, when considering respiratory cycle duration of 6 s. CONCLUSION: PSO-based 4D-IMRT represents an attractive technique to further improve OAR sparing in lung SBRT. Efficient delivery of a large number of sparse apertures (control points) introduces a challenge in 4D-IMRT treatment planning and delivery. Through judicious optimization of the aperture sequence across all phases, such delivery can be performed on a clinically feasible time scale.
Assuntos
Tomografia Computadorizada Quadridimensional , Neoplasias Pulmonares/diagnóstico por imagem , Neoplasias Pulmonares/radioterapia , Radiocirurgia , Planejamento da Radioterapia Assistida por Computador/métodos , Humanos , Neoplasias Pulmonares/fisiopatologia , Movimento , Dosagem Radioterapêutica , Respiração , Estudos RetrospectivosRESUMO
OBJECTIVE: Evolutionary stochastic global optimization algorithms are widely used in large-scale, nonconvex problems. However, enhancing the search efficiency and repeatability of these techniques often requires well-customized approaches. This study investigates one such approach. METHODS: We use particle swarm optimization (PSO) algorithm to solve a 4D radiation therapy (RT) inverse planning problem, where the key idea is to use respiratory motion as an additional degree of freedom in lung cancer RT. The primary goal is to administer a lethal dose to the tumor target while sparing surrounding healthy tissue. Our optimization iteratively adjusts radiation fluence-weights for all beam apertures across all respiratory phases. We implement three PSO-based approaches: conventionally used unconstrained, hard-constrained, and our proposed virtual search. As proof of concept, five lung cancer patient cases are optimized over ten runs using each PSO approach. For comparison, a dynamically penalized likelihood (DPL) algorithm-a popular RT optimization technique is also implemented and used. RESULTS: The proposed technique significantly improves the robustness to random initialization while requiring fewer iteration cycles to converge across all cases. DPL manages to find the global optimum in 2 out of 5 RT cases over significantly more iterations. CONCLUSION: The proposed virtual search approach boosts the swarm search efficiency, and consequently, improves the optimization convergence rate and robustness for PSO. SIGNIFICANCE: RT planning is a large-scale, nonconvex optimization problem, where finding optimal solutions in a clinically practical time is critical. Our proposed approach can potentially improve the optimization efficiency in similar time-sensitive problems.