AJR Women's Imaging Online
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Matsuoka, S.
Right arrow Articles by Nakajima, Y.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Matsuoka, S.
Right arrow Articles by Nakajima, Y.
Social Bookmarking
 Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Reddit   Add to Technorati  
What's this?
DOI:10.2214/AJR.05.0047
AJR 2006; 186:1272-1279
© American Roentgen Ray Society


Original Research

Quantification of Thin-Section CT Lung Attenuation in Acute Pulmonary Embolism: Correlations with Arterial Blood Gas Levels and CT Angiography

Shin Matsuoka1, Yasuyuki Kurihara1, Kunihiro Yagihashi1, Hiroshi Niimi1 and Yasuo Nakajima1

1 All authors: Department of Radiology, St. Marianna University School of Medicine, 2-16-1 Sugao, Miyamae-Ku, Kawasaki, Kanagawa, 216-8511, Japan.

Received January 10, 2005; accepted after revision March 8, 2005.

 
Address correspondence to S. Matsuoka (shinma{at}d9.dion.ne.jp).


Abstract
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References
 
OBJECTIVE. The purposes of this study were to investigate the frequency histogram of lung attenuation on CT angiography (CTA) in patients with and without acute pulmonary embolism (PE) and to evaluate the relation of the frequency histogram of lung attenuation and hypoxemia.

MATERIALS AND METHODS. Twenty-six patients with PE and 11 patients without PE who underwent CTA were evaluated with frequency histograms. We obtained quantitative parameters such as mean lung attenuation, median lung attenuation, SD, skewness, kurtosis, and the proportion of lung attenuation except for the median ± 50 H (P ± 50 H). Lung attenuation was also assessed visually and scored. The relationship between those histogram parameters, or visual score, and PaO2 was evaluated. CTA scores for evaluation of the degree of pulmonary artery obstruction were obtained, and the relation with PaO2 was assessed.

RESULTS. No significant differences were found in mean lung attenuation and median lung attenuation between patients with and without PE. Meanwhile, SD, skewness, kurtosis, and P ± 50 H were significantly different between patients with and without PE (p = 0.0003, 0.0071, 0.0047, and 0.0028, respectively) and significantly correlated with PaO2 (r = -0.770, 0.797, 0.786, -0.871, respectively). Significant differences were found in visual scores between patients with and without PE (p < 0.0001). There were significant but relatively low correlations between CTA score and arterial blood gas levels (r = -0.442, p = 0.03).

CONCLUSION. In patients with acute PE, heterogeneity in lung attenuation is more prominent than in patients without PE.

Keywords: CT • CT angiography • embolism • lung


Introduction
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References
 
The importance of CT in the diagnosis of pulmonary embolism (PE) has grown in recent years, mainly because of the advent of fast CT image acquisition techniques. However, most previous investigations have focused on the ability of CT scans to reveal PE and the widely described features of arterial emboli. Only a few reports have addressed the use of CT angiography (CTA) to evaluate the severity of acute PE [1, 2].

In patients with acute PE, various pathophysiologic changes such as regional alterations of pulmonary blood flow [3, 4], bronchoconstriction [5-7], edema [8], and infarction [9, 10] are seen. Lung attenuation values on CT images are affected by reduced pulmonary perfusion [11-14], bronchoconstriction [15], and increased intravascular volume or edema [16, 17]. These pathophysiologic changes are also recognized in nonembolic regions [18] and result in ventilation-perfusion (V/Q) ratio inequality [3, 19-22]. Thus, theoretically, lung attenuation in patients with acute PE may become heterogeneous. However, several investigators have reported that heterogeneous lung attenuation, termed "mosaic perfusion," is a less frequent and uncommon feature of acute PE [23, 24]. Unfortunately, in those previous studies, mosaic perfusion was assessed with visual and subjective evaluations. To our knowledge, quantitative analysis of lung attenuation in patients with acute PE has not been assessed.

Frequency histograms of CT scans of the lung provide quantitative indexes—such as mean lung attenuation, skewness, and kurtosis—for the evaluation of the distribution of lung attenuation values [25, 26]. Skewness is a measure of the degree of asymmetry of a distribution. Kurtosis is a measure of whether the data are peaked or flat relative to a normal distribution. This quantitative analysis has been used for the evaluation of interstitial pulmonary fibrosis, and the quantitative indexes have been shown to correlate with pulmonary function [26]. Considering the pathophysiologic characteristics of acute PE, we hypothesized that the frequency histograms of lung attenuation in patients with PE differ from those in patients without PE (referred to hereafter as "non-PE patients"). In addition, many patients with acute PE develop clinically significant hypoxemia. Although the exact mechanism by which this hypoxemia occurs is controversial, the potential mechanism of hypoxemia in acute PE can primarily be explained by V/Q ratio inequality [3, 19-22], which causes heterogeneous regional alterations of ventilation and perfusion. Therefore, we also hypothesized that there may be a relation between distribution of lung attenuation values and the degree of hypoxemia in patients with acute PE.

