Pinhole SPECT Research Articles

Efficient Monte-Carlo based system modelling for image reconstruction in preclinical pinhole SPECT

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, 111In, and 123I, and simulated scans with 67Ga and 90Y 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. 67Ga and 111In) or with a continuous energy spectrum (e.g. 90Y). To conclude, FMC is an efficient, accurate, and versatile tool for creating system matrices for ultra-high-resolution pinhole SPECT.

Super Spatial Resolution (SSR) method for small animal SPECT imaging: A Monte Carlo study

This work used Monte Carlo (MC) simulations in order to investigate the feasibility of the application of a super-resolution method, namely Super Spatial Resolution (SSR), to a preclinical SPECT system. The SSR method is used to obtain a resolution enhancement by combining several images containing slightly different perspectives of the same scene.MC simulations are a useful tool in the development of new imaging devices, as well as in the assessment of new reconstruction methods. In this respect, GATE (Geant4 Application for Emission Tomography) is an open source platform dedicated to numerical MC simulations that has become a well validated toolkit for the simulation of SPECT and PET systems.We modelled a four-headed preclinical SPECT scanner capable of the proper movements in order to obtain the SSR acquisition sequences. Each detector is based on a CsI(Tl) pixelated scintillator coupled to a low-energy tungsten collimator with parallel square holes.Several phantoms and radioactive sources were simulated, in order to assess the performance of the proposed SPECT scanner. In addition, a complex voxelized phantom of a mouse was used to obtain realistic images with the aim of evaluating the system as a suitable tool for small animal imaging.The whole system sensitivity was 164.1 cps/MBq. In order to assess the impact of the SSR this value must be divided by the number of the images acquired to perform this method. The average value of the trans-axial resolution improves from 2.4 mm to 1.54 mm, 1.21 mm and 1.03 mm by respectively applying the SSR based on two, three or four images. While the average axial resolution changes from 1.69 mm to 1.49 mm, 1.15 mm and 0.98 mm by respectively using the SSR based on two, three or four images. The mouse images obtained by using the voxel phantom have demonstrated the good capability of the system as a suitable tool for small animal imaging.Finally, a comparison with commercial preclinical scanners has proved that the presented SSR scanner provides an alternative to pinhole SPECT systems for many preclinical research studies.

Geometric Sensitivity of Variable Pinhole SPECT with a Keel-edge Pinhole Model

A variable pinhole (VP) collimator is formed by stacking several tungsten plates with different apertures. The configuration of the VP can be optimized for target region-of-interests at different rotation angles. Hence, both the resolution and sensitivity of the system can be improved. However, it is difficult to analytically implement the exact system model that simulates the step shapes between adjacent layers. In this study, we compared the conventional keel-edge and knife-edge pinhole models with the VP model using simulations. The keel-edge and knife-edge models were evaluated based on various specifications to identify the optimum system model. The simulation results showed that the keel-edge model with 80% channel length (keel-edge 80%) elicits the best performance with a fixed acceptance angle and pinhole plate thickness. The results in this study are expected to be applied to full-scale SPECT with a VP collimator which is currently under development.

Optimized image acquisition for dopamine transporter imaging with ultra-high resolution clinical pinhole SPECT

