Determining clinical photon beam spectra from measured depth dose with the Cimmino algorithm

Comparison between transmission and scattering spectrum reconstruction methods based on EPID images

Numerous improved physics-based methods for Linac photon spectra reconstruction have been published; some of them are based on transmission data analysis and others on scattering data. In this work, the two spectrum unfolding approaches are compared in order to experimentally validate its robustness and to determine which is the optimal methodology for application on a clinical quality assurance routine. Both studied methods are based on EPID images generated when the incident photon beam impinges onto plastic blocks. The distribution of transmitted/scatter radiation produced by this object centered at the beam field size was measured. Measurements were performed using a 6 MeV photon beam produced by the linear accelerator. The same radiation distribution conditions were also simulated with Monte Carlo code for a series of monoenergetic identical geometry photon beams for both cases. Two systems of linear equations were generated to combine the polyenergetic EPID measurements with the monoenergetic simulation results. Regularization techniques were applied to solve the systems for obtaining the incident photon spectrum. We present a comparison between the well-known photon Spectral Reconstruction based on Transmission Data (Trans-based) technology and the Spectral Reconstruction based on Scattering Data (Scatt-based), which we both developed using EPID images. It is shown that Trans-based reconstruction results display much better agreement with photon spectrum theoretical predictions.

Enhanced bremsstrahlung spectrum reconstruction from depth–dose gradients

We present a technique to significantly improve the accuracy of bremsstrahlung spectra reconstructed from central-axis depth-dose data using inverse methods. While typical approaches directly use the measured depth-dose data, we show the advantage of using the gradient of these data for reconstruction. The inverse problem in terms of gradients is shown to be markedly less ill-conditioned than the usual inverse problem. In each case, a regularization is introduced to alleviate the effects of noise due to measurement and computation. The error in taking derivatives of depth-dose data is demonstrated to be sufficiently controllable as not to overtake the improvements in the conditioning of the inverse problem in a pair of simulated examples.

Enhanced spectral discrimination through the exploitation of interface effects in photon dose data

The convolution/superposition algorithm for computing dose from photon beams in radiation therapy planning requires knowledge of the energy spectrum. The algorithm can compute the dose for a polyenergetic beam as the weighted sum of the individual dose contributions from monoenergetic beams. In this study we exploit interface effects apparent in the dose distributions to discriminate among spectra of high energy photon beams. We have studied the sensitivity of the depth dose distribution to the energy components using a hypothetical beam for various field sizes and depths in water and water-lung-water media. Six theoretical spectra were simulated. We compared depth dose data from these spectra using three quantitative measures which are inherently free of normalization ambiguities: for homogeneous water, the ratio D20/D10 and a logarithmic derivative in the buildup region LD(build-up) and for inhomogeneous lung/water, the lung correction factor (CF). It was found that the ability of both the CF and the LD(build-up) tests to discriminate between the various theoretical spectra were superior to that of the D20/D10 test. This discriminating power of the CF test decreases with increasing field size due to restored electronic equilibrium. The CF test, though, has some advantages over the LD(build-up) test since it is less prone to electron contamination issues and numerical errors. A practical example with a 15 MV photon beam illustrates the process. Consequently, we suggest that as part of a beam-commissioning methodology, designated electronic disequilibrium test cases be implemented in unambiguously determining the correct energy spectrum to be used.

Fast treatment plan modification with an over-relaxed Cimmino algorithm

A method to quickly modify a treatment plan in adaptive radiotherapy was proposed and studied. The method is based on a Cimmino-type algorithm in linear programming. The fast convergence speed is achieved by over-relaxing the algorithm relaxation parameter from its sufficient convergence range of (0, 2) to (0, infinity). The algorithm parameters are selected so that the over-relaxed Cimmino (ORC) algorithm can effectively approximate an unconstrained re-optimization process in adaptive radiotherapy. To demonstrate the effectiveness and flexibility of the proposed method in adaptive radiotherapy, two scenarios with different organ motion/deformation of one nasopharyngeal case were presented with comparisons made between this method and the re-optimization method. In both scenarios, the ORC algorithm modified treatment plans have dose distributions that are similar to those given by the re-optimized treatment plans. It takes us using the ORC algorithm to finish a treatment plan modification at least three times faster than the re-optimization procedure compared.