The purpose of this study was to investigate the distribution of lung attenuation values on CTA with quantitative histogram analysis in patients with acute PE in comparison with non-PE patients. In addition, we assessed the relation between the degree of hypoxemia and the distribution of lung attenuation values compared with the degree of arterial obstruction scored using CTA.


Materials and Methods
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References
 
Subjects
Our institutional review board required neither its approval nor patient informed consent for this retrospective study. One radiologist retrospectively reviewed 186 consecutive patients who had undergone CTA on a 4-MDCT scanner for clinically suspected acute PE between July 2001 and August 2004. Of those patients, 56 patients had an endoluminal filling defect on CTA. Of those positive studies, we selected 48 patients who underwent arterial blood gas (ABG) measurements while breathing room air within 24 hr of CTA. Twenty-two patients were excluded because of prior cardiopulmonary disease (n = 3), obvious abnormal parenchymal lesions (consolidation, atelectasis, or emphysema) (n = 11), extensive pleural effusions (n = 3), nondiagnostic CT image due to poor opacification of contrast material or major motion artifacts (n = 2), and CTA performed 48 hr after clinical presentation (n = 3). Thus, 26 patients (seven men and 19 women; mean age, 63 years; age range, 32-83 years) were included in the PE group. No patient had a history of smoking in the PE group. None of the patients in the PE group received therapeutic anticoagulation before CTA and ABG measurements.

In the 130 patients who had no sign of PE on CTA, 58 patients underwent CT venography of the lower extremities at the same time. Of those patients, we selected 15 patients who also had a normal D-dimer level. We excluded three patients who had prior cardiopulmonary disease (n = 1) or obvious abnormal parenchymal lesions (n = 2). One patient had a history of smoking, so that patient was also excluded. Thus, 11 patients (five men and six women; mean age, 56 years; age range, 39-87 years) were included in the non-PE group.

CT Scanning
CTA was performed on one of two MDCT scanners: Aquilion (Toshiba Medical Systems), with a collimation of 4 x 1 mm, a table speed of 8 mm per rotation, and a 0.5-sec rotation time (PE group, n = 16; non-PE group, n = 7), or Asteion (Toshiba Medical Systems), with a collimation of 4 x 2 mm, a table speed of 16 mm per rotation, and a 0.75-sec rotation time (PE group, n = 10; non-PE group, n = 5), 120 kV, and 200-250 mA using a 512 x 512 matrix. CTA was performed from 1 cm above the top of the aortic arch to the diaphragm at deep inspiration during a single breath-hold. A bolus of 100 mL of iodinated contrast material was injected at a rate of 3 mL/sec with an automatic injector. The bolus-tracking method was used for optimizing pulmonary artery opacification. A region of interest was placed on the main pulmonary artery, and scanning was initiated at a preset level of opacification. With this protocol, we consistently achieved high and uniform contrast enhancement throughout the thorax in all patients. For each patient, the raw data were retrospectively reconstructed with a section thicknesses of 2 mm.

Evaluation of CT Findings
Quantitative assessment of lung attenuation—CT images were transferred to a PC for quantitative analysis of lung attenuation. The CT images with 2-mm section thickness and 10-mm spacing were selected from among all of the CT images of each patient, and these CT images were segmented using a semiautomatic image-processing program (ImageJ [version 1.31], a public domain Java image-processing program inspired by the National Institutes of Health program Image for the Macintosh [Apple Computer]). This program uses a semiautomatic threshold technique to isolate the lungs from other tissues and structures and selects all pixels between -500 and -1,024 H. Minimal user intervention by one radiologist was required to exclude nonlung structures that satisfy the threshold criteria, such as the trachea, blood vessels, and large bronchi near the hilum. CT slices with a significant interlobar fissure area or extensive artifacts caused by the inflow of high concentrations of contrast medium in the superior vena cava were removed from the analysis.

We used a 16-bit histogram plug-in, the Java programming language, which is available on the ImageJ Web site [27], to obtain frequency histograms. Gray-scale frequency histograms for the combined lung parenchyma of each subject were plotted. Because lung size varied from one patient to another, the number of pixels was preferentially expressed as a percentage of the total number of pixels. Values for the mean lung attenuation, median lung attenuation, SD, skewness, and kurtosis of the gray-scale frequency histogram for each subject were calculated using statistical software (JMP software [version 5.0.1], SAS Institute). Skewness describes the degree of asymmetry of a histogram; a histogram with a long tail to the right has a positive skewness value, and a completely symmetric distribution has a skewness value of zero. Kurtosis describes how sharply peaked a histogram is; a histogram that is more peaked than a normal distribution has a positive kurtosis value, and a normal distribution has a kurtosis of zero. Skewness is given by the following equation:

Formula
and kurtosis is given by the following equation:

Formula
where ni is the number of pixels at lung attenuation value xi, x is the mean value of lung attenuation value, s is the SD, and N is the total number of pixels. In addition, the proportion of pixels that existed except for ± 50 H of the median value of segmented lung attenuation was calculated as P ± 50 H (%).

