Quantification of Blood Velocity with 4D Digital Subtraction Angiography Using the Shifted Least-Squares Method

Y. Wu; G. Shaughnessy; C.A. Hoffman; E.L. Oberstar; S. Schafer; T. Schubert; K.L. Ruedinger; B.J. Davis; C.A. Mistretta; C.M. Strother; M.A. Speidel

doi:10.3174/ajnr.A5793

Abstract

BACKGROUND AND PURPOSE: 4D-DSA provides time-resolved 3D-DSA volumes with high temporal and spatial resolutions. The purpose of this study is to investigate a shifted least squares method to estimate the blood velocity from the 4D DSA images. Quantitative validation was performed using a flow phantom with an ultrasonic flow probe as ground truth. Quantification of blood velocity in human internal carotid arteries was compared with measurements generated from 3D phase-contrast MR imaging.

MATERIALS AND METHODS: The centerlines of selected vascular segments and the time concentration curves of each voxel along the centerlines were determined from the 4D-DSA dataset. The temporal shift required to achieve a minimum difference between any point and other points along the centerline of a segment was calculated. The temporal shift as a function of centerline point position was fit to a straight line to generate the velocity. The proposed shifted least-squares method was first validated using a flow phantom study. Blood velocities were also estimated in the 14 ICAs of human subjects who had both 4D-DSA and phase-contrast MR imaging studies. Linear regression and correlation analysis were performed on both the phantom study and clinical study, respectively.

RESULTS: Mean velocities of the flow phantom calculated from 4D-DSA matched very well with ultrasonic flow probe measurements with 11% relative root mean square error. Mean blood velocities of ICAs calculated from 4D-DSA correlated well with phase-contrast MR imaging measurements with Pearson correlation coefficient r = 0.835.

CONCLUSIONS: The availability of 4D-DSA provides the opportunity to use the shifted least-squares method to estimate velocity in vessels within a 3D volume.

ABBREVIATIONS:

PC: phase-contrast
SBR: sideband ratio
TCC: time concentration curve
VIPR: vastly undersampled isotropic projection reconstruction

2D digital subtraction angiography and a 3D rotational angiography acquisition followed by 3D-DSA reconstruction are the standards for vascular morphology assessment. There is an increasing demand for hemodynamic information, including blood flow rate and velocity for diagnosis, treatment planning, and evaluation.^1⇓–3 Many algorithms have been proposed to estimate blood flow in arteries from 2D-DSA; these are divided into 2 major classes: bolus-tracking and computational methods.^{4⇓⇓⇓⇓⇓⇓–11} A thorough review of these techniques can be found in Shpilfoygel et al.¹² However, 2D-DSA based methods are challenged by vessel overlap and vessels being aligned with the primary x-ray beam direction.

Recently, a new reconstruction method^{13⇓⇓–16} has been applied to a single 3D rotational angiography acquisition to generate time-resolved 3D-DSA at frame rates up to 30 frames/second (4D DSA). The availability of the geometric and temporal data in a 4D-DSA reconstruction provides the opportunity to estimate velocity and flow more accurately than has previously been possible with 2D DSA. In this article, a shifted least-squares based technique was used with the 4D-DSA data to estimate blood velocity. The proposed algorithm was first validated using flow phantom studies, in which the velocity could be documented using an ultrasonic flow probe. Blood velocity was also determined in 4D-DSA datasets from 14 ICAs of volunteers who also underwent 3D fast phase-contrast with vastly undersampled isotropic projection reconstruction (PC VIPR) studies.^17⇓–19

Materials and Methods

4D DSA Acquisition and Data Preparation

Noncontrast (mask) and contrast-enhanced (fill) rotational acquisitions were performed using a conventional flat panel detector angiographic system (Artis zee; Siemens, Erlangen, Germany) with the following settings: 70 kV(peak), 0.36 μGy/frame, 304 images, 260° arc. Contrast injection (iohexol, Omnipaque 300; GE Healthcare, Piscataway, New Jersey) into the subject was started at the beginning of the fill run (3 mL/second, 7-second injection duration). The rotational acquisition protocol yields x-ray projection images, which are used to automatically reconstruct both 3D-DSA and 4D-DSA with an isotropic spatial resolution of 0.46 × 0.46 × 0.46 mm³.

