4D digital subtraction angiography: implementation and demonstration of feasibility

Simulation-informed learning for time-resolved angiographic contrast agent concentration reconstruction

Three-dimensional Digital Subtraction Angiography (3D-DSA) is a well-established X-ray-based technique for visualizing vascular anatomy. Recently, four-dimensional DSA (4D-DSA) reconstruction algorithms have been developed to enable the visualization of volumetric contrast flow dynamics through time-series of volumes. This reconstruction problem is ill-posed mainly due to vessel overlap in the projection direction and geometric vessel foreshortening, which leads to information loss in the recorded projection images. However, knowledge about the underlying fluid dynamics can be leveraged to constrain the solution space. In our work, we implicitly include this information in a neural network-based model that is trained on a dataset of image-based blood flow simulations. The model predicts the spatially averaged contrast agent concentration for each centerline point of the vasculature over time, lowering the overall computational demand. The trained network enables the reconstruction of relative contrast agent concentrations with a mean absolute error of 0.02±0.02 and a mean absolute percentage error of 5.31±9.25 %. Moreover, the network is robust to varying degrees of vessel overlap and vessel foreshortening. Our approach demonstrates the potential of the integration of machine learning and blood flow simulations in time-resolved angiographic contrast agent concentration reconstruction.

Improved temporal resolution and acceleration on 4D-MR angiography based on superselective pseudo-continuous arterial spin labeling combined with CENTRA-keyhole and view-sharing (4D-S-PACK) using an interpolation algorithm on the temporal axis and compressed sensing-sensitivity encoding (CS-SENSE)

PURPOSE

Two major drawbacks of 4D-MR angiography based on superselective pseudo-continuous arterial spin labeling combined with CENTRA-keyhole and view-sharing (4D-S-PACK) are the low temporal resolution and long scanning time. We investigated the feasibility of increasing the temporal resolution and accelerating the scanning time on 4D-S-PACK by using CS-SENSE and PhyZiodynamics, a novel image-processing program that interpolates images between phases to generate new phases and reduces image noise.

METHODS

Seven healthy volunteers were scanned with a 3.0 T MR scanner to visualize the internal carotid artery (ICA) system. PhyZiodynamics is a novel image-processing that interpolates images between phases to generate new phases and reduces image noise, and by increasing temporal resolution using PhyZiodynamics, inflow dynamic data (reference) were acquired by changing the labeling durations (100-2000 msec, 31 phases) in 4D-S-PACK. From this set of data, we selected seven time intervals to calculate interpolated time points with up to 61 intervals using ×10 for the generation of interpolated phases with PhyZiodynamics. In the denoising process of PhyZiodynamics, we processed the none, low, medium, high noise reduction dataset images. The time intensity curve (TIC), the contrast-to-noise ratio (CNR) were evaluated. In accelerating with CS-SENSE for 4D-S-PACK, 4D-S-PACK were scanned different SENSE or CS-SENSE acceleration factors: SENSE3, CS3-6. Signal intensity (SI), CNR, were evaluated for accelerating the 4D-S-PACK. With regard to arterial vascular visualization, we evaluated the middle cerebral artery (MCA: M1-4 segments).

RESULTS

In increasing temporal resolution, the TIC showed a similar trend between the reference dataset and the interpolated dataset. As the noise reduction weight increased, the CNR of the interpolated dataset were increased compared to that of the reference dataset. In accelerating 4D-S-PACK, the SI values of the SENSE3 dataset and CS dataset with CS3-6 were no significant differences. The image noise increased with the increase of acceleration factor, and the CNR decreased with the increase of acceleration factor. Significant differences in CNR were observed between acceleration factor of SENSE3 and CS6 for the M1-4 (P < 0.05). Visualization of small arteries (M4) became less reliable in CS5 or CS6 images. Significant differences were found for the scores of M2, M3 and M4 segments between SENSE3 and CS6.

CONCLUSION

With PhyZiodynamics and CS-SENSE in 4D-S-PACK, we were able to shorten the scan time while improving the temporal resolution.

Uterine Arteriovenous Malformation: Diagnostic and Therapeutic Challenges

Uterine arteriovenous malformations are a rare cause of puerperal haemorrhage, but their incidence is increasing due to both improved diagnosis and the more frequent use of uterine surgery in recent years. The use of ultrasound, both B-mode and Doppler, is recommended for diagnosis and follow-up, as it has been shown to be the simplest and most cost-effective method. Endometrial thickening associated with an anechoic and vascular intramiometrial structure is very useful for diagnosis and can help to exclude other causes of dysfunctional bleeding. Pulsed Doppler shows low-resistance vessels and high pulsatility indices with a high peak systolic velocity (PSV). In a healthy myometrium, the vessels have a peak systolic velocity of 9-40 cm/s and a resistance index between 0.6 and 0.8, whereas in the case of AVMs, the systolic and diastolic velocities are 4-6 times higher (PSV 25-110 cm/s with a mean of 60 cm/s and a resistance index of 0.27-0.75 with a mean of 0.41). For treatment, we must individualise each case, taking into account haemodynamic stability, the patient's reproductive wishes, and the severity of the AVM as assessed by its size and PSV.

Enhancing cerebral vasculature analysis with pathlength-corrected 2D angiographic parametric imaging: A feasibility study

BACKGROUND

2D angiographic parametric imaging (API) quantitatively extracts imaging biomarkers related to contrast flow and is conventionally applied to 2D digitally subtracted angiograms (DSA's). In the interventional suite, API is typically performed using 1-2 projection views and is limited by vessel overlap, foreshortening, and depth-integration of contrast motion.

PURPOSE

This work explores the use of a pathlength-correction metric to overcome the limitations of 2D-API: the primary objective was to study the effect of converting 3D contrast flow to projected contrast flow using a simulated angiographic framework created with computational fluid dynamics (CFD) simulations, thereby removing acquisition variability.

METHODS

The pathlength-correction framework was applied to in-silico angiograms, generating a reference (i.e., ground-truth) volumetric contrast distribution in four patient-specific intracranial aneurysm geometries. Biplane projections of contrast flow were created from the reference volumetric contrast distributions, assuming a cone-beam geometry. A Parker-weighted reconstruction was performed to obtain a binary representation of the vessel structure in 3D. Standard ray tracing techniques were then used to track the intersection of a ray from the focal spot with each voxel of the reconstructed vessel wall to a pixel in the detector plane. The lengths of each ray through the 3D vessel lumen were then projected along each ray-path to create a pathlength-correction map, where the pixel intensity in the detector plane corresponds to the vessel width along each source-detector ray. By dividing the projection sequences with this correction map, 2D pathlength-corrected in-silico angiograms were obtained. We then performed voxel-wise (3D) API on the ground-truth contrast distribution and compared it to pixel-wise (2D) API, both with and without pathlength correction for each biplane view. The percentage difference (PD) between the resultant API biomarkers in each dataset were calculated within the aneurysm region of interest (ROI).

RESULTS

Intensity-based API parameters, such as the area under the curve (AUC) and peak height (PH), exhibited notable changes in magnitude and spatial distribution following pathlength correction: these now accurately represent conservation of mass of injected contrast media within each arterial geometry and accurately reflect regions of stagnation and recirculation in each aneurysm ROI. Improved agreement was observed between these biomarkers in the pathlength-corrected biplane maps: the maximum PD within the aneurysm ROI is 3.3% with pathlength correction and 47.7% without pathlength correction. As expected, improved agreement with ROI-averaged ground-truth 3D counterparts was observed for all aneurysm geometries, particularly large aneurysms: the maximum PD for both AUC and PH was 5.8%. Temporal parameters (mean transit time, MTT, time-to-peak, TTP, time-to-arrival, TTA) remained unaffected after pathlength correction.

CONCLUSIONS

This study indicates that the values of intensity-based API parameters obtained with conventional 2D-API, without pathlength correction, are highly dependent on the projection orientation, and uncorrected API should be avoided for hemodynamic analysis. The proposed metric can standardize 2D API-derived biomarkers independent of projection orientation, potentially improving the diagnostic value of all acquired 2D-DSA's. Integration of a pathlength correction map into the imaging process can allow for improved interpretation of biomarkers in 2D space, which may lead to improved diagnostic accuracy during procedures involving the cerebral vasculature.

Interleaved x-ray imaging: A method for simultaneous acquisition of quantitative and diagnostic digital subtraction angiography

BACKGROUND

Flow altering angiographic procedures suffer from ill-defined, qualitative endpoints. Quantitative digital subtraction angiography (qDSA) is an emerging technology that aims to address this issue by providing intra-procedural blood velocity measurements from time-resolved, 2D angiograms. To date, qDSA has used 30 frame/s DSA imaging, which is associated with high radiation dose rate compared to clinical diagnostic DSA (up to 4 frame/s).

PURPOSE

The purpose of this study is to demonstrate an interleaved x-ray imaging method which decreases the radiation dose rate associated with high frame rate qDSA while simultaneously providing low frame rate diagnostic DSA images, enabling the acquisition of both datasets in a single image sequence with a single injection of contrast agent.

METHODS

Interleaved x-ray imaging combines low radiation dose image frames acquired at a high rate with high radiation dose image frames acquired at a low rate. The feasibility of this approach was evaluated on an x-ray system equipped with research prototype software for x-ray tube control. qDSA blood velocity quantification was evaluated in a flow phantom study for two lower dose interleaving protocols (LD1:3.7 ± 0.02 mGy / s $3.7 \pm 0.02\ {\mathrm{mGy}}/{\mathrm{s}}$ and LD2:1.7 ± 0.04 mGy / s $1.7 \pm 0.04{\mathrm{\ mGy}}/{\mathrm{s}}$ ) and one conventional (full dose) protocol (11.4 ± 0.04 mGy / s ) $11.4 \pm 0.04{\mathrm{\ mGy}}/{\mathrm{s}})$ . Dose was measured at the interventional reference point. Fluid velocities ranging from 24 to 45 cm/s were investigated. Gold standard velocities were measured using an ultrasound flow probe. Linear regression and Bland-Altman analysis were used to compare ultrasound and qDSA.