Visual assessment of lung attenuation—The CT images obtained with 2-mm section thickness and 10-mm spacing were selected from among all of the CT images of each patient. The same CT images were viewed as those for quantitative assessment. These CT images were evaluated by two experienced radiologists independently without knowledge of clinical data, and final decisions about the findings were reached by consensus. Regional attenuation differences between normal and abnormal lung parenchyma are often small [28]; thus, window settings for the detection of subtle lung attenuation differences (width, 250 H; center, median lung density) were applied. The patterns of lung attenuation on single-detector images were divided into two categories: 1, homogeneously normal or gravity-dependent gradient lung attenuation with no obvious difference in attenuation values or patterns between the right and left lungs; and 2, geographic regions of heterogeneous lung attenuation or attenuation values or patterns that differ between the right and left lung parenchyma. Category 1 was assigned 0 points, and category 2 was assigned 1 point. In each patient, the visual score was calculated by adding scores for each slice and then dividing by slice number.

CTA—CTA was evaluated on a monitor in a scroll-through or cine mode with the standard window settings for viewing CTA (width, 500 H; center, 80 H) at our institution. Two experienced radiologists without knowledge of the clinical data independently evaluated the images, and final decisions about the findings were reached by consensus. PE was diagnosed when either complete or partial filling defects within the main, lobar, or segmental arteries were identified. CTA scores described in the literature [1] were applied to evaluate the degree of pulmonary artery obstruction.

To define the CT obstruction index, we regarded the arterial tree of the lung as having 18 segmental arteries: three to the right upper lobes, two to the left upper lobes, two to the middle lobe and to the lingula, five to the right lower lobes, and four to the left lower lobes. The presence of an embolus in a segmental artery was scored as 1 point, and emboli in the most proximal arterial level were scored a value equal to the number of segmental arteries arising distally. To provide additional information about the residual perfusion distal to the embolus, a weighting factor was assigned to each value, depending on the degree of vascular obstruction. This factor was equal to 0 when no thrombus was observed; 1, when partially occlusive thrombus was observed; or 2, with total occlusion. Thus, the maximal CT obstruction index was 36 per patient. Isolated subsegmental embolus was considered as a partially occluded segmental artery and was assigned a value of 1. The percentage of vascular obstruction was calculated by dividing the patient score by the maximal total score and by multiplying the result by 100. Therefore, the CT obstruction index can be expressed as follows:

Formula
where n is the value of the proximal thrombus in the pulmonary arterial tree equal to the number of segmental branches arising distally (minimum, 1; maximum, 18), and d is the degree of obstruction (minimum, 0; maximum, 2).

ABG Measurements
In the present study, we evaluated patients with PE who had ABG levels measured while breathing room air. All measurements of ABG were obtained with a blood gas analyzer (ABL520, Radiometer) during the 24 hr before or after CTA. The following physiologic variable was measured: partial pressure of O2 (PaO2).

Statistical Analysis
Comparisons of the mean for mean lung attenuation, median lung attenuation, SD, skewness, kurtosis, and P ± 50 H between the PE group and the non-PE group were performed using the Mann-Whitney U test. Comparison of the visual score of lung attenuation between the PE group and the non-PE group was also performed using the Mann-Whitney U test. The linear regression analysis and Spearman's rank correlation analysis were used to evaluate the relationship between, first, PaO2 and the mean lung attenuation, median lung attenuation, SD, skewness, kurtosis, or P ± 50 H; second, PaO2 and the visual lung attenuation score; third, PaO2 and the CTA score; and, fourth, the CTA score and the mean lung attenuation, median lung attenuation, SD, skewness, kurtosis, or P ± 50 H. All statistical analyses were performed using JMP 5.0.1 software. Data are expressed as mean ± SD. For all statistical analyses, a p value of less than 0.05 was considered significant.


Results
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References
 
Comparisons of the Measured Parameters of Histograms Between PE and Non-PE Groups
A total of 416 slices were obtained from the 26 patients with acute PE for quantitative analyses. Thirty-seven slices that had a significant low-attenuation area due to interlobar fissure or extensive artifacts caused by the inflow of high concentrations of contrast medium in the superior vena cava were excluded, and 379 slices were included in the analysis of the frequency histogram of lung attenuation. A total of 180 slices were obtained from the 11 non-PE patients, and 21 slices that had significant extent of interlobar fissure or extensive artifacts caused by contrast medium in the superior vena cava were excluded, leaving 159 slices for quantitative analysis.

The results of quantitative analysis of lung attenuation are shown in Table 1. The mean values of the mean lung attenuation and median lung attenuation for the PE group and non-PE group were not significantly different (p = 0.273 and p = 0.158, respectively). Other measured parameters—SD, skewness, kurtosis, and P ± 50 H—of the frequency histogram were significantly different between patients with PE and those without PE. The mean value of the SD of the PE group (87.1 ± 8.4) was significantly higher than that of non-PE group (78.9 ± 2.9) (p = 0.0003). The mean value of skewness of the PE group (1.36 ± 0.5) was significantly lower than that of the non-PE group (1.82 ± 0.2) (p = 0.0071). The mean value of kurtosis of the PE group (2.04 ± 1.7) was significantly lower than that of the non-PE group (3.78 ± 1.1) (p = 0.0047). The mean value of P ± 50 H in the PE group (42.2% ± 10.2%) was significantly greater than that of the non-PE group (31.7% ± 4.0%) (p = 0.0028). Typical frequency histograms are shown in a patient with PE (Figs. 1A and 1B) and one without PE (Figs. 2A and 2B). All PE and non-PE patients had a unimodal distribution histogram. Non-PE patients had the histogram with a skewed long tail to the right, toward increasing gray-scale density. The patients with PE, however, had a distribution that was more symmetric and flatter (lower kurtosis) than that of the non-PE patients.