The centerline of each vascular segment was determined from the static 3D volume using an efficient 3D parallel-thinning algorithm.²⁰ The contrast waveform map, which is the time concentration curves (TCCs) along the centerline path (curvilinear length z), was extracted from the 4D-DSA and stored for each vascular segment for velocity calculation.

Shifted Least-Squares Algorithm

In an arterial contrast injection, there is a temporal oscillation in iodine concentration that arises from the mixing of contrast medium, which is injected at a fixed rate, and nonopacified blood, which flows at a variable rate driven by the cardiac cycle. This temporal variation in contrast, referred to as pulsatility, appears at points downstream of the injection with a time delay. The time delay is related to the distance downstream and the blood velocity. Therefore, velocity can be estimated by measuring the distance along the vessel centerline and the time delay.

For any 2 points, i and j, along the centerline, the shifted least-squares difference between their TCCs, c_i(t) and c_j(t), can be calculated as The value τ_ij⁰, τ_ij⁰ ∈ [0, T] that minimizes the shifted least-squares difference, ε (τ_ij), is regarded as the time of bolus transport between these 2 points. The spatial distance z between the 2 points is calculated from the 3D path length along the vessel centerline. The time-shift τ⁰ as a function of the spatial distance z can then be fit to a linear relation: where the slope a is the inverse of the velocity v.

Optimizing Waveform Selection

It has been found that the shifted least-squares algorithm provides the most reliable velocity calculation in vessel segments where pulsatility is strong and consistent. To automatically select the waveform regions with strong and consistent pulsatility, we generated a sideband ratio (SBR) map.²¹

For each point in the waveform map, a short-time Fourier transform was applied to generate the local power spectrum PS_fn, where f_n is the n^th characteristic frequency. The SBR at this point can be calculated as where the denominator summation includes the neighboring frequency points and excludes the fundamental frequency f₀. From this SBR map, a low-pass filter was applied and a threshold of one-fourth of the median SBR was enforced to generate a mask. This mask excludes low pulsatile signals from consideration in the velocity calculation.

Phantom Study

To validate the proposed method for velocity estimation, we used a flow phantom containing a loop of plastic tubing to perform 4D-DSA studies. Pulsatile water flow with different mean flow rates was generated using a perfusion pump with a controllable flow setting (S3; Storkcert, Freiburg, Germany). A 5F catheter was introduced into the tubing for injection of iodinated contrast medium (Omnipaque 300 diluted with 25% water) with a power injector. 4D-DSA-derived flow was compared with flow measurements made with an ultrasonic flow probe (400-Series Multi-Channel Flowmeter; Transonic, Ithaca, New York) attached to the tubing 5 inches away from the catheter tip.

Rotational projection images were acquired using the mask and iodine-filled scans with 12-second acquisition over a 260° arc. The frame rate was 30 frames per second. The 4D-DSA volume data were reconstructed using the vendor's prototype software (Siemens Healthineers, Forchheim, Germany).

A time-resolved flow profile was recorded during each fill scan using the ultrasonic probe. The reference flow measurement was taken as the average of the flow profile over the same time window that was used for 4D-DSA-derived flow.

Five different flow settings were selected to provide average flow rates ranging from 285 to 1140 mL/min (corresponding to the mean velocity of 15–60 cm/s for the tubing with a 0.64-cm diameter). Three 4D-DSA scans and reference flow measurements were made for each flow setting.

Clinical Study

Nine subjects who had both DSA and PC VIPR examinations were selected from a data base generated under University of Wisconsin review board approval. Because 5 of 9 subjects had bilateral studies, a total of 14 ICAs were available for evaluation. Scanning parameters for postcontrast PC VIPR were the following: FOV = 22 × 22 × 22 cm², TR/TE = 12.5/4.8 ms, velocity-encoded gradient-echo imaging = 80~100 cm/s, bandwidth = 83.3 kHz, readout matrix = 320 points per projections, spatial resolution for the complex difference image = 0.68 × 0.68 × 0.68 mm³. The complex difference images from the PC VIPR scan were segmented for vessel content using Materialise Mimics (Materialise NV, Leuven, Belgium). The segmented MR imaging data were then registered to the 3D-DSA using a correlated point registration with manual point placement, followed by an affine rigid registration. Centerline coordinates of the ICAs generated from the 3D-DSA volumes were used for velocity estimation. An automated workflow for the PC VIPR analysis was developed in Matlab (MathWorks, Natick, Massachusetts). This software performed cross-sectional plane placement for every voxel along the 3D centerline path, segmentation of the vessel boundaries in each plane, and averaging of the normal velocities for all points inside the segmented boundaries. The reported velocity value for a vessel was then the average result over all the centerline points that were present in both the PC VIPR and 3D-DSA volumes.

