[en] Iterative reconstruction from single photon emission computed tomography (SPECT) data requires regularization to avoid noise amplification and edge artefacts in the reconstructed image. This is often accomplished by stopping the iteration process at a relatively low number of iterations or by post-filtering the reconstructed image. The aim of this paper is to develop a method to automatically select an optimal combination of stopping iteration number and filters for a particular imaging situation. To this end different error measures between the distribution of a phantom and a corresponding filtered SPECT image are minimized for different iteration numbers. As a study example, simulated data representing a brain study are used. For post-reconstruction filtering, the performance of 3D linear diffusion (Gaussian filtering) and edge preserving 3D nonlinear diffusion (Catte scheme) is investigated. For reconstruction methods which model the image formation process accurately, error measures between the phantom and the filtered reconstruction are significantly reduced by performing a high number of iterations followed by optimal filtering compared with stopping the iterative process early. Furthermore, this error reduction can be obtained over a wide range of iteration numbers. Only a negligibly small additional reduction of the errors is obtained by including spatial variance in the filter kernel. Compared with Gaussian filtering, Catte diffusion can further reduce the error in some cases. For the examples considered, using accurate image formation models during iterative reconstruction is far more important than the choice of the filter. (author)

[en] Block-iterative image reconstruction methods, such as ordered subset expectation maximization (OSEM), are commonly used to accelerate image reconstruction. In OSEM, the speed-up factor over maximum likelihood expectation maximization (MLEM) is approximately equal to the number of subsets in which the projection data are divided. Traditionally, each subset consists of a couple of projection views, and the more subsets are used, the more the solution deviates from MLEM solutions. We found for multi-pinhole single photon emission computed tomography (SPECT) that even moderate acceleration factors in OSEM lead to inaccurate reconstructions. Therefore, we introduce pixel-based ordered subset expectation maximization (POSEM), which is based on an alternative subset choice. Pixels in each subset are spread out regularly over projections and are spatially separated as much as possible. We validated POSEM for data acquired with a focusing multi-pinhole SPECT system. Performance was compared with traditional OSEM and MLEM for a rat total body bone scan, a gated mouse myocardial perfusion scan and a Defrise phantom scan. We found that POSEM can be operated at acceleration factors that are often an order of magnitude higher than in traditional OSEM.

[en] A rotation-based Monte Carlo (MC) simulation method (RMC) has been developed, designed for rapid calculation of downscatter through non-uniform media in SPECT. A possible application is downscatter correction in dual isotope SPECT. With RMC, only a fraction of all projections of a SPECT study have to be MC simulated in a standard manner. The other projections can be estimated rapidly using the results of these standard MC calculations. For efficiency, approximations have to be made in RMC with regard to the final scatter angle of the detected photons. Further speed-up is obtained by combining RMC with convolution-based forced detection (CFD) instead of forced detection (FD), which is a more common variance reduction technique for MC. The RMC method was compared with standard MC for ^99mTc downscatter in a ²⁰¹Tl window (72 keV±10%) using a digital thorax phantom. The resulting scatter projections are in good agreement (maximum bias a few per cent of the largest value in the projection), but RMC with CFD is about three orders in magnitude faster than standard MC with FD and up to 25 times faster than standard MC with CFD. Using RMC combined with CFD, the generation of 64 almost noise-free downscatter projections (64x64) takes only a couple of minutes on a 500 MHz Pentium processor. Therefore, rotation-based Monte Carlo could serve as a practical tool for downscatter correction schemes in dual isotope SPECT. (author)