SPECT can be used to image dopamine transporter (DaT) availability in the human striatum, e.g. for diagnosis of Parkinson’s disease (PD). As traditional SPECT provides limited resolution and sensitivity, we proposed a full ring focusing multi-pinhole SPECT system (G-SPECT-I (Beekman 2015 Eur. J. Nucl. Med. Mol. Imaging 42 S209)) which demonstrated a 2.5 mm reconstructed resolution in phantom scans. G-SPECT-I achieves data completeness in the scan region of interest by translating the patient bed with an xyz-stage and combining projections from all bed positions into image reconstruction using a scanning focus method (SFM). This paper aims to develop dedicated SFM parameters for performing a DaTscan with high effective sensitivity and appropriate sampling. To this end, the axial scanning length was restricted and transaxial bed trajectories with a reduced number of positions based on a convex hull data-completeness model were tested. Quantitative accuracy was assessed using full G-SPECT-I simulations of an Alderson phantom based on measured system matrices. For each sampling strategy, the specific binding ratio (SBR) and asymmetry index (AI) in the left and right striatum, as well as the Localized SBR (L-SBR) and the Localized AI (L-AI) in eight striatal sub-regions were calculated and compared to those of the reference scan which performs full brain oversampling using 112 bed positions. Results show that structures essential for PD diagnosis were visually and quantitatively barely affected even when using the lowest number of bed translations (i.e. 4). The maximum deviation from the reference was only 1.5%, 1.5%, 5.5% and 7.0% for the SBR, AI, L-SBR and L-AI, respectively, when 4 positions were used. Thus, it is possible to perform an accurate DaTscan with a confined axial scan region and a limited number of focused bed positions. This enables protocols for extremely fast dynamic SPECT scans with less than half-minute time frames, which can be useful for motion correction.

Design and evaluation of two multi-pinhole collimators for brain SPECT

SPECT is a powerful tool for diagnosing or staging brain diseases such as Alzheimer's disease (AD) and Parkinson's disease (PD) but is limited by its inferior resolution and sensitivity. At the same time, pinhole SPECT provides superior resolution and detection efficiency trade-off as compared to the conventional parallel-hole collimator for imaging small field-of-view (FOV), which fits for the case of brain imaging. In this study, we propose to develop and evaluate two multi-pinhole (MPH) collimator designsto improve the imaging of cerebral blood flow and striatum. We set the target resolutions to be 12 and 8mm, respectively, and the FOV at 200mm which is large enough to cover the whole brain. The constraints for system optimization include maximum and minimum detector-to-center-of-FOV (CFOV) distances of 344 and 294mm, respectively, and minimal radius-of-rotation (ROR) of 135mm to accommodate patients' shoulder. According to the targeted FOV, resolutions, and constraints, we determined the pinhole number, ROR, focal length, apertureacceptance angle, and aperturediameter which maximized the system sensitivity. We then assessed the imaging performance of the proposedMPH and standard low-energy high-resolution (LEHR) collimators using analytical simulations of a digital NCAT brain phantom with 99mTc-HMPAO/99mTc-TRODAT-1 distributions; Monte Carlo simulations of a hot-rod phantom; and a Defrise phantom using GATE v6.1. Projections were generated over 360° and reconstructed using the 3D MPH/LEHR OS-EM methods with up to 720 updates. The normalized mean square error (NMSE) was calculated over the cerebral and striatal regions extracted from the reconstructed images for 99mTc-HMPAO and 99mTc-TRODAT-1 simulations, respectively, and average normalized standard deviation (NSD) based on 20 noise realizations was assessed on selected uniform 3D regions as the noise index. Visual assessment and image profiles were applied to the results of Monte Carlo simulations. The optimized design parameters of the MPH collimators were 9 pinholes with 4.7 and 2.8mm pinhole diameter, 73° acceptance angle, 127mm focal length, 167mm ROR for 12mm and 8mm target resolution, respectively. According to the optimization results, the detection efficiencies of the proposed collimators were 270 and 40% more as compared to LEHR. The Monte Carlo simulations showed that 7.9 and 6.4mm rods can be discriminated for the MPH collimators with target resolutions of 12 and 8mm, respectively. The eight 12mm-thick discs of the Defrise phantom can also be resolved clearly in the axial plane as demonstrated by the image profiles generated with the MPH collimators. The two collimator designs provide superior image quality as compared to the conventional LEHR, and shows potential to improve current brain SPECT imaging based on a conventional SPECT scanner.

GATE simulation of a new design of pinhole SPECT system for small animal brain imaging