Statistical Analysis

The validation process was performed using linear regression fit between flow probe (V_fp) and V_dsa (calculated using 4D-DSA) for the phantom study. Correlation analysis and Bland-Altman plots were used to compare the velocity measurements using PC VIPR versus 4D-DSA for the clinical study. The Pearson r was used for correlation analysis, and a P value < .05 for the correlation coefficient was considered statistically significant.

Results

Phantom Study

Figure 1 shows sample data from the phantom study. First, the centerline was generated from the static 3D volume (Fig 1A). Figure 1B is the contrast waveform map (TCCs along the centerline path) of the selected centerline. Each point in the map M(t,z) represents the signal intensity of the voxel at distance z along the centerline and at timeframe t [C_z(t)]. The intensity variations along the vertical axis reflect the concentration of iodine contrast medium versus the distance from the injection site for a given time point. Variations along the horizontal axis reflect variations in concentration versus time for a given point along the centerline. The slope of the intensity ridges reflects the blood velocity. In the displayed format, higher blood velocities produce ridges that are closer to vertical, whereas lower blood velocities produce ridges closer to horizontal. The highlighted region is the optimized waveform region for the flow calculation determined using the SBR method. TCCs of 2 selected voxels along the centerline (marked as blue and red stars on Fig 1A and blue and red dashed lines on Fig 1B) are shown in Fig 1C. Figure 1D shows the least-squares differences of these 2 signals as a function of the time-shift, where the minimal appears at τ⁰ = 5 is considered as the time of bolus transport from the blue voxel to the red voxel. Figure 1F demonstrates the linear fitting between the time-shift of each voxel along the centerline and its corresponding centerline position. Figure 1E is the flow profile recorded from the flow probe (downsampled to 1/30 second to correspond to the 4D-DSA timeframes). Flow measurement was taken as an average over the same temporal window used in the 4D-DSA flow calculation.

Fig 1.

A representative phantom study demonstrates the flow calculation from 4D-DSA. A, MIP image of the reconstructed phantom with centerline (red) overlaid. B, The contrast waveform map of the voxels along the centerline. The horizontal direction is the timeframe in the 4D-DSA scan, and the vertical direction is the position along the vessel centerline. Each point in the map M(t,z) represents the signal intensity of the voxel at distance z along the centerline and at time frame t. The highlighted area is the optimized waveform region where the pulsatility is strong and consistent. TCCs of 2 selected voxels along the centerline (marked as blue and red stars on A and blue and red curves on B) are shown in C. D, The least-squares differences of these 2 signals as a function of the time-shift, in which the minimal appears at τ⁰ = 5 shown as red arrow on D is considered as the time of bolus transport from the blue voxel to the red voxel. The time-shift as a function of the centerline position was fit to a linear relation (F), where the slope is the inverse of the velocity. E, The flow profile recorded from the flow probe (downsampled to 1/30 second to correspond to the 4D-DSA timeframes). Flow measurement is the average taken between the red lines, which correspond to the optimized window shown in B. AUI indicates arbitrary unit of intensity.

The linear regression fit between the flow probe measurement and the 4D-DSA calculation is shown in Fig 2. The linear regression equation is as follows: V_dsa = 1.0069V_fp + 2.8482, where V_dsa is the velocity calculated from 4D-DSA and V_fp is the measurement using the flow probe. The determination coefficient R² of the fit was 0.9677. The 95% confidence interval slope was 0.8966–1.1172 and the 95% CI intercept was −1.4917–7.1881. These evaluations have been obtained with significance levels of P = 4.521E-11 and P = 0.18, respectively. The relative root mean square error (normalized by mean velocity) was 11.3%.

Fig 2.