RESULTS

The LD1 and LD2 interleaved protocols resulted in dose rate reductions of -67.7% and -85.5%, compared to the full dose qDSA scan. For the full dose protocol, the Bland-Altman limits of agreement (LOA) between qDSA and ultrasound velocities were [0.7, 6.7] cm/s with a mean difference of 3.7 cm/s. The LD1 interleaved protocol results were similar (LOA: [0.3, 6.9] cm/s, bias: 3.6 cm/s). The LD2 interleaved protocol resulted in slightly larger LOA: [-2.5, 5.5] cm/s with a decrease in the bias: 1.5 cm/s. Linear regression analysis showed a strong correlation between ultrasound and qDSA derived velocities using the LD1 protocol, with aR 2 ${R}^2$ of0.96 $0.96$ , a slope of1.05 $1.05$ and an offset of1.9 $1.9$  cm/s. Similar values were also found for the LD2 protocol, with aR 2 ${R}^2$ of0.93 $0.93$ , a slope of0.98 $0.98$ and an offset of2.0 $2.0$  cm/s.

CONCLUSIONS

The interleaved method enables simultaneous acquisition of low-dose high-rate images for intra-procedural blood velocity quantification (qDSA) and high-dose low-rate images for vessel morphology evaluation (diagnostic DSA).

AIRPORT: A Data Consistency Constrained Deep Temporal Extrapolation Method To Improve Temporal Resolution In Contrast Enhanced CT Imaging

Typical tomographic image reconstruction methods require that the imaged object is static and stationary during the time window to acquire a minimally complete data set. The violation of this requirement leads to temporal-averaging errors in the reconstructed images. For a fixed gantry rotation speed, to reduce the errors, it is desired to reconstruct images using data acquired over a narrower angular range, i.e., with a higher temporal resolution. However, image reconstruction with a narrower angular range violates the data sufficiency condition, resulting in severe data-insufficiency-induced errors. The purpose of this work is to decouple the trade-off between these two types of errors in contrast-enhanced computed tomography (CT) imaging. We demonstrated that using the developed data consistency constrained deep temporal extrapolation method (AIRPORT), the entire time-varying imaged object can be accurately reconstructed with 40 frames-per-second temporal resolution, the time window needed to acquire a single projection view data using a typical C-arm cone-beam CT system. AIRPORT is applicable to general non-sparse imaging tasks using a single short-scan data acquisition.

Spatiotemporal frequency domain analysis for blood velocity measurement during embolization procedures

BACKGROUND

Currently, determining procedural endpoints and treatment efficacy of vascular interventions is largely qualitative and relies on subjective visual assessment of digital subtraction angiography (DSA) images leading to large interobserver variabilities and poor reproducibility. Quantitative metrics such as the residual blood velocity in embolized vessel branches could help establish objective and reproducible endpoints. Recently, velocity quantification techniques based on a contrast enhanced X-ray sequence such as qDSA and 4D DSA have been proposed. These techniques must be robust, and, to avoid radiation dose concerns, they should be compatible with low dose per frame image acquisition.

PURPOSE

To develop and evaluate a technique for robust blood velocity quantification from low dose contrast enhanced X-ray image sequences that leverages the oscillating signal created by pulsatile blood flow.

METHODS

The proposed spatiotemporal frequency domain (STF) approach quantifies velocities from time attenuation maps (TAMs) representing the oscillating signal over time for all points along a vessel centerline. Due to the time it takes a contrast bolus to travel along the vessel centerline, the resulting TAM resembles a sheared sine wave. The shear angle is related to the velocity and can be determined in the spatiotemporal frequency domain after applying the 2D Fourier transform to the TAM. The approach was evaluated in a straight tube phantom using three different radiation dose levels and compared to ultrasound transit-time-based measurements. The STF velocity results were also compared to previously published approaches for the measurement of blood velocity from contrast enhanced X-ray sequences including shifted least squared (SLS) and phase shift (PHS). Additionally, an in vivo porcine study (n = 8) was performed where increasing amounts of embolic particles were injected into a hepatic or splenic artery with intermittent velocity measurements after each injection to monitor the resulting reduction in velocity.

RESULTS

At the lowest evaluated dose level (average air kerma rate 1.3 mGy/s at the interventional reference point), the Pearson correlation between ultrasound and STF velocity measurements was99 % $99\%$ . This was significantly higher (p < 0.0001 $p < 0.0001$ ) than corresponding correlation results between ultrasound and the previously published SLS and PHS approaches (91 $\hskip.001pt 91$ and93 % $93\%$ , respectively). In the in vivo study, a reduction in velocity was observed in85.7 % $85.7\%$ of cases after injection of 1 mL,96.4 % $96.4\%$ after 3 mL, and100.0 % $100.0\%$ after 4 mL of embolic particles.

CONCLUSIONS

The results show good agreement of the spatiotemporal frequency domain approach with ultrasound even in low dose per frame image sequences. Additionally, the in vivo study demonstrates the ability to monitor the physiological changes due to embolization. This could provide quantitative metrics during vascular procedures to establish objective and reproducible endpoints.

Experimental study using phantom models of cerebral aneurysms and 4D-DSA to measure blood flow on 3D-color-coded images

BACKGROUND

The current 3D-iFlow application can only measure the arrival time of contrast media through intensity values. If the flow rate could be estimated by 3D-iFlow, patient-specific hemodynamics could be determined within the scope of normal diagnostic management, eliminating the need for additional resources for blood flow rate estimation.

OBJECTIVE

The aim of this study is to develop and validate a method for measuring the flow rate by data obtained from 3D-iFlow images - a prototype application in Four-dimensional digital subtraction angiography (4D-DSA).

METHODS

Using phantom model and experimental circuit with circulating glycerin solution, an equation for the relationship between contrast media intensity and flow rate was developed. Applying the equation to the aneurysm phantom models, the derived flow rate was evaluated.

RESULTS

The average errors between the derived flow rate and setting flow rate became larger when the glycerin flow and the X-rays from the X-ray tube of the angiography system were parallel to each other or when the measurement point included overlaps with other contrast enhanced areas.

CONCLUSION

