Influence of convolution filtering on coronary plaque attenuation values: observations in an ex vivo model of multislice computed tomography coronary angiography.

Attenuation variability (measured in Hounsfield Units, HU) of human coronary plaques using multislice computed tomography (MSCT) was evaluated in an ex vivo model with increasing convolution kernels. MSCT was performed in seven ex vivo left coronary arteries sunk into oil followingthe instillation of saline (1/infinity) and a 1/50 solution of contrast material (400 mgI/ml iomeprol). Scan parameters were: slices/collimation, 16/0.75 mm; rotation time, 375 ms. Four convolution kernels were used: b30f-smooth, b36f-medium smooth, b46f-medium and b60f-sharp. An experienced radiologist scored for the presence of plaques and measured the attenuation in lumen, calcified and noncalcified plaques and the surrounding oil. The results were compared by the ANOVA test and correlated with Pearson's test. The signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR) were calculated. The mean attenuation values were significantly different between the four filters (p < 0.0001) in each structure with both solutions. After clustering for the filter, all of the noncalcified plaque values (20.8 +/- 39.1, 14.2 +/- 35.8, 14.0 +/- 32.0, 3.2 +/- 32.4 HU with saline; 74.7 +/- 66.6, 68.2 +/- 63.3, 66.3 +/- 66.5, 48.5 +/- 60.0 HU in contrast solution) were significantly different, with the exception of the pair b36f-b46f, for which a moderate-high correlation was generally found. Improved SNRs and CNRs were achieved by b30f and b46f. The use of different convolution filters significantly modifief the attenuation values, while sharper filtering increased the calcified plaque attenuation and reduced the noncalcified plaque attenuation.

Introduction

Multislice computed tomography (MSCT) of the coronary arteries is becoming an established method to exclude the presence of significant stenosis because of the high negative predictive value [1-6]. Stenosis, however, occurs in a later stage of coronary artery disease (CAD); consequently, the atherosclerotic plaque composition and morphology may represent a better predictor of plaque stability [7-9]. Until recently, only invasive techniques, such as intracoronary ultrasound (ICUS), have been able to provide consistently reproducible quantitative and qualitative information on plaque composition [10]. ICUS, however, cannot be used for routine evaluation of plaque characteristics because of its invasiveness and related increased risk, additional time and cost.

Magnetic resonance imaging may be used to characterize in vivo the atherosclerotic plaque morphology and composition in carotid [11] and coronary arteries [12, 13].

Several recent in vivo studies have stressed that MSCT can be used to characterize coronary atherosclerotic plaque and has the potential to become the noninvasive modality of choice for plaque imaging [14-17]. This technique is able to provide information on noncalcified coronary plaques, which may be useful to identify a high-risk plaque [18-20].

MSCT plaque assessment is based upon the variable X-ray attenuation of the tissue components. However, several parameters, such as lumen attenuation, convolution filtering, body mass index of the patient and contrast-to-noise ratio (CNR) of the images, are able to modify the attenuation values that are being used to define the composition of coronary plaques.

The aim of the present study was to investigate the effect of increasing convolution filtering on plaque components with varying coronary attenuation in an ex vivo coronary model.

Material and methods

Specimens

Seven ex vivo left anterior descending (LAD) coronary arteries were obtained at autopsy (Fig. 1) from seven deceased patients: four male and three female, with a mean age of 68 ± 10 years: five (four males and one female with an age between 62 and 71 years) jad died of noncardiovascular diseases, and two (one male and one female with an age of 65 and 73 years, respectively) had died of ischemic heart disease. The Institutional Review Board approved the study protocol.

Fig. 1

Example of proximal left coronary artery specimen. CX Left circumflex, LAD left anterior descending artery, LM left main coronary artery

The specimens were prepared and examined separately. Each coronary artery was fitted with two cannulas fixed with surgical thread in the proximal (in the left main) and distal ends (in the LAD). The circumflex artery had been closed earlier at its end with thread. Only the LAD was used because of the major length (segments considered with a mean length of 7.57 cm) and higher prevalence of atherosclerotic disease and for technical reasons (ability to fix a surgical thread in the left main, without hampering the evaluation of the proximal LAD).

Contrast material

