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

A deep-learning-based scatter correction with water equivalent path length map for digital radiography

We proposed a new deep learning (DL) model for accurate scatter correction in digital radiography. The proposed network featured a pixel-wise water equivalent path length (WEPL) map of subjects with diverse sizes and 3D inner structures. The proposed U-Net model comprises two concatenated modules: one for generating a WEPL map and the other for predicting scatter using the WEPL map as auxiliary information. First, 3D CT images were used as numerical phantoms for training and validation, generating observed and scattered images by Monte Carlo simulation, and WEPL maps using Siddon's algorithm. Then, we optimised the model without overfitting. Next, we validated the proposed model's performance by comparing it with other DL models. The proposed model obtained scatter-corrected images with a peak signal-to-noise ratio of 44.24 ± 2.89 dB and a structural similarity index measure of 0.9987 ± 0.0004, which were higher than other DL models. Finally, scatter fractions (SFs) were compared with other DL models using an actual phantom to confirm practicality. Among DL models, the proposed model showed the smallest deviation from measured SF values. Furthermore, using an actual radiograph containing an acrylic object, the contrast-to-noise ratio (CNR) of the proposed model and the anti-scatter grid were compared. The CNR of the images corrected using the proposed model are 16% and 82% higher than those of the raw and grid-applied images, respectively. The advantage of the proposed method is that no actual radiography system is required for collecting training dataset, as the dataset is created from CT images using Monte Carlo simulation.

Determination of contrast-detail curves in mammography image quality assessment by a parametric model observer

A novel approach is proposed for the determination of contrast-detail curves in mammography image quality assessment. The approach is compared with current practice using virtual mammography. A binary parametric model observer is applied to images of the CDMAM phantom. The observer accounts for the simple disc shaped objects in the phantom and is applied separately to each cell of the phantom. For each of these applications, the area under the ROC curve (AUC) of the model observer is determined. The different AUCs, calculated from different applications of the parametric model observer, are then combined to a single contrast-detail curve quantifying the ability of the observer to detect details in the images. Virtual mammography is developed as a tool to simulate X-ray images of single CDMAM cells and to quantitatively assess the approach in comparison with current practice. It is shown that the proposed approach can lead to similar contrast-detail curves as current practice. The precision of the estimated contrast-detail curves is increased, i.e. using 5 images yields about the same precision for the proposed approach as 16 images when applying current practice. We conclude that contrast-detail curves in mammography image quality assessment can also be determined through the AUC of a binary parametric model observer. Since the proposed approach has higher precision than current practice, it is a promising candidate for contrast-detail analysis in mammography image quality assessment.

Detection and quantification of coronary calcium from dual energy chest x-rays: Phantom feasibility study

PURPOSE

We have demonstrated the ability to identify coronary calcium, a reliable biomarker of coronary artery disease, using nongated, 2-shot, dual energy (DE) chest x-ray imaging. Here we will use digital simulations, backed up by measurements, to characterize DE calcium signals and the role of potential confounds such as beam hardening, x-ray scatter, cardiac motion, and pulmonary artery pulsation. For the DE calcium signal, we will consider quantification, as compared to CT calcium score, and visualization.

METHODS

We created stylized and anatomical digital 3D phantoms including heart, lung, coronary calcium, spine, ribs, pulmonary artery, and adipose. We simulated high and low kVp x-ray acquisitions with x-ray spectra, energy dependent attenuation, scatter, ideal detector, and automatic exposure control (AEC). Phantoms allowed us to vary adipose thickness, cardiac motion, etc. We used specialized dual energy coronary calcium (DECC) processing that includes corrections for scatter and beam hardening.

RESULTS

Beam hardening over a wide range of adipose thickness (0-30 cm) reduced the change in intensity of a coronary artery calcification (ΔI_CAC ) by < 3% in DECC images. Scatter correction errors of ±50% affected the calcium signal (ΔI_CAC ) in DECC images ±9%. If a simulated pulmonary artery fills with blood between exposures, it can give rise to a residual signal in DECC images, explaining pulmonary artery visibility in some clinical images. Residual misregistration can be mostly compensated by integrating signals in an enlarged region encompassing registration artifacts. DECC calcium score compared favorably to CT mass and volume scores over a number of phantom perturbations.