Small animal SPECT imaging has gained an increased interest over the past decade since it is an excellent tool for developing new drugs and tracers. Therefore, there is a huge effort on the development of cost-effective SPECT detectors with high capabilities. The aim of this study is to simulate the performance characteristics of new designs for a cost effective, stationary SPECT system dedicated to small animal imaging with a focus on mice brain. The conceptual design of this SPECT system platform, Stationary Small Animal SSA-SPECT, is to use many pixelated CsI:TI detector modules with 0.4 mm × 0.4 mm pixels in order to achieve excellent intrinsic detector resolution where each module is backed by a single pinhole collimator with 0.3 mm hole diameter. In this work, we present the simulation results of four variations of the SSA-SPECT platform where the number of detector modules and FOV size is varied while keeping the detector size and collimator hole size constant. Using the NEMA NU-4 protocol, we performed spatial resolution, sensitivity, image quality simulations followed by a Derenzo-like phantom evaluation. The results suggest that all four SSA-SPECT systems can provide better than 0.063% system sensitivity and < 1.5 mm FWHM spatial resolution without resolution recovery or other correction techniques. Specifically, SSA-SPECT-1 showed a system sensitivity of 0.09% in combination with 1.1 mm FWHM spatial resolution.

Modified Gaussian Interpolations with Geometric Calibrations for Rapid Constructions of Circular-Orbit Pinhole SPECT System Matrices

For high-spatial-resolution reconstructions in SPECT imaging, an accurate and fine system matrix (H matrix) should be constructed. In this study, the modified Gaussian interpolation method (GIM) and the modified GIM combined with geometric parameter estimations (GIMGPE) are developed to rapidly construct the complete H matrices of circular-orbit 1-pinhole/4-pinhole SPECT systems. The multi-module multi-resolution (M3R) pinhole SPECT system (in the Center for Gamma-Ray Imaging, University of Arizona) is utilized to verify the usefulness of the proposed interpolation methods. The modified GIM and GIMGPE apply a geometric calibration, a rough grid-scan experiment, and point response parameterization based on Gaussian fitting to construct the full H matrices with a 0.5-mm voxel size and 60 projection angles. The geometric calibration uses a 3-point calibration phantom projected at 60 projection angles to construct the geometric projection model. In the rough grid-scan experiment, a point source is stepped on a regular grid pattern with a 1-mm grid spacing to measure the point response functions (PRFs) at the 0° projection angle. Sequentially, the measured PRFs are parameterized into two-dimensional Gaussians. The PRFs of full H matrices are interpolated by the relations between the Gaussian coefficients and the geometric parameters. The processing time of constructing the complete H matrices with the modified GIM and GIMGPE are respectively reduced by about 347× and 336×, compared to that of the full three-dimensional grid-scan experiment. According to the results of geometric calibration, the magnification of M3R with 1- and 4-pinhole patterns are 1.9× and 2.6×, respectively. The reconstructions of a hot-rod phantom show that the spatial resolutions of M3R with 1- and 4-pinhole patterns are 1.4 mm and 1.0 mm, respectively.

TH‐EF‐207A‐01: Respiratory Motion Correction in 4D‐Multi Pinhole Small Animal SPECT

Purpose:Respiration causes lesion smearing, image blurring and quality degradation, affecting lesion contrast and the ability to define correct lesion size. The spatial resolution of current multi pinhole SPECT (MPHS) scanners is sub‐millimeter. Therefore, the effect of motion is more noticeable in comparison to conventional SPECT scanner. Gated imaging aims to reduce motion artifacts. A major issue in gating is the lack of statistics and individual reconstructed frames are noisy. The increased noise in each frame, deteriorates the quantitative accuracy of the MPHS Images. The objective of this work, is to enhance the image quality in 4D‐MPHS imaging, by 4D image reconstruction.Methods:The new algorithm requires deformation vector fields (DVFs) that are calculated by non‐rigid Demons registration. The algorithm is based on the motion‐incorporated version of ordered subset expectation maximization (OSEM) algorithm. This iterative algorithm is capable to make full use of all projections to reconstruct each individual frame. To evaluate the performance of the proposed algorithm a simulation study was conducted. A fast ray tracing method was used to generate MPHS projections of a 4D digital mouse phantom with a small tumor in liver in eight different respiratory phases. To evaluate the 4D‐OSEM algorithm potential, tumor to liver activity ratio was compared with other image reconstruction methods including 3D‐MPHS and post reconstruction registered with Demons‐derived DVFs.Results:Image quality of 4D‐MPHS is greatly improved by the 4D‐OSEM algorithm. When all projections are used to reconstruct a 3D‐MPHS, motion blurring artifacts are present, leading to overestimation of the tumor size and 24% tumor contrast underestimation. This error reduced to 16% and 10% for post reconstruction registration methods and 4D‐OSEM respectively.Conclusion:4D‐OSEM method can be used for motion correction in 4D‐MPHS. The statistics and quantification are improved since all projection data are combined together to update the image.