Although the error increases dependent on the imaging direction and overlap of contrast enhanced area, the developed equation can estimate the flow rate using the image intensity value measured on 3D-iFlow based on 4D-DSA.

Updates on Neuronavigation: Emerging tools for tumor resection

Multiple studies have been conducted to properly elucidate the various tools available to help enhance the resection of tumor tissue, aneurysms, and arteriovenous malformations (AVM). Diffusion tensor imaging (DTI) tractography is useful in providing a map of the tumor borders, allowing the optimal preservation of function and structure of specific regions of the brain. During neurosurgery, especially craniotomies, the possibility of the brain shifting due to swelling or gravity is high. Thus, tools for intraoperative imaging such as high-frequency linear array ultrasound transducers and doppler ultrasonography are utilized for high resolution images and detecting frequency shifts. 4D-digital subtraction angiography (DSA) is another technique used to create spatial resolutions and 3D maps for aneurysms. These similar techniques can also be utilized to assess the integrity of white matter in AVM. By implementing effective evaluation strategies, healthcare professionals can make informed decisions regarding treatment options, preventive measures, and long-term care plans tailored to individual patients.

4D-DSA for Assessment of the Angioarchitecture and Grading of Cranial Dural AVF

BACKGROUND AND PURPOSE

Time-resolved 3D rotational angiography (4D-DSA) has been used to demonstrate details of the angioarchitecture of AVM, whereas it has rarely been used to describe features of dural AVF. In this exploratory study, we analyzed dural AVFs with a novel 4D software prototype, developed and provided by Siemens, to determine whether identification of the location of the fistulous point, grading, and treatment planning were feasible.

MATERIALS AND METHODS

4D-DSA volumes were calculated from existing 3D rotational angiography data sets of patients with dural AVFs. The 4D-DSA volumes were displayed in a virtual DSA mode and MPR or MIP in 3 orthogonal planes and compared with 2D-DSA by 2 experienced neuroradiologists. Fusions with unenhanced CT or MR images were used to improve visualization of adjacent anatomic structures.

RESULTS

Comparison with 2D-DSA showed that evaluation of the fistulous point and grading according to the classification of Borden, Cognard, or Barrow was feasible in 26 of 27 cases. In 8 of 27 cases, 4D-DSA was considered advantageous for determining the fistulous point and the course of the draining vein in the dural AVF with cortical venous drainage, especially in the frontoethmoidal and frontoparietal regions. In 6 cases, the display of angioarchitecture was considered inferior to that of 2D-DSA due to motion artifacts, suboptimal selection of the injected vessel, and lack of temporal resolution.

CONCLUSIONS

Detailed analysis of dural AVFs according to the standardized display of 4D-DSA volumes was feasible and helpful in understanding the angioarchitecture in selected cases. Further improvement and validation of the 4D software should solidify the complementary role of 4D-DSA to conventional 2D-DSA series.

Combined cone-beam CT imaging and microsurgical dissection of cadaver specimens to study cerebral venous anatomy: a technical note

PURPOSE

Cadaver dissections and X-ray based 3D angiography are considered gold standards for studying neurovascular anatomy. We sought to develop a model that utilize the combination of both these techniques to improve current tools for anatomical research, teaching and preoperative surgical planning, particularly addressing the venous system of the brain.

MATERIALS AND METHODS

Seven ethanol-fixed human cadaveric heads and one arm were injected with a latex-barium mixture into the internal jugular veins and the brachial artery. After the ethanol-based fixation, specimens were scanned by high-resolution cone-beam CT and images were post-processed on a 3D-workstation. Subsequent, microsurgical dissections were performed by an experienced neurosurgeon and venous anatomy was compared with relevant 3D venograms.

RESULTS

Latex-barium mixtures resulted in a homogenous cast with filling of the cerebral venous structures down to 150 μm in diameter. The ethanol-based preparation of the cadaveric brains allowed for near-realistic microsurgical maneuverability during dissection. The model improves assessment of the venous system for anatomical education and hands-on surgical training.

CONCLUSION

To our knowledge we describe the first preparation method which combines near-realistic microsurgical dissection of human heads with high-resolution 3D imaging of the cerebral venous system in the same specimens.

Feasibility of 4D HDR brachytherapy source tracking using x-ray tomosynthesis: Monte Carlo investigation

PURPOSE

High dose rate (HDR) brachytherapy rapidly delivers dose to targets with steep dose gradients. This treatment method must adhere to prescribed treatment plans with high spatiotemporal accuracy and precision, as failure to do so may degrade clinical outcomes. One approach to achieving this goal is to develop imaging techniques to track HDR sources in vivo in reference to surrounding anatomy. This work investigates the feasibility of using an isocentric C-arm x-ray imager and tomosynthesis methods to track Ir-192 HDR brachytherapy sources in vivo over time (4D).

METHODS

A tomosynthesis imaging workflow was proposed and its achievable source detectability, localization accuracy, and spatiotemporal resolution were investigated in silico. An anthropomorphic female XCAT phantom was modified to include a vaginal cylinder applicator and Ir-192 HDR source (0.5 × 0.5 × 5.0 mm3 ), and the workflow was carried out using the MC-GPU Monte Carlo image simulation platform. Source detectability was characterized using the reconstructed source signal-difference-to-noise-ratio (SDNR), localization accuracy by the absolute 3D error in its measured centroid location, and spatiotemporal resolution by the full-width-at-half-maximum (FWHM) of line profiles through the source in each spatial dimension considering a maximum C-arm angular velocity of 30° per second. The dependence of these parameters on acquisition angular range (θtot = 0°-90°), number of views, angular increment between views (Δθ = 0°-15°), and volumetric constraints imposed in reconstruction was evaluated. Organ voxel doses were tallied to derive the workflow's attributable effective dose.

RESULTS

The HDR source was readily detected and its centroid was accurately localized with the proposed workflow and method (SDNR: 10-40, 3D error: 0-0.144 mm). Tradeoffs were demonstrated for various combinations of image acquisition parameters; namely, increasing the tomosynthesis acquisition angular range improved resolution in the depth-encoded direction, for example from 2.5 mm to 1.2 mm between θtot = 30o and θtot = 90o , at the cost of increasing acquisition time from 1 to 3 s. The best-performing acquisition parameters (θtot = 90o , Δθ = 1°) yielded no centroid localization error, and achieved submillimeter source resolution (0.57 × 1.21 × 5.04 mm3 apparent source dimensions, FWHM). The total effective dose for the workflow was 263 µSv for its required pre-treatment imaging component and 7.59 µSv per mid-treatment acquisition thereafter, which is comparable to common diagnostic radiology exams.

CONCLUSIONS

A system and method for tracking HDR brachytherapy sources in vivo using C-arm tomosynthesis was proposed and its performance investigated in silico. Tradeoffs in source conspicuity, localization accuracy, spatiotemporal resolution, and dose were determined. The results suggest this approach is feasible for localizing an Ir-192 HDR source in vivo with submillimeter spatial resolution, 1-3 second temporal resolution and minimal additional dose burden.

Estimating Pulsatile Flow Velocity using Four-Dimensional Digital Subtraction Angiography. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2023;2023:1-4. [PMID: 38082691 DOI: 10.1109/embc40787.2023.10340843] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Digital subtraction angiography (DSA) is a X-ray based imaging modality for interventional procedures, and remains criterion standard for diagnosing vascular diseases. The imaging protocol of DSA involves administration of a foreign contrast medium into the blood vessel that opacifies the vasculature during the imaging. Using two-dimensional (2D) DSA with high temporal resolution, it was recently demonstrated that the pulsatile velocity can be estimated by evaluating the temporal and spatial variations of the contrast medium distributions in the blood vessel. In this paper, we evaluate the feasibility to estimate the pulsatile flow velocity using the four-dimensional (4D) DSA. To overcome the noise and artefacts of 4D-DSA data, a empirical mode decomposition plus autocorrelation based method is proposed to estimate the pulsatile velocities, and the pulsatile velocities estimated using 2D-DSA are used as reference for comparison.Clinical Relevance-4D-DSA encompasses both structural and temporal information; it theoretically reduces the need of multiple scans, hence reducing the radiation doses. The estimated pulsatile flow velocities open up a new parameter for hemodynamic studies and potential for real-time diagnostic and therapeutic monitoring during interventional procedures.

