David A. Chesler

Noise due to photon counting statistics in computed x-ray tomography

A general expression is derived for the noise due to photon counting statistics in computed X-ray tomography. The variance is inversely proportional to the cube of the resolution distance. For scanners using a water box, the noise in the reconstructed image depends inversely on the number of detected primary photons, summed over all angles, that have passed through a resolution element. Predictions of this formula agree well with the results of computer simulations. It is shown how this formula can be used to determine such parameters as required X-ray flux, detector counting rate, and dose, with special emphasis on tradeoffs between these parameters and resolution. It is also shown that to determine the X-ray attenuation coefficient of a resolution element to a given precision, the number of photons required by computed X-ray tomography is close to a theoretical limit.

Modeling Cerebral Blood Flow and Flow Heterogeneity from Magnetic Resonance Residue Data

Existing model-free approaches to determine cerebral blood flow by external residue detection show a marked dependence of flow estimates on tracer arrival delays and dispersion. In theory, this dependence can be circumvented by applying a specific model of vascular transport and tissue flow heterogeneity. The authors present a method to determine flow heterogeneity by magnetic resonance residue detection of a plasma marker. Probability density functions of relative flows measured in six healthy volunteers were similar among tissue types and volunteers, and were in qualitative agreement with literature measurements of capillary red blood cell and plasma velocities. Combining the measured flow distribution with a model of vascular transport yielded excellent model fits to experimental residue data. Fitted gray-to-white flow-rate ratios were in good agreement with PET literature values, as well as a model-free singular value decomposition (SVD) method in the same subjects. The vascular model was found somewhat sensitive to data noise, but showed far less dependence on vascular delay and dispersion than the model-free SVD approach.

Estimation of the Local Statistical Noise in Emission Computed Tomography

A simple modification of the filtered backprojection algorithm is presented for the computation of the local statistical noise in emission computed tomography. The technique is general in that any distribution of radioactivity may be accommodated. When applied to positron emission tomography, it is shown that the effects of photon absorption, random coincidences, radioactive decay, and detector nonuniformity may be included. Calculations have shown the effects of resolution, object size, and photon absorption on the statistical noise of disk-shaped emitters. Comparison of calculation and experiment show close agreement both in magnitude and spatial variation. Measurements of the noise level in tomograms of the brain obtained during continuous inhalation of 150-CO2 demonstrate that estimates of radioactivity concentration with a precision of a few percent are readily attainable.

Density of organic matrix of native mineralized bone measured by water- and fat-suppressed proton projection MRI.

Water‐ and fat‐suppressed projection MR imaging (WASPI) utilizes the large difference between the proton T  2* s of the solid organic matrix and the fluid constituents of bone to suppress the fluid signals while preserving solid matrix signals. The solid constituents include collagen and some molecularly immobile water and exhibit very short T  2* . The fluid constituents include mobile water and fat, with long T  2* . In WASPI, chemical shift selective low‐power π/2 pulses excite mobile water and fat magnetization which is subsequently dephased by gradient pulses, while the magnetization of collagen and immobile water remains mostly in the z‐direction. Additional selective π pulses in alternate scans further cancel the residual water and fat magnetization. Following water and fat suppression, the matrix signal is excited by a short hard pulse and the free induction decay acquired in the presence of a gradient in a 3D projection method. WASPI was implemented on a 4.7 T MR imaging system and tested on phantoms and bone specimens, enabling excellent visualization of bone matrix. The bone matrix signal per unit volume of bovine trabecular specimens was measured by this MR technique and compared with that determined by chemical analysis. This method could be used in combination with bone mineral density measurement by solid state 31P projection MRI to determine the degree of bone mineralization. Magn Reson Med 50:59–68, 2003.

Magnetic Resonance Imaging Mapping of Brain Function Human Visual Cortex

&NA; Belliveau JW, Kwong KK, Kennedy DN, Baker JR, Stern CE, Benson R, Chesler DA, Weisskoff RM, Cohen MS, Tootell RBH, Fox PT, Brady TJ, Rosen BR. Magnetic resonance imaging mapping of brain function: human visual cortex. Invest Radiol 1992;27:S59‐S65. Magnetic resonance imaging (MRI) studies of human brain activity are described. Task‐induced changes in brain cognitive state were measured using high‐speed MRI techniques sensitive to changes in cerebral blood volume (CBV), blood flow (CBF), and blood oxygenation. These techniques were used to generate the first functional MRI maps of human task activation, by using a visual stimulus paradigm. The methodology of MRI brain mapping and results from the investigation of the functional organization and frequency response of human primary visual cortex (V1) are presented.

Mr Velocity Imaging by Phase Display