View this table:
[in this window]
[in a new window]

 
TABLE 1: Quantitative Analysis of Lung Attenuation in Patients With Pulmonary Embolism (PE) and Those Without PE (Non-PE)

 

Figure 1
View larger version (14K):
[in this window]
[in a new window]
[as a PowerPoint slide]
 
Fig. 1A —48-year-old woman with acute pulmonary embolism (PE). Graph shows frequency histogram for lung parenchyma. Note reduced skewness (1.4) and kurtosis (2.3) compared with those in non-PE patient. SD and P ± 50 H are 88.3% and 43.3%, respectively.

 

Figure 2
View larger version (86K):
[in this window]
[in a new window]
[as a PowerPoint slide]
 
Fig. 1B —48-year-old woman with acute pulmonary embolism (PE). Transverse thin-section (2-mm collimation) CT scan (window level, -873 H; window width, 250 H) shows inhomogeneous lung attenuation.

 

Figure 3
View larger version (18K):
[in this window]
[in a new window]
[as a PowerPoint slide]
 
Fig. 2A —61-year-old woman without acute pulmonary embolism. Graph shows frequency histogram for lung parenchyma. SD, skewness, kurtosis, and P ± 50 H are 77.1, 2.3, 6.2, and 28.5%, respectively.

 

Figure 4
View larger version (87K):
[in this window]
[in a new window]
[as a PowerPoint slide]
 
Fig. 2B —61-year-old woman without acute pulmonary embolism. Transverse thin-section (2-mm collimation) CT scan (window level, -903 H; window width, 250 H) shows relatively homogeneous lung attenuation.

 

The Relations Between Measured Histogram Parameters and ABG Measurements
The correlation between the measured parameters of the frequency histogram and measurements of PaO2 is shown in Table 2. The mean PaO2 was 67.7 ± 10.8 mm Hg. No significant correlation was found between the mean lung attenuation and PaO2 (r = 0.280, p = 0.39). Similarly, there was no significant correlation between the median lung attenuation and PaO2 (r = 0.309, p = 0.23). In contrast, significant correlations were found between the other measured parameters (SD, skewness, kurtosis, and P ± 50 H) of the frequency histogram and ABG levels (Figs. 3A, 3B, 3C, and 3D). A significant inverse relationship between the SD and PaO2 (r = -0.770, p = 0.001) was found. Significant positive relationships were found between PaO2 and both skewness (r = 0.797, p < 0.0001) and kurtosis (r = 0.786, p < 0.0001). The values of the P ± 50 H were highly correlated with PaO2. The values of P ± 50 H were increased with decreasing PaO2 levels (r = -0.871, p < 0.0001).


View this table:
[in this window]
[in a new window]

 
TABLE 2: Correlation Between the Histogram Parameters and PaO2

 

Figure 5
View larger version (10K):
[in this window]
[in a new window]
[as a PowerPoint slide]
 
Fig. 3A —Correlation between PaO2 and SD, skewness, kurtosis, and P ± 50 H in patients with pulmonary embolism. Scatterplots show significant correlation between PaO2 and SD (A) (r = -0.770, p = 0.001), PaO2 and skewness (B) (r = 0.797, p < 0.0001), PaO2 and kurtosis (C) (r = 0.786, p < 0.0001), and PaO2 and P ± 50 H (D) (r = -0.871, p < 0.0001).

 

Figure 6
View larger version (10K):
[in this window]
[in a new window]
[as a PowerPoint slide]
 
Fig. 3B —Correlation between PaO2 and SD, skewness, kurtosis, and P ± 50 H in patients with pulmonary embolism. Scatterplots show significant correlation between PaO2 and SD (A) (r = -0.770, p = 0.001), PaO2 and skewness (B) (r = 0.797, p < 0.0001), PaO2 and kurtosis (C) (r = 0.786, p < 0.0001), and PaO2 and P ± 50 H (D) (r = -0.871, p < 0.0001).

 

Figure 7
View larger version (10K):
[in this window]
[in a new window]
[as a PowerPoint slide]
 
Fig. 3C —Correlation between PaO2 and SD, skewness, kurtosis, and P ± 50 H in patients with pulmonary embolism. Scatterplots show significant correlation between PaO2 and SD (A) (r = -0.770, p = 0.001), PaO2 and skewness (B) (r = 0.797, p < 0.0001), PaO2 and kurtosis (C) (r = 0.786, p < 0.0001), and PaO2 and P ± 50 H (D) (r = -0.871, p < 0.0001).

 