Importance of Anatomical Variation of the Hepatic Artery for Complicated Liver and Pancreatic Surgeries: A Review Emphasizing Origin and Branching

Knowledge of anatomical variations of the hepatic artery from its origin to intrahepatic segmentation is of utmost importance for planning upper abdominal surgeries including liver transplantation, pancreatoduodenectomy, and biliary reconstruction. The origin and branching pattern of the hepatic artery was thoroughly described by the classification of Michels and Hiatt. Some rare variations of the hepatic artery were classified by Kobayashi and Koops. By the use of the multidetector computed tomography (MDCT) technique, the branching pattern of the hepatic artery can be visualized quite accurately. Unawareness of these arterial variations may lead to intraoperative injuries such as necrosis, abscess, and failure of the liver and pancreas. The origin and course of the aberrant hepatic arteries are crucial in the surgical planning of carcinoma of the head of the pancreas and hepatobiliary surgeries. In liver transplant surgeries, to minimize intraoperative bleeding complications and postoperative thrombosis, exact anatomy of the branching of the hepatic artery, its variations and intrahepatic course is of utmost importance. This review discusses variations in the anatomy of the hepatic artery from its origin to branching by the use of advanced imaging techniques and its effect on the liver, pancreatic, biliary and gastric surgeries.

Imaging of intracranial aneurysms in animals: a systematic review of modalities

Intracranial aneurysm (IA) animal models are paramount to study IA pathophysiology and to test new endovascular treatments. A number of in vivo imaging modalities are available to characterize IAs at different stages of development in these animal models. This review describes existing in vivo imaging techniques used so far to visualize IAs in animal models. We systematically searched for studies containing in vivo imaging of induced IAs in animal models in PubMed and SPIE Digital library databases between 1 January 1945 and 13 July 2022. A total of 170 studies were retrieved and reviewed in detail, and information on the IA animal model, the objective of the study, and the imaging modality used was collected. A variety of methods to surgically construct or endogenously induce IAs in animals were identified, and 88% of the reviewed studies used surgical methods. The large majority of IA imaging in animals was performed for 4 reasons: basic research for IA models, testing of new IA treatment modalities, research on IA in vivo imaging of IAs, and research on IA pathophysiology. Six different imaging techniques were identified: conventional catheter angiography, computed tomography angiography, magnetic resonance angiography, hemodynamic imaging, optical coherence tomography, and fluorescence imaging. This review presents and discusses the advantages and disadvantages of all in vivo IA imaging techniques used in animal models to help future IA studies finding the most appropriate IA imaging modality and animal model to answer their research question.

Four-dimensional digital subtraction angiography to assess cerebral arteriovenous malformations

BACKGROUND AND PURPOSE

The performance of a novel prototype four-dimensional (4D) digital subtraction angiography (DSA) for cerebral arteriovenous malformation (AVM) diagnosis was evaluated and compared with that of two-dimensional (2D) and three-dimensional (3D) DSA.

METHODS

In this retrospective study, 37 consecutive cerebral AVM patients were included. The standard diagnostic results were concluded from the 2D and 3D DSA. Two 4D DSA volumes were reconstructed for each patient by a commercial and a prototype software, then evaluated by two independent experienced neurosurgeons, who were blinded to the diagnosis and treatment process. The evaluation results were compared with the diagnostic results on Spetzler-Martin (SM) Grading Scale, number of feeding arteries, number of draining veins, and intranidal aneurysms.

RESULTS

Complete agreement was achieved between 4D DSA and 2D and 3D DSA in SM Grading Scale and intracranial aneurysm identification (agreement coefficient: 1) for both reviewers. The agreement coefficients were .888 and .917 for both reviewers in feeding artery number determination using 4D DSA product and 4D DSA prototype, respectively. The agreement coefficients in draining vein number determination were all larger than .94 for both reviewers using both 4D DSA volumes.

CONCLUSIONS

The performance of this prototype 4D DSA in cerebral AVMs diagnosis was largely equivalent to that of 2D and 3D DSA combination. Four-dimensional DSA can be regarded as a very good complement for 2D DSA in cerebral AVM diagnosis.

Validation of SART 3.5D algorithm for cerebrovascular dynamics and artery versus vein classification in presurgical 3D digital subtraction angiographies

Classification of arteries and veins in cerebral angiograms can increase the safety of neurosurgical procedures, such as StereoElectroEncephaloGraphy, and aid the diagnosis of vascular pathologies, as arterovenous malformations. We propose a new method for vessel classification using the contrast medium dynamics in rotational digital subtraction angiography (DSA). After 3D DSA and angiogram segmentation, contrast enhanced projections are processed to suppress soft tissue and bone structures attenuation effect and further enhance the CM flow. For each voxel labelled as vessel, a time intensity curve (TIC) is obtained as a linear combination of temporal basis functions whose weights are addressed by simultaneous algebraic reconstruction technique (SART 3.5D), expanded to include dynamics. Each TIC is classified by comparing the areas under the curve in the arterial and venous phases. Clustering is applied to optimize the classification thresholds. On a dataset of 60 patients, a median value of sensitivity (90%), specificity (91%), and accuracy (92%) were obtained with respect to annotated arterial and venous voxels up to branching order 4–5. Qualitative results are also presented about CM arrival time mapping and its distribution in arteries and veins respectively. In conclusion, this study shows a valuable impact, at no protocol extra-cost or invasiveness, concerning surgical planning related to the enhancement of arteries as major organs at risk. Also, it opens a new scope on the pathophysiology of cerebrovascular dynamics and its anatomical relationships.

Using temporal and structural data to reconstruct 3D cerebral vasculature from a pair of 2D digital subtraction angiography sequences

PURPOSE

The purpose of this work is to present a new method for reconstructing patient-specific three-dimensional (3D) vasculature of the brain from a pair of digital subtraction angiography (DSA) image sequences from different viewpoints, e.g., from bi-plane angiography. Our long-term goal is to provide high resolution visualization of 3D vasculature with dynamic flow of contrast agent from limited data that is readily available during surgical procedures. The proposed method is the second of a three-stage process composed of 1) augmenting vessel segmentation with vessel radii and timing of the arrival of a bolus of contrast agent, 2) reconstructing a volumetric representation of the augmented vessel data from the augmented 2D segmentations, and 3) generating a 3D model of vessels and flow of contrast agent from the volumetric reconstruction. Unlike previous methods, which are either limited to relatively simple vessel structures or rely on multiple views and/or prior models of the vasculature, our method requires only a single pair of 2D DSA sequences taken from different view directions.

METHODS

We developed a new mathematical algorithm that augments vessel centerlines with vessel radii and bolus arrival times derived directly from the 2D DSA sequences to constrain the 3D reconstruction. We validated this method on digital phantoms derived from clinical data and from fractal models of branching tree structures.

RESULTS

In standard reconstruction methods, reconstruction by projection of two views into 3D space results in 'ghosting' artifacts, i.e., false 3D structure that occurs where vessels or vessel segments overlap in the 2D images. For the complex vascular of the brain, this ghosting is severe and is a major hurdle for methods that attempt to generate 3D structure from 2D images. We show that our approach reduces ghosting by up to 99% in digital phantoms derived from clinical data.

CONCLUSION

Our dramatic reduction in ghosting artifacts in 3D reconstructions from a pair of 2D image sequences is an important step towards generating high resolution 3D vasculature with dynamic flow information from a single DSA sequence acquired using bi-plane angiography.

Extraction of patient-specific boundary conditions from 4D-DSA and their influence on CFD simulations of cerebral aneurysms

We developed a new technique for extracting patient-specific inflow conditions, such as the pulse cycle duration and blood flow velocity, from four-dimensional digital subtraction angiography images and experimentally examined its validity. The maximum error between the values extracted by the technique and measured values was 14.3%. We performed blood flow simulations and calculated representative haemodynamic parameters. The maximum differences between the parameters obtained using general and patient-specific inflow conditions were approximately 400%, 150%, and 50% for the velocity, normalised wall shear stress, and pressure loss coefficient, respectively. These results indicate that patient-specific conditions are critical for accurately reproducing aneurysmal haemodynamics.

