Direct reconstruction of single-photon emission computed tomography images using retained matrix elements

In clinical applications, two methods of single-photon emission computed tomography (SPECT) reconstruction are widely used. These are filtered backprojection and iterative reconstruction. Filtered backprojection is fast and produces acceptable images. Iterative reconstruction is slow, but produces images of greater accuracy than backprojection. The authors sought to develop a method of SPECT reconstruction that would have the advantages of both established methods: close in speed to backprojection and with the accuracy of iterative reconstruction. This was accomplished by computing a direct solution to the set of linear equations governing SPECT reconstruction. We tested this method of SPECT reconstruction using a set of projections from a cold rod and sphere phantom. Direct reconstruction produced images having equivalent resolution to backprojected images, but with double the contrast ratio. The direct method required 10 seconds of computation per slice on a Macintosh Quadra 950 (Apple Computer; Cupertin, CA), significantly faster than most iterative methods.     

Iterative reconstruction techniques in emission computed tomography

In emission tomography statistically based iterative methods can improve image quality relative to analytic image reconstruction through more accurate physical and statistical modelling of high-energy photon production and detection processes. Continued exponential improvements in computing power, coupled with the development of fast algorithms, have made routine use of iterative techniques practical, resulting in their increasing popularity in both clinical and research environments. Here we review recent progress in developing statistically based iterative techniques for emission computed tomography. We describe the different formulations of the emission image reconstruction problem and their properties. We then describe the numerical algorithms that are used for optimizing these functions and illustrate their behaviour using small scale simulations.     

Critical examination of the uniformity requirements for single-photon emission computed tomography

It is generally recognized that single-photon emission computed tomography (SPECT) imposes very stringent requirements on gamma camera uniformity to prevent the occurrence of ring artifacts. The purpose of this study was to examine the relationship between nonuniformities in the planar data and the magnitude of the consequential ring artifacts in the transaxial data, and how the perception of these artifacts is influenced by factors such as reconstruction matrix size, reconstruction filter, and image noise. The study indicates that the relationship between ring artifact magnitude and image noise is essentially independent of the acquisition or reconstruction matrix sizes, but is strongly dependent upon the type of smoothing filter applied during the reconstruction process. Furthermore, the degree to which a ring artifact can be perceived above image noise is dependent on the size and location of the nonuniformity in the planar data, with small nonuniformities (1-2 pixels wide) close to the center of rotation being less perceptible than those further out (8-20 pixels). Small defects or nonuniformities close to the center of rotation are thought to cause the greatest potential corruption to tomographic data. The study indicates that such may not be the case. Hence the uniformity requirements for SPECT may be less demanding than was previously thought.     

Iterative concurrent reconstruction algorithms for emission computed tomography

Direct reconstruction techniques, such as those based on filtered backprojection, are typically used for emission computed tomography (ECT), even though it has been argued that iterative reconstruction methods may produce better clinical images. The major disadvantage of iterative reconstruction algorithms, and a significant reason for their lack of clinical acceptance, is their computational burden. We outline a new class of 'concurrent' iterative reconstruction techniques for ECT in which the reconstruction process is reorganized such that a significant fraction of the computational processing occurs concurrently with the acquisition of ECT projection data. These new algorithms use the 10-30 min required for acquisition of a typical SPECT scan to iteratively process the available projection data, significantly reducing the requirements for post-acquisition processing. These algorithms are tested on SPECT projection data from a Hoffman brain phantom acquired with a 2 x 10(5) counts in 64 views each having 64 projections. The SPECT images are reconstructed as 64 x 64 tomograms, starting with six angular views. Other angular views are added to the reconstruction process sequentially, in a manner that reflects their availability for a typical acquisition protocol. The results suggest that if T s of concurrent processing are used, the reconstruction processing time required after completion of the data acquisition can be reduced by at least 1/3T s.     

Whole-body single-photon emission computed tomography using dual, large-field-of-view scintillation cameras.

