Adaptive Angular Sampling for SPECT Imaging

Joint Regional Uptake Quantification of Thorium-227 and Radium-223 Using a Multiple-Energy-Window Projection-Domain Quantitative SPECT Method

Thorium-227 ( )-based -particle radiopharmaceutical therapies ( -RPTs) are currently being investigated in several clinical and pre-clinical studies. After administration, decays to , another -particle-emitting isotope, which redistributes within the patient. Reliable dose quantification of both and is clinically important, and SPECT may perform this quantification as these isotopes also emit X- and -ray photons. However, reliable quantification is challenging for several reasons: the orders-of-magnitude lower activity compared to conventional SPECT, resulting in a very low number of detected counts, the presence of multiple photopeaks, substantial overlap in the emission spectra of these isotopes, and the image-degrading effects in SPECT. To address these issues, we propose a multiple-energy-window projection-domain quantification (MEW-PDQ) method that jointly estimates the regional activity uptake of both and directly using the SPECT projection data from multiple energy windows. We evaluated the method with realistic simulation studies conducted with anthropomorphic digital phantoms, including a virtual imaging trial, in the context of imaging patients with bone metastases of prostate cancer who were treated with -based -RPTs. The proposed method yielded reliable (accurate and precise) regional uptake estimates of both isotopes and outperformed state-of-the-art methods across different lesion sizes and contrasts, as well as in the virtual imaging trial. This reliable performance was also observed with moderate levels of intra-regional heterogeneous uptake as well as when there were moderate inaccuracies in the definitions of the support of various regions. Additionally, we demonstrated the effectiveness of using multiple energy windows and the variance of the estimated uptake using the proposed method approached the Cramér-Rao-lower-bound-defined theoretical limit. These results provide strong evidence in support of this method for reliable uptake quantification in -based -RPTs.

Pareto optimization of SPECT acquisition and reconstruction settings for 177Lu activity quantification

BACKGROUND

The aim was to investigate the noise and bias properties of quantitative 177Lu-SPECT with respect to the number of projection angles, and the number of subsets and iterations in the OS-EM reconstruction, for different total acquisition times.

METHODS

Experimental SPECT acquisition of six spheres in a NEMA body phantom filled with 177Lu was performed, using medium-energy collimators and 120 projections with 180 s per projection. Bootstrapping was applied to generate data sets representing acquisitions with 20 to 120 projections for 10 min, 20 min, and 40 min, with 32 noise realizations per setting. Monte Carlo simulations were performed of 177Lu-DOTA-TATE in an anthropomorphic computer phantom with three tumours (2.8 mL to 40.0 mL). Projections representing 24 h and 168 h post administration were simulated, each with 32 noise realizations. Images were reconstructed using OS-EM with compensation for attenuation, scatter, and distance-dependent resolution. The number of subsets and iterations were varied within a constrained range of the product number of iterations × number of projections ≤ 2400 . Volumes-of-interest were defined following the physical size of the spheres and tumours, the mean activity-concentrations estimated, and the absolute mean relative error and coefficient of variation (CV) over noise realizations calculated. Pareto fronts were established by analysis of CV versus mean relative error.

RESULTS

Points at the Pareto fronts with low CV and high mean error resulted from using a low number of subsets, whilst points at the Pareto fronts associated with high CV but low mean error resulted from reconstructions with a high number of subsets. The number of projection angles had limited impact.

CONCLUSIONS

For accurate estimation of the 177Lu activity-concentration from SPECT images, the number of projection angles has limited importance, whilst the total acquisition time and the number of subsets and iterations are parameters of importance.

Inclusion of quasi-vertex views in a brain-dedicated multi-pinhole SPECT system for improved imaging performance