Imaging Modalities for Intracranial Aneurysm: More Than Meets the Eye

Intracranial aneurysms (IA) are often asymptomatic and have a prevalence of 3 to 5% in the adult population. The risk of IA rupture is low, however when it occurs half of the patients dies from subarachnoid hemorrhage (SAH). To avoid this fatal evolution, the main treatment is an invasive surgical procedure, which is considered to be at high risk of rupture. This risk score of IA rupture is evaluated mainly according to its size and location. Therefore, angiography and anatomic imaging of the intracranial aneurysm are crucial for its diagnosis. Moreover, it has become obvious in recent years that several other factors are implied in this complication, such as the blood flow complexity or inflammation. These recent findings lead to the development of new IA imaging tools such as vessel wall imaging, 4D-MRI, or molecular MRI to visualize inflammation at the site of IA in human and animal models. In this review, we will summarize IA imaging techniques used for the patients and those currently in development.

Flat-panel conebeam CT in the clinic: history and current state

Research into conebeam CT concepts began as soon as the first clinical single-slice CT scanner was conceived. Early implementations of conebeam CT in the 1980s focused on high-contrast applications where concurrent high resolution ( < 200    μ m ), for visualization of small contrast-filled vessels, bones, or teeth, was an imaging requirement that could not be met by the contemporaneous CT scanners. However, the use of nonlinear imagers, e.g., x-ray image intensifiers, limited the clinical utility of the earliest diagnostic conebeam CT systems. The development of consumer-electronics large-area displays provided a technical foundation that was leveraged in the 1990s to first produce large-area digital x-ray detectors for use in radiography and then compact flat panels suitable for high-resolution and high-frame-rate conebeam CT. In this review, we show the concurrent evolution of digital flat panel (DFP) technology and clinical conebeam CT. We give a brief summary of conebeam CT reconstruction, followed by a brief review of the correction approaches for DFP-specific artifacts. The historical development and current status of flat-panel conebeam CT in four clinical areas-breast, fixed C-arm, image-guided radiation therapy, and extremity/head-is presented. Advances in DFP technology over the past two decades have led to improved visualization of high-contrast, high-resolution clinical tasks, and image quality now approaches the soft-tissue contrast resolution that is the standard in clinical CT. Future technical developments in DFPs will enable an even broader range of clinical applications; research in the arena of flat-panel CT shows no signs of slowing down.

Uterine Arteriovenous Malformation: A Pictorial Review of Diagnosis and Management

Uterine arteriovenous malformation (UAVM) is a rare condition and is classified as either congenital or acquired UAVM. Patients with UAVMs usually experience miscarriages or recurrent menorrhagia. Ultrasound is used for the initial estimation of UAVMs. Computed tomography and magnetic resonance imaging are noninvasive and valuable methods that provide good compatibility with digital subtraction angiography to support the diagnosis and treatment of UAVM. Timely diagnosis is crucial to provide appropriate treatment for alleviating complications. This article presents a pictorial and literature review of the current evidence of the diagnosis and management of UAVM.

Estimating blood flow velocity using time-resolved 3D angiography and a derived physical law of contrast media

OBJECTIVE

Four-dimensional (4D) digital subtraction angiography (DSA) offers a method for evaluating hemodynamics. It is, however, unclear how the delivered contrast medium interacts with the physiological blood flow, and how hemodynamic information may be inferred from the mixture of the contrast medium and blood. In this study, we present a theoretical explanation of contrast dynamics, and an accompanying algorithm for estimating blood flow velocity.

APPROACH

We retrospectively recruited 23 patients who underwent both 4D DSA and magnetic resonance (MR) phase-contrast imaging. The 4D DSA-reconstructed contrast dynamics were first studied for the internal carotid arteries. Using physical laws governing fluid motion within a curved tube, we showed that the reconstructed contrast dynamics obeyed a simple advection equation. We then proposed an algorithm for estimating the contrast dynamics using angiographic data, and subsequently estimated the axial blood flow velocity using an advection equation.

MAIN RESULTS

The estimated velocities were compared using three techniques: the Fourier technique, Lin's method, and MR phase contrast. Testing with noise-corrupted artificial data showed that the proposed algorithm was noise resistant. The velocities of 23 patients computed by 4D DSA using the proposed algorithm showed a moderate correlation with the MR phase contrast (r = 0.61), and good correlations with the other two techniques (r = 0.75 and r = 0.72).

SIGNIFICANCE

The proposed algorithm and has been applied to blood vessel segments with poor signal-to-noise ratios and axial lengths of less than 3 cm, and has a physical basis for computing axial flow velocities using an advection equation. The results of the proposed algorithm are consistent with existing methods.

Four-dimensional digital subtraction angiography for exploration of intraosseous arteriovenous fistula in the sphenoid bone

Background

Intraosseous arteriovenous fistula (AVF) is a rare clinical entity that typically presents with symptoms from their effect on surrounding structures. Here, we report a case of intraosseous AVF in the sphenoid bone that presented with bilateral abducens palsy.

Case Description

A previously healthy man presented with tinnitus for 1 month, and initial imaging suspected dural AVF of the cavernous sinus. Four-dimensional digital subtraction angiography (4D-DSA) imaging and a three-dimensional (3D) fused image from the bilateral external carotid arteries revealed that the shunt was in a large venous pouch within the sphenoid bone that was treated through transvenous coil embolization. His symptoms improved the day after surgery.

Conclusion

This is a case presentation of intraosseous AVF in the sphenoid bone and highlights the importance of 4D-DSA and 3D fused images for planning the treatment strategy.

4D-DSA: Development and Current Neurovascular Applications

Originally described by Davis et al in 2013, 4D-Digital Subtraction Angiography (4D-DSA) has developed into a commercially available application of DSA in the angiography suite. 4D-DSA provides the user with 3D time-resolved images, allowing observation of a contrast bolus at any desired viewing angle through the vasculature and at any time point during the acquisition (any view at any time). 4D-DSA mitigates some limitations that are intrinsic to both 2D- and 3D-DSA images. The clinical applications for 4D-DSA include evaluations of AVMs and AVFs, intracranial aneurysms, and atherosclerotic occlusive disease. Recent advances in blood flow quantification using 4D-DSA indicate that these data provide both the velocity and geometric information necessary for the quantification of blood flow. In this review, we will discuss the development, acquisition, reconstruction, and current neurovascular applications of 4D-DSA volumes.

4D Digital Subtraction Angiography for the Temporal Flow Visualization of Intracranial Aneurysms and Vascular Malformations

PURPOSE

To assess the benefit and radiation dose of four-dimensional (4D) digital subtraction angiography (DSA) - a time resolved three-dimensional (3D) DSA application - to evaluate the flow and architecture of aneurysms and vascular malformations.

METHODS

All patients with cerebrovascular disease were considered who underwent 4D-DSA at our institution between January 2015 and February 2016. For the aneurysm patients, we evaluated the image quality in terms of the visualization of contrast flow in the aneurysm on a 3-point scale (excellent, fair and poor). Interrater agreement between two raters was estimated using Cohen's Kappa statistics. For the shunt disease patients, the additional information obtained from the 4D-DSA was described for each disease. The median radiation dose and volume of contrast medium required for the acquisitions were estimated.

RESULTS

A total of 173 patients underwent 4D-DSA; 126 intracranial aneurysms, 10 arteriovenous malformations (AVM), 15 dural arteriovenous fistula (dAVF) and 22 other diseases. For aneurysm patients, excellent and fair visualization of the intra-aneurysmal flow was observed in 27.7%, 72.3%, and excellent (κ = 0.9) agreement between the raters was found. For AVM and dAVF patients, 4D-DSA clarified the complex vasculature by viewing the discrete time phase of contrast filling. Median radiation dose for intracranial lesions was 79.6 mGy for 6s 4D-DSA, and 175 mGy for 12s 4D-DSA. The median amount of contrast medium used was 18.0 ml for 6s 4D-DSA and 21.0 ml for 12s 4D-DSA.

