Search PubMed⌕ Search

Biomedical subjects

B F Hutton

Publications and source records attributed to B F Hutton.

At least 19 recordsLinked to original sources
PubMed2004

Improved tolerance to missing data in myocardial perfusion SPET using OSEM reconstruction.

When projection data are incomplete for various technical reasons, artefacts may occur in the reconstructed images. This study examines whether an iterative reconstruction method, the ordered subsets implementation of the EM algorithm (OSEM), can improve reconstruction and minimise the artefacts compared to filtered back-projection (FBP). We varied the number and location of projections removed to investigate when significant artefacts occur, and whether diagnosis is affected. Phantom studies were analysed with sequential orthogonal pairs of projection angles removed (as would typically occur when either data loss or severe motion is detected during acquisition with a right-angled, dual-head cardiac single-photon emission tomography system) and reconstructed with both FBP and OSEM. Twelve normal myocardial perfusion studies were also assessed to study the effect of missing projections on clinical diagnosis. Differences between reconstructions with intact versus missing data were measured. Also, reconstructed images were clinically assessed and scored on a five-point scale based on whether the artefacts would alter clinical interpretation. Although both reconstruction methods showed artefacts, the absolute differences between reconstructed phantom data with intact and missing projection sets were significantly greater (P<0.005) for FBP than for OSEM for all numbers of missing projections. The clinical data showed similar differences between FBP and OSEM reconstructions. The three observers noted superiority of OSEM compared to FBP, with reduced incidence of clinically significant artefacts. However, neither reconstruction method could tolerate six or more missing pairs from 32 projections. There was no significant dependence on the angular location of missing projections. In the absence of any attempt to correct for missing projections, OSEM reduced the influence of artefacts compared to FBP.

Algorithms↗
PubMed2001

Choice of collimator for cardiac SPET when resolution compensation is included in iterative reconstruction.

In clinical cardiac single-photon emission tomography (SPET) studies, collimators of different spatial resolution and geometric efficiency are available for imaging. In selecting the appropriate collimator for clinical use, there is a trade-off between spatial resolution, which can limit the contrast of the reconstructed image, and detection efficiency, which determines the noise in the image. Our objective was to assess which collimator is best suited for cardiac SPET when reconstruction is performed with and without compensation for distance-dependent resolution (CDR). The dynamic MCAT thorax phantom was used to simulate 180 degree technetium-99m cardiac data, acquired using either a general-purpose (GP) or high-resolution (HR) collimator. For GP and HR, the resolution at 15 cm was 11.5 mm and 9.5 mm respectively, and the corresponding relative efficiency was 1.0 and 0.52 respectively. Distance-dependent resolution, attenuation and noise were included in the projection data; scatter was not included. Ordered subsets expectation maximisation reconstruction (subset size 4) was performed with and without CDR. Results were evaluated by comparing the myocardial recovery coefficient and contrast between myocardium and ventricle relative to the original phantom, each plotted for different noise levels corresponding to increasing iteration number. The study demonstrated that, without CDR, HR gave the best results. However, for any given noise level with CDR, GP gave superior recovery and contrast. These findings were confirmed in a physical phantom study. Results suggest that improved reconstruction can be achieved using a GP collimator in combination with resolution compensation.

Algorithms↗
PubMed2001

Non-rigid image registration using a median-filtered coarse-to-fine displacement field and a symmetric correlation ratio.

