Development of accelerated angiography techniques.

This article summarizes a progression of techniques designed to provide higher spatial and temporal resolution for angiographic acquisition and, in some cases, significant dose reduction. These methods were developed over a time period from 1976 to the present.

Optimization of quantitative time-resolved 3D (4D) digital subtraction angiography in a porcine liver model.

BACKGROUND: Time-resolved three-dimensional digital subtraction angiography (4D-DSA) can be used to quantify blood velocity. Contrast pulsatility, a major discriminant on 4D-DSA, is yet to be optimized. We investigated the effects of different imaging and injection parameters on sideband ratio (SBR), a measure of contrast pulsatile strength, within the hepatic vasculature of an in vivo porcine model. METHODS: Fifty-nine hepatic 4D-DSA procedures were performed in three female domestic swine (mean weight 54 kg). Contrast injections were performed in the common hepatic artery with different combinations of imaging duration (6 s or 12 s), injection rates (from 1.0 to 2.5 mL/s), contrast concentration (50% or 100%), and catheter size (4 Fr or 5 Fr). Reflux was recorded. SBR and vessel cross-sectional areas were calculated in 289 arterial segments. Multiple linear mixed-effects models were estimated to determine the effects of parameters on SBR and cross-sectional vessel area. RESULTS: Twelve-second acquisitions yielded a SBR higher than 6 s (p < 0.001). No significant differences in SBR were seen between different catheter sizes (p = 0.063) or contrast concentration (p = 0.907). For higher injection rates (2.5 mL/s), SBR was lower (p = 0.007) and cross-sectional area was higher (p < 0.001). Reflux of contrast does not significantly affect SBR (p = 0.087). CONCLUSIONS: The strength of contrast pulsatility used for flow quantitation with 4D-DSA can be increased by adjusting injection rates and using longer acquisition times. Reduction of contrast concentration to 50% is feasible and reflux of contrast does not significantly hinder contrast pulsatility.

Quantitative 4D-Digital Subtraction Angiography to Assess Changes in Hepatic Arterial Flow during Transarterial Embolization: A Feasibility Study in a Swine Model.

PURPOSE: To determine the feasibility of using time-resolved 3D-digital subtraction angiography (4D-DSA) for quantifying changes in hepatic arterial blood flow and velocity during transarterial embolization. MATERIALS AND METHODS: Hepatic arteriography and selective transarterial embolization were performed in 4 female domestic swine (mean weight, 54 kg) using 100-300-µm microspheres. Conventional 2D and 4D-DSA were performed before, during, and after each embolization. From the 4D-DSA reconstructions, blood flow and velocity values were calculated for hepatic arterial branches using a pulsatility-based algorithm. 4D-DSA velocity values were compared to those measured using an intravascular Doppler wire with a linear regression analysis. Paired t-tests were used to compare data before and after embolization. RESULTS: There was a weak-to-moderate but statistically significant correlation of flow velocities measured with 4D-DSA and the Doppler wire (r = 0.35, n = 39, P = .012). For vessels with high pulsatility, the correlation was higher (r = 0.64, n = 11, P = .034), and the relationship between 4D-DSA and the Doppler wire fit a linear model with a positive bias toward the Doppler wire (failed to reject at 95% confidence level, P = .208). 4D-DSA performed after partial embolization showed a reduction in velocity in the embolized hepatic arteries compared to pre-embolization (mean, 3.96 ± 0.74 vs 11.8 2± 2.15 cm/s, P = .006). CONCLUSION: Quantitative 4D-DSA can depict changes in hepatic arterial blood velocity during transarterial embolization in a swine model. Further work is needed to optimize 4D-DSA acquisitions and to investigate its applicability in humans.

Ultra-Low Radiation Dose CT Fluoroscopy for Percutaneous Interventions: A Porcine Feasibility Study.

