PubMed Health⌕ Search

Biomedical subjects

B M McCurdy

Publications and source records attributed to B M McCurdy.

5 recordsLinked to original sources

Dosimetric investigation and portal dose image prediction using an amorphous silicon electronic portal imaging device.

A two step algorithm to predict portal dose images in arbitrary detector systems has been developed recently. The current work provides a validation of this algorithm on a clinically available, amorphous silicon flat panel imager. The high-atomic number, indirect amorphous silicon detector incorporates a gadolinium oxysulfide phosphor scintillating screen to convert deposited radiation energy to optical photons which form the portal image. A water equivalent solid slab phantom and an anthropomorphic phantom were examined at beam energies of 6 and 18 MV and over a range of air gaps (approximately 20-50 cm). In the many examples presented here, portal dose images in the phosphor were predicted to within 5% in low-dose gradient regions, and to within 5 mm (isodose line shift) in high-dose gradient regions. Other basic dosimetric characteristics of the amorphous silicon detector were investigated, such as linearity with dose rate (+/- 0.5%), repeatability (+/- 2%), and response with variations in gantry rotation and source to detector distance. The latter investigation revealed a significant contribution to the image from optical photon spread in the phosphor layer of the detector. This phenomenon is generally known as "glare," and has been characterized and modeled here as a radially symmetric blurring kernel. This kernel is applied to the calculated dose images as a convolution, and is successfully demonstrated to account for the optical photon spread. This work demonstrates the flexibility and accuracy of the two step algorithm for a high-atomic number detector. The algorithm may be applied to improve performance of dosimetric treatment verification applications, such as direct image comparison, backprojected patient dose calculation, and scatter correction in megavoltage computed tomography. The algorithm allows for dosimetric applications of the new, flat panel portal imager technology in the indirect configuration, taking advantage of a greater than tenfold increase in detector sensitivity over a direct configuration.

Algorithms↗

A two-step algorithm for predicting portal dose images in arbitrary detectors.

Recently, portal imaging systems have been successfully demonstrated in dosimetric treatment verification applications, where measured and predicted images are quantitatively compared. To advance this approach to dosimetric verification, a two-step model which predicts dose deposition in arbitrary portal image detectors is presented. The algorithm requires patient CT data, source-detector distance, and knowledge of the incident beam fluence. The first step predicts the fluence entering a portal imaging detector located behind the patient. Primary fluence is obtained through ray-tracing techniques, while scatter fluence prediction requires a library of Monte Carlo-generated scatter fluence kernels. These kernels allow prediction of basic radiation transport parameters characterizing the scattered photons, including fluence and mean energy. The second step of the algorithm involves a superposition of Monte Carlo-generated pencil beam kernels, describing dose deposition in a specific detector, with the predicted incident fluence. This process is performed separately for primary and scatter fluence, and yields a predicted dose image. A small but noticeable improvement in prediction is obtained by explicitly modeling the off-axis energy spectrum softening due to the flattening filter. The algorithm is tested on a slab phantom and a simple lung phantom (6 MV). Furthermore, an anthropomorphic phantom is utilized for a simulated lung treatment (6 MV), and simulated pelvis treatment (23 MV). Data were collected over a range of air gaps (10-80 cm). Detectors incorporating both low and high atomic number buildup are used to measure portal image profiles. Agreement between predicted and measured portal dose is better than 3% in areas of low dose gradient (<30%/cm) for all phantoms, air gaps, beam energies, and detector configurations tested here. It is concluded that this portal dose prediction algorithm is fast, accurate, allows separation of primary and scatter dose, and can model arbitrary detectors.

Air↗

Photon scatter in portal images: physical characteristics of pencil beam kernels generated using the EGS Monte Carlo code.

