Browsing by Subject "Digital Tomosynthesis"

Imaging respiratory induced tumor motion in the radiation therapy treatment room could eliminate the necessity for large motion encompassing margins that result in excessive irradiation of healthy tissues. Currently available image guidance technologies are ill-suited for this task. Two-dimensional fluoroscopic images are acquired with sufficient speed to image respiratory motion. However, volume information is not present, and soft tissue structures are often not visible because a large volume is projected onto a single plane. Currently available volumetric imaging modalities are not acquired with sufficient speed to capture full motion trajectory information. Four-dimensional cone-beam computed tomography (4D CBCT) using a gantry mounted 2D flat panel imaging device has been proposed but has been limited by high doses, long scan times and severe under-sampling artifacts. The focus of the work completed in this thesis was to find ways to improve 4D imaging using a gantry mounted 2D kV imaging system. Specifically, the goals were to investigate methods for minimizing imaging dose and scan time while achieving consistent, controllable, high quality 4D images.

First, we introduced four-dimensional digital tomosynthesis (4D DTS) and characterized its potential for 3D motion analysis using a motion phantom. The motion phantom was programmed to exhibit motion profiles with various known amplitudes in all three dimensions and scanned using a 2D kV imaging system mounted on a linear accelerator. Two arcs of projection data centered about the anterior-posterior and lateral axes were used to reconstruct phase resolved DTS coronal and sagittal images. Respiratory signals were obtained by analyzing projection data, and these signals were used to derive phases for each of the projection images. Projection images were sorted according to phase, and DTS phase images were reconstructed for each phase bin. 4D DTS target location accuracies for peak inhalation and peak exhalation in all three dimensions were limited only by the 0.5 mm pixel resolution for all DTS scan angles. The average localization errors for all phases of an assymetric motion profile with a 2 cm peak-to-peak amplitude were 0.68, 0.67 and 1.85 mm for 60 o 4D DTS, 360 o CBCT and 4DCT, respectively. Motion artifacts for 4D DTS were found to be substantially less than those seen in 4DCT, which is the current clinical standard in 4D imaging.

We then developed a comprehensive framework for relating patient respiratory parameters with acquisition and reconstruction parameters for slow gantry rotation 4D DTS and 4D CBCT imaging. This framework was validated and optimized with phantom and lung patient studies. The framework facilitates calculation of optimal frame rates and gantry rotation speeds based on patient specific respiratory parameters and required temporal resolution (task dependent). We also conducted lung patient studies to investigate required scan angles for 4D DTS and achievable dose and scan times for 4D DTS and 4D CBCT using the optimized framework. This explicit and comprehensive framework of relationships allowed us to demonstrate that under-sampling artifacts can be controlled, and 4D CBCT images can be acquired using lower doses than previously reported. We reconstructed 4D CBCT images of three patients with accumulated doses of 4.8 to 5.7 cGy. These doses are three times less than the doses used for the only previously reported 4D CBCT investigation that did not report images characterized by severe under-sampling artifacts.

We found that scan times for 200 o 4D CBCT imaging using acquisition sequences optimized for reduction of imaging dose and under-sampling artifacts were necessarily between 4 and 7 minutes (depending on patient respiration). The results from lung patient studies concluded that scan times could be reduced using 4D DTS. Patient 4D DTS studies demonstrated that tumor visibility for the lung patients we studied could be achieved using 30 o scan angles for coronal views. Scan times for those cases were between 41 and 50 seconds. Additional dose reductions were also demonstrated. Image doses were between 1.56 and 2.13 cGy. These doses are well below doses for standard CBCT scans. The techniques developed and reported in this thesis demonstrate how respiratory motion can be imaged in the radiotherapy treatment room using clinically feasible imaging doses and scan times.

The developments of highly conformal and precise radiation therapy techniques promote the necessity of more accurate treatment target localization and tracking. On-board imaging techniques, especially the x-ray based techniques, have found a great popularity nowadays for on-board target localization and tracking. With an objective to improve the accuracy of on-board imaging for lung cancer patients, the dissertation work focuses on the investigations of using limited-angle on-board x-ray projections for image guidance. The limited-angle acquisition enables scan time and imaging dose reduction and improves the mechanical clearance of imaging.

