Extending the dynamic range of biomedical micro-computed tomography for application to geomaterials.

BACKGROUND: X-ray computed tomography (CT) can non-destructively examine objects by producing three-dimensional images of their internal structure. Although the availability of biomedical micro-CT offers the increased access to scanners, CT images of dense objects are susceptible to artifacts particularly due to beam hardening. OBJECTIVE: This study proposes and evaluates a simple semi-empirical correction method for beam hardening and scatter that can be applied to biomedical scanners. METHODS: Novel calibration phantoms of varying diameters were designed and built from aluminum and poly[methyl-methacrylate]. They were imaged using two biomedical micro-CT scanners. Absorbance measurements made through different phantom sections were fit to polynomial and inversely exponential functions and used to determine linearization parameters. Corrections based on the linearization equations were applied to the projection data before reconstruction. RESULTS: Correction for beam hardening was achieved when applying both scanners with the correction methods to all test objects. Among them, applying polynomial correction method based on the aluminum phantom provided the best improvement. Correction of sample data demonstrated a high agreement of percent-volume composition of dense metallic inclusions between using the Bassikounou meteorite from the micro-CT images (13.7%) and previously published results using the petrographic thin sections (14.6% 8% metal and 6.6% troilite). CONCLUSIONS: Semi-empirical linearization of X-ray projection data with custom calibration phantoms allows accurate measurements to be obtained on the radiodense samples after applying the proposed correction method on biomedical micro-CT images.

Quantifying lung morphology with respiratory-gated micro-CT in a murine model of emphysema.

Non-invasive micro-CT imaging techniques have been developed to investigate lung structure in free-breathing rodents. In this study, we investigate the utility of retrospectively respiratory-gated micro-CT imaging in an emphysema model to determine if anatomical changes could be observed in the image-derived quantitative analysis at two respiratory phases. The emphysema model chosen was a well-characterized, genetically altered model (TIMP-3 knockout mice) that exhibits a homogeneous phenotype. Micro-CT scans of the free-breathing, anaesthetized mice were obtained in 50 s and retrospectively respiratory sorted and reconstructed, providing 3D images representing peak inspiration and end expiration with 0.15 mm isotropic voxel spacing. Anatomical measurements included the volume and CT density of the lungs and the volume of the major airways, along with the diameters of the trachea, left bronchus and right bronchus. From these measurements, functional parameters such as functional residual capacity and tidal volume were calculated. Significant differences between the wild-type and TIMP-3 knockout groups were observed for measurements of CT density over the entire lung, indicating increased air content in the lungs of TIMP-3 knockout mice. These results demonstrate retrospective respiratory-gated micro-CT, providing images at multiple respiratory phases that can be analyzed quantitatively to investigate anatomical changes in murine models of emphysema.

Rigid-body motion correction in hybrid PET/MRI using spherical navigator echoes.

Integrated positron emission tomography and magnetic resonance imaging (PET/MRI) is an imaging technology that provides complementary anatomical and functional information for medical diagnostics. Both PET and MRI are highly susceptible to motion artifacts due, in part, to long acquisition times. The simultaneous acquisition of the two modalities presents the opportunity to use MRI navigator techniques for motion correction of both PET and MRI data. For this task, we propose spherical navigator echoes (SNAVs)-3D k-space navigators that can accurately and rapidly measure rigid body motion in all six degrees of freedom. SNAVs were incorporated into turbo FLASH (tfl)-a product fast gradient echo sequence-to create the tfl-SNAV pulse sequence. Acquiring in vivo brain images from a healthy volunteer with both sequences first compared the tfl-SNAV and product tfl sequences. It was observed that incorporation of the SNAVs into the image sequence did not have any detrimental impact on the image quality. The SNAV motion correction technique was evaluated using an anthropomorphic brain phantom. Following a stationary reference image where the tfl-SNAV sequence was acquired along with simultaneous list-mode PET, three identical PET/MRI scans were performed where the phantom was moved several times throughout each acquisition. This motion-up to 11° and 14 mm-resulted in motion artifacts in both PET and MR images. Following SNAV motion correction of the MRI and PET list-mode data, artifact reduction was achieved for both the PET and MR images in all three motion trials. The corrected images have improved image quality and are quantitatively more similar to the ground truth reference images.