Two solutions were used: a saline and a 1/50 dilution of contrast material (400 mgI/ml Iomeprol; Bracco, Italy). The attenuation values (Hounsfield Units, HU) of the two solutions, measured in a 10-ml syringe after dilution, were 15.9 ± 1.8 HU (1/∞; defined as saline; no contrast material was diluted) and 414.8 ± 5.6 HU (1/50; defined as contrast solution). We preferred to use a 1/50 dilution to achieve a mean coronary in-lumen density >300 HU – which is not typically achieved in in vivo cardiac computed tomography angiography (CTA) studies – in order to emphasize the difference with the saline solution. Our study benefits from two consecutive scans of the same vessel (e.g. with and without contrast media) for a comparative evaluation of plaques, which is not truly feasible during in vivo studies due to high radiation exposure.

Experimental settings

A box was filled with olive oil. Prior to positioning the specimen in the oil, saline was injected through the sheaths to wash out the air bubbles in the lumen as much as possible.

Once the specimen had sunk into the oil, which simulated epicardial fat, the solution was instilled through the sheaths using a 10-ml syringe from the proximal end of the specimen. The injection was finished when the solution was observed leaking out of the distal end of specimen. The leaking solution was removed from the specimen using a syringe. The same procedure was used to fill the LAD with both solutions. The specimen was kept for all investigations in the same longitudinal position (head-LM-to feet-distal end of LAD ) to obtain direct cross-section images.

Scan parameters

A MSCT scan (Somatom Sensation 16; Siemens, Germany) was performed following the intra-coronary injection of two solutions.

Scans were performed at the following parameters: slices/collimation, 16/0.75 mm; rotation time, 375 ms; feed per rotation 3.0 mm (pitch: 0.25 mm); 120 kv; 400 mAs; effective slice thickness, 1 mm; reconstruction increment, 0.5 mm; field of view (FOV), 100 mm. Four convolution filters were used: b30f (smooth), b36f (medium smooth), b46f (medium) and b60f (sharp). Based on our clinical experience in routine evaluation, b30f was used instead of b20f because of a better delineation of the coronary arteries and atherosclerotic lesions.

The scan geometry was based on a retrospective electrocardiogram (ECG)-gated protocol (the same used for the in vivo examination). This protocol is based on a low pitch that allows a retrospective reconstruction of multiple phases within the cardiac cycle. A demo-ECG was switched on and a heart rate of 71 occurred during the scan. The half-rotation reconstruction algorithm brought the effective temporal resolution down to 187 ms.

Data collection and analysis

A radiologist with a level 3 expertise in cardiac CT [21] loaded the data sets into a dedicated workstation (Leonardo; Siemens, Germany) and performed orthogonal views of all specimens for each convolution filter with two solutions.

One operator performed all of the measurements. Each specimen was evaluated for the presence of atherosclerotic plaque. The coronary wall thickening of plaque (minimum size observed: 0.3 mm2) was clearly visible from the surrounding oil with low attenuation and from the lumen after injection of contrast material. Plaques were targeted regardless of their size. The operator loaded the two data sets for each solution of the same specimen into a workstation screen, scrolling the data sets in parallel with standard soft-tissue window settings (window width:700 HU; window center:140 HU).

Once a plaque was detected in the two solutions by means of a comparative evaluation, the operator drew four regions of interest (ROIs) for the contrast material: in the lumen of the vessel (defined as lumen), the soft tissue of the coronary plaque (defined as noncalcified plaque), the calcification within the coronary artery wall (defined as calcified plaque) and the oil surrounding the plaque (defined as surrounding). The ROIs (minimum size: 0.1 mm2) were drawn as large as possible while avoiding the borders of each structure in order to limit the partial volume effect. Once the four ROIs were drawn, the operator could copy and paste these into other stacks of images. The two data sets for each solution were then evaluated using four convolution filters (b30f, b36f, b46f and b60f). For each ROI the mean and the standard deviation of the attenuation value were collected. Overall, 3200 attenuation values were obtained.

Statistical analysis

The attenuation values are presented as means and standard deviations. The signal-to-noise ratio (SNR) was calculated as the mean attenuation value of the ROI/standard deviation of the surrounding oil attenuation value. The standard deviation HU value of the surrounding air is often used for the estimation of noise. However, we did not include the air in our images. The standard deviation of the surrounding oil can estimate the noise because of its reliable signal. SNRs are presented as means for each kernel.

The contrast-to-noise ratio (CNR) was calculated as (mean attenuation A – mean attenuation B)/image noise (defined as the standard deviation of the surrounding oil). CNRs of (lumen – noncalcified plaque), (noncalcified plaque – surrounding), (calcified plaque – noncalcified plaque) and (calcified plaque – lumen) were calculated and presented as average for each kernel and solution.