A whole-body single-photon emission computed tomography system (SPECT) consisting of two large-field-of-view scintillation cameras mounted on a rotatable gantry, a minicomputer and a display station has been designed, constructed and evaluated. In its usual mode of operation, eleven contiguous transverse sections, each 12.5 or 25 mm thick, are reconstructed from projection data acquired during a single, continuous 360 degree rotation lasting from 2 to 22 min. A generalised filtered and weighted backprojection algorithm is used to reconstruct data obtained with conventional parallel-hole collimators in the case of body scanning, or with specially designed fan beam collimators in the case of centrally positioned organs. A simple, yet effective, correction is used to compensate for the effects of gamma ray attenuation within the patient. In addition to providing transverse section images, the system is capable of simultaneous acquisition of opposed conventional scintigrams, the reconstruction of longitudinal section images, and the acquisition of gated cardiac transverse sections. Resolutions in the reconstructed images are typically 15 mm for body scans and 11 mm for brain scans, with only slight variations in sensitivity and resolution within the image. Phantoms and clinical data demonstrate that the SPECT system generates high quality section images while maintaining most of the flexibility of normal scintillation cameras, with the added advantage of dual heads.     

Simultaneous iterative reconstruction for emission and attenuation images in positron emission tomography

The quality of the attenuation correction strongly influences the outcome of the reconstructed emission scan in positron emission tomography. Usually the attenuation correction factors are calculated from the transmission and blank scan and thereafter applied during the reconstruction on the emission data. However, this is not an optimal treatment of the available data, because the emission data themselves contain additional information about attenuation: The optimal treatment must use this information for the determination of the attenuation correction factors. Therefore, our purpose is to investigate a simultaneous emission and attenuation image reconstruction using a maximum likelihood estimator, which takes the attenuation information in the emission data into account. The total maximum likelihood function for emission and transmission is used to derive a one-dimensional Newton-like algorithm for the calculation of the emission and attenuation image. Log-likelihood convergence, mean differences, and the mean of squared differences for the emission image and the attenuation correction factors of a mathematical thorax phantom were determined and compared. As a result we obtain images improved with respect to log likelihood in all cases and with respect to our figures of merit in most cases. We conclude that the simultaneous reconstruction can improve the performance of image reconstruction.     

Convolutional image reconstruction for quantitative single photon emission computed tomography

Two image reconstruction algorithms have been investigated. They are based on filtered backprojection, and are useful when the tissue attenuation is considered to be uniform in the object. The first method uses a weighted backprojection, the weighting factor being determined in such a way that the photon attenuation is compensated with low noise propagation. The parameters involved in the convolution kernel and the correction function were determined by a computer iteration program. The second method, which is a simplified version of the first, uses conventional backprojection, and takes a shorter computation time than the first method. The statistical noise of an image can be minimised by suitable positioning of the coordinate origin for the reconstruction. The theory of the two methods, their performance on statistical noise and some results of mathematical and experimental phantom studies are described.     

Accelerated median root prior reconstruction for pinhole single-photon emission tomography (SPET)

Pinhole collimation can be used to improve spatial resolution in SPET. However, the resolution improvement is achieved at the cost of reduced sensitivity, which leads to projection images with poor statistics. Images reconstructed from these projections using the maximum likelihood expectation maximization (ML-EM) algorithms, which have been used to reduce the artefacts generated by the filtered backprojection (FBP) based reconstruction, suffer from noise/bias trade-off: noise contaminates the images at high iteration numbers, whereas early abortion of the algorithm produces images that are excessively smooth and biased towards the initial estimate of the algorithm. To limit the noise accumulation we propose the use of the pinhole median root prior (PH-MRP) reconstruction algorithm. MRP is a Bayesian reconstruction method that has already been used in PET imaging and shown to possess good noise reduction and edge preservation properties. In this study the PH-MRP algorithm was accelerated with the ordered subsets (OS) procedure and compared to the FBP, OS-EM and conventional Bayesian reconstruction methods in terms of noise reduction, quantitative accuracy, edge preservation and visual quality. The results showed that the accelerated PH-MRP algorithm was very robust. It provided visually pleasing images with lower noise level than the FBP or OS-EM and with smaller bias and sharper edges than the conventional Bayesian methods.     

The peanut orbit: a modified elliptical orbit for single-photon emission computed tomography imaging

The front collimator surface of a conventional single-photon emission computed tomography (SPECT) detector system transcribing an elliptical orbit to approximate body contour could, under certain circumstances, penetrate the ellipse and make patient contact. The problem is associated both with the large front surface dimensions of tomographic cameras and the need to maintain the camera perpendicular to a radius through the axis of rotation. To aid in the development of an improved body contour orbit, software has been developed to simulate SPECT imaging systems. A major feature of the algorithms is a spatially calibrated graphic representation of the pallet, patient and orbit of the camera head. Based on computer simulations performed with this software, a modified elliptical orbit has been proposed for patient contour SPECT scanning on two different types of tomographic acquisition systems: a dual camera, rotating ring system and a single camera, multimotion stand system. This peanut orbit has the advantage of minimizing collimator patient distance, as does the elliptical orbit, but also compensates for camera motion that could potentially produce patient contact. Versions of the peanut orbit algorithms have been installed and are operational on clinical systems similar to those that were simulated. The ring-based gantry system (with cameras mounted on cantilevered arms) requires additional software to correct for translational shift in the field of view as a function of radial distance from the center of the system. This is done during uniformity correction. Standard unmodified backprojection software is used for reconstruction. The peanut orbit may represent an improved approach to body contour imaging in rotational SPECT.     