Implementation of dual- and triple-energy cone-beam micro-CT for postreconstruction material decomposition.

Micro-CT has become a powerful tool for small animal research, having the ability to obtain high-resolution in vivo and ex vivo images for analyzing bone mineral content, organ vasculature, and bone microarchitecture extraction. The use of exogenous contrast agents further extends the use of micro-CT techniques, but despite advancements in contrast agents, single-energy micro-CT is still limited in cases where two different materials share similar grey-scale intensity values. This study specifically addresses the development of multiple-energy cone-beam micro-CT, for applications where bone must be separated from blood vessels filled with a Pb-based contrast material (Microfil) in ex vivo studies of rodents and tissue specimens. The authors report the implementation of dual- and triple-energy CT algorithms for material-specific imaging using postreconstruction decomposition of micro-CT data; the algorithms were implemented on a volumetric cone-beam micro-CT scanner (GE Locus Ultra). For the dual-energy approach, extrinsic filtration was applied to the x-ray beam to produce spectra with different proportions of x rays above the K edge of Pb. The optimum x-ray tube energies (140 kVp filtered with 1.45 mm Cu and 96 kVp filtered with 0.3 mm Pb) that maximize the contrast between bone and Microfil were determined through numerical simulation. For the triple-energy decomposition, an additional low-energy spectrum (70 kVp, no added filtration) was used. The accuracy of decomposition was evaluated through simulations and experimental verification of a phantom containing a cortical bone simulating material (SB3), Microfil, and acrylic. Using simulations and phantom experiments, an accuracy greater than 95% was achieved in decompositions of bone and Microfil (for noise levels lower than 11 HU), while soft tissue was separated with accuracy better than 99%. The triple-energy technique demonstrated a slightly higher, but not significantly different, decomposition accuracy than the dual-energy technique for the same achieved noise level in the micro-CT images acquired at the multiple energies. The dual-energy technique was applied to the decomposition of an ex vivo rat specimen perfused with Microfil; successful decomposition of the bone and Microfil was achieved, enabling the visualization and characterization of the vasculature both in areas where the vessels traverse soft tissue and when they are surrounded by bone. In comparison, in single energy micro-CT, vessels surrounded by bone could not be distinguished from the cortical bone, based on grey-scale intensity alone. This work represents the first postreconstruction application of material-specific decomposition that directly takes advantage of the K edge characteristics of a contrast material injected into an animal specimen; the application of the technique resulted in automatic, accurate segmentation of 3D micro-CT images into bone, vessel, and tissue components. The algorithm uses only reconstructed images, rather than projection data, and is calibrated by an operator with signal values in regions identified as being comprised entirely of either cortical bone, contrast-enhanced vessel, or soft tissue; these required calibration values are observed directly within reconstructed CT images acquired at the multiple energies. These features facilitate future implementation on existing research micro-CT systems.

In vivo small-animal imaging using micro-CT and digital subtraction angiography.

Small-animal imaging has a critical role in phenotyping, drug discovery and in providing a basic understanding of mechanisms of disease. Translating imaging methods from humans to small animals is not an easy task. The purpose of this work is to review in vivo x-ray based small-animal imaging, with a focus on in vivo micro-computed tomography (micro-CT) and digital subtraction angiography (DSA). We present the principles, technologies, image quality parameters and types of applications. We show that both methods can be used not only to provide morphological, but also functional information, such as cardiac function estimation or perfusion. Compared to other modalities, x-ray based imaging is usually regarded as being able to provide higher throughput at lower cost and adequate resolution. The limitations are usually associated with the relatively poor contrast mechanisms and potential radiation damage due to ionizing radiation, although the use of contrast agents and careful design of studies can address these limitations. We hope that the information will effectively address how x-ray based imaging can be exploited for successful in vivo preclinical imaging.

