Feature Guided Motion Artifact Reduction with Structure-Awareness in 4D CT Images

Rapid estimation of 4DCT motion-artifact severity based on 1D breathing-surrogate periodicity

PURPOSE

Motion artifacts are common in patient four-dimensional computed tomography (4DCT) images, leading to an ill-defined tumor volume with large variations for radiotherapy treatment and a poor foundation with low imaging fidelity for studying respiratory motion. The authors developed a method to estimate 4DCT image quality by establishing a correlation between the severity of motion artifacts in 4DCT images and the periodicity of the corresponding 1D respiratory waveform (1DRW) used for phase binning in 4DCT reconstruction.

METHODS

Discrete Fourier transformation (DFT) was applied to analyze 1DRW periodicity. The breathing periodicity index (BPI) was defined as the sum of the largest five Fourier coefficients, ranging from 0 to 1. Distortional motion artifacts (excluding blurring) of cine-scan 4DCT at the junctions of adjacent couch positions around the diaphragm were classified in three categories: incomplete, overlapping, and duplicate anatomies. To quantify these artifacts, discontinuity of the diaphragm at the junctions was measured in distance and averaged along six directions in three orthogonal views. Artifacts per junction (APJ) across the entire diaphragm were calculated in each breathing phase and phase-averaged APJ¯, defined as motion-artifact severity (MAS), was obtained for each patient. To make MAS independent of patient-specific motion amplitude, two new MAS quantities were defined: MAS(D) is normalized to the maximum diaphragmatic displacement and MAS(V) is normalized to the mean diaphragmatic velocity (the breathing period was obtained from DFT analysis of 1DRW). Twenty-six patients' free-breathing 4DCT images and corresponding 1DRW data were studied.

RESULTS

Higher APJ values were found around midventilation and full inhalation while the lowest APJ values were around full exhalation. The distribution of MAS is close to Poisson distribution with a mean of 2.2 mm. The BPI among the 26 patients was calculated with a value ranging from 0.25 to 0.93. The DFT calculation was within 3 s per 1DRW. Correlations were found between 1DRW periodicity and 4DCT artifact severity: -0.71 for MAS(D) and -0.73 for MAS(V). A BPI greater than 0.85 in a 1DRW suggests minimal motion artifacts in the corresponding 4DCT images.

CONCLUSIONS

The breathing periodicity index and motion-artifact severity index are introduced to assess the relationship between 1DRW and 4DCT. A correlation between 1DRW periodicity and 4DCT artifact severity has been established. The 1DRW periodicity provides a rapid means to estimate 4DCT image quality. The rapid 1DRW analysis and the correlative relationship can be applied prospectively to identify irregular breathers as candidates for breath coaching prior to 4DCT scan and retrospectively to select high-quality 4DCT images for clinical motion-management research.

Reconstruction of four-dimensional computed tomography lung images by applying spatial and temporal anatomical constraints using a Bayesian model

Current four-dimensional computed tomography (4-D CT) lung image reconstruction methods rely on respiratory gating, such as surrogate, to sort the large number of axial images captured during multiple breathing cycles into serial three-dimensional CT images of different respiratory phases. Such sorting methods may be subject to external surrogate signal noises due to poor reproducibility of breathing cycles. New image-matching-based reconstruction algorithms refine the 4-D CT reconstruction by matching neighboring image slices, and they generally work better for the cine mode of 4-D CT acquisition than the helical mode due to different table positions of axial images in the helical mode. We propose a Bayesian model (BM) based automated 4-D CT lung image reconstruction for helical mode scans. BM allows for applying new spatial and temporal anatomical constraints in the optimization procedure. Using an iterative optimization procedure, each axial image is assigned to a respiratory phase to make sure the anatomical structures are spatially and temporally smooth based on the BM framework. In experiments, we visually and quantitatively compared the results of the proposed BM-based 4-D CT reconstruction with the respiratory surrogate and the normalized cross-correlation based image matching method using both simulated and actual 4-D patient scans. The results indicated that the proposed algorithm yielded more accurate reconstruction and fewer artifacts in the 4-D CT image series.

Assessment of a quantitative metric for 4D CT artifact evaluation by observer consensus

The benefits of four-dimensional computed tomography (4D CT) are limited by the presence of artifacts that remain difficult to quantify. A correlation-based metric previously proposed for ciné 4D CT artifact identification was further validated as an independent artifact evaluator by using a novel qualitative assessment featuring a group of observers reaching a consensus decision on artifact location and magnitude. The consensus group evaluated ten ciné 4D CT scans for artifacts over each breathing phase of coronal lung views assuming one artifact per couch location. Each artifact was assigned a magnitude score of 1-5, 1 indicating lowest severity and 5 indicating highest severity. Consensus group results served as the ground truth for assessment of the correlation metric. The ten patients were split into two cohorts; cohort 1 generated an artifact identification threshold derived from receiver operating characteristic analysis using the Youden Index, while cohort 2 generated sensitivity and specificity values from application of the artifact threshold. The Pearson correlation coefficient was calculated between the correlation metric values and the consensus group scores for both cohorts. The average sensitivity and specificity values found with application of the artifact threshold were 0.703 and 0.476, respectively. The correlation coefficients of artifact magnitudes for cohort 1 and 2 were 0.80 and 0.61, respectively, (p < 0.001 for both); these correlation coefficients included a few scans with only two of the five possible magnitude scores. Artifact incidence was associated with breathing phase (p < 0.002), with presentation less likely near maximum exhale. Overall, the correlation metric allowed accurate and automated artifact identification. The consensus group evaluation resulted in efficient qualitative scoring, reduced interobserver variation, and provided consistent identification of artifact location and magnitudes.

A Parameterization Based Correspondence Method for PDM Building

Place-march of corresponding landmarks is one of the major factors influencing 3D Points Distribution Model (PDM) quality. In this study, we propose a semi-automatic correspondence method based on surface parameterization theory. All the training sets are mapped into a spherical domain previously. The rotation transformation of training samples is regarded as spherical rotation of their maps. We solve it by comparing the density distribution of surface map of training sample with respect to the reference model. Simultaneously, the corresponding landmarks across the whole training set are marketed depending on the spherical coordinates on parameter domain. In this paper, we also compared the corresponding results with two constraint conditions of spherical conformal mapping: 3 datum points constrain and zero-mass constrain. Experimental results are given for left lung training sets of 3D shapes. The mean result with the 3 datum points constraint and the zero mass-center constraint was 21.65 mm and 20.19 mm respectively.

Helical mode lung 4D-CT reconstruction using Bayesian model

4D computed tomography (CT) has been widely used for treatment planning of thoracic and abdominal cancer radiotherapy. Current 4D-CT lung image reconstruction methods rely on respiratory gating to rearrange the large number of axial images into different phases, which may be subject to external surrogate errors due to poor reproducibility of breathing cycles. New image-matching-based reconstruction works better for the cine mode of 4D-CT acquisition than the helical mode because the table position of each axial image is different in helical mode and image matching might suffer from bigger errors. In helical mode, not only the phases but also the un-uniform table positions of images need to be considered. We propose a Bayesian method for automated 4D-CT lung image reconstruction in helical mode 4D scans. Each axial image is assigned to a respiratory phase based on the Bayesian framework that ensures spatial and temporal smoothness of surfaces of anatomical structures. Iterative optimization is used to reconstruct a series of 3D-CT images for subjects undergoing 4D scans. In experiments, we compared visually and quantitatively the results of the proposed Bayesian 4D-CT reconstruction algorithm with the respiratory surrogate and the image matching-based method. The results showed that the proposed algorithm yielded better 4D-CT for helical scans.