Thomas E Marchant

Shading correction algorithm for improvement of cone-beam CT images in radiotherapy

Cone-beam CT (CBCT) images have recently become an established modality for treatment verification in radiotherapy. However, identification of soft-tissue structures and the calculation of dose distributions based on CBCT images is often obstructed by image artefacts and poor consistency of density calibration. A robust method for voxel-by-voxel enhancement of CBCT images using a priori knowledge from the planning CT scan has been developed and implemented. CBCT scans were enhanced using a low spatial frequency grey scale shading function generated with the aid of a planning CT scan from the same patient. This circumvents the need for exact correspondence between CBCT and CT and the process is robust to the appearance of unshared features such as gas pockets. Enhancement was validated using patient CBCT images. CT numbers in regions of fat and muscle tissue in the processed CBCT were both within 1% of the values in the planning CT, as opposed to 10-20% different for the original CBCT. Visual assessment of processed CBCT images showed improvement in soft-tissue visibility, although some cases of artefact introduction were observed.

Measurement of inter and intra fraction organ motion in radiotherapy using cone beam CT projection images

A method is presented for extraction of intra and inter fraction motion of seeds/markers within the patient from cone beam CT (CBCT) projection images. The position of the marker is determined on each projection image and fitted to a function describing the projection of a fixed point onto the imaging panel at different gantry angles. The fitted parameters provide the mean marker position with respect to the isocentre. Differences between the theoretical function and the actual projected marker positions are used to estimate the range of intra fraction motion and the principal motion axis in the transverse plane. The method was validated using CBCT projection images of a static marker at known locations and of a marker moving with known amplitude. The mean difference between actual and measured motion range was less than 1 mm in all directions, although errors of up to 5 mm were observed when large amplitude motion was present in an orthogonal direction. In these cases it was possible to calculate the range of motion magnitudes consistent with the observed marker trajectory. The method was shown to be feasible using clinical CBCT projections of a pancreas cancer patient.

Inter-fraction motion and dosimetric consequences during breast intensity-modulated radiotherapy (IMRT)

BACKGROUND AND PURPOSE Intensity-modulated radiotherapy (IMRT) can improve dose homogeneity within the breast planned target volume (PTV), but may be more susceptible to patient/organ motion than standard tangential radiotherapy (RT). We used daily cone-beam CT (CBCT) imaging to assess inter-fraction motion during breast IMRT and its subsequent impact on IMRT and standard RT dose homogeneity. MATERIALS AND METHODS Ten breast cancer patients selected for IMRT were studied. CBCT images were acquired immediately after daily treatment. Automatic image co-registration was used to determine patient positioning variations. Daily PTV contours were used to calculate PTV variations and daily delivered IMRT and theoretically planned tangential RT dose. RESULTS Group systematic (and random) setup errors detected by CBCT were 5.7 (3.9)mm laterally, 2.8 (3.5)mm vertically and 2.3 (3.2)mm longitudinally. Rotations >2 degrees in any axis occurred on 53/106 (50%) occasions. Daily PTV volume varied up to 23%. IMRT dose homogeneity was superior at planning and throughout the treatment compared with standard RT (1.8% vs. 15.8% PTV received >105% planned mean dose), despite increased motion sensitivity. CONCLUSIONS CBCT revealed inadequacies of current patient positioning and verification procedures during breast RT and confirmed improved dose homogeneity using IMRT for the patients studied.

Monitoring Dosimetric Impact of Weight Loss With Kilovoltage (KV) Cone Beam CT (CBCT) During Parotid-Sparing IMRT and Concurrent Chemotherapy

PURPOSE Parotid-sparing head-and-neck intensity-modulated radiotherapy (IMRT) can reduce long-term xerostomia. However, patients frequently experience weight loss and tumor shrinkage during treatment. We evaluate the use of kilovoltage (kV) cone beam computed tomography (CBCT) for dose monitoring and examine if the dosimetric impact of such changes on the parotid and critical neural structures warrants replanning during treatment. METHODS AND MATERIALS Ten patients with locally advanced oropharyngeal cancer were treated with contralateral parotid-sparing IMRT concurrently with platinum-based chemotherapy. Mean doses of 65 Gy and 54 Gy were delivered to clinical target volume (CTV)1 and CTV2, respectively, in 30 daily fractions. CBCT was prospectively acquired weekly. Each CBCT was coregistered with the planned isocenter. The spinal cord, brainstem, parotids, larynx, and oral cavity were outlined on each CBCT. Dose distributions were recalculated on the CBCT after correcting the gray scale to provide accurate Hounsfield calibration, using the original IMRT plan configuration. RESULTS Planned contralateral parotid mean doses were not significantly different to those delivered during treatment (p > 0.1). Ipsilateral and contralateral parotids showed a mean reduction in volume of 29.7% and 28.4%, respectively. There was no significant difference between planned and delivered maximum dose to the brainstem (p = 0.6) or spinal cord (p = 0.2), mean dose to larynx (p = 0.5) and oral cavity (p = 0.8). End-of-treatment mean weight loss was 7.5 kg (8.8% of baseline weight). Despite a ≥10% weight loss in 5 patients, there was no significant dosimetric change affecting the contralateral parotid and neural structures. CONCLUSIONS Although patient weight loss and parotid volume shrinkage was observed, overall, there was no significant excess dose to the organs at risk. No replanning was felt necessary for this patient cohort, but a larger patient sample will be investigated to further confirm these results. Nevertheless, kilovoltage CBCT is a valuable tool for patient setup verification and monitoring of dosimetric variation during radiotherapy.