Morphology-Specific Discrimination between MS White Matter Lesions and Benign White Matter Hyperintensities Using Ultra-High-Field MRI.

BACKGROUND AND PURPOSE: Recently published North American Imaging in Multiple Sclerosis guidelines call for derivation of a specific radiologic definition of MS WM lesions and mimics. The purpose of this study was to use SWI and magnetization-prepared FLAIR images for sensitive differentiation of MS from benign WM lesions using the morphologic characteristics of WM lesions. MATERIALS AND METHODS: Seventeen patients with relapsing-remitting MS and 18 healthy control subjects were enrolled retrospectively. For each subject, FLAIR and multiecho gradient-echo images were acquired using 7T MR imaging. Optimized postprocessing was used to generate single-slice SWI of cerebral veins. SWI/FLAIR images were registered, and 3 trained readers performed lesion assessment. Morphology, location of lesions, and the time required for assessment were recorded. Analyses were performed on 3 different pools: 1) lesions of >3 mm, 2) nonconfluent lesions of >3 mm, and 3) nonconfluent lesions of >3 mm with no or a single central vein. RESULTS: The SWI/FLAIR acquisition and processing protocol enabled effective assessment of central veins and hypointense rims in WM lesions. Assessment of nonconfluent lesions with ≥1 central vein enabled the most specific and sensitive differentiation of patients with MS from controls. A threshold of 67% perivenous WM lesions separated patients with MS from controls with a sensitivity of 94% and specificity of 100%. Lesion assessment took an average of 12 minutes 10 seconds and 4 minutes 33 seconds for patients with MS and control subjects, respectively. CONCLUSIONS: Nonconfluent lesions of >3 mm with ≥1 central vein were the most sensitive and specific differentiators between patients with MS and control subjects.

In vivo characterization of lung morphology and function in anesthetized free-breathing mice using micro-computed tomography.

Lung morphology and function in human subjects can be monitored with computed tomography (CT). Because many human respiratory diseases are routinely modeled in rodents, a means of monitoring the changes in the structure and function of the rodent lung is desired. High-resolution images of the rodent lung can be attained with specialized micro-CT equipment, which provides a means of monitoring rodent models of lung disease noninvasively with a clinically relevant method. Previous studies have shown respiratory-gated images of intubated and respirated mice. Although the image quality and resolution are sufficient in these studies to make quantitative measurements, these measurements of lung structure will depend on the settings of the ventilator and not on the respiratory mechanics of the individual animals. In addition, intubation and ventilation can have unnatural effects on the respiratory dynamics of the animal, because the airway pressure, tidal volume, and respiratory rate are selected by the operator. In these experiments, important information about the symptoms of the respiratory disease being studied may be missed because the respiration is forced to conform to the ventilator settings. In this study, we implement a method of respiratory-gated micro-CT for use with anesthetized free-breathing rodents. From the micro-CT images, quantitative analysis of the structure of the lungs of healthy unconscious mice was performed to obtain airway diameters, lung and airway volumes, and CT densities at end expiration and during inspiration. Because the animals were free breathing, we were able to calculate tidal volume (0.09 +/- 0.03 ml) and functional residual capacity (0.16 +/- 0.03 ml).

Magnetic resonance imaging compatible remote catheter navigation system with 3 degrees of freedom.