Linear regression between the flow probe measurements and the 4D-DSA calculations from 15 phantom studies. Evaluations have been obtained with significance levels P = 4.521E-11 for the slope and P = 0.18 for the intercept.

Clinical Study

The average ICA flow rates were calculated from 13 of the 14 ICAs that had 4D-DSA examinations. One subject was excluded due to the poor pulsatility in the TCCs. Figure 3 shows an example of the ICA studies. Average ICA flow rates were also measured with the PC VIPR studies. Correlation between 2 measurements using 4D-DSA and PC VIPR is shown in Fig 4. Average ICA flow rates calculated using 4D-DSA and measured using PC VIPR are well-correlated, with the Pearson r = 0.835 and a significance level P = .0002. Figure 5 is the Bland-Altman analysis of 4D-DSA compared with PC VIPR, which shows all points within 2 SDs, with a bias of 6.5 cm/s and 2-SD limits of agreement of −3.56 cm/s and 16.5 cm/s.

Fig 3.

A representative ICA study. A, The centerline positions of the ICA. B, The contrast waveform map of the voxels along the centerline. The highlighted area is the optimized waveform region for the flow calculation. C, The time-shift as a function of the centerline position was fit to a linear relation in which the slope is the inverse of the velocity.

Fig 4.

Correlation analysis between velocity estimated from 4D-DSA and measured using PC VIPR. The significance level is P = .0002.

Fig 5.

Bland-Altman analysis of 4D-DSA compared with PC VIPR.

Discussion

The shifted least-squares algorithm is a simple-but-robust method to calculate flow rates by finding the bolus transport time between any 2 points along the centerline and performing linear regression with the corresponding positions. The combination of the spatial and temporal data available in the 4D-DSA reconstruction provides the opportunity to apply this technique to measurement of blood flow in routine 4D-DSA acquisitions. The flow phantom study demonstrated that by using this method, we were able to calculate the flow rates with 11% relative root mean square error compared with ultrasonic flow probe measurements. Results from the 13 human ICAs showed that the calculated velocities from 4D-DSA studies correlated well with the measurements using the 3D phase-contrast MR imaging method. We believe that the reason for the higher velocities in the 4D-DSA studies compared with the PC VIPR measurements (6.5 cm/s bias) is the increment in flow caused by the intra-arterial power injection of contrast medium.

The proposed method calculates the mean flow rate within a selected window defined by the SBR map. This truncated average of pulsatile flow varies depending on the window size and location. This variation can be demonstrated in Fig 6A and B, where the average of the flow profile was performed with a gradually expanding average window. The truncated average oscillates depending on where the truncation window ends and eventually converges to the mean signal as the window expands. Similar oscillation has been observed in the flow calculation using 4D-DSA. Figure 6D shows the estimated average velocity changes as the right side of the window gradually expands as shown in Fig 6C. The current algorithm automatically selects a window purely on the basis of the SBR map so that the selected window range may fall anywhere between 4 and 6 pulsatile cycles. Further criteria could be included to trim the window to an integer number of pulsatile cycles to minimize this error.

Fig 6.

Fluctuation of the velocity estimation with a different average range. A, The flow profile recorded from the flow probe. Red lines define the starting average window. This window was gradually expanded to the purple line. The corresponding mean velocity was shown as the blue curve. A zoom-in of the green window is shown in B. Similarly, velocity calculations have been performed by gradually expanding the window coverage from yellow to purple as shown in the waveform map (C). Estimated velocity varies (D) as the window edge slides from position A toward B. AU indicates arbitrary units.

The proposed algorithm worked best when there was good and consistent pulsatility in the 4D-DSA TCCs. With each cardiac cycle, as an intra-arterial contrast bolus is injected, there is an oscillation in the ratio of contrast medium to nonopacified blood. This pulsatility can be identified and extracted from the contrast waveform map by thresholding the SBR map. For vessels more distal to the contrast injection site where the pulsatility can be low, injection protocols should be modified to enhance the pulsatility strength or contrast kinetics. The present work focused on the evaluation of blood velocity within the ICAs. Investigation of more downstream segments, including the MCAs and basilar arteries, should be performed in the future.

Conclusions

4D DSA is a novel technique that provides both the geometric and temporal information required to measure blood flow velocity and flow. In this study, calculated velocity values from 4D-DSA TCCs correlated well with human ICA measurements derived from phase-contrast MR imaging and flow phantom measurements using an ultrasonic flow probe.

