Evaluation of electronic portal imaging device for missing tissue compensator design and verification

Correction for dose-response variations in a scanning liquid ion chamber EPID as a function of linac gantry angle

The response of electronic portal imaging devices (EPIDs) of the scanning liquid ionization chamber (SLIC) type is known to vary with linear accelerator gantry angle. This work considered several contributing factors, quantified the artefacts, monitored their reproducibility and investigated the effects of repeated gantry rotations. Unflatness of up to 5% was found. A correction technique was devised using nonlinear regression of a three-variable sinusoidal modulation. Comparison with two existing techniques found our method to be the most effective, providing a flatness well within 2%. This improved accuracy is expected to benefit more accurate dosimetric studies in particular. The post-acquisition correction process required no change in imaging protocols. Applicability of the new technique was demonstrated on images acquired on different days and with different beam sizes. Since the artefacts compromise both accurate dosimetry and image quality, their successful removal should benefit a broad range of SLIC EPID applications.

Monte Carlo simulation of portal dosimetry on a rectilinear voxel geometry: a variable gantry angle solution

A software solution has been developed to carry out Monte Carlo simulations of portal dosimetry using the BEAMnrc/DOSXYZnrc code at oblique gantry angles. The solution is based on an integrated phantom, whereby the effect of incident beam obliquity was included using geometric transformations. Geometric transformations are accurate within +/- 1 mm and +/- 1 degrees with respect to exact values calculated using trigonometry. An application in portal image prediction of an inhomogeneous phantom demonstrated good agreement with measured data, where the root-mean-square of the difference was under 2% within the field. Thus, we achieved a dose model framework capable of handling arbitrary gantry angles, voxel-by-voxel phantom description and realistic particle transport throughout the geometry.

Compensator quality control with an amorphous silicon EPID

The calibration and quality control of compensators is conventionally performed with an ion chamber in a water-equivalent phantom. In our center, the compensator factor and four off-axis fluence ratios are measured to verify the central axis beam modulation and orientation of the compensator. Here we report the investigation of an alternative technique for compensator quality control using an amorphous silicon electronic portal imaging device (a-Si EPID). Preliminary experiments were performed to identify appropriate EPID operating parameters for this relative dosimetric study and also to quantify EPID operation. The pixel value versus energy fluence response of the EPID for both open and compensated fields was then determined, and expressed via calibration curves. For open fields the response was seen to be linear, whereas for compensated fields it exhibited a small quadratic component. To account for field size effects, we measured EPID scatter factors. These exhibited small but non-negligible dependencies on compensator thickness and source-detector distance. Finally, a number of test and clinical compensators were evaluated to assess the suitability of the EPID for compensator quality control. Our results indicate that the a-Si EPID can measure clinical compensator factors and off-axis energy fluence ratios to within 2% of values measured by a Farmer chamber on average, and so is a suitable ion chamber replacement.

An investigation of a new amorphous silicon electronic portal imaging device for transit dosimetry

The relationship between the pixel value and exit dose was investigated for a new commercially available amorphous silicon electronic portal imaging device. The pixel to dose mapping function was established to be linear for detector distances between 116.5 cm to 150 cm from the source, radiation field sizes from 5 x 5 cm2 to 20 x 20 cm2 and beam energies of 6 to 18 MV. Coefficients in the mapping function were found to be dependent on beam energy and field size. Open and wedged field profiles measured with the device showed agreement to a maximum of 5% and 8%, respectively, as compared to film. A comparison of relative transmission measurements between the EPID and ion chamber indicate a maximum deviation of 6% and 2% at 6 and 18 MV, respectively, for an attenuator thickness of 21 cm and SDD > or = 130 cm. It was found that accuracies of better than 1% could be obtained if detector position and field size specific fitting parameters were used to generate unique mapping functions for each configuration.

Assessment of flatness and symmetry of megavoltage x-ray beam with an electronic portal imaging device (EPID)

