Search

In this paper, we propose a coupled level set (LS) framework for segmentation of bladder wall using T1-weighted magnetic resonance (MR) images with clinical applications to virtual cystoscopy (i.e., MR cystography). The framework uses two collaborative LS functions and a regional adaptive clustering algorithm to delineate the bladder wall for the wall thickness measurement on a voxel-by-voxel basis. It is significantly different from most of the pre-existing bladder segmentation work in four aspects. First of all, while most previous work only segments the inner border of the wall or at most manually segments the outer border, our framework extracts both the inner and outer borders automatically except that the initial seed point is given by manual selection. Secondly, it is adaptive to T1-weighted images with decreased intensities in urine, as opposed to enhanced intensities in T1-weighted scenario and computed tomography. Thirdly, by considering the image global intensity distribution and local intensity contrast, the defined image energy function in the framework is more immune to inhomogeneity effect, motion artifacts and image noise. Finally, the bladder wall thickness is measured by the length of integral path between the two borders which mimic the electric field line between two iso-potential surfaces. The framework was tested on six datasets with comparison to the well-known Chan-Vese (C-V) LS model. Five experts blindly scored the segmented inner and outer borders of the presented framework and the C-V model. The scores demonstrated statistically the improvement in detecting the inner and outer borders. [ABSTRACT FROM AUTHOR]

This work presents an iterative expectation-maximization (EM) approach to the maximum a posteriori (MAP) solution of segmenting tissue mixtures inside each image voxel. Each tissue type is assumed to follow a normal distribution across the field-of-view (FOV). Furthermore, all tissue types are assumed to be independent from each other. Under these assumptions, the summation of all tissue mixtures inside each voxel leads to the image density mean value at that voxel. The summation of all the tissue mixtures' unobservable random processes leads to the observed image density at that voxel, and the observed image density value also follows a normal distribution (image data are observed to follow a normal distribution in many applications). By modeling the underlying tissue distributions as a Markov random field across the FOV, the conditional expectation of the posteriori distribution of the tissue mixtures inside each voxel is determined, given the observed image data and the current-iteration estimation of the tissue mixtures. Estimation of the tissue mixtures at next iteration is computed by maximizing the conditional expectation. The iterative EM approach to a MAP solution is achieved by a finite number of iterations and reasonable initial estimate. This MAP-EM framework provides a theoretical solution to the partial volume effect, which has been a major cause of quantitative imprecision in medical image processing. Numerical analysis demonstrated its potential to estimate tissue mixtures accurately and efficiently. [ABSTRACT FROM AUTHOR]

Computed tomography colonography (CTC) or CT-based virtual colonoscopy (VC) is an emerging tool for detection of colonic polyps. Compared to the conventional fiber-optic colonoscopy, VC has demonstrated the potential to become a mass screening modality in terms of safety, cost, and patient compliance. However, current CTC delivers excessive X-ray radiation to the patient during data acquisition. The radiation is a major concern for screening application of CTC. In this work, we performed a simulation study to demonstrate a possible ultra low-dose CT technique for VC. The ultra low-dose abdominal CT images were simulated by adding noise to the sinograms of the patient CTC images acquired with normal dose scans at 100 mAs levels. The simulated noisy sinogram or projection data were first processed by a Karhunen-Loëve domain penalized weighted least-squares (KL-PWLS) restoration method and then reconstructed by a filtered backprojection algorithm for the ultra low-dose CT images. The patient-specific virtual colon lumen was constructed and navigated by a VC system after electronic colon cleansing of the orally-tagged residue stool and fluid. By the KL-PWLS noise reduction, the colon lumen can successfully be constructed and the colonic polyp can be detected in an ultra low-dose level below 50 mAs. Polyp detection can be found more easily by the KL-PWLS noise reduction compared to the results using the conventional noise filters, such as Hanning filter. These promising results indicate the feasibility of an ultra low-dose CTC pipeline for colon screening with less-stressful bowel preparation by fecal tagging with oral contrast. [ABSTRACT FROM AUTHOR]

In this paper, we present a new hardware acceleration method to speedup the ordered-subsets expectation-maximization (OS-EM) algorithm for quantitative single photon emission computed tomography (SPECT) image reconstruction with varying focal-length fan-beam (VFF) collimation. By utilizing the geometrical symmetry of VFF point-spread function (PSF), compensation for object-specific attenuation and system-specific PSF are accelerated using currently available PC video/graphics card technologies. A tenfold acceleration of quantitative SPECT reconstruction is achieved. [ABSTRACT FROM AUTHOR]

When excessive quantum noise is present in extremely low dose X-ray CT imaging, statistical properties of the data has to be considered to achieve a satisfactory image reconstruction. Statistical Iterative reconstruction with accurate modeling of the noise, rather than a filtered back-projection (FBP) with low-pass filtering, is one way to deal with the problem. Estimating a noise-free sinogram to satisfy the FBP reconstruction for the Radon transform is another way. The benefits of the latter include a higher computation efficiency, more uniform spatial resolution in the reconstructed image, and less modification of the current machine configurations. In a clinic X-ray CT system, the acquired raw data must be calibrated, in addition to the logarithmic transform, to achieve the high diagnostic image quality. The calibrated projection data or sinogram no longer follow a compound Poisson distribution in general, but are close to a Gaussian distribution with signal-dependent variance. In this paper, we first investigated a relatively accurate statistical model for the sinogram data, based on several phantom experiments. Then we developed a penalized likelihood method to smooth the sinogram, which led to a set of nonlinear equations that can be solved by iterated conditional mode (ICM) algorithm within a reasonable computing time. The method was applied to several experimental datasets acquired at 120 kVp, 10 mA/20 mA/50 mA protocols with a GE HiSpeed multi-slice detector CT scanner and demonstrated a significant noise suppression without noticeable sacrifice of the spatial resolution. [ABSTRACT FROM AUTHOR]

We present an automatic mixture-based algorithm for segmentation of brain tissues (white and gray matters--WM and GM), cerebral spinal fluid (CSF), and brain lesions to quantitatively analyze multiple sclerosis. The method performs intensity-based tissue classification using multispectral magnetic resonance (MR) images based on a stochastic model. With the existence of white Gaussian noise and spatially invariant blurring in acquired MR images, a Karhunen-Loéve (K-L) domain Wiener filter is applied for accurate noise reduction and resolution restoration on blurred and noisy images to minimize the partial volume effect (PVE), which is a major limiting factor for the quantitative analysis. Following that, we utilize a Markov random field Gibbs model to integrate the local spatial information into the well-established expectation-maximization model-fitting algorithm. Each voxel is then classified by a maximum a posterior (MAP) criterion, indicating its probabilities of belonging to each class, i.e., each voxel is labeled as a mixel with different tissue percentages, leading to further minimization of the PVE. The volumes of WM, GM, CSF, and brain lesions are extracted from the mixture-based segmentation and the corresponding brain atrophies are computed. In this study, we have investigated the accuracy and repeatability of the algorithm with inclusion of noise analysis and point spread function for image resolution enhancement. Experimental results on phantom, healthy volunteer, and patient studies are presented. [ABSTRACT FROM AUTHOR]