Norbert J. Pelc

Generalized image combinations in dual KVP digital radiography

Dual energy basis decomposition techniques apply to single projection radiographic imaging. The high and low energy images are non-linearly transformed to generate two energy-independent images characterizing the integrated Compton/photoelectric attenuation components. Characteristic linear combinations of these two basis images identify unknown materials, cancel known materials, and generate synthesized monoenergetic images. The problems of intervening materials and material displacement are solved in general for a wide class of clinical imaging tasks. The basis projection angle identifies one from a family of energy selective imaging tasks, and such performance measures as the contrast enhancement factor (CEF) and signal to noise ratio (SNR) are expressed as functions of this angle. Algorithms for the decomposition of high and low energy measurements are compared and experimental images are included.

Iterative Decomposition of Water and Fat With Echo Asymmetry and Least-Squares Estimation (IDEAL): Application With Fast Spin-Echo Imaging

Chemical shift based methods are often used to achieve uniform water–fat separation that is insensitive to Bo inhomogeneities. Many spin‐echo (SE) or fast SE (FSE) approaches acquire three echoes shifted symmetrically about the SE, creating time‐dependent phase shifts caused by water–fat chemical shift. This work demonstrates that symmetrically acquired echoes cause artifacts that degrade image quality. According to theory, the noise performance of any water–fat separation method is dependent on the proportion of water and fat within a voxel, and the position of echoes relative to the SE. To address this problem, we propose a method termed “iterative decomposition of water and fat with echo asymmetric and least‐squares estimation” (IDEAL). This technique combines asymmetrically acquired echoes with an iterative least‐squares decomposition algorithm to maximize noise performance. Theoretical calculations predict that the optimal echo combination occurs when the relative phase of the echoes is separated by 2π/3, with the middle echo centered at π/2+πk (k = any integer), i.e., (–π/6+πk, π/2+πk, 7π/6+πk). Only with these echo combinations can noise performance reach the maximum possible and be independent of the proportion of water and fat. Close agreement between theoretical and experimental results obtained from an oil–water phantom was observed, demonstrating that the iterative least‐squares decomposition method is an efficient estimator. Magn Reson Med, 2005.

Unaliasing by Fourier-Encoding the Overlaps Using the Temporal Dimension (UNFOLD), Applied to Cardiac Imaging and fMRI

In several applications, MRI is used to monitor the time behavior of the signal in an organ of interest; e.g., signal evolution because of physiological motion, activation, or contrast‐agent accumulation. Dynamic applications involve acquiring data in a k–t space, which contains both temporal and spatial information. It is shown here that in some dynamic applications, the t axis of k–t space is not densely filled with information. A method is introduced that can transfer information from the k axes to the t axis, allowing a denser, smaller k–t space to be acquired, and leading to significant reductions in the acquisition time of the temporal frames.

Multicoil Dixon Chemical Species Separation with an Iterative Least-Squares Estimation Method

This work describes a new approach to multipoint Dixon fat–water separation that is amenable to pulse sequences that require short echo time (TE) increments, such as steady‐state free precession (SSFP) and fast spin‐echo (FSE) imaging. Using an iterative linear least‐squares method that decomposes water and fat images from source images acquired at short TE increments, images with a high signal‐to‐noise ratio (SNR) and uniform separation of water and fat are obtained. This algorithm extends to multicoil reconstruction with minimal additional complexity. Examples of single‐ and multicoil fat–water decompositions are shown from source images acquired at both 1.5T and 3.0T. Examples in the knee, ankle, pelvis, abdomen, and heart are shown, using FSE, SSFP, and spoiled gradient‐echo (SPGR) pulse sequences. The algorithm was applied to systems with multiple chemical species, and an example of water–fat–silicone separation is shown. An analysis of the noise performance of this method is described, and methods to improve noise performance through multicoil acquisition and field map smoothing are discussed. Magn Reson Med 51:35–45, 2004.

The anatomy of the posterior communicating artery as a risk factor for ischemic cerebral infarction.

BACKGROUND After the occlusion of an internal carotid artery the principal source of collateral flow is through the arteries of the circle of Willis, but the size and patency of these arteries are quite variable. Study of the anatomy of the collateral pathways in patients with internal-carotid-artery occlusion with or without infarction in the watershed area of the deep white matter may identify patterns that afford protection from ischemic infarction. METHODS Using conventional magnetic resonance imaging and three-dimensional phase-contrast magnetic resonance angiography, we evaluated 29 consecutive patients (32 hemispheres at risk) with angiographically proved occlusion of the internal carotid artery. Four collateral pathways to the occluded vessel were evaluated: the proximal segment of the anterior cerebral artery, the posterior communicating artery, the ophthalmic artery, and leptomeningeal collateral vessels from the posterior cerebral artery. RESULTS Only features of the ipsilateral posterior communicating artery were related to the risk of watershed infarction. The presence of posterior communicating arteries measuring at least 1 mm in diameter was associated with the absence of watershed infarction (13 hemispheres, no infarcts; P < 0.001). Conversely, there were 4 watershed infarcts in the 6 hemispheres with posterior communicating arteries measuring less than 1 mm in diameter and 10 infarcts in the 13 hemispheres with no detectable flow in the ipsilateral posterior communicating artery. CONCLUSIONS A small (< 1 mm in diameter) or absent ipsilateral posterior communicating artery is a risk factor for ischemic cerebral infarction in patients with internal-carotid-artery occlusion.

