Leaf sequencing via difference matrices reduces aperture number and beam-on time in IMRT
Robert R. Meyer,Athula Gunawardena,Warren D. D'Souza,Wensha Yang,Shahid A. Naqvi,Leyuan Shi
DOI: https://doi.org/10.1016/s0360-3016(03)01109-x
2003-01-01
Abstract:Purpose/Objective: Leaf-sequencing in the Corvus (Nomos, PA) planning system takes place after the optimization of the treatment plan to be delivered in step-and-shoot IMRT. The Corvus leaf-sequencing algorithm generates a large number of apertures, which can result in significant treatment delivery time. We present a new approach for the intensity map decomposition problem that is solved to obtain leaf sequencing.Materials/Methods: Mathematically, the version of the problem that we consider may be modeled as the optimal decomposition of a matrix M of nonnegative integer intensities into a weighted sum: M = σ ci Si where the intensities ci and the binary matrices Si are to be determined by the algorithm so as to make the number of matrices (segments) and sum of the intensities σ ci in the decomposition as small as possible. For certain linacs (e.g., Elekta), the transition time between segments can dominate the actual radiation delivery time. In order to define the difference matrices used in our approach, we assume that we are given an m × n intensity map M = (Mij) and we add bordering columns on the left and right via ghost entries Mi,-1 = Mi,n+1 = 0 for i=1,…,m. After this augmentation, the difference operator Δ and initial difference matrix D are defined as follows: D=ΔM=(Dij) is the m × (n+1) matrix such that Dij=Mi,j-Mi,j-1. The algorithm generates a finite sequence of intensity maps and corresponding difference matrices, at each stage extracting from the current intensity map a shape matrix and its corresponding intensity multiplier, followed by a corresponding update of the intensity map and difference matrices. The extraction process is guided by examination of the current difference matrix. The chosen extraction is the one that zeros out the maximum number of entries in the difference matrix. This approach was applied to clinical test cases: head and neck and prostate. Plans were generated with 5, 10 and 100 intensity levels. The number of apertures (segments) and monitor units (beam-on time) resulting from our approach was compared with the output from Corvus and a previously published method by Boland et al (BHL).Results: A 5-field and 6-field plan was generated for the head and neck and prostate cases, respectively. The prescription doses were 45 Gy and 75.6 Gy for the head and neck and prostate cases, respectively. The reduction in the number of segments with our approach was 72%, 67% and 32% with 5, 10 and 100 intensity levels, respectively compared with Corvus for the head and neck case. The corresponding reduction compared to the BHL method was 74%, 62% and 13%, respectively. For the prostate case, the reduction in segments was 65%, 62% and 21% with 5, 10 and 100 intensity levels, respectively between our method and Corvus. The corresponding reduction in comparison with the BHL method was 60%, 51% and 6%. The number of monitor units is reduced by > 20% compared with the Corvus plans. We have processed the aperture shapes from our approach and the Corvus leaf sequence through a Monte Carlo calculation. The intensity profiles outside the beam’s-eye-view (BEV) projection of the tumor display leakage radiation with the Corvus leaf sequence which, can be attributed to the use of the partial transmission block used by Corvus in conjunction with Elekta linear accelerators. The leakage radiation outside the BEV from our approach is significantly reduced. However, tongue-and-groove effects are more prominent in our approach. The impact of tongue-and-groove has been studied by calculating a full 3-D dose distribution of the Corvus leaf sequence and the output from our method. Differences in the dose distribution resulting from our sequencing approach and the Corvus plan are negligible.Conclusions: Our results demonstrate that we have developed a fast heuristic for leaf sequencing which reduces the number of aperture shapes and beam-on time significantly compared with the output from Corvus and existing methods. We estimate that the overall delivery time will be 15 to 20 minutes with 5 to 10 intensity levels, a significant reduction compared to current delivery times. Purpose/Objective: Leaf-sequencing in the Corvus (Nomos, PA) planning system takes place after the optimization of the treatment plan to be delivered in step-and-shoot IMRT. The Corvus leaf-sequencing algorithm generates a large number of apertures, which can result in significant treatment delivery time. We present a new approach for the intensity map decomposition problem that is solved to obtain leaf sequencing. Materials/Methods: Mathematically, the version of the problem that we consider may be modeled as the optimal decomposition of a matrix M of nonnegative integer intensities into a weighted sum: M = σ ci Si where the intensities ci and the binary matrices Si are to be determined by the algorithm so as to make the number of matrices (segments) and sum of the intensities σ ci in the decomposition as small as possible. For certain linacs (e.g., Elekta), the transition time between segments can dominate the actual radiation delivery time. In order to define the difference matrices used in our approach, we assume that we are given an m × n intensity map M = (Mij) and we add bordering columns on the left and right via ghost entries Mi,-1 = Mi,n+1 = 0 for i=1,…,m. After this augmentation, the difference operator Δ and initial difference matrix D are defined as follows: D=ΔM=(Dij) is the m × (n+1) matrix such that Dij=Mi,j-Mi,j-1. The algorithm generates a finite sequence of intensity maps and corresponding difference matrices, at each stage extracting from the current intensity map a shape matrix and its corresponding intensity multiplier, followed by a corresponding update of the intensity map and difference matrices. The extraction process is guided by examination of the current difference matrix. The chosen extraction is the one that zeros out the maximum number of entries in the difference matrix. This approach was applied to clinical test cases: head and neck and prostate. Plans were generated with 5, 10 and 100 intensity levels. The number of apertures (segments) and monitor units (beam-on time) resulting from our approach was compared with the output from Corvus and a previously published method by Boland et al (BHL). Results: A 5-field and 6-field plan was generated for the head and neck and prostate cases, respectively. The prescription doses were 45 Gy and 75.6 Gy for the head and neck and prostate cases, respectively. The reduction in the number of segments with our approach was 72%, 67% and 32% with 5, 10 and 100 intensity levels, respectively compared with Corvus for the head and neck case. The corresponding reduction compared to the BHL method was 74%, 62% and 13%, respectively. For the prostate case, the reduction in segments was 65%, 62% and 21% with 5, 10 and 100 intensity levels, respectively between our method and Corvus. The corresponding reduction in comparison with the BHL method was 60%, 51% and 6%. The number of monitor units is reduced by > 20% compared with the Corvus plans. We have processed the aperture shapes from our approach and the Corvus leaf sequence through a Monte Carlo calculation. The intensity profiles outside the beam’s-eye-view (BEV) projection of the tumor display leakage radiation with the Corvus leaf sequence which, can be attributed to the use of the partial transmission block used by Corvus in conjunction with Elekta linear accelerators. The leakage radiation outside the BEV from our approach is significantly reduced. However, tongue-and-groove effects are more prominent in our approach. The impact of tongue-and-groove has been studied by calculating a full 3-D dose distribution of the Corvus leaf sequence and the output from our method. Differences in the dose distribution resulting from our sequencing approach and the Corvus plan are negligible. Conclusions: Our results demonstrate that we have developed a fast heuristic for leaf sequencing which reduces the number of aperture shapes and beam-on time significantly compared with the output from Corvus and existing methods. We estimate that the overall delivery time will be 15 to 20 minutes with 5 to 10 intensity levels, a significant reduction compared to current delivery times.