Quantifying motion for pancreatic radiotherapy margin calculation.

BACKGROUND AND PURPOSE Pancreatic radiotherapy (RT) is limited by uncertain target motion. We quantified 3D patient/organ motion during pancreatic RT and calculated required treatment margins. MATERIALS AND METHODS Cone-beam computed tomography (CBCT) and orthogonal fluoroscopy images were acquired post-RT delivery from 13 patients with locally advanced pancreatic cancer. Bony setup errors were calculated from CBCT. Inter- and intra-fraction fiducial (clip/seed/stent) motion was determined from CBCT projections and orthogonal fluoroscopy. RESULTS Using an off-line CBCT correction protocol, systematic (random) setup errors were 2.4 (3.2), 2.0 (1.7) and 3.2 (3.6)mm laterally (left-right), vertically (anterior-posterior) and longitudinally (cranio-caudal), respectively. Fiducial motion varied substantially. Random inter-fractional changes in mean fiducial position were 2.0, 1.6 and 2.6mm; 95% of intra-fractional peak-to-peak fiducial motion was up to 6.7, 10.1 and 20.6mm, respectively. Calculated clinical to planning target volume (CTV-PTV) margins were 1.4 cm laterally, 1.4 cm vertically and 3.0 cm longitudinally for 3D conformal RT, reduced to 0.9, 1.0 and 1.8 cm, respectively, if using 4D planning and online setup correction. CONCLUSIONS Commonly used CTV-PTV margins may inadequately account for target motion during pancreatic RT. Our results indicate better immobilisation, individualised allowance for respiratory motion, online setup error correction and 4D planning would improve targeting.

Reduction of motion artefacts in on-board cone beam CT by warping of projection images.

OBJECTIVE We describe the development and testing of a motion correction method for flat panel imager-based cone beam CT (CBCT) based on warping of projection images. METHODS Markers within or on the surface of the patient were tracked and their mean three-dimensional (3D) position calculated. The two-dimensional (2D) cone beam projection images were then warped before reconstruction to place each marker at the projection from its mean 3D position. The motion correction method was tested using simulated cone beam projection images of a deforming virtual phantom, real CBCT images of a moving breast phantom and clinical CBCT images of a patient with breast cancer and another with pancreatic cancer undergoing radiotherapy. RESULTS In phantom studies, the method was shown to greatly reduce motion artefacts in the locality of the radiotherapy target and allowed the true surface shape to be accurately recovered. The breast phantom motion-compensated surface was within 1 mm of the true surface shape for 90% of surface points and greater than 2 mm from the true surface at only 2% of points. Clinical CBCT images showed improved image quality in the locality of the radiotherapy target after motion correction. CONCLUSION The proposed method is effective in reducing motion artefacts in CBCT images.

Automatic tracking of implanted fiducial markers in cone beam CT projection images

PURPOSE This paper describes a novel method for simultaneous intrafraction tracking of multiple fiducial markers. Although the proposed method is generic and can be adopted for a number of applications including fluoroscopy based patient position monitoring and gated radiotherapy, the tracking results presented in this paper are specific to tracking fiducial markers in a sequence of cone beam CT projection images. METHODS The proposed method is accurate and robust thanks to utilizing the mean shift and random sampling principles, respectively. The performance of the proposed method was evaluated with qualitative and quantitative methods, using data from two pancreatic and one prostate cancer patients and a moving phantom. The ground truth, for quantitative evaluation, was calculated based on manual tracking preformed by three observers. RESULTS The average dispersion of marker position error calculated from the tracking results for pancreas data (six markers tracked over 640 frames, 3840 marker identifications) was 0.25 mm (at iscoenter), compared with an average dispersion for the manual ground truth estimated at 0.22 mm. For prostate data (three markers tracked over 366 frames, 1098 marker identifications), the average error was 0.34 mm. The estimated tracking error in the pancreas data was < 1 mm (2 pixels) in 97.6% of cases where nearby image clutter was detected and in 100.0% of cases with no nearby image clutter. CONCLUSIONS The proposed method has accuracy comparable to that of manual tracking and, in combination with the proposed batch postprocessing, superior robustness. Marker tracking in cone beam CT (CBCT) projections is useful for a variety of purposes, such as providing data for assessment of intrafraction motion, target tracking during rotational treatment delivery, motion correction of CBCT, and phase sorting for 4D CBCT.