Conventional approaches to image registration are generally limited to image-wide rigid transformations. However, the body and its internal organs are non-rigid structures that change shape due to changes in the body's posture during image acquisition, and due to normal, pathological and treatment-related variations. Inter-subject matching also constitutes a non-rigid registration problem. In this paper, we present a fully automated non-rigid image registration method that maximizes a local voxel-based similarity metric. Overlapping image blocks are defined on a 3D grid. The transformation vector field representing image deformation is found by translating each block so as to maximize the local similarity measure. The resulting sparsely sampled vector field is median filtered and interpolated by a Gaussian function to ensure a locally smooth transformation. A hierarchical strategy is adopted to progressively establish local registration associated with image structures at diminishing scale. Simulation studies were carried out to evaluate the proposed algorithm and to determine the robustness of various voxel-based cost functions. Mutual information, normalized mutual information, correlation ratio (CR) and a new symmetric version of CR were evaluated and compared. A T1-weighted magnetic resonance (MR) image was used to test intra-modality registration. Proton density and T2-weighted MR images of the same subject were used to evaluate inter-modality registration. The proposed algorithm was tested on the 2D MR images distorted by known deformations and 3D images simulating inter-subject distortions. We studied the robustness of cost functions with respect to image sampling. Results indicate that the symmetric CR gives comparable registration to mutual information in intra- and inter-modality tasks at full sampling and is superior to mutual information in registering sparsely sampled images.

Algorithms↗
PubMed2000

Improved efficiency for MRI-SPET registration based on mutual information.

Mutual information has been proposed as a criterion for image registration. The criterion is calculated from a two-dimensional grey-scale histogram of the image pair being registered. In this paper we study how sparse sampling can be used to increase speed performance using the registration algorithm of Maes et al. (IEEE Trans Med Imaging 1997; 16: 187-198) with a focus on registration of MRI-SPET brain images. In particular we investigate how sparse sampling and parameters such as the number of bins used for the grey-scale histograms and smoothing of the data prior to registration affect accuracy and robustness of the registration. The method was validated using both simulated and human data. Our results show that sparse sampling introduced local maxima into the mutual information similarity function when the number of bins used for the histograms was large. To speed up registration while retaining robustness, smoothing of the data prior to registration was used and a coarse to fine subsampling protocol, where the number of bins in the histograms were dependent on the subsampling factor, was employed, For the simulated data, the method was able to recover known transformations with an accuracy of about 1 mm. Using the human data, there were no significant differences in the recovered transformation parameters when the suggested subsampling scheme was used compared with when no subsampling was used, but there was a more than tenfold increase in speed. Our results show that, with the appropriate choice of parameters, the method can accurately register MRI-SPET brain images even when very efficient sampling protocols are used.

Algorithms↗
PubMed1999

A distance-assisted training program for nuclear medicine technologists.

We have developed training materials for nuclear medicine technologists to be used in distance-assisted training programs. We have completed our first pilot project in Asia and report that there will be nearly 500 students around the world, in Asia, Africa, Central America and South America, using our materials during the coming year.

Allied Health Personnel↗
PubMed1998

Efficient scatter modelling for incorporation in maximum likelihood reconstruction.

Definition of a simplified model of scatter which can be incorporated in maximum likelihood reconstruction for single-photon emission tomography (SPET) continues to be appealing; however, implementation must be efficient for it to be clinically applicable. In this paper an efficient algorithm for scatter estimation is described in which the spatial scatter distribution is implemented as a spatially invariant convolution for points of constant depth in tissue. The scatter estimate is weighted by a space-dependent build-up factor based on the measured attenuation in tissue. Monte Carlo simulation of a realistic thorax phantom was used to validate this approach. Further efficiency was introduced by estimating scatter once after a small number of iterations using the ordered subsets expectation maximisation (OSEM) reconstruction algorithm. The scatter estimate was incorporated as a constant term in subsequent iterations rather than modifying the scatter estimate each iteration. Monte Carlo simulation was used to demonstrate that the scatter estimate does not change significantly provided at least two iterations OSEM reconstruction, subset size 8, is used. Complete scatter-corrected reconstruction of 64 projections of 40¿128 pixels was achieved in 38 min using a Sun Sparc20 computer.

Algorithms↗
PubMed1998

Application of distance-dependent resolution compensation and post-reconstruction filtering for myocardial SPECT.