PURPOSE: To facilitate MRI-guided catheterization procedures, we present an MRI-compatible remote catheter navigation system that allows remote navigation of steerable catheters with 3 degrees of freedom. METHODS: The system consists of a user interface (master), a robot (slave), and an ultrasonic motor control servomechanism. The interventionalist applies conventional motions (axial, radial and plunger manipulations) on an input catheter in the master unit; this user input is measured and used by the servomechanism to control a compact catheter manipulating robot, such that it replicates the interventionalist's input motion on the patient catheter. The performance of the system was evaluated in terms of MRI compatibility (SNR and artifact), feasibility of remote navigation under real-time MRI guidance, and motion replication accuracy. RESULTS: Real-time MRI experiments demonstrated that catheter was successfully navigated remotely to desired target references in all 3 degrees of freedom. The system had an absolute value error of [Formula: see text]1 mm in axial catheter motion replication over 30 mm of travel and [Formula: see text] for radial catheter motion replication over [Formula: see text]. The worst case SNR drop was observed to be [Formula: see text]3 %; the robot did not introduce any artifacts in the MR images. CONCLUSION: An MRI-compatible compact remote catheter navigation system has been developed that allows remote navigation of steerable catheters with 3 degrees of freedom. The proposed system allows for safe and accurate remote catheter navigation, within conventional closed-bore scanners, without degrading MR image quality.

A laboratory CT scanner for dynamic imaging.

A high-resolution laboratory CT scanner has been developed for imaging objects undergoing periodic motion. The scanner comprises an x-ray image intensifier, optically coupled to a linear photodiode array. Gated time-evolved projections of a single slice of the moving object are acquired, reformatted, and reconstructed. The resulting series of CT images shows the object at different phases of its motion cycle. The scanner has an adjustable field of view (FOV) and the resolution can be as high as 3.2 mm-1 (for the 40-mm FOV). The spatial resolution depends on the inherent resolution of the scanner and on the object's velocity. For objects moving at 1 cm s-1, the spatial resolution is reduced by 9% in the direction of motion. The signal intensity in the reconstructed image is linear for materials with attenuation coefficients as high as 1.5 cm-1 (for a 90-kVp x-ray beam), with an average accuracy of +/- 0.02 cm-1. The average accuracy of circumference measurements made from the CT images is +/- 0.3 mm. Lastly, an application of this dynamic CT scanner to imaging excised human arterial specimens under simulated physiological pressure conditions is presented as an example.

Measurement of the spatial Wiener spectrum of nonstorage imaging devices.

All previous methods for measuring image noise spectra require a noise realization, a static image, typified as a photograph which can be scanned to create the Wiener spectrum. We wished to analyze the spatial noise power spectrum at the output phosphor of a continuously irradiated imaging device, an x-ray image intensifier (XRII), which is incapable of image storage and thus the image is continually changing as a function of both time and space. Our new method utilizes a pair of slits to measure the relative Wiener spectrum of the temporally changing components of the image (i.e., x-ray quantum and XRII gain noises). By measuring the modulation transfer function and the Wiener spectrum of the same XRII on the same apparatus it was possible to demonstrate the spatial frequency dependence of the detective quantum efficiency. Adaptations of the method should permit the measurement of Wiener spectra of fluoroscopic television systems directly from the TV monitor.

Optical factors affecting the detective quantum efficiency of radiographic screens.

Parameters related to the detective quantum efficiency (DQE) of several representative screens of different thicknesses, phosphor grain sizes, and optical properties were measured by the scintillation spectrum method, using monoenergetic x rays produced from x-ray fluorescence. The experimental results, including those for spectral shape and average light energies (EA) emitted, are compared with conventional theories of the operation of screens. It was hoped that this would vindicate the theory of the effect of optical properties and so permit the simple calculation of all parameters related to DQE from standard x-ray attenuation tables. Rather more substantial energy-dependent deviations of EA are found than was previously realized, which preliminary analysis suggests are due to both optical effects and photoelectron escape. We conclude that although DQE for a single energy can be calculated by simplified methods to within +/- 10%, the effective DQE when polyenergetic beams are used is much less accurately estimated and requires a fuller theoretical treatment.

A high-resolution XRII-based quantitative volume CT scanner.