First of all, the dissertation developed a phase-matched digital tomosynthesis (DTS) technique using limited-angle (<=30 deg) projections for lung tumor localization. This technique acquires the same traditional motion-blurred on-board DTS image as the 3D-DTS technique, but uses the planning 4D computed tomography (CT) to synthesize a phase-matched reference DTS to register with the on-board DTS for tumor localization. Of the 324 different scenarios simulated using the extended cardiac torso (XCAT) digital phantom, the phase-matched DTS technique localizes the 3D target position with an localization error of 1.07 mm (± 0.57 mm) (average ± standard deviation (S.D.)). Similarly, for the total 60 scenarios evaluated using the computerized imaging reference system (CIRS) 008A physical phantom, the phase-matched DTS technique localizes the 3D target position with an average localization error of 1.24 mm (± 0.87 mm). In addition to the phantom studies, preliminary clinical cases were also studied using imaging data from three lung cancer patients. Using the localization results of 4D cone beam computed tomography (CBCT) as `gold-standard', the phase-matched DTS techniques localized the tumor to an average localization error of 1.5 mm (± 0.5 mm).

The phantom and patient study results show that the phase-matched DTS technique substantially improved the accuracy of moving lung target localization, as compared to the 3D-DTS technique. The phase-matched DTS technique can provide accurate lung target localizations like 4D-DTS, but with much reduced imaging dose and scan time. The phase-matched DTS technique is also found more robust, being minimally affected by variations of respiratory cycle lengths, fractions of respiration cycle contained within the DTS scan and the scan directions, which potentially enables quasi-instantaneous (within a sub-breathing cycle) moving target verification during radiation therapy, preferably arc therapy.

Though the phase-matched DTS technique can provide accurate target localization under normal scenarios, its accuracy is limited when the patient on-board breathing experiences large variations in motion amplitudes. In addition, the limited-angle based acquisition leads to severe structural distortions in DTS images reconstructed by the current clinical gold-standard Feldkamp-Davis-Kress (FDK) reconstruction algorithm, which prohibit accurate target deformation tracking, delineation and dose calculation.

To solve the above issues, the dissertation further developed a prior knowledge based image estimation technique to fundamentally change the landscape of limited-angle based imaging. The developed motion modeling and free-form deformation (MM-FD) method estimates high quality on-board 4D-CBCT images through applying deformation field maps to existing prior planning 4D-CT images. The deformation field maps are solved using two steps: first, a principal component analysis based motion model is built using the planning 4D-CT (motion modeling). The deformation field map is constructed as an optimized linear combination of the extracted motion modes. Second, with the coarse deformation field maps obtained from motion modeling, a further fine-tuning process called free-form deformation is applied to further correct the residual errors from motion modeling. Using the XCAT phantom, a lung patient with a 30 mm diameter tumor was simulated to have various anatomical and respirational variations from the planning 4D-CT to on-board 4D-CBCTs, including respiration amplitude variations, tumor size variations, tumor average position variations, and phase shift between tumor and body respiratory cycles. The tumors were contoured in both the estimated and the `ground-truth' on-board 4D-CBCTs for comparison. 3D volume percentage error (VPE) and center-of-mass error (COME) were calculated to evaluate the estimation accuracy of the MM-FD technique. For all simulated patient scenarios, the average (± S.D.) VPE / COME of the tumor in the prior image without image estimation was 136.11% (± 42.76%) / 15.5 mm (± 3.9 mm). Using orthogonal-view 30 deg scan angle, the average VPE/COME of the tumors in the MM-FD estimated on-board images was substantially reduced to 5.22% (± 2.12%) / 0.5 mm (± 0.4 mm).

In addition to XCAT simulation, CIRS phantom measurements and actual patient studies were also performed. For these clinical studies, we used the normalized cross-correlation (NCC) as a new similarity metric and developed an updated MMFD-NCC method, to improve the robustness of the image estimation technique to the intensity mismatches between CT and CBCT imaging systems. Using 4D-CBCT reconstructed from fully-sampled on-board projections as `gold-standard', for the CIRS phantom study, the average (± S.D.) VPE / COME of the tumor in the prior image and the tumors in the MMFD-NCC estimated images was 257.1% (± 60.2%) / 10.1 mm (± 4.5 mm) and 7.7% (± 1.2%) / 1.2 mm (± 0.2mm), respectively. For three patient cases, the average (± S.D.) VPE / COME of tumors in the prior images and tumors in the MMFD-NCC estimated images was 55.6% (± 45.9%) / 3.8 mm (± 1.9 mm) and 9.6% (± 6.1%) / 1.1 mm (± 0.5 mm), respectively. With the combined benefits of motion modeling and free-form deformation, the MMFD-NCC method has achieved highly accurate image estimation under different scenarios.

Another potential benefit of on-board 4D-CBCT imaging is the on-board dose calculation and verification. Since the MMFD-NCC estimates the on-board 4D-CBCT through deforming prior 4D-CT images, the 4D-CBCT inherently has the same image quality and Hounsfield unit (HU) accuracy as 4D-CT and therefore can potentially improve the accuracy of on-board dose verification. Both XCAT and CIRS phantom studies were performed for the dosimetric study. Various inter-fractional variations featuring patient motion pattern change, tumor size change and tumor average position change were simulated from planning CT to on-board images. The doses calculated on the on-board CBCTs estimated by MMFD-NCC (MMFD-NCC doses) were compared to the doses calculated on the `gold-standard' on-board images (gold-standard doses). The absolute deviations of minimum dose (DDmin), maximum dose (DDmax), mean dose (DDmean) and prescription dose coverage (DV100%) of the planning target volume (PTV) were evaluated. In addition, 4D on-board treatment dose accumulations were performed using 4D-CBCT images estimated by MMFD-NCC in the CIRS phantom study. The accumulated doses were compared to those measured using optically stimulated luminescence (OSL) detectors and radiochromic films.