Figure 8
View larger version (11K):
[in this window]
[in a new window]
[as a PowerPoint slide]
 
Fig. 3D —Correlation between PaO2 and SD, skewness, kurtosis, and P ± 50 H in patients with pulmonary embolism. Scatterplots show significant correlation between PaO2 and SD (A) (r = -0.770, p = 0.001), PaO2 and skewness (B) (r = 0.797, p < 0.0001), PaO2 and kurtosis (C) (r = 0.786, p < 0.0001), and PaO2 and P ± 50 H (D) (r = -0.871, p < 0.0001).

 

Visual Scores of Lung Parenchyma
The same CT slices that were evaluated for the analysis of distribution histogram of lung attenuation were used for the visual score assessment. The mean visual scores for lung attenuation were 0.60 ± 0.16 (range, 0.28-0.88) in the PE group and 0.22 ± 0.16 (range, 0-0.60) in the non-PE group. There was a significant difference between the PE and non-PE groups (p < 0.0001). Significant correlations were found between visual score and PaO2 (r = -0.709, p = 0.004) (Fig. 4).


Figure 9
View larger version (11K):
[in this window]
[in a new window]
[as a PowerPoint slide]
 
Fig. 4 —Correlation between visual scores and PaO2. Significant correlation was found between visual scores and PaO2 (r = -0.709, p = 0.0004).

 
CTA Score
The mean CTA score was 35.4% ± 19.4% (range, 5.6-69.4%). There was significant inverse correlation between the CTA score and PaO2 (r = -0.442, p = 0.03) (Fig. 5). No significant correlation was found between the CTA score and mean lung attenuation, median lung attenuation, SD, skewness, kurtosis, or P ± 50 H (r = 0.01, 0.05, 0.27, -0.226, -0.25, and -0.26, respectively).


Figure 10
View larger version (11K):
[in this window]
[in a new window]
[as a PowerPoint slide]
 
Fig. 5 —Correlation between CT angiography (CTA) score and PaO2. Mean CTA score was 35.4% ± 19.4% (range, 5.6-69.4%). There was significant correlation between CTA score and PaO2 (r = -0.442, p = 0.03).

 

Discussion
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References
 
In the present study, we found that the frequency histogram of the lung attenuation on CTA was different between patients with and without acute PE. There was heterogeneous lung attenuation in patients with acute PE. In patients with acute PE with less skewness and less kurtosis, distribution histograms of lung attenuation were more symmetric and flatter than those of patients without PE. In addition, we found significant correlations between the parameters of the frequency histogram and hypoxia. The variation of the lung attenuation increased with the deterioration of PaO2. Although we cannot explain the correlation seen in this study, some speculation is possible.

Recently, visualization and quantification of pulmonary perfusion based on various CT techniques have been introduced [11-14]. Contrast-enhanced dynamic helical CT densitometry allows the assessment of at least segmental perfusion defects in patients with acute PE [11-14]. The segments with reduced perfusion showed a significantly lower mean CT attenuation value on the enhanced scans than the segments with normal perfusion. Im et al. [15] reported that they did not find lung attenuation changes in pigs during the 28 days immediately after experimentally induced PE. However, they evaluated lung attenuation with unenhanced CT. IV contrast administration produced a rise in the overall mean lung attenuation [29]. Thus, attenuation differences between PE and normal areas could not be detected with unenhanced CT. Although the change of the attenuation value by the perfusion defect is subtle, it may be related to some extent to the heterogeneity of the lung attenuation in patients with PE.

In patients with acute PE, generalized or local bronchoconstriction has been invoked [5-7]. The proposed mechanism includes bronchoactive amines, such as serotonin and prostaglandins, that are released during platelet aggregation in the thrombus, as has been shown in canine lungs [30], or a change in parasympathetic nervous system tension, which controls bronchial smooth-muscle tension [18]. In an experimental study, bronchial obstruction resulted in lobular areas of low attenuation on high-resolution CT [15]. Arakawa et al. [31] studied 15 patients with PE and confirmed that clots were identified on CTA in 90.6% of the arteries proximal to the areas of decreased attenuation. They also reported that 71.9% of the areas with decreased attenuation on inspiratory scans showed air-trapping on expiratory scans. Therefore, decreased blood flow and bronchoconstriction may lead to the existence of low-attenuation areas in patients with acute PE.

Lung parenchyma distal to a pulmonary thromboembolus is pathologically normal or shows only mild atelectasis and minimal intraalveolar hemorrhage or edema. Especially in the early stage, histologic examination shows only intraalveolar hemorrhage and edema, with intact alveolar wall. Grossman et al. [32] found increased lung attenuation after arterial occlusion in dogs. In addition, pathophysiologically the blood flow increase may be generated in non-PE regions [4], and increased intravascular volume increases lung attenuation [16]. These pathologic and physiologic changes might lead, in part, to increases in lung attenuation.

