David A. Jaffray

Patient dose from kilovoltage cone beam computed tomography imaging in radiation therapy

Kilovoltage cone-beam computerized tomography (kV-CBCT) systems integrated into the gantry of linear accelerators can be used to acquire high-resolution volumetric images of the patient in the treatment position. Using on-line software and hardware, patient position can be determined accurately with a high degree of precision and, subsequently, set-up parameters can be adjusted to deliver the intended treatment. While the patient dose due to a single volumetric imaging acquisition is small compared to the therapy dose, repeated and daily image guidance procedures can lead to substantial dose to normal tissue. The dosimetric properties of a clinical CBCT system have been studied on an Elekta linear accelerator (Synergy RP, XVI system) and additional measurements performed on a laboratory system with identical geometry. Dose measurements were performed with an ion chamber and MOSFET detectors at the center, periphery, and surface of 30 and 16-cm-diam cylindrical shaped water phantoms, as a function of x-ray energy and longitudinal field-of-view (FOV) settings of 5,10,15, and 26 cm. The measurements were performed for full 360 degrees CBCT acquisition as well as for half-rotation scans for 120 kVp beams using the 30-cm-diam phantom. The dose at the center and surface of the body phantom were determined to be 1.6 and 2.3 cGy for a typical imaging protocol, using full rotation scan, with a technique setting of 120 kVp and 660 mAs. The results of our measurements have been presented in terms of a dose conversion factor fCBCT, expressed in cGy/R. These factors depend on beam quality and phantom size as well as on scan geometry and can be utilized to estimate dose for any arbitrary mAs setting and reference exposure rate of the x-ray tube at standard distance. The results demonstrate the opportunity to manipulate the scanning parameters to reduce the dose to the patient by employing lower energy (kVp) beams, smaller FOV, or by using half-rotation scan.

Volume CT with a flat-panel detector on a mobile, isocentric C-arm: pre-clinical investigation in guidance of minimally invasive surgery.

A mobile isocentric C-arm (Siemens PowerMobil) has been modified in our laboratory to include a large area flat-panel detector (in place of the x-ray image intensifier), providing multi-mode fluoroscopy and cone-beam computed tomography (CT) imaging capability. This platform represents a promising technology for minimally invasive, image-guided surgical procedures where precision in the placement of interventional tools with respect to bony and soft-tissue structures is critical. The image quality and performance in surgical guidance was investigated in pre-clinical evaluation in image-guided spinal surgery. The control, acquisition, and reconstruction system are described. The reproducibility of geometric calibration, essential to achieving high three-dimensional (3D) image quality, is tested over extended time scales (7 months) and across a broad range in C-arm angulation (up to 45 degrees), quantifying the effect of improper calibration on spatial resolution, soft-tissue visibility, and image artifacts. Phantom studies were performed to investigate the precision of 3D localization (viz., fiber optic probes within a vertebral body) and effect of lateral projection truncation (limited field of view) on soft-tissue detectability in image reconstructions. Pre-clinical investigation was undertaken in a specific spinal procedure (photodynamic therapy of spinal metastases) in five animal subjects (pigs). In each procedure, placement of fiber optic catheters in two vertebrae (L1 and L2) was guided by fluoroscopy and cone-beam CT. Experience across five procedures is reported, focusing on 3D image quality, the effects of respiratory motion, limited field of view, reconstruction filter, and imaging dose. Overall, the intraoperative cone-beam CT images were sufficient for guidance of needles and catheters with respect to bony anatomy and improved surgical performance and confidence through 3D visualization and verification of transpedicular trajectories and tool placement. Future investigation includes improvement in image quality, particularly regarding x-ray scatter, motion artifacts and field of view, and integration with optical tracking and navigation systems.

A simple, direct method for x-ray scatter estimation and correction in digital radiography and cone-beam CT