Purpose To determine the feasibility of ultra-low-dose (ULD) CT fluoroscopy for performing percutaneous CT-guided interventions in an in vivo porcine model and to compare radiation dose, spatial accuracy, and metal artifact for conventional CT versus CT fluoroscopy. Materials and Methods An in vivo swine model was used (n = 4, ∼50 kg) for 20 procedures guided by 246 incremental conventional CT scans (mean, 12.5 scans per procedure). The procedures were approved by the Institutional Animal Care and Use Committee and performed by two experienced radiologists from September 7, 2017, to August 8, 2018. ULD CT fluoroscopic acquisitions were simulated by using only two of 984 conventional CT projections to locate and reconstruct the needle, which was superimposed on a previously acquired and motion-compensated CT scan. The authors (medical physicists) compared the ULD CT fluoroscopy results to those of conventional CT for needle location, radiation dose, and metal artifacts using Deming regression and generalized mixed models. Results The average distance between the needle tip reconstructed using conventional CT and ULD CT fluoroscopy was 1.06 mm. Compared with CT fluoroscopy, the estimated dose for a percutaneous procedure, including planning acquisitions, was 0.99 mSv (21% reduction) for patients (effective dose) and 0.015 µGy (97% reduction) for physicians (scattered dose in air). Metal artifacts were statistically significantly reduced (P < .001, bootstrapping), and the average registration error of the motion compensation was within 1-3 mm. Conclusion Ultra-low-dose CT fluoroscopy has the potential to reduce radiation exposure for intraprocedural scans to patients and staff by a factor of approximately 500 times compared with conventional CT acquisition, while maintaining image quality without metal artifacts. © RSNA, 2019.

Measuring blood velocity using 4D-DSA: A feasibility study.

PURPOSE: Four-dimensional (4D) DSA reconstruction provides three-dimensional (3D) time-resolved visualization of contrast bolus passage through arterial vasculature in the interventional setting. The purpose of this study was to evaluate the feasibility of using these data in measuring blood velocity and flow. METHODS: The pulsatile signals in the time concentration curves (TCCs) measured at different points along a vessel are markers of the movement of a contrast bolus and thus of blood flow. When combined with the spatial content, that is, geometry of the vasculature, this information then provides the data required to determine blood velocity. A Fourier-based algorithm was used to identify and follow the pulsatility signal. A Side Band Ratio (SBR) metric was used to reduce uncertainty in identifying the pulsatility in regions where the signal was weak. We tested this method using 4D-DSA reconstructions from vascular phantoms as well as from human studies. RESULTS: In five studies using 3D printed patient-specific cerebrovascular phantoms, velocities calculated from the 4D-DSAs were found to be within 10% of velocities measured with a flow meter. Calculated velocity and flow values from three human studies were within the range of those reported in the literature. CONCLUSIONS: 4D-DSA provides temporal and spatial information about blood flow and vascular geometry. This information is obtained using conventional rotational angiographic systems. In this small feasibility study, these data allowed calculations of velocity values that correlated well with measured values. The availability of velocity and blood flow information in the interventional setting would support a more quantitative approach to diagnosis, treatment planning and post-treatment evaluations of a variety of cerebrovascular diseases.

4D DSA reconstruction using tomosynthesis projections.

We investigate the use of tomosynthesis in 4D DSA to improve the accuracy of reconstructed vessel time-attenuation curves (TACs). It is hypothesized that a narrow-angle tomosynthesis dataset for each time point can be exploited to reduce artifacts caused by vessel overlap in individual projections. 4D DSA reconstructs time-resolved 3D angiographic volumes from a typical 3D DSA scan consisting of mask and iodine-enhanced C-arm rotations. Tomosynthesis projections are obtained either from a conventional C-arm rotation, or from an inverse geometry scanning-beam digital x-ray (SBDX) system. In the proposed method, rays of the tomosynthesis dataset which pass through multiple vessels can be ignored, allowing the non-overlapped rays to impart temporal information to the 4D DSA. The technique was tested in simulated scans of 2 mm diameter vessels separated by 2 to 5 cm, with TACs following either early or late enhancement. In standard 4D DSA, overlap artifacts were clearly present. Use of tomosynthesis projections in 4D DSA reduced TAC artifacts caused by vessel overlap, when a sufficient fraction of non-overlapped rays was available in each time frame. In cases where full overlap between vessels occurred, information could be recovered via a proposed image space interpolation technique. SBDX provides a tomosynthesis scan for each frame period in a rotational acquisition, whereas a standard C-arm geometry requires the grouping of multiple frames.

Feasibility of reduced-dose three-dimensional/four-dimensional-digital subtraction angiogram using a weighted edge preserving filter.