CONCLUSIONS

4D-DSA provided additional information regarding intra-aneurysmal flow and contributed to detect different component of nidus or shunt points.

Virtual 2D angiography from four-dimensional digital subtraction angiography (4D-DSA): A feasibility study

BACKGROUND

Digital subtraction angiography (DSA) remains the gold standard for angiographic evaluation of cerebrovascular pathology, however, multiple acquisitions requiring additional time and radiation are often needed. In contrast, 3D-DSA provides volumetric information from a single injection but neglects temporal information. Four-dimensional-DSA (4D-DSA) combines temporal information of 2D-DSA with volumetric information of 3D-DSA to provide time-resolved tomographic 3D reconstructions, potentially reducing procedure time and radiation. This work evaluates the diagnostic quality of virtual single-frame 4D-DSA relative to 2D-DSA images by assessing clinicians' ability to evaluate cerebrovascular pathology.

METHODS

Single-frame images of four projections from 4D-DSA and their corresponding 2D-DSA images (n = 15) were rated by two neurointerventional radiologists. Images were graded based on diagnostic quality (0 = non-diagnostic, 1 = poor, 2 = acceptable, 3 = good). Dose area product (DAP) for each case was recorded for all 2D-DSA, 4D-DSA acquisitions, and the overall procedure.

RESULTS

The mean diagnostic quality of all four 4D-DSA projections from both raters was 1.75 while the mean of 2D-DSA projections was 2.8. Student's t-test revealed significant difference in diagnostic quality between 4D-DSA and 2D-DSA at all four projections (p < 0.001). On average 4D-DSA acquisitions accounted for 30% dose compared to the overall average aggregated dose per procedure.

CONCLUSIONS

The difference in image quality between virtual single-frame 4D-DSA and their respective 2D-DSA images is statistically significant. Furthermore, 4D-DSA acquisitions require less radiation dose than conventional procedures with 2D-DSA acquisitions.

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.

Four-dimensional digital subtraction angiography for exploration of spinal cord vascular malformations: preliminary experience

BACKGROUND

The precise understanding of the angioarchitecture of spinal vascular malformations (SVMs) is often difficult to reach with conventional digital subtraction angiography (DSA). The purpose of our study was to evaluate the potential of four-dimensional DSA (4D-DSA) (Siemens Healthcare) in the exploration of SVMs.

METHODS

We retrospectively studied all patients who underwent spinal DSA, including 4D-DSA acquisition, from July 2018 to June 2019 at a single institution. All spinal DSA acquisitions were performed under general anesthesia. 4D-DSA acquisitions were acquired with the protocol '12 s DSA Dyna4D Neuro'. 12 mL of iodixanol 320 mg iodine/mL were injected via a 5 F catheter (1 mL/s during the 12 s 4D-DSA acquisition). Inter-rater (three independent reviewers) and intermodality agreements were assessed.

RESULTS

Nine consecutive patients (six men, three women, mean age 55.3±19.8 years) with 10 SVMs (spinal dural arteriovenous fistulas n=3, spinal epidural arteriovenous fistulas n=2, spinal pial arteriovenous fistulas n=2, and spinal arteriovenous malformations n=2; one patient had two synchronous pial fistulas) had spinal DSA, including 4D-DSA acquisition. Inter-rater agreement was good and moderate for the venous drainage pattern and the SVM subtype, respectively. In 9 of 10 cases, the quality of the acquisition was graded as good. Satisfactory concordance between 4D-DSA and the selective microcatheterization was observed in 90% of cases for the location of the shunt point.

CONCLUSION

4D-DSA acquisition may be helpful for a better understanding of the angioarchitecture of SVMs. Larger series are warranted to confirm these preliminary results.

Optimizing the Quality of 4D-DSA Temporal Information

BACKGROUND AND PURPOSE

Quantification of blood flow using a 4D-DSA would be useful in the diagnosis and treatment of cerebrovascular diseases. A protocol optimizing identification of density variations in the time-density curves of a 4D-DSA has not been defined. Our purpose was to determine the contrast injection protocol most likely to result in the optimal pulsatility signal strength.

MATERIALS AND METHODS

Two 3D-printed patient-specific models were used and connected to a pulsatile pump and flow system, which delivered 250-260 mL/min to the model. Contrast medium (Isovue, 370 mg I/mL, 75% dilution) was injected through a 6F catheter positioned upstream from the inlet of the model. 4D-DSA acquisitions were performed for the following injection rates: 1.5, 2.0, 2.5, 3.0 and 3.5 mL/s for 8 seconds. To determine pulsatility, we analyzed the time-density curve at the inlets using the oscillation amplitude and a previously described numeric metric, the sideband ratio. Vascular geometry from 4D-DSA reconstructions was compared with ground truth and micro-CT measurements of the model. Dimensionless numbers that characterize hemodynamics, Reynolds and Craya-Curtet, were calculated for each injection rate.

RESULTS

The strongest pulsatility signal occurred with the 2.5 mL/s injections. The largest oscillation amplitudes were found with 2.0- and 2.5-mL/s injections. Geometric accuracy was best preserved with injection rates of >1.5 mL/s.

CONCLUSIONS

An injection rate of 2.5 mL/s provided the strongest pulsatility signal in the 4D-DSA time-density curve. Geometric accuracy was best preserved with injection rates above 1.5 mL/s. These results may be useful in future in vivo studies of blood flow quantification.

Quantitative and Qualitative Comparison of 4D-DSA with 3D-DSA Using Computational Fluid Dynamics Simulations in Cerebral Aneurysms

BACKGROUND AND PURPOSE

4D-DSA allows time-resolved 3D imaging of the cerebral vasculature. The aim of our study was to evaluate this method in comparison with the current criterion standard 3D-DSA by qualitative and quantitative means using computational fluid dynamics.

MATERIALS AND METHODS

3D- and 4D-DSA datasets were acquired in patients with cerebral aneurysms. Computational fluid dynamics analysis was performed for all datasets. Using computational fluid dynamics, we compared 4D-DSA with 3D-DSA in terms of both aneurysmal geometry (quantitative: maximum diameter, ostium size [OZ1/2], volume) and hemodynamic parameters (qualitative: flow stability, flow complexity, inflow concentration; quantitative: average/maximum wall shear stress, impingement zone, low-stress zone, intra-aneurysmal pressure, and flow velocity). Qualitative parameters were descriptively analyzed. Correlation coefficients (r, P value) were calculated for quantitative parameters.

RESULTS

3D- and 4D-DSA datasets of 10 cerebral aneurysms in 10 patients were postprocessed. Evaluation of aneurysmal geometry with 4D-DSA (r _{maximum diameter} = 0.98, P _{maximum diameter} <.001; r _OZ1/OZ2 = 0.98/0.86, P _OZ1/OZ2 < .001/.002; r _volume = 0.98, P _volume <.001) correlated highly with 3D-DSA. Evaluation of qualitative hemodynamic parameters (flow stability, flow complexity, inflow concentration) did show complete accordance, and evaluation of quantitative hemodynamic parameters (r _{average/maximum wall shear stress diastole} = 0.92/0.88, P _{average/maximum wall shear stress diastole} < .001/.001; r _{average/maximum wall shear stress systole} = 0.94/0.93, P _{average/maximum wall shear stress systole} < .001/.001; r _{impingement zone} = 0.96, P _{impingement zone} < .001; r _{low-stress zone} = 1.00, P _{low-stress zone} = .01; r _{pressure diastole} = 0.84, P _{pressure diastole} = .002; r _{pressure systole} = 0.9, P _{pressure systole} < .001; r _{flow velocity diastole} = 0.95, P _{flow velocity diastole} < .001; r _{flow velocity systole} = 0.93, P _{flow velocity systole} < .001) did show nearly complete accordance between 4D- and 3D-DSA.

CONCLUSIONS

Despite a different injection protocol, 4D-DSA is a reliable basis for computational fluid dynamics analysis of the intracranial vasculature and provides equivalent visualization of aneurysm geometry compared with 3D-DSA.

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.

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

BACKGROUND AND PURPOSE

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

MATERIALS AND METHODS

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