The input/output characteristics of the Wellhofer BIS 710 electronic portal imaging device (EPID) have been investigated to establish its efficacy for periodic quality assurance (QA) applications. Calibration curves have been determined for the energy fluence incident on the detector versus the pixel values. The effect of the charge coupled device (CCD) camera sampling time and beam parameters (such as beam field size, dose rate, photon energy) on the calibration have been investigated for a region of interest (ROI) around the central beam axis. The results demonstrate that the pixel output is a linear function of the incident exposure, as expected for a video-based electronic portal imaging system. The field size effects of the BIS 710 are similar to that of an ion chamber for smaller field sizes up to 10 x 10 cm2. However, for larger field sizes the pixel value increases more rapidly. Furthermore, the system is slightly sensitive to dose rate and is also energy dependent The BIS 710 has been used in the current study to develop a QA procedure for measurements of flatness and symmetry of a linac x-ray beam. As a two-dimensional image of the radiation field is obtained from a single exposure of the BIS 710, a technique has been developed to calculate flatness and symmetry from a defined radiation area. The flatness and symmetry values obtained are different from those calculated conventionally from major axes only (inplane, crossplane). This demonstrates that the technique can pick up the "cold" and "hot" spots in the analysed area, providing thus more information about the radiation beam. When calibrated against the water tank measurements, the BIS 710 can be used as a secondary device to monitor the x-ray beam flatness and symmetry.

Dosimetric properties of the Theraview fluoroscopic electronic portal imaging device

Electronic portal imaging devices (EPIDs) can be used for non-imaging applications in radiotherapy such as patient dosimetry. Of the systems available, the fluoroscopic camera-based EPID Theraview (InfiMed Inc.) has not been studied to date, and a review of the dosimetric properties of the system is presented here. In the "single set-up" mode of image acquisition, pixel intensity increases sublinearly with applied dose. The response was dependent on the system's video signal gain and showed a threshold dose to the detector in the range 0.05-0.35 cGy, and pixel saturation at detector doses in the range 1.2-1.6 cGy. Repeated exposures of the EPID were observed to be extremely reproducible (standard deviation 0.5%). The sensitivity of the system showed a linear decline of 0.04% day-1 over a 68-day period, during which time the relative off-axis response within 10 x 10 cm2 field was constant to within a standard deviation of 0.56%. The system shows spatial non-uniformity, which requires correction for application to dose measurements in two-dimensions. Warm-up of the camera control unit required a period of at least 40 min and was associated with an enhancement in pixel intensity of up to 12%. A radiation dose history effect was observed at doses as low as 0.2 Gy. Camera dark current was shown to be negligible at normal accelerator operation. No discernible image distortion was found. Mechanical stability on gantry rotation was also assessed and image displacement of up to 5 mm at the isocentre was observed. It was concluded that the device could be used for dosimetry provided necessary precautions were observed and corrections made.

Direct measurement and analytical modeling of scatter in portal imaging

In this study a direct measurement of scatter in portal imaging for various air gaps and scatterer thicknesses at a beam energy of 6 MV is presented. The experimental data are compared with results from a Monte Carlo (MC) scatter model. In the regime where the air gap is larger than 9.3 cm the MC and the experiment agree. Based on this MC model an analytical model is developed, which takes all important interaction processes into account. It comprises a rigorous treatment of first order scattering and an estimation of photons scattered more than once within the phantom. This estimation is based on the assumption that higher order scattering can be considered as isotropically distributed around a certain scatter origin located in the midplane of the phantom. It is found that relative deviations between the MC model and the analytical model are of 2% to 3% in regions where scattering is very large.

A method to estimate the transit dose on the beam axis for verification of dose delivery with portal images

PURPOSE AND BACKGROUND

A feasibility study is performed to evaluate the possibility of using the transit dose of portal images on the beam axis to measure the accuracy in dose delivery. The algorithm and the method are tested on a breast phantom and on patients with a breast disease.

MATERIALS AND METHODS

To estimate the transit dose at various air gaps behind the patient, a method is proposed which applies, for a given air gap, the inverse square law to the primary component of the exit dose and an experimentally determined function for the scatter component of the exit dose. It is assumed that the primary component and the scattered component of the exit dose are given by the treatment planning system. The experimental function for the variation of the scattered component with the air gap, determined by phantom measurements, is modelled by an analytical function which contains only field size, air gap and one energy-dependent parameter.

