Four Dimensional Intravenous Cone-Beam Computed Tomographic Subtraction Angiography: In Vitro Study of Feasibility

Feasibility of texture analysis for the assessment of biochemical changes in meniscal tissue on T1 maps calculated from delayed gadolinium-enhanced magnetic resonance imaging of cartilage data: comparison with conventional relaxation time measurements

OBJECTIVES

To (1) establish the feasibility of texture analysis for the in vivo assessment of biochemical changes in meniscal tissue on delayed gadolinium-enhanced magnetic resonance imaging of cartilage (dGEMRIC), and (2) compare textural with conventional T1 relaxation time measurements calculated from dGEMRIC data ("T1(Gd) relaxation times").

MATERIALS AND METHODS

We enrolled 10 asymptomatic volunteers (7 men and 3 women; mean age, 27.2 +/- 4.5 years), without a history of meniscus damage, in our study. MRI of the right knee was performed at 3.0 T. An isotropic, 3-dimensional (3D), double-echo steady-state sequences was used for morphologic evaluation, and a dual flip angle 3D gradient echo sequence was used for T1(Gd) mapping. All MRI scans were performed 90 minutes after injection of 0.2 mmol/kg of Gd-diethylenetriamine pentaacetic acid (DTPA), and subsequently, during application of a compressive force (50% of the body weight) in the axial direction. Regions of interest, covering the central portions of the posterior horn of the medial meniscus, were defined on 3 adjacent sagittal sections. Based on the relaxation time maps, mean T1(Gd), as well as the T1(Gd) texture features derived from the co-occurrence matrix (COC: Angular Second Moment, Entropy, Inverse Difference Moment) and wavelet transform (WAV: WavEnLL, WavEnHL, WavEnHH, WavEnLH), were calculated. Paired t tests were used to assess differences between baseline and compression, and intraclass correlation coefficients (ICC) were calculated to establish the intrarater reliability of the measurements.

RESULTS

Mean T1(Gd) (-67.3 ms, P = 0.011), Angular Second Moment (-0.0002, P = 0.009), Entropy (+0.033, P = 0.025), WavEnLL (+1011.16, P = 0.002), WavEnHL (+18.64, P = 0.012), and WavEnLH (+72.74, P = 0.035) differed significantly between baseline and compression. Intrarater reliability was substantial for mean T1(Gd) relaxation times (ICC = 0.99-1.0), and also for T1(Gd) co-occurrence matrix (ICC = 0.63-0.92) and WAV (ICC = 0.86-0.98) features.

CONCLUSIONS

Texture features extracted from T1 maps calculated from dGEMRIC data are feasible for the in vivo assessment of biochemical changes in the menisci, such as might be induced by mechanical loading. Thus, T1(Gd) texture features complement conventional relaxation time measurements. Further studies are necessary to determine whether the mechanical compression, or a prolonged Gd-DTPA uptake, or both, are responsible for the observed decrease in mean T1(Gd) relaxation times in the menisci.

Sub-Nyquist acquisition and constrained reconstruction in time resolved angiography

In 1980 DSA provided a real time series of digitally processed angiographic images that facilitated and reduced the risk of angiographic procedures. This technique has become an enabling technology for interventional radiology. Initially it was hoped that intravenous DSA could eliminate the need for arterial injections. However the 2D nature of the images resulted in overlap of vessels and repeat injections were often required. Ultimately the use of smaller arterial catheters and reduced iodine injections resulted in significant reduction in complications. During the next two decades time resolved MR DSA angiographic methods were developed that produced time series of 3D images. These 4D displays were initially limited by tradeoffs in temporal and spatial resolution when acquisitions obeying the Nyquist criteria were employed. Then substantial progress was made in the implementation of undersampled non-Cartesian acquisitions such as VIPR and constrained reconstruction methods such as HYPR, which removed this tradeoff and restored SNR usually lost by accelerated techniques. Recently, undersampled acquisition and constrained reconstruction have been applied to generate time series of 3D x-ray DSA volumes reconstructed using rotational C-arm acquisition completing a 30 year evolution from DSA to 4D DSA. These 4D DSA volumes provide a flexible series of roadmaps for interventional procedures and solve the problem of vessel overlap for intravenous angiography. Full time-dependent behavior can be visualized in three dimensions. When a biplane system is used, 4D fluoroscopy is also possible, enabling the interventionalist to track devices in vascular structures from any angle without moving the C-arm gantrys. Constrained reconstruction methods have proved useful in a broad range of medical imaging applications, where substantial acquisition accelerations and dose reductions have been reported.

Intravascular functional maps of common neurovascular lesions derived from volumetric 4D CT data

PURPOSE

Current computed tomography angiography (CTA) postprocessing tools do not support quantitative assessment of intravascular physiology. Dynamic volumetric CT, acquired at a sufficiently high temporal resolution, is ideal for such analysis. Following preliminary experiments in flow phantoms, we examine the segmentation of blood vessels from 4D CT angiography by curve fit and encoding of functional blood flow information into the resulting functional intravascular maps.

MATERIALS AND METHODS

Flow phantoms were constructed consisting of a single pipe input and 4 simultaneous outputs of varying flow rates. Two outflow pipe diameters were tested. Bolus transit time (TT), time to peak (TTP), and time of arrival (TOA) were analyzed using contrast bolus profiles generated from 4D volumetric CT examinations on a 320 detector scanner in regions of interest placed 10 cm apart in all outflow pipes. Six subjects with various neurovascular lesions were next examined using a volumetric contrast-enhanced 4D CT angiography protocol. Segmentation was performed by quadratic curve fit after comparative analysis and optimization of the segmentation technique using quadratic curves, the gamma variate function, and a simplified formulation of the gamma variate function. After segmentation, quantitative analysis of spatially congruent intravascular voxels including TTP, rise, TT, and slope of the contrast upstroke was employed to encode physiologic information into the segmentations and produce intravascular functional maps. Comparison was made in each case to the patient's routine imaging.

RESULTS

Increasing volumetric flow rates correspond to reduction of bolus TT in flow phantoms. TT elongation was observed as the contrast bolus moved distally in all pipes, with greater elongation seen at slower flow rates and larger pipe diameters. A greater difference was observed between TTP proximally and distally in pipes compared with TOA, an effect most prominent at slower flow rates and larger pipe lumens, and thus TTP was chosen for functional encoding into segmentations of the clinical series. In vivo, the quadratic function demonstrated the lowest coefficient of variation when fit to intravascular time density series and outperformed 2 formulations of the gamma variate function. After segmentation with quadratic curves, Gaussian distributions were chosen over gamma variate functions to characterize contrast bolus profiles while neglecting recirculation and to calculate functional parameters for spatial encoding. Intravascular functional maps free of bone artifacts were created in every case that demonstrated all appropriate vessels and showed agreement with conventional imaging modalities in terms of vessel delineation and the diagnosis of vascular pathology. The most useful and interesting functional maps are discussed in each case.

CONCLUSIONS

The above approach to quantitative CT angiography provides a method of evaluating dynamic CTA data by means of intravascular functional maps. The techniques are broadly applicable in the clinical assessment of a variety of vascular diseases.