Design and assessment of a novel SPECT system for desktop open-gantry imaging of small animals: A simulation study.

Given increasing efforts in biomedical research utilizing molecular imaging methods, development of dedicated high-performance small-animal SPECT systems has been growing rapidly in the last decade. In the present work, we propose and assess an alternative concept for SPECT imaging enabling desktop open-gantry imaging of small animals. The system, PERSPECT, consists of an imaging desk, with a set of tilted detector and pinhole collimator placed beneath it. The object to be imaged is simply placed on the desk. Monte Carlo (MC) and analytical simulations were utilized to accurately model and evaluate the proposed concept and design. Furthermore, a dedicated image reconstruction algorithm, finite-aperture-based circular projections (FABCP), was developed and validated for the system, enabling more accurate modeling of the system and higher quality reconstructed images. Image quality was quantified as a function of different tilt angles in the acquisition and number of iterations in the reconstruction algorithm. Furthermore, more complex phantoms including Derenzo, Defrise, and mouse whole body were simulated and studied. The sensitivity of the PERSPECT was 207 cps/MBq. It was quantitatively demonstrated that for a tilt angle of 30°, comparable image qualities were obtained in terms of normalized squared error, contrast, uniformity, noise, and spatial resolution measurements, the latter at ∼0.6 mm. Furthermore, quantitative analyses demonstrated that 3 iterations of FABCP image reconstruction (16 subsets/iteration) led to optimally reconstructed images. The PERSPECT, using a novel imaging protocol, can achieve comparable image quality performance in comparison with a conventional pinhole SPECT with the same configuration. The dedicated FABCP algorithm, which was developed for reconstruction of data from the PERSPECT system, can produce high quality images for small-animal imaging via accurate modeling of the system as incorporated in the forward- and back-projection steps. Meanwhile, the developed MC model and the analytical simulator of the system can be applied for further studies on development and evaluation of the system.

Enhanced PET resolution by combining pinhole collimation and coincidence detection

Spatial resolution of clinical PET scanners is limited by detector design and photon non-colinearity. Although dedicated small animal PET scanners using specialized high-resolution detectors have been developed, enhancing the spatial resolution of clinical PET scanners is of interest as a more available alternative. Multi-pinhole 511 keV SPECT is capable of high spatial resolution but requires heavily shielded collimators to avoid significant background counts. A practical approach with clinical PET detectors is to combine multi-pinhole collimation with coincidence detection. In this new hybrid modality, there are three locations associated with each event, namely those of the two detected photons and the pinhole aperture. These three locations over-determine the line of response and provide redundant information that is superior to coincidence detection or pinhole collimation alone. Multi-pinhole collimation provides high resolution and avoids non-colinearity error but is subject to collimator penetration and artifacts from overlapping projections. However the coincidence information, though at lower resolution, is valuable for determining whether the photon passed near a pinhole within the cone acceptance angle and for identifying through which pinhole the photon passed. This information allows most photons penetrating through the collimator to be rejected and avoids overlapping projections. With much improved event rejection, a collimator with minimal shielding may be used, and a lightweight add-on collimator for high resolution imaging is feasible for use with a clinical PET scanner. Monte Carlo simulations were performed of a 18F hot rods phantom and a 54-pinhole unfocused whole-body mouse collimator with a clinical PET scanner. Based on coincidence information and pinhole geometry, events were accepted or rejected, and pinhole-specific crystal-map projections were generated. Tomographic images then were reconstructed using a conventional pinhole SPECT algorithm. Hot rods of 1.4 mm diameter were resolved easily in a simulated phantom. System sensitivity was 0.09% for a simulated 70-mm line source corresponding to the NEMA NU-4 mouse phantom. Higher resolution is expected with further optimization of pinhole design, and higher sensitivity is expected with a focused and denser pinhole configuration. The simulations demonstrate high spatial resolution and feasibility of small animal imaging with an add-on multi-pinhole collimator for a clinical PET scanner. Further work is needed to develop geometric calibration and quantitative data corrections and, eventually, to construct a prototype device and produce images with physical phantoms.

