Characterization of cardiac quiescence from retrospective cardiac computed tomography using a correlation-based phase-to-phase deviation measure

Coronary Computed Tomography Angiography Using an Optimal Acquisition Time Window Based on Heart Rate Determined During Breath-Holding Following Free Breathing

OBJECTIVES

To compare the image quality and radiation dose in coronary computed tomography angiography (CCTA) based on different acquisition time windows corresponding to the heart rate of breath-holding after free breathing.

METHODS

Two hundred patients who underwent CCTA with a basal heart rate between 70 and 85 beats/min were divided into groups A and B, with 100 patients in each group. Patients in groups A and B were scanned with the acquisition time window corresponding to the heart rate determined during a breath hold obtained after free breathing and the basal heart rate during free breathing, respectively. Computed tomography (CT) attenuation values of the coronary artery, signal-to-noise ratio (SNR), and contrast-to-noise ratio (CNR) were calculated. The subjective image scores of the groups were assessed blindly by 2 experienced physicians using a 4-point system, and score consistency was compared using the κ test. The volume CT dose index and dose-length product were recorded for each patient, and the effective dose (ED) was calculated. The Kruskal-Wallis H test was performed to evaluate differences in age, heart rate, and body mass index. A χ2 test was used to evaluate sex differences. An independent-sample t test was employed to compare objective and subjective data such as dose-length product, volume CT dose index, ED, SNR, CNR, and averaged subjective assessment scores. Statistical significance was set at P  < 0.05.

RESULTS

No statistically significant differences occurred in sex, age, or body mass index between patients in group A and group B (all P  > 0.05). No significant differences occurred in the mean CT values, mean SNR values, mean CNR values, or mean subjective scores of CCTA images between the patients in groups A and B ( P  > 0.05). The ED values of the patients in group A were 52.93% lower than those in group B ( P  < 0.001).

CONCLUSION

The radiation dose in CCTA examinations can be significantly reduced while maintaining image quality by narrowing the acquisition time window for breath-holding after free breathing.

Near Real-Time Implementation of An Adaptive Seismocardiography – ECG Multimodal Framework for Cardiac Gating

Objective: Accurate gating for data acquisition of computed tomography (CT) is crucial to obtaining high quality images for diagnosing cardiovascular diseases. To illustrate the feasibility of an optimized cardiac gating strategy, we present a near real-time implementation based on fusing seismocardiography (SCG) and ECG. Methods: The implementation was achieved via integrating commercial hardware and software platforms. Testing was performed on five healthy subjects (age: 24–27; m/f: 4/1) and three cardiac patients (age: 41–71; m/f: 2/1), and compared with baseline quiescence derived from echocardiography. Results: The average latency introduced by computerized processing was 5.1 ms, well within a 100 ms tolerance bounded by data accumulation time for quiescence prediction. The average prediction error associated with conventional ECG-only versus SCG-ECG-based method over all subjects were 59.58 ms and 27.24 ms, respectively. Discussion: The results demonstrate that the multimodal framework can achieve improved quiescence prediction accuracy over the ECG-only-based method in near real-time.

Automated quantification and evaluation of motion artifact on coronary CT angiography images

PURPOSE

This study developed and validated a Motion Artifact Quantification algorithm to automatically quantify the severity of motion artifacts on coronary computed tomography angiography (CCTA) images. The algorithm was then used to develop a Motion IQ Decision method to automatically identify whether a CCTA dataset is of sufficient diagnostic image quality or requires further correction.

METHOD

The developed Motion Artifact Quantification algorithm includes steps to identify the right coronary artery (RCA) regions of interest (ROIs), segment vessel and shading artifacts, and to calculate the motion artifact score (MAS) metric. The segmentation algorithms were verified against ground-truth manual segmentations. The segmentation algorithms were also verified by comparing and analyzing the MAS calculated from ground-truth segmentations and the algorithm-generated segmentations. The Motion IQ Decision algorithm first identifies slices with unsatisfactory image quality using a MAS threshold. The algorithm then uses an artifact-length threshold to determine whether the degraded vessel segment is large enough to cause the dataset to be nondiagnostic. An observer study on 30 clinical CCTA datasets was performed to obtain the ground-truth decisions of whether the datasets were of sufficient image quality. A five-fold cross-validation was used to identify the thresholds and to evaluate the Motion IQ Decision algorithm.

RESULTS