[en] The use of multi-pinhole collimation has enabled ultra-high-resolution imaging of SPECT and PET tracers in small animals. Key for obtaining high-quality images is the use of statistical iterative image reconstruction with accurate energy-dependent photon transport modelling through collimator and detector. This can be incorporated in a system matrix that contains the probabilities that a photon emitted from a certain voxel is detected at a specific detector pixel. Here we introduce a fast Monte-Carlo based (FMC-based) matrix generation method for pinhole imaging that is easy to apply to various radionuclides. The method is based on accelerated point source simulations combined with model-based interpolation to straightforwardly change or combine photon energies of the radionuclide of interest. The proposed method was evaluated for a VECTor PET-SPECT system with (i) a HE-UHR-M collimator and (ii) an EXIRAD-3D 3D autoradiography collimator. Both experimental scans with ^99mTc, ¹¹¹In, and ¹²³I, and simulated scans with ⁶⁷Ga and ⁹⁰Y were performed for evaluation. FMC was compared with two currently used approaches, one based on a set of point source measurements with ^99mTc (dubbed traditional method), and the other based on an energy-dependent ray-tracing simulation (ray-tracing method). The reconstruction results show better image quality when using FMC-based matrices than when applying the traditional or ray-tracing matrices in various cases. FMC-based matrices generalise better than the traditional matrices when imaging radionuclides with energies deviating too much from the energy used in the calibration and are computationally more efficient for very-high-resolution imaging than the ray-tracing matrices. In addition, FMC has the advantage of easily combining energies in a single matrix which is relevant when imaging radionuclides with multiple photopeak energies (e.g. ⁶⁷Ga and ¹¹¹In) or with a continuous energy spectrum (e.g. ⁹⁰Y). To conclude, FMC is an efficient, accurate, and versatile tool for creating system matrices for ultra-high-resolution pinhole SPECT. (paper)

[en] Statistical reconstruction (SR) methods provide a general and flexible framework for obtaining tomographic images from projections. For several applications SR has been shown to outperform analytical algorithms in terms of resolution-noise trade-off achieved in the reconstructions. A disadvantage of SR is the long computational time required to obtain the reconstructions, in particular when large data sets characteristic for x-ray computer tomography (CT) are involved. As was shown recently, by combining statistical methods with block iterative acceleration schemes [e.g., like in the ordered subsets convex (OSC) algorithm], the reconstruction time for x-ray CT applications can be reduced by about two orders of magnitude. There are, however, some factors lengthening the reconstruction process that hamper both accelerated and standard statistical algorithms to similar degree. In this simulation study based on monoenergetic and scatter-free projection data, we demonstrate that one of these factors is the extremely high number of iterations needed to remove artifacts that can appear around high-contrast structures. We also show (using the OSC method) that these artifacts can be adequately suppressed if statistical reconstruction is initialized with images generated by means of Radon inversion algorithms like filtered back projection (FBP). This allows the reconstruction time to be shortened by even as much as one order of magnitude. Although the initialization of the statistical algorithm with FBP image introduces some additional noise into the first iteration of OSC reconstruction, the resolution-noise trade-off and the contrast-to-noise ratio of final images are not markedly compromised

[en] Imaging of ^99mTc-labelled tracers is gaining popularity for detecting breast tumours. Recently, we proposed a novel design for molecular breast tomosynthesis (MBT) based on two sliding focusing multi-pinhole collimators that scan a modestly compressed breast. Simulation studies indicate that MBT has the potential to improve the tumour-to-background contrast-to-noise ratio significantly over state-of-the-art planar molecular breast imaging. The aim of the present paper is to optimize the collimator-detector geometry of MBT. Using analytical models, we first optimized sensitivity at different fixed system resolutions (ranging from 5 to 12 mm) by tuning the pinhole diameters and the distance between breast and detector for a whole series of automatically generated multi-pinhole designs. We evaluated both MBT with a conventional continuous crystal detector with 3.2 mm intrinsic resolution and with a pixelated detector with 1.6 mm pixels. Subsequently, full system simulations of a breast phantom containing several lesions were performed for the optimized geometry at each system resolution for both types of detector. From these simulations, we found that tumour-to-background contrast-to-noise ratio was highest for systems in the 7 mm–10 mm system resolution range over which it hardly varied. No significant differences between the two detector types were found. (paper)

[en] SPECT quantitation and image contrast are degraded by photon scatter. Water equivalent depths (WEDs) have been used by several investigators to model scatter responses in non-uniform attenuators. The drawback of this approach is the occurrence of undesired fluctuations in the shape of the scatter responses, as is shown by measurements. An improvement of the WED method is presented, based on the assumption that only a part of the scattering object (the region in the 'scatter cone') contributes significantly to the detected scatter events. The remaining part of the object is treated as a uniform medium. The extension of the WED method with extra-conical invariance (WED_ECI) is evaluated by projection measurements of a phantom with a ^99mTc source. Shapes of scatter responses predicted by the WED_ECI method are found to agree better with the measurements than t hose predicted by conventional WEDs. (author)

