Investigating CT to CBCT image registration for head and neck proton therapy as a tool for daily dose recalculation

PURPOSE

Intensity modulated proton therapy (IMPT) of head and neck (H&N) cancer patients may be improved by plan adaptation. The decision to adapt the treatment plan based on a dose recalculation on the current anatomy requires a diagnostic quality computed tomography (CT) scan of the patient. As gantry-mounted cone beam CT (CBCT) scanners are currently being offered by vendors, they may offer daily or weekly updates of patient anatomy. CBCT image quality may not be sufficient for accurate proton dose calculation and it is likely necessary to perform CBCT CT number correction. In this work, the authors investigated deformable image registration (DIR) of the planning CT (pCT) to the CBCT to generate a virtual CT (vCT) to be used for proton dose recalculation.

METHODS

Datasets of six H&N cancer patients undergoing photon intensity modulated radiation therapy were used in this study to validate the vCT approach. Each dataset contained a CBCT acquired within 3 days of a replanning CT (rpCT), in addition to a pCT. The pCT and rpCT were delineated by a physician. A Morphons algorithm was employed in this work to perform DIR of the pCT to CBCT following a rigid registration of the two images. The contours from the pCT were deformed using the vector field resulting from DIR to yield a contoured vCT. The DIR accuracy was evaluated with a scale invariant feature transform (SIFT) algorithm comparing automatically identified matching features between vCT and CBCT. The rpCT was used as reference for evaluation of the vCT. The vCT and rpCT CT numbers were converted to stopping power ratio and the water equivalent thickness (WET) was calculated. IMPT dose distributions from treatment plans optimized on the pCT were recalculated with a Monte Carlo algorithm on the rpCT and vCT for comparison in terms of gamma index, dose volume histogram (DVH) statistics as well as proton range. The DIR generated contours on the vCT were compared to physician-drawn contours on the rpCT.

RESULTS

The DIR accuracy was better than 1.4 mm according to the SIFT evaluation. The mean WET differences between vCT (pCT) and rpCT were below 1 mm (2.6 mm). The amount of voxels passing 3%/3 mm gamma criteria were above 95% for the vCT vs rpCT. When using the rpCT contour set to derive DVH statistics from dose distributions calculated on the rpCT and vCT the differences, expressed in terms of 30 fractions of 2 Gy, were within [-4, 2 Gy] for parotid glands (D(mean)), spinal cord (D(2%)), brainstem (D(2%)), and CTV (D(95%)). When using DIR generated contours for the vCT, those differences ranged within [-8, 11 Gy].

CONCLUSIONS

In this work, the authors generated CBCT based stopping power distributions using DIR of the pCT to a CBCT scan. DIR accuracy was below 1.4 mm as evaluated by the SIFT algorithm. Dose distributions calculated on the vCT agreed well to those calculated on the rpCT when using gamma index evaluation as well as DVH statistics based on the same contours. The use of DIR generated contours introduced variability in DVH statistics.

Using CBCT for pretreatment range check in proton therapy: a phantom study for prostate treatment by anterior-posterior beam

This study explores the potential of cone‐beam computed tomography (CBCT) for monitoring relative beam range variations due to daily changes in patient anatomy for prostate treatment by anterior proton beams. CBCT was used to image an anthropomorphic pelvic phantom, in eight sessions on eight different days. In each session, the phantom was scanned twice, first at a standard position as determined by the room lasers, and then after it was shifted by 10 mm translation randomly along one of the X, Y, or Z directions. The filling of the phantom bladder with water was not refreshed from day to day, inducing gradual change of the water‐equivalent path length (WEPL) across the bladder. MIMvista (MIM) software was used to perform image registration and re‐alignment of all the scans with the scan from the first session. The XiO treatment planning system was used to perform data analysis. It was found that, although the Hounsfield unit numbers in CBCT have substantially larger fluctuations than those in diagnostic CT, CBCT datasets taken for daily patient positioning could potentially be used to monitor changes in patient anatomy. The reproducibility of the WEPL, computed using CBCT along anterior–posterior (AP) paths across and around the phantom prostate, over a volume of 360 cc, is sufficient for detecting daily WEPL variations that are equal to or larger than 3 mm. This result also applies to CBCT scans of the phantom after it is randomly shifted from the treatment position by 10 mm. limiting the interest to WEPL variation over a specific path within the same CBCT slice, one can detect WEPL variation smaller than 1 mm. That is the case when using CBCT for tracking daily change of the WEPL across the phantom bladder that was induced by spontaneous change in the bladder filling due to evaporation. In summary, the phantom study suggests that CBCT can be used for monitoring day to day WEPL variations in a patient. The method can detect WEPL variation equal to or greater than 3 mm. The study calls for further investigation using the CBCT data from real patients. If confirmed with real patients' data, CBCT could become, in addition to patient setup, a standard tool for proton therapy pretreatment beam range check.

PACS number: 87.55.Tm

Phantom based evaluation of CT to CBCT image registration for proton therapy dose recalculation

The ability to perform dose recalculation on the anatomy of the day is important in the context of adaptive proton therapy. The objective of this study was to investigate the use of deformable image registration (DIR) and cone beam CT (CBCT) imaging to generate the daily stopping power distribution of the patient. We investigated the deformation of the planning CT scan (pCT) onto daily CBCT images to generate a virtual CT (vCT) using a deformable phantom designed for the head and neck (H & N) region. The phantom was imaged at a planning CT scanner in planning configuration, yielding a pCT and in deformed, treatment day configuration, yielding a reference CT (refCT). The treatment day configuration was additionally scanned at a CBCT scanner. A Morphons DIR algorithm was used to generate a vCT. The accuracy of the vCT was evaluated by comparison to the refCT in terms of corresponding features as identified by an adaptive scale invariant feature transform (aSIFT) algorithm. Additionally, the vCT CT numbers were compared to those of the refCT using both profiles and regions of interest and the volumes and overlap (DICE coefficients) of various phantom structures were compared. The water equivalent thickness (WET) of the vCT, refCT and pCT were also compared to evaluate proton range differences. Proton dose distributions from the same initial fluence were calculated on the refCT, vCT and pCT and compared in terms of proton range. The method was tested on a clinical dataset using a replanning CT scan acquired close in time to a CBCT scan as reference using the WET evaluation. Results from the aSIFT investigation suggest a deformation accuracy of 2-3 mm. The use of the Morphon algorithm did not distort CT number intensity in uniform regions and WET differences between vCT and refCT were of the order of 2% of the proton range. This result was confirmed by proton dose calculations. The patient results were consistent with phantom observations. In conclusion, our phantom study suggests the vCT approach is adequate for proton dose recalculation on the basis of CBCT imaging.