Footnotes

Disclosures: Gabe Shaughnessy—RELATED: Grant: National Institutes of Health.* Erick L. Oberstar—RELATED: Grant: University of Wisconsin, Madison, Comments: National Institutes of Health Grant R01HL116567*. Sebastian Schafer—UNRELATED: Employment: Siemens Healthineers, Comments: full-time employee. Michael A. Speidel—RELATED: Grant: Siemens Healthineers, Comments: Partial support was received from a research grant*; UNRELATED: Patents (Planned, Pending or Issued): Wisconsin Alumni Research Foundation, Comments: I am a coinventor on a Wisconsin Alumni Research Foundation patent (pending) on resolving artifacts in 4D angiographic data. Money of the patents was paid to Wisconsin Alumni Research Foundation. Charles A. Mistretta—RELATED: Grant: National Institute of Biomedical Imaging and Bioengineering, Comments: Strother/Mistretta R01 on 4DDSA*;UNRELATED: Patents (Planned, Pending or Issued): Wisconsin Alumni Research Foundation. Charles M. Strother—RELATED: Grant: University of Wisconsin School of Medicine and Public Health*; UNRELATED: Patents (Planned, Pending or Issued): Wisconsin Alumni Research Foundation, Comments: I have a pending IP for flow quantification; the IP is assigned to the Wisconsin Alumni Research Foundation. *Money paid to the institution.
This work was supported by the National Institutes of Health Center for Scientific Review 2R01HL116567 and Siemens Healthineers.
The concepts and information presented in this article are based on research and are not commercially available.

Indicates open access to non-subscribers at www.ajnr.org

References

