Verification of field match lines in whole breast radiation therapy using Cherenkov imaging

PURPOSE

In mono-isocentric radiation therapy treatment plans designed to treat the whole breast and supraclavicular lymph nodes, the fields meet at isocenter, forming the match line. Insufficient coverage at the match line can lead to recurrence, and overlap over weeks of treatment can lead to increased risk of healthy tissue toxicity. Cherenkov imaging was used to assess the accuracy of delivery at the match line and identify potential incidents during patient treatments.

METHODS AND MATERIALS

A controlled calibration was constructed from the deconvolved Cherenkov images from the delivery of a modified patient treatment plan to an anthropomorphic phantom with introduced separation and overlap. The trend from this calibration was then used to evaluate the field match line for accuracy and inter-fraction consistency for two patients.

RESULTS

The intersection point between matching field profiles was directly correlated to the distance (gap/overlap) between the fields (anthropomorphic phantom R² = 0.994 "breath hold" and R² = 0.990 "free breathing"). The profile intersection points from two patients' imaging sessions yielded an average of +1.40 mm offset (overlap) and -1.32 mm offset (gap), thereby introducing roughly a 25.0% over-dose and a -23.6% under-dose (R² = 0.994).

CONCLUSIONS

This study shows that field match regions can be detected and quantified by taking deconvolved Cherenkov images and using their product image to create steep intensity gradients, causing match lines to stand out. These regions can then be quantitatively translated into a dose consequence. This approach offers a high sensitivity detection method which can quantify match line variability and errors in vivo.

Analysis of corrected Cerenkov emission during electron radiotherapy by Monte Carlo method

Cerenkov emission during electron radiotherapy had been emerging as a new dose assessment approach for clinical radiotherapy and could be imaged through a standard commercial camera. The purpose of this work aimed to study the accuracy of corrected Cerenkov emission method during electron radiotherapy by Monte Carlo (MC) method. GAMOS MC software was used to model the physics of electron therapy and calculated dose and Cerenkov photon distribution in water phantom. Compared to ionization chamber and diode measurement, MC simulated dose discrepancy was less than 1% in percentage depth dose (PDD) curves and less than. 2% in crossline profile curves, which was acceptable for clinical criterion. Compared to ionization chamber dose measurement, MC simulated Cerenkov discrepancy was less than 2% in crossline profile distribution, which was acceptable for clinical criterion. However, the Cerenkov PDD curves tended to overestimate the dose at the build-up region and underestimate the dose at the remaining attenuation region. After increasing the Cerenkov distribution depth to 2-3 mm, the discrepancy became well within 1% at the remaining attenuation region, which was acceptable for clinical criterion. Therefore, corrected Cerenkov emission could be used to assess PDD accuracy and crossline profile accuracy during electron radiotherapy.

Space-variant deconvolution of Cerenkov light images acquired from a curved surface

PURPOSE

Cerenkov photons are generated by high-energy radiation used in external beam radiation therapy (EBRT). This study expands upon the Cerenkov light dosimetry formula previously developed to relate an image of Cerenkov photons to the primary beam fluence. Extension of this formulation allows for deconvolution to be performed on images acquired from curved geometries.

METHODS

The integral equation, which represented the formation of Cerenkov photon image from an incident high-energy photon beam, was expanded to allow for space-variance of the convolution kernel called as the Cerenkov scatter function (CSF). The GAMOS (Geant4-based Architecture for Medicine-Oriented Simulations) Monte Carlo (MC) particle simulation software was used to obtain the CSF for different incident beam angles. The image of a curved surface was first projected to a flat plane by using a perspective correction method. Then, the planar image was partitioned into small segments (or blocks), where a CSF corresponding to a specific beam incident angle was applied for deconvolution. The block size and the margin around the block were optimized by studying the effects of those parameters on the deconvolution accuracy for a test image. We evaluated three deconvolution techniques: Richardson-Lucy, Blind, and Total Variation minimization (TV/L2) algorithms, to select the most accurate method for the current applications.

RESULTS