CONCLUSION

Simulations indicate that proper DECC processing can faithfully recover coronary calcium signals. Beam hardening, errors in scatter estimation, cardiac motion, calcium residual misregistration etc., are all manageable. Simulations are valuable as we continue to optimize DE coronary calcium image processing and quantitative analysis.

Quantification of fractional flow reserve based on angiographic image data

Coronary angiography provides excellent visualization of coronary arteries, but has limitations in assessing the clinical significance of a coronary stenosis. Fractional flow reserve (FFR) has been shown to be reliable in discerning stenoses responsible for inducible ischemia. The purpose of this study is to validate a technique for FFR quantification using angiographic image data. The study was carried out on 10 anesthetized, closed-chest swine using angioplasty balloon catheters to produce partial occlusion. Angiography based FFR was calculated from an angiographically measured ratio of coronary blood flow to arterial lumen volume. Pressure-based FFR was measured from a ratio of distal coronary pressure to aortic pressure. Pressure-wire measurements of FFR (FFR_P) correlated linearly with angiographic volume-derived measurements of FFR (FFR_V) according to the equation: FFR_P = 0.41 FFR_V + 0.52 (P-value < 0.001). The correlation coefficient and standard error of estimate were 0.85 and 0.07, respectively. This is the first study to provide an angiographic method to quantify FFR in swine. Angiographic FFR can potentially provide an assessment of the physiological severity of a coronary stenosis during routine diagnostic cardiac catheterization without a need to cross a stenosis with a pressure-wire.

Estimation of coronary artery hyperemic blood flow based on arterial lumen volume using angiographic images

The purpose of this study is to develop a method to estimate the hyperemic blood flow in a coronary artery using the sum of the distal lumen volumes in a swine animal model. The limitations of visually assessing coronary artery disease are well known. These limitations are particularly important in intermediate coronary lesions where it is difficult to determine whether a particular lesion is the cause of ischemia. Therefore, a functional measure of stenosis severity is needed using angiographic image data. Coronary arteriography was performed in 10 swine (Yorkshire, 25-35 kg) after power injection of contrast material into the left main coronary artery. A densitometry technique was used to quantify regional flow and lumen volume in vivo after inducing hyperemia. Additionally, 3 swine hearts were casted and imaged post-mortem using cone-beam CT to obtain the lumen volume and the arterial length of corresponding coronary arteries. Using densitometry, the results showed that the stem hyperemic flow (Q) and the associated crown lumen volume (V) were related by Q = 159.08 V(3/4) (r = 0.98, SEE = 10.59 ml/min). The stem hyperemic flow and the associated crown length (L) using cone-beam CT were related by Q = 2.89 L (r = 0.99, SEE = 8.72 ml/min). These results indicate that measured arterial branch lengths or lumen volumes can potentially be used to predict the expected hyperemic flow in an arterial tree. This, in conjunction with measured hyperemic flow in the presence of a stenosis, could be used to predict fractional flow reserve based entirely on angiographic data.

Quantification of breast density with dual energy mammography: a simulation study

Breast density, the percentage of glandular breast tissue, has been identified as an important yet underutilized risk factor in the development of breast cancer. A quantitative method to measure breast density with dual energy imaging was investigated using a computer simulation model. Two configurations to measure breast density were evaluated: the usage of monoenergetic beams and an ideal detector, and the usage of polyenergetic beams with spectra from a tungsten anode x-ray tube with a detector modeled after a digital mammography system. The simulation model calculated the mean glandular dose necessary to quantify the variability of breast density to within 1/3%. The breast was modeled as a semicircle 10 cm in radius with equal homogenous thicknesses of adipose and glandular tissues. Breast thicknesses were considered in the range of 2-10 cm and energies in the range of 10-150 keV for the two monoenergetic beams, and 20-150 kVp for spectra with a tungsten anode x-ray tube. For a 4.2 cm breast thickness, the required mean glandular doses were 0.183 microGy for two monoenergetic beams at 19 and 71 keV, and 9.85 microGy for two polyenergetic spectra from a tungsten anode at 32 and 96 kVp with beam filtrations of 50 microm Rh and 300 microm Cu for the low and high energy beams, respectively. The results suggest that for either configuration, breast density can be precisely measured with dual energy imaging requiring only a small amount of additional dose to the breast. The possibility of using a standard screening mammogram as the low energy image is also discussed.