1.↵
2. Murayama Y,
3. Usami S,
4. Abe T, et al
. Transvenous Doppler guidewire sonographic monitoring during treatment of a complex vertebral arteriovenous fistula associated with neurofibromatosis type 1. Neuroradiology 1999;41:328–33 doi:10.1007/s002340050758 pmid:10379588
CrossRef PubMed
2.↵
2. Mast H,
3. Mohr J,
4. Thompson J, et al
. Transcranial Doppler ultrasonography in cerebral arteriovenous malformations: diagnostic sensitivity and association of flow velocity with spontaneous hemorrhage and focal neurological deficit. Stroke 1995;26:1024–27 doi:10.1161/01.STR.26.6.1024 pmid:7762018
Abstract/FREE Full Text
3.↵
2. Turski P,
3. Levine R,
4. Turnipseed W, et al
. MR angiography flow analysis: neurovascular applications. Magn Reson Imaging Clin N Am 1995;3:541–55 pmid:7584256
PubMed
4.↵
2. Rosen L,
3. Silverman NR
, Videodensitometric measurements of blood flow using crosscorrelation techniques. Radiology 1973;109:305–10 doi:10.1148/109.2.305 pmid:4742317
CrossRef PubMed
5.↵
2. Silverman NR,
3. Rosen L
. Arterial blood flow measurement: assessment of velocity estimation methods. Invest Radiol 1977:12:319–24 doi:10.1097/00004424-197707000-00004 pmid:885690
CrossRef PubMed
6.↵
2. Seifalian AM,
3. Hawkes DJ,
4. Colchester AC, et al
. A new algorithm for deriving pulsatile blood flow waveforms tested using simulated dynamic angiographic data, Neuroradiology 1989;31:263–69 pmid:2779778
CrossRef PubMed
7.↵
2. Hoffmann KR,
3. Doi K,
4. Fencil LE
. Determination of instantaneous and average blood flow rates from digital angiograms of vessel phantoms using distance-density curves. Invest Radiol 1991;26:207–12 pmid:2055725
CrossRef PubMed
8.↵
2. Imbert B,
3. Meunier J,
4. Mongrain R, et al
. Stenosis parameter assessment from contrast medium tracking in cineangiography with an optical flow method. Proc SPIE Int Soc Opt Eng 1997;l3034:631–40 doi:10.1117/12.274149
CrossRef
9.↵
2. Pereira VM,
3. Bonnefous O,
4. Ouared R, et al
. A DSA-based method using contrast-motion estimation for the assessment of the intra-aneurysmal flow changes induced by flow-diverter stents. AJNR Am J Neuroradiol 2013;34:808–15 doi:10.3174/ajnr.A3322 pmid:23124641
Abstract/FREE Full Text
10.↵
2. Rhode KS,
3. Lambrou T,
4. Hawkes DJ, et al
. Novel approaches to the measurement of arterial blood flow from dynamic digital X-ray images. IEEE Trans Med Imaging 2005;24:500–13 doi:10.1109/TMI.2004.843202 pmid:15822808
CrossRef PubMed
11.↵
2. Pereira VM,
3. Ouared R,
4. Brina O, et al
. Quantification of internal carotid artery flow with digital subtraction angiography: validation of an optical flow approach with Doppler ultrasound. AJNR Am J Neuroradiol 2014;35:156–63 doi:10.3174/ajnr.A3662 pmid:23928145
Abstract/FREE Full Text
12.↵
2. Shpilfoygel SD,
3. Close RA,
4. Valentino DJ, et al
. X-ray videodensitometric methods for blood flow and velocity measurement: a critical review of literature. Med Phys 2000;27:2008–23 doi:10.1118/1.1288669 pmid:11011728
CrossRef PubMed
13.↵
2. Mistretta CA,
3. Oberstar EL,
4. Davis JB, et al
. 4D DSA and 4D fluoroscopy: preliminary implementation. Proc SPIE Int Soc Opt Eng 2010;7622 doi:10.1117/12.844561
CrossRef
14.↵
2. Strother CM,
3. Bender F,
4. Deuerling-Zheng Y, et al
. Parametric color coding of digital subtraction angiography. AJNR Am J Neuroradiol 2010;31:919–24 doi:10.3174/ajnr.A2020 pmid:20167651
Abstract/FREE Full Text
15.↵
2. Benndorf G
. Color-coded digital subtraction angiography: the end of a monochromatic era? AJNR Am J Neuroradiol 2010;31:925–27 doi:10.3174/ajnr.A2077 pmid:20395396
FREE Full Text
16.↵
2. Davis B,
3. Royalty K,
4. Kowarschik M, et al
. 4D digital subtraction angiography: implementation and demonstration of feasibility. AJNR Am J Neuroradiol 2013;34:1914–21 doi:10.3174/ajnr.A3529 pmid:23620072
Abstract/FREE Full Text
17.↵
2. Gu T,
3. Korosec FR,
4. Block W, et al
. PC VIPR: a high-speed 3D phase-contrast method for flow quantification and high-resolution angiography. AJNR Am J Neuroradiol 2005;26:743–49 pmid:15814915
Abstract/FREE Full Text
18.↵
2. Chang W,
3. Landgraf B,
4. Johnson KM, et al
. Velocity measurements in the middle cerebral arteries of healthy volunteers using 3D radial phase-contrast HYPRFlow: comparison with transcranial Doppler sonography and 2D phase-contrast MR imaging. AJNR Am J Neuroradiol 2011;32:54–59 doi:10.3174/ajnr.A2240 pmid:20947642
Abstract/FREE Full Text
19.↵
2. Schneiders JJ,
3. Ferns SP,
4. van Ooij P, et al
. Comparison of phase-contrast MR imaging and endovascular sonography for intracranial blood flow velocity measurements. AJNR Am J Neuroradiol 2012;33:1786–90 doi:10.3174/ajnr.A3142 pmid:22576898
Abstract/FREE Full Text
20.↵
2. Lee TC,
3. Kashyap RL,
4. Chu CN
. Building skeleton models via 3-D medial surface axis thinning algorithms. Graphical Models and Image Processing 1994;56:462–78 doi:10.1006/cgip.1994.1042
CrossRef
21.↵
2. Shaughnessy G,
3. Hoffman C,
4. Schafer S, et al
. Quantitative flow and velocity measurements of pulsatile blood flow with 4D DSA. Proc SPIE 10132, Medical Imaging 2017: Physics of Medical Imaging, 101325R (19 April 2017); doi: 10.1117/12.2254143 doi:10.1117/12.2254143
CrossRef

Received February 6, 2018.
Accepted after revision June 11, 2018.

View Abstract

In this issue

Download PDF

Email Article

Citation Tools