A conventional three-dimensional/four-dimensional (3D/4D) digital subtraction angiogram (DSA) requires two rotational acquisitions (mask and fill) to compute the log-subtracted projections that are used to reconstruct a 3D/4D volume. Since all of the vascular information is contained in the fill acquisition, it is hypothesized that it is possible to reduce the x-ray dose of the mask acquisition substantially and still obtain subtracted projections adequate to reconstruct a 3D/4D volume with noise level comparable to a full-dose acquisition. A full-dose mask and fill acquisition were acquired from a clinical study to provide a known full-dose reference reconstruction. Gaussian noise was added to the mask acquisition to simulate a mask acquisition acquired at 10% relative dose. Noise in the low-dose mask projections was reduced with a weighted edge preserving filter designed to preserve bony edges while suppressing noise. Two-dimensional (2D) log-subtracted projections were computed from the filtered low-dose mask and full-dose fill projections, and then 3D/4D-DSA reconstruction algorithms were applied. Additional bilateral filtering was applied to the 3D volumes. The signal-to-noise ratio measured in the filtered 3D/4D-DSA volumes was compared to the full-dose case. The average ratio of filtered low-dose SNR to full-dose SNR was 0.856 for the 3D-DSA and 0.849 for the 4D-DSA, indicating that the method is a feasible approach to restoring SNR in DSA scans acquired with a low-dose mask. The method was also tested in a phantom study with full-dose fill and 22%-dose mask.

Time-of-Arrival Parametric Maps and Virtual Bolus Images Derived From Contrast-Enhanced Time-Resolved Radial Magnetic Resonance Angiography Improve the Display of Brain Arteriovenous Malformation Vascular Anatomy.

OBJECTIVES: Time-of-arrival (TOA) maps can be derived from high-resolution 4-dimensional (4D) contrast-enhanced magnetic resonance angiography (MRA) data sets to provide a quantitative description of contrast material arrival time in each voxel. This information can further be processed to create a compressed time evolution curve that virtually shortens the contrast bolus (virtual bolus [VB]). The purpose of this project was to determine whether TOA-enhanced 4D MRA and/or VB imaging improve the display of contrast kinetics in patients with vascular disease. METHODS: High-resolution whole-brain contrast-enhanced 4D MRA examinations with 1.2-second temporal reconstruction were acquired by using radial acquisition and highly constrained projection reconstruction (radial 4D contrast-enhanced HYPRFlow, abbreviated as HFMRA in this article) in 10 patients (8 patients with arteriovenous malformations [AVM], 1 patient with an arteriovenous fistula, and 1 patient with a high-grade intracranial stenosis). The TOA for each voxel was defined as the time point when the signal intensity reached 20% of its maximum. In the first method, TOA maps were generated, color-encoded, and then multiplied with the time-resolved contrast-enhanced MRA images at each time frame to form new 4D MRA images (TOA-enhanced HFMRA), which contains the contrast arrival times with defined color encoding. In the second method, each time frame was weighted by a Gaussian distribution in the time domain to form a virtual 4D bolus map. This 4D bolus map was then color-coded and multiplied with the HFMRA images to form a digital subtraction angiography (DSA)-like VB, where at each time frame, only vessels with certain TOA values within the defined bolus length appear. HFMRA, TOA maps, and VB images were scored qualitatively with regard to delineation of arteries, veins, and nidus, as well as artifacts. Furthermore, diagnostic confidence and arteriovenous overlap were evaluated and compared between techniques. A comparison with DSA was performed where DSA served as the reference standard in terms of number of arterial feeders, draining veins, and Spetzler-Martin score of AVMs. In addition, TOA maps were evaluated quantitatively. RESULTS: Overall, diagnostic confidence score of TOA was significantly higher compared with that of HFMRA (P = 0.03). Virtual bolus showed significantly higher scores for overall diagnostic confidence (P = 0.02) and reduced arteriovenous overlap (0.01) compared with HFMRA. Furthermore, VB-reduced arteriovenous overlap scores were significantly higher compared with TOA (P = 0.04). Agreement regarding AVM draining veins was lower between DSA and HFMRA (κ = 0.3) compared with TOA and VB (κ = 0.56). Agreement regarding Spetzler-Martin score was lower between DSA and HFMRA (κ = 0.56) compared with TOA and VB (κ = 0.74). CONCLUSIONS: TOA-enhanced HFMRA provides serial images and time of arrival maps in one inclusive display. In this study, TOA mapping combined with Virtual Bolus imaging improved diagnostic confidence in AVM patients and facilitated arteriovenous separation. The VB method further reduced overlap of arterial and venous structures.

