4D Cone-beam CT reconstruction using a motion model based on principal component analysis

• Published in: Medical physics (Lancaster) 2011-12, Vol.38 (12), p.6697-6709
• Main Authors: Staub, David, Docef, Alen, Brock, Robert S., Vaman, Constantin, Murphy, Martin J.
• Format: Article
• Language: English
• Subjects: 4D cone-beam CT; Algorithms; Computed tomography; computerised tomography; Cone beam computed tomography; Cone-Beam Computed Tomography - instrumentation; Cone-Beam Computed Tomography - methods; Digital radiography; Eigenvalues; eigenvalues and eigenfunctions; image reconstruction; image registration; Imaging, Three-Dimensional - methods; medical image processing; Medical image reconstruction; Medical imaging; Motion; motion modeling; Normal modes; numerical analysis; Numerical approximation and analysis; Numerical modeling; optimisation; Phantoms, Imaging; Pneumodyamics, respiration; pneumodynamics; Principal Component Analysis; Probability theory, stochastic processes, and statistics; projection matching; Radiation Imaging Physics; Radiographic Image Enhancement - methods; Radiographic Image Interpretation, Computer-Assisted - methods; Reconstruction; Registration; Reproducibility of Results; Respiratory-Gated Imaging Techniques - methods; Sensitivity and Specificity; time series; Time series analysis; Trajectory models
• ISSN: 0094-2405; 2473-4209; 0094-2405
• DOI: 10.1118/1.3662895

Purpose: To provide a proof of concept validation of a novel 4D cone-beam CT (4DCBCT) reconstruction algorithm and to determine the best methods to train and optimize the algorithm. Methods: The algorithm animates a patient fan-beam CT (FBCT) with a patient specific parametric motion model in order to generate a time series of deformed CTs (the reconstructed 4DCBCT) that track the motion of the patient anatomy on a voxel by voxel scale. The motion model is constrained by requiring that projections cast through the deformed CT time series match the projections of the raw patient 4DCBCT. The motion model uses a basis of eigenvectors that are generated via principal component analysis (PCA) of a training set of displacement vector fields (DVFs) that approximate patient motion. The eigenvectors are weighted by a parameterized function of the patient breathing trace recorded during 4DCBCT. The algorithm is demonstrated and tested via numerical simulation. Results: The algorithm is shown to produce accurate reconstruction results for the most complicated simulated motion, in which voxels move with a pseudo-periodic pattern and relative phase shifts exist between voxels. The tests show that principal component eigenvectors trained on DVFs from a novel 2D/3D registration method give substantially better results than eigenvectors trained on DVFs obtained by conventionally registering 4DCBCT phases reconstructed via filtered backprojection. Conclusions: Proof of concept testing has validated the 4DCBCT reconstruction approach for the types of simulated data considered. In addition, the authors found the 2D/3D registration approach to be our best choice for generating the DVF training set, and the Nelder-Mead simplex algorithm the most robust optimization routine.