With brain-dedicated multi-detector systems employing pinhole apertures the usage of detectors facing the top of the patient's head (i.e. quasi-vertex (QV) views) can provide the advantage of additional viewing from close to the brain for improved detector coverage. In this paper, we report the results of simulation and reconstruction studies to investigate the impact of the QV views on the imaging performance of AdaptiSPECT-C, a brain-dedicated stationary SPECT system under development. In this design, both primary and scatter photons from regions located inferior to the brain can contribute to SPECT projections acquired by the QV views, and thus degrade AdaptiSPECT-C imaging performance. In this work, we determined the proportion, origin, and nature (i.e. primary, scatter, and multiple-scatter) of counts emitted from structures within the head and throughout the body contributing to projections from the different AdaptiSPECT-C detector rings, as well as from a true vertex view detector. We simulated phantoms used to assess different aspects of image quality (i.e. uniform activity concentration sphere, and Derenzo), as well as anthropomorphic phantoms with different count levels emulating clinical ¹²³I activity distributions (i.e. DaTscan and perfusion). We determined that attenuation and scatter in the patient's body greatly diminish the probability of the photons emitted outside the volume of interest reaching to detectors and being recorded within the 15% photopeak energy window. In addition, we demonstrated that the inclusion of the residual of such counts in the system acquisition does not degrade visual interpretation or quantitative analysis. The addition of the QV detectors improves volumetric sensitivity, angular sampling, and spatial resolution leading to significant enhancement in image quality, especially in the striato-thalamic and superior regions of the brain. Besides, the use of QV detectors improves the recovery of clinically relevant metrics such as the striatal binding ratio and mean activity in selected cerebral structures. Our findings proving the usefulness of the QV ring for brain imaging with ¹²³I agents can be generalized to other commonly used SPECT imaging agents labelled with isotopes, such as ^99mTc and likely ¹¹¹In.

Investigation of Axial and Angular Sampling in Multi-Detector Pinhole-SPECT Brain Imaging

We designed a dedicated multi-detector multi-pinhole brain SPECT scanner to generate images of higher quality compared to general-purpose systems. The system, AdaptiSPECT-C, is intended to adapt its sensitivity-resolution trade-off by varying its aperture configurations allowing both high-sensitivity dynamic and high-spatial-resolution static imaging. The current system design consists of 23 detector heads arranged in a truncated spherical geometry. In this work, we investigated the axial and angular sampling capability of the current stationary system design. Two data acquisition schemes using limited rotation of the gantry and two others using axial translation of the imaging bed were also evaluated concerning their impact on image quality through improved sampling. Increasing both angular and axial sampling in the current prototype system resulted in quantitative improvements in image quality metrics and qualitative appearance of the images as determined in studies with specifically selected phantoms. Visual improvements for the brain phantoms with clinical distributions were less pronounced but presented quantitative improvements in the fidelity (normalized root-mean-square error (NRMSE)) and striatal specific binding ratio (SBR) for a dopamine transporter (DAT) distribution, and in NRMSE and activity recovery for a brain perfusion distribution. More pronounced improvements with increased sampling were seen in contrast recovery coefficient, bias, and coefficient of variation for a lesion in the brain perfusion distribution. The negligible impact of the most cranial ring of detectors on axial sampling, but its significant impact on sensitivity and angular sampling in the cranial portion of the imaging volume-of-interest were also determined.

Fast quantitative reconstruction with focusing collimators for liver SPECT

BACKGROUND

Generation of a SPECT scan during procedure may aid in the optimization of treatments as liver radioembolization by offering image-guided dosimetry. This, however, requires both shortened acquisition times and fast quantitative reconstruction. Focusing collimators increase sensitivity and thus may speed up imaging. Monte Carlo-based iterative reconstruction has shown to provide quantitative results for parallel hole collimators but may be slow. The purpose of this work is to develop fast Monte Carlo-based reconstruction for focusing collimators and to evaluate the impact of reconstruction and collimator choice on quantitative accuracy of liver dosimetry by means of simulations.

RESULTS

The developed fast Monte Carlo simulator was found to accurately generate projections compared to a full Monte Carlo simulation, providing projections in several seconds instead of several days. Monte Carlo-based scatter correction was superior to other scatter correction methods in describing recovered activity and reached similar noise levels as dual-energy window scatter correction. Although truncation artifacts were present in the cone beam collimator (50 cm), the region inside the field of view (FOV) could be reconstructed without loss of accuracy. Provided the object to image is inside the FOV, the focusing collimator with 50 cm focal distance could retrieve the same noise levels as a parallel hole collimator in 68% of the total scanning time, the multifocal collimator in 73% of the time, and the 100-cm focal distance collimator in 84% of the time.

CONCLUSION

