Spatially Variant Resolution Modelling for Iterative List-Mode PET Reconstruction

Foveated panoramic ghost imaging

Panoramic ghost imaging (PGI) is a novel method by only using a curved mirror to enlarge the field of view (FOV) of ghost imaging (GI) to 360°, making GI a breakthrough in the applications with a wide FOV. However, high-resolution PGI with high efficiency is a serious challenge because of the large amount of data. Therefore, inspired by the variant-resolution retina structure of human eye, a foveated panoramic ghost imaging (FPGI) is proposed to achieve the coexistence of a wide FOV, high resolution and high efficiency on GI by reducing the resolution redundancy, and further to promote the practical applications of GI with a wide FOV. In FPGI system, a flexible variant-resolution annular pattern structure via log-rectilinear transformation and log-polar mapping is proposed to be used for projection, which can allocate the resolution of the region of interest (ROI) and the other region of non-interest (NROI) by setting related parameters in the radial and poloidal directions independently to meet different imaging requirements. In addition, in order to reasonably reduce the resolution redundancy and avoid the loss of the necessary resolution on NROI, the variant-resolution annular pattern structure with a real fovea is further optimized to keep the ROI at any position in the center of 360° FOV by flexibly changing the initial position of the start-stop boundary on the annular pattern structure. The experimental results of the FPGI with one fovea and multiple foveae demonstrate that, compared to the traditional PGI, the proposed FPGI not only can improve the imaging quality on the ROIs with a high resolution and flexibly remain a lower-resolution imaging on the NROI with different required resolution reduction; but also reduce the reconstruction time to improve the imaging efficiency due to the reduction of the resolution redundancy.

Hybrid model of photon propagation based on the analytical and Monte Carlo methods for a dual-head PET system

The construction of photon propagation has a close relationship with the quality of reconstructed images. The classical Monte Carlo (MC) based method can model the photon propagation precisely, but it is time-consuming. The analytical method can often quickly construct a model, but its precision is a problem. How to fully exploit the advantages of the MC simulation and analytical model is an open problem. Inspired by the characteristics of the depth of interaction (DOI) detectors, which can help confirm the deposited position of a photon with DOI-encoding technology, we virtually discretize each crystal into several subcrystals to obtain the statistical distribution by MC-based simulation. Then, the statistical distribution is combined with a spatially variant solid-angle model. This combination strategy provides a hybrid model to describe photon propagation with relatively high accuracy and low computational cost. Three discretization schemes are compared to optimize the constructed photon propagation model. Several experiments are carried out to evaluate the performance of the proposed hybrid method. The metrics of full width at half maximum (FWHM), contrast recovery (CR), and coefficient of variation (COV) are adopted to quantitate the imaging results. The classical MC-based method is compared as a gold-standard reference. When a crystal is divided into two discretized positions, the convergent tendencies of CRs and COVs are consistent with that based on MC simulation method, respectively. In terms of FWHMs, the resolutions of using the MC-based model and the proposed hybrid model are 0.71 mm and 0.68 mm in the direction parallel to the detector head, respectively. This indicates the potential of the proposed method in positron emission tomography imaging.

Motion Dependent and Spatially Variant Resolution Modeling for PET Rigid Motion Correction

Recent advances in positron emission tomography (PET) have allowed to perform brain scans of freely moving animals by using rigid motion correction. One of the current challenges in these scans is that, due to the PET scanner spatially variant point spread function (SVPSF), motion corrected images have a motion dependent blurring since animals can move throughout the entire field of view (FOV). We developed a method to calculate the image-based resolution kernels of the motion dependent and spatially variant PSF (MD-SVPSF) to correct the loss of spatial resolution in motion corrected reconstructions. The resolution kernels are calculated for each voxel by sampling and averaging the SVPSF at all positions in the scanner FOV where the moving object was measured. In resolution phantom scans, the use of the MD-SVPSF resolution model improved the spatial resolution in motion corrected reconstructions and corrected the image deformation caused by the parallax effect consistently for all motion patterns, outperforming the use of a motion independent SVPSF or Gaussian kernels. Compared to motion correction in which the SVPSF is applied independently for every pose, our method performed similarly, but with more than two orders of magnitude faster computation time. Importantly, in scans of freely moving mice, brain regional quantification in motion-free and motion corrected images was better correlated when using the MD-SVPSF in comparison with motion independent SVPSF and a Gaussian kernel. The method developed here allows to obtain consistent spatial resolution and quantification in motion corrected images, independently of the motion pattern of the subject.

Open-field PET: Simultaneous brain functional imaging and behavioural response measurements in freely moving small animals

A comprehensive understanding of how the brain responds to a changing environment requires techniques capable of recording functional outputs at the whole-brain level in response to external stimuli. Positron emission tomography (PET) is an exquisitely sensitive technique for imaging brain function but the need for anaesthesia to avoid motion artefacts precludes concurrent behavioural response studies. Here, we report a technique that combines motion-compensated PET with a robotically-controlled animal enclosure to enable simultaneous brain imaging and behavioural recordings in unrestrained small animals. The technique was used to measure in vivo displacement of [¹¹C]raclopride from dopamine D2 receptors (D2R) concurrently with changes in the behaviour of awake, freely moving rats following administration of unlabelled raclopride or amphetamine. The timing and magnitude of [¹¹C]raclopride displacement from D2R were reliably estimated and, in the case of amphetamine, these changes coincided with a marked increase in stereotyped behaviours and hyper-locomotion. The technique, therefore, allows simultaneous measurement of changes in brain function and behavioural responses to external stimuli in conscious unrestrained animals, giving rise to important applications in behavioural neuroscience.

Sinogram Blurring Matrix Estimation From Point Sources Measurements With Rank-One Approximation for Fully 3-D PET

An accurate system matrix is essential in positron emission tomography (PET) for reconstructing high quality images. To reduce storage size and image reconstruction time, we factor the system matrix into a product of a geometry projection matrix and a sinogram blurring matrix. The geometric projection matrix is computed analytically and the sinogram blurring matrix is estimated from point source measurements. Previously, we have estimated a 2-D blurring matrix for a preclinical PET scanner. The 2-D blurring matrix only considers blurring effects within a transaxial sinogram and does not compensate for inter-sinogram blurring effects. For PET scanners with a long axial field of view, inter-sinogram blurring can be a major problem influencing the image quality in the axial direction. Hence, the estimation of a 4-D blurring matrix is desirable to further improve the image quality. The 4-D blurring matrix estimation is an ill-conditioned problem due to the large number of unknowns. Here, we propose a rank-one approximation for each blurring kernel image formed by a row vector of the sinogram blurring matrix to improve the stability of the 4-D blurring matrix estimation. The proposed method is applied to the simulated data as well as the real data obtained from an Inveon microPET scanner. The results show that the newly estimated 4-D blurring matrix can improve the image quality over those obtained with a 2-D blurring matrix and requires less point source scans to achieve similar image quality compared with an unconstrained 4-D blurring matrix estimation.