Fast iterative reconstruction for helical pinhole SPECT imaging.

Pinhole SPECT for small animal has become a routine procedure in many applications of molecular biology and pharmaceutical development. There is an increasing demand in the whole body imaging of lab animals. A simple and direct solution is to scan the object along a helical trajectory, similar to a helical CT scan. The corresponding acquisition time can be greatly reduced, while the over-lapping and gap between consecutive bed positions can be avoided. However, helical pinhole SPECT inevitably leads to the tremendous increase in computational complexity when the iterative reconstruction algorithms are applied. We suggest a novel voxel-driven (VD) system model which can be integrated with geometric symmetries from helical trajectory for fast iterative image reconstruction. Such a model construction can also achieve faster calculation and lower storage requirement of the system matrix. Due to the independence among various symmetries, it permits parallel coding to further boost computation efficiency of forward/backward projection. From phantom study, the results also indicate that the proposed VD model can adequately model the helical pinhole SPECT scanner with manageable storage size of system matrix and clinically acceptable computation loading of reconstruction.

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.

Design optimization of multi-pinhole micro-SPECT configurations by signal detection tasks and system performance evaluations for mouse cardiac imaging

An optimized configuration of multi-pinhole aperture can improve the spatial resolution and the sensitivity of pinhole SPECT simultaneously. In this study, an optimization strategy of the multi-pinhole configuration with a small detector is proposed for mouse cardiac imaging. A 14 mm-diameter spherical field-of-view (FOV) is used to accommodate the mouse heart. To accelerate the optimization process, the analytic models are applied to rapidly obtain the projection areas of the FOV, the sensitivities and the spatial resolutions of numerous system designs. The candidates of optimal multi-pinhole configuration are then decided by the preliminary evaluations with the analytic models. Subsequently, the pinhole SPECT systems equipped with the designed multi-pinhole apertures are modeled in GATE to generate the imaging system matrices (H matrices) for the system performance assessments. The area under the ROC curves (AUC) of the designed systems is evaluated by signal-known-exactly/background-known-statistically detection tasks with their corresponding H matrices. In addition, the spatial resolutions are estimated by the Fourier crosstalk approach, and the sensitivities are calculated with the H matrices of designed systems, respectively. Furthermore, a series of OSEM reconstruction images of synthetic phantoms, including the hot-rod phantom, mouse heart phantom and Defrise phantom, are reconstructed with the H matrices of designed systems. To quantify the sensitivity and resolution competition in the optimization process, the AUC from the detection tasks and the resolution estimated by the Fourier crosstalk are used as the figure of merits. A trade-off function of AUC and resolution is introduced to find the optimal multi-pinhole configuration. According to the examining results, a 22.5° rotated detector plus a 4-pinhole aperture with 22.5° rotation, 20% multiplexing and 1.52X magnification is the optimized multi-pinhole configuration for the micro pinhole-SPECT applied to mouse cardiac imaging with a camera of 49 × 49 mm2 active area.

Efficient simulation of voxelized phantom in GATE with embedded SimSET multiple photon history generator