The MMFD-NCC doses matched very well with the gold-standard doses. For the XCAT phantom study, the average (± S.D.) DDmin, DDmax, DDmean and DV100% (values normalized by the prescription dose or the total PTV volume) between the MMFD-NCC PTV doses and the gold-standard PTV doses were 0.3% (± 0.2%), 0.9% (± 0.6%), 0.6% (± 0.4%) and 1.0% (± 0.8%), respectively. Similarly, for the CIRS phantom study, the corresponding values between the MMFD-NCC PTV doses and the gold-standard PTV doses were 0.4% (± 0.8%), 0.8% (± 1.0%), 0.5% (± 0.4%) and 0.8% (± 0.8%), respectively. For the 4D dose accumulation study, the average (± S.D.) absolute dose deviation (normalized by local doses) between the accumulated doses and the OSL measured doses was 3.0% (± 2.4%). The average gamma index (3%/3mm) between the accumulated doses and the radiochromic film measured doses was 96.1%. The MMFD-NCC estimated 4D-CBCT enables accurate on-board dose calculation and accumulation for lung radiation therapy under different scenarios. It can potentially be valuable for treatment quality assessment and adaptive radiation therapy.

However, a major limitation of the estimated 4D-CBCTs above is that they can only capture inter-fractional patient variations as they were acquired prior to each treatment. The intra-treatment patient variations cannot be captured, which can also affect the treatment accuracy. In light of this issue, an aggregated kilo-voltage (kV) and mega-voltage (MV) imaging scheme was developed to enable intra-treatment imaging. Through using the simultaneously acquired kV and MV projections during the treatment, the MMFD-NCC method enabled 4D-CBCT estimation using combined kV and MV projections.

For all XCAT-simulated patient scenarios, the average (± S.D.) VPE / COME of the tumor in the prior image and tumors in the MMFD-NCC estimated images (using kV + open field MV) was 136.11% (± 42.76%) / 15.5 mm (± 3.9 mm) and 4.5% (± 1.9%) / 0.3 mm (± 0.4 mm), respectively. In contrast, the MMFD-NCC estimation using kV + beam's eye view (BEV) MV projections yielded results of 4.3% (± 1.5%) / 0.3 mm (± 0.3 mm). The kV + BEV MV aggregation can estimate the target as accurately as the kV + open field MV aggregation. The impact of this study is threefold: 1. the kV and MV projections can be acquired at the same time. The imaging time will be cut to half as compared to the cases which use kV projections only. 2. The kV and MV aggregation enables intra-treatment imaging and target tracking, since the MV projections can be the side products of the treatment beams (BEV MV). 3. As the BEV MV projections originate from the treatment beams, there will be no extra MV imaging dose to the patient.

The above introduced 4D-CBCT estimation techniques were all based on limited-angle acquisition. Though limited-angle acquisition enables substantial scan time and dose reduction as compared to the full-angle scan, it is still not real-time and cannot provide `cine' imaging, which refers to the instantaneous imaging with negligible scan time and imaging dose. Cine imaging is important in image guided radiation therapy practice, considering the respirational variations may occur quickly and frequently during the treatment. For instance, the patient may experience a breathing baseline shift after every respiratory cycle. The limited-angle 4D-CBCT approach still requires a scan time of multiple respiratory cycles, which will not be able to capture the baseline shift in a timely manner.

In light of this issue, based on the previously developed MMFD-NCC method, an AI-FD-NCC method was further developed to enable quasi-cine CBCT imaging using extremely limited-angle (<=6 deg) projections. Using pre-treatment 4D-CBCTs acquired just before the treatment as prior information, AI-FD-NCC enforces an additional prior adaptive constraint to estimate high quality `quasi-cine' CBCT images. Two on-board patient scenarios: tumor baseline shift and continuous motion amplitude change were simulated through the XCAT phantom. Using orthogonal-view 6 deg projections, for the baseline shift scenario, the average (± S.D.) VPE / COME of the tumors in the AI-FD-NCC estimated images was 1.3% (± 0.5%) / 0.4 mm (± 0.1 mm). For the amplitude variation scenario, the average (± S.D.) VPE / COME of the tumors in the AI-FD-NCC estimated images was 1.9% (± 1.1%) / 0.5 mm (± 0.2 mm). The impact of this study is three-fold: first, the quasi-cine CBCT technique enables actual real-time volumetric tracking of tumor and normal tissues. Second, the method enables real-time tumor and normal tissues dose calculation and accumulation. Third, the high-quality volumetric images obtained can potentially be used for real-time adaptive radiation therapy.