X-ray scatter poses a significant limitation to image quality in cone-beam CT (CBCT), resulting in contrast reduction, image artifacts, and lack of CT number accuracy. We report the performance of a simple scatter correction method in which scatter fluence is estimated directly in each projection from pixel values near the edge of the detector behind the collimator leaves. The algorithm operates on the simple assumption that signal in the collimator shadow is attributable to x-ray scatter, and the 2D scatter fluence is estimated by interpolating between pixel values measured along the top and bottom edges of the detector behind the collimator leaves. The resulting scatter fluence estimate is subtracted from each projection to yield an estimate of the primary-only images for CBCT reconstruction. Performance was investigated in phantom experiments on an experimental CBCT bench-top, and the effect on image quality was demonstrated in patient images (head, abdomen, and pelvis sites) obtained on a preclinical system for CBCT-guided radiation therapy. The algorithm provides significant reduction in scatter artifacts without compromise in contrast-to-noise ratio (CNR). For example, in a head phantom, cupping artifact was essentially eliminated, CT number accuracy was restored to within 3%, and CNR (breast-to-water) was improved by up to 50%. Similarly in a body phantom, cupping artifact was reduced by at least a factor of 2 without loss in CNR. Patient images demonstrate significantly increased uniformity, accuracy, and contrast, with an overall improvement in image quality in all sites investigated. Qualitative evaluation illustrates that soft-tissue structures that are otherwise undetectable are clearly delineated in scatter-corrected reconstructions. Since scatter is estimated directly in each projection, the algorithm is robust with respect to system geometry, patient size and heterogeneity, patient motion, etc. Operating without prior information, analytical modeling, or Monte Carlo, the technique is easily incorporated as a preprocessing step in CBCT reconstruction to provide significant scatter reduction.

Accurate technique for complete geometric calibration of cone-beam computed tomography systems

Cone-beam computed tomography systems have been developed to provide in situ imaging for the purpose of guiding radiation therapy. Clinical systems have been constructed using this approach, a clinical linear accelerator (Elekta Synergy RP) and an iso-centric C-arm. Geometric calibration involves the estimation of a set of parameters that describes the geometry of such systems, and is essential for accurate image reconstruction. We have developed a general analytic algorithm and corresponding calibration phantom for estimating these geometric parameters in cone-beam computed tomography (CT) systems. The performance of the calibration algorithm is evaluated and its application is discussed. The algorithm makes use of a calibration phantom to estimate the geometric parameters of the system. The phantom consists of 24 steel ball bearings (BBs) in a known geometry. Twelve BBs are spaced evenly at 30 deg in two plane-parallel circles separated by a given distance along the tube axis. The detector (e.g., a flat panel detector) is assumed to have no spatial distortion. The method estimates geometric parameters including the position of the x-ray source, position, and rotation of the detector, and gantry angle, and can describe complex source-detector trajectories. The accuracy and sensitivity of the calibration algorithm was analyzed. The calibration algorithm estimates geometric parameters in a high level of accuracy such that the quality of CT reconstruction is not degraded by the error of estimation. Sensitivity analysis shows uncertainty of 0.01 degrees (around beam direction) to 0.3 degrees (normal to the beam direction) in rotation, and 0.2 mm (orthogonal to the beam direction) to 4.9 mm (beam direction) in position for the medical linear accelerator geometry. Experimental measurements using a laboratory bench Cone-beam CT system of known geometry demonstrate the sensitivity of the method in detecting small changes in the imaging geometry with an uncertainty of 0.1 mm in transverse and vertical (perpendicular to the beam direction) and 1.0 mm in the longitudinal (beam axis) directions. The calibration algorithm was compared to a previously reported method, which uses one ball bearing at the isocenter of the system, to investigate the impact of more precise calibration on the image quality of cone-beam CT reconstruction. A thin steel wire located inside the calibration phantom was imaged on the conebeam CT lab bench with and without perturbations in source and detector position during the scan. The described calibration method improved the quality of the image and the geometric accuracy of the object reconstructed, improving the full width at half maximum of the wire by 27.5% and increasing contrast of the wire by 52.8%. The proposed method is not limited to the geometric calibration of cone-beam CT systems but can be used for many other systems, which consist of one or more point sources and area detectors such as calibration of megavoltage (MV) treatment system (focal spot movement during the beam delivery, MV source trajectory versus gantry angle, the axis of collimator rotation, and couch motion), cross calibration between Kilovolt imaging and MV treatment system, and cross calibration between multiple imaging systems. Using the complete information of the system geometry, it was demonstrated that high image quality in CT reconstructions is possible even in systems with large geometric nonidealities.