The automated segmentation algorithms in the Motion Artifact Quantification algorithm resulted in Dice coefficients of 0.84 for the segmented vessel regions and 0.75 for the segmented shading artifact regions. The MAS calculated using the automated algorithm was within 10% of the values obtained using ground-truth segmentations. The MAS threshold and artifact-length thresholds were determined by the ROC analysis to be 0.6 and 6.25 mm by all folds. The Motion IQ Decision algorithm demonstrated 100% sensitivity, 66.7% ± 27.9% specificity, and a total accuracy of 86.7% ± 12.5% for identifying datasets in which the RCA required correction. The Motion IQ Decision algorithm demonstrated 91.3% sensitivity, 71.4% specificity, and a total accuracy of 86.7% for identifying CCTA datasets that need correction for any of the three main vessels.

CONCLUSION

The Motion Artifact Quantification algorithm calculated accurate (<10% error) motion artifact scores using the automated segmentation methods. The developed algorithms demonstrated high sensitivity (91.3%) and specificity (71.4%) in identifying datasets of insufficient image quality. The developed algorithms for automatically quantifying motion artifact severity may be useful for comparing acquisition techniques, improving best-phase selection algorithms, and evaluating motion compensation techniques.

An Adaptive Seismocardiography (SCG)-ECG Multimodal Framework for Cardiac Gating Using Artificial Neural Networks

To more accurately trigger data acquisition and reduce radiation exposure of coronary computed tomography angiography (CCTA), a multimodal framework utilizing both electrocardiography (ECG) and seismocardiography (SCG) for CCTA prospective gating is presented. Relying upon a three-layer artificial neural network that adaptively fuses individual ECG- and SCG-based quiescence predictions on a beat-by-beat basis, this framework yields a personalized quiescence prediction for each cardiac cycle. This framework was tested on seven healthy subjects (age: 22-48; m/f: 4/3) and eleven cardiac patients (age: 31-78; m/f: 6/5). Seventeen out of 18 benefited from the fusion-based prediction as compared to the ECG-only-based prediction, the traditional prospective gating method. Only one patient whose SCG was compromised by noise was more suitable for ECG-only-based prediction. On average, our fused ECG-SCG-based method improves cardiac quiescence prediction by 47% over ECG-only-based method; with both compared against the gold standard, B-mode echocardiography. Fusion-based prediction is also more resistant to heart rate variability than ECG-only- or SCG-only-based prediction. To assess the clinical value, the diagnostic quality of the CCTA reconstructed volumes from the quiescence derived from ECG-, SCG- and fusion-based predictions were graded by a board-certified radiologist using a Likert response format. Grading results indicated the fusion-based prediction improved diagnostic quality. ECG may be a sub-optimal modality for quiescence prediction and can be enhanced by the multimodal framework. The combination of ECG and SCG signals for quiescence prediction bears promise for a more personalized and reliable approach than ECG-only-based method to predict cardiac quiescence for prospective CCTA gating.

Evaluation of motion artifact metrics for coronary CT angiography

PURPOSE

This study quantified the performance of coronary artery motion artifact metrics relative to human observer ratings. Motion artifact metrics have been used as part of motion correction and best-phase selection algorithms for Coronary Computed Tomography Angiography (CCTA). However, the lack of ground truth makes it difficult to validate how well the metrics quantify the level of motion artifact. This study investigated five motion artifact metrics, including two novel metrics, using a dynamic phantom, clinical CCTA images, and an observer study that provided ground-truth motion artifact scores from a series of pairwise comparisons.

METHOD

Five motion artifact metrics were calculated for the coronary artery regions on both phantom and clinical CCTA images: positivity, entropy, normalized circularity, Fold Overlap Ratio (FOR), and Low-Intensity Region Score (LIRS). CT images were acquired of a dynamic cardiac phantom that simulated cardiac motion and contained six iodine-filled vessels of varying diameter and with regions of soft plaque and calcifications. Scans were repeated with different gantry start angles. Images were reconstructed at five phases of the motion cycle. Clinical images were acquired from 14 CCTA exams with patient heart rates ranging from 52 to 82 bpm. The vessel and shading artifacts were manually segmented by three readers and combined to create ground-truth artifact regions. Motion artifact levels were also assessed by readers using a pairwise comparison method to establish a ground-truth reader score. The Kendall's Tau coefficients were calculated to evaluate the statistical agreement in ranking between the motion artifacts metrics and reader scores. Linear regression between the reader scores and the metrics was also performed.

RESULTS

On phantom images, the Kendall's Tau coefficients of the five motion artifact metrics were 0.50 (normalized circularity), 0.35 (entropy), 0.82 (positivity), 0.77 (FOR), 0.77(LIRS), where higher Kendall's Tau signifies higher agreement. The FOR, LIRS, and transformed positivity (the fourth root of the positivity) were further evaluated in the study of clinical images. The Kendall's Tau coefficients of the selected metrics were 0.59 (FOR), 0.53 (LIRS), and 0.21 (Transformed positivity). In the study of clinical data, a Motion Artifact Score, defined as the product of FOR and LIRS metrics, further improved agreement with reader scores, with a Kendall's Tau coefficient of 0.65.