RESULTS

The measurements on the breast phantom yield a maximum deviation between measured and estimated transit doses of 4.5%. The mean deviation is 0.9% with a standard deviation of the distribution of 2.3%. In vivo diode measurements on the same phantom yield a maximum deviation of 2.7%. Transit dose measurements on the beam axis for 45 portal images of breast patients show a mean deviation of 0.0% between the measured transit dose and the estimated transit dose. The standard deviation of the distribution is 4.4%. The method seems to be very sensitive to patient positioning and to discrepancies in breast thicknesses used for treatment planning.

CONCLUSION

Preliminary results on breast patients show that the method proposed to evaluate transit doses on the beam axis from portal images may be a valuable alternative to conventional in vivo exit dosimetry. The method can be implemented in a simple way and does not require additional time during the irradiation session, as exit dosimetry with diodes does. The transit dose is only considered in one point. Nevertheless, in the framework of quality assurance of treatment delivery, this study is an example of the possibilities of monitoring at the same time the visual evaluation of the irradiated volume as well as the dosimetric control (i.e. in Gy) of treatment delivery with portal images.

Variation of relative transit dose profiles with patient-detector distance

BACKGROUND AND PURPOSE

In view of using portal images for exit dosimetry, an experimental study is performed of relative transit dose profiles at different distances behind patients (and phantoms) and of their relation to the exit dose profile.

MATERIALS AND METHODS

Irregular, homogeneous polystyrene phantoms with a variable thickness to simulate head and neck (H&N) treatments (6-MV photon beam) are investigated by ionization chamber measurements performed close to the exit surface and at various distances behind the phantom (10, 20 and 30 cm). Similar measurements are performed for a rectangular phantom with large inhomogeneities (A1 and air). For one irregular homogeneous phantom and an irregular phantom containing an A1 inhomogeneity, ionization chamber measurements are performed at the exit surface, and a portal film image is taken at 30 cm behind the phantom. Portal films of a patient treated for a head and neck malignancy are evaluated for different air gaps behind the patient.

RESULTS

For the irregular phantoms, deviations up to 15% and more are observed between the exit dose profile (along the shaped surface of the phantom) and the transit profile close to the phantom (perpendicular to the beam axis). There is, however, a good agreement--within 3%--between the exit profile and the transit profile at 30 cm. For the rectangular, inhomogeneous phantom, the deviation between the exit profile and the transit dose profile at 30 cm does not exceed 5%; transit dose profiles overestimate the exit dose for the air cavity and underestimate the dose for the A1 inhomogeneity. Measurements on portal films of a H&N patient for different air gaps confirm the order of magnitude of the difference observed between transit dose profiles close to the patient and transit dose profiles at some distance behind the patient.

CONCLUSIONS

For 6-MV photon beam treatments with significant thickness variations (H&N), large variations (> 10%) are observed in transit dose profiles as a function of the air gap between the patient and the portal film. For this energy, a good agreement is found between the exit profile and the transit profile at about 30 cm behind the patient.

Verification of compensator thicknesses using a fluoroscopic electronic portal imaging device

A method is presented for verification of compensator thicknesses using a fluoroscopic electronic portal imaging device (EPID). The method is based on the measured transmission through the compensator, defined by the ratio of the portal dose with the compensator in the beam and the portal dose without the compensator in the beam. The transmission is determined with the EPID by dividing two images, acquired with and without compensator inserted, which are only corrected for the nonlinear response of the fluoroscopic system. The transmission has a primary and a scatter component. The primary component is derived from the measured transmission by subtracting the predicted scatter component. The primary component for each point is only related to the radiological thickness of the compensator along the ray line between the focus and that point. Compensator thicknesses are derived from the primary components taking into account off-axis variations in beam quality. The developed method has been tested for various compensators made of a granulate of stainless steel. The compensator thicknesses could be determined with an accuracy of 0.5 mm (1 s.d.), corresponding to a change in the transmitted dose of about 1% for a 10 MV beam. The method is fast, accurate, and insensitive to long-term output and beam profile fluctuations of the linear accelerator.