Quantitative coronary angiography using image recovery techniques for background estimation in unsubtracted images

Densitometry measurements have been performed previously using subtracted images. However, digital subtraction angiography (DSA) in coronary angiography is highly susceptible to misregistration artifacts due to the temporal separation of background and target images. Misregistration artifacts due to respiration and patient motion occur frequently, and organ motion is unavoidable. Quantitative densitometric techniques would be more clinically feasible if they could be implemented using unsubtracted images. The goal of this study is to evaluate image recovery techniques for densitometry measurements using unsubtracted images. A humanoid phantom and eight swine (25-35 kg) were used to evaluate the accuracy and precision of the following image recovery techniques: Local averaging (LA), morphological filtering (MF), linear interpolation (LI), and curvature-driven diffusion image inpainting (CDD). Images of iodinated vessel phantoms placed over the heart of the humanoid phantom or swine were acquired. In addition, coronary angiograms were obtained after power injections of a nonionic iodinated contrast solution in an in vivo swine study. Background signals were estimated and removed with LA, MF, LI, and CDD. Iodine masses in the vessel phantoms were quantified and compared to known amounts. Moreover, the total iodine in left anterior descending arteries was measured and compared with DSA measurements. In the humanoid phantom study, the average root mean square errors associated with quantifying iodine mass using LA and MF were approximately 6% and 9%, respectively. The corresponding average root mean square errors associated with quantifying iodine mass using LI and CDD were both approximately 3%. In the in vivo swine study, the root mean square errors associated with quantifying iodine in the vessel phantoms with LA and MF were approximately 5% and 12%, respectively. The corresponding average root mean square errors using LI and CDD were both 3%. The standard deviations in the differences between measured iodine mass in left anterior descending arteries using DSA and LA, MF, LI, or CDD were calculated. The standard deviations in the DSA-LA and DSA-MF differences (both approximately 21 mg) were approximately a factor of 3 greater than that of the DSA-LI and DSA-CDD differences (both approximately 7 mg). Local averaging and morphological filtering were considered inadequate for use in quantitative densitometry. Linear interpolation and curvature-driven diffusion image inpainting were found to be effective techniques for use with densitometry in quantifying iodine mass in vitro and in vivo. They can be used with unsubtracted images to estimate background anatomical signals and obtain accurate densitometry results. The high level of accuracy and precision in quantification associated with using LI and CDD suggests the potential of these techniques in applications where background mask images are difficult to obtain, such as lumen volume and blood flow quantification using coronary arteriography.

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.

Effect of vessel orientation on videodensitometry quantitative coronary arteriography

Cross-sectional area measurement using videodensitometry coronary arteriography is a function of the out-of-plane angle of the vessel segment with respect to the face of the image intensifier. Therefore, absolute cross-sectional area measurement using videodensitometry requires correction for the out-of-plane angle. In this study, a simple technique to measure the out-of-plane angle of a vessel segment using images from two different projections is presented. The technique was tested using vessel phantoms imaged at different out-of-plane angles. A vessel segment with an out-of-plane angle of less than 20 degrees introduces less than 5% error in the cross-sectional area measurements. However, an out-of-plane angle of 60 degrees can introduce more than 90% error. The results of the measurements show that correction for the out-of-plane angle reduces the error in cross-sectional area measurement to less that 5%. The correction for the out-of-plane angle makes it possible to measure the absolute cross-sectional area using videodensitometery, which can be used to assess a diseased vessel segment with any complex geometry.

Regional volumetric coronary blood flow measurement by digital angiography: in vivo validation

RATIONALE AND OBJECTIVES