CONCLUSION

The metrics of FOR, LIRS, and the product of the two metrics provided the highest agreement in motion artifact ranking when compared to the readers, and the highest linear correlation to the reader scores. The validated motion artifact metrics may be useful for developing and evaluating methods to reduce motion in Coronary Computed Tomography Angiography (CCTA) images.

Seismocardiography-Based Cardiac Computed Tomography Gating Using Patient-Specific Template Identification and Detection

To more accurately trigger cardiac computed tomography angiography (CTA) than electrocardiography (ECG) alone, a sub-system is proposed as an intermediate step toward fusing ECG with seismocardiography (SCG). Accurate prediction of quiescent phases is crucial to prospectively gating CTA, which is susceptible to cardiac motion and, thus, can affect the diagnostic quality of images. The key innovation of this sub-system is that it identifies the SCG waveform corresponding to heart sounds and determines their phases within the cardiac cycles. Furthermore, this relationship is modeled as a linear function with respect to heart rate. For this paper, B-mode echocardiography is used as the gold standard for identifying the quiescent phases. We analyzed synchronous ECG, SCG, and echocardiography data acquired from seven healthy subjects (mean age: 31; age range: 22–48; males: 4) and 11 cardiac patients (mean age: 56; age range: 31–78; males: 6). On average, the proposed algorithm was able to successfully identify 79% of the SCG waveforms in systole and 68% in diastole. The simulated results show that SCG-based prediction produced less average phase error than that of ECG. It was found that the accuracy of ECG-based gating is more susceptible to increases in heart rate variability, while SCG-based gating is susceptible to high cycle to cycle variability in morphology. This pilot work of prediction using SCG waveforms enriches the framework of a comprehensive system with multiple modalities that could potentially, in real time, improve the image quality of CTA.

Automated selection of the optimal cardiac phase for single-beat coronary CT angiography reconstruction

PURPOSE

Reconstructing a low-motion cardiac phase is expected to improve coronary artery visualization in coronary computed tomography angiography (CCTA) exams. This study developed an automated algorithm for selecting the optimal cardiac phase for CCTA reconstruction. The algorithm uses prospectively gated, single-beat, multiphase data made possible by wide cone-beam imaging. The proposed algorithm differs from previous approaches because the optimal phase is identified based on vessel image quality (IQ) directly, compared to previous approaches that included motion estimation and interphase processing. Because there is no processing of interphase information, the algorithm can be applied to any sampling of image phases, making it suited for prospectively gated studies where only a subset of phases are available.

METHODS

An automated algorithm was developed to select the optimal phase based on quantitative IQ metrics. For each reconstructed slice at each reconstructed phase, an image quality metric was calculated based on measures of circularity and edge strength of through-plane vessels. The image quality metric was aggregated across slices, while a metric of vessel-location consistency was used to ignore slices that did not contain through-plane vessels. The algorithm performance was evaluated using two observer studies. Fourteen single-beat cardiac CT exams (Revolution CT, GE Healthcare, Chalfont St. Giles, UK) reconstructed at 2% intervals were evaluated for best systolic (1), diastolic (6), or systolic and diastolic phases (7) by three readers and the algorithm. Pairwise inter-reader and reader-algorithm agreement was evaluated using the mean absolute difference (MAD) and concordance correlation coefficient (CCC) between the reader and algorithm-selected phases. A reader-consensus best phase was determined and compared to the algorithm selected phase. In cases where the algorithm and consensus best phases differed by more than 2%, IQ was scored by three readers using a five point Likert scale.

RESULTS

There was no statistically significant difference between inter-reader and reader-algorithm agreement for either MAD or CCC metrics (p > 0.1). The algorithm phase was within 2% of the consensus phase in 15/21 of cases. The average absolute difference between consensus and algorithm best phases was 2.29% ± 2.47%, with a maximum difference of 8%. Average image quality scores for the algorithm chosen best phase were 4.01 ± 0.65 overall, 3.33 ± 1.27 for right coronary artery (RCA), 4.50 ± 0.35 for left anterior descending (LAD) artery, and 4.50 ± 0.35 for left circumflex artery (LCX). Average image quality scores for the consensus best phase were 4.11 ± 0.54 overall, 3.44 ± 1.03 for RCA, 4.39 ± 0.39 for LAD, and 4.50 ± 0.18 for LCX. There was no statistically significant difference (p > 0.1) between the image quality scores of the algorithm phase and the consensus phase.