In this study, the heterogeneity of lung attenuation in patients with acute PE was confirmed even by visual analysis. However, several investigators have reported that mosaic perfusion or the presence of low-attenuation area is not a common feature of acute PE [23, 24, 33]. Coche et al. [23] compared lung parenchymal findings of 88 patients with suspected acute PE who underwent CTA. Those researchers found mosaic perfusion in 12% of patients with acute PE and in 10% of patients without this finding. Shah et al. [24] evaluated parenchymal findings in 28 patients with acute PE and found only 7% with mosaic perfusion. However, Arakawa et al. [31] found mosaic perfusion in four (44.4%) of nine patients with acute PE.

The difference of frequency in these studies may result from variability in patient selection; these studies may differ from our study because we used a special setting for window level and window width. Groell et al. [11] reported that the segment with reduced perfusion showed a lower attenuation value on enhanced CT than segments with normal perfusion. However, their attempts to detect parenchymal perfusion defects by visual assessment have been unsuccessful because of subtle differences between defects and normal perfusion. We analyzed lung parenchyma with enhanced contrast imaging using a narrow window setting. Thus, a subtle attenuation difference could be detected in our study. The existence of mosaic perfusion in patients with acute PE may be missed with the usual window setting for lung parenchyma. However, there are various problems in the visual assessment of lung attenuation. We used the median value of lung attenuation as the center level, but we cannot assert that this method was the most appropriate one. Future investigation, including the reproducibility of our method, is required.

In our PE group, significant correlations were found between PaO2 and SD, skewness, kurtosis, or P ± 50 H. Therefore, the degree of heterogeneity in lung attenuation may relate to hypoxia. Hypoxemia is commonly reported in patients with acute PE. The exact mechanisms by which hypoxemia occurs remain in dispute. In general, hypoxemia is mainly explained by V/Q ratio inequality [3, 19-22]. The high proportion of the V/Q ratio is explained by the preferential reduction in perfusion within the region of the emboli. The low V/Q ratio is mainly due to a shifting of blood flow away from the region of the emboli toward regions with normal ventilated units [3, 19-21]. The heterogeneity of lung attenuation on CTA might reflect these pathophysiologic changes in acute PE.

Quantitative CT parameters, such as mean lung attenuation, skewness, and kurtosis, can be obtained from frequency histograms of the lung [25, 26]. This method has been applied to the quantitative evaluation for interstitial pulmonary fibrosis, and those parameters correlate with results of pulmonary function tests [26]. In our PE patients, these parameters of frequency histograms of lung attenuation differed from those of non-PE patients. The distribution of the lung attenuation values in the PE patients was characterized by a simultaneous higher proportion of lung attenuation at both extremities of the histogram, with subsequent widening of the SD and flattening of the middle peak. Quantitative analysis of lung scintigraphy scans based on the frequency histograms of the V/Q ratio was applied in patients with PE, and the histograms of the V/Q distribution showed widening of the second moment of the distribution and flattening of the middle peak, comparable to the results of our histogram analysis [34].

With regard to the distribution of the V/Q ratio using multiple perfusions of inert gases in patients with acute PE, data have also shown an abnormal distribution of lung V/Q ratios to alveolar units, with low V/Q ratios, and have helped explain the abnormalities of blood gas exchange [3, 19]. Undoubtedly, lung attenuation values and the V/Q ratio are based on the different pathophysiologic mechanisms. Low-attenuation areas on enhanced CT in patients with acute PE may relate mainly to decreased blood flow and bronchoconstriction. However, histogram parameters in our PE patients correlated with the ABG levels and might relate to abnormalities in blood gas exchange.

In previous studies, CTA scores were used for the quantitative evaluation of acute PE severity, and those scores correlated with clinical data such as oxygen saturation and echocardiographic findings [1, 2]. With regard to the importance of the size of embolic clots, experimental studies have shown that embolus size significantly influences gas exchange characteristics [35]. We adopted the method of Qanadli et al. [1] considering the size of the embolic clots. However, in our PE patients, there was relatively a low correlation between the CTA score and ABG levels. The method of quantitative assessment of lung attenuation can be used to evaluate pulmonary blood flow without influencing the existence or size of the thrombus. Furthermore, there were no significant correlations between the CTA score and the SD, skewness, kurtosis, or P ± 50 H. These parameters may reflect the reserve capacity of pulmonary function in patients with PE. By combining CTA and quantitative analysis of lung attenuation, a comprehensive and noninvasive diagnosis of thoracic structure and function is feasible with a single technique in patients with acute PE.

There are various problems with quantitative analysis of lung attenuation in patients with acute PE. First, we could not clarify the relation of lung attenuation and ventilation or blood flow in this study. Low-attenuation areas might reflect hypoperfusion, bronchoconstriction, or both. In our PE patients, few patients underwent V/Q scintigraphy because the diagnosis of PE was based on CTA examinations in most cases. To clarify the relation of bronchoconstriction, the combined use of expiration CT will be considered in the future. Second, there was considerable variability in the lung density map even among healthy individuals [29]. In addition, an attenuation gradient lung is normally present, with the most dependent lung regions being densest. However, indexes of the quantitative attenuation analysis were significantly different between the PE and non-PE groups. Third, quantitative measurements of lung attenuation are known to be affected by several factors, such as variation in CT techniques [36]. These variations in our study might have influenced the results to some extent. Although several methods were applied to correct lung attenuation for the quantitative histogram analysis of lung disease [25, 26], to our knowledge, standardization of CT densitometry has not been achieved completely. Thus, we did not correct the gray-scale values measured in current study. Further work needs to be done to develop an appropriate correction factor for lung attenuation.