GEANT4 Application for Tomographic Emission (GATE) is a powerful Monte Carlo simulator that combines the advantages of the general-purpose GEANT4 simulation code and the specific software tool implementations dedicated to emission tomography. However, the detailed physical modelling of GEANT4 is highly computationally demanding, especially when tracking particles through voxelized phantoms. To circumvent the relatively slow simulation of voxelized phantoms in GATE, another efficient Monte Carlo code can be used to simulate photon interactions and transport inside a voxelized phantom. The simulation system for emission tomography (SimSET), a dedicated Monte Carlo code for PET/SPECT systems, is well-known for its efficiency in simulation of voxel-based objects. An efficient Monte Carlo workflow integrating GATE and SimSET for simulating pinhole SPECT has been proposed to improve voxelized phantom simulation. Although the workflow achieves a desirable increase in speed, it sacrifices the ability to simulate decaying radioactive sources such as non-pure positron emitters or multiple emission isotopes with complex decay schemes and lacks the modelling of time-dependent processes due to the inherent limitations of the SimSET photon history generator (PHG). Moreover, a large volume of disk storage is needed to store the huge temporal photon history file produced by SimSET that must be transported to GATE. In this work, we developed a multiple photon emission history generator (MPHG) based on SimSET/PHG to support a majority of the medically important positron emitters. We incorporated the new generator codes inside GATE to improve the simulation efficiency of voxelized phantoms in GATE, while eliminating the need for the temporal photon history file. The validation of this new code based on a MicroPET R4 system was conducted for 124I and 18F with mouse-like and rat-like phantoms. Comparison of GATE/MPHG with GATE/GEANT4 indicated there is a slight difference in energy spectra for energy below 50 keV due to the lack of x-ray simulation from 124I decay in the new code. The spatial resolution, scatter fraction and count rate performance are in good agreement between the two codes. For the case studies of 18F-NaF (124I-IAZG) using MOBY phantom with 1 × 1 × 1 mm3 voxel sizes, the results show that GATE/MPHG can achieve acceleration factors of approximately 3.1 × (4.5 ×), 6.5 × (10.7 ×) and 9.5 × (31.0 ×) compared with GATE using the regular navigation method, the compressed voxel method and the parameterized tracking technique, respectively. In conclusion, the implementation of MPHG in GATE allows for improved efficiency of voxelized phantom simulations and is suitable for studying clinical and preclinical imaging.

Fabrication of the pinhole aperture for AdaptiSPECT.

AdaptiSPECT is a pre-clinical pinhole SPECT imaging system under final construction at the Center for Gamma-Ray Imaging. The system is designed to be able to autonomously change its imaging configuration. The system comprises 16 detectors mounted on translational stages to move radially away and towards the center of the field-of-view. The system also possesses an adaptive pinhole aperture with multiple collimator diameters and pinhole sizes, as well as the possibility to switch between multiplexed and non-multiplexed imaging configurations. In this paper, we describe the fabrication of the AdaptiSPECT pinhole aperture and its controllers.

Rapid construction of pinhole SPECT system matrices by distance-weighted Gaussian interpolation method combined with geometric parameter estimations

In pinhole SPECT applied to small-animal studies, it is essential to have an accurate imaging system matrix, called H matrix, for high-spatial-resolution image reconstructions. Generally, an H matrix can be obtained by various methods, such as measurements, simulations or some combinations of both methods. In this study, a distance-weighted Gaussian interpolation method combined with geometric parameter estimations (DW-GIMGPE) is proposed. It utilizes a simplified grid-scan experiment on selected voxels and parameterizes the measured point response functions (PRFs) into 2D Gaussians. The PRFs of missing voxels are interpolated by the relations between the Gaussian coefficients and the geometric parameters of the imaging system with distance-weighting factors. The weighting factors are related to the projected centroids of voxels on the detector plane. A full H matrix is constructed by combining the measured and interpolated PRFs of all voxels. The PRFs estimated by DW-GIMGPE showed similar profiles as the measured PRFs. OSEM reconstructed images of a hot-rod phantom and normal rat myocardium demonstrated the effectiveness of the proposed method. The detectability of a SKE/BKE task on a synthetic spherical test object verified that the constructed H matrix provided comparable detectability to that of the H matrix acquired by a full 3D grid-scan experiment. The reduction in the acquisition time of a full 1.0-mm grid H matrix was about 15.2 and 62.2 times with the simplified grid pattern on 2.0-mm and 4.0-mm grid, respectively. A finer-grid H matrix down to 0.5-mm spacing interpolated by the proposed method would shorten the acquisition time by 8 times, additionally.