Pencil beam kernels describing scattered photon fluence behind homogeneous water slabs at various air gap distances were generated using the EGS Monte Carlo code. Photon scatter fluence was scored in separate bins based on the particle's history: singly scattered, multiply scattered, and bremsstrahlung and positron annihilation photons. Simultaneously, the mean energy and mean angle with respect to the incident photon pencil beam were tallied. Kernels were generated for incident photon pencil beams exhibiting monoenergetic spectra of 2.0 and 10.0 MeV, and polyenergetic spectra representative of 6 and 24 MV beams. Reciprocity was used to generate scatter fractions on the central axis for various field sizes, phantom thicknesses, and air gaps. The scatter kernels were further characterized by full width at half-maximum estimates. Modulation transfer functions were calculated, providing theoretical estimates of the limit of performance of portal imaging systems due to the intrinsic scattering of photon radiation through the patient.

Humans↗

Photon scatter in portal images: accuracy of a fluence based pencil beam superposition algorithm.

The accuracy of a pencil beam algorithm to predict scattered photon fluence into portal imaging systems was studied. A data base of pencil beam kernels describing scattered photon fluence behind homogeneous water slabs (1-50 cm thick) at various air gap distances (0-100 cm) was generated using the EGS Monte Carlo code. Scatter kernels were partitioned according to particle history: singly-scattered, multiply-scattered, and bremsstrahlung and positron annihilation photons. Mean energy and mean angle with respect to the incident photon pencil beam were also scored. This data allows fluence, mean energy, and mean angular data for each history type to be predicted using the pencil beam algorithm. Pencil beam algorithm predictions for 6 and 24 MV incident photon beams were compared against full Monte Carlo simulations for several inhomogeneous phantoms, including approximations to a lateral neck, and a mediastinum treatment. The accuracy of predicted scattered photon fluence, mean energy, and mean angle was investigated as a function of air gap, field size, photon history, incident beam resolution, and phantom geometry. Maximum errors in mean energies were 0.65 and 0.25 MeV for the higher and lower energy spectra, respectively, and 15 degrees for mean angles. The ability of the pencil beam algorithm to predict scatter fluence decreases with decreasing air gap, with the largest error for each phantom occurring at the exit surface. The maximum predictive error was found to be 6.9% with respect to the total fluence on the central axis. By maintaining even a small air gap (approximately 10 cm), the error in predicted scatter fluence may be kept under 3% for the phantoms and beam energies studied here. It is concluded that this pencil beam algorithm is sufficiently accurate (using International Commission on Radiation Units and Measurements Report No. 24 guidelines for absorbed dose) over the majority of clinically relevant air gaps, for further investigation in a portal dose prediction algorithm.

Air↗

Determination of equivalent photon fields through integrated 1D convolution kernels.

The equivalent fields concept has been studied using a variation of the convolution dose calculation method. This allows a better understanding of this classic, widely-used concept and its relation to modern dose calculation techniques. Total scatter energy contribution as a function of field size was computed by radially integrating three-dimensional point interaction scatter kernels representing cobalt-60, 6 MV and 24 MV photon beam spectra and utilizing reciprocity. For arbitrarily shaped treatment fields, equivalent square or circular fields are chosen based on a conservation of energy approach. Associated depth dose curves are generated via convolution of primary fluence with pre-integrated scatter kernels. For the energy spectra used, the resulting equivalent circle to square relationships agree with the empirical equation recommended in the British Journal of Radiology Supplement 25 (BJR 25) (to within 2%). Tables converting rectangular fields to equivalent squares and circles have been generated and most (approximately 63%) of the equivalent square (and circular) field sizes are within 0.5 cm of the BJR 25 tables. This agreement provides an explanation for the successful use of a single relationship for all therapeutic energies. Fields of different shape will be approximately equivalent only when the scatter kernels are similar in shape and total energy content. The equivalent field concept inevitably results in an optimal match of depth dose curves at a single depth only. These ideas are demonstrated for a highly elongated (2 x 30 cm2) field and its equivalent square field. The invalidity of the equivalent fields concept off-axis is demonstrated by examining predicted dose deposition for several points in and around a 15 x 15 cm2 cobalt-60 beam. Modern convolution techniques are shown to offer a unique approach to the classic equivalent fields method. The appropriate use of the BJR 25 equivalent fields tables is recommended for all photon energies.

Monte Carlo Method↗