Analysis of deconvolution algorithms showed that the TV/L2 method provided the most accurate solution to the deconvolution problem for Cerenkov imaging. Optimization of space-variant deconvolution parameters showed that including a margin that is at least 42.9% of the image width provided the most accurate product image. There was no optimal size for the deconvolution area and should be chosen based on the presence of unique CSF kernels within an image. Space-variant deconvolution improved measured field size in Cerenkov photon images by 7.37%, as compared with 1.74% by space-invariant deconvolution. Space-variant deconvolution improved measured penumbra by 99.3%, as compared with 76.7% by space-invariant deconvolution. Space-variant deconvolution introduced artifacts in flat regions of the beam. Artifacts were avoided through selective space-variant deconvolution in only the penumbra region.

CONCLUSIONS

Primary photon fluence distributions of a curved surface can be obtained by using space-variant deconvolution methods in Cerenkov light dosimetry. The TV/L2 algorithm is the best method for deconvolution of Cerenkov photon images from an open-field beam derived from either a flat or curved surface. The partition size chosen for space-variant deconvolution should be at least six times the full width at half maximum (FWHM) of the corresponding scatter kernel used in deconvolution. Space-variant deconvolution is necessary if the incident beam angle difference is larger than 6 ∘ between regions of an image.

Characterization of the Cerenkov scatter function: a convolution kernel for Cerenkov light dosimetry

Cerenkov light is created in clinical applications involving high-energy radiation such as in radiation therapy. There is considerable interest in using Cerenkov light as a means to perform in vivo dosimetry during radiation therapy; however, a better understanding of the light-to-dose relationship is needed. One such method to solve this relationship is that of a deconvolution formulation, which relies on the Cerenkov scatter function (CSF). The CSF describes the creation of Cerenkov photons by a pencil beam of high-energy radiation, and the subsequent scattering that occurs before emission from the irradiated medium surface. This study investigated the dependence of the CSF on common radiation beam parameters (beam energy and incident angle) and the type of irradiated medium. An analytical equation with fitting coefficients of the CSF was obtained for common beam energies in a stratified skin model and optical phantom. Perturbation analysis was performed to investigate the dependence of the deconvolved Cerenkov images on the full-width at half-maximum and amplitude of the CSF. The irradiated material and beam angle had a large impact on the deconvolution process, whereas the beam energy had little effect.

A mathematical deconvolution formulation for superficial dose distribution measurement by Cerenkov light dosimetry

PURPOSE

Cerenkov photons are created by high-energy radiation beams used for radiation therapy. In this study, we developed a Cerenkov light dosimetry technique to obtain a two-dimensional dose distribution in a superficial region of medium from the images of Cerenkov photons by using a deconvolution method.

METHODS

An integral equation was derived to represent the Cerenkov photon image acquired by a camera for a given incident high-energy photon beam by using convolution kernels. Subsequently, an equation relating the planar dose at a depth to a Cerenkov photon image using the well-known relationship between the incident beam fluence and the dose distribution in a medium was obtained. The final equation contained a convolution kernel called the Cerenkov dose scatter function (CDSF). The CDSF function was obtained by deconvolving the Cerenkov scatter function (CSF) with the dose scatter function (DSF). The GAMOS (Geant4-based Architecture for Medicine-Oriented Simulations) Monte Carlo particle simulation software was used to obtain the CSF and DSF. The dose distribution was calculated from the Cerenkov photon intensity data using an iterative deconvolution method with the CDSF. The theoretical formulation was experimentally evaluated by using an optical phantom irradiated by high-energy photon beams.

RESULTS

The intensity of the deconvolved Cerenkov photon image showed linear dependence on the dose rate and the photon beam energy. The relative intensity showed a field size dependence similar to the beam output factor. Deconvolved Cerenkov images showed improvement in dose profiles compared with the raw image data. In particular, the deconvolution significantly improved the agreement in the high dose gradient region, such as in the penumbra. Deconvolution with a single iteration was found to provide the most accurate solution of the dose. Two-dimensional dose distributions of the deconvolved Cerenkov images agreed well with the reference distributions for both square fields and a multileaf collimator (MLC) defined, irregularly shaped field.

CONCLUSIONS

The proposed technique improved the accuracy of the Cerenkov photon dosimetry in the penumbra region. The results of this study showed initial validation of the deconvolution method for beam profile measurements in a homogeneous media. The new formulation accounted for the physical processes of Cerenkov photon transport in the medium more accurately than previously published methods.