Compensation for distance-dependent resolution can be directly incorporated in maximum likelihood reconstruction. Our objective was to examine the effectiveness of this compensation using either the standard expectation maximization (EM) algorithm or an accelerated algorithm based on use of ordered subsets (OSEM). We also investigated the application of post-reconstruction filtering in combination with resolution compensation. Using the MCAT phantom, projections were simulated for 360 degrees data, including attenuation and distance-dependent resolution. Projection data were reconstructed using conventional EM and OSEM with subset size 2 and 4, with/without 3D compensation for detector response (CDR). Also post-reconstruction filtering (PRF) was performed using a 3D Butterworth filter of order 5 with various cutoff frequencies (0.2-1.2 cycles cm(-1)). Image quality and reconstruction accuracy were improved when CDR was included. Image noise was lower with CDR for a given iteration number. PRF with cutoff frequency greater than 0.6 cycles cm(-1) improved noise with no reduction in recovery coefficient for myocardium but the effect was less when CDR was incorporated in the reconstruction. CDR alone provided better results than use of PRF without CDR. Results suggest that using CDR without PRF, and stopping at a small number of iterations, may provide sufficiently good results for myocardial SPECT. Similar behaviour was demonstrated for OSEM.

Algorithms↗
PubMed1998

Dynamic imaging and tracer kinetic modeling for emission tomography using rotating detectors.

When performing dynamic studies using emission tomography the tracer distribution changes during acquisition of a single set of projections. This is particularly true for some positron emission tomography (PET) systems which, like single photon emission computed tomography (SPECT), acquire data over a limited angle at any time, with full projections obtained by rotation of the detectors. In this paper, an approach is proposed for processing data from these systems, applicable to either PET or SPECT. A method of interpolation, based on overlapped parabolas, is used to obtain an estimate of the total counts in each pixel of the projections for each required frame-interval, which is the total time to acquire a single complete set of projections necessary for reconstruction. The resultant projections are reconstructed using traditional filtered backprojection (FBP) and tracer kinetic parameters are estimated using a method which relies on counts integrated over the frame-interval rather than instantaneous values. Simulated data were used to illustrate the technique's capabilities with noise levels typical of those encountered in either PET or SPECT. Dynamic datasets were constructed, based on kinetic parameters for fluoro-deoxy-glucose (FDG) and use of either a full ring detector or rotating detector acquisition. For the rotating detector, use of the interpolation scheme provided reconstructed dynamic images with reduced artefacts compared to unprocessed data or use of linear interpolation. Estimates for the metabolic rate of glucose had similar bias to those obtained from a full ring detector.

Algorithms↗
PubMed1998

Half-fanbeam collimators combined with scanning point sources for simultaneous emission-transmission imaging.

UNLABELLED: One type of SPECT system often used for simultaneous emission-transmission tomography is equipped with parallel-hole collimators, moving line sources (MLS) and electronic windows that move in synchrony with the sources. Although downscatter from the emission distribution is reduced by the use of the electronic window, this still can represent a sizable fraction of the transmitted counts. These systems have relatively poor spatial resolution and use costly transmission sources. METHODS: Using a two-head SPECT system, with heads at right angles, two 153Gd line sources (5800 MBq each) were replaced by two 153Gd point sources of only 750 MBq each and positioned to move along the focal lines of two half-fanbeam collimators. A suitable acquisition protocol for a moving point source (MPS) system was selected by considering the results of a simulation study. With this protocol, physical phantom experiments were conducted. RESULTS: Simulations showed that by using two half-fanbeam collimators, a gantry rotation of 90 degrees, such as used for 180 degrees acquisition with parallel-beam collimators for cardiac imaging, was insufficient. A gantry rotation of 180 degrees resulted in attenuation maps where only an area to the posterior of a 400-mm wide thorax phantom was affected by truncation. The MPS system had a 14.7 times higher sensitivity for transmission counts than the MLS system. Despite the smaller sources in the MPS system, the number of acquired transmission counts was a factor 1.91 times higher compared with the MLS system, resulting in reduced noise. The relative downscatter contribution from 99mTc (140 keV) in the 153Gd moving electronic window (100 keV) was reduced by a factor of 1.81. Transmission images of a rod phantom with segments containing acrylic rods of different diameters showed an improvement of resolution in favor of the MPS system from about 11 mm to about 6 mm (five instead of two segments of rods were clearly visible). In addition, the noise level in the MPS thorax transmission images was significantly lower. CONCLUSION: The MPS system has important advantages when compared with the MLS system. The use of low-activity point sources is economically beneficial when compared with line sources and reduces radiation exposure to staff and patients.