4D interventional device reconstruction from biplane fluoroscopy.

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.

Volumetric limiting spatial resolution analysis of four-dimensional digital subtraction angiography.

C-Arm CT three-dimensional (3-D) digital subtraction angiography (DSA) reconstructions cannot provide temporal information to radiologists. Four-dimensional (4-D) DSA provides a time series of 3-D volumes utilizing temporal dynamics in the two-dimensional (2-D) projections using a constraining image reconstruction approach. Volumetric limiting spatial resolution (VLSR) of 4-D DSA is quantified and compared to a 3-D DSA. The effects of varying 4-D DSA parameters of 2-D projection blurring kernel size and threshold of the 3-D DSA (constraining image) of an in silico phantom (ISPH) and physical phantom (PPH) were investigated. The PPH consisted of a 76-micron tungsten wire. An [Formula: see text] scan protocol acquired the projection data. VLSR was determined from MTF curves generated from each 2-D transverse slice of every (248) 4-D temporal frame. 4-D DSA results for PPH and ISPH were compared to the 3-D DSA. 3-D DSA analysis resulted in a VLSR of 2.28 and [Formula: see text] for ISPH and PPH, respectively. Kernel sizes of either [Formula: see text] or [Formula: see text] with a 3-D DSA constraining image threshold of 10% provided 4-D DSA VLSR nearest to the 3-D DSA. 4-D DSA yielded 2.21 and [Formula: see text] with a percent error of 3.1 and 1.2% for ISPH and PPH, respectively, as compared to 3-D DSA. This research indicates 4-D DSA is capable of retaining the resolution of 3-D DSA.

Accelerated Time-Resolved Contrast-Enhanced Magnetic Resonance Angiography of Dural Arteriovenous Fistulas Using Highly Constrained Reconstruction of Sparse Cerebrovascular Data Sets.

OBJECTIVE: Time-resolved contrast-enhanced magnetic resonance angiography (MRA) is commonly used to noninvasively characterize vascular malformations. However, the spatial and temporal resolution of current methods often compromises the clinical value of the examinations. Constrained reconstruction is a temporal spatial correlation strategy that exploits the relative sparsity of vessels in space to dramatically reduce the amount of data required to generate fast high-resolution time-resolved contrast-enhanced MRA studies. In this report, we use a novel temporal spatial acceleration method termed HYPRFlow to diagnose and classify dural arteriovenous fistulas (DAVFs). Our hypothesis is that HYPRFlow images are of adequate diagnostic image quality to delineate the arterial and venous components of DAVFs and allow correct classification using the Cognard system. SUBJECTS AND METHODS: Eight patients with known DAVFs underwent HYPRFlow imaging with isotropic resolution of 0.68 mm and temporal resolution of 0.75 second and 3-dimensional time-of-flight (3DTOF) MRA. The 3DTOF images and HYPRFlow images were evaluated by 2 readers and scored for arterial anatomic image quality. Digital subtraction angiography (DSA) was available for comparison in 7 subjects, and for these patients, each DAVF was classified according to the Cognard system using HYPRFlow and DSA examinations. Digital subtraction angiography was considered the reference examination or criterion standard. RESULTS: HYPRFlow imaging classification was concordant with DSA in all but 1 case. There was no difference in the arterial image quality scores between HYPRFlow and 3DTOF MRA (95% confidence interval). Arterial-to-venous separation was rated excellent (n = 3), good (n = 4), or poor (n = 1), and arteriovenous shunting was easily appreciated. Undersampling artifacts were reduced by using a low pass filter and did not interfere with the diagnostic quality of the examinations. CONCLUSIONS: HYPRFlow is a novel acquisition and reconstruction technique that exploits the relative sparsity of intracranial vessels in space to increase temporal and spatial resolution and provides accurate delineation of DAVF vasculature.

Noise reduction for curve-linear structures in real time fluoroscopy applications using directional binary masks.

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.

Simultaneous static and cine nonenhanced MR angiography using radial sampling and highly constrained back projection reconstruction.

