PubMed Health⌕ Search

Biomedical subjects

M Defrise

Publications and source records attributed to M Defrise.

29 records · Page 2Linked to original sources

Three-dimensional multimodality medical image registration using a parameter accumulation approach.

A parameter accumulation method based on the Hough transformation is proposed to register three-dimensional (3-D) multimodality medical images. The estimation of registration parameters is decomposed into separate estimations of rotation, using directional vectors, and translation, using positional vectors. Similarly, the rotation parameters are decomposed into the rotation axis and angle, which are then estimated separately. This kind of decomposition reduces the parametric dimension and improves the computing efficiency which has been a major concern in implementing the Hough transformation. When 3-D rotation is involved, evaluating registration error is not straightforward. This paper introduces an equivalent error angle as a criterion to evaluate the performance of 3-D registration methods. Experimental results indicate that a least-squares fitting is superior to the parameter accumulation with data contaminated by additive noise only. When mismatched feature points (outliers) exist in the data set, however, the parameter accumulation approach is more accurate. The application of the proposed approach to the registration of 3-D PET and CT images is demonstrated.

Algorithms↗

Filtered backprojection reconstruction of combined parallel beam and cone beam SPECT data.

The reconstruction problem for a combined parallel beam (PB) and cone beam (CB) imaging geometry has been addressed. The general algorithm (CB-FBP) of Defrise and Clack has been applied to this geometry and shown to provide accurate images as expected. A second algorithm specifically tailored to the PB-CB geometry was developed. It uses the general principles of the CB-FBP method to combine a shift-variant filtering of the PB data with a standard reconstruction of the CB data using the algorithm of Feldkamp, Davis and Kress (FDK). This 'mixed' algorithm has the advantage of fewer interpolation steps, thereby reducing reconstruction time and providing more accurate reconstructions. The algorithms were applied to noiseless data from a computer-generated 3D Shepp phantom. Both the CB-FBP algorithm and the mixed algorithm successfully combine the CB and the PB data to correct the well known artefacts observed when reconstructing CB data acquired with a circular orbit. The mixed algorithm is about twice as fast as CB-FBP and results in better image quality, due to decreased discretization errors. However, the two algorithms yield comparable image quality when applied to a disc phantom measured with a two-headed SPECT system.

Algorithms↗

A performance study of 3D reconstruction algorithms for positron emission tomography.

This paper investigates the statistical and systematic accuracy of five three-dimensional reconstruction algorithms for multi-ring PET scanners operated without septa: the reprojection method, the direct Fourier reconstruction, the FAVOR algorithm, and the single-slice and multi-slice rebinning algorithms. Simulated data of a uniform cylinder, of Gaussian sources, and of spherical sources are used to compare respectively the noise properties, the modulation transfer function, and the recovery coefficients of the algorithms. Brain scans reconstructed with the different algorithms are compared by calculating the linear regression of the mean values within regions of interest. The most significant observations are a slight loss of transaxial resolution with the reprojection algorithm in the external slices of the scanner, and increased noise in the images reconstructed using multi-slice rebinning.

Algorithms↗

A rotating PET scanner using BGO block detectors: design, performance and applications.

Recent advances in fully three-dimensional reconstruction for multi-ring PET scanners have led us to explore the potential of a prototype scanner based on the rotation of two opposing arrays of BGO block detectors. The prototype contains only one-third of the number of detectors in the equivalent full ring scanner, resulting in reduced cost. With a lower energy threshold at 250 keV, the absolute efficiency of the scanner is 0.5% and the scatter fraction is 35% for a 20-cm cylinder. Transaxial and axial spatial resolution is about 6 mm. The maximum noise equivalent count rate estimated for a 15-cm diameter cylinder is 36,000 cps at a concentration of 26 kBq/ml. The minimum scan time for a 18F-fluoro-2-deoxyglucose (FDG) brain study is 55 sec. The camera has been validated for clinical applications using both FDG and 82Rb.

Algorithms↗

A normalization technique for 3D PET data.

