Absolute volumetric blood flow measurements using dual-energy digital subtraction angiography

Dynamic CT perfusion measurement in a cardiac phantom

Widespread clinical implementation of dynamic CT myocardial perfusion has been hampered by its limited accuracy and high radiation dose. The purpose of this study was to evaluate the accuracy and radiation dose reduction of a dynamic CT myocardial perfusion technique based on first pass analysis (FPA). To test the FPA technique, a pulsatile pump was used to generate known perfusion rates in a range of 0.96-2.49 mL/min/g. All the known perfusion rates were determined using an ultrasonic flow probe and the known mass of the perfusion volume. FPA and maximum slope model (MSM) perfusion rates were measured using volume scans acquired from a 320-slice CT scanner, and then compared to the known perfusion rates. The measured perfusion using FPA (P(FPA)), with two volume scans, and the maximum slope model (P(MSM)) were related to known perfusion (P(K)) by P(FPA) = 0.91P(K) + 0.06 (r = 0.98) and P(MSM) = 0.25P(K) - 0.02 (r = 0.96), respectively. The standard error of estimate for the FPA technique, using two volume scans, and the MSM was 0.14 and 0.30 mL/min/g, respectively. The estimated radiation dose required for the FPA technique with two volume scans and the MSM was 2.6 and 11.7-17.5 mSv, respectively. Therefore, the FPA technique can yield accurate perfusion measurements using as few as two volume scans, corresponding to approximately a factor of four reductions in radiation dose as compared with the currently available MSM. In conclusion, the results of the study indicate that the FPA technique can make accurate dynamic CT perfusion measurements over a range of clinically relevant perfusion rates, while substantially reducing radiation dose, as compared to currently available dynamic CT perfusion techniques.

Estimation of regional myocardial mass at risk based on distal arterial lumen volume and length using 3D micro-CT images

The determination of regional myocardial mass at risk distal to a coronary occlusion provides valuable prognostic information for a patient with coronary artery disease. The coronary arterial system follows a design rule which allows for the use of arterial branch length and lumen volume to estimate regional myocardial mass at risk. Image processing techniques, such as segmentation, skeletonization and arterial network tracking, are presented for extracting anatomical details of the coronary arterial system using micro-computed tomography (micro-CT). Moreover, a method of assigning tissue voxels to their corresponding arterial branches is presented to determine the dependent myocardial region. The proposed micro-CT technique was utilized to investigate the relationship between the sum of the distal coronary arterial branch lengths and volumes to the dependent regional myocardial mass using a polymer cast of a porcine heart. The correlations of the logarithm of the total distal arterial lengths (L) to the logarithm of the regional myocardial mass (M) for the left anterior descending (LAD), left circumflex (LCX) and right coronary (RCA) arteries were log(L)=0.73log(M)+0.09 (R=0.78), log(L)=0.82log(M)+0.05 (R=0.77) and log(L)=0.85log(M)+0.05 (R=0.87), respectively. The correlation of the logarithm of the total distal arterial lumen volumes (V) to the logarithm of the regional myocardial mass for the LAD, LCX and RCA were log(V)=0.93log(M)-1.65 (R=0.81), log(V)=1.02log(M)-1.79 (R=0.78) and log(V)=1.17log(M)-2.10 (R=0.82), respectively. These morphological relations did not change appreciably for diameter truncations of 600-1400microm. The results indicate that the image processing procedures successfully extracted information from a large 3D dataset of the coronary arterial tree to provide prognostic indications in the form of arterial tree parameters and anatomical area at risk.

Feasibility of real time dual-energy imaging based on a flat panel detector for coronary artery calcium quantification