There are well-known limitations to the use of visual estimation to assess the severity of coronary artery disease and luminal stenosis. This is especially true in the case of an intermediate coronary lesion (30%-70% diameter stenosis), where coronary arteriography is very limited in distinguishing ischemia-producing intermediate coronary lesions from non-ischemia-producing ones. For this reason, a functional measure of stenosis severity is desirable. The goal of this study is to validate a video densitometry technique for quantitative assessment of regional volumetric coronary blood flow.

MATERIALS AND METHODS

Coronary arteriography was performed in eight swine (body weight, 25-50 kg) after power injection of contrast material into the left main coronary artery. Phase-matched subtracted images were used to quantify regional coronary blood flow using a video densitometry technique. The in vivo regional flow measurements were validated using a transit time ultrasound flow probe.

RESULTS

In 44 measurements, the ultrasound (Q(US)) and video densitometry (Q(VD)) regional flow measurements were related by Q(VD) = 0.98 Q(US) + 0.11 mL/min (r = 0.98). The results of mean regional coronary blood flow measurements for repeated coronary arteriograms of the first (Q(VD1)) and second (Q(VD2)) measured flows were related by Q(VD1) = 1.04 Q(VD2) + 0.05 mL/min (r = 0.97).

CONCLUSIONS

A video densitometry technique for quantification of regional coronary blood flow was validated using a swine animal model. The results demonstrated the feasibility and potential utility of the video densitometry technique for accurate measurement of regional coronary blood flow, in vivo. This study provides an angiographic method that can potentially be used to evaluate intermediate coronary lesions during routine coronary arteriography.

Assessment of vasoreactivity using videodensitometry coronary angiography

BACKGROUND

Previous studies demonstrated that the dysfunction of vasomotor tone (VT) is closely linked to the development of atherosclerosis and it is considered important in the very early stages of atherogenesis. Currently, the evaluation of VT relies on lumen changes in response to vasoactive stimuli using quantitative coronary angiography (QCA) based on geometric edge detection (ED). However, using ED for measuring lumen diameters is inherently associated with large uncertainties. Videodensitometry (VD) methods have important advantages over ED for QCA. The objective of this study was to investigate the reliability of VD and ED techniques in determining the effect of nitroglycerin (NTG) on cross-sectional area (CSA) and volume changes in a swine animal model for evaluating coronary vasoreactivity.

METHODS AND RESULTS

Coronary angiography was performed on four anesthetized swine. CSA and volume were measured in the left anterior descending (LAD) coronary artery using VD before and after intracoronary injection of 0.3 mg of NTG. CSA was also calculated using standard QCA based on ED. The average CSA changes in the proximal, middle and distal branches measured using VD were 23.83% (+/-10.76%), 30.78% (+/-18.39%), and 27.34% (+/-36.53%), respectively. Similarly, the average CSA changes in the proximal, middle, and distal branches measured using ED were 15.02% (+/-36.38%), 22.02% (+/-26.12), and 38.00% (+/-48.31%), respectively. The average lumen volume change measured using VD was 29.79% (+/-14.79%). In order to evaluate the relative reliability of the techniques. the significance of deviation (SOD) was calculated, which is the ratio of the change after NTG and the measurement error. The average SOD for CSA for all the branches based on VD and ED were 1.86 and 0.69, respectively. The SOD for volume measurement was 2.78.

CONCLUSIONS

Lumen changes measured by VD showed substantial improvement in reliability when compared to the ED. Moreover, VD can be used to measure substantially smaller changes in lumen dimension in response to vasoactive stimuli than the standard QCA based on ED. Finally, VD allows the measurement of arterial volume, which is not possible with ED.

Effect of passive myocardium on the compliance of porcine coronary arteries

