Computerized assessment of motion-contaminated calcified plaques in cardiac multidetector CT

Can nontriggered thoracic CT be used for coronary artery calcium scoring? A phantom study

PURPOSE

Coronary artery calcium score, traditionally based on electrocardiography (ECG)-triggered computed tomography (CT), predicts cardiovascular risk. However, nontriggered CT is extensively utilized. The study-purpose is to evaluate the in vitro agreement in coronary calcium score between nontriggered thoracic CT and ECG-triggered cardiac CT.

METHODS

Three artificial coronary arteries containing calcifications of different densities (high, medium, and low), and sizes (large, medium, and small), were studied in a moving cardiac phantom. Two 64-detector CT systems were used. The phantom moved at 0-90 mm∕s in nontriggered low-dose CT as index test, and at 0-30 mm∕s in ECG-triggered CT as reference. Differences in calcium scores between nontriggered and ECG-triggered CT were analyzed by t-test and 95% confidence interval. The sensitivity to detect calcification was calculated as the percentage of positive calcium scores.

RESULTS

Overall, calcium scores in nontriggered CT were not significantly different to those in ECG-triggered CT (p>0.05). Calcium scores in nontriggered CT were within the 95% confidence interval of calcium scores in ECG-triggered CT, except predominantly at higher velocities (≥50 mm∕s) for the high-density and large-size calcifications. The sensitivity for a nonzero calcium score was 100% for large calcifications, but 46%±11% for small calcifications in nontriggered CT.

CONCLUSIONS

When performing multiple measurements, good agreement in positive calcium scores is found between nontriggered thoracic and ECG-triggered cardiac CT. Agreement decreases with increasing coronary velocity. From this phantom study, it can be concluded that a high calcium score can be detected by nontriggered CT, and thus, that nontriggered CT likely can identify individuals at high risk of cardiovascular disease. On the other hand, a zero calcium score in nontriggered CT does not reliably exclude coronary calcification.

Automated tracing of the adventitial contour of aortoiliac and peripheral arterial walls in CT angiography (CTA) to allow calculation of non-calcified plaque burden

Aortoiliac and lower extremity arterial atherosclerotic plaque burden is a risk factor for the development of visceral and peripheral ischemic and aneurismal vascular disease. While prior research allows automated quantification of calcified plaque in these body regions using CT angiograms, no automated method exists to quantify soft plaque. We developed an automatic algorithm that defines the outer wall contour and wall thickness of vessels to quantify non-calcified plaque in CT angiograms of the chest, abdomen, pelvis, and lower extremities. The algorithm encodes the search space as a constrained graph and calculates the outer wall contour by deriving a minimum cost path through the graph, following the visible outer wall contour while minimizing path tortuosity. Our algorithm was statistically equivalent to a reference standard made by two reviewers. Absolute error was 1.9 ± 2.3% compared to the inter-observer variability of 3.9 ± 3.6%. Wall thickness in vessels with atherosclerosis was 3.4 ± 1.6 mm compared to 1.2 ± 0.4 mm in normal vessels. The algorithm shows promise as a tool for quantification of non-calcified plaque in CT angiography. When combined with previous research, our method has the potential to quantify both non-calcified and calcified plaque in all clinically significant systemic arteries, from the thoracic aorta to the arteries of the calf, over a wide range of diameters. This algorithm has the potential to enable risk stratification of patients and facilitate investigations into the relationships between asymptomatic atherosclerosis and a variety of behavioral, physiologic, pathologic, and genotypic conditions.

Computerized method for evaluating diagnostic image quality of calcified plaque images in cardiac CT: validation on a physical dynamic cardiac phantom

PURPOSE

In cardiac computed tomography (CT), important clinical indices, such as the coronary calcium score and the percentage of coronary artery stenosis, are often adversely affected by motion artifacts. As a result, the expert observer must decide whether or not to use these indices during image interpretation. Computerized methods potentially can be used to assist in these decisions. In a previous study, an artificial neural network (ANN) regression model provided assessability (image quality) indices of calcified plaque images from the software NCAT phantom that were highly agreeable with those provided by expert observers. The method predicted assessability indices based on computer-extracted features of the plaque. In the current study, the ANN-predicted assessability indices were used to identify calcified plaque images with diagnostic calcium scores (based on mass) from a physical dynamic cardiac phantom. The basic assumption was that better quality images were associated with more accurate calcium scores.

METHODS

A 64-channel CT scanner was used to obtain 500 calcified plaque images from a physical dynamic cardiac phantom at different heart rates, cardiac phases, and plaque locations. Two expert observers independently provided separate sets of assessability indices for each of these images. Separate sets of ANN-predicted assessability indices tailored to each observer were then generated within the framework of a bootstrap resampling scheme. For each resampling iteration, the absolute calcium score error between the calcium scores of the motion-contaminated plaque image and its corresponding stationary image served as the ground truth in terms of indicating images with diagnostic calcium scores. The performances of the ANN-predicted and observer-assigned indices in identifying images with diagnostic calcium scores were then evaluated using ROC analysis.

RESULTS