Equipment Design↗
PubMed1997

A clinical perspective of accelerated statistical reconstruction.

Although the potential benefits of maximum likelihood reconstruction have been recognised for many years, the technique has only recently found widespread popularity in clinical practice. Factors which have contributed to the wider acceptance include improved models for the emission process, better understanding of the properties of the algorithm and, not least, the practicality of application with the development of acceleration schemes and the improved speed of computers. The objective in this article is to present a framework for applying maximum likelihood reconstruction for a wide range of clinically based problems. The article draws particularly on the experience of the three authors in applying an acceleration scheme involving use of ordered subsets to a range of applications. The potential advantages of statistical reconstruction techniques include: (a) the ability to better model the emission and detection process, in order to make the reconstruction converge to a quantitative image, (b) the inclusion of a statistical noise model which results in better noise characteristics, and (c) the possibility to incorporate prior knowledge about the distribution being imaged. The great flexibility in adapting the reconstruction for a specific model results in these techniques having wide applicability to problems in clinical nuclear medicine.

Algorithms↗
PubMed1996

Correlation between post-ejection shortening and improvement in regional wall motion after revascularization in patients with coronary artery disease.

This study investigates the relationship between regional post-ejection shortening and improvement in wall motion after revascularization in patients with severe regional left ventricular contractile dysfunction. Canine studies have documented the phenomenon of post-systolic shortening during acute ischemia, and this predicted recovery of contractile function when ischemia was relieved. This delayed shortening from the ischemic segment would cause regional emptying while the other normal portion of the left ventricle starts to have diastolic relaxation. This was detectable by standard phase analysis of the radionuclide ventriculogram, the abnormal region with delayed emptying being distinguished from the normal portion of left ventricle as an area with a homogeneous phase delay. Twelve patients with successful revascularization and a pre-operative study demonstrating a single homogeneous area of phase delay were identified. The area was demarcated and the regional volume-time curve constructed. The volume-time curve of the abnormal region is analogous to the myocardial segment length versus time relationship in the canine model. We quantitated the amount of delayed (post-systolic) emptying in the demarcated region as the difference between end-systole counts and post-systolic nadir counts, and this was normalised to left ventricular stroke count. After revascularization, regional ejection fraction improved from 44 +/- 10% to 62 +/- 14% (P < 0.001), representing a 47 +/- 50% improvement over baseline. The percentage improvement in regional ejection fraction correlated with post-systolic emptying (r = 0.74, P < 0.05) but not with initial regional ejection fraction. In conclusion, post-ejection shortening causes regional post-systolic emptying and this correlates with post-revascularization improvement in regional wall motion.

Adult↗
PubMed1996

Transmission-based scatter correction of 180 degrees myocardial single-photon emission tomographic studies.