Statistical evaluation was performed with dedicated software (SPSS ver. 10.1; SPSS, Chicago, Ill.). A one-way ANOVA test was performed to analyze the statistical difference between kernels in each structure for two solutions. A p < 0.0001 was considered for statistical significance. The Tukey HSD test evaluated the statistical absolute difference between any two groups of values. The Pearson’s correlation was also tested between all the different convolution filters.

After clustering for kernel and structure, the attenuation values obtained with the two solutions were compared with Student’s t-test and correlated with Pearson’s test.

Results

A total of 100 sections containing mixed plaques (n = 20) were evaluated in the seven coronary specimens. Two solutions and four convolution kernels were available (800 slices) at each level, and four ROIs (lumen, surrounding, calcified plaque, noncalcified plaque) were sampled in each slice. A total of 3200 measurements were performed. The results are summarized in Table 1. The attenuation values were all significantly different between the convolution filters in each structure with both solutions (p < 0.0001). The results of the Tukey HSD to test the difference between pairs of four filters showed that the noncalcified plaque attenuation values were all significantly different using both solutions (p < 0.01), with the exception of pair b36f–b46f. The attenuation values between b30f and b46f were significantly different at a high concentration in calcified plaques, whereas for the same pair of filters the attenuation values of the lumen and the surrounding were not significantly different. After clustering the paired attenuation values of the four structures for the convolution filter, a moderate-high correlation was generally found, with the exception of pairs b30f–b46f in surrounding (r = 0.150), b30f–b60f in lumen (r = 0.206) and surrounding (r = 0.311), b36f–b60f in lumen (r = 0.132) and b46f-b60f (r = 0.201) in lumen in the saline solution.

Table 1

Summary of the attenuation values measured in each structure in the two solutions (e.g. saline and contrast solutions) with increasing convolution filter

		b30f			b36f			b46f			b60f
		Mean	SD	SNR	Mean	SD	SNR	Mean	SD	SNR	Mean	SD	SNR
Lumen	Saline	58.1	40.4	29.9	52.1	40.5	9.4	48.2	36.0	19.1	23.5	6.0	1.8
	Contrast	329.1	93.7	209.8	312.6	85.5	44.0	334.1	91.7	169.7	331.3	138.3	17.8
Noncalcified Plaque	Saline	20.8	39.1	8.0	14.2	35.8	2.4	14.0	32.0	4.0	3.2	32.4	0.3
	Contrast	74.7	66.6	39.3	68.2	63.3	9.2	66.3	66.5	25.9	48.5	60.0	2.4
Calcified Plaque	Saline	740.5	392.0	327.0	758.5	360.3	133.6	785.7	388.5	294.9	1145.8	517.4	87.8
	Contrast	795.4	333.8	498.0	838.7	364.2	116.0	885.6	382.6	488.7	1194.8	520.7	61.0
Surrounding	Saline	127.1	2.7	−54.7	−128.0	3.1	−22.6	−128.7	2.7	−47.6	−134.8	11.3	−10.5
	Contrast	−128.2	2.6	−84.5	−126.8	4.3	−17.6	−128.5	2.4	−7.2	−123.7	10.0	−6.4

Whereas the convolution filter did not significantly affect the attenuation values of the lumen and the surrounding, sharper filters decreased the attenuation of the plaque and, conversely, increased the attenuation of the calcification with both solutions (Fig. 2).

Fig. 2

Mean attenuation of the four structures in the two solutions (e.g. saline and contrast) with increasing convolution filters. Whereas the convolution filter does not significantly affect the attenuation values of the lumen and the surrounding, sharper filters decrease the attenuation of the noncalcified plaque and, conversely, increase the attenuation of the calcified plaque with both solutions

The SNR values appear to have a common trend with increasing convolution filtering and both solutions. However, the best SNR was achieved with the b30f and b46f filters. Plaque values follow the lumen values, while the surrounding and calcification patterns are similar (Fig. 3).

Fig. 3

The SNR values appear to have a common trend with increasing convolution filtering and both solutions (e. g. saline and contrast). However, the best SNR is achieved with b30 and b46 filters. Noncalcified plaque values follow the lumen values, while the surrounding and the calcified plaque patterns are similar