The ability of the nuclear magnetic resonance signal to encode information about macroscopic motion has been recognized since the works of Hahn and Carr and Purcell. In the medical imaging setting this ability has led to a variety of ingenious magnetic resonance flow imaging schemes that ultimately may become competitive with X-ray angiography in sensitivity and specificity while remaining radically noninvasive. This work demonstrates that conventional spin-echo Fourier transform image acquisitions naturally encode a component of flow velocity that lies within the image plane. By displaying just the real part of the complex image data (phase display), the velocity distribution within the subject is revealed in image form. This method of flow imaging requires neither special pulse sequences nor image reconstruction and format software for its implementation. Further, images that intersect a flow channel longitudinally, demonstrating in-plane flow, yield an unusually large quantity of physiologic information per image. Phantom and in vivo flow images are presented. Also described is a phantom based on a rotating disk that enables calibration of the velocity/phase-shift constant for an untested pulse sequence.

Combined Diffusion-Weighted and Perfusion-Weighted Flow Heterogeneity Magnetic Resonance Imaging in Acute Stroke

BACKGROUND AND PURPOSE The heterogeneity of microvascular flows is known to be an important determinant of the efficacy of oxygen delivery to tissue. Studies in animals have demonstrated decreased flow heterogeneity (FH) in states of decreased perfusion pressure. The purpose of the present study was to assess microvascular FH changes in acute stroke with use of a novel perfusion-weighted MRI technique and to evaluate the ability of combined diffusion-weighted MRI and FH measurements to predict final infarct size. METHODS Cerebral blood flow, FH, and plasma mean transit time (MTT) were measured in 11 patients who presented with acute (<12 hours after symptom onset) stroke. Final infarct size was determined with follow-up MRI or CT scanning. RESULTS In normal brain tissue, the distribution of relative flows was markedly skewed toward high capillary flow velocities. Within regions of decreased cerebral blood flow, plasma MTT was prolonged. Furthermore, subregions were identified with significant loss of the high-flow component of the flow distribution, thereby causing increased homogeneity of flow velocities. In parametric maps that quantify the acute deviation of FH from that of normal tissue, areas of extreme homogenization of capillary flows predicted final infarct size on follow-up scans of 10 of 11 patients. CONCLUSIONS Flow heterogeneity and MTT can be rapidly assessed as part of a routine clinical MR examination and may provide a tool for planning of individual stroke treatment, as well as in targeting and evaluation of emerging therapeutic strategies.

Determination of the quantitative thallium imaging variables that optimize detection of coronary artery disease

Although quantification of exercise thallium images has been previously reported, the relative value of different imaging variables for detection of coronary artery disease has not been analyzed in a large group of patients with cardiac catheterization data. Regional initial thallium uptake, redistribution and clearance on thallium study were measured in 325 patients also undergoing cardiac catheterization (281 patients with and 44 patients without coronary artery disease). Normal values were defined in 55 other clinically normal subjects. When five myocardial segments were analyzed in each view, the respective values for sensitivity and specificity were 95 and 50% for initial thallium uptake, 60 and 87% for redistribution and 74 and 66% for clearance. Initial thallium uptake was the most sensitive but least specific (p less than 0.001), whereas redistribution was the least sensitive and most specific (p less than 0.001). Using stepwise logistic regression analysis, the best correlate of coronary artery disease was initial thallium uptake. Addition of redistribution to a mathematical model of the probability of coronary artery disease did not alter sensitivity, but increased specificity from 50 to 70% (p less than 0.001). Once initial uptake and redistribution were considered, myocardial thallium clearance provided no additional improvement in the correlation. Excluding the two basal segments in each view from the analysis increased the specificity from 70 to 80% (p less than 0.001) without affecting sensitivity. Of the 15 patients (5%) with coronary disease not detected using this approach, none had left main disease and 10 (67%) had one vessel disease. A combination of variables derived from quantification of exercise thallium images provides a superior sensitivity and specificity for the detection of coronary artery disease compared with the use of a single variable.

Attenuation correction in gamma emission computed tomography

Abstract: A method of computing tomographic images from single photon radionuclide emission data is presented. The method takes into account attenuation of gamma rays inside the source and makes use of an iterative technique, based on the difference between the projection data obtained from the source and computed projections, called reprojections, from successive reconstructions of the sources. The method has been tested both by computer simulations and reconstruction of plastic phantoms imaged with 99mTc radionuclides. Substantial improvement in reconstruction accuracy over algorithms uncorrected for internal attenuation is demonstrated. Since the technique is iterative, it can be used with a variety of reconstruction algorithms or combined with other first approximation techniques of attenuation correction.

Ripple suppression during reconstruction in transverse tomography

In many transverse tomographic reconstruction techniques using filtered projections, ringing or overshoot is found in regions where the attenuation suddenly changes. Such artifacts are obtained with iterative techniques as well. A modification of the filter often used in corrected tomography which produces very small ringing in such regions is described. A transmission device such as the EMI Scanner was computer simulated. Projections were calculated through a phantom and the phantom then reconstructed using different filters. The phantom consisted of a disc with diameter 32 arbitrary units centred at the origin of the plane.