Seven-pinhole tomography of the thyroid

In this paper we introduce a new application of the sevenpinhole (7P) collimator: tomographic imaging of the thyroid. The collimator design has been reoptimized for this particular application by diminishing the distance from the collimator plate to the crystal and by choosing a smaller pinhole diameter. To reconstruct thyroid images from the two-dimensional projection data we use a method which we developed for 7P tomographic imaging of the heart [1]. Phantom experiments and patient studies demonstrate that this new device is capable of producing tomographic images of good quality and high resolution. Therefore, it seems to offer a promising alternative to conventional planar imaging of the thyroid (using a single-pinhole collimator).     

Time-coded aperture tomography: experimental results

In this paper, we discuss the applicability of a time-coded aperture system especially designed for thyroid tomography on the basis of phantom experiments. Our studies show that 1) the quality of the reconstructions is high (e.g., a cold spot of 6 mm diameter in a thyroid phantom can easily be detected), and 2) the reconstruction can be carried out in less than 11 min on a standard 16 bit minicomputer (HP1000). It is therefore concluded that the clinical potentiality of the device is good.     

Fast SPECT simulation including object shape dependent scatter

A fast simulator of SPECT projection data taking into account attenuation, distance dependent detector response, and scatter has been developed, based on an analytical point spread function model. The parameters of the scatter response are obtained from a single line source measurement with a triangular phantom. The simulator is able to include effects of object curvature on the scatter response to a high accuracy. The simulator has been evaluated for homogeneous media by measurements of (99m)Tc point sources placed at different locations in a water-filled cylinder at energy windows of 15% and 20%. The asymmetrical shapes of measured projections of point sources are In excellent agreement with simulations for both energy windows. Scatter-to-primary ratio (SPR) calculations of point sources at different positions in a cylindrical phantom differ not more than a few percent from measurements. The simulator uses just a few megabytes of memory for storing the tables representing the forward model; furthermore, simulation of 60 SPECT projections from a three-dimensional digital brain phantom with 6-mm cubic voxels takes only ten minutes on a standard workstation. Therefore, the simulator could serve as a projector in iterative true 3-D SPECT reconstruction.     

Improved SPECT quantitation using fully three-dimensional iterative spatially variant scatter response compensation

The quality and quantitative accuracy of iteratively reconstructed SPECT images improves when better point spread function (PSF) models of the gamma camera are used during reconstruction. Here, inclusion in the PSF model of photon crosstalk between different slices caused by limited gamma camera resolution and scatter is examined. A three-dimensional (3-D) projector back-projector (proback) has been developed which models both the distance dependent detector point spread function and the object shape-dependent scatter point spread function of single photon emission computed tomography (SPECT). A table occupying only a few megabytes of memory is sufficient to represent this scatter model. The contents of this table are obtained by evaluating an analytical expression for object shape-dependent scatter. The proposed approach avoids the huge memory requirements of storing the full transition matrix needed for 3-D reconstruction including object shape-dependent scatter. In addition, the method avoids the need for lengthy Monte Carlo simulations to generate such a matrix. In order to assess the quantitative accuracy of the method, reconstructions of a water filled cylinder containing regions of different activity levels and of simulated 3-D brain projection data have been evaluated for technetium-99m. It is shown that fully 3-D reconstruction including complete detector response and object shape-dependent scatter modeling clearly outperforms simpler methods that lack a complete detector response and/or a complete scatter response model. Fully 3-D scatter correction yields the best quantitation of volumes of interest and the best contrast-to-noise curves.     

The Gaussian scale-space paradigm and the multiscale local jet

A representation of local image structure is proposed which takes into account both the image's spatial structure at a given location, as well as its deep structure, that is, its local behaviour as a function of scale or resolution (scale-space). This is of interest for several low-level image tasks. The proposed basis of scale-space, for example, enables a precise local study of interactions of neighbouring image intensities in the course of the blurring process. It also provides an extrapolation scheme for local image data, obtained at a given spatial location and resolution, to a finite scale-space neighbourhood. This is especially useful for the determination of sampling rates and for interpolation algorithms in a multilocal context. Another, particularly straightforward application is image enhancement or deblurring, which is an instance of data extrapolation in the high-resolution direction.A potentially interesting feature of the proposed local image parametrisation is that it captures a trade-off between spatial and scale extrapolations from a given interior point that do not exceed a given tolerance. This (rade-off suggests the possibility of a fairly coarse scale sampling at the expense of a dense spatial sampling large relative spatial overlap of scale-space kernels).The central concept developed in this paper is an equivalence class called the multiscale local jet, which is a hierarchical, local characterisation of the image in a full scale-space neighbourhood. For this local jet, a basis of fundamental polynomials is constructed that captures the scale-space paradigm at the local level up to any given order.     