Evaluation of an electronic portal imaging device for transit dosimetry

The possibility of using a commercially available electronic portal imaging device for transit dosimetry was investigated. The detection unit of the device comprises a metal plate/fluorescence screen and a camera. Basic parameters of this system were investigated: stability, detector uniformity, dose-response curve, field-size dependence and phantom-thickness dependence. It was found that in terms of relative dosimetry, portal images corrected for detector non-uniformity are in good agreement (within 5%) with transit dose distributions measured by film dosimetry. For dose determination, it was found that the use of the device is hampered by an important field-size dependence and phantom-thickness dependence. Both correction factors should be applied if the device is to be used for this purpose.

Radiological thickness measurement using a liquid ionization chamber electronic portal imaging device

We present a method of calibrating the Portal Vision electronic portal imaging device to obtain radiological thickness maps for compensator design. In this method, coefficients are derived to describe the relationship between intensity and thickness for a set of water-equivalent blocks. The effects of four parameters were studied: (a) The dose response of the system was measured and found to be describable by a square-root function. (b) The calibration data and images were taken with a wedge in situ. The effects of using different wedges and different wedge orientations were investigated. The intrinsic accuracy of the accelerator/imager system was found to be 1.9 mm, for both 15 degrees and 30 degrees wedges. Changing the wedge orientation between calibration and imaging and rotating the calibration coefficients accordingly led to an error of 3.5 mm. (c) The variation in detector response with gantry angle was measured and corrected. The residual error in this process was 2.4 mm. (d) The use of a model to correct the effects of imaging with different field sizes was investigated and found to yield a residual error of 2.9 mm. The overall error in image calibrations was approximately 4 mm or 2% in dose. This is considered to be sufficiently small for the intended use of designing compensators for tangential breast irradiation.

A large-area ionization chamber for portal image calibration

Methods of removing the effects of linear accelerator (linac) output fluctuation from electronic portal images are described and compared. The output of the linac is measured using a specially constructed large-area ionization chamber during imaging and recorded with the image. The use of a dose-rate signal directly from the linac monitor chamber is discussed. Various versions of a quadratic thickness calibration scheme are tested, incorporating linac output data measured by the ionization chamber. Experimental results are presented showing that the incorporation of data from the ionization chamber gives improved absolute calibration accuracy and flatness. Immediately after calibration, the mean systematic thickness error in calibration of a uniform 136.8 mm water-equivalent slab was shown to be no more than 0.6 mm with a thickness variation within each image also of no more than +/-0.8 mm. This was true even when imaging with an unstable linac beam giving mean thickness errors between images of 8.8 mm and variations within each image of +/-4.9 mm without the ionization chamber correction. Up to one month after calibration, use of the ionization chamber to remove short-term linac fluctuations is shown to still keep mean thickness errors to less than 1.6 mm with variations within each image of no more than +/-1.4 mm.

Assessment of dose inhomogeneities in clinical practice by film dosimetry

AIM

To use portal images acquired in routine circumstances for assessment of midplane dose variations in the patient.

MATERIAL AND METHODS

Optical density readings are performed on routinely acquired Verification films of breast and ear-nose-throat (ENT) cancer patients and these readings are converted into relative doses with the sensitometric curve. ( 1 ) The impact of redistribution is evaluated on films taken close to the patient exit surface and at routine focus film distance. (2) Midplane doses are estimated from film readings to assess dose variations in the patient. The influence of wedges is evaluated. Film measurements doses are compared with calculated exit doses.

RESULTS

(1) In regions with large variations in the distance between the patient exit surface and the film but without inhomogeneities in tissue density, the relative doses distributions read on films acquired at large focus-film-distance (FFD) are proportional to exit doses. In regions with flat exit surfaces but with inhomogeneities in tissue density, the redistribution has only a small impact. (2) Large variations in relative midplane doses were found in both breast (85-115%) and ENT (-3.6 to +15%) patients. The application of a wedge was shown to increase dose homogeneity in the midplane. A good agreement (differences < 3%) was found between exit doses obtained from film readings and exit doses calculated with the treatment planning system (TPS).

