Publications by Year: 2003

A new method for mapping diffusivity profiles in tissue is presented. The Bloch-Torrey equation is modified to include a diffusion term with an arbitrary rank Cartesian tensor. This equation is solved to give the expression for the generalized Stejskal-Tanner formula quantifying diffusive attenuation in complicated geometries. This makes it possible to calculate the components of higher-rank tensors without using the computationally-difficult spherical harmonic transform. General theoretical relations between the diffusion tensor (DT) components measured by traditional (rank-2) DT imaging (DTI) and 3D distribution of diffusivities, as measured by high angular resolution diffusion imaging (HARDI) methods, are derived. Also, the spherical tensor components from HARDI are related to the rank-2 DT. The relationships between higher- and lower-rank Cartesian DTs are also presented. The inadequacy of the traditional rank-2 tensor model is demonstrated with simulations, and the method is applied to excised rat brain data collected in a spin-echo HARDI experiment.

Several techniques have been described in the literature in recent years for the reconstruction of a regular volume out of a series of ultrasound (US) slices with arbitrary orientations, typically scanned by means of US freehand systems. However, a systematic approach to such a problem is still missing. This paper focuses on proposing a theoretical framework for the 3-D US volume reconstruction problem. We introduce a statistical method for the construction and trimming of the sampling grid where the reconstruction will be carried out. The results using in vivo US data demonstrate that the computed reconstruction grid that encloses the region-of-interest (ROI) is smaller than those obtained from other reconstruction methods in those cases where the scanning trajectory deviates from a pure straight line. In addition, an adaptive Gaussian interpolation technique is studied and compared with well-known interpolation methods that have been applied to the reconstruction problem in the past. We find that the proposed method numerically outperforms former proposals in several control studies; subjective visual results also support this conclusion and highlight some potential deficiencies of methods previously proposed.

Diffusion Tensor MRI (DT-MRI) can provide important in vivo information for the detection of brain abnormalities in diseases characterized by compromised neural connectivity. To quantify diffusion tensor abnormalities based on voxel-based statistical analysis, spatial normalization is required to minimize the anatomical variability between studied brain structures. In this article, we used a multiple input channel registration algorithm based on a demons algorithm and evaluated the spatial normalization of diffusion tensor image in terms of the input information used for registration. Registration was performed on 16 DT-MRI data sets using different combinations of the channels, including a channel of T2-weighted intensity, a channel of the fractional anisotropy, a channel of the difference of the first and second eigenvalues, two channels of the fractional anisotropy and the trace of tensor, three channels of the eigenvalues of the tensor, and the six channel tensor components. To evaluate the registration of tensor data, we defined two similarity measures, i.e., the endpoint divergence and the mean square error, which we applied to the fiber bundles of target images and registered images at the same seed points in white matter segmentation. We also evaluated the tensor registration by examining the voxel-by-voxel alignment of tensors in a sample of 15 normalized DT-MRIs. In all evaluations, nonlinear warping using six independent tensor components as input channels showed the best performance in effectively normalizing the tract morphology and tensor orientation. We also present a nonlinear method for creating a group diffusion tensor atlas using the average tensor field and the average deformation field, which we believe is a better approach than a strict linear one for representing both tensor distribution and morphological distribution of the population.

RATIONALE AND OBJECTIVES: A novel method to correct flow data from magnetic resonance phase contrast (MR-PC) angiography, based on combining computational fluid dynamics and segmentation in a level set framework, was developed and tested in this study. MATERIALS AND METHODS: The MR-PC velocity data was used in a partial differential equation-based level set method for vessel segmentation. The results were supplied as the quantitative description of the vessel wall to the flow field solver using computational fluid dynamics, based on the level set method, to obtain a physically meaningful flow. The most significant characteristic of our novel approach is that it requires light computational loads, especially insofar as it avoids generation of complex computational grid system. The integration of segmentation and computational fluid dynamics in a level set framework is shown to be both robust and economic, and yet yields a physically correct velocity field and optimal vessel geometry. RESULTS: The application to the flow field in a straight tube with circular cross section of constant radius demonstrated the validity of out new approach, especially the treatment of the velocity boundary conditions on the solid wall. Simulation of the velocity field in both common carotid artery and bifurcation of basilar and vertebral arteries, based on clinical MR-PC data, provided with smooth and stable results. CONCLUSION: Applying this procedure to both synthetic and clinical data, significant improvement of the blood velocity field, such as a smooth velocity distribution aligned along the vessels and removal of burst or error vectors, could be observed. This procedure also offers possibilities for improved vessel segmentation.

Evidence suggests that a disruption in limbic system network integrity and, in particular, the cingulate gyrus (CG), may play a role in the pathophysiology of schizophrenia; however, the cingulum bundle (CB), the white matter tract furnishing both input and output to CG, and the most prominent white matter fiber tract in the limbic system, has not been evaluated in schizophrenia using the new technology of diffusion tensor imaging (DTI). We used line scan DTI to evaluate diffusion in the CB in 16 male schizophrenia patients and 18 male control subjects, group-matched for age, parental socioeconomic status, and handedness. We acquired 4-mm-thick coronal slices through the entire brain. Maps of fractional anisotropy (FA) were generated to quantify diffusion within the left and right CB on eight slices that included the central portion of the CB. Results showed group differences, bilaterally, in area and mean FA for CB, where patients showed smaller area and less anisotropy than controls. For patients, decreased left CB correlated significantly with attention and working memory measures as assessed by the Wisconsin Card Sorting Test. These data provide strong evidence for CB disruptions in schizophrenia, which may be related to disease-related attention and working memory abnormalities.