An analysis of breast motion using high-frequency, dense surface points captured by an optical sensor during radiotherapy treatment delivery

Patient motion is an important factor affecting the quality of external beam radiotherapy in breast patients. We analyse the motion of a dense set of surface points on breast patients throughout their treatment schedule to assess the magnitude and stability of motion, in particular, with respect to breast volume. We use an optical sensor to measure the surface motion of 13 breast cancer patients. Patients were divided into two cohorts dependent upon breast volume. Measurements were made during radiotherapy treatment beam delivery for an average of 12 fractions per patient (total 158 datasets). The motion of each surface point is parameterized in terms of its period, amplitude and relative phase. Inter-comparison of the motion parameters across treatment schedules and between patients is made through the creation of corresponding regions on the breast surfaces. The motion period is spatially uniform and is similar in both patient groups (mean 4 s), with the small volume cohort exhibiting greater inter-fraction period variability. The mean motion amplitude is also similar in both groups with a range between 2 mm and 4 mm and an inter-fraction variability generally less than 1 mm. There is a phase lag of up to 0.4 s across the breast, led by the sternum. Breast patient motion is reasonably stable between and during treatment fractions, with the large volume cohort exhibiting greater repeatability than the small volume one.

A megavoltage scatter correction technique for cone-beam CT images acquired during VMAT delivery.

Kilovoltage cone-beam CT (kV CBCT) can be acquired during the delivery of volumetric modulated arc therapy (VMAT), in order to obtain an image of the patient during treatment. However, the quality of such CBCTs is degraded by megavoltage (MV) scatter from the treatment beam onto the imaging panel. The objective of this paper is to introduce a novel MV scatter correction method for simultaneous CBCT during VMAT, and to investigate its effectiveness when compared to other techniques. The correction requires the acquisition of a separate set of images taken during VMAT delivery, while the kV beam is off. These images--which contain only the MV scatter contribution on the imaging panel--are then used to correct the corresponding kV/MV projections. To test this method, CBCTs were taken of an image quality phantom during VMAT delivery and measurements of contrast to noise ratio were made. Additionally, the correction was applied to the datasets of three VMAT prostate patients, who also received simultaneous CBCTs. The clinical image quality was assessed using a validated scoring system, comparing standard CBCTs to the uncorrected simultaneous CBCTs and a variety of correction methods. Results show that the correction is able to recover some of the low and high-contrast signal to noise ratio lost due to MV scatter. From the patient study, the corrected CBCT scored significantly higher than the uncorrected images in terms of the ability to identify the boundary between the prostate and surrounding soft tissue. In summary, a simple MV scatter correction method has been developed and, using both phantom and patient data, is shown to improve the image quality of simultaneous CBCTs taken during VMAT delivery.

Comprehensive Monte Carlo study of patient doses from cone-beam CT imaging in radiotherapy.

Accurate knowledge of ionizing radiation dose from cone-beam CT (CBCT) imaging in radiotherapy is important to allow concomitant risks to be estimated and for justification of imaging exposures. This study uses a Monte Carlo CBCT model to calculate imaging dose for a wide range of imaging protocols for male and female patients. The Elekta XVI CBCT system was modeled using GATE and simulated doses were validated against measurements in a water tank and thorax phantom. Imaging dose was simulated in the male and female ICRP voxel phantoms for a variety of anatomical sites and imager settings (different collimators, filters, full and partial rotation). The resulting dose distributions were used to calculate effective doses for each scan protocol. The Monte Carlo simulated doses agree with validation measurements within 5% and 10% for water tank and thorax phantom respectively. Effective dose for head CBCT scans was generally lower for scans centred on the pituitary than the larynx (0.03 mSv versus 0.06 mSv for male ICRP phantom). Pelvis CBCT scan effective dose was higher for the female than male phantom (5.11 mSv versus 2.80 mSv for M15 collimator scan), principally due to the higher dose received by gonads for the female scan. Medium field of view thorax scan effective doses ranged from 1.38-3.19 mSv depending on scan length and phantom sex. Effective dose for half rotation thorax scans with offset isocentre varied by almost a factor of three depending on laterality of the isocentre, patient sex and imaged field length. The CBCT imaging doses simulated here reveal large variations in dose depending on imaging isocentre location, patient sex and partial rotation angles. This information may be used to estimate risks from CBCT and to optimize CBCT imaging protocols.