In summary, the dissertation work uses limited-angle on-board x-ray projections to reconstruct/estimate volumetric images for lung tumor localization, delineation and dose calculation. Limited-angle acquisition reduces imaging dose, scan time and improves imaging mechanical clearance. Using limited-angle projections enables continuous, sub respiratory-cycle tumor localization, as validated in the phase-matched DTS study. The combination of prior information, motion modeling, free-form deformation and limited-angle on-board projections enables high-quality on-board 4D-CBCT estimation, as validated by the MM-FD / MMFD-NCC techniques. The high-quality 4D-CBCT not only can be applied for accurate target localization and delineation, but also can be used for accurate treatment dose verification, as validated in the dosimetric study. Through aggregating the kV and MV projections for image estimation, intra-treatment 4D-CBCT imaging was also proposed and validated for its feasibility. At last, the introduction of more accurate prior information and additional adaptive prior knowledge constraints also enables quasi-cine CBCT imaging using extremely-limited angle projections. The dissertation work contributes to lung on-board imaging in many aspects with various approaches, which can be beneficial to the future lung image guided radiation therapy practice.

Digital tomosynthesis (DTS) is a quasi-three-dimensional (3D) imaging technique which reconstructs images from a limited angle of cone-beam projections with shorter acquisition time, lower imaging dose, and less mechanical constraint than full cone-beam CT (CBCT). However, DTS images reconstructed by the conventional filtered back projection method have low plane-to-plane resolution, and they do not provide full volumetric information for target localization due to the limited angle of the DTS acquisition.

This dissertation presents the optimization and clinical implementation of image guided radiation therapy using limited-angle projections.

A hybrid multiresolution rigid-body registration technique was developed to automatically register reference DTS images with on-board DTS images to guide patient positioning in radiation therapy. This hybrid registration technique uses a faster but less accurate static method to achieve an initial registration, followed by a slower but more accurate adaptive method to fine tune the registration. A multiresolution scheme is employed in the registration to further improve the registration accuracy, robustness and efficiency. Normalized mutual information is selected as the criterion for the similarity measure, and the downhill simplex method is used as the search engine. This technique was tested using image data both from an anthropomorphic chest phantom and from head-and-neck cancer patients. The effects of the scan angle and the region-of-interest size on the registration accuracy and robustness were investigated. The average capture ranges in single-axis simulations with a 44° scan angle and a large ROI covering the entire DTS volume were between -31 and +34 deg for rotations and between -89 and +78 mm for translations in the phantom study, and between -38 and +38 deg for rotations and between -58 and +65 mm for translations in the patient study.

Additionally, a novel limited-angle CBCT estimation method using a deformation field map was developed to optimally estimate volumetric information of organ deformation for soft tissue alignment in image guided radiation therapy. The deformation field map is solved by using prior information, a deformation model, and new projection data. Patients' previous CBCT data are used as the prior information, and the new patient volume to be estimated is considered as a deformation of the prior patient volume. The deformation field is solved by minimizing bending energy and maintaining new projection data fidelity using a nonlinear conjugate gradient method. The new patient CBCT volume is then obtained by deforming the prior patient CBCT volume according to the solution to the deformation field. The method was tested for different scan angles in 2D and 3D cases using simulated and real projections of a Shepp-Logan phantom, liver, prostate and head-and-neck patient data. Hardware acceleration and multiresolution scheme are used to accelerate the 3D estimation process. The accuracy of the estimation was evaluated by comparing organ volume, similarity and pixel value differences between limited-angle CBCT and full-rotation CBCT images. Results showed that the respiratory motion in the liver patient, rectum volume change in the prostate patient, and the weight loss and airway volume change in the head-and-neck patient were accurately estimated in the 60° CBCT images. This new estimation method is able to optimally estimate the volumetric information using 60-degree projection images. It is both technically and clinically feasible for image-guidance in radiation therapy.