A laboratory volume CT scanner has been developed, with high spatial resolution in all three dimensions, which can be used for quantitative analysis of excised tissue samples in vitro. The system incorporates an x-ray image intensifier, optically coupled to a time-delay integration (TDI) CCD to obtain low-noise and low-scatter projections of the sample volume. A water bath surrounds the sample to equalize the exposure to the image intensifier, thereby reducing the dynamic range of the input signal. The scanner operates in two modes, producing either a single, transverse image through the sample or a three-dimensional image of the sample volume. Spatial resolution is adjustable over the range of 1.2 to 2.8 mm-1. System response is linear over the range -1000 to 3500 Houndsfield units (HU), with an average precision of +/- 80 HU. The precision of geometric measurements in the transverse plane allows circumference measurements to within +/- 0.1 mm. Finally, applications of this technique of nondestructive analysis in biomedical research are discussed.

Quantitative angiographic blood-flow measurement using pulsed intra-arterial injection.

A technique for quantitative blood-flow measurement using a novel pulsed injection of radiographic contrast agent is reported. A pressurized source of contrast agent is interrupted by a rotary valve at rates ranging from 1 to 30 Hz, producing well-defined boli at the end of a catheter. The position of these boli can be recorded by a digital radiographic system and analyzed by one of several previously reported techniques, to produce quantitative measurements of blood velocity and flow rate throughout the cardiac cycle. The contrast-agent flow wave form produced by the pulsed injector has been measured with an electromagnetic flow meter, for driving pressures ranging from 600 to 1500 kPa. Excellent modulation of the contrast agent is observed for injection frequencies up to 20 Hz, through catheters up to 100 cm in length. Preliminary in vitro angiographic flow measurements have been performed using an x-ray image intensifier, coupled to a linear photodiode array as the digital detector. Both constant flow and pulsatile human blood-flow wave forms were simulated within a 6.4-mm-diam straight tube and monitored with an electromagnetic flow meter. These experiments indicate that the pulsed injector can be used to provide estimates of arterial blood flow over the entire cardiac cycle (including reverse flow), to within about +/-11%, following injection of less than 10 ml of iodinated contrast agent.

Experimental and theoretical x-ray imaging performance comparison of iodine and lanthanide contrast agents.

Contrast agents based on the lanthanide elements gadolinium and holmium have recently been developed for magnetic resonance imaging (MRI). Because of the increased atomic number of these elements relative to iodine, these new compounds, used as x-ray contrast agents, may yield higher radiographic contrast, and hence improved x-ray image quality, relative to conventional iodinated compounds, for clinically useful x-ray spectra. This possibility has been investigated, in independent experimental and theoretical studies, for two x-ray imaging systems: a digital radiographic system, using an x-ray image intensifier (XRII) and charge-coupled device (CCD) detector; and a conventional screen/film system, using a Lanex Regular screen. Iodine, gadolinium, and holmium contrast agents were investigated over a wide range of concentration-thickness products (0.1-0.6 M cm) and diagnostic x-ray spectra (60-120 kVp). A simple theoretical model of x-ray detector response predicts the experimental radiographic contrast measurements with a mean absolute error of 8.0% for the XRII/CCD system and 5.9% for the screen/film system, and shows that the radiographic contrast for these two systems is representative of all XRII and screen/film systems. An index of image quality is defined, and its dependence on radiographic contrast, x-ray fluence per unit dose, and detective quantum efficiency (DQE) is shown. Theoretical values of the index, predicted by our model, are then used to compare the performance of the three contrast agents for the two systems investigated. In general, iodine performance decreases steadily with increasing kVp, gadolinium performance has a broad maximum near 85 kVp, and gadolinium outperforms holmium. Gadolinium outperforms iodine for spectra above (and vice versa below) about 72 kVp, depending slightly on spectrum filtration, object thickness, and detector type. Thus, raising the kVp to shorten exposure times or reduce x-ray tube heat loading results in a loss of image quality with iodine, but not with gadolinium. Similarly, beam-hardening artifacts in performing video densitometry with iodine would be reduced with gadolinium. Gadolinium-based contrast agents are thus shown to offer several practical advantages over conventional iodinated contrast agents.