Adaptive Stochastic Gradient Descent Optimisation for Image Registration

We present a stochastic gradient descent optimisation method for image registration with adaptive step size prediction. The method is based on the theoretical work by Plakhov and Cruz (J. Math. Sci. 120(1):964–973, 2004). Our main methodological contribution is the derivation of an image-driven mechanism to select proper values for the most important free parameters of the method. The selection mechanism employs general characteristics of the cost functions that commonly occur in intensity-based image registration. Also, the theoretical convergence conditions of the optimisation method are taken into account. The proposed adaptive stochastic gradient descent (ASGD) method is compared to a standard, non-adaptive Robbins-Monro (RM) algorithm. Both ASGD and RM employ a stochastic subsampling technique to accelerate the optimisation process. Registration experiments were performed on 3D CT and MR data of the head, lungs, and prostate, using various similarity measures and transformation models. The results indicate that ASGD is robust to these variations in the registration framework and is less sensitive to the settings of the user-defined parameters than RM. The main disadvantage of RM is the need for a predetermined step size function. The ASGD method provides a solution for that issue.     

Three-dimensional simulation of time-coded emission tomography

In this paper we discuss the potentialities of a time-coded single photon emissive imaging system for thyroid tomography. We have performed three-dimensional simulation studies in order to answer two questions: 1) does this coded aperture device produce good quality reconstructions, and 2) can the reconstruction be carried out sufficiently fast on a microcomputer. Our study leads to the conclusion that both questions can be answered affirmatively. Hence, time-coded emission tomography is a potentially useful imaging technique for diagnostic clinical practice.     

Objective quantification of the motion of soft tissues in the orbit

Orbital soft-tissue motion analysis aids in the localization and diagnosis of orbital disorders. A technique has been developed to objectively quantify and visualize motion in the orbit during gaze. T1-weighted MR volume sequences are acquired during gaze and soft-tissue motion is quantified using optical flow techniques. The flow field is visualized using color-coding: orientation of the flow vector is coded by hue and magnitude by saturation of the pixel. Current clinical circumstances limit MR image acquisition to short sequences and short acquisition times. The effect of these limitations on the performance of optical flow computation has been studied for four representative optical flow algorithms: on short (nine frames) and long (21 frames) simulated sequences of rotation of a magnetic resonance (MR) imaged object, on short measured MR sequences of controlled rotation of the same object and on short MR sequences of motion in the orbit. On the short simulated and motion-controlled sequences, the Lucas and Kanade algorithm showed the best performance with respect to both accuracy and robustness. These motion estimates were accurate to within 20%. Motion in the orbit ranged between 0.05 and 0.25 mm/degree gaze. Color-coding was found to be attractive as a visualization technique, because it shows both magnitude and orientation of all flow vectors without cluttering.     

Medical image matching-a review with classification

A classification scheme for multimodal image matching is considered. The scope of the classification is restricted to methods that register data after acquisitions. The classification scheme may be used for any modality; not only for (2-D) projection images and (3-D) tomographic images, but also for other signal modalities that provide spatial insight into function or anatomy, e.g., EEG (electroencephalography) or MEG (magnetoencephalography) and for the real physical patient. The available literature on image matching is discussed and classified     

The stretching nonlinearity of the basilar membrane in a cochlear model

Recently it has been suggested that the stretching nonlinearity of the basilar membrane might be responsible for the observed nonlinear behaviour of basilar membrane motion. In the present study this type of nonlinearity is investigated, both by estimating its influence in an analytical manner, and by calculating its effect numerically, using a regular pertubation method. The conclusion reads that the stretching nonlinearity does not explain the observed nonlinear phenomena; not only is stretching negligible at normal sound levels, but it also fails to fit the data qualitatively, because of its typical hard-spring effect. In consequence, the origin of cochlear nonlinearity is not to be sought in the macromechanics of the inner ear, but in the more detailed processes in th organ of Corti-tectorial membrane complex.