This paper presents a new variational method for the segmentation of a moving object against a still background, over a sequence of [two-dimensional or three-dimensional (3-D)] image frames. The method is illustrated in application to myocardial gated single photon emission computed tomography (SPECT) data, and incorporates a level set framework to handle topological changes while providing closed boundaries. The key innovation is the introduction of a geometrical constraint into the derivation of the Euler-Lagrange equations, such that the segmentation of each individual frame can be interpreted as a closed boundary of an object (an isolevel of a set of hyper-surfaces) while integrating information over the entire sequence. This results in the definition of an evolution velocity normal to the object boundary. Applying this method to 3-D myocardial gated SPECT sequences, the left ventricle endocardial and epicardial limits can be computed in each frame. This space-time segmentation method was tested on simulated and clinical 3-D myocardial gated SPECT sequences and the corresponding ejection fractions were computed.

The qualitative and quantitative accuracy of SPECT images is degraded by physical factors of attenuation, Compton scatter and spatially varying collimator geometric response. This paper presents a 3D ray-tracing technique for modelling attenuation, scatter and geometric response for SPECT imaging in an inhomogeneous attenuating medium. The model is incorporated into a three-dimensional projector-backprojector and used with the maximum-likelihood expectation-maximization algorithm for reconstruction of parallel-beam data. A transmission map is used to define the inhomogeneous attenuating and scattering object being imaged. The attenuation map defines the probability of photon attenuation between the source and the scattering site, the scattering angle at the scattering site and the probability of attenuation of the scattered photon between the scattering site and the detector. The probability of a photon being scattered through a given angle and being detected in the emission energy window is approximated using a Gaussian function. The parameters of this Gaussian function are determined using physical measurements of parallel-beam scatter line spread functions from a non-uniformly attenuating phantom. The 3D ray-tracing scatter projector-backprojector produces the scatter and primary components. Then, a 3D ray-tracing projector-backprojector is used to model the geometric response of the collimator. From Monte Carlo and physical phantom experiments, it is shown that the best results are obtained by simultaneously correcting attenuation, scatter and geometric response, compared with results obtained with only one or two of the three corrections. It is also shown that a 3D scatter model is more accurate than a 2D model. A transmission map is useful for obtaining measurements of attenuation and scatter in SPECT data, which can be used together with a model of the geometric response of the collimator to obtain corrected images with quantitative and diagnostically accurate information.

No abstract

This paper presents a new variational method for the segmentation of a moving object against a still background, over a sequence of (2-D or 3-D) image frames. The method is illustrated in application to myocardial gated SPECT data, and incorporates a level set framework to handle topological changes while providing closed boundaries. The key innovation is the introduction of a geometrical constraint into the derivation of the Euler-Lagrange equations, such that the segmentation of each individual frame can be interpreted as a closed boundary of an object (an isolevel of a set of hyper-surfaces) while integrating information over the entire sequence. This results in the definition of an evolution velocity normal to the object boundary. Applying this method to 3-D myocardial gated SPECT sequences, the left ventricle endocardial and epicardial limits can be computed in each frame. This space-time segmentation method was tested on simulated and clinical 3-D myocardial gated SPECT sequences and the corresponding ejection fractions were computed.

SPECT systems incorporate the use of one or more rotating gamma cameras which can be equipped with cone-beam collimators to improve the trade-off between spatial resolution and sensitivity. The geometry of the cone-beam collimators implies that a specific 3D reconstruction algorithm must be applied. Algebraic methods provide the possibility of including the physical characteristics, such as attenuation, Compton scatter and detector response, in the reconstruction process. However, the reconstruction problem is an ill-posed problem which should be regularized. This paper presents a 3D algebraic method that combines three regularizing constraints. These constraints deal respectively with penalizing negative voxels, local noise smoothing and missing data compensation. The results presented were obtained from imaging simulations, phantom data and from a thyroid clinical study of a normal volunteer.