CONCLUSION

Films acquired in routine circumstances at large FFD can be used to obtain information on exit doses and to assess midplane doses in breast and ENT, without the use of a TPS. Film dosimetry can also provide a quality assurance tool to check actually delivered doses in patients by comparing exit doses estimated on film to expected exit doses calculated by the TPS.

Linear accelerator output variations and their consequences for megavoltage imaging

An experimental study of radiation output intensity fluctuations of a Philips SL25 linear accelerator is presented. Measurements are obtained using an electronic portal imaging device, and the consequences of the measured fluctuations for various different applications of megavoltage imaging including portal imaging, transit dosimetry and megavoltage computed tomography (MVCT) are discussed with examples. Fluctuations in output of +/- 0.7% (1 SD) are seen on every radiation pulse after photon noise and uncertainties caused by the detection system have been accounted for. Large fluctuations are also seen during the initial beam stabilization period (15%), during normal accelerator operation after the beam has been on for more than 1 min (4.5%) and during are therapy as a repeatable function of gantry angle (9%). Such output intensity fluctuations are shown to produce image artifacts in portal imaging devices with scanned detector readout and can also produce systematic errors in detector calibration that would lead to uncertainty in transit dose calculations. The propagation of these intensity fluctuations through MVCT image reconstruction is shown to produce ring artifacts in the reconstructed image. Sample portal and MVCT images are presented. All observed fluctuations in accelerator output are well within the manufacturer's specifications and do not affect the total dose delivered during normal treatment. Finally, megavoltage imaging is shown to be a powerful tool for accelerator quality assurance and treatment verification.

Measurements of exit dose profiles in 60Co beams with a conventional portal film system

An important step in the verification of the reliability of portal films as in vivo dosemeters is the evaluation of the agreement between exit dose profiles and optical density profiles measured on the portal film. To test the possibilities of a conventional portal film system in 60Co beams suitable for head and neck irradiation, we verified the agreement between relative exit doses (measured by ionization chamber) and relative optical densities, on cubic homogeneous phantoms, on an homogeneous "step" phantom and on a cubic phantom including air and aluminium inhomogeneities. The optical density profiles were corrected with the appropriate sensitometric curves. For an homogeneous phantom 10.8 cm thick and with the film in contact with the phantom, the agreement was found to be excellent with a mean deviation of 0.8% and a maximum deviation of 1.5%. The agreement was worse when the air gap between the exit surface of the phantom and the portal film was increased (with an air gap equal to 15 cm the maximum deviation was 4%), and when the thickness of the phantom was increased (for a thickness of 14.4 cm the maximum deviation was 3.1%). The agreement was found to be acceptable for the "step" phantom too, with a mean deviation around 1% and a maximum deviation within 2% (air gap equal to zero). When air and aluminium inhomogeneities were incorporated into the phantom a maximum deviation of 6% and a mean deviation less than 3% were found. Furthermore, the relative optical density profiles show an underestimate of measured off-axis exit dose values under a high density inhomogeneity and a small overestimate under a low density inhomogeneity. Results suggest the possibility of using conventional portal films for exit relative dosimetry in head and neck irradiation with 60Co beams if the air gap is kept as small as possible.

An electronic portal imaging device for transit dosimetry

An electronic portal imaging device has been designed and constructed. It consists of an array of 128 CsI scintillation crystals coupled to photodiodes which is scanned across the field in 4 seconds. The linac is operated at a dose rate of 400 cGy min-1 and the dose delivered for image acquisition is approximately 27 cGy. The data acquisition controller is a stand-alone STE computer located within the scan arm. Sample images are presented showing contrast and spatial resolution of the system together with a humanoid phantom image and a clinical image of a breast cancer patient. The phantom images show the detector has a contrast resolution of 0.3% (at 15 mm diameter) and a spatial resolution of 2.5-3.2 mm. Images of uniform Perspex blocks have also been calibrated for thickness, indicating the system can measure radiological thickness to an accuracy of 2-3 mm of water. These results indicate the detector may be used for transit dosimetry applications including compensator design.