RESULTS

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

CONCLUSIONS

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

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.

A new time-resolved 3D angiographic technique (4D DSA): Description, and assessment of its reliability in Spetzler-Martin grading of cerebral arteriovenous malformations

BACKGROUND AND PURPOSE

The Spetzler and Martin (SM) cerebral arteriovenous malformation (AVM) classification is a widely used 5-tier classification. This common language allows specialists to exchange about AVMs and must be reliably characterized by the imaging methods. We presented an agreement study on a new method of digital subtracted 3D rotational angiography resolved in time (four-dimensional DSA: 4D DSA) compared to the gold standard (two-dimensional DSA: 2D DSA) in AVM grading using the SM classification.

METHODS

Ten patients with AVMs were included during one year, they had an angiographic exploration with both 4D DSA and 2D DSA. Three readers assessed the SM classification. One reader conducted a second reading. The inter-, intra-observer and intermodality agreements were calculated by Kappas. Dose to patient was reported.

RESULTS

Considering the SM grade, the inter-observer agreement between 4D DSA and 2D DSA was equivalent (κ=0.45 and 0.46), and calculated as substantial κ=0.76 between the 2 methods. The agreement between 4D DSA and 2D DSA was calculated as moderate κ=0.46 assessing the size of the nidus, slight κ=0.18 analyzing the drainage and almost perfect κ=0.95 depicting the localization. 4D DSA performed during a standard initial angiographic assessment of AVM represented approximately 6% of the total dose.

CONCLUSION

The addition of this new technique 4D DSA could be performed regularly in addition to the 2D DSA if available, to assess SM grading, with an acceptable exposure to ionizing radiation.

4D DSA for Dynamic Visualization of Cerebral Vasculature: A Single-Center Experience in 26 Cases

BACKGROUND AND PURPOSE

4D DSA allows acquisition of time-resolved 3D reconstructions of cerebral vessels by using C-arm conebeam CT systems. The aim of our study was to evaluate this new method by qualitative and quantitative means.

MATERIALS AND METHODS

2D and 4D DSA datasets were acquired in patients presenting with AVMs, dural arteriovenous fistulas, and cerebral aneurysms. 4D DSA was compared with 2D DSA in a consensus reading of qualitative and quantitative parameters of AVMs (eg, location, feeder, associated aneurysms, nidus size, drainage, Martin-Spetzler Score), dural arteriovenous fistulas (eg, fistulous point, main feeder, diameter of the main feeder, drainage), and cerebral aneurysms (location, neck configuration, aneurysmal size). Identifiability of perforators and diameters of the injection vessel (ICA, vertebral artery) were analyzed in 2D and 4D DSA. Correlation coefficients and a paired t test were calculated for quantitative parameters. The effective patient dose of the 4D DSA protocol was evaluated with an anthropomorphic phantom.

RESULTS

In 26 patients, datasets were acquired successfully (AVM = 10, cerebral aneurysm = 10, dural arteriovenous fistula = 6). Qualitative and quantitative evaluations of 4D DSA in AVMs (nidus size: r = 0.99, P = .001), dural arteriovenous fistulas (diameter of the main feeder: r = 0.954, P = .03), and cerebral aneurysms (aneurysmal size: r = 1, P = .001) revealed nearly complete accordance with 2D DSA. Perforators were comparably visualized with 4D DSA. Measurement of the diameter of the injection vessel in 4D DSA was equivalent to that in 2D DSA (P = .039). The effective patient dose of 4D DSA was 1.2 mSv.

CONCLUSIONS

4D DSA is feasible for imaging of AVMs, dural arteriovenous fistulas, and cerebral aneurysms. 4D DSA offers reliable visualization of the cerebral vasculature and may improve the understanding and treatment of AVMs and dural arteriovenous fistulas. The number of 2D DSA acquisitions required for an examination may be reduced through 4D DSA.

Application of Time-Resolved 3D Digital Subtraction Angiography to Plan Cerebral Arteriovenous Malformation Radiosurgery

BACKGROUND AND PURPOSE

Time-resolved 3D-DSA (4D-DSA) enables viewing vasculature from any desired angle and time frame. We investigated whether these advantages may facilitate treatment planning and the feasibility of using 4D-DSA as a single imaging technique in AVM/dural arteriovenous fistula radiosurgery.

MATERIALS AND METHODS

Twenty consecutive patients (8 dural arteriovenous fistulas and 12 AVMs; 13 men and 7 women; mean age, 45 years; range, 18-64 years) who were scheduled for gamma knife radiosurgery were recruited (November 2014 to October 2015). An optimal volume of reconstructed time-resolved 3D volumes that defines the AVM nidus/dural arteriovenous fistula was sliced into 2D-CT-like images. The original radiosurgery treatment plan was overlaid retrospectively. The registration errors of stereotactic 4D-DSA were compared with those of integrated stereotactic imaging. AVM/dural arteriovenous fistula volumes were contoured, and disjoint and conjoint components were identified. The Wilcoxon signed rank test and the Wilcoxon rank sum test were adopted to evaluate registration errors and contoured volumes of stereotactic 4D-DSA and integration of stereotactic MR imaging and stereotactic 2D-DSA.

RESULTS

Sixteen of 20 patients were successfully registered in Advanced Leksell GammaPlan Program. The registration error of stereotactic 4D-DSA was smaller than that of integrated stereotactic imaging (P = .0009). The contoured AVM volume of 4D-DSA was smaller than that contoured on the integration of MR imaging and 2D-DSA, while major inconsistencies existed in cases of dural arteriovenous fistula (P = .042 and 0.039, respectively, for measurements conducted by 2 authors).

CONCLUSIONS

Implementation of stereotactic 4D-DSA data for gamma knife radiosurgery for brain AVM/dural arteriovenous fistula is feasible. The ability of 4D-DSA to demonstrate vascular morphology and hemodynamics in 4 dimensions potentially reduces the target volumes of irradiation in vascular radiosurgery.

Time-resolved 3D Rotational Angiography (4D DSA) of the Lenticulostriate Arteries: Display of Normal Anatomic Variants and Collaterals in Cases with Chronic Obstruction of the MCA

PURPOSE

The lenticulostriate arteries (LSA) and other perforators may play a role for collateral supply in cases with ischemia due to stenosis or occlusions of the middle cerebral artery (MCA). Purpose of this case series was to evaluate the potential of time-resolved 3D rotational angiography data sets (4D DSA) for detailed visualization of anatomic variants of LSA feeders and for display of local collaterals involving the LSA in cases with chronic MCA obstruction.

METHODS

Multiplanar and volume rendering reconstructions of 4D DSA data were computed in addition to standard postprocessing in 24 patients who had indications for 3D rotational angiography (3DRA) of the internal carotid artery (ICA) without pathologies of the ICA, middle cerebral artery (MCA) and anterior cerebral artery (ACA) main stems (n = 18) or with stenosis or chronic occlusion of the MCA (n = 6). For acquisition of 3DRA, we used a modified digital subtraction angiography (DSA) image acquisition protocol with an extended rotation angle of 260° and a prolonged scan time of 12 s on a Siemens Axiom Artis Zee biplane neuroangiography equipment. The 4D reconstructions of existing 3DRA data were computed on a dedicated workstation. Origin and course of LSA and other perforators were analyzed according to coronal multiplanar reconstructions (MPRs) with slice thicknesses between 6 and 28 mm.

RESULTS

In all cases 4D reconstructions of the LSA were technically feasible and evaluable. As expected, origin and course of LSA showed a wide range of variations: The most common pattern was a common trunk dividing into multiple ascending branches originating from the proximal M1 (n = 5) or the proximal A1 segment (n = 4). Alternatively, 8 patients showed several individual branches that directly originated from the proximal M1 segment of the MCA and occasionally from the A1 segment of the ACA. In patients with M1 stenosis or occlusion, 4 out of 6 cases had local collaterals with involvement of proximal LSA trunks and a network parallel to the obstructed vessel segment. The 4D reconstructions were found to be equivalent (n = 16) or superior to 3D reconstructions (n = 8).

CONCLUSION