PURPOSE: To describe a pulse sequence for simultaneous static and cine nonenhanced magnetic resonance angiography (NEMRA) of the peripheral arteries. METHODS: The peripheral arteries of 10 volunteers and 6 patients with peripheral arterial disease (PAD) were imaged with the proposed cine NEMRA sequence on a 1.5 Tesla (T) system. The impact of multi-shot imaging and highly constrained back projection (HYPR) reconstruction was examined. The propagation rate of signal along the length of the arterial tree in the cine nonenhanced MR angiograms was quantified. RESULTS: The cine NEMRA sequence simultaneously provided a static MR angiogram showing vascular anatomy as well as a cine display of arterial pulse wave propagation along the entire length of the peripheral arteries. Multi-shot cine NEMRA improved temporal resolution and reduced image artifacts. HYPR reconstruction improved image quality when temporal reconstruction footprints shorter than 100 ms were used (P < 0.001). Pulse wave propagation within the arterial tree as displayed by cine NEMRA was slower in patients with PAD than in volunteers. CONCLUSION: Simultaneous static and cine NEMRA of the peripheral arteries is feasible. Multi-shot acquisition and HYPR reconstruction can be used to improve arterial conspicuity and temporal resolution.

Noncontrast dynamic 3D intracranial MR angiography using pseudo-continuous arterial spin labeling (PCASL) and accelerated 3D radial acquisition.

PURPOSE: To develop a novel dynamic 3D noncontrast magnetic resonance angiography (MRA) technique that combines dynamic pseudo-continuous arterial spin labeling (dynamic PCASL), accelerated 3D radial sampling (VIPR), and time-of-arrival (TOA) mapping to provide quantitative assessment of arterial flow. MATERIALS AND METHODS: Digital simulations were performed to investigate the effects of acquisition scheme and sequence parameters on image quality and TOA mapping fidelity. Five patients with vascular malformations (arteriovenous malformation [AVM] = 3, dural arteriovenous fistula [DAVF] = 2) were scanned and the images were compared to digital subtraction angiography (DSA) for the ability to identify the arterial supply, AVM location, nidus size, and venous drainage. RESULTS: Digital simulations demonstrated reduced image artifacts and improved TOA accuracy using radial acquisition over Cartesian. TOA mapping accuracy is more sensitive to sampling window length than time spacing. Dynamic PCASL MRA depicted seven of eight arterial pedicles, and accurately measured the AVM nidus size when the nidus was compact. The venous drainage in the AVM patients was not consistently visualized. CONCLUSION: Dynamic 3D PCASL-VIPR with TOA mapping is able to acquire both high temporal and spatial resolution inflow dynamics that could improve diagnosis of high-flow intracranial vascular diseases.

Design of a digital beam attenuation system for computed tomography: part I. System design and simulation framework.

PURPOSE: The purpose of this work is to introduce a new device that allows for patient-specific imaging-dose modulation in conventional and cone-beam CT. The device is called a digital beam attenuator (DBA). The DBA modulates an x-ray beam by varying the attenuation of a set of attenuating wedge filters across the fan angle. The ability to modulate the imaging dose across the fan beam represents another stride in the direction of personalized medicine. With the DBA, imaging dose can be tailored for a given patient anatomy, or even tailored to provide signal-to-noise ratio enhancement within a region of interest. This modulation enables decreases in: dose, scatter, detector dynamic range requirements, and noise nonuniformities. In addition to introducing the DBA, the simulation framework used to study the DBA under different configurations is presented. Finally, a detailed study on the choice of the material used to build the DBA is presented. METHODS: To change the attenuator thickness, the authors propose to use an overlapping wedge design. In this design, for each wedge pair, one wedge is held stationary and another wedge is moved over the stationary wedge. The composite thickness of the two wedges changes as a function of the amount of overlap between the wedges. To validate the DBA concept and study design changes, a simulation environment was constructed. The environment allows for changes to system geometry, different source spectra, DBA wedge design modifications, and supports both voxelized and analytic phantom models. A study of all the elements from atomic number 1 to 92 were evaluated for use as DBA filter material. The amount of dynamic range and tube loading for each element were calculated for various DBA designs. Tube loading was calculated by comparing the attenuation of the DBA at its minimum attenuation position to a filtered non-DBA acquisition. RESULTS: The design and parametrization of DBA implemented FFMCT has been introduced. A simulation framework was presented with which DBA-FFMCT, bowtie filter CT acquisitions, and unmodulated CT acquisitions can be simulated. The study on wedge filter design concluded that the ideal filter material should have an atomic number in the range of 21-34. Iron was chosen for an experimental relative-tube-loading measurement and showed that DBA-FFMCT scans could be acquired with negligible increases in tube power demands. CONCLUSIONS: The basic idea of DBA implemented fluence field modulated CT, a simulation framework to verify the concept, and a filter selection study have been presented. The use of a DBA represents another step toward the ultimate in patient specific CT dose delivery as patient dose can be delivered uniquely as a function of view and fan angle using this device.

Design of a digital beam attenuation system for computed tomography. Part II. Performance study and initial results.

PURPOSE: The purpose of this work is to present a performance study of the digital beam attenuator (DBA) for implementing fluence field modulated CT (FFMCT) using a simulation framework developed to model the incorporation of the DBA into an existing CT system. Additionally, initial results will be presented using a prototype DBA and the realization of the prototype will be described. To our knowledge, this study represents the first experimental use of a device capable of modulating x-ray fluence as a function of fan angle using a CT geometry. METHODS: To realize FFMCT, the authors propose to use a wedge design in which one wedge is held stationary and another wedge is moved over the stationary wedge. Due to the wedge shape, the composite thickness of the two wedges changes as a function of the amount of overlap between the wedges. This design allows for the wedges to modulate the photon fluence incident onto a patient. Using a simulation environment, the effect of changing the number of wedges has on dose, scatter, detector dynamic range, and noise uniformity is explored. Experimental results are presented using a prototype DBA having ten Fe wedges and a c-arm CT system geometry. The experimental DBA results are compared to non-DBA scans using scatter and detector dynamic range as metrics. Both flat field and bowtie filtered CT acquisitions were simulated for comparison with the DBA. RESULTS: Numerical results suggest that substantial gains in noise uniformity and scatter-to-primary ratio (SPR) can be obtained using only seven wedges. After seven wedges, the decrease in noise ununiformity and SPR falls off at a lower rate. Simulations comparing CT acquisitions between flat field, bowtie enabled, and DBA CT acquisitions suggest DBA-FFMCT can reduce dose relative to flat field CT by ≈3 times. A bowtie filter under the same imaging conditions was shown to only allow a dose reduction of 1.65 times. Experimentally, a 10 wedge DBA prototype result showed a SPR reduction of ≈4 times relative to flat field CT. The dynamic range for the DBA prototype was 3.7 compared to 84.2 for the flat field scan. CONCLUSIONS: Based on the results presented in this paper and the companion paper [T. Szczykutowicz and C. Mistretta, "Design of a digital beam attenuation system for computed tomography. Part I. System design and simulation framework," Med. Phys. 40, 021905 (2013)], FFMCT implemented via the DBA device seems feasible and should result in both a dose reduction and an improvement in image quality as judged by noise uniformity and scatter reduction. In addition, the dynamic range reduction achievable using the DBA may allow photon counting imaging to become a clinical reality. This study may allow for yet another step to be taken in the field of patient specific dose modulation.

Noncontrast-enhanced three-dimensional (3D) intracranial MR angiography using pseudocontinuous arterial spin labeling and accelerated 3D radial acquisition.

Pseudocontinuous arterial spin labeling (PCASL) can be used to generate noncontrast magnetic resonance angiograms of the cerebrovascular structures. Previously described PCASL-based angiography techniques were limited to two-dimensional projection images or relatively low-resolution three-dimensional (3D) imaging due to long acquisition time. This work proposes a new PCASL-based 3D magnetic resonance angiography method that uses an accelerated 3D radial acquisition technique (VIPR, spoiled gradient echo) as the readout. Benefiting from the sparsity provided by PCASL and noise-like artifacts of VIPR, this new method is able to obtain submillimeter 3D isotropic resolution and whole head coverage with a 8-min scan. Intracranial angiography feasibility studies in healthy (N = 5) and diseased (N = 5) subjects show reduced saturation artifacts in PCASL-VIPR compared with a standard time-of-flight protocol. These initial results show great promise for PCASL-VIPR for static, dynamic, and vessel selective 3D intracranial angiography.

Improved kinetic analysis of dynamic PET data with optimized HYPR-LR.

PURPOSE: Highly constrained backprojection-local reconstruction (HYPR-LR) has made a dramatic impact on magnetic resonance angiography (MRA) and shows promise for positron emission tomography (PET) because of the improvements in the signal-to-noise ratio (SNR) it provides dynamic images. For PET in particular, HYPR-LR could improve kinetic analysis methods that are sensitive to noise. In this work, the authors closely examine the performance of HYPR-LR in the context of kinetic analysis, they develop an implementation of the algorithm that can be tailored to specific PET imaging tasks to minimize bias and maximize improvement in variance, and they provide a framework for validating the use of HYPR-LR processing for a particular imaging task. METHODS: HYPR-LR can introduce errors into non sparse PET studies that might bias kinetic parameter estimates. An implementation of HYPR-LR is proposed that uses multiple temporally summed composite images that are formed based on the kinetics of the tracer being studied (HYPR-LR-MC). The effects of HYPR-LR-MC and of HYPR-LR using a full composite formed with all the frames in the study (HYPR-LR-FC) on the kinetic analysis of Pittsburgh compound-B ([11C]-PIB) are studied. HYPR-LR processing is compared to spatial smoothing. HYPR-LR processing was evaluated using both simulated and human studies. Nondisplaceable binding potential (BP(ND)) parametric images were generated from fifty noise realizations of the same numerical phantom and eight [(11)C]-PIB positive human scans before and after HYPR-LR processing or smoothing using the reference region Logan graphical method and receptor parametric mapping (RPM2). The bias and coefficient of variation in the frontal and parietal cortex in the simulated parametric images were calculated to evaluate the absolute performance of HYPR-LR processing. Bias in the human data was evaluated by comparing parametric image BP(ND) values averaged over large regions of interest (ROIs) to Logan estimates of the BP(ND) from TACs averaged over the same ROIs. Variance was assessed qualitatively in the parametric images and semiquantitatively by studying the correlation between voxel BP(ND) estimates from Logan analysis and RPM2. RESULTS: Both the simulated and human data show that HYPR-LR-FC overestimates BP(ND) values in regions of high [(11)C]-PIB uptake. HYPR-LR-MC virtually eliminates this bias. Both implementations of HYPR-LR reduce variance in the parametric images generated with both Logan analysis and RPM2, and HYPR-LR-FC provides a greater reduction in variance. This reduction in variance nearly eliminates the noise-dependent Logan bias. The variance reduction is greater for the Logan method, particularly for HYPR-LR-MC, and the variance in the resulting Logan images is comparable to that in the RPM2 images. HYPR-LR processing compares favorably with spatial smoothing, particularly when the data are analyzed with the Logan method, as it provides a reduction in variance with no loss of spatial resolution. CONCLUSIONS: HYPR-LR processing shows significant potential for reducing variance in parametric images, and can eliminate the noise-dependent Logan bias. HYPR-LR-FC processing provides the greatest reduction in variance but introduces a positive bias into the BP(ND) of high-uptake border regions. The proposed method for forming HYPR composite images, HYPR-LR-MC, eliminates this bias at the cost of less variance reduction.

Time-resolved angiography: Past, present, and future.

The introduction of digital subtraction angiography (DSA) in 1980 provided a method for real time 2D subtraction imaging. Later, 4D magnetic resonance (MR) angiography emerged beginning with techniques like Keyhole and time-resolved imaging of contrast kinetics (TRICKS) that provided frame rates of one every 5 seconds with limited spatial resolution. Undersampled radial acquisition was subsequently developed. The 3D vastly undersampled isotropic projection (VIPR) technique allowed undersampling factors of 30-40. Its combination with phase contrast displays time-resolved flow dynamics within the cardiac cycle and has enabled the measurement of pressure gradients in small vessels. Meanwhile similar accelerations were achieved using Cartesian acquisition with projection reconstruction (CAPR), a Cartesian acquisition with 2D parallel imaging. Further acceleration is provided by constrained reconstruction techniques such as highly constrained back-projection reconstruction (HYPR) and its derivatives, which permit acceleration factors approaching 1000. Hybrid MRA combines a separate phase contrast, time-of flight, or contrast-enhanced acquisition to constrain the reconstruction of contrast-enhanced time frames providing exceptional spatial and temporal resolution and signal-to-noise ratio (SNR). This can be extended to x-ray imaging where a 3D DSA examination can be used to constrain the reconstruction of time-resolved 3D volumes. Each 4D DSA (time-resolved 3D DSA) frame provides spatial resolution and SNR comparable to 3D DSA, thus removing a major limitation of intravenous DSA. Similar techniques have provided the ability to do 4D fluoroscopy.