Purchase

Cited By...

More in this TOC Section

Show more INTERVENTIONAL

[1] 1.↵

Murayama Y,
Usami S,
Abe T, et al
. Transvenous Doppler guidewire sonographic monitoring during treatment of a complex vertebral arteriovenous fistula associated with neurofibromatosis type 1. Neuroradiology 1999;41:328–33 doi:10.1007/s002340050758 pmid:10379588
CrossRef PubMed

[3] Murayama Y,

[4] Usami S,

[5] Abe T, et al

[6] 2.↵

Mast H,
Mohr J,
Thompson J, et al
. Transcranial Doppler ultrasonography in cerebral arteriovenous malformations: diagnostic sensitivity and association of flow velocity with spontaneous hemorrhage and focal neurological deficit. Stroke 1995;26:1024–27 doi:10.1161/01.STR.26.6.1024 pmid:7762018
Abstract/FREE Full Text

[8] Mast H,

[9] Mohr J,

[10] Thompson J, et al

[11] 3.↵

Turski P,
Levine R,
Turnipseed W, et al
. MR angiography flow analysis: neurovascular applications. Magn Reson Imaging Clin N Am 1995;3:541–55 pmid:7584256
PubMed

[13] Turski P,

[14] Levine R,

[15] Turnipseed W, et al

[16] 4.↵

Rosen L,
Silverman NR
, Videodensitometric measurements of blood flow using crosscorrelation techniques. Radiology 1973;109:305–10 doi:10.1148/109.2.305 pmid:4742317
CrossRef PubMed

[18] Rosen L,

[19] Silverman NR

[20] 5.↵

Silverman NR,
Rosen L
. Arterial blood flow measurement: assessment of velocity estimation methods. Invest Radiol 1977:12:319–24 doi:10.1097/00004424-197707000-00004 pmid:885690
CrossRef PubMed

[22] Silverman NR,

[23] Rosen L

[24] 6.↵

Seifalian AM,
Hawkes DJ,
Colchester AC, et al
. A new algorithm for deriving pulsatile blood flow waveforms tested using simulated dynamic angiographic data, Neuroradiology 1989;31:263–69 pmid:2779778
CrossRef PubMed

[26] Seifalian AM,

[27] Hawkes DJ,

[28] Colchester AC, et al

[29] 7.↵

Hoffmann KR,
Doi K,
Fencil LE
. Determination of instantaneous and average blood flow rates from digital angiograms of vessel phantoms using distance-density curves. Invest Radiol 1991;26:207–12 pmid:2055725
CrossRef PubMed

[31] Hoffmann KR,

[32] Doi K,

[33] Fencil LE

[34] 8.↵

Imbert B,
Meunier J,
Mongrain R, et al
. Stenosis parameter assessment from contrast medium tracking in cineangiography with an optical flow method. Proc SPIE Int Soc Opt Eng 1997;l3034:631–40 doi:10.1117/12.274149
CrossRef

[36] Imbert B,

[37] Meunier J,

[38] Mongrain R, et al

[39] 9.↵

Pereira VM,
Bonnefous O,
Ouared R, et al
. A DSA-based method using contrast-motion estimation for the assessment of the intra-aneurysmal flow changes induced by flow-diverter stents. AJNR Am J Neuroradiol 2013;34:808–15 doi:10.3174/ajnr.A3322 pmid:23124641
Abstract/FREE Full Text

[41] Pereira VM,

[42] Bonnefous O,

[43] Ouared R, et al

[44] 10.↵

Rhode KS,
Lambrou T,
Hawkes DJ, et al
. Novel approaches to the measurement of arterial blood flow from dynamic digital X-ray images. IEEE Trans Med Imaging 2005;24:500–13 doi:10.1109/TMI.2004.843202 pmid:15822808
CrossRef PubMed

[46] Rhode KS,

[47] Lambrou T,

[48] Hawkes DJ, et al

[49] 11.↵

Pereira VM,
Ouared R,
Brina O, et al
. Quantification of internal carotid artery flow with digital subtraction angiography: validation of an optical flow approach with Doppler ultrasound. AJNR Am J Neuroradiol 2014;35:156–63 doi:10.3174/ajnr.A3662 pmid:23928145
Abstract/FREE Full Text