Focusing collimators combined with Monte Carlo-based reconstruction have the ability to enable quantitative imaging of the FOV in a significantly shorter timeframe. The proposed approach to the forward projector will additionally make it possible to reconstruct within minutes. These are crucial steps in moving toward real-time dosimetry during interventions.

Simulation study of the second-generation MR-compatible SPECT system based on the inverted compound-eye gamma camera design

In this paper, we present simulation studies for the second-generation MRI compatible SPECT system, MRC-SPECT-II, based on an inverted compound eye (ICE) gamma camera concept. The MRC-SPECT-II system consists of a total of 1536 independent micro-pinhole-camera-elements (MCEs) distributed in a ring with an inner diameter of 6 cm. This system provides a FOV of 1 cm diameter and a peak geometrical efficiency of approximately 1.3% (the typical levels of 0.1%-0.01% found in modern pre-clinical SPECT instrumentations), while maintaining a sub-500 μm spatial resolution. Compared to the first-generation MRC-SPECT system (MRC-SPECT-I) (Cai 2014 Nucl. Instrum. Methods Phys. Res. A 734 147-51) developed in our lab, the MRC-SPECT-II system offers a similar resolution with dramatically improved sensitivity and greatly reduced physical dimension. The latter should allow the system to be placed inside most clinical and pre-clinical MRI scanners for high-performance simultaneous MRI and SPECT imaging.

Optimization of an Adaptive SPECT System with the Scanning Linear Estimator

A method for optimization of an adaptive Single Photon Emission Computed Tomography (SPECT) system is presented. Adaptive imaging systems can quickly change their hardware configuration in response to data being generated in order to improve image quality for a specific task. In this work we simulate an adaptive SPECT system and propose a method for finding the adaptation that maximizes the performance on a signal estimation task. To start with, a simulated object model containing a spherical signal is imaged with a scout configuration. A Markov-Chain Monte Carlo (MCMC) technique utilizes the scout data to generate an ensemble of possible objects consistent with the scout data. This object ensemble is imaged by numerous simulated hardware configurations and for each system estimates of signal activity, size and location are calculated via the Scanning Linear Estimator (SLE). A figure of merit, based on a Modified Dice Index (MDI), quantifies the performance of each imaging configuration and it allows for optimization of the adaptive SPECT. This figure of merit is calculated by multiplying two terms: the first term uses the definition of the Dice similarity index to determine the percent of overlap between the actual and the estimated spherical signal, the second term utilizes an exponential function that measures the squared error for the activity estimate. The MDI combines the error in estimates of activity, size, and location, in one convenient metric and it allows for simultaneous optimization of the SPECT system with respect to all the estimated signal parameters. The results of our optimizations indicate that the adaptive system performs better than a non-adaptive one in conditions where the diagnostic scan has a low photon count - on the order of thousand photons per projection. In a statistical study, we optimized the SPECT system for one hundred unique objects and demonstrated that the average MDI on an estimation task is 0.84 for the adaptive system and 0.65 for the non-adaptive system.

Optimal channels for channelized quadratic estimators

We present a new method for computing optimized channels for estimation tasks that is feasible for high-dimensional image data. Maximum-likelihood (ML) parameter estimates are challenging to compute from high-dimensional likelihoods. The dimensionality reduction from M measurements to L channels is a critical advantage of channelized quadratic estimators (CQEs), since estimating likelihood moments from channelized data requires smaller sample sizes and inverting a smaller covariance matrix is easier. The channelized likelihood is then used to form ML estimates of the parameter(s). In this work we choose an imaging example in which the second-order statistics of the image data depend upon the parameter of interest: the correlation length. Correlation lengths are used to approximate background textures in many imaging applications, and in these cases an estimate of the correlation length is useful for pre-whitening. In a simulation study we compare the estimation performance, as measured by the root-mean-squared error (RMSE), of correlation length estimates from CQE and power spectral density (PSD) distribution fitting. To abide by the assumptions of the PSD method we simulate an ergodic, isotropic, stationary, and zero-mean random process. These assumptions are not part of the CQE formalism. The CQE method assumes a Gaussian channelized likelihood that can be a valid for non-Gaussian image data, since the channel outputs are formed from weighted sums of the image elements. We have shown that, for three or more channels, the RMSE of CQE estimates of correlation length is lower than conventional PSD estimates. We also show that computing CQE by using a standard nonlinear optimization method produces channels that yield RMSE within 2% of the analytic optimum. CQE estimates of anisotropic correlation length estimation are reported to demonstrate this technique on a two-parameter estimation problem.

Evaluation of Fisher Information Matrix-Based Methods for Fast Assessment of Image Quality in Pinhole SPECT

The accurate determination of the local impulse response and the covariance in voxels from penalized maximum likelihood reconstructed images requires performing reconstructions from many noise realizations of the projection data. As this is usually a very time-consuming process, efficient analytical approximations based on the Fisher information matrix (FIM) have been extensively used in PET and SPECT to estimate these quantities. For 3D imaging, however, additional approximations need to be made to the FIM in order to speed up the calculations. The most common approach is to use the local shift-invariant (LSI) approximation of the FIM, but this assumes specific conditions which are not always necessarily valid. In this paper we take a single-pinhole SPECT system and compare the accuracy of the LSI approximation against two other methods that have been more recently put forward: the non-uniform object-space pixelation (NUOP) and the subsampled FIM. These methods do not assume such restrictive conditions while still increasing the speed of the calculations considerably. Our results indicate that in pinhole SPECT the NUOP and subsampled FIM approaches could be more reliable than the LSI approximation, especially when a high accuracy is required.

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.

Comparison of FBP and Iterative Algorithms with Non-Uniform Angular Sampling

Some people may believe that the filtered backprojection (FBP) algorithm does not work if the projection data are measured non-uniformly. Some may also believe that iterative algorithms can automatically handle the non-uniformly sampled data in the projector/backprojector. This paper claims that the FBP algorithm can effectively handle the situation where the angular sampling is not uniform. This paper compares the images that are reconstructed by both the FBP and the iterative Landweber algorithms when the angular sampling is nonuniform. When the iteration number is low, the iterative algorithms do not handle the non-uniform sampling properly. A weighting strategy is then suggested and it makes the image resolution more isotropic. In few-view tomography, the FBP and iterative algorithms both perform poorly if no other prior information is used. We have made the following observations: 1) When using an iterative algorithm, one must use early solutions due to noise amplification. 2) An early solution can have anisotropic spatial resolution if the angular sampling is not uniform. 3) The anisotropic resolution problem can be solved by introducing angle dependent weighting, which is not noise dependent. 4) The weighting is not effective when the iteration number is large. The weighting only affects the early solutions, and does not affect the converged solution. 5) When the iteration number is large, the model-mismatch errors are amplified and cause artifacts in the image. 6) The FBP algorithm is not sensitive to the model-mismatch errors, and does not have the "early solution" problems. 7) In few-view tomography, both FBP and iterative algorithms perform poorly, while the FBP algorithm gives a sharper image than the iterative algorithm does.

Efficient determination of the uncertainty for the optimization of SPECT system design: a subsampled fisher information matrix

System designs in single photon emission tomography (SPECT) can be evaluated based on the fundamental trade-off between bias and variance that can be achieved in the reconstruction of emission tomograms. This trade off can be derived analytically using the Cramer-Rao type bounds, which imply the calculation and the inversion of the Fisher information matrix (FIM). The inverse of the FIM expresses the uncertainty associated to the tomogram, enabling the comparison of system designs. However, computing, storing and inverting the FIM is not practical with 3-D imaging systems. In order to tackle the problem of the computational load in calculating the inverse of the FIM, a method based on the calculation of the local impulse response and the variance, in a single point, from a single row of the FIM, has been previously proposed for system design. However this approximation (circulant approximation) does not capture the global interdependence between the variables in shift-variant systems such as SPECT, and cannot account e.g., for data truncation or missing data. Our new formulation relies on subsampling the FIM. The FIM is calculated over a subset of voxels arranged in a grid that covers the whole volume. Every element of the FIM at the grid points is calculated exactly, accounting for the acquisition geometry and for the object. This new formulation reduces the computational complexity in estimating the uncertainty, but nevertheless accounts for the global interdependence between the variables, enabling the exploration of design spaces hindered by the circulant approximation. The graphics processing unit accelerated implementation of the algorithm reduces further the computation times, making the algorithm a good candidate for real-time optimization of adaptive imaging systems. This paper describes the subsampled FIM formulation and implementation details. The advantages and limitations of the new approximation are explored, in comparison with the circulant approximation, in the context of design optimization of a parallel-hole collimator SPECT system and of an adaptive imaging system (similar to the commercially available D-SPECT).

Performance assessment of the single photon emission microscope: high spatial resolution SPECT imaging of small animal organs

The single photon emission microscope (SPEM) is an instrument developed to obtain high spatial resolution single photon emission computed tomography (SPECT) images of small structures inside the mouse brain. SPEM consists of two independent imaging devices, which combine a multipinhole collimator, a high-resolution, thallium-doped cesium iodide [CsI(Tl)] columnar scintillator, a demagnifying/intensifier tube, and an electron-multiplying charge-coupling device (CCD). Collimators have 300- and 450-µm diameter pinholes on tungsten slabs, in hexagonal arrays of 19 and 7 holes. Projection data are acquired in a photon-counting strategy, where CCD frames are stored at 50 frames per second, with a radius of rotation of 35 mm and magnification factor of one. The image reconstruction software tool is based on the maximum likelihood algorithm. Our aim was to evaluate the spatial resolution and sensitivity attainable with the seven-pinhole imaging device, together with the linearity for quantification on the tomographic images, and to test the instrument in obtaining tomographic images of different mouse organs. A spatial resolution better than 500 µm and a sensitivity of 21.6 counts·s^-1·MBq^-1 were reached, as well as a correlation coefficient between activity and intensity better than 0.99, when imaging ^99mTc sources. Images of the thyroid, heart, lungs, and bones of mice were registered using ^99mTc-labeled radiopharmaceuticals in times appropriate for routine preclinical experimentation of <1 h per projection data set. Detailed experimental protocols and images of the aforementioned organs are shown. We plan to extend the instrument's field of view to fix larger animals and to combine data from both detectors to reduce the acquisition time or applied activity.

SmartCAM: An adaptive clinical SPECT camera

An adaptive pinhole aperture that fits a GE MaxiCam Single-Photon-Emission Computed Tomography (SPECT) system has been designed, built, and is undergoing testing. The purpose of an adaptive aperture is to allow the imaging system to make adjustments to the aperture while imaging data are being acquired. Our adaptive pinhole aperture can alter several imaging parameters, including field of view, resolution, sensitivity, and magnification. The dynamic nature of such an aperture allows for imaging of specific regions of interest based on initial measurements of the patient. Ideally, this mode of data collection will improve the understanding of a patient's condition, and will facilitate better diagnosis and treatment. The aperture was constructed using aluminum and a low melting point, high-stopping-power metal alloy called Cerrobend. The aperture utilizes a rotating disk for the selection of a pinhole configuration; as the aluminum disk rotates, different pinholes move into view of the camera face and allow the passage of gamma rays through that particular pinhole. By controlling the angular position of the disk, the optical characteristics of the aperture can be modified, allowing the system to acquire data from controlled regions of interest. First testing was performed with a small radioactive source to prove the functionality of the aperture.

SPECT system optimization against a discrete parameter space

In this paper, we present an analytical approach for optimizing the design of a static SPECT system or optimizing the sampling strategy with a variable/adaptive SPECT imaging hardware against an arbitrarily given set of system parameters. This approach has three key aspects. First, it is designed to operate over a discretized system parameter space. Second, we have introduced an artificial concept of virtual detector as the basic building block of an imaging system. With a SPECT system described as a collection of the virtual detectors, one can convert the task of system optimization into a process of finding the optimum imaging time distribution (ITD) across all virtual detectors. Thirdly, the optimization problem (finding the optimum ITD) could be solved with a block-iterative approach or other nonlinear optimization algorithms. In essence, the resultant optimum ITD could provide a quantitative measure of the relative importance (or effectiveness) of the virtual detectors and help to identify the system configuration or sampling strategy that leads to an optimum imaging performance. Although we are using SPECT imaging as a platform to demonstrate the system optimization strategy, this development also provides a useful framework for system optimization problems in other modalities, such as positron emission tomography and x-ray computed tomography (Moore et al (2009 IEEE Nucl. Sci. Symp. Conf. Rec. pp 4154-7), Freed et al (2008 Med. Phys. 35 1912-25)).