Time-resolved three-dimensional phase-contrast MRI.

To demonstrate the feasibility of a four‐dimensional phase contrast (PC) technique that permits spatial and temporal coverage of an entire three‐dimensional volume, to quantitatively validate its accuracy against an established time resolved two‐dimensional PC technique to explore advantages of the approach with regard to the four‐dimensional nature of the data.

Comparison of linear and circular polarization for magnetic resonance imaging

Abstract A comparison of experimental imaging results obtained with linearly polarized and circularly polarized radiofrequency excitation and reception is presented. Simulation images in good agreement with the experimental scans are described. The simulations are calculated with a model in which a homogeneous, isotropic cylinder of lossy dielectric material and infinite axial extent is immersed in a uniform rf magnetic field perpendicular to the axis. It is found that with the usual linear polarization, reconstructions of uniform objects have regions of decreased intensity. These artifacts are shown to arise from dielectric standing wave effects and eddy currents. The effects become more severe as the frequency or object size is increased, and depend upon the complex conductivity of the object. Results indicate that a significant reduction in the artifact intensity is achieved when circular polarization is employed for both transmission and reception. The expected benefits of circular polarization over linear polarization in reduction of excitation power (up to 50% reduction) and signal-to-noise advantage (√2) have been realized in practice with cylindrical objects and human subjects.

An algorithm for the reduction of metal clip artifacts in CT reconstructions.

Implanted surgical metal clips often produce objectionable artifacts in CT reconstructions. The artifacts appear as streaks which emanate radially from the site of the clip. It is shown in this paper that these artifacts stem primarily from motion of the clip during the scan. An algorithm is described which reduces the intensity of these artifacts. The procedure attempts to remove the metal object entirely from the scan data by replacing the measured projection values of rays that passed through a neighborhood of the clip with calculated values consistent with an object whose density is an average of the surround. Examples are given for head and body scans as well as for computer simulations which show substantial reduction of the streak intensity. Key words: computed tomography, artifacts, reconstruction algorithms, surgical clip streaks, starburst artifacts

Field map estimation with a region growing scheme for iterative 3-point water-fat decomposition

Robust fat suppression techniques are required for many clinical applications. Multi‐echo water‐fat separation methods are relatively insensitive to B0 field inhomogeneity compared to the fat saturation method. Estimation of this field inhomogeneity, or field map, is an essential and important step, which is well known to have ambiguity. For an iterative water‐fat decomposition method recently proposed, ambiguities still exist, but are more complex in nature. They were studied by analytical expressions and simulations. To avoid convergence to incorrect field map solutions, an initial guess closer to the true field map is necessary. This can be achieved using a region growing process, which correlates the estimation among neighboring pixels. Further improvement in stability is achieved using a low‐resolution reconstruction to guide the selection of the starting pixels for the region growing. The proposed method was implemented and shown to significantly improve the algorithms immunity to field inhomogeneity. Magn Reson Med, 2005.

Time-resolved 3-dimensional Velocity Mapping in the Thoracic Aorta: Visualization of 3-directional Blood Flow Patterns in Healthy Volunteers and Patients

Objective: An analysis of thoracic aortic blood flow in normal subjects and patients with aortic pathologic findings is presented. Various visualization tools were used to analyze blood flow patterns within a single 3-component velocity volumetric acquisition of the entire thoracic aorta Methods: Time-resolved, 3-dimensional phase-contrast magnetic resonance imaging (3D CINE PC MRI) was employed to obtain complete spatial and temporal coverage of the entire thoracic aorta combined with spatially registered 3-directional pulsatile blood flow velocities. Three-dimensional visualization tools, including time-resolved velocity vector fields reformatted to arbitrary 2-dimensional cut planes, 3D streamlines, and time-resolved 3D particle traces, were applied in a study with 10 normal volunteers. Results from 4 patient examinations with similar scan prescriptions to those of the volunteer scans are presented to illustrate flow features associated with common pathologic findings in the thoracic aorta. Results: Previously reported blood flow patterns in the thoracic aorta, including right-handed helical outflow, late systolic retrograde flow, and accelerated passage through the aortic valve plane, were visualized in all volunteers. The effects of thoracic aortic disease on spatial and temporal blood flow patterns are illustrated in clinical cases, including ascending aortic aneurysms, aortic regurgitation, and aortic dissection. Conclusion: Time-resolved 3D velocity mapping was successfully applied in a study of 10 healthy volunteers and 4 patients with documented aortic pathologic findings and has proven to be a reliable tool for analysis and visualization of normal characteristic as well as pathologic flow features within the entire thoracic aorta.