The feasibility of a real-time dual-energy imaging technique with dynamic filtration using a flat panel detector for quantifying coronary arterial calcium was evaluated. In this technique, the x-ray beam was switched at 15 Hz between 60 kVp and 120 kVp with the 120 kVp beam having an additional 0.8 mm silver filter. The performance of the dynamic filtration technique was compared with a static filtration technique (4 mm Al+0.2 mm Cu for both beams). The ability to quantify calcium mass was evaluated using calcified arterial vessel phantoms with 20-230 mg of hydroxylapatite. The vessel phantoms were imaged over a Lucite phantom and then an anthropomorphic chest phantom. The total thickness of Lucite phantom ranges from 13.5-26.5 cm to simulate patient thickness of 16-32 cm. The calcium mass was measured using a densitometric technique. The effective dose to patient was estimated from the measured entrance exposure. The effects of patient thickness on contrast-to-noise ratio (CNR), effective dose, and the precision of calcium mass quantification (i.e., the frame to frame variability) were studied. The effects of misregistration artifacts were also measured by shifting the vessel phantoms manually between low- and high-energy images. The results show that, with the same detector signal level, the dynamic filtration technique produced 70% higher calcium contrast-to-noise ratio with only 4% increase in patient dose as compared to the static filtration technique. At the same time, x-ray tube loading increased by 30% with dynamic filtration. The minimum detectability of calcium with anatomical background was measured to be 34 mg of hydroxyapatite. The precision in calcium mass measurement, determined from 16 repeated dual-energy images, ranges from 13 mg to 41 mg when the patient thickness increased from 16 to 32 cm. The CNR was found to decrease with the patient thickness linearly at a rate of (-7%/cm). The anatomic background produced measurement root-mean-square (RMS) errors of 13 mg and 18 mg when the vessel phantoms were imaged over a uniform (over the rib) and nonuniform (across the edge of rib) bone background, respectively. Misregistration artifacts due to motions of up to 1.0 mm between the low- and high-energy images introduce RMS error of less than 4.3 mg, which is much smaller than the frame to frame variability and the measurement error due to anatomic background. The effective dose ranged from 1.1 to 6.6 microSv for each dual-energy image, depending on patient thickness. The study shows that real-time dual-energy imaging can potentially be used as a low dose technique for quantifying coronary arterial calcium.

Automated technique for angiographic determination of coronary blood flow and lumen volume

RATIONALE AND OBJECTIVES

Visual interpretation of angiographic images has been shown to be inadequate for assessing the severity of intermediate coronary stenoses. An approach for evaluating both the anatomic and functional impact of a stenosis is needed. An automated technique for determining both coronary blood flow and lumen volume based on a first-pass analysis (FPA) of coronary angiograms and a template matching algorithm was evaluated.

MATERIALS AND METHODS

Coronary angiograms of a swine animal model were obtained during power injections of contrast material into the left coronary ostium. Background anatomy was subtracted with an automated phase matching program. A template matching algorithm and first-pass analysis were then used to quantify coronary blood flow and lumen volume. Coronary blood flow and lumen volume measurements were validated with a transit-time ultrasound flow probe and a polymer cast of the coronary arteries, respectively.

RESULTS

In 14 independent comparisons, the mean coronary blood flow measured with FPA showed strong correlation with the mean flow measured with the ultrasound flow probe (Q(FPA) = 0.88Q(probe) - 1.99; r = 0.977; standard error of estimate = 3.23 mL/minute). The lumen volumes determined with FPA and cast measurements demonstrated excellent correlation and can be related to each other by V(FPA) = 0.95V(C) - 0.01 (r = 0.997; standard error of estimate = 0.01 mL).

CONCLUSIONS

The proposed automated method for accurate determination of coronary blood flow and lumen volume can supplement visual evaluation of coronary anatomy with quantitative physiologic data. This automated technique potentially offers a clinically feasible method of quantifying coronary blood flow and lumen volume in conjunction with routine cardiac catheterization.

Novel approaches to the measurement of arterial blood flow from dynamic digital X-ray images

We have developed two new algorithms for the measurement of blood flow from dynamic X-ray angiographic images. Both algorithms aim to improve on existing techniques. First, a model-based (MB) algorithm is used to constrain the concentration-distance curve matching approach. Second, a weighted optical flow algorithm (OP) is used to improve on point-based optical flow methods by averaging velocity estimates along a vessel with weighting based on the magnitude of the spatial derivative. The OP algorithm was validated using a computer simulation of pulsatile blood flow. Both the OP and the MB algorithms were validated using a physiological blood flow circuit. Dynamic biplane digital X-ray images were acquired following injection of iodine contrast medium into a variety of simulated arterial vessels. The image data were analyzed using our integrated angiographic analysis software SARA to give blood flow waveforms using the MB and OP algorithms. These waveforms were compared to flow measured using an electromagnetic flow meter (EMF). In total 4935 instantaneous measurements of flow were made and compared to the EMF recordings. It was found that the new algorithms showed low measurement bias and narrow limits of agreement and also out-performed the concentration-distance curve matching algorithm (ORG) and a modification of this algorithm (PA) in all studies.