A modified x-ray image intensifier with continuously variable field of view: resolution considerations.

A conventional x-ray image intensifier (XRII) has been modified to enable the field of view (FOV) to be varied continuously, by adjusting the potentials at the focusing electrodes. The benefit, to system resolution, from decreasing the FOV has been characterized by measuring the modulation transfer function (MTF) of the XRII coupled to a high-resolution photo-diode array (PDA), at a number of different FOVs achieved either by electronic or optical zooming. Electronic zooming of the XRII from FOV = 24 cm to FOV = 10 cm led to an increase in f0.1 (the frequency at which MTF = 0.1) from 1.41 to 3.05 mm-1, while optical zooming increased f0.1 from 1.41 mm-1 only to 1.88 mm-1. It is proposed that the advantage, with respect to resolution gain, of electronic zooming over optical zooming was realized only when the XRII limits system resolution. The MTF of the XRII coupled to a video camera, with lower resolving power than the PDA, was measured at different FOVs to show that using electronic zooming is only marginally beneficial when the optical detector and the XRII contribute equally to the resolution degradation. However, when a higher-resolution optical detector is used, electronic zooming always yields a greater gain in resolution.

A geometrically accurate vascular phantom for comparative studies of x-ray, ultrasound, and magnetic resonance vascular imaging: construction and geometrical verification.

A technique for producing accurate models of vascular segments for use in experiments that assess vessel geometry and flow has been developed and evaluated. The models are compatible with x-ray, ultrasound, and magnetic resonance (MR) imaging systems. In this paper, a model of the human carotid artery bifurcation, is evaluated that has been built using this technique. The phantom consists of a thin-walled polyester-resin replica of the bifurcation through which a blood-mimicking fluid may be circulated. The phantom is surrounded by an agar tissue-mimicking material and a series of fiducial markers. The blood- and tissue-mimicking materials have x-ray, ultrasound, and MR properties similar to blood and tissue; fiducial markers provide a means of aligning images acquired by different modalities. The root-mean-square difference between the inner wall geometry of the constructed model and the desired dimensions was 0.33 mm. Static images were successfully acquired using x-ray, ultrasound, and MR imaging systems, and are free of significant artifacts. Flow images acquired with ultrasound and MR agree qualitatively with each other, and with previously published flow patterns. Volume-flow measurements obtained with ultrasound and MR were within 4.4% of the actual values.

Optical detection of triggered atherosclerotic plaque disruption by fluorescence emission analysis.

Fluorescence emission analysis (FEA) has proven to be very sensitive for the detection of elastin, collagen and lipids, which are recognized as the major sources of autofluorescence in vascular tissues. FEA has also been reported to detect venous thromboemboli. In this paper we have tested the hypothesis that FEA can reproducibly detect in vivo and in vitro triggered plaque disruption and thrombosis in a rabbit model. Fluorescence emission (FE) spectra, recorded in vivo, detected Russell's viper venom (RVV)-induced transformation of atherosclerotic plaque. FE intensity at 410-490 nm 4 weeks after angioplasty was significantly lower (P < 0.0033 by analysis of variance) in RVV-treated rabbits when compared to control animals with stable plaque. FE spectral profile analyses also demonstrated a significant change in curve shape as demonstrated by polynomial regression analysis (R2 from 0.980 to 0.997). We have also demonstrated an excellent correlation between changes in FE intensity and the structural characteristics detected at different stages of "unstable atherosclerotic plaque" development using multiple regression analysis (R2 = 0.989). Thus, FEA applied in vivo is a sensitive and highly informative diagnostic technique for detection of triggered atherosclerotic plaque disruption and related structural changes, associated with plaque transformation, in a rabbit model.