OBJECTIVE: To investigate the geometry of the glenohumeral joint using three-dimensional (3D) models; define landmarks, planes, angles and regions of interest; and analyze the exact morphology of the scapula. MATERIALS AND METHODS: We reconstructed 3D scapula models based on computed tomography (CT) data sets of 12 healthy controls. Three-dimensional models were reconstructed using the 3D Slicer2 (Surgical Planning Lab, Boston, MA), which provides interactive measurement. The 3D model and measuring tools can be freely shifted and rotated in all planes.

RATIONALE AND OBJECTIVES: The accuracy of diagnostic test and imaging segmentation is important in clinical practice because it has a direct impact on therapeutic planning. Statistical validations of classification accuracy was conducted based on parametric receiver operating characteristic analysis, illustrated on three radiologic examples, MATERIALS AND METHODS: Two parametric models were developed for diagnostic or imaging data. Example 1: A semi-automated fractional segmentation algorithm was applied to magnetic resonance imaging of nine cases of brain tumors. The tumor and background pixel data were assumed to have bi-beta distributions. Fractional segmentation was validated against an estimated composite pixel-wise gold standard based on multi-reader manual segmentations. Example 2: The predictive value of 100 cases of spiral computed tomography of ureteral stone sizes, distributed as bi-normal after a non-linear transformation, under two treatment options received. Example 3: One hundred eighty cases had prostate-specific antigen levels measured in a prospective clinical trial. Radical prostatectomy was performed in all to provide a binary gold standard of local and advanced cancer stages. Prostate-specific antigen level was transformed and modeled by bi-normal distributions. In all examples, areas under the receiver operating characteristic curves were computed. RESULTS. The areas under the receiver operating characteristic curves were: Example 1: Fractional segmentation of magnetic resonance imaging of brain tumors: meningiomas (0.924-0.984); astrocytomas (0.786-0.986); and other low-grade gliomas (0.896-0.983). Example 3: Ureteral stone size for treatment planning (0.813). Example 2: Prostate-specific antigen for staging prostate cancer (0.768). CONCLUSION: All clinical examples yielded fair to excellent accuracy. The validation metric area under the receiver operating characteristic curves may be generalized to evaluating the performances of several continuous classifiers related to imaging.

OBJECTIVE: An abnormal superior temporal gyrus has figured prominently in schizophrenia research, and left superior temporal gyrus volume has been shown to be smaller in male subjects with schizotypal personality disorder. This is the first structural magnetic resonance imaging study to examine a group of female subjects with schizotypal personality disorder.

Schizophrenia has long been associated with abnormal patterns of arousal that are thought to reflect disturbances in the reticular-activating system of the brain. Psychophysiological investigations of sensory responsivity have repeatedly demonstrated reduced reactivity and habituation to moderately intense stimuli in patients with schizophrenia. While not traditionally used as a measure of physiological arousal, the startle reflex represents an alternative method for studying reactivity and habituation in schizophrenia. This study examined eye blink responsivity to a repeatedly presented intense acoustic startle probe in men with chronic schizophrenia and healthy normal controls. Subjects’ personality profiles were also measured, as increased reactivity and arousal have been traditionally implicated as a physiological component to the personality trait of neuroticism. Results indicated that schizophrenic subjects did demonstrate significantly reduced rates of habituation to the acoustic startle probe and higher scores on measures of neuroticism in comparison to normal controls. However, no correlation between habituation rate and neuroticism emerged. These studies replicate previous findings of habituation in schizophrenia and provide further evidence for sensory reactivity disturbances in schizophrenia. The relationship of these findings to cognitive disturbances in schizophrenia is considered and directions for future research are discussed.

BACKGROUND: The Heschl gyrus and planum temporale have crucial roles in auditory perception and language processing. Our previous investigation using magnetic resonance imaging (MRI) indicated smaller gray matter volumes bilaterally in the Heschl gyrus and in left planum temporale in patients with first-episode schizophrenia but not in patients with first-episode affective psychosis. We sought to determine whether there are progressive decreases in anatomically defined MRI gray matter volumes of the Heschl gyrus and planum temporale in patients with first-episode schizophrenia and also in patients with first-episode affective psychosis. METHODS: At a private psychiatric hospital, we conducted a prospective high spatial resolution MRI study that included initial scans of 28 patients at their first hospitalization (13 with schizophrenia and 15 with affective psychosis, 13 of whom had a manic psychosis) and 22 healthy control subjects. Follow-up scans occurred, on average, 1.5 years after the initial scan. RESULTS: Patients with first-episode schizophrenia showed significant decreases in gray matter volume over time in the left Heschl gyrus (6.9%) and left planum temporale (7.2%) compared with patients with first-episode affective psychosis or control subjects. CONCLUSIONS: These findings demonstrate a left-biased progressive volume reduction in the Heschl gyrus and planum temporale gray matter in patients with first-episode schizophrenia in contrast to patients with first-episode affective psychosis and control subjects. Schizophrenia but not affective psychosis seems to be characterized by a postonset progression of neocortical gray matter volume loss in the left superior temporal gyrus and thus may not be developmentally fixed.