Assessability indices provided by the first observer and the corresponding ANN performed similarly (AUC(OBS1) = 0.80 [0.73, 0.86] vs AUC(ANN1) = 0.88 [0.82, 0.92]) as that of the second observer and the corresponding ANN (AUC(OBS2) = 0.87 [0.83,0.91] vs. AUC(ANN2) = 0.90 [0.85, 0.94]). Moreover, the ANN-predicted indices were generated in a fraction of the time required to obtain the observer-assigned indices.

CONCLUSIONS

ANN-predicted assessability indices performed similar to observer-assigned assessability indices in identifying images with diagnostic calcium scores from the physical dynamic cardiac phantom. The results of this study demonstrate the potential of using computerized methods for identifying images with diagnostic clinical indices in cardiac CT images.

Anniversary paper: History and status of CAD and quantitative image analysis: the role of Medical Physics and AAPM

The roles of physicists in medical imaging have expanded over the years, from the study of imaging systems (sources and detectors) and dose to the assessment of image quality and perception, the development of image processing techniques, and the development of image analysis methods to assist in detection and diagnosis. The latter is a natural extension of medical physicists' goals in developing imaging techniques to help physicians acquire diagnostic information and improve clinical decisions. Studies indicate that radiologists do not detect all abnormalities on images that are visible on retrospective review, and they do not always correctly characterize abnormalities that are found. Since the 1950s, the potential use of computers had been considered for analysis of radiographic abnormalities. In the mid-1980s, however, medical physicists and radiologists began major research efforts for computer-aided detection or computer-aided diagnosis (CAD), that is, using the computer output as an aid to radiologists-as opposed to a completely automatic computer interpretation-focusing initially on methods for the detection of lesions on chest radiographs and mammograms. Since then, extensive investigations of computerized image analysis for detection or diagnosis of abnormalities in a variety of 2D and 3D medical images have been conducted. The growth of CAD over the past 20 years has been tremendous-from the early days of time-consuming film digitization and CPU-intensive computations on a limited number of cases to its current status in which developed CAD approaches are evaluated rigorously on large clinically relevant databases. CAD research by medical physicists includes many aspects-collecting relevant normal and pathological cases; developing computer algorithms appropriate for the medical interpretation task including those for segmentation, feature extraction, and classifier design; developing methodology for assessing CAD performance; validating the algorithms using appropriate cases to measure performance and robustness; conducting observer studies with which to evaluate radiologists in the diagnostic task without and with the use of the computer aid; and ultimately assessing performance with a clinical trial. Medical physicists also have an important role in quantitative imaging, by validating the quantitative integrity of scanners and developing imaging techniques, and image analysis tools that extract quantitative data in a more accurate and automated fashion. As imaging systems become more complex and the need for better quantitative information from images grows, the future includes the combined research efforts from physicists working in CAD with those working on quantitative imaging systems to readily yield information on morphology, function, molecular structure, and more-from animal imaging research to clinical patient care. A historical review of CAD and a discussion of challenges for the future are presented here, along with the extension to quantitative image analysis.

Computerized assessment of motion-contaminated calcified plaques in cardiac multidetector CT

An automated method for evaluating the image quality of calcified plaques with respect to motion artifacts in noncontrast-enhanced cardiac computed tomography (CT) images is introduced. This method involves using linear regression (LR) and artificial neural network (ANN) regression models for predicting two patient-specific, region-of-interest-specific, reconstruction-specific and temporal phase-specific image quality indices. The first is a plaque motion index, which is derived from the actual trajectory of the calcified plaque and is represented on a continuous scale. The second is an assessability index, which reflects the degree to which a calcified plaque is affected by motion artifacts, and is represented on an ordinal five-point scale. Two sets of assessability indices were provided independently by two radiologists experienced in evaluating cardiac CT images. Inputs for the regression models were selected from 12 features characterizing the dynamic, morphological, and intensity-based properties of the calcified plaques. Whereas LR-velocity (LR-V) used only a single feature (three-dimensional velocity), the LR-multiple (LR-M) and ANN regression models used the same subset of these 12 features selected through stepwise regression. The regression models were parameterized and evaluated using a database of simulated calcified plaque images from the dynamic NCAT phantom involving nine heart rate/multi-sector gating combinations and 40 cardiac phases covering two cardiac cycles. Six calcified plaques were used for the plaque motion indices and three calcified plaques were used for both sets of assessability indices. In one configuration, images from the second cardiac cycle were used for feature selection and regression model parameterization, whereas images from the first cardiac cycle were used for testing. With this configuration, repeated measures concordance correlation coefficients (CCCs) and associated 95% confidence intervals for the LR-V, LR-M, and ANN were 0.817 [0.785, 0.848], 0.894 [0.869, 0.916], and 0.917 [0.892, 0.936] for the plaque motion indices. For the two sets of assess-ability indices, CCC values for the ANN model were 0.843 [0.791, 0.877] and 0.793 [0.747, 0.828]. These two CCC values were statistically greater than the CCC value of 0.689 [0.648, 0.727], which was obtained by comparing the two sets of assessability indices with each other. These preliminary results suggest that the variabilities of assessability indices provided by regression models can lie within the variabilities of the indices assigned by independent observers. Thus, the potential exists for using regression models and assessability indices for determining optimal phases for cardiac CT image interpretation.