[51] Pereira VM,

[52] Ouared R,

[53] Brina O, et al

[54] 12.↵

Shpilfoygel SD,
Close RA,
Valentino DJ, et al
. X-ray videodensitometric methods for blood flow and velocity measurement: a critical review of literature. Med Phys 2000;27:2008–23 doi:10.1118/1.1288669 pmid:11011728
CrossRef PubMed

[56] Shpilfoygel SD,

[57] Close RA,

[58] Valentino DJ, et al

[59] 13.↵

Mistretta CA,
Oberstar EL,
Davis JB, et al
. 4D DSA and 4D fluoroscopy: preliminary implementation. Proc SPIE Int Soc Opt Eng 2010;7622 doi:10.1117/12.844561
CrossRef

[61] Mistretta CA,

[62] Oberstar EL,

[63] Davis JB, et al

[64] 14.↵

Strother CM,
Bender F,
Deuerling-Zheng Y, et al
. Parametric color coding of digital subtraction angiography. AJNR Am J Neuroradiol 2010;31:919–24 doi:10.3174/ajnr.A2020 pmid:20167651
Abstract/FREE Full Text

[66] Strother CM,

[67] Bender F,

[68] Deuerling-Zheng Y, et al

[69] 15.↵

Benndorf G
. Color-coded digital subtraction angiography: the end of a monochromatic era? AJNR Am J Neuroradiol 2010;31:925–27 doi:10.3174/ajnr.A2077 pmid:20395396
FREE Full Text

[71] Benndorf G

[72] 16.↵

Davis B,
Royalty K,
Kowarschik M, et al
. 4D digital subtraction angiography: implementation and demonstration of feasibility. AJNR Am J Neuroradiol 2013;34:1914–21 doi:10.3174/ajnr.A3529 pmid:23620072
Abstract/FREE Full Text

[74] Davis B,

[75] Royalty K,

[76] Kowarschik M, et al

[77] 17.↵

Gu T,
Korosec FR,
Block W, et al
. PC VIPR: a high-speed 3D phase-contrast method for flow quantification and high-resolution angiography. AJNR Am J Neuroradiol 2005;26:743–49 pmid:15814915
Abstract/FREE Full Text

[79] Gu T,

[80] Korosec FR,

[81] Block W, et al

[82] 18.↵

Chang W,
Landgraf B,
Johnson KM, et al
. Velocity measurements in the middle cerebral arteries of healthy volunteers using 3D radial phase-contrast HYPRFlow: comparison with transcranial Doppler sonography and 2D phase-contrast MR imaging. AJNR Am J Neuroradiol 2011;32:54–59 doi:10.3174/ajnr.A2240 pmid:20947642
Abstract/FREE Full Text

[84] Chang W,

[85] Landgraf B,

[86] Johnson KM, et al

[87] 19.↵

Schneiders JJ,
Ferns SP,
van Ooij P, et al
. Comparison of phase-contrast MR imaging and endovascular sonography for intracranial blood flow velocity measurements. AJNR Am J Neuroradiol 2012;33:1786–90 doi:10.3174/ajnr.A3142 pmid:22576898
Abstract/FREE Full Text

[89] Schneiders JJ,

[90] Ferns SP,

[91] van Ooij P, et al

[92] 20.↵

Lee TC,
Kashyap RL,
Chu CN
. Building skeleton models via 3-D medial surface axis thinning algorithms. Graphical Models and Image Processing 1994;56:462–78 doi:10.1006/cgip.1994.1042
CrossRef

[94] Lee TC,

[95] Kashyap RL,

[96] Chu CN

[97] 21.↵

Shaughnessy G,
Hoffman C,
Schafer S, et al
. Quantitative flow and velocity measurements of pulsatile blood flow with 4D DSA. Proc SPIE 10132, Medical Imaging 2017: Physics of Medical Imaging, 101325R (19 April 2017); doi: 10.1117/12.2254143 doi:10.1117/12.2254143
CrossRef

[99] Shaughnessy G,

[100] Hoffman C,

[101] Schafer S, et al

Main menu

User menu

Search

American Journal of Neuroradiology

Quantification of Blood Velocity with 4D Digital Subtraction Angiography Using the Shifted Least-Squares Method

Abstract

ABBREVIATIONS: