Roel Van Holen

The effect of errors in segmented attenuation maps on PET quantification

PURPOSE Accurate attenuation correction is important for PET quantification. Often, a segmented attenuation map is used, especially in MRI-based attenuation correction. As deriving the attenuation map from MRI images is difficult, different errors can be present in the segmented attenuation map. The goal of this paper is to determine the effect of these errors on quantification. METHODS The authors simulated the digital XCAT phantom using the GATE Monte Carlo simulation framework and a model of the Philips Gemini TF. A whole body scan was simulated, spanning an axial field of view of 70 cm. A total of fifteen lesions were placed in the lung, liver, spine, colon, prostate, and femur. The acquired data were reconstructed with a reference attenuation map and with different attenuation maps that were modified to reflect common segmentation errors. The quantitative difference between reconstructed images was evaluated. RESULTS Segmentation into five tissue classes, namely cortical bone, spongeous bone, soft tissue, lung, and air yielded errors below 5%. Large errors were caused by ignoring lung tissue (up to 45%) or cortical bone (up to 17%). The interpatient variability of lung attenuation coefficients can lead to errors of 10% and more. Up to 20% tissue misclassification from bone to soft tissue yielded errors below 5%. The same applies for up to 10% misclassification from lung to air. CONCLUSIONS When using a segmented attenuation map, at least five different tissue types should be considered: cortical bone, spongeous bone, soft tissue, lung, and air. Furthermore, the interpatient variability of lung attenuation coefficients should be taken into account. Limited misclassification from bone to soft tissue and from lung to air is acceptable, as these do not lead to relevant errors.

Rapid additive manufacturing of MR compatible multipinhole collimators with selective laser melting of tungsten powder

PURPOSE The construction of complex collimators with a high number of oblique pinholes is very labor intensive, expensive or is sometimes impossible with the current available techniques (drilling, milling or electric discharge machining). All these techniques are subtractive: one starts from solid plates and the material at the position of the pinholes is removed. The authors used a novel technique for collimator construction, called metal additive manufacturing. This process starts with a solid piece of tungsten on which a first layer of tungsten powder is melted. Each subsequent layer is then melted on the previous layer. This melting is done by selective laser melting at the locations where the CAD design file defines solid material. METHODS A complex collimator with 20 loftholes with 500 μm diameter pinhole opening was designed and produced (16 mm thick and 70 × 52 mm(2) transverse size). The density was determined, the production accuracy was measured (GOM ATOS II Triple Scan, Nikon AZ100M microscope, Olympus IMT200 microscope). Point source measurements were done by mounting the collimator on a SPECT detector. Because there is increasing interest in dual-modality SPECT-MR imaging, the collimator was also positioned in a 7T MRI scanner (Bruker Pharmascan). A uniform phantom was acquired using T1, T2, and T2* sequences to check for artifacts or distortion of the phantom images due to the collimator presence. Additionally, three tungsten sample pieces (250, 500, and 750 μm thick) were produced. The density, attenuation (140 keV beam), and uniformity (GE eXplore Locus SP micro-CT) of these samples were measured. RESULTS The density of the collimator was equal to 17.31 ± 0.10 g∕cm(3) (89.92% of pure tungsten). The production accuracy ranges from -260 to +650 μm. The aperture positions have a mean deviation of 5 μm, the maximum deviation was 174 μm and the minimum deviation was -122 μm. The mean aperture diameter is 464 ± 19 μm. The calculated and measured sensitivity and resolution of point sources at different positions in the field-of-view agree well. The measured and expected attenuation of the three sample pieces are in a good agreement. There was no influence of the 7T magnetic field on the collimator (which is paramagnetic) and minimal distortion was noticed on the MR scan of the uniform phantom. CONCLUSIONS Additive manufacturing is a very promising technique for the production of complex multipinhole collimators and may also be used for producing other complex collimators. The cost of this technique is only related to the amount of powder needed and the time it takes to have the collimator built. The timeframe from design to collimator production is significantly reduced.

Iterative CT Reconstruction Using Shearlet-Based Regularization

Total variation (TV) methods have been proposed to improve the image quality in count-reduced images, by reducing the variation between neighboring pixels. Although very easy to implement and fast to compute, TV-based methods may lead to a loss of texture information when applied to images with complex textures, such as high-resolution abdominal CT images. Here, we investigate the use of another regularization approach in the context of medical images based on multiresolution transformations. One such transformation is the shearlet transform, which is optimally sparse for images that are C2 except for discontinuities along C2 curves, and has better directional sensitivity than most other, related, wavelet transform approaches. We propose to solve the convex problem using the split-Bregman (augmented Lagrangian) approach. One of the primary advantages of the split-Bregman approach, is that the shearlet transform can easily be incorporated into the sparse-view CT reconstruction. The required sparsity prior is the l1 norm of the shearlet coefficients. Results are shown for this method in comparison to the same framework with TV as the regularization term on simulated data. The noise-resolution performance is investigated at different contrast levels. At equal image noise, TV-based regularization outperforms shearlet-based regularization. However, when image texture is analyzed on measured mouse data, shearlets outperform TV, which suffers from staircasing effects. Our results show that there are benefits in using shearlets in CT imaging: texture is reconstructed more accurately compared to when TV is used, without biasing the image towards a piecewise constant image model. However, due to the larger support of the basis functions, our results suggest that uncareful usage of shearlets may lead to wavy artifacts, which can be equally unwanted as staircasing effects.

Fast and memory‐efficient Monte Carlo‐based image reconstruction for whole‐body PET

PURPOSE Several studies have shown the benefit of an accurate system modeling using Monte Carlo techniques. For state-of-the-art whole-body positron emission tomography (PET) scanners, Monte Carlo-based image reconstruction is associated with a significant computational cost to calculate the system matrix as well as a large memory capacity to store it. In this article, the authors present a simulation-reconstruction framework to solve these problems on the Philips Gemini GS PET scanner. METHODS A fast, realistic system matrix simulation module was developed using egs_pet, which is an efficient PET simulation code based on EGSnrc. The generated system matrix was then used in a rotator-based ordered subset expectation maximization (OS-EM) algorithm, which exploits the rotational symmetry of a cylindrical PET scanner. The system matrix was further compressed by using sparse storage techniques. RESULTS The system matrix simulation took five days on 50 cores of Xeon 2.66 GHz, resulting in a system matrix of 2.01 GB. The entire system matrix could be stored in the main memory of a standard personal computer. The image quality in terms of contrast-noise trade-offs was considerably improved compared to a standard OS-EM algorithm. The image quality was also compared to the clinical software on the scanner using routine parameter settings. The contrast recovery coefficient of small hot spheres and cold spheres was significantly improved. CONCLUSIONS The results indicated that the proposed framework could be used for this PET scanner with improved image quality. This method could also be applied to other state-of-the-art whole-body PET scanners and preclinical PET scanners with a similar shape.

Fast simulation of yttrium-90 bremsstrahlung photons with GATE

PURPOSE Multiple investigators have recently reported the use of yttrium-90 (Y90) bremsstrahlung single photon emission computed tomography (SPECT) imaging for the dosimetry of targeted radionuclide therapies. Because Monte Carlo (MC) simulations are useful for studying SPECT imaging, this study investigates the MC simulation of Y90 bremsstrahlung photons in SPECT. To overcome the computationally expensive simulation of electrons, the authors propose a fast way to simulate the emission of Y90 bremsstrahlung photons based on prerecorded bremsstrahlung photon probability density functions (PDFs). METHODS The accuracy of bremsstrahlung photon simulation is evaluated in two steps. First, the validity of the fast bremsstrahlung photon generator is checked. To that end, fast and analog simulations of photons emitted from a Y90 point source in a water phantom are compared. The same setup is then used to verify the accuracy of the bremsstrahlung photon simulations, comparing the results obtained with PDFs generated from both simulated and measured data to measurements. In both cases, the energy spectra and point spread functions of the photons detected in a scintillation camera are used. RESULTS Results show that the fast simulation method is responsible for a 5% overestimation of the low-energy fluence (below 75 keV) of the bremsstrahlung photons detected using a scintillation camera. The spatial distribution of the detected photons is, however, accurately reproduced with the fast method and a computational acceleration of ∼17-fold is achieved. When measured PDFs are used in the simulations, the simulated energy spectrum of photons emitted from a point source of Y90 in a water phantom and detected in a scintillation camera closely approximates the measured spectrum. The PSF of the photons imaged in the 50-300 keV energy window is also accurately estimated with a 12.4% underestimation of the full width at half maximum and 4.5% underestimation of the full width at tenth maximum. CONCLUSIONS Despite its limited accuracy, the fast bremsstrahlung photon generator is well suited for the simulation of bremsstrahlung photons emitted in large homogeneous organs, such as the liver, and detected in a scintillation camera. The computational acceleration makes it very useful for future investigations of Y90 bremsstrahlung SPECT imaging.