Several other limitations of this study exist. First, we excluded patients with severe PE such as those requiring mechanical ventilation or having pulmonary infarction. In the complete obstruction of a unilateral pulmonary artery, lung attenuation may lower uniformly. Second, we could not conclude whether the threshold value that evaluates the low-attenuation area was appropriate. Third, although all subjects were nonsmokers, we could not exclude the possibility of lung attenuation being affected in subjects with subclinical small airways disease.

In conclusion, the parameters of frequency histograms of lung attenuation on contrast-enhanced CT were significantly different between patients with and without PE. In patients with acute PE, heterogeneity in lung attenuation was more prominent, which may be due to alteration of perfusion and ventilation distribution. In addition, significant correlations were found between those parameters and hypoxia. The degree of heterogeneity in lung attenuation may relate to abnormal gas change. By combining CTA and quantitative analyses of lung attenuation, diagnosis and severity assessment of acute PE are feasible with a single technique. Moreover, quantitative CT analyses of lung attenuation may be helpful for pathophysiologic analysis of patients with acute PE.


References
Top
Abstract
Introduction
Materials and Methods
Results
Discussion
References
 

  1. Qanadli SD, El Hajjam M, Vieillard-Baron A, et al. New CT index to quantify arterial obstruction in pulmonary embolism: comparison with angiographic index and echocardiography. AJR2001; 176:1415 -1420[Abstract/Free Full Text]
  2. Bankier AA, Janata K, Fleischmann D, et al. Severity assessment of acute pulmonary embolism with spiral CT: evaluation of two modified angiographic scores and comparison with clinical data. J Thorac Imaging 1997; 12:150 -158[Medline]
  3. Dantzker DR, Wagner PD, Tornabene VW, Alazraki NP, West JB. Gas exchange after pulmonary thromboembolization in dogs. Circ Res 1978; 42:92 -103[Abstract/Free Full Text]
  4. Altemeier WA, Robertson HT, McKinney S, Glenny RW. Pulmonary embolization causes hypoxemia by redistributing regional blood flow without changing ventilation. J Appl Physiol1998; 85:2337 -2343[Abstract/Free Full Text]
  5. Gurewich V, Thomas D, Stein M, Wessler S. Bronchoconstriction in the presence of pulmonary embolism. Circulation1963; 27:339 -345[Abstract/Free Full Text]
  6. Robinson AE, Puckett CL, Green JD, Silver D. In vivo demonstration of small-airway bronchoconstriction following pulmonary embolism. Radiology 1973;109 : 283-286[Medline]
  7. Sasahara AA, Cannilla JE, Morse RL, Sidd JJ, Tremblay GM. Clinical and physiologic studies in pulmonary thromboembolism. Am J Cardiol 1967; 20:10 -20[CrossRef][Medline]
  8. Hyers TM, Fowler AA, Wicks AB. Focal pulmonary edema after massive pulmonary embolism. Am Rev Respir Dis1981; 123:232 -233[Medline]
  9. Horvath G, Noma S, Moskowitz GW, Khan A, Balazs G, Herman PG. Experimental pulmonary infarction in a pig model. Invest Radiol 1992; 27:829 -835[CrossRef][Medline]
  10. Dalen JE, Haffajee CI, Alpert JS 3rd, Howe JP, Ockene IS, Paraskos JA. Pulmonary embolism, pulmonary hemorrhage and pulmonary infarction. N Engl J Med 1977;296 : 1431-1435[Abstract]
  11. Groell R, Peichel KH, Uggowitzer MM, Schmid F, Hartwagner K. Computed tomography densitometry of the lung: a method to assess perfusion defects in acute pulmonary embolism. Eur J Radiol1999; 32:192 -196[CrossRef][Medline]
  12. Geraghty JJ, Stanford W, Landas SK, Galvin JR. Ultrafast computed tomography in experimental pulmonary embolism. Invest Radiol 1992; 27:60 -63[Medline]
  13. Screaton NJ, Coxson HO, Kalloger SE, et al. Detection of lung perfusion abnormalities using computed tomography in a porcine model of pulmonary embolism. J Thorac Imaging2003; 18:14 -20[CrossRef][Medline]
  14. Schoepf UJ, Bruening R, Konschitzky H, et al. Pulmonary embolism: comprehensive diagnosis by using electron-beam CT for detection of emboli and assessment of pulmonary blood flow. Radiology2000; 217:693 -700[Abstract/Free Full Text]
  15. Im JG, Choi YW, Kim HD, Jeong YK, Han MC. Thin-section CT findings of the lungs: experimentally induced bronchial and pulmonary artery obstruction in pigs. AJR 1996;167 : 631-636[Abstract/Free Full Text]
  16. Herold CJ, Wetzel RC, Robotham JL, Herold SM, Zerhouni EA. Acute effects of increased intravascular volume and hypoxia on the pulmonary circulation: assessment with high-resolution CT. Radiology 1992;183 : 655-662[Abstract/Free Full Text]
  17. Slutsky RA, Long S, Peck WW, Higgins CB, Mattrey R. Pulmonary density distribution in experimental noncardiac canine pulmonary edema evaluated by computed transmission tomography. Invest Radiol 1984; 19:168 -173[CrossRef][Medline]
  18. Austin JH, Sagel SS. Alterations of airway caliber after pulmonary embolization in the dog. Invest Radiol1972; 7:135 -139[CrossRef][Medline]
  19. Santolicandro A, Prediletto R, Fornai E, et al. Mechanisms of hypoxemia and hypocapnia in pulmonary embolism. Am J Respir Crit Care Med 1995; 152:336 -347[Abstract]
  20. Manier G, Castaing Y, Guenard H. Determinants of hypoxemia during the acute phase of pulmonary embolism in humans. Am Rev Respir Dis 1985; 132:332 -338[Medline]
  21. Huet Y, Lemaire F, Brun-Buisson C, et al. Hypoxemia in acute pulmonary embolism. Chest 1985;88 : 829-836[Abstract/Free Full Text]
  22. Levy SE, Simmons DH. Mechanism of arterial hypoxemia following pulmonary thromboembolism in dogs. J Appl Physiol1975; 39:41 -46[Abstract/Free Full Text]
  23. Coche EE, Müller NL, Kim KI, Wiggs BR, Mayo JR. Acute pulmonary embolism: ancillary findings at spiral CT. Radiology 1998;207 : 753-758[Abstract/Free Full Text]
  24. Shah AA, Davis SD, Gamsu G, Intriere L. Parenchymal and pleural findings in patients with and patients without acute pulmonary embolism detected at spiral CT. Radiology 1999;211 : 147-153[Abstract/Free Full Text]
  25. Hartley PG, Galvin JR, Hunninghake GW, et al. High-resolution CT-derived measures of lung density are valid indexes of interstitial lung disease. J Appl Physiol 1994;76 : 271-277[Abstract/Free Full Text]
  26. Best AC, Lynch AM, Bozic CM, Miller D, Grunwald GK, Lynch DA. Quantitative CT indexes in idiopathic pulmonary fibrosis: relationship with physiologic impairment. Radiology 2003;228 : 407-414[Abstract/Free Full Text]
  27. ImageJ Web site. Available at: http://rsb.info.nih.gov/ij/. Accessed February 13, 2006
  28. Yang GZ, Chabat F, Hansell DM. Enhancement of subtle density differences of the lung parenchyma on CT. Br J Radiol1998; 71:686 -690[Abstract]
  29. Adams H, Bernard MS, McConnochie K. An appraisal of CT pulmonary density mapping in normal subjects. Clin Radiol1991; 43:238 -242[CrossRef][Medline]
  30. Lindsey HE, Wyllie JH. Release of prostaglandins from embolized lungs. Br J Surg 1970;57 : 738-741[Medline]
  31. Arakawa H, Kurihara Y, Sasaka K, Nakajima Y, Webb WR. Air trapping on CT of patients with pulmonary embolism. AJR2002; 178:1201 -1207[Abstract/Free Full Text]
  32. Grossman ZD, Thomas FD, Gagne G, et al. Transmission computed tomographic diagnosis of experimentally produced acute pulmonary vascular occlusion in the dog. Radiology 1979;131 : 767-769[Abstract]
  33. Johnson PT, Wechsler RJ, Salazar AM, Fisher AM, Nazarian LN, Steiner RM. Spiral CT of acute pulmonary thromboembolism: evaluation of pleuroparenchymal abnormalities. J Comput Assist Tomogr 1999; 23:369 -373[CrossRef][Medline]
  34. Itti E, Nguyen S, Robin F, et al. Distribution of ventilation/perfusion ratios in pulmonary embolism: an adjunct to the interpretation of ventilation/perfusion lung scans. J Nucl Med 2002; 43:1596 -1602[Abstract/Free Full Text]
  35. Delcroix M, Melot C, Vanderhoeft P, Naeije R. Embolus size affects gas exchange in canine autologous blood clot pulmonary embolism. J Appl Physiol 1993; 74:1140 -1148[Abstract/Free Full Text]
  36. Thurlbeck WM, Müller NL. Emphysema: definition, imaging, quantification. AJR 1994;163 : 1017-1025[Abstract/Free Full Text]

Add to CiteULike CiteULike   Add to Complore Complore   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us   Add to Digg Digg   Add to Reddit Reddit   Add to Technorati Technorati    What's this?



This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Matsuoka, S.
Right arrow Articles by Nakajima, Y.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Matsuoka, S.
Right arrow Articles by Nakajima, Y.
Social Bookmarking
 Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Reddit   Add to Technorati  
What's this?


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS