Search

Purpose: Biplane angiography systems provide time resolved 2D fluoroscopic images from two different angles, which can be used for the positioning of interventional devices such as guidewires and catheters. The purpose of this work is to provide a novel algorithm framework, which allows the 3D reconstruction of these curvilinear devices from the 2D projection images for each time frame. This would allow creating virtual projection images from arbitrary view angles without changing the position of the gantries, as well as virtual endoscopic 3D renderings., Methods: The first frame of each time sequence is registered to and subtracted from the following frame using an elastic grid registration technique. The images are then preprocessed by a noise reduction algorithm using directional adaptive filter kernels and a ridgeness filter that emphasizes curvilinear structures. A threshold based segmentation of the device is then performed, followed by a flux driven topology preserving thinning algorithm to extract the segments of the device centerline. The exact device path is determined using Dijkstra's algorithm to minimize the curvature and distance between adjacent segments as well as the difference to the device path of the previous frame. The 3D device centerline is then reconstructed using epipolar geometry., Results: The accuracy of the reconstruction was measured in a vascular head phantom as well as in a cadaver head and a canine study. The device reconstructions are compared to rotational 3D acquisitions. In the phantom experiments, an average device tip accuracy of 0.35 ± 0.09 mm, a Hausdorff distance of 0.65 ± 0.32 mm, and a mean device distance of 0.54 ± 0.33 mm were achieved. In the cadaver head and canine experiments, the device tip was reconstructed with an average accuracy of 0.26 ± 0.20 mm, a Hausdorff distance of 0.62 ± 0.08 mm, and a mean device distance of 0.41 ± 0.08 mm. Additionally, retrospective reconstruction results of real patient data are presented., Conclusions: The presented algorithm is a novel approach for the time resolved 3D reconstruction of interventional devices from biplane fluoroscopic images, thus allowing the creation of virtual projection images from arbitrary view angles as well as virtual endoscopic 3D renderings. Availability of this technique would enhance the ability to accurately position devices in minimally invasive endovascular procedures.

Purpose: Recent efforts in the reconstruction of interventional devices from two distinct views require the segmentation of the object in both fluoroscopic images. Noise might decrease the quality of the segmentation and cause artifacts in the reconstruction. The noise level depends on the x-ray dose the patient is exposed to. The proposed algorithm reduces the noise and enhances the separability of curvilinear devices in background subtracted fluoroscopic images to allow a more accurate segmentation., Methods: The algorithm uses a set of binary masks to estimate a line conformity measure that determines the best direction for a directional filter kernel. If the calculated value exceeds a certain threshold, the directional kernel is used to obtain the filtered value. Otherwise, an isotropic filter kernel is used., Results: The evaluation was performed on a set of 36 fluoroscopic images using a vascular head phantom with three different guidewires and nine different x-ray dosages from 6 nGy/pulse to 45 nGy/pulse as well as a clinical data set containing ten images. Compared with wavelet shrinkage and the bilateral filter, the proposed algorithm increased the average contrast to noise ratio by at least 17.8% for the phantom and 68.9% for the clinical images. The accuracy of the device segmentation was improved on average by at least 17.3% and 14.0%, respectively., Conclusions: The proposed algorithm was able to significantly reduce the amount of noise in the images and therefore increase the quality of the device segmentations compared to both the bilateral filter and the wavelet thresholding approach for all acquired noise levels using rotating directional filter kernels near line structures and isotropic kernels for the background. The application of the proposed algorithm for the 3D reconstruction of curvilinear devices from two views would allow a more accurate reconstruction of the device.

During the past decade the use of undersampled acquisition and constrained reconstruction have led to significant increases in data acquisition speed, SNR, spatial resolution and temporal resolution in MR imaging. When a separately acquired constraining image is combined with an angiographic time series the traditional tradeoff between spatial and temporal resolution is greatly reduced. Artifacts and limited resolution that would normally be associated with a rapid highly undersampled temporal image series are mitigated by the constrained reconstruction process which transfers the SNR and spatial resolution of the constraining image to the individual time frames. In rotational C-Arm DSA a 3D image volume is formed from all the projections acquired during the C-Arm rotation. Although the individual projections contain temporal information, the reconstructed 3D image has no temporal information and represents a composite of the vascular filling that has occurred during the iodine injection. However, the 3D cone beam CT reconstruction can be used to constrain the reconstruction of one 3D volume for each of the rotational projections. This extends the traditional DSA time series of 2D images to a series of 3D volumes at rates up to 30 per second. Similar techniques can be used to provide fluoroscopy that can be embedded in the 3D space of the constraining volume and viewed from arbitrary angles without gantry motion. This overcomes the problem of forbidden views and guarantees that an intervention can be done without having to send patients to surgery. Unlike 4D DSA which requires only one source and receptor, 4D Fluoroscopy requires a bi-plane fluoroscopy system., Learning Objectives: 1. To understand the application of under sampling and constrained reconstruction to 4D DSA and Fluoroscopy., (© 2012 American Association of Physicists in Medicine.)

Computed radiography (CR) has shown promise in digital mammographic screening due to its good low spatial frequency MTF and its relatively wide exposure latitude. The CR image format has not gained acceptance clinically because of reduced high spatial frequency resolution as compared to film-screen images. X-ray capillary optics, aligned between the breast and CR phosphor imaging plate, will capture primary x-ray photons almost exclusively. Due to the very small angle of acceptance, scattered photons angled more than about 1.6 x 10(-3) radians from primary trajectory will not be accepted at the capillary optic entrance. The virtual elimination of detected scatter means almost 100% of the possible primary contrast should be visible in the image. In addition, the image can be magnified without focal spot blurring. Effective resolution of CR images can be increased by a factor equal to that magnification. Clinical implementation of future capillary optics are expected to be either in the form of a large, stationary, post-patient optic that accepts primary from the entire breast or a fan-shaped optic that is scanned across the breast. Measurements of a test capillary optic showed a reduction of scatter fraction to 0.018. Images of a lucite contrast detail phantom revealed a corresponding increase in image contrast when compared to anti-scatter grid and no grid methods. Spectral transmission measurements using a high-purity germanium detector showed good primary transmission (45%-50%) in the mammographic energy range. The MTF measurements of both stationary and scanned capillary optics showed improvement at the 5% MTF level to 8.4 mm-1 for scanned optics and 9.2 mm-1 for stationary optics representing a 68% and 84% respective increase over the CR MTF without magnification or capillary optics.

Capillary optic arrays are bundles of hollow glass capillaries which guide x rays in a manner similar to the way fiber optics guide light. Focused postpatient capillary optic arrays have the potential to significantly improve both contrast and resolution of mammographic images compared to conventional antiscatter grids. Contrast can be improved by the nearly total scatter rejection of the optic. Effective resolution can be improved by geometric magnification without increased focal spot blurring. The best results were found for borosilicate glasses, with transmissions in excess of 60% for 22-cm-long fibers. To evaluate the scatter rejection properties, the transmission of off-axis radiation was measured. Transmission drops to < 1% at an angular displacement of 2.7 mrad. Transmission of a bulk capillary array dropped to near zero if the source was at an angle of 2.5 mrad. This implies excellent scatter rejection capabilities. To evaluate whether unchanneled photons might still reach the detector, absorption measurements were also performed on fibers and arrays. Absorption was found to be adequate for scatter rejection. All of the data agreed well with numerical simulations. Performance calculations for two potential optics geometries gave promising results.

Single kernel scatter correction algorithms are based on the model that the scatter field can be predicted by convolution of the primary intensity (Iprim) with a spatially invariant scatter point-spread function (PSF). Practical limitations (Iprim unknown) suggest the substitution of the total detected intensity (Idet) for Iprim as the source image in the convolution. In regions of high scatter fraction (SF), Idet is a poor approximation of Iprim, thereby causing an overestimation of scatter originating in the region. This contributes to errors in estimating detected scatter in the mediastinum and neighboring regions. A technique using a regionally variable point-spread function that significantly reduces RMS error in estimation of the primary image as compared to the single PSF method is investigated. The regionally variable convolution method employs a larger PSF in the mediastinum and a smaller PSF in the lungs to reduce the error in estimating the scatter throughout the image. The method to allow for patient differences has also been expanded and various implementations of these methods have been compared. Results show that the dual-kernel algorithm is always more effective than an equivalent single-kernel algorithm. The dual-kernel algorithm using a predicted scatter fraction curve gives an overall RMS error in the primary of as low as 20.8% which is equivalent to 8.7% RMS error in the scatter. The dual-kernel method using a predicted scatter fraction curve approaches the accuracy of the single-kernel method using patient specific scatter measurements. Because using individual scatter measurements is a less desirable method for clinical use, we feel that the dual-kernel algorithm which uses two regions specific convolution kernels and a variable scatter fraction curve is the preferable method.

The generation of monoenergetic or quasimonoenergetic x-ray spectra has been accomplished by several methods including the use of K-edge filtration, characteristic radiation, crystal monochromators, and multilayer mirrors. In this paper some practical questions are discussed regarding the possibility of obtaining narrow-band spectra using x-ray reflection in glass capillary optics which have been reported recently in connection with the focusing and generation of parallel x-ray beams. Derivation of formulas for double and triple reflection with and without additional K-edge filtration imply that these methods are superior to the use of filtration alone. It is shown that the double reflection scheme is most sensitive to the angular divergence of the beam emanating from the capillary array used to generate the incident parallel beam. Simulations using three reflections predict output spectra which are relatively insensitive to blurring due to angular divergence. A small amount of K-edge filtration in combination with three reflections provides considerable sharpening of the x-ray spectrum. Aside from the spectral selectivity provided by the capillary array, the reduced divergence of the primary transmission may have advantages for scatter rejection in situations where the air gap between the patient and detector is increased by an amount consistent with the resolution requirements of the application. In parallel beam situations, the use of capillary arrays between the patient and detector may also be useful for the rejection of divergent scattered radiation.

The use of characteristic-line radiation from rare-earth targets bombarded by high-energy (up to 1 MeV) electron beams has been evaluated as an x-ray source for dual energy K-edge subtraction imaging of the human coronary arteries. Two characteristic-line x-ray sources, one using the split K alpha 1 and K alpha 2 lines of lanthanum excited by a high-energy electron beam and the other using the K alpha lines of barium and cerium, were studied. A Monte Carlo electron-photon simulation was used to calculate x-ray spectra and energy deposition profiles from targets of these elements bombarded by electrons in the energy range 140 keV to 1 MeV. A general dual-energy imaging model was developed that used these calculated source spectra to numerically investigate the dependence of the subtraction image signal-to-noise ratio on such factors as the ratio of K-line to x-ray continuum yield, continuum spectral shape, x-ray filtering, and detector response. A signal averaging technique for enhancing the signal-to-noise ratio was also evaluated. The results of these calculations were used to identify an optimum electron beam, target, filter, and detector configuration. A compact electron accelerator capable of providing the required electron beam parameters was designed. Calculations indicate that under ideal conditions the optimized system would be capable of imaging 2 mg/cm2 of iodine contrast agent in 20 g/cm2 of tissue with a signal-to-noise ratio of 5, a detector pixel size of 0.25 mm2, and a total image acquisition time of 10 ms. These parameters are consistent with those needed to image the human coronary arteries after an intravenous injection of iodine contrast agent. These capabilities, along with the relatively modest hardware requirements of this system, make it attractive as an x-ray source for dual energy transvenous coronary angiography.

Due to the well-documented problems associated with visual interpretation of coronary angiograms, more physiologic means of assessing coronary artery stenosis are being investigated. One physiologic parameter that has been suggested is coronary flow reserve (CFR). A digital subtraction angiographic technique based on first pass distribution analysis (FPA) is proposed as a means of measuring CFR and absolute coronary flow. The theory of the FPA method is first outlined, and the implementation of a preliminary version of the FPA algorithm is described. Experiments verifying the utility of this algorithm for measuring absolute flow through a flow phantom, and through the canine circumflex artery are reported. It was determined that the preliminary FPA algorithm is capable of measuring canine coronary flow ratios (R) with accuracy and precision characteristics meeting or exceeding those reported for the parametric imaging technique (RFPA = 0.933.Rtrue, SEE = 0.16, r = 0.984). Accurate absolute flow (Q) measurements were obtained in all of the phantom experiments (QFPA = 1.054.Qtrue, r = 0.993), and in one of the three dogs that were studied (QFPA = 0.977.Qtrue, r = 0.935). The difficulty encountered in the other two dog experiments is attributed to the effects of system temporal lag, and would likely be corrected through the use of improved cameras. The feasibility of the general FPA method for measuring relative flow is established, and the potential for routine, absolute flow measurement is demonstrated.

A digital beam attenuator system has been developed to automatically generate patient-specific compensating filters for chest radiography. An initial low-dose test image is used to generate the attenuator, which is fabricated by overprinting multiple layers of a heavy-metal material onto a nonattenuating substrate. The attenuator is subsequently inserted into the x-ray beam for a final compensated radiograph. The effects of focal spot blurring and limited attenuator resolution result in the final compensated image containing only high-spatial frequency information. The frequency response of the process is not strictly describable by a modulation transfer function, but an approximation of the frequencies remaining in the compensated image is obtained for low-contrast conditions. It is found that a 4 X 4 blurring function on the original 64 X 64 test image is required for the attenuator to give appropriate compensated image appearance. A proposed attenuator printing scheme prints the attenuator in a 16 X 16 matrix, staggering successively printed layers to achieve the required 64 X 64 sampling with appropriate blurring. The resulting compensated image has good anatomical definition and contains a frequency response similar to that obtained by compensation techniques being investigated by Plewes and Sorenson.

Our lab has previously generated selective iodine images with an image intensifier fluoroscopic system using a three-beam K-edge approach. Logarithmically amplified video images Li were linearly combined to yield the final image k1L1 + k2L2 + k3L3. This paper discusses refinements of the K-edge technique. A study is made of the manner in which contrast-reducing effects such as x-ray scatter and image intensifier veiling glare enter into the final image. If such biases can be approximated as multiplicative and independent of the x-ray spectrum, and if the sum of the ki is zero, then the biases are canceled. Experimental data is presented which demonstrates that the inaccuracy due to such biases can be reduced by a factor as large as 10. The theorem that K1 + K2 + K3 is approximately equal to 0 is proven rigorously and discussed. Because the ki add to zero, the final image can be expressed as a linear combination of two of the differences between the Li. A difference-based algorithm which reduces biases and make allowance for nonlinearities such as beam hardening is proposed and discussed.

In dual-beam selective iodine imaging, images of an object are made with each of two spectrally different x-ray beams. The mean beam energies may either straddle the 33 keV iodine K-edge or both lie above the K-edge. Both patient exposure considerations and the availability of sufficient x-ray flux make the latter approach favorable for tissue thicknesses exceeding 5 cm. Consider such an approach in which image contrast from tissue is suppressed in the difference image. It is proven theoretically that the residual bone-to-iodine contrast is a constant independent of the two mean beam energies used. This invariance principle is demonstrated experimentally by comparing images made from different pairs of x-ray spectra. Observed contrast ratios match the predicted value very well. In dual-beam imaging, contrast from only one material may be suppressed. Other substances yield residual signals which compete with the iodine. Subtleties of this incomplete cancellation are demonstrated, discussed, and quantitated. A contrast enhancement factor (CEF) is defined as the factor by which iodine contrast is enhanced in a multiple beam subtraction technique relative to monoenergetic imaging at 40 keV. CEFs are determined for tissue and bone cancellation separately and their limits are discussed. Images of a simulated artery containing iodine superimposed over a Rando head and neck phantom show that the CEF limitation for dual beam imaging is quite severe compared to a time dependent mask mode imaging approach. Finally, optimum, energies for dual beam images are discussed.

We have previously reported a motion immune dual-energy subtraction technique in which x-ray tube voltage and x-ray beam filtration are switched at 30 Hz between 60 kVp (2.0-mm Al filter) and 120 kVp (2.0-mm Al + 2.5-mm Cu filtration). In this paper we consider the suitability of these dual-energy images for quantitative measurements of iodine thickness and volume. Optimized iodine signal-to-noise ratio (S/N) was measured as a function of phantom thickness. Using a fixed mAs, the S/N of the dual-energy images was found to decrease by sevenfold as lucite thickness increased from 10 to 25 cm. For the same increase in lucite thickness S/N for time subtraction images decreased by fivefold. Image quality in two human volunteers was subjectively judged to be good. In order to quantitate physiological parameters such as ejection fraction and left ventricular volume, energy dependent corrections for scatter and veiling glare, beam hardening, detector nonuniformity, heel effect, and uncanceled bone signals were developed. Since the dual-energy technique does not completely cancel bone, a preinjection dual-energy subtraction image was used to estimate integrated bone contributions to iodine volume measurements. In a phantom measurement simulating exercise ventriculography, the known (Vk) and videodensitometrically measured (Vm) volumes of 19 mg/cm3 solution of iodine were related by Vm = 0.95 Vk + 1.50 cm3 (r greater than 0.99).

The properties of tomography, K-edge imaging, and K-edge tomography are discussed in relation to the imaging of small concentrations of elements such as iodine and xenon and are compared by means of phantom images. It is demonstrated that the complementary selectivities provided by depth and energy subtraction are combined in K-edge tomography. Using a three-spectrum subtraction technique, the iodine difference signal predicted by computer calculations is on the order of 8000 times that of an equal concentration fo bone. The corresponding ratio in tomography without energy subtraction is 20:1. It is argued that K-edge tomography can successfully eliminate artifacts due to tissue inhomogeneities which presently enable 0.6% variations in tissue attenuation to mimic minimum detectable iodine signals in conventional computed tomography. Various instrumentation possibilities and energy subtraction techniques are discussed.

In our previous reports on absorption-edge fluoroscopy, it was not possible to relate fully the subtleties involved in the selection of spectral parameters. This paper is intended as an overview of this important aspect of the technique. It is shown that, by using the 1-kVp, 2-filter technique, it is possible to image certain elements (e.g., iodine and xenon) in the presence of tissue variations of +/-2 cm about the thickness at which perfect tissue cancellation takes place. Use of logarithmic signal processing extends this range, but bone thickness variations may not be accomodated because only two x-ray energies are involved in the imaging process. Use of a 3-kVp, 3-filter technique with logarithmic signal processing is shown to solve this problem. Computer simulations show that 1-mg/cm2 iodine may be imaged in the presence of 10 cm or more tissue variations and 2000-mg/cm2 bone variations.

It has long been recognized that the problems of motion artifacts in conventional time subtraction digital subtraction angiography (DSA) may be overcome using energy subtraction techniques. Of the variety of energy subtraction techniques investigated, non-k-edge dual-energy subtraction offers the best signal-to-noise ratio (SNR). However, this technique achieves only 55% of the temporal DSA SNR. Noise reduction techniques that average the noisier high-energy image produce various degrees of noise improvement while minimally affecting iodine contrast and resolution. A more significant improvement in dual-energy DSA iodine SNR, however, results when the correlated noise that exists in material specific images is appropriately cancelled. The correlated noise reduction (CNR) algorithm presented here follows directly from the dual-energy computed tomography work of Kalender who made explicit use of noise correlations in material specific images to reduce noise. The results are identical to those achieved using a linear version of the two-stage filtering process described by Macovski in which the selective image is filtered to reduce high-frequency noise and added to a weighted, high SNR, nonselective image which has been processed with a high-frequency bandpass filter. The dual-energy DSA CNR algorithm presented here combines selective tissue and iodine images to produce a significant increase in the iodine SNR while fully preserving iodine spatial resolution. Theoretical calculations predict a factor of 2-4 improvement in SNR compared to conventional dual-energy images. The improvement factor achieved is dependent upon the x-ray beam spectra and the size of blurring kernel used in the algorithm.(ABSTRACT TRUNCATED AT 250 WORDS)

Several investigators have approximated the scatter component of digital radiographic images by applying a constant weighting factor to a convolved version of the detected image. The scatter approximation is then subtracted from the detected image. When a constant weighting factor is used, the scatter in areas of low transmission is underestimated. We have extended this technique by allowing for a spatially variable weighting factor based on the local image intensity. This technique improves the scatter estimate and should provide better results for general videodensitometric applications.

Iodine is commonly used as a contrast material in computerized x-ray tomography. In some cases the determination of the iodine distribution in the image may be prevented by the presence of bone or tissue variations within the tomographic slice. This paper describes a method for quantitative selective imaging of the iodine concentration in the slice. The method employs scans using three heavily filtered x-ray beams, two having mean energies which straddle the iodine K edge (33 keV) and another at a slightly higher energy. The results are independent of tissue and bone over a broad range of projection path lengths. It is shown that, for separation of iodine from one other material, a two-beam K-edge approach requires less integral dose than a two-beam technique at conventional CT energies for slice diameters up to 30 cm. For selective iodine imaging in the presence of more than one other material, the three-spectrum K-edge technique is a necessity. Exposure requirements and beam-hardening corrections are discussed in detail and a computer-simulated CT image generated by the proposed scheme is presented.

In an earlier article we discussed the rationale for using differences between video images in three-beam selective iodine K-edge imaging. Rather than combining three initial images Li linearly to yield the final image k1L1 + k2L2 + k3L3, differences between the Li were first generated and then combined either to linear or quadratic order. This approach was motivated by the desire to suppress the large multiplicative biases of fluoroscopic imaging and justified by theoretically proving that k1 + k2 + k3 is approximately equal to 0. In this paper we discuss the instrumentation and experimental results obtained from this difference-based technique. A specially-constructed apparatus is described which automatically selects the optimum combination coefficients and combines the difference images up to quadratic order at realtime video rates. Three methods for generating K-edge subtraction images are compared: the former approach in which the Li are linearly combined and combination of differences to linear and quadratic order. In imaging phantoms in which the iodine distribution is known, the resultant subtraction images from all three methods appear similar. Inspection of signal sizes shows that the quadratic difference-based approach provides superior bone and tissue residual suppression by about a factor of 2. In imaging phantoms in which the iodine distribution is unknown, incomplete suppression of x-ray scatter and image intensifier veiling glare prevent a quantitative comparison of performance of the three algorithms. An experiment verification is provided of the theorem which states that k1 + k2 + k3 is approximately equal to 0.

Previous attempts to use time subtraction intravenous digital subtraction angiography for ventricular imaging have been hampered by artifacts due to cardiac and respiratory motion. We have previously reported a motion-immune dual-energy technique in which kVp is switched between 60 and 120, at 300-500 mA, 30 times/s. In order to quantitate parameters such as ejection fraction and left ventricular volume, it is necessary to correct for scatter and veiling glare (SVG), which are the major sources of nonlinearities in videodensitometric digital subtraction angiography (DSA). In this report, a convolution filtering method has been investigated to estimate SVG in DSA images. In the first step, a grey level transformation of the detected image is utilized to get an estimated SVG image. In the second step this image is convolved to produce an image with appropriate spatial frequency content. Estimates of SVG in several Humanoid chest phantom images were obtained using Gaussian convolution kernels with a full width at half-maximum (FWHM) of 51-125 pixels. The root-mean-square (rms) percentage error of these estimates was obtained by comparison with direct SVG measurement. A convolution kernel with a FWHM of 75 pixels in each dimension applied to 16 Humanoid phantom images with various projections, thicknesses, and beam energies resulted in an average rms percentage error of 9.7% in the SVG estimate, for the 16 cases studied. The SVG estimation consisting of grey scale-to-SVG fraction lookup table (LUT) is made based on previous measurements. The x-ray settings required for each patient are utilized to alter the LUT in order to account for patient thickness variations.(ABSTRACT TRUNCATED AT 250 WORDS)