The stability of mechanical calibration for a kV cone beam computed tomography system integrated with linear accelerator

The geometric accuracy and precision of an image-guided treatment system were assessed. Image guidance is performed using an x-ray volume imaging (XVI) system integrated with a linear accelerator and treatment planning system. Using an amorphous silicon detector and x-ray tube, volumetric computed tomography images are reconstructed from kilovoltage radiographs by filtered backprojection. Image fusion and assessment of geometric targeting are supported by the treatment planning system. To assess the limiting accuracy and precision of image-guided treatment delivery, a rigid spherical target embedded in an opaque phantom was subjected to 21 treatment sessions over a three-month period. For each session, a volumetric data set was acquired and loaded directly into an active treatment planning session. Image fusion was used to ascertain the couch correction required to position the target at the prescribed iso-center. Corrections were validated independently using megavoltage electronic portal imaging to record the target position with respect to symmetric treatment beam apertures. An initial calibration cycle followed by repeated image-guidance sessions demonstrated the XVI system could be used to relocate an unambiguous object to within less than 1 mm of the prescribed location. Treatment could then proceed within the mechanical accuracy and precision of the delivery system. The calibration procedure maintained excellent spatial resolution and delivery precision over the duration of this study, while the linear accelerator was in routine clinical use. Based on these results, the mechanical accuracy and precision of the system are ideal for supporting high-precision localization and treatment of soft-tissue targets.

Respiration correlated cone-beam computed tomography and 4DCT for evaluating target motion in Stereotactic Lung Radiation Therapy

An image-guidance process for using cone-beam computed tomography (CBCT) for stereotactic body radiation therapy (SBRT) of peripheral lung lesions is presented. Respiration correlated CBCT on the treatment unit and four dimensional computed tomography (4DCT) from planning are evaluated for assessing respiration-induced target motion during planning and treatment fractions. Image-guided SBRT was performed for 12 patients (13 lesions) with inoperable early stage non-small cell lung carcinoma. Kilovoltage (kV) projections were acquired over a 360 degree gantry rotation and sorted based on the pixel value of an image-based aperture located at the air-tissue interface of the diaphragm. The sorted projections were reconstructed to provide volumetric respiration correlated CBCT image datasets at different phases of the respiratory cycle. The 4D volumetric datasets were directly compared with 4DCT datasets acquired at the time of planning. For ten of 12 patients treated, the lung tumour motion, as measured by respiration correlated CBCT on the treatment unit, was consistent with the tumour motion measured by 4DCT at the time of planning. However, in two patients, maximum discrepancies observed were 6 and 10 mm in the anterior-posterior and superior-inferior directions, respectively. Respiration correlated CBCT acquired on the treatment unit allows target motion to be assessed for each treatment fraction, allows target localization based on different phases on the breathing cycle, and provides the facility for adaptive margin design in radiation therapy of lung malignancies. The current study has shown that the relative motion and position of the tumour at the time of treatment may not match that of the planning 4DCT scan. Therefore, application of breathing motion data acquired at simulation for tracking or gating radiation therapy may not be suitable for all patients — even those receiving short course treatment techniques such as SBRT.

Compensators for dose and scatter management in cone‐beam computed tomography

The ability of compensators (e.g., bow-tie filters) designed for kV cone-beam computed tomography (CT) to reduce both scatter reaching the detector and dose to the patient is investigated. Scattered x rays reaching the detector are widely recognized as one of the most significant challenges to cone-beam CT imaging performance. With cone-beam CT gaining popularity as a method of guiding treatments in radiation therapy, any methods that have the potential to reduce the dose to patients and/or improve image quality should be investigated. Simple compensators with a design that could realistically be implemented on a cone-beam CT imaging system have been constructed to determine the magnitude of reduction of scatter and/or dose for various cone-beam CT imaging conditions. Depending on the situation, the compensators were shown to reduce x-ray scatter at the detector and dose to the patient by more than a factor of 2. Further optimization of the compensators is a possibility to achieve greater reductions in both scatter and dose.

Characterization and real-time optical measurements of the ionizing radiation dose response for a new radiochromic medium