The objective of this study was to determine the effect of passive myocardium on the coronary arteries under distension and compression. To simulate distension and compression, we placed a diastolic-arrested heart in a Lucite box, where both the intravascular pressure and external (box) pressure were varied independently and expressed as a pressure difference (DeltaP = intravascular pressure - box pressure). The DeltaP-cross-sectional area relationship of the first several generations of porcine coronary arteries and the DeltaP-volume relationship of the coronary arterial tree (vessels >0.5 mm in diameter) were determined using a video densitometric technique in the range of +150 to -150 mmHg. The vasodilated left anterior descending (LAD) coronary artery of six KCl-arrested hearts were perfused with iodine and 3% Cab-O-Sil. The intravascular pressure was varied in a triangular pattern, whereas the absolute cross-sectional area of each vessel and the total arterial volume were calculated using video densitometry under different box pressures (0, 50, 100, and 150 mmHg). In the range of positive DeltaP, we found that the compliance of the proximal LAD artery in situ (4.85 +/- 3.8 x 10-3 mm2/mmHg) is smaller than that of the same artery in vitro (16.5 +/- 6 x 10-3 mm2/mmHg; P = 0.009). Hence, the myocardium restricts the compliance of the epicardial artery under distension. In the negative DeltaP range, the LAD artery does not collapse, whereas the same vessel readily collapses when tested in vitro. Hence, we conclude that myocardial tethering prevents collapse of large blood vessel under compression.

X-ray beam equalization: feasibility and performance of an automated prototype system in a phantom and swine

PURPOSE

To investigate the feasibility of an automated implementation of a beam equalization technique and to evaluate the experimental performance of the prototype system.

MATERIALS AND METHODS

X-ray beam equalization involved the process of low-dose image acquisition, attenuator thickness calculation, mask generation, mask positioning, equalized image acquisition, and mask reshaping. The entire equalization process was performed in approximately 7 seconds. The equalized images were assessed both qualitatively and quantitatively by using a humanoid phantom and a swine animal model. The general image quality was assessed for the ability to visualize arterial branches and other anatomic structures. The level of equalization was quantitatively assessed by segmenting the images into an 8 x 8 matrix of square regions.

RESULTS

The ratio of the root-mean-squared variance for the equalized and unequalized images from humanoid phantom and swine animal studies was 0.49 and 0.59, respectively. Furthermore, qualitative assessment of the images showed substantial improvement in image quality and visualization of arterial branches after beam equalization in phantom and animal studies.

CONCLUSION

Automated area beam equalization is feasible and improves image quality in previously underpenetrated regions of the image.

Quantification of coronary artery lumen volume by digital angiography: in vivo validation

BACKGROUND

Coronary artery lumen volume may potentially have several advantages over the commonly used variables, such as percent stenosis or minimal lumen diameter, in the assessment of coronary artery disease. The goal of this study is to validate a quantitative assessment of lumen volume using a video densitometry technique.

METHODS AND RESULTS

Coronary arteriography was performed in 9 swine (body weight 20 to 55 kg) after power injection of contrast material (2 mL/s for 3 seconds) into the left main coronary artery. Phase-matched subtracted images were used to quantify regional lumen volume by a video densitometry technique. The in vivo volume measurements were validated by a polymer cast of the coronary arterial tree made at physiological pressure. The measured cast volume (V(C)) and video densitometric regional lumen volume (V(VD)) were related by V(VD)=1.06 V(C)-0.01 mL (r=0.99). The root mean square and systematic errors for these measurements were 17% and -3%, respectively.

CONCLUSIONS

A video densitometry technique for quantification of coronary lumen volume was validated both in vitro and in vivo in a swine animal model. The present results demonstrated the feasibility and potential utility of the video densitometry technique for accurate measurement of regional lumen volume in vivo. This study contributes to the understanding of the angiographic methods used for the assessment of coronary artery disease and indicates that this technique can potentially be used for quantification of diffuse coronary artery disease during routine coronary arteriography.

Cross-sectional area and volume compliance of porcine left coronary arteries

We have determined the cross-sectional area (CSA) compliance of the first several generations of pig coronary arteries and the volume compliance of the coronary arterial tree (vessels >0.5 mm in diameter) using a videodensitometric technique. The coronary arteries of four KCl-arrested maximally vasodilated pig hearts were perfused with iodine and 3% Cab-O-Sil. Because Cab-O-Sil occludes small arteries, the flow can be stopped and the pressure can be maintained while the trunk of the coronary artery and its subbranches are imaged using digital angiography. The coronary arteries were preconditioned several times with cyclic changes in pressure from 0 to 160 mmHg. The pressure was then varied in a triangular pattern, and the absolute CSA of each vessel and the total arterial volume were calculated using videodensitometry in conjunction with digital subtraction angiography. Our results have shown that the pressure-diameter and pressure-volume relationships are linear in the 60-140 mmHg pressure range. Furthermore, the compliance of the coronary arteries is small; i.e., the diameter of the coronary artery changes by <15% in the 80-mmHg pressure range. The compliance data couples the mechanics of the blood vessel wall to the mechanics of blood flow to yield a pressure-flow relationship for each coronary arterial segment.