[en] Planar molecular breast imaging (MBI) is rapidly gaining in popularity in diagnostic oncology. To add 3D capabilities, we introduce a novel molecular breast tomosynthesis (MBT) scanner concept based on multi-pinhole collimation. In our design, the patient lies prone with the pendant breast lightly compressed between transparent plates. Integrated webcams view the breast through these plates and allow the operator to designate the scan volume (e.g. a whole breast or a suspected region). The breast is then scanned by translating focusing multi-pinhole plates and NaI(Tl) gamma detectors together in a sequence that optimizes count yield from the volume-of-interest. With simulations, we compared MBT with existing planar MBI. In a breast phantom containing different lesions, MBT improved tumour-to-background contrast-to-noise ratio (CNR) over planar MBI by 12% and 111% for 4.0 and 6.0 mm lesions respectively in case of whole breast scanning. For the same lesions, much larger CNR improvements of 92% and 241% over planar MBI were found in a scan that focused on a breast region containing several lesions. MBT resolved 3.0 mm rods in a Derenzo resolution phantom in the transverse plane compared to 2.5 mm rods distinguished by planar MBI. While planar MBI cannot provide depth information, MBT offered 4.0 mm depth resolution. Our simulations indicate that besides offering 3D localization of increased tracer uptake, multi-pinhole MBT can significantly increase tumour-to-background CNR compared to planar MBI. These properties could be promising for better estimating the position, extend and shape of lesions and distinguishing between single and multiple lesions. (paper)

[en] For the purpose of obtaining x-ray tomographic images, statistical reconstruction (SR) provides a general framework with possible advantages over analytical algorithms such as filtered backprojection (FBP) in terms of flexibility, resolution, contrast and image noise. However, SR images may be seriously affected by some artefacts that are not present in FBP images. These artefacts appear as aliasing patterns and as severe overshoots in the areas of sharp intensity transitions ('edge artefacts'). We characterize this inherent property of iterative reconstructions and hypothesize how discretization errors during reconstruction contribute to the formation of the artefacts. An adequate solution to the problem is to perform the reconstructions on an image grid that is finer than that typically employed for FBP reconstruction, followed by a downsampling of the resulting image to a granularity normally used for display. Furthermore, it is shown that such a procedure is much more effective than post-filtering of the reconstructions. Resulting SR images have superior noise-resolution trade-off compared to FBP, which may facilitate dose reduction during CT examinations

[en] We recently developed a dedicated focusing multi-pinhole collimator for a stationary SPECT system that offers down to 120 µm (or 1.7 nL) spatial resolution SPECT images of cryo-cooled tissue samples (EXIRAD-3D). This collimator is suitable for imaging isotopes that are often used in small animal and diagnostic SPECT such as ¹²⁵I (27 keV), ²⁰¹Tl (71 keV), ^99mTc (140 keV), and ¹¹¹In (171 and 245 keV). The goal of the present work is to develop high-resolution pinhole imaging of tissue samples containing isotopes with high-energy photon emissions, for example, therapeutic alpha and beta emitters that co-emit high energy gammas (e.g. ²¹³Bi (440 keV) and ¹³¹I (364 keV)) or 511 keV annihilation photons from PET isotopes. To this end, we optimise and evaluate a new high energy small-bore multi-pinhole collimator through simulations. The collimator-geometry was first optimised by simulating a Derenzo phantom scan with a biologically realistic activity concentration of ¹⁸F at two system sensitivities (0.30% and 0.60%) by varying pinhole placements. Subsequently, the wall thickness was selected based on reconstructions of a Derenzo phantom and a uniform phantom. The obtained collimators were then evaluated for ¹³¹I (364 keV), ²¹³Bi (440 keV), ⁶⁴Cu (511 keV), and ¹²⁴I (511 + 603 keV) with biologically realistic activity concentrations, and also for some high activity concentrations of ¹⁸F, using digital resolution, mouse knee joint, and xenograft phantoms. Our results show that placing pinhole centres at a distance of 8 mm from the collimator inner wall yields good image quality, while a wall thickness of 43 mm resulted in sufficient shielding. The collimators offer resolutions down to 0.35 mm, 0.6 mm, 0.5 mm, 0.6 mm, and 0.5 mm when imaging ¹³¹I, ²¹³Bi, ¹⁸F, ⁶⁴Cu, and ¹²⁴I, respectively, contained in tissue samples at biologically achievable activity concentrations. (paper)