A line‐source method for aligning on‐board and other pinhole SPECT systems

In order to achieve functional and molecular imaging as patients are in position for radiation therapy, a robotic multipinhole SPECT system is being developed. Alignment of the SPECT system-to the linear accelerator (LINAC) coordinate frame and to the coordinate frames of other on-board imaging systems such as cone-beam CT (CBCT)-is essential for target localization and image reconstruction. An alignment method that utilizes line sources and one pinhole projection is proposed and investigated to achieve this goal. Potentially, this method could also be applied to the calibration of the other pinhole SPECT systems. An alignment model consisting of multiple alignment parameters was developed which maps line sources in three-dimensional (3D) space to their two-dimensional (2D) projections on the SPECT detector. In a computer-simulation study, 3D coordinates of line-sources were defined in a reference room coordinate frame, such as the LINAC coordinate frame. Corresponding 2D line-source projections were generated by computer simulation that included SPECT blurring and noise effects. The Radon transform was utilized to detect angles (α) and offsets (ρ) of the line-source projections. Alignment parameters were then estimated by a nonlinear least squares method, based on the α and ρ values and the alignment model. Alignment performance was evaluated as a function of number of line sources, Radon transform accuracy, finite line-source width, intrinsic camera resolution, Poisson noise, and acquisition geometry. Experimental evaluations were performed using a physical line-source phantom and a pinhole-collimated gamma camera attached to a robot. In computer-simulation studies, when there was no error in determining angles (α) and offsets (ρ) of the measured projections, six alignment parameters (three translational and three rotational) were estimated perfectly using three line sources. When angles (α) and offsets (ρ) were provided by the Radon transform, estimation accuracy was reduced. The estimation error was associated with rounding errors of Radon transform, finite line-source width, Poisson noise, number of line sources, intrinsic camera resolution, and detector acquisition geometry. Statistically, the estimation accuracy was significantly improved by using four line sources rather than three and by thinner line-source projections (obtained by better intrinsic detector resolution). With five line sources, median errors were 0.2 mm for the detector translations, 0.7 mm for the detector radius of rotation, and less than 0.5° for detector rotation, tilt, and twist. In experimental evaluations, average errors relative to a different, independent registration technique were about 1.8 mm for detector translations, 1.1 mm for the detector radius of rotation (ROR), 0.5° and 0.4° for detector rotation and tilt, respectively, and 1.2° for detector twist. Alignment parameters can be estimated using one pinhole projection of line sources. Alignment errors are largely associated with limited accuracy of the Radon transform in determining angles (α) and offsets (ρ) of the line-source projections. This alignment method may be important for multipinhole SPECT, where relative pinhole alignment may vary during rotation. For pinhole and multipinhole SPECT imaging on-board radiation therapy machines, the method could provide alignment of SPECT coordinates with those of CBCT and the LINAC.

An Analytical Geometric Calibration Method for Circular Cone-Beam Geometry

