Abstract
BACKGROUND AND PURPOSE: From the time-attenuation curves of DSA flow parameters, maximal intensity, maximal slope, and full width at half maximum of selected vascular points are defined. The study explores the reliability of defining the flow parameters by the time-attenuation curves of DSA.
MATERIALS AND METHODS: Seventy patients with unilateral carotid artery stenosis (group A) and 56 healthy controls (group B) were retrospectively enrolled. Fixed contrast injection protocols and DSA acquisition parameters were used with all patients. The M1, sigmoid sinus, and internal jugular vein on anteroposterior view DSA and the M2, parietal vein, and superior sagittal sinus on lateral view DSA were chosen as ROI targets for measuring flow parameters. The difference of time of maximal intensity between 2 target points was defined as the circulation time between the target points.
RESULTS: The maximal intensity difference of 2 selected points from the ICA to the M1, sigmoid sinus, internal jugular vein, M2, parietal vein, and superior sagittal sinus was significantly longer in group A than in group B. The maximum slope of M1, M2, and the superior sagittal sinus was significantly lower in group A than in group B. The full width at half maximum of M1 and M2 was significantly larger in group A than in group B. The maximal slope of M1 demonstrated the best diagnostic performance.
CONCLUSIONS: The maximal intensity difference of 2 selected points derived from DSA can be used as a definitive alternative flow parameter for intracranial circulation time measurement. Maximal slope and full width at half maximum complement the maximal intensity difference of 2 selected points in defining flow characteristics of healthy subjects and patients with carotid stenosis.
ABBREVIATIONS:
- FWHM
- full width at half maximum
- MS
- maximal slope
- PV
- parietal vein
- SSS
- superior sagittal sinus
- TDC
- time-attenuation curve
- Tmax
- time of maximal intensity
- rTmax
- Tmax difference of 2 selected points
DSA is the standard reference for diagnosing cerebrovascular diseases with its superior temporal and spatial resolution compared with other imaging methods.1⇓–3 Using different approaches, a number of recent studies have demonstrated the feasibility of quantitative flow parameter measurement by using flat detector DSA.4⇓⇓–7 Compared with optical flow methods and computer fluid dynamics simulations, flow parameter analysis with the time-attenuation curve (TDC) of DSA images is less demanding in terms of computer power and processing time.8,9 Clinically, TDCs are used to assess the “real-time” peritherapeutic hemodynamics of various vascular disorders in an angiography suite.10⇓–12 The TDC represents the dynamic intensity changes of a contrast bolus passing an ROI. It is affected by the bolus characteristics and physiologic and anatomic conditions (eg, arterial stenosis or arteriovenous shunts).11,13 From the TDC, we may measure the time of maximal intensity (Tmax), maximal slope (MS), and full width at half maximum (FWHM) of any selected vascular point on the DSA image. The time difference to reach maximum intensity (rTmax) of 2 selected vascular points indicates the circulation time between these 2 points. Accordingly, cerebral circulation time is defined as the rTmax between the internal carotid artery and the parietal vein (PV). The PV is closer to the brain parenchyma compared with the transverse sinus or jugular vein and thus better represents the time for blood flow to travel through the brain parenchyma.14 Circulation time is an objective flow parameter for various vascular disorders (eg, carotid stenosis, carotid cavernous fistula, and peritherapeutic assessment).10,12,15,16 The aim of the present study was to compare the diagnostic accuracy of rTmax, MS, and FWHM for detecting blood flow property changes by using a large sample of healthy subjects and patients with carotid stenosis.
Materials and Methods
Patient Selection
The institutional review board of the hospital approved this retrospective study, and patient consent forms were waived. From October 2011 to December 2013, 546 patients were consecutively referred to our institution for DSA. Patients with poor heart function present with polymorphic TDC waveforms and thus make identification of Tmax and subsequent measurement of MS and FWHM difficult. After we excluded patients with poor renal and/or heart functions and previous large territorial infarct and those whose imaging did not follow the standard DSA acquisition protocol, 126 patients eligible for analysis were retrospectively recruited for the current study. Seventy of 126 patients (mean age, 73.6 years; 59 men and 11 women) with unilateral extracranial internal carotid artery stenosis (>70%, based on the NASCET criteria) were classified as group A. The other 56 patients (mean age, 65.6 years; 26 men and 30 women) referred for DSA with nonarterial occlusive disorders (namely, post-aneurysm clipping/coiling follow-up and suspected vascular lesions with negative angiography findings) were classified as group B.
Imaging Protocol and Data Analysis
DSA acquisitions with a standard, clinically routine protocol were performed in all 126 cases. A power injector (Liebel-Flarsheim Angiomat; Illumena, San Diego, California) was used to create a contrast bolus after placing a 4F angiocatheter in the common carotid artery at the C4 vertebral body level. A bolus of 12 mL of 60% diluted contrast medium (340 mg I/mL) was administered within 1.5 seconds. Neither extra contrast medium nor extra radiation was used. The acquisition parameters were 7.5 frames/s for the first 5 seconds, followed by 4 frames/s for 3 seconds, 3 frames/s for 2 seconds, and finally 2 frames/s for 2 seconds. The entire DSA acquisition time was 12 seconds. However, it might be manually tailored to be shortened or prolonged for optimized internal jugular vein opacification.10 The mean and range of irradiation parameters and dosage were the following: 91.4 (88–97) kV, 287.8 (268–318) mA, 29.1 (29.1–29.2) ms, and 1.21 (0.7–1.78) mGy/frame for anteroposterior views and 73.4 (70–92) kV, 402.5 (283–426) mA, 29.1 (29.1–29.2) ms, and 1.30 (0.68–1.86) mGy/frame for lateral views.
DSAs of contralateral (normal) sides in group A were not evaluated, mainly due to different FOVs routinely used in our institution for the contralateral (normal) side of patients with unilateral internal carotid artery stenosis. There were 24 patients with hypoplasia or aplasia of the anterior communicating artery in group A. The same biplane angiography suite (Axiom Artis dBA; Siemens, Erlangen, Germany) was used for DSA throughout the entire study. The degree of arterial stenosis was determined by the more severe degree on anteroposterior lateral views according to the NASCET criteria.17 All DSA analyses were performed on a workstation equipped with prototype software, DSA Analyzer (Siemens). DSA Analyzer is a software package for analyzing 2D DSA image series. It represents an extended version of the commercially available software tool syngo iFlow (Siemens). On the basis of the TDC, DSA Analyzer extracts flow parameters (eg, Tmax, MS, and FWHM) of user-selected vascular points or ROIs on DSA.
Selection of ROIs
Parametric color-coding of DSA according to the Tmax of individual pixels was displayed instantly.4 On the basis of previous research, we placed ROIs on the first segment of the middle cerebral artery, sigmoid sinus, and ipsilateral internal jugular vein on anteroposterior view DSA and on the second branch of the middle cerebral artery, parietal vein, and superior sagittal sinus (SSS) on lateral view DSA for flow parameter analyses (Fig 1).10,14 The ROI placement was standardized to avoid overlapping anatomic structures and inhomogeneous areas and to use the caliber of the target vessel as the diameter of an ROI.10 One neuroradiologist with 8 years' experience and 1 angiographic technician performed the ROI placements by consensus. Both were unaware of the clinical conditions of the studied subjects.
Anteroposterior (A) and lateral (B) views of parametric color-coding of quantitative DSA. A, The ROI of the ICA is located at the midpoint of the cervical portion of the ICA. The ROI of M1 is located at the midpoint of the first segment of the middle cerebral artery. The sigmoid sinus ROI is located at the midpoint of the ipsilateral sigmoid sinus. The jugular vein ROI is located in the internal jugular vein at the same level as the ICA ROI. B, The ICA ROI is located in the cavernous portion of the ICA. The M2 ROI is located in the insular branch of the MCA. The PV ROI is located in the outlet of the parietal vein. The SSS ROI is located 2 cm above the confluence of the SSS.
Definition of Flow Parameters
After ROI placements, 3 flow parameters (Tmax, MS, and FWHM) were immediately extracted from the respective TDCs. “Tmax” was defined as the time point at which the pixel reached its maximum concentration during the angiographic series. “MS” was defined as the maximum tangential slope between the arrival time of the contrast medium and the Tmax on the TDC. “FWHM” was defined as the width of the 2 time points on the TDC when the concentration reached half of the maximum concentration (Fig 2). Tmax of the ROIs was normalized by subtraction of the Tmax of the cervical ICA on the anteroposterior view or cavernous ICA on the lateral view. The normalized time difference was defined as rTmax—that is, the circulation time between 2 selected vascular points. Accordingly, rTmax ICA-PV was referred to as “cerebral circulation time,” the time for blood flow to travel through the brain parenchyma.14
Time-attenuation curves of lateral view DSA in a healthy subject. The blue curve represents the TDC of the ROI in the cavernous portion of the ICA. Tmax is the time point at which the ROI reaches maximal intensity. The MS of an ROI is defined by the maximal tangential slope located between arrival time and Tmax. FWHM is the width of the waveform at the level of half maximum concentration.
Statistical Analysis
All statistical analyses were performed by using SPSS 20 (2010; IBM, Armonk, New York). The correlation between the stenotic degree and TDC flow parameters (rTmax, MS, and FWHM) was explored by using the Pearson correlations. The differences of rTmax, MS, and FWHM between groups A and B were compared by using a Student t test. We used receiver operating characteristic curves to evaluate the diagnostic performance of the above-mentioned flow parameters in carotid stenosis. The maximum area under the curve of the receiver operating characteristic was used to determine the optimal cutoff value. Significance was set to P < .05 for all statistical tests.
Results
No subjects in group A had acute stroke peritherapeutically evidenced by MR imaging. The patient demographic data are listed in Table 1. The average age in group A (73.6 ± 11.6 years) was older than that for group B (65.6 ± 10.2 years). Stenotic degree had a mild positive correlation with rTmax in M2 (r = 0.380) and M1 (r = 0.486) and FWHM in M2 (r = 0.34) and a mild negative correlation with MS in PV (−0.448) (Table 2). MS of the PV was the only venous ROI showing a negative correlation with stenotic degree.
Patient characteristics for groups A and B
Correlation of degree of stenosis with rTmax, MS, and FWHMa
The comparison of rTmax, MS, and FWHM in 6 ROIs between groups A and B is shown in Table 3. All rTmax values were significantly longer in group A than in group B. For MS, only the M1, M2, and SSS of group A demonstrated flatter slopes than in group B. For FWHM, only M1 and M2 demonstrated significantly wider waveforms in group A than in group B. None of the venous ROIs demonstrated significantly wider waveforms in group A.
Comparison of group A and B values for rTmax, MS, and FWHM for 6 different ROIs
M1 and M2 are the only ROIs with all 3 flow parameters demonstrating significant differences between 2 groups. SSS is the only venous ROI with 2 flow parameters (rTmax and MS) showing significant differences between the 2 groups.
Figure 3 illustrates the receiver operating characteristic curves of the flow parameters (6 rTmax, 3 MS, and 2 FWHM) to compare their diagnostic performance in detecting carotid stenosis. The 4 best diagnostic flow parameters were MS in M1 (cutoff value, 378.3; sensitivity, 66.7%; specificity, 60.0%), followed by FWHM in M1 (cutoff value, 3.40; sensitivity, 57.1%; specificity, 78.4%), rTmax of M1 (cutoff value, 0.45; sensitivity, 81%; specificity, 58.2%), and rTmax of SSS (cutoff value, 6.895; sensitivity, 55.2%; specificity, 73.7%) (Table 4).
Receiver operating characteristic curves of all significantly different rTmax values (A), all significantly different MS values (B), and all significantly different FWHM values (C) in differentiating patients with stenosis from control groups.
Cutoff values of 10 significant TDC parameters for detecting stenosis flow with optimized sensitivity and specificity
Discussion
Within normal physiologic conditions, the degree of carotid stenosis has greater impact on rTmax than on MS and FWHM. Carotid stenosis makes the waveform of the poststenotic arterial segment delayed (prolonged rTmax in group A) and wider (larger FWHM values in group A) (Fig 4). This effect has also been observed as delay and dispersion effects in MR imaging and CT perfusion studies.18,19 Doppler sonography illustrates the phenomenon, known as “tardus parvus,” due to the widened and flatter poststenotic waveform.20,21
TDC of lateral views of angiography are from a healthy individual (A) and a patient with 80% carotid stenosis (B). Both green lines are TDCs of M2. The rTmax of M2 in the patient with stenosis is 0.34 seconds, which was more prolonged than that (0.27 seconds) in the healthy individual. The FWHM of M2 (2.9 seconds) in the patient with stenosis is longer than that (1.8 seconds) in the healthy individual. The ICA (dark blue curves), sigmoid sinus (red curves), and internal jugular vein (light blue curves) all demonstrate prolonged rTmax and are right-shifted in the patient with stenosis. These delayed and dispersed phenomena of angiographic TDCs are consistent with those observed in MR imaging and CT perfusion imaging.
rTmax values of arteries and veins (artery-vein) were prolonged in group A (defined by the presence of carotid stenosis). Proximal arterial stenosis may transit the slowdown effects on blood flow to the downstream areas of the brain parenchyma. Poststenting normalization of rTmax artery-vein in carotid occlusive arterial disease indicates restoration of both arterial and venous blood flows.14,22 Patients who have ipsilateral stenotic and hypoplastic venous outlets are prone to develop hyperperfusion syndrome after stent placement.22⇓–24
MS was considered a surrogate marker for intravascular velocity because it represents the maximum gradient of x-ray attenuation and therefore reflects the speed of contrast medium (ie, the blood flow velocity).11,25 Theoretically, MS should be more sensitive to flow changes compared with rTmax. Nevertheless, only 2 arterial ROIs (M1, M2) and 1 venous ROI (SSS) demonstrated a significant decrease in MS in group A. As for FWHM, values were wider only in M2 and M1 in group A but not in any venous ROIs. If we take the receiver operating characteristic curves into consideration, M1 was the best location for detecting stenosis because its 3 flow parameters all showed reasonably high diagnostic performance. rTmax is still the best flow parameter among the 3 because the rTmax of all ROIs was prolonged and therefore demonstrated delayed reach of wave peaks in the TDC chart in group A (Fig 4). MS and FWHM illustrated less impact on venous flows. One plausible explanation would be that the smaller amplitude of the waveform and the encountering of more anatomic variation in the intracranial venous system make venous flow more complicated and MS and FWHM less sensitive. Additionally, the unpredictable inflow of blood without contrast from other vascular territories changes the attenuation of TDC and makes analysis of venous ROIs challenging.
Optical flow methods and computer fluid dynamics simulations quantify flows with an alternative approach, and both are highly demanding in terms of computer power and processing time. Until now, they have been available mainly in research institutes, and their roles in clinical applications are still limited. Nevertheless, computer fluid dynamics measures exclusively in vivo flow information in aneurysmal and juxta-aneurysmal areas.26,27 Its role in stenotic cerebrovascular disease is not clear. The current study results showed that MS and FWHM were only clinically significant in limited ROIs and were complementary to rTmax in describing flows.
Several limitations exist in our study. The DSA series was 2D; therefore, the overlapping of anatomic structures mixed and distorted the TDCs. Future research by using 3D DSA datasets would resolve this issue.5,28⇓⇓–31 Curve fitting or using the average slopes of ascending parts of the TDC would improve the reliability of MS and FWHM. Other potential flow parameters (eg, arrival time, kurtosis, and skewness) can describe bolus characteristics. The size and location of previous infarction might affect the measurements of the TDC but were not explored in the current study. According to our initial experience, those with old small lacunar infarcts or subterritorial infarct usually develop collateral circulation. Therefore, chronic ischemic insults might influence the rTmax or circulation, but stenotic degree still dominates in terms of TDC changes.
Our study offers reassurance about the reliability of evaluating intravascular flows by TDC and offers a real-time clinically diagnostic tool to tailor endovascular treatment within the angiography suite.26,32 Currently, the ROIs are placed manually. Future work on automation in ROI placement based on selecting the best-described waveform of flows may further shorten the assessment time and broaden the clinical application of color-coded quantitative DSA.
Conclusions
Our study confirmed that rTmax remains the preferred measurement for all ROIs when a single-flow parameter is used. MS and FWHM demonstrated equivalent diagnostic performance in differentiating carotid stenosis and thus can serve as complementary diagnostic parameters because neither alone achieves excellent sensitivity or specificity.
Acknowledgments
We thank Josephe Lavalle for manuscript preparation and. C.S. Lin for statistical support.
Footnotes
Disclosures: Wan-Yuo Guo—UNRELATED: Grant: collaboration between Institution and Siemens*; Support for Travel to Meetings for the Study or Other Purposes: Taipei Veterans General Hospital,* Comments: to present the study results; Provision of Writing Assistance, Medicines, Equipment, or Administrative Support: software provided by Siemens under the collaboration contract. Janina Beilner—OTHER RELATIONSHIPS: In this study, Siemens software was used. I am a Siemens employee. Markus Kowarschik—OTHER RELATIONSHIPS: Employee of Siemens. *Money paid to the institution.
This work was cosponsored by the Taipei Veterans General Hospital and Siemens (grant No: T1100200).
References
- Received August 10, 2014.
- Accepted after revision December 1, 2014.
- © 2015 by American Journal of Neuroradiology