Prior to reconstruction, emission data from a multi-ring PET camera must be corrected (normalized) for variations in detector sensitivity. The appropriate correction coefficients are obtained by measuring the response of all coincidence lines to a calibrated source of activity (a blank scan). State-of-the-art cameras may contain up to a million such lines of response (LORs), and therefore around 400 million counts will be required to calibrate each LOR to a statistical accuracy of 5%. Alternatively, by modelling the LOR sensitivity as the product of the individual detector efficiencies and a geometrical factor, a calibration procedure has been proposed which requires the determination of only 6000 parameters from this same data set. A significant improvement in the statistical accuracy of the coefficients can therefore be expected. Recently, multi-ring scanners have been operated with the septa retracted, increasing the number of measured LORs by a factor of eight. The acquisition of the calibration data necessary to achieve adequate statistical accuracy then becomes prohibitive. We show that, by modelling the LOR sensitivity, it is possible, with certain approximations, to normalize a septa-retracted emission data set with good accuracy. The input to the model is a high statistics blank scan acquired with the septa extended, which offers a number of practical advantages.

Calibration↗

Normalisation and reconstruction of PET data acquired by a multi-ring camera with septa retracted.

Emission scan data acquired by a multi-ring PET camera operated with septa retracted must be corrected for (1) geometrical and detector sensitivity variations between the different lines of response (normalisation), (2) photon attenuation, and (3) mispositioned events due to photon scattering. These corrections must be applied to the full 3-D set of lines of response before reconstruction. The standard normalisation and attenuation correction procedures for 2-D scans increase the statistical noise in the emission scan, a problem which becomes even more serious in 3-D because of the large number of LORs involved (approximately 8 million). This paper will describe a fully 3-D reconstruction algorithm for multi-angle PET data incorporating a practical normalisation and attenuation correction procedure which minimises the increase in emission scan statistical noise. The correction factors are derived from 2-D, septa extended scans. The algorithm is currently used to reconstruct 3-D emission data from an ECAT 953B, a sixteen-ring PET camera with retractable septa.

Algorithms↗

Implementation of three-dimensional image reconstruction for multi-ring positron tomographs.

In view of the number of PET studies involving low count rate acquisitions, there has been increasing interest recently in the development of positron cameras capable of fully three-dimensional acquisition and reconstruction. This interest has given impetus to the study of algorithms for 3D reconstruction, including those algorithms suitable for application to multi-ring PET scanners. While 2D reconstruction methods can often be generalised to 3D, a number of implementation problems arise which are unique to the 3D approach. This paper examines some of the difficulties associated with the generalisation of the filtered backprojection algorithm to 3D, paying particular attention to the approximations and variable transformations required for application to data from a multi-ring scanner.

Humans↗

Aspects of three dimensional reconstruction for a multi ring positron tomograph.

An important feature of multi ring positron tomographs is the inter plane septa, the purpose of which is to reduce random and scattered coincidences. In general, such septa also eliminate the coincidence lines of response between pairs of detectors more than one ring apart. The operation of a camera without septa must result in an increase not only in the true coincidence rate, but also in the singles, and therefore in the dead time and randoms rate, and in the scattered coincidences. A configuration option in the coincidence hardware of the 8 ring, 15 slice ECAT 931/08-12 enables a full set of 64 sinograms to be acquired when the septa are removed. The detector normalisation and transmission data for studies with the septa out can be obtained using a rotating pin source. To take maximum advantage of the additional signal, the emission data must be reconstructed using a fully three dimensional reconstruction algorithm. This paper presents an analysis of some phantom studies acquired without septa and reconstructed in three dimensions. The results are compared with data acquired with septa for the same phantoms imaged under similar conditions. It is found that, with the septa removed, the signal to noise for a uniform, 20 cm diameter cylinder improves by a factor of 2.8 in the centre of the field of view, whereas in regions distant from the centre in the axial direction, the signal to noise decreases due to the increase in scatter and randoms. An improvement in signal to noise is observed in 6 cm of the 10 cm axial length of the tomograph.

Image Processing, Computer-Assisted↗

Three-dimensional image reconstruction from complete projections.

Three-dimensional medical image reconstruction for both transmission and emission tomography has traditionally decomposed the problem into a set of two-dimensional reconstructions on parallel transverse sections. There is, however, increasing interest in reconstructing projection data directly in three dimensions. For emission tomography in particular, such a reconstruction procedure would clearly make more efficient use of the available photon flux. In the past few years, a number of authors have studied the problems associated with full three-dimensional reconstruction, especially in the case of positron tomography where three-dimensional reconstruction is likely to offer the greatest benefits. While most approaches follow that of filtered backprojection, the relationship between the various filters that have been proposed is far from evident. This paper clarifies this relationship by analysing and generalising the different classes of published filters and establishes the properties and characteristics of a general solution to the three-dimensional reconstruction problem. Some guidelines are suggested for the choice of an appropriate filter in a given situation.

Algorithms↗