This work is a continuation of our previous work on geometric calibration in the circular cone-beam geometry. It is well known that seven parameters completely describe such a geometry in either flat-panel X-ray computed tomography or single pinhole SPECT imaging. Previously we developed a graphical procedure to determine the detector in-plane rotation angle independently of the other six parameters. Using the discovered geometrical relationships, in this paper we determine the remaining six parameters using the cone-beam projections of a minimum of three point objects. Our method is analytical. It makes use of the parameters of the fitted ellipse from the calibration data. The parameter estimation is accurate in the noise-free case or when there is moderate projection data truncation or shorter calibration scan range ( ≤ 360°). We perform numerical evaluations to study the robustness of the proposed method under different projection noise levels and using different data acquisition ranges. Using a full 360° scan range, the estimation accuracy and precision of our method are comparable or superior to previous methods. Using a shorter acquisition range, there may be bias in the ellipse parameters obtained by simple algebraic fitting methods. This bias will propagate to the estimated geometric parameters. Such bias can be mostly eliminated by using a more sophisticated fitting algorithm. At the same noise level, the geometric parameter estimation accuracies are comparable, but the estimation precision degrades, as the acquisition range becomes shorter.

PET-based geometrical calibration of a pinhole SPECT add-on for an animal PET scanner

We developed SPECT imaging capability on an animal PET scanner to provide a cost effective option for animal SPECT imaging. The SPECT add-on sub-system was enabled by mechanically integrating a multiple-pinhole collimator in the PET detector ring. This study introduces a method to calibrate the geometrical parameters of the SPECT add-on using the PET imaging capability of the scanner. The proposed PET imaging-based calibration method consists of two steps: (1) paint the pinhole apertures of the collimator with a positron emitting radioactive solution; and (2) image the collimator inside the scanner in PET mode. The geometrical parameters of the multi-pinhole SPECT add-on can then be derived directly from a set of PET images by simple linear calculation and used in defining the SPECT system. The method was compared to our implementation of a SPECT calibration approach with model-based fitting of SPECT projection data. The procedure for carrying out the PET imaging-based calibration method is simpler and faster than that of our implementation of the SPECT model-based calibration method. Since it does not require model fitting, the uniqueness of the calibration result is warranted. Better quality SPECT images were reconstructed using the PET-derived calibration parameters rather than our implementation of the SPECT model-based calibration parameters. We conclude that the proposed PET imaging-based calibration method provides a highly effective means for enabling SPECT imaging on a PET scanner.

Murine cardiac images obtained with focusing pinhole SPECT are barely influenced by extra-cardiac activity

Ultra-high-resolution SPECT images can be obtained with focused multipinhole collimators. Here we investigate the influence of unwanted high tracer uptake outside the scan volume on reconstructed tracer distributions inside the scan volume, for 99mTc-tetrofosmin myocardial perfusion scanning in mice. Simulated projections of a digital mouse phantom (MOBY) in a focusing multipinhole SPECT system (U-SPECT-II, MILabs, The Netherlands) were generated. With this system differently sized user-defined scan volumes can be selected, by translating the animal in 3D through the focusing collimators. Scan volume selections were set to (i) a minimal volume containing just the heart, acquired without translating the animal during scanning, (ii) a slightly larger scan volume as is typically applied for the heart, requiring only small XYZ translations during scanning, (iii) same as (ii), but extended further transaxially, and (iv) same as (ii), but extended transaxially to cover the full thorax width (gold standard). Despite an overall negative bias that is significant for the minimal scan volume, all selected volumes resulted in visually similar images. Quantitative differences in the reconstructed myocardium between gold standard and the results from the smaller scan volume selections were small; the 17 standardized myocardial segments of a bull's eye plot, normalized to the myocardial mean of the gold standard, deviated on average 6.0%, 2.5% and 1.9% for respectively the minimal, the typical and the extended scan volume, while maximum absolute deviations were respectively 18.6%, 9.0% and 5.2%. Averaged over ten low-count noisy simulations, the mean absolute deviations were respectively 7.9%, 3.2% and 1.9%. In low-count noisy simulations, the mean and maximum absolute deviations for the minimal scan volume could be reduced to respectively 4.2% and 12.5% by performing a short survey scan of the exterior activity and focusing the remaining scan time at the organ of interest. We conclude that reconstructed tracer distribution in the myocardium can be influenced by activity in surrounding organs when a too narrow scan volume is used. With slightly larger scan volumes this problem is adequately suppressed. This approach produced a smaller mean deviation and may be more effective than employing a narrow scan volume with an additional survey scan.