The 4D DSA reconstructions provide a reliable display of normal LSA variants and connections to local collateral networks in cases with chronic MCA obstruction. The possibility to select a correct angiographic phase is advantageous compared to 3D DSA.

Comparison of the Diagnostic Utility of 4D-DSA with Conventional 2D- and 3D-DSA in the Diagnosis of Cerebrovascular Abnormalities

BACKGROUND AND PURPOSE

4D-DSA is a time-resolved technique that allows viewing of a contrast bolus at any time and from any desired viewing angle. Our hypothesis was that the information content in a 4D-DSA reconstruction was essentially equivalent to that in a combination of 2D acquisitions and a 3D-DSA reconstruction.

MATERIALS AND METHODS

Twenty-six consecutive patients who had both 2D- and 3D-DSA acquisitions were included in the study. The angiography report was used to obtain diagnoses and characteristics of abnormalities. Diagnoses included AVM/AVFs, aneurysms, stenosis, and healthy individuals. 4D-DSA reconstructions were independently reviewed by 3 experienced observers who had no part in the clinical care. Using an electronic evaluation form, these observers recorded their assessments based only on the 4D reconstructions. The clinical evaluations were then compared with the 4D evaluations for diagnosis and lesion characteristics.

RESULTS

Results showed both interrater and interclass agreements (κ = 0.813 and 0.858). Comparing the 4D diagnosis with the clinical diagnosis for the 3 observers yielded κ values of 0.906, 0.912, and 0.906. The κ values for agreement among the 3 observers for the type of abnormality were 0.949, 0.845, and 0.895. There was complete agreement on the presence of an abnormality between the clinical and 4D-DSA in 23/26 cases. In 2 cases, there were conflicting opinions.

CONCLUSIONS

In this study, the information content of 4D-DSA reconstructions was largely equivalent to that of the combined 2D/3D studies. The availability of 4D-DSA should reduce the requirement for 2D-DSA acquisitions.

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.

Brain Vascular Imaging Techniques

Recent major improvements in a number of imaging techniques now allow for the study of the brain in ways that could not be considered previously. Researchers today have well-developed tools to specifically examine the dynamic nature of the blood vessels in the brain during development and adulthood; as well as to observe the vascular responses in disease situations in vivo. This review offers a concise summary and brief historical reference of different imaging techniques and how these tools can be applied to study the brain vasculature and the blood-brain barrier integrity in both healthy and disease states. Moreover, it offers an overview on available transgenic animal models to study vascular biology and a description of useful online brain atlases.

Time-resolved 3D rotational angiography: display of detailed neurovascular anatomy in patients with intracranial vascular malformations

PURPOSE

The purpose of this pilot study was to demonstrate the applicability of time-resolved three-dimensional (3D) reconstructions from 3D digital subtraction angiography (DSA) rotational angiography (RA) datasets (four-dimensional (4D) DSA) to provide a more detailed display of the architecture of intracranial vascular malformations.

METHODS

The experimental reconstruction software was applied to the existing 3D DSA datasets obtained with Siemens Artis zee biplane neuroangiography equipment. We included 27 patients with clinical indications for 3DRA for preinterventional or preoperative evaluation of intracranial dural arteriovenous fistulas (dAVFs, n=8) or arteriovenous malformations (AVMs, n=19). A modified DSA acquisition protocol covering an extended rotation angle of the C-arm of 260° during a scan time of 12 s was used. 4D volumes were displayed with up to 30 frames/s in a transparent volume rendering (VRT) mode and time-resolved multiplanar reconstructions (MPRs). Arterial feeders, fistulous points, or the shunt zone within the AVM nidus and venous drainage patterns as well as associated aneurysms were assessed after definition of a standardized evaluation procedure by consensus of two reviewers in comparison with 2D DSA and conventional 3D reconstructions.

RESULTS

In all cases calculation of 4D reconstructions were technically feasible and evaluable. In two cases image quality was slightly compromised by movement artifacts. Compared with standard DSA projection images and 3D reconstructions, 4D VRTs and MPRs were rated significantly superior to define a proper projection and display of the shunt zone. In 12 out of 27 cases 4D reconstructions showed details of the angioarchitecture at the fistulous point or the nidus better than the other modalities and came close to the quality of superselective angiography. The efficacy of 3D and 4D applications was equal in the detection of pre- and intranidal aneurysms. The course of long arterial feeders and draining veins was difficult to assess on VRTs and MPRs. Especially for dAVFs, 2D DSA was clearly superior in identifying meningeal feeders. For detecting smaller vessels and for distinction between angiographic phases, 2D DSA is still considered to be superior to 4D imaging. Venous drainage was slightly better displayed in 4D reconstructions.

CONCLUSIONS

Time-resolved 3DRA with 4D VRTs and MPRs is technically feasible and provides a detailed display of the angioarchitecture at the fistulous point or the nidus. Visualization of all angiographic features demands additional post-processing. Further standardization of evaluation tools and studies with blinded independent reviewers are necessary before the new technique can replace conventional neuroangiographic approaches.

Application of 4-Dimensional Digital Subtraction Angiography for Dural Arteriovenous Fistulas

BACKGROUND

Three-dimensional reconstruction of digital subtraction angiography (3D-DSA) is a useful imaging modality for assessing complex cerebrovascular lesions. However, due to the importance of flow over time in certain vascular lesions, 3D-DSA is of less value as it lacks the temporal resolution. Dural arteriovenous fistulas (AVFs) are complex lesions in which an arteriovenous shunt exists between meningeal arteries and a dural venous sinus or cortical vein. Traditional 2D-DSA, especially with superselective injections of feeding arteries, is currently the gold standard for assessment, but overlapping of opacified vessels can complicate interpretation. A novel imaging technique, 4D-DSA, merges 3D reconstructions of multiple temporal series. It offers a unique perspective on complex cerebrovascular lesions and may offer several advantages in the assessment of dural AVF.

METHODS

4D-DSA images were acquired in 5 patients who presented with dural AVFs. All relevant clinical data, imaging, and procedural/operative reports were reviewed retrospectively. 4D-DSA images were reconstructed on a separate 3D workstation and compared to 2D and 3D-DSA images in an offline fashion.

RESULTS

In all 5 cases, 4D-DSA proved to be useful in lesion assessment and treatment planning. This included observation (n = 2), microsurgery (n = 1), and endovascular embolization (n = 2).

CONCLUSIONS

In the small series of patients in which it has been evaluated, 4D-DSA offers several advantages in assessing dural AVFs. The ability to see and manipulate feeding arteries in 3D combined with temporal resolution was useful in assessment and treatment planning. Continued experience with this imaging technique will be needed to identify its optimal use.

Finding the optimal deconvolution algorithm for MR perfusion in carotid stenosis: Correlations with angiographic cerebral circulation time

PURPOSE

The aim of our study is to explore the impacts of different deconvolution algorithms on correlations between CBF, MTT, CBV, TTP, Tmax from MR perfusion (MRP) and angiography cerebral circulation time (CCT).

METHODS

Retrospectively, 30 patients with unilateral carotid stenosis, and available pre-stenting MRP and angiography were included for analysis. All MRPs were conducted in a 1.5-T MR scanner. Standard singular value decomposition, block-circulant, and two delay-corrected algorithms were used as the deconvolution methods. All angiographies were obtained in the same bi-plane flat-detector angiographic machine. A contrast bolus of 12mL was administrated via angiocatheter at a rate of 8mL/s. The acquisition protocols were the same for all cases. CCT was defined as the difference between time to peak from the cavernous ICA and the parietal vein in lateral view. Pearson correlations were calculated for CCT and CBF, MTT, CBV, TTP, Tmax.

RESULTS

The correlation between CCT and MTT was highest with Tmax (r=0.65), followed by MTT (r=0.60), CBF (r=-0.57), and TTP (r=0.33) when standard singular value decomposition was used. No correlation with CBV was noted.

CONCLUSIONS

MRP using a singular value decomposition algorithm confirmed the feasibility of quantifying cerebral blood flow deficit in steno-occlusive disease within the angio-room. This approach might further improve patient safety by providing immediate cerebral hemodynamics without extraradiation and iodine contrast.