Estimation of deformation gradient and strain from cine-PC velocity data.

Phase contrast magnetic resonance imaging (MRI) can provide in vivo myocardial velocity field measurements. These data allow densely spaced material points to be tracked throughout the whole heart cycle using, for example, the Fourier tracking algorithm. To process the tracking results for myocardial deformation and strain quantification, we developed a method that is based on fitting the tracking results to an appropriate local deformation model. We further analyzed the accuracy and precision of the method and provided performance predictions for several local models. In order to validate the method and the theoretical performance analysis, we conducted controlled computer simulations and a phantom study. The results agreed well with expectations. Human heart data were also acquired and analyzed, and provided encouraging results. At the signal-to-noise ratio (SNR) level and spatial resolution expected in clinical settings, the study predicts strain quantification accuracy and precision that may allow the technique to become a practical and powerful noninvasive approach for the study of cardiac function, although clinically acceptable data acquisition strategies for three-dimensional (3-D) data are still a challenge.

Three-dimensional echocardiography: assessment of inter- and intra-operator variability and accuracy in the measurement of left ventricular cavity volume and myocardial mass.

Accurate left ventricular (LV) volume and mass estimation is a strong predictor of cardiovascular morbidity and mortality. We propose that our technique of 3D echocardiography provides an accurate quantification of LV volume and mass by the reconstruction of 2D images into 3D volumes, thus avoiding the need for geometric assumptions. We compared the accuracy and variability in LV volume and mass measurement using 3D echocardiography with 2D echocardiography, using in vitro studies. Six operators measured the LV volume and mass of seven porcine hearts, using both 3D and 2D techniques. Regression analysis was used to test the accuracy of results and an ANOVA test was used to compute variability in measurement. LV volume measurement accuracy was 9.8% (3D) and 18.4% (2D); LV mass measurement accuracy was 5% (3D) and 9.2% (2D). Variability in LV volume quantification with 3D echocardiography was %SEMinter = 13.5%, %SEMintra = 11.4%, and for 2D echocardiography was %SEMinter = 21.5%, %SEMintra = 19.1%. We derived an equation to predict uncertainty in measurement of LV volume and mass using 3D echocardiography, the results of which agreed with our experimental results to within 13%. 3D echocardiography provided twice the accuracy for LV volume and mass measurement and half the variability for LV volume measurement as compared with 2D echocardiography.

In vitro verification of myocardial motion tracking from phase-contrast velocity data.

The ability to track motion from cine phase-contrast (PC) magnetic resonance (MR) velocity measurements was investigated using an in vitro model. A computer-controlled deformable phantom was used for the characterization of the accuracy and precision of the forward-backward and the compensated Fourier integration techniques. Trajectory accuracy is limited by temporal resolution when the forward-backward technique is used. With this technique the extent of the calculated trajectories is underestimated by an amount related to the motion period and the sequence repetition time, because of the band-limiting caused in the cine interpolation step. When the compensated Fourier integration technique is used, trajectory accuracy is independent of temporal resolution and is better than 1 mm for excursions of less than 15 mm, which are comparable to those observed in the myocardium. Measurement precision is dominated by the artifact level in the phase-contrast images. If no artifacts are present precision is limited by the inherent signal-to-noise ratio of the images. In the presence of artifacts, similar in magnitude to those observed in vivo, the reproducibility of tracking a 2.2 x 2.2 mm2 region of interest is better than 0.5 mm. When the Fourier integration technique is used, the improved accuracy is accompanied by a reduction in precision. We verified that tracking three-dimensional (3D) motion from velocity measurements of a single slice can lead to underestimations of the trajectory if there is a through-plane component of the motion that is not truly represented by the measured velocities. This underestimation can be overcome if volumetric cine phase-contrast velocity data are acquired and full three-dimensional analysis is performed.