CONCLUSIONS

The proposed algorithm was statistically equivalent to a reader in selecting an optimal cardiac phase for CCTA exams. When reader and algorithm phases differed by >2%, image quality as rated by blinded readers was statistically equivalent. By detecting the optimal phase for CCTA reconstruction, the proposed algorithm is expected to improve coronary artery visualization in CCTA exams.

Echocardiography as an indication of continuous-time cardiac quiescence

Cardiac computed tomography (CT) angiography using prospective gating requires that data be acquired during intervals of minimal cardiac motion to obtain diagnostic images of the coronary vessels free of motion artifacts. This work is intended to assess B-mode echocardiography as a continuous-time indication of these quiescent periods to determine if echocardiography can be used as a cost-efficient, non-ionizing modality to develop new prospective gating techniques for cardiac CT. These new prospective gating approaches will not be based on echocardiography itself but on CT-compatible modalities derived from the mechanics of the heart (e.g. seismocardiography and impedance cardiography), unlike the current standard electrocardiogram. To this end, echocardiography and retrospectively-gated CT data were obtained from ten patients with varied cardiac conditions. CT reconstructions were made throughout the cardiac cycle. Motion of the interventricular septum (IVS) was calculated from both echocardiography and CT reconstructions using correlation-based, deviation techniques. The IVS was chosen because it (1) is visible in echocardiography images, whereas the coronary vessels generally are not, and (2) has been shown to be a suitable indicator of cardiac quiescence. Quiescent phases were calculated as the minima of IVS motion and CT volumes were reconstructed for these phases. The diagnostic quality of the CT reconstructions from phases calculated from echocardiography and CT data was graded on a four-point Likert scale by a board-certified radiologist fellowship-trained in cardiothoracic radiology. Using a Wilcoxon signed-rank test, no significant difference in the diagnostic quality of the coronary vessels was found between CT volumes reconstructed from echocardiography- and CT-selected phases. Additionally, there was a correlation of 0.956 between the echocardiography- and CT-selected phases. This initial work suggests that B-mode echocardiography can be used as a tool to develop CT-compatible gating techniques based on modalities derived from cardiac mechanics rather than relying on the ECG alone.

Relationship between cardiac quiescent periods derived from seismocardiography and echocardiography

The seismocardiogram (SCG) is a measure of chest wall acceleration due to cardiac motion that could potentially supplement the electrocardiogram (ECG) to more reliably predict cardiac quiescence. Accurate prediction is critical for modalities requiring minimal motion during imaging data acquisition, such as cardiac computed tomography (CT) and magnetic resonance imaging (MRI). For seven healthy subjects, SCG and B-mode echocardiography were used to identify quiescent periods on a beat-by-beat basis. Quiescent periods were detected as time intervals when the magnitude of the velocity signals calculated from SCG and echocardiography were less than a specified threshold. The quiescent periods detected from SCG were compared to those detected from B-mode echocardiography. The quiescent periods of the SCG were found to occur before those detected by echocardiography. A linear relationship between the delay from SCG- to echocardiography-detected phases with respect to heart rate was found. This delay could potentially be used to predict cardiac quiescence from SCG-observed quiescence for use with cardiac imaging modalities such as CT and MRI.

Seismocardiography-Based Detection of Cardiac Quiescence

Cardiac-computed tomography angiography (CTA) is a minimally invasive imaging technology for characterizing coronary arteries. A fundamental limitation of CTA imaging is cardiac movement, which can cause artifacts and reduce the quality of the obtained images. To mitigate this problem, current approaches involve gating the image based on the electrocardiogram (ECG) to predict the timing of quiescent periods of the cardiac cycle. This paper focuses on developing a foundation for using a mechanical alternative to the ECG for finding these quiescent periods: the seismocardiogram (SCG). SCG was used to determine beat-by-beat systolic and diastolic quiescent periods of the cardiac cycle for nine healthy subjects, and 11 subjects with various cardiovascular diseases. To reduce noise in the SCG, and quantify these quiescent periods, a Kalman filter was designed to extract the velocity of chest wall movement from the recorded SCG signals. The average systolic and diastolic quiescent periods were centered at 29% and 76% for the healthy subjects, and 33% and 79% for subjects with cardiovascular disease. Both inter and intrasubject variability in the quiescent phases were observed compared to ECG-predicted phases, suggesting that the ECG may be a suboptimal modality for predicting quiescence, and that the SCG provides complementary data to the ECG.