A new radiochromic film, GafChromic EBT, was investigated for use in a real-time radiation dosimetry system. It was found to be approximately eight times more sensitive to ionizing radiation dose, exhibited less postexposure development and achieved stable readout faster than one of its predecessors, GafChromic MD-55. A clear distinction in change in optical density between exposure and postexposure was observed, but the measurements obtained during exposure were not linear with time or dose. This could not be explained by a shift in wavelength of maximum change in absorbance, as it was stable at approximately 636 nm during the entire exposure range (up to 9.52 Gy). Increasing the spectral window of interest over which calculations were performed did little to correct the nonlinearity. The radiochromic film exhibited small dose rate dependence in real-time measurements, with an increase in standard deviation of change in optical density measurements from 0.9% to 1.8% over a sixfold variation in dose rate. Overall, GafChromic EBT has increased sensitivity and decreased postexposure darkening, and this bodes well for its potential role as a radiation dosimeter, including real-time applications.

Geometric calibration of a mobile C-arm for intraoperative cone-beam CT.

A geometric calibration method that determines a complete description of source-detector geometry was adapted to a mobile C-arm for cone-beam computed tomography (CBCT). The non-iterative calibration algorithm calculates a unique solution for the positions of the source (X(s), Y(s), Z(s)), detector (X(d), Y(d), Z(d)), piercing point (U(o), V(o)), and detector rotation angles (phi, theta, eta) based on projections of a phantom consisting of two plane-parallel circles of ball bearings encased in a cylindrical acrylic tube. The prototype C-arm system was based on a Siemens PowerMobil modified to provide flat-panel CBCT for image-guided interventions. The magnitude of geometric nonidealities in the source-detector orbit was measured, and the short-term (approximately 4 h) and long-term (approximately 6 months) reproducibility of the calibration was evaluated. The C-arm exhibits large geometric nonidealities due to mechanical flex, with maximum departures from the average semicircular orbit of deltaU(o) = 15.8 mm and deltaV(o) = 9.8 mm (for the piercing point), deltaX and deltaY = 6-8 mm and deltaZ = 1 mm (for the source and detector), and deltaphi approximately 2.9 degrees, deltatheta approximately 1.9 degrees, and delta eta approximately 0.8 degrees (for the detector tilt/rotation). Despite such significant departures from a semicircular orbit, these system parameters were found to be reproducible, and therefore correctable by geometric calibration. Short-term reproducibility was < 0.16 mm (subpixel) for the piercing point coordinates, < 0.25 mm for the source-detector X and Y, < 0.035 mm for the source-detector Z, and < 0.02 degrees for the detector angles. Long-term reproducibility was similarly high, demonstrated by image quality and spatial resolution measurements over a period of 6 months. For example, the full-width at half-maximum (FWHM) in axial images of a thin steel wire increased slightly as a function of the time (delta) between calibration and image acquisition: FWHM=0.62, 0.63, 0.66, 0.71, and 0.72 mm at delta = 0 s, 1 h, 1 day, 1 month, and 6 months, respectively. For ongoing clinical trials in CBCT-guided surgery at our institution, geometric calibration is conducted monthly to provide sufficient three-dimensional (3D) image quality while managing time and workflow considerations of the calibration and quality assurance process. The sensitivity of 3D image quality to each of the system parameters was investigated, as was the tolerance to systematic and random errors in the geometric parameters, showing the most sensitive parameters to be the piercing point coordinates (U(o), V(o)) and in-plane positions of the source (X(s), Y(s)) and detector (X(d), Y(d)). Errors in the out-of-plane position of the source (Z(s)) and detector (Z(d)) and the detector angles (phi, theta, eta) were shown to have subtler effects on 3D image quality.

Energy dependence "75 kVp to 18 MV… of radiochromic films assessed using a real-time optical dosimeter

The response of radiochromic film, GafChromic EBT, was investigated for dependence on x-ray beam energy using a previously reported real-time optical readout approach. X-ray beams of energy from 75 kVp to 18 MV were employed. The dose-induced change in optical density for the EBT film was compared to values obtained for GafChromic HS and MD-55 films, exposed under the same conditions. All responses were normalized to that obtained for 60Co irradiation. While change in optical density for 1 Gy of applied dose as measured with HS and MD-55 films decreased by approximately 40% at low energies, the mean change in optical density of EBT film remained within 3% of that in the 60Co beam over the entire energy range.