Review of SPECT collimator selection, optimization, and fabrication for clinical and preclinical imaging

In single photon emission computed tomography, the choice of the collimator has a major impact on the sensitivity and resolution of the system. Traditional parallel-hole and fan-beam collimators used in clinical practice, for example, have a relatively poor sensitivity and subcentimeter spatial resolution, while in small-animal imaging, pinhole collimators are used to obtain submillimeter resolution and multiple pinholes are often combined to increase sensitivity. This paper reviews methods for production, sensitivity maximization, and task-based optimization of collimation for both clinical and preclinical imaging applications. New opportunities for improved collimation are now arising primarily because of (i) new collimator-production techniques and (ii) detectors with improved intrinsic spatial resolution that have recently become available. These new technologies are expected to impact the design of collimators in the future. The authors also discuss concepts like septal penetration, high-resolution applications, multiplexing, sampling completeness, and adaptive systems, and the authors conclude with an example of an optimization study for a parallel-hole, fan-beam, cone-beam, and multiple-pinhole collimator for different applications.

Comparing planar image quality of rotating slat and parallel hole collimation: influence of system modeling

The main remaining challenge for a gamma camera is to overcome the existing trade-off between collimator spatial resolution and system sensitivity. This problem, strongly limiting the performance of parallel hole collimated gamma cameras, can be overcome by applying new collimator designs such as rotating slat (RS) collimators which have a much higher photon collection efficiency. The drawback of a RS collimated gamma camera is that, even for obtaining planar images, image reconstruction is needed, resulting in noise accumulation. However, nowadays iterative reconstruction techniques with accurate system modeling can provide better image quality. Because the impact of this modeling on image quality differs from one system to another, an objective assessment of the image quality obtained with a RS collimator is needed in comparison to classical projection images obtained using a parallel hole (PH) collimator. In this paper, a comparative study of image quality, achieved with system modeling, is presented. RS data are reconstructed to planar images using maximum likelihood expectation maximization (MLEM) with an accurate Monte Carlo derived system matrix while PH projections are deconvolved using a Monte Carlo derived point-spread function. Contrast-to-noise characteristics are used to show image quality for cold and hot spots of varying size. Influence of the object size and contrast is investigated using the optimal contrast-to-noise ratio (CNR(o)). For a typical phantom setup, results show that cold spot imaging is slightly better for a PH collimator. For hot spot imaging, the CNR(o) of the RS images is found to increase with increasing lesion diameter and lesion contrast while it decreases when background dimensions become larger. Only for very large background dimensions in combination with low contrast lesions, the use of a PH collimator could be beneficial for hot spot imaging. In all other cases, the RS collimator scores better. Finally, the simulation of a planar bone scan on a RS collimator revealed a hot spot contrast improvement up to 54% compared to a classical PH bone scan.

System characteristics of SPECT with a slat collimated strip detector.

In classical SPECT with parallel hole collimation, the sensitivity is constant over the field of view (FOV). This is no longer the case if a rotating slat collimator with planar photon collection is used: there will be a significant variation of the sensitivity within the FOV. Since not compensating for this inhomogeneous sensitivity distribution would result in non-quantitative images, an accurate knowledge of the sensitivity is mandatory to account for it during reconstruction. On the other hand, the spatial resolution versus distance dependency remains unaltered compared to parallel hole collimation. For deriving the sensitivity, different factors have to be taken into account: a first factor concerns the intrinsic detector properties and will be incorporated into the calculations as a detection efficiency term depending on the incident angle. The calculations are based on a second and more pronounced factor: the collimator and detector geometry. Several assumptions will be made for the calculation of the sensitivity formulae and it will be proven that these calculations deliver a valid prediction of the sensitivity at points far enough from the collimator. To derive a close field model which also accounts for points close to the collimator surface, a modified calculation method is used. After calculating the sensitivity in one plane it is easy to obtain the tomographic sensitivity. This is done by rotating the sensitivity maps for spin and camera rotation. The results derived from the calculations are then compared to simulation results and both show good agreement after including the aforementioned detection efficiency term. The validity of the calculations is also proven by measuring the sensitivity in the FOV of a prototype rotating slat gamma camera. An expression for the resolution of these planar collimation systems is obtained. It is shown that for equal collimator dimensions the same resolution-distance relationship is obtained as for parallel hole collimators. Although, a better spatial resolution can be obtained with our prototype camera due to the smaller pitch of the slats. This can be achieved without a major drop in system sensitivity due to the fact that the slats consist of less collimator material compared to a parallel hole collimator. The accuracy of the calculated resolution is proven by comparison with Monte Carlo simulation and measurement resolution values.

Characterization of a SPECT pinhole collimator for optimal detector usage (the lofthole)

In single-photon emission computed tomography (SPECT), multi-pinhole collimation is often employed nowadays. Most multi-pinhole collimators avoid overlap (multiplexing) of the projections on the detector. This can be done by using additional shielding or by spacing the pinholes far enough apart. Using additional shielding has the drawback that it increases weight, design complexity and cost. Spacing the pinholes far enough apart results in sub-optimal detector usage, the valuable detector area is not entirely used. This is due to the circular projections of pinholes on the detector; these ellipses can not be tiled with high detector coverage. To overcome this we designed a new pinhole geometry, the lofthole, that has a rectangular projection on the detector. The lofthole has a circular aperture and a rectangular entrance/exit opening. Sensitivity formulae have been derived for pinholes and loftholes. These formulae take the penumbra effect into account; the proposed formulae do not take penetration into account. The derived formulae are valid for geometries where the field-of-view and the sensitivity of the aperture are solely limited by the exit window. A flood map measurement was performed to compare the rectangular projection of a lofthole with the circular projection of a pinhole. Finally, measurements were done to compare the amount of penetration of pinholes with the amount of penetration of a lofthole. A square lofthole collimator has less penetration than a knife-edge pinhole collimator that irradiates the same rectangular detector area with full coverage. A multi-lofthole collimator allows high detector coverage without using additional shielding. An additional advantage is the lower amount of penetration.

Design and performance of a compact and stationary microSPECT system.

PURPOSE Over the last ten years, there has been an extensive growth in the development of microSPECT imagers. Most of the systems are based on the combination of conventional, relatively large gamma cameras with poor intrinsic spatial resolution and multipinhole collimators working in large magnification mode. Spatial resolutions range from 0.58 to 0.76 mm while peak sensitivities vary from 0.06% to 0.4%. While pushing the limits of performance is of major importance, the authors believe that there is a need for smaller and less complex systems that bring along a reduced cost. While low footprint and low-cost systems can make microSPECT available to more researchers, the ease of operation and calibration and low maintenance cost are additional factors that can facilitate the use of microSPECT in molecular imaging. In this paper, the authors simulate the performance of a microSPECT imager that combines high space-bandwidth detectors and pinholes with truncated projection, resulting in a small and stationary system. METHODS A system optimization algorithm is used to determine the optimal SPECT systems, given our high resolutions detectors and a fixed field-of-view. These optimal system geometries are then used to simulate a Defrise disk phantom and a hot rod phantom. Finally, a MOBY mouse phantom, with realistic concentrations of Tc99m-tetrofosmin is simulated. RESULTS Results show that the authors can successfully reconstruct a Defrise disk phantom of 24 mm in diameter without any rotating system components or translation of the object. Reconstructed spatial resolution is approximately 800 μm while the peak sensitivity is 0.23%. Finally, the simulation of the MOBY mouse phantom shows that the authors can accurately reconstruct mouse images. CONCLUSIONS These results show that pinholes with truncated projections can be used in small magnification or minification mode to obtain a compact and stationary microSPECT system. The authors showed that they can reach state-of-the-art system performance and can successfully reconstruct images with realistic noise levels in a preclinical context. Such a system can be useful for dynamic SPECT imaging.