Area x-ray beam equalization for digital angiography

An area beam equalization technique has been investigated in order to generate patient-specific compensating filters for digital angiography. An initial image was used to generate the compensating filter, which was fabricated using a deformable compensating material, containing CeO2, and an array of square pistons. The CeO2 attenuator thicknesses were calculated using the gray level information from the initial unequalized image. The array of pistons was pressed against a uniform thickness of attenuating material to generate a filter for x-ray beam equalization. The filter was subsequently inserted into the x-ray beam for the final equalized radiograph. It was positioned close to the focal spot (magnification of 8.0) in order to minimize edge artifacts from the filter. The equalization of x-ray transmission across the field exiting from the object significantly improved the image quality by preserving local contrast throughout the image. The contrast-to-noise ratio (CNR) in the equalized images was increased-by up to fivefold. Phantom studies indicate that equalized images using a relatively small array of pistons (e.g., 8 x 8) produce significant improvement in image quality with negligible perceptible artifacts. Animal studies showed that beam equalization significantly improved fluoroscopic and angiographic image quality. X-ray beam equalization produced an image with a relatively uniform scatter-glare intensity and it reduced the scatter-glare fraction in the previously underpenetrated region of the image from 0.65 to 0.50. Also, x-ray tube loading due to the mask assembly itself was negligible. In conclusion, area beam equalization reduces the scatter-glare fraction and significantly improves CNR in the previously underpenetrated region of the image.

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.

Scatter-glare estimation for digital radiographic systems: comparison of digital filtration and sampling techniques

The scatter and veiling glare distribution in images acquired with a digital subtraction angiography imaging system was estimated using a digital filtration and a beam-stop technique. The digital filtration technique utilizes exposure parameters and image gray levels to estimate scatter-glare intensity based on previous phantom measurements. The beam-stop technique uses an array of lead discs in order to sample scatter-glare intensity for each patient. To test the ability of digital filtration and beam-stop techniques to estimate the scatter-glare intensity, they were applied to images of postmortem swine animal models at different projections and beam energies. The systematic and root-mean-square (rms) percentage errors of these estimates were obtained by comparison to directly measured scatter-glare images using a scanning lead strip technique. The average rms percentage error for the digital filtration and beam-stop techniques were 8.07% and 6.67%, respectively. The changes in scatter-glare intensity due to contrast injection during coronary arteriography and ventriculography were also measured using the beam-stop technique. The maximum changes in scatter-glare intensities during coronary arteriography and ventriculography were 19 and 88%, respectively. The results indicate that the digital filtration technique is more suited for applications such as coronary arteriography and ventriculography where the iodinated contrast material significantly changes the scatter-glare intensity.

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.

Measurement of a cross-sectional area of normal and stenotic arteries with videodensitometric quantitative arteriography and intravascular ultrasound

RATIONALE AND OBJECTIVES

A videodensitometric technique that allows measurement of absolute cross-sectional area of any complex lesion was compared with an intravascular ultrasound (US) technique.

METHODS

Stenotic devices (10-15 mm long) with cross sections of different shapes were placed in the distal aortas of five anesthetized pigs (weight, 40-50 kg). The stenotic devices were imaged by using an intravascular US probe after power injection of contrast material.

RESULTS

A comparison of actual areas and measured cross-sectional areas of the stenotic devices showed that videodensitometry and intravascular US produced better results than edge-detection techniques for both unsubtracted and temporal subtraction images.

CONCLUSION

These data suggest that the videodensitometric technique can be used to measure absolute cross-sectional areas of arteries with different shapes.

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.