X-ray videodensitometric methods for blood flow and velocity measurement: a critical review of literature

Blood flow rate and velocity are important parameters for the study of vascular systems, and for the diagnosis, monitoring and evaluation of treatment of cerebro- and cardiovascular disease. For rapid imaging of cerebral and cardiac blood vessels, digital x-ray subtraction angiography has numerous advantages over other modalities. Roentgen-videodensitometric techniques measure blood flow and velocity from changes of contrast material density in x-ray angiograms. Many roentgen-videodensitometric flow measurement methods can also be applied to CT, MR and rotational angiography images. Hence, roentgen-videodensitometric blood flow and velocity measurement from digital x-ray angiograms represents an important research topic. This work contains a critical review and bibliography surveying current and old developments in the field. We present an extensive survey of English-language publications on the subject and a classification of published algorithms. We also present descriptions and critical reviews of these algorithms. The algorithms are reviewed with requirements imposed by neuro- and cardiovascular clinical environments in mind.

Scatter and veiling glare estimation based on sampled primary intensity

Scatter and veiling glare are predominant sources of error in videodensitometric iodine quantification. Standard beam stop techniques such as lead strips or an array of lead discs, placed before the patients, have previously been used to measure scatter and veiling glare in digital radiographic images. However, these techniques significantly increase patient x-ray exposure. In order to overcome this limitation, a scatter measurement technique based on sampled primary intensity has been investigated. This technique uses an array of apertures in a lead sheet to sample the primary x-ray intensity. The scatter-glare intensity in these locations is calculated by subtracting the sampled primary intensity from an open field image which contains both primary and scatter-glare. The calculated scatter-glare values can be interpolated or combined with digital filtration to estimate the scatter-glare intensity on a pixel by pixel basis. The technique was evaluated using a Lucite step phantom and an anthropomorphic chest phantom. The average rms percentage errors of scatter and veiling glare estimation using bi-cubic interpolation and digital filtration techniques were 8.02% and 7.53%, respectively. The average rms percentage errors of primary intensity estimation using bi-cubic interpolation and digital filtration techniques were 10.01% and 8.91%, respectively. The x-ray exposure-area product (EAP) from the aperture array was only 4.38% of the EAP from the open field. These results indicate that the scatter-glare intensity can be accurately estimated with minimal x-ray exposure using sampled primary intensity.

Comparison of methods for instantaneous angiographic blood flow measurement

Several different algorithms have been reported for measurement of blood flow rates and velocities from digital x-ray angiograms. We compare four videodensitometric methods: (1) distance-density curve matching (DDCM), (2) distance-density curve matching with curve-fitting (DDCM-F), (3) bolus mass tracking with curve-fitting (BMT-F) and (4) fluid continuity method (FCM). We tested the flow algorithms with simulated angiograms and with images obtained from a programmable flow phantom under clinically realistic flow and contrast injection conditions including imperfect mixing. All methods perform well for simulated angiograms. On phantom angiograms with constant flow, all methods tended to underestimate flow velocities by at least 7% and demonstrate high variability between consecutive measurements. The FCM demonstrated relatively low variability, but a large negative bias. The DDCM method was moderately biased and had the highest variability. The BMT-F method demonstrated the lowest bias (-7.1%) and the lowest variability both within (27%) and between (27%) studies. No method yields reliable measurements near the peak contrast opacification, when little or no gradient of contrast is present. The extrapolating version of the BMT-F method was also the most robust for estimation of interframe displacements longer than the field of view.

Absolute volumetric coronary blood flow measurement with digital subtraction angiography

The problems associated with visual interpretation of coronary arteriograms have been well-documented. There is a need for more physiologic means of assessing coronary artery stenosis during routine coronary arteriography. Volumetric coronary blood flow assessed as a function of time can be a valuable aid in the analysis of functional significance of arterial obstruction. A volumetric coronary blood flow measurement technique was investigated in a swine animal model using phase matched temporal subtraction images. The left anterior descending (LAD) coronary artery of swine animal models were instrumented with an ultrasound flow probe (US) and a vascular occluder producing stenosis. Contrast material injections (2-4 ml/sec for 3 sec) were made into the left coronary ostium during image acquisition. Phase-matched temporal subtraction (DSA) images were used to measure volumetric coronary blood flow in the LAD. In addition, a technique for automatic estimation of iodine calibration slope was also investigated. In 49 independent comparisons, the mean coronary blood flow (FPA) correlated extremely well with the mean US flow (FPA = 0.92US + 1.42 ml/min, r = 0.99, standard error of estimate (SEE) = 4.32 ml/min). Further more, the automatic iodine calibration technique was shown to be very accurate. In conclusion, accurate volumetric coronary blood flow measurements can be made before the onset of significant changes in coronary blood flow due to contrast injection.

