NMR Object Boundaries: B-Spline Modeling and Estimator Performance

We give estimation error bounds and specify optimal estimators for continuous, closed boundary curves in an NMR image. The boundary is parameterized using periodic B-splines. A Cramer-Rao lower bound on mean-square-estimate error in the presence of system smoothing and Gaussian noise is derived, and the performance of maximum likelihood and penalized maximum likelihood estimators is compared to this bound. Finally, we comment on the usefulness of estimates of the boundary for providing anatomical side information in the reconstruction of functional tomographic images like those of a PET or SPECT system.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/85835/1/Fessler133.pd

Minimax Emission Computed Tomography using High-Resolution Anatomical Side Information and B-Spline Models

In this paper a minimax methodology is presented for combining information from two imaging modalities having different intrinsic spatial resolutions. The focus application is emission computed tomography (ECT), a low-resolution modality for reconstruction of radionuclide tracer density, when supplemented by high-resolution anatomical boundary information extracted from a magnetic resonance image (MRI) of the same imaging volume. The MRI boundary within the two-dimensional (2-D) slice of interest is parameterized by a closed planar curve. The Cramer-Rao (CR) lower bound is used to analyze estimation errors for different boundary shapes. Under a spatially inhomogeneous Gibbs field model for the tracer density a representation for the minimax MRI-enhanced tracer density estimator is obtained. It is shown that the estimator is asymptotically equivalent to a penalized maximum likelihood (PML) estimator with resolution-selective Gibbs penalty. Quantitative comparisons are presented using the iterative space alternating generalized expectation maximization (SAGE-FM) algorithm to implement the PML estimator with and without minimax weight averaging.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/85822/1/Fessler86.pd

Adaptive morphological filters based on a multiple orientation vector field dependent on image local features

This paper addresses the formulation of adaptive morphological filters based on spatially-variant structuring elements. The adaptivity of these filters is achieved by modifying the shape and orientation of the structuring elements according to a multiple orientation vector field. This vector field is provided by means of a bank of directional openings which can take into account the possible multiple orientations of the contours in the image. After reviewing and formalizing the definition of the spatially-variant dilation, erosion, opening and closing, the proposed structuring elements are described. These spatially-variant structuring elements are based on ellipses which vary over the image domain adapting locally their orientation according to the multiple orientation vector field and their shape (the eccentricity of the ellipses) according to the distance to relevant contours of the objects. The proposed adaptive morphological filters are used on gray-level images and are compared with spatially-invariant filters, with spatially-variant filters based on a single orientation vector field, and with adaptive morphological bilateral filters. Results show that the morphological filters based on a multiple orientation vector field are more adept at enhancing and preserving structures which contains more than one orientation

Blind deconvolution of medical ultrasound images: parametric inverse filtering approach

©2007 IEEE. Personal use of this material is permitted. However, permission to reprint/republish this material for advertising or promotional purposes or for creating new collective works for resale or distribution to servers or lists, or to reuse any copyrighted component of this work in other works must be obtained from the IEEE. This material is presented to ensure timely dissemination of scholarly and technical work. Copyright and all rights therein are retained by authors or by other copyright holders. All persons copying this information are expected to adhere to the terms and constraints invoked by each author's copyright. In most cases, these works may not be reposted without the explicit permission of the copyright holder.DOI: 10.1109/TIP.2007.910179The problem of reconstruction of ultrasound images by means of blind deconvolution has long been recognized as one of the central problems in medical ultrasound imaging. In this paper, this problem is addressed via proposing a blind deconvolution method which is innovative in several ways. In particular, the method is based on parametric inverse filtering, whose parameters are optimized using two-stage processing. At the first stage, some partial information on the point spread function is recovered. Subsequently, this information is used to explicitly constrain the spectral shape of the inverse filter. From this perspective, the proposed methodology can be viewed as a ldquohybridizationrdquo of two standard strategies in blind deconvolution, which are based on either concurrent or successive estimation of the point spread function and the image of interest. Moreover, evidence is provided that the ldquohybridrdquo approach can outperform the standard ones in a number of important practical cases. Additionally, the present study introduces a different approach to parameterizing the inverse filter. Specifically, we propose to model the inverse transfer function as a member of a principal shift-invariant subspace. It is shown that such a parameterization results in considerably more stable reconstructions as compared to standard parameterization methods. Finally, it is shown how the inverse filters designed in this way can be used to deconvolve the images in a nonblind manner so as to further improve their quality. The usefulness and practicability of all the introduced innovations are proven in a series of both in silico and in vivo experiments. Finally, it is shown that the proposed deconvolution algorithms are capable of improving the resolution of ultrasound images by factors of 2.24 or 6.52 (as judged by the autocorrelation criterion) depending on the type of regularization method used