Meaningful comparison of single-photon emission tomographic (SPET) reconstructions for data acquired over 180 degrees or 360 degrees can only be performed if both attenuation and scatter correction are applied. Convolution subtraction has appeal as a practical method for scatter correction; however, it is limited to data acquired over 360 degrees. A new algorithm is proposed which can be applied equally well to data acquired over 180 degrees or 360 degrees. The method involves estimating scatter based on knowledge of reconstructed transmission data in combination with a reconstructed estimate of the activity distribution, obtained using attenuation correction with broad beam attenuation coefficients. Processing is implemented for planes of activity parallel to the projection images for which a simplified model for the scatter distribution may be applied, based on the measured attenuation. The appropriate broad beam (effective) attenuation coefficients were determined by considering the scatter buildup equation. It was demonstrated that narrow beam attenuation coefficients should be scaled by 0.75 and 0.65 to provide broad beam attenuation coefficients for technetium-99m and thallium-201 respectively. Using a thorax phantom, quantitative accuracy of the new algorithm was compared with conventional transmission-based convolution subtraction (TDCS) for 360 degrees data. Similar heart to lung contrasts were achieved and correction of 180 degrees data yielded a 10.4% error for cardiac activity compared to 5.2% for TDCS. Contrast for myocardium to ventricular cavity was similarly good for scatter-corrected 180 degrees and 360 degrees data, in contrast to attenuation-corrected data, where contrast was significantly reduced. The new algorithm provides a practical method for correction of scatter applicable to 180 degrees myocardial SPET.

Algorithms↗
PubMed1996

Monte Carlo and experimental evaluation of accuracy and noise properties of two scatter correction methods for SPECT.

Scatter correction is a prerequisite for quantitative SPECT, but potentially increases noise. Monte Carlo simulations (EGS4) and physical phantom measurements were used to compare accuracy and noise properties of two scatter correction techniques: the triple-energy window (TEW), and the transmission dependent convolution subtraction (TDCS) techniques. Two scatter functions were investigated for TDCS: (i) the originally proposed mono-exponential function (TDCSmono) and (ii) an exponential plus Gaussian scatter function (TDCSGauss) demonstrated to be superior from our Monte Carlo simulations. Signal to noise ratio (S/N) and accuracy were investigated in cylindrical phantoms and a chest phantom. Results from each method were compared to the true primary counts (simulations), or known activity concentrations (phantom studies). 99mTc was used in all cases. The optimized TDCS(Gauss) method overall performed best, with an accuracy of better than 4% for all simulations and physical phantom studies. Maximum errors for TEW and TDCS(mono) of -30 and -22%, respectively, were observed in the heart chamber of the simulated chest phantom. TEW had the worst S/N ratio of the three techniques. The S/N ratios of the two TDCS methods were similar and only slightly lower than those of simulated true primary data. Thus, accurate quantitation can be obtained with TDCS(Gauss), with a relatively small reduction in S/N ratio.

Air↗
PubMed1996

Minimum cross-entropy reconstruction of PET images using prior anatomical information.

An algorithm is presented for the reconstruction of PET images using prior anatomical information derived from MR images of the same subject. The cross-entropy or Kullback-Leiber distance is a measure of dissimilarity between two images. We propose to reconstruct PET images by minimizing a weighted sum of two cross-entropy terms. The first is the cross-entropy between the measured emission data and the forward projection of the current estimate of the PET image. Minimizing this term alone is equivalent to the ML-EM reconstruction. The second term is the cross-entropy between the current estimate of the PET image and a prior image model which incorporates anatomical information derived from registered MR images. A weighting parameter determines the relative emphasis given to the emission data and the prior model in the reconstruction. Details of this algorithm are presented as well as test reconstructions for real and simulated data. The performance of the algorithm was evaluated with respect to errors in prior anatomical information. The algorithm provided significant improvement in the quality of reconstructed images as compared with the ML-EM reconstruction technique. The reconstructed images had higher resolution as compared with the images obtained from MAP-like reconstructions which do not utilize anatomical information. The algorithm displayed robustness with respect to errors in prior anatomical information.

Algorithms↗
PubMed1996

Automated interstudy image registration technique for SPECT and PET.