Quantification of volumetric coronary blood flow with dual-energy digital subtraction angiography

BACKGROUND

As a solution to the well-documented problems associated with visual interpretation of coronary arteriograms, more physiological methods of assessing coronary artery stenosis are being investigated. Volumetric coronary blood flow (BF) can be a valuable aid in the analysis of functional significance of arterial obstruction.

METHODS AND RESULTS

The left anterior descending coronary artery (LAD) of 15 anesthetized pigs (40 to 50 kg) was dissected free from the epicardium in its proximal portion, and a transmit-time ultrasound flow probe of the appropriate size was applied. A vascular occluder was positioned distal to the flow probe for flow adjustments. Contrast injections (2 to 4 mL/s for 3 seconds) were made into the left main coronary artery during image acquisition with a motion-immune dual-energy digital subtraction angiography (DE DSA) system. Tissue-suppressed energy-subtracted images were used to generate time-density curves. BF measurements were made in the LAD vascular bed with use of the time-density curve, with consideration that blood was momentarily replaced with contrast during the injection. In 19 comparisons, the mean BF, measured with the use of DE DSA, correlated extremely well with the mean ultrasound flow (DE DSA=0.90 ultrasound+3.10 mL/min, r=.96). Also, contrast injection increased the BF by an average of only 15% during the image-acquisition time interval.

CONCLUSIONS

Accurate BF measurements can be made with motion-immune DE DSA. The BF measurements can be completed before the onset of significant changes in BF due to contrast injection. Furthermore, it is possible to make the BF measurements during routine coronary arteriography.

A digital filtration technique for scatter-glare correction based on thickness estimation

In order to quantitate anatomical and physiological parameters such as vessel dimensions and volumetric blood flow, it is necessary to make corrections for scatter and veiling glare, which are the major sources of nonlinearities in videodensitometric digital subtraction angiography (DSA). A convolution filtering technique has been investigated to estimate scatter-glare distribution in DSA images without the need to sample the scatter-glare intensity for each patient. This technique utilizes exposure parameters and image gray levels to assign equivalent Lucite thickness for every pixel in the image. The thickness information is then used to estimate scatter-glare intensity on a pixel-by-pixel basis. To test its ability to estimate scatter-glare intensity, the correction technique was applied to images of a Lucite step phantom, anthropomorphic chest phantom, head phantom, and animal models at different thicknesses, projections, and beam energies. The root-mean-square (rms) percentage error of these estimates was obtained by comparison with direct scatter-glare measurements made behind a lead strip. The average rms percentage errors in the scatter-glare estimate for the 25 phantom studies and the 17 animal studies were 6.44% and 7.96%, respectively. These results indicate that the scatter-glare intensity can be estimated with adequate accuracy for a wide range of thicknesses, projections, and beam energies using exposure parameters and gray level information.

CCD camera for dual-energy digital subtraction angiography

A motion immune dual-energy subtraction technique in which X-ray tube voltage and beam filtration were switched at 30 Hz between 60 kVp (2.0 mm Al filter) and 120 kVp (2.00 mm Al+2.5 mm Cu filter) was previously reported. In this study the effects of camera lag on the dual-energy iodine signal is investigated. The temporal lag of the lead oxide vidicon tested reduced the dual-energy iodine signal by a factor of 2.3, as compared to a mode that included 4 scrub frames between low- and high-energy images, for an iodine phantom with thicknesses of 0-86.0 mg/cm(2), imaged over a 15 cm thick Lucite phantom. On the other hand, the Charge-Coupled Device (CCD) camera has inherently no temporal lag and its versatile scanning characteristics make it near ideal for dual-energy DSA. The CCD camera eliminates the reduction of dual-energy iodine signal, since it does not mix low- and high-energy image data. Another benefit of the CCD camera is that the separation time between low and high-energy images is not limited to the frame period, as is the lead oxide vidicon; and as small as a 5-msec time difference is possible. The short time interval between low and high-energy images minimizes motion misregistration artifacts. Due to these advantages, the CCD camera significantly improves the utility of dual-energy DSA.