Principal component analysis-based imaging angle determination for 3D motion monitoring using single-slice on-board imaging.

Abstract

Through-plane motion introduces uncertainty in three-dimensional (3D) motion monitoring when using single-slice on-board imaging (OBI) modalities such as cine MRI. We propose a principal component analysis (PCA)-based framework to determine the optimal imaging plane to minimize the through-plane motion for single-slice imaging-based motion monitoring. Four-dimensional computed tomography (4DCT) images of eight thoracic cancer patients were retrospectively analyzed. The target volumes were manually delineated at different respiratory phases of 4DCT. We performed automated image registration to establish the 4D respiratory target motion trajectories for all patients. PCA was conducted using the motion information to define the three principal components of the respiratory motion trajectories. Two imaging planes were determined perpendicular to the second and third principal component, respectively, to avoid imaging with the primary principal component of the through-plane motion. Single-slice images were reconstructed from 4DCT in the PCA-derived orthogonal imaging planes and were compared against the traditional AP/Lateral image pairs on through-plane motion, residual error in motion monitoring, absolute motion amplitude error and the similarity between target segmentations at different phases. We evaluated the significance of the proposed motion monitoring improvement using paired t test analysis. The PCA-determined imaging planes had overall less through-plane motion compared against the AP/Lateral image pairs. For all patients, the average through-plane motion was 3.6mm (range: 1.6-5.6mm) for the AP view and 1.7mm (range: 0.6-2.7mm) for the Lateral view. With PCA optimization, the average through-plane motion was 2.5mm (range: 1.3-3.9mm) and 0.6mm (range: 0.2-1.5mm) for the two imaging planes, respectively. The absolute residual error of the reconstructed max-exhale-to-inhale motion averaged 0.7mm (range: 0.4-1.3mm, 95% CI: 0.4-1.1mm) using optimized imaging planes, averaged 0.5mm (range: 0.3-1.0mm, 95% CI: 0.2-0.8mm) using an imaging plane perpendicular to the minimal motion component only and averaged 1.3mm (range: 0.4-2.8mm, 95% CI: 0.4-2.3mm) in AP/Lateral orthogonal image pairs. The root-mean-square error of reconstructed displacement was 0.8mm for optimized imaging planes, 0.6mm for imaging plane perpendicular to the minimal motion component only, and 1.6mm for AP/Lateral orthogonal image pairs. When using the optimized imaging planes for motion monitoring, there was no significant absolute amplitude error of the reconstructed motion (P=0.0988), while AP/Lateral images had significant error (P=0.0097) with a paired t test. The average surface distance (ASD) between overlaid two-dimensional (2D) tumor segmentation at end-of-inhale and end-of-exhale for all eight patients was 0.6±0.2mm in optimized imaging planes and 1.4±0.8mm in AP/Lateral images. The Dice similarity coefficient (DSC) between overlaid 2D tumor segmentation at end-of-inhale and end-of-exhale for all eight patients was 0.96±0.03 in optimized imaging planes and 0.89±0.05 in AP/Lateral images. Both ASD (P=0.034) and DSC (P=0.022) were significantly improved in the optimized imaging planes. Motion monitoring using imaging planes determined by the proposed PCA-based framework had significantly improved performance. Single-slice image-based motion tracking can be used for clinical implementations such as MR image-guided radiation therapy (MR-IGRT).