The CNR values were calculated between structures with increasing convolution filters and different solutions (Table 2). Improved CNRs were achieved with the b30f and b46f filters. The lumen-noncalcified plaque CNR was particularly improved with the b30f and b46f filters using a high solution of contrast media. The noncalcified plaque-surrounding CNR was improved with the b30f and b46f filters in contrast solution. Calcified plaque-noncalcified plaque CNR increased with the contrast solution with the b30f and b46f filters, while it decreased using the b36f and b60f filters. The b46f filter was found to be ideal for calcified plaque-lumen contrast since it improved the CNR at contrast solution.

Table 2

Summary of CNR values calculated between structures with increasing convolution filters and different solutions

		b30f	b36f	b46f	b60f
Saline	Lumen – noncalcified plaque	21.8	7.1	15.1	1.5
	Non – calcified plaque-surrounding	62.7	25.0	51.6	10.8
	Calcified plaque – noncalcified plaque	318.9	131.3	290.8	87.4
	Calcified plaque – lumen	297.1	124.2	275.7	85.9
Contrast solution	Lumen – noncalcified plaque	170.6	34.8	143.8	15.4
	Noncalcified plaque – surrounding	123.7	26.9	97.1	8.8
	Calcified plaque – noncalcified plaque	458.7	106.7	462.8	58.5
	Calcified plaque – lumen	288.2	71.9	319.0	43.1

After clustering the paired attenuation values obtained for the four structures in order to test the statistical difference between saline and contrast solution for each convolution filter, values were all significantly different (p < 0.05) apart from surrounding with b46f and calcification with b60f. Pearson’s test has shown good correlation only between the attenuation values of calcification (r = 0.738, r = 0.826, r= 0.759, r = 0.723 with b30f, b36f, b46f, b60f, respectively).

Discussion

The evaluation of the coronary artery wall is a data set that has become available since the introduction of coronary MSCT. Several studies have reported the ability of MSCT to detect coronary plaques and subsequently characterize their tissue composition (i.e. calcified, mixed, noncalcified) [14, 20]. In an in vivo study, Kopp et al. showed an excellent correspondence between ICUS criteria (soft, intermediate and calcified plaques) and MSCT attenuation values [14]. Schroeder et al. subsequently found a strong correlation between tissue density measurements within the plaque and the qualitative ICUS criteria. The CT attenuation results were 14 ± 26 HU for plaques classified by ICUS as soft (i.e. predominantly lipid) plaques, 91 ± 21 HU for plaques classified by ICUS as intermediate (i.e. predominantly fibrous) plaques and 419 ± 194 HU for plaques classified by ICUS as calcified. These researchers reported no overlap in the mean attenuation among the three groups [15]. Leber et al. showed that lesion echogenity correlated well with MSCT attenuation measurements in coronary plaque, with MSCT correctly classifying 78% of sections containing hypoechoic plaque areas (soft plaques), 78% of sections containing hyperechoic plaque areas (fibrous plaques) and 95% of sections containing calcified plaque tissue. The MSCT density measurements within coronary lesions revealed significantly different values for hypoechoic, hyperechoic, and calcified plaques [17]. The close correlation between ICUS and MSCT in terms of the detection of calcified and noncalcified coronary plaques was also shown by Achenbach et al., although this group found a sensitivity of only 53% for exclusively noncalcified plaque [16].

Becker et al. compared the atherosclerotic lesions of 11 human cadaver heart specimens detected by MSCT with the histopathological macroscopic findings according to American Heart Association criteria. They concluded that MSCT is a promising tool for the characterization of atherosclerotic coronary plaques [22].

Such information is of great value since the risk of acute coronary syndromes caused by plaque disruption and thrombosis depends on plaque composition (i.e. noncalcified, predominantly lipid plaques with positive remodeling) rather than stenosis severity [18-20]. Reliable and reproducible non-invasive methods for the assessment of plaque constitution could be important in risk stratification of patients with CAD.

Nevertheless, several aspects (i.e. the impact of lumen attenuation, convolution filters, body mass index of the patient, CNR of the images and coronary calcification) of the methodology need to be better addressed in order to validate the accuracy of attenuation values measured within the coronary plaques using MSCT.

It is controversial if the absolute HU values are valid measurements for determining the plaque composition and identifying the corresponding tissue. A phantom study revealed that intracoronary attenuation can significantly affect the measured plaque attenuation so that accurate plaque tissue differentiation may be influenced [23]. Furthermore, Cademartiri et al. recently reported that the intravascular attenuation modifies significantly the attenuation of the coronary atherosclerotic plaques in an ex vivo coronary specimen [24]. Therefore, the characterization of the plaque on the basis of absolute attenuation values should be reported with caution, and the intraluminal attenuation should also be reported. The calibration of tissue attenuation values with contrast in-lumen values may improve accurate tissue component identification.

New possibilities for coronary artery plaque imaging may be offered by dual-source CT, which provides robust diagnostic image quality independent of the heart rate. Potential applications of dual-energy CT include tissue characterization and Ca quantification of atherosclerotic plaques [25, 26].

In our experimental study we addressed the issue of the variability of coronary plaque attenuation by varying the convolution filters (Fig. 4).

Fig. 4

Example of plaque with increasing filtering using a saline and contrast solution

The direct back projection of CT attenuation profiles results in unsharp images. To prevent unsharp and approximate images, each projection has to be convolved with a predetermined mathematical function, the convolution kernel; this kernel represents a filtering procedure that is able to enhance or reduce object boundaries. Convolution influences image characteristics by the choice and design of the kernel (i.e. smooth, sharp). A smooth filter reduces spatial resolution as well as image noise; a sharp filter increases spatial resolution as well as image noise. The noise is caused by the fluctuations in the number of X-ray quanta registered by the detector and is affected by convolution and dose. The SNR evaluates the ratio between the intensity of the signal and its statistical fluctuation. Suzuki et al. have demonstrated that the accuracy of a vascular diameter measurement with an automated software is affected by the choice of convolution filtering [27]. Seifarth et al. recently compared the effect of different convolution filtering on the visualization of coronary artery stents. Using the edge-enhancing b46f, they observed a significant improvement in lumen assessment compared to the conventional b30f [28].

Our study shows how the use of different convolution filters significantly affects the measured coronary plaque attenuation (Fig. 3). The use of sharper convolution filters produces an increase in the mean attenuation value of the calcified component of plaque. Increasingly sharper filtering produces a reduction in the mean attenuation values of the noncalcified plaque component. On the basis of our results, it is still challenging to rely on absolute HU values since the protocol used in coronary MSCT has not been standardized. Furthermore, absolute plaque attenuation values determined in the various studies are not comparable. It is reasonable to use a proper convolution filtering with improved SNR and CNR to better focus on plaque components. The b36f intermediate filter obtained a slight decrease of noncalcified plaque attenuation and an increase of calcified plaque attenuation. However, improved SNRs and CNRs were also achieved by the b30f and b46f filters. The b36f underestimated the noncalcified plaque attenuation and overestimated the calcified plaque attenuation in the likely attempt to weigh both components.

There are several limitations to our study. The first one is inherent to the lack of an evaluation of intra- and inter-observer variability. A second limitation is related to the fact that the coronary specimens were collected from elderly patients with a high prevalence of calcification. The blooming artifacts of calcification may impair coronary noncalcified plaque visualization and characterization based on HU values [29]. Furthermore, we did not investigate the influence of mAs and voltage (kV) on the blooming effect. Nevertheless, it has been reported that the voltage setting employed (i.e.120 kV) might be reasonably chosen during in vivo studies [29]. In vivo studies should also better address the issue of body mass, which affects coronary plaque measurements [30]. Moreover, according to data mined from the literature, it is reasonable to expect that plaque size (especially in calcified plaques) and lumen area would be influenced by the use of different convolution filters due to variable partial volume effect. Edge-enhancing sharper reconstruction filters are reported to increase the visible in-stent lumen diameter [28, 31].

The pathological correlation has not been provided since the aim of our study was not to compare pathology with MSCT, but to assess the influence of convolution filtering on plaque attenuation measurement.

In conclusion, the use of different convolution filtering significantly modified plaque attenuation values. Therefore, convolution filtering should be reported when attenuation values are shown in order to provide a standardization of the methodology. Sharper convolution kernels increased the attenuation of the calcium within the coronary plaques and reduced the attenuation of soft plaque tissues. Improved SNR and CNR were achieved by using the b30f and b46f filters. A smooth filter (e.g. b30f) may be used for clinical routine evaluation, while a medium filter (e.g. b46f) may be considered the best choice for the assessment of highly calcified or stented vessels. The use of a proper filtering according to plaque type could give a more reliable assessment of plaque attenuation values in terms of HU.