UNLABELLED: We report the extended application of an automated computer technique for three-dimensional spatial registration of SPECT and PET studies. METHODS: The technique iteratively reslices a misaligned data set until the sum of the absolute differences (SAD) from a reference data set is minimized. The registration accuracy was assessed in Hoffman brain phantom studies collected with known misalignments and transmission studies of a thorax phantom with fiducial markers. The SAD was compared with three other cost functions: stochastic sign change criterion, sum of products and standard deviation (s.d.) of ratios. In clinical neurological and myocardial perfusion studies, registration accuracy was estimated from the relative locations of landmarks in the reference and registered data sets. RESULTS: Registration accuracy in the Hoffman brain phantom studies was -0.07 +/- 0.46 mm (mean +/- s.d.) for translations and -0.01 +/- 0.20 degrees for rotations, with maximum translation and rotation errors of 1.2 mm and 0.8 degree, respectively. The SAD was the most accurate and reliable cost function. Registration errors in the thorax phantom were 3.1 +/- 1.7 mm. Mean accuracy in the neurological studies, estimated from landmark pairs, was 2.0 +/- 1.1 mm for SPECT to SPECT and 1.8 +/- 1.1 mm for PET to SPECT registrations. Average registration accuracy in 201Tl myocardial perfusion studies was 2.1 +/- 1.2 mm. CONCLUSION: Our registration method (a) provided accurate registrations for phantom and clinical SPECT and PET studies, (b) is fully automated, (c) simplifies comparison of data sets obtained at different times and with different modalities, and (d) can be applied retrospectively.

Algorithms↗
PubMed1995

Simultaneous emission and transmission measurements for attenuation correction in whole-body PET.

UNLABELLED: We describe a methodology for measuring and correcting for attenuation in whole-body PET using simultaneous emission and transmission (SET) measurements. METHODS: The main components of the methodology are: (a) sinogram windowing of low activity (< or = 50 MBq) rotating 68Ge/Ga rod sources, (b) segmented attenuation correction (SAC) and (c) maximum likelihood reconstruction using the ordered subsets EM (OS-EM) algorithm. The methods were implemented on a whole-body positron emission tomograph. Quantitative accuracy and the signal-to-noise ratio (SNR) were measured for a thorax-tumor phantom as functions of acquisition time (range: 2-20 min per position). RESULTS: When a typical rod source activity (200 MBq 68Ge/Ga) was used, emission SNR was 60% lower in simultaneous than in separate measurements. The difference was only 14% when the rods contained 45 MBq 68Ge/Ga. The SNR was further improved by SAC in conjunction with OS-EM reconstruction and the relative gain increased with increasing acquisition time. Quantitative estimates of tumor, liver and lung radioactivity agreed with values obtained from a separate high count measurement to within 8%, independent of acquisition time. CONCLUSION: Attenuation correction of whole-body PET images is feasible using SET measurements. There is good quantitative agreement with conventional methods and increased noise is offset by the use of SAC and OS-EM reconstruction.

Adult↗
PubMed1994

Accelerated EM reconstruction in total-body PET: potential for improving tumour detectability.

Total-body positron emission tomography (PET) is a useful diagnostic tool for evaluating malignant disease. However, tumour detection is limited by image artefacts due to the lack of attenuation correction and noise. Attenuation correction may be possible using transmission data acquired after or simultaneously with emission data. Despite the elimination of attenuation artefacts, however, tumour detection is still hampered by noise, which is amplified during image reconstruction by filtered backprojection (FBP). We have investigated, as an alternative to FBP, an accelerated expectation maximization (EM) algorithm for its potential to improve tumour detectability in total-body PET. Signal to noise ratio (SNR), calculated for a tumour with respect to the surrounding background, is used as a figure of merit. A software tumour phantom, with conditions typical of those encountered in a total-body PET study using simultaneous acquisition, is used to optimize and compare various reconstruction approaches. Accelerated EM reconstruction followed by two-dimensional filtering is shown to yield significantly higher SNR than FBP for a range of tumour sizes, concentrations and counting statistics (deltaSNR = 6.3 +/- 3.9, p < 0.001). The methods developed are illustrated by examples derived from physical phantom and patient data.

Algorithms↗