Convergence of the maximum likelihood reconstruction algorithm for emission computed tomography

Convergence properties of the maximum likelihood estimator (MLE) for emission computed tomographic (ECT) image reconstruction are evaluated as a function of Poisson noise, precision of the assumed system resolution model and iteration number up to 10,000 iterations. In the ECT reconstruction problem, the photon-emitting source distribution is to be estimated from measurements of projections of the emitted photon flux. The MLE algorithm seeks a source distribution which will maximise the maximum likelihood function relating the estimated and the measured projections. A Monte Carlo model of the system transfer function of a single photon emission computed tomographic (SPECT) system allowed realistic projection data to be simulated from a known source distribution. Poisson noise was added to the Monte Carlo simulations. By using projection data from a known source distribution generated through a known system transfer function, we were able to simultaneously evaluate the convergence of both the projection estimations as well as the source distribution estimations. As predicted by theory, the estimates of the projections did continue to improve (or remain the same) for all combinations of Poisson noise (up to 10% RMS) and system resolution (+/- 10% of true value) tested. Convergence of source distribution estimates to the true value was found for up to 10,000 iterations only for low noise (0.1% RMS) with the correct resolution function. For all other combinations, there was some optimum iteration (between 30 and 400) after which the source estimate was degraded even though the estimate of the projections was improved.     

Simultaneous iterative reconstruction of emission and attenuation images in positron emission tomography from emission data only

For quantitative image reconstruction in positron emission tomography attenuation correction is mandatory. In case that no data are available for the calculation of the attenuation correction factors one can try to determine them from the emission data alone. However, it is not clear if the information content is sufficient to yield an adequate attenuation correction together with a satisfactory activity distribution. Therefore, we determined the log likelihood distribution for a thorax phantom depending on the choice of attenuation and activity pixel values to measure the crosstalk between both. In addition an iterative image reconstruction (one-dimensional Newton-type algorithm with a maximum likelihood estimator), which simultaneously reconstructs the images of the activity distribution and the attenuation coefficients is used to demonstrate the problems and possibilities of such a reconstruction. As result we show that for a change of the log likelihood in the range of statistical noise, the associated change in the activity value of a structure is between 6% and 263%. In addition, we show that it is not possible to choose the best maximum on the basis of the log likelihood when a regularization is used, because the coupling between different structures mediated by the (smoothing) regularization prevents an adequate solution due to crosstalk. We conclude that taking into account the attenuation information in the emission data improves the performance of image reconstruction with respect to the bias of the activities, however, the reconstruction still is not quantitative.     

Increased single-photon emission computed tomography image processing speed achieved in personal computers with memory-intensive algorithms

Recent dramatic reductions in the cost of computer random access memory (RAM) and the ability of newer microprocessors and associated personal computer operating systems to address large amounts of memory make novel strategies for high-speed image processing possible. We developed image processing algorithms that use this newly available memory to achieve increases in effective processing speed. These algorithms rely on the use of precomputed lookup tables to avoid repeated use of relatively expensive machine instructions, such as multiplications and divisions. Programs using this strategy to perform single photon emission computer tomography (SPECT) analysis were written in C and assembly language and tested on a Macintosh Quadra 950 (Apple Computer, Cupertino, CA) having 64 megabytes of RAM. The measured processing times are competitive with most dedicated nuclear medicine computers. A general implementation of such programs will allow personal computers to compete with dedicated imaging systems, at a substantial reduction in cost.     

An approximate reconstruction method for helical cone-beam differential phase-contrast computed tomography images

Differential phase-contrast computed tomography (DPC-CT) is a novel x-ray inspection method. Currently, DPC-CT reconstruction problems are solved by using parallel-beam, fan-beam and cone-beam algorithms. The above algorithms cannot show the internal structures of rod-shaped objects well enough for only reconstructing a few slices. The helical cone-beam algorithms have significant advantages for rod-shaped objects over other algorithms. Along with our numerical evaluation and verification, we report a PI-line-based approximate algorithm for helical cone-beam DPC-CT, which can be applied to reconstruct the refractive index decrement distribution of the samples directly from phase-contrast projection images. Simulations using a 3D Shepp-Logan phantom are carried out to verify the proposed algorithm. Reconstruction results show that the proposed algorithm can provide higher quality performance compared with the existing interpolation-based reconstruction algorithm.     

Use of fuzzy edge single-photon emission computed tomography analysis in definite Alzheimer's disease - a retrospective study

Background  

Definite Alzheimer's disease (AD) requires neuropathological confirmation. Single-photon emission computed tomography (SPECT) may enhance diagnostic accuracy, but due to restricted sensitivity and specificity, the role of SPECT is largely limited with regard to this purpose.     

Modifying constrained least-squares restoration for application to single photon emission computed tomography projection images

Image restoration methods have been shown to increase the contrast of nuclear medicine images by decreasing the effects of scatter and septal penetration. Image restoration can also reduce the high-frequency noise in the image. This study applies constrained least-squares (CLS) restoration to the projection images of single photon emission computed tomography (SPECT). In a previous study, it was noted that CLS restoration has the potential advantage of automatically adapting to the blurred object. This potential is confirmed using planar images. CLS restoration is then modified to improve its performance when applied to SPECT projection image sets. The modification was necessary because the Poisson noise in low count SPECT images causes considerable variation in the CLS filter. On phantom studies, count-dependent Metz restoration was slightly better than the modified CLS restoration method, according to measures of contrast and noise. However, CLS restoration was generally judged as yielding the best results when applied to clinical studies, apparently because of its ability to adapt to the image being restored.     

正电子断层扫描-PET的图像重建

随着计算机技术的迅速发展，ＰＥＴ图像重建方法得以很大的改进，其研究工作越来越受到人们的重视。特别是对能压制噪声和提高空间分辨率的图像快速重建算法的研究是目前一个热门课题。本文简介了ＰＥＴ图像重建的基本原理与基本方法及其研究与应用现状。此外，对ＰＥＴ的三维重建也做了一些简介。     

99Tcm-MIBI SPECT/CT显像联合超声检查在甲状旁腺功能亢进症合并甲状腺癌诊断中的应用

目的 探讨⁹⁹Tc^m-甲氧基异丁基异腈(MIBI)单光子发射计算机断层成像/计算机断层扫描(SPECT/CT)显像联合超声对甲状旁腺功能亢进症(HPT)合并甲状腺癌(TC)的术前诊断价值。方法 回顾性研究。纳入2013年9月—2019年12月经东南大学医学院附属江阴医院乳甲科手术及病理证实的HPT合并TC患者15例,其中男3例、女12例,年龄35～57岁(中位年龄49岁);继发性HPT 10例,原发性HPT 5例。统计术中检出并经术后病理证实的甲状旁腺病灶及TC病灶数;由多学科团队分析术前⁹⁹Tc^m-MIBI SPECT/CT显像及超声检查,比较⁹⁹Tc^m-MIBI SPECT/CT显像及超声对病变甲状旁腺及TC的检出率有无差异。结果 15例HPT合并TC患者术中检出并经术后病理证实的甲状旁腺病灶44枚,⁹⁹Tc^m-MIBI SPECT/CT显像的检出率为86.4%(38/44),高于超声的检出率65.9%(29/44),差异有统计学意义(P＜0.05)。10例继发性HPT合并TC患者手术检出39枚甲状旁腺病灶,⁹⁹Tc^m-MIBI SPECT/CT显像及超声则分别检出33枚(84.6%)和27枚(69.2%),两者比较差异无统计学意义(P＞0.05);5例原发性HPT合并TC患者均为单发病灶,⁹⁹Tc^m-MIBI SPECT/CT显像检出5枚,超声检出2枚、漏检3枚。15例患者术中检出并经术后病理证实的TC病灶共16枚,其中单发14例、多发1例(2枚),⁹⁹Tc^m-MIBI SPECT/CT显像对TC病灶的检出率为5/16低于超声的12/16,两者比较差异有统计学意义(P＜0.05)。结论 ⁹⁹Tc^m-MIBI SPECT/CT显像对HPT合并TC的甲状旁腺病灶的检出率优于超声,而超声对合并的TC诊断优于⁹⁹Tc^m-MIBI SPECT/CT显像,两者联合应用有助于HPT合并TC的诊断。