Deep learning-based positron emission tomography (PET) image denoising offers the potential to reduce radiation exposure and scanning time by transforming low-count images into high-count equivalents. However, existing methods typically blur crucial details, leading to inaccurate lesion quantification. This paper proposes a lesion-perceived and quantification-consistent modulation (LeqMod) strategy for enhanced PET image denoising, via employing downstream lesion quantification analysis as auxiliary tools. The LeqMod is a plug-and-play design adaptable to a wide range of model architectures, modulating the sampling and optimization procedures of model training without adding any computational burden to the inference phase. Specifically, the LeqMod consists of two components, the lesion-perceived modulation (LeMod) and the multiscale quantification-consistent modulation (QuMod). The LeMod enhances lesion contrast and visibility by allocating higher sampling weights and stricter loss criteria to lesion-present samples determined by an auxiliary segmentation network than lesion-absent ones. The QuMod further emphasizes quantification accuracy for both the mean and maximum standardized uptake value ( ${\mathrm {SUV}}_{{\textit {mean}}}$ and ${\mathrm {SUV}}_{{\textit {max}}}$ ) across multiscale sub-regions throughout the entire image, thereby reducing biases of denoised results relative to high-count references. Experiments conducted on large PET datasets from multiple centers and vendors, and varying noise levels demonstrated the LeqMod efficacy across various denoising frameworks. Compared to frameworks without LeqMod, the integration of LeqMod reduces the lesion ${\mathrm {SUV}}_{{\textit {max}}}$ bias by 5.92% on average and increases the peak signal-to-noise ratio (PSNR) by 0.36 on average, when denoising images across participating sites. (Code is available at https://github.com/mhxiaaa/LeqMod_PET_denoising).

Quantitative PET attenuation correction (AC) for cardiac PET/CT and PET/MR is a challenging problem. We propose and evaluate an AC approach that uses coincidences from a relatively weak and physically fixed sparse external source, in combination with that from the patient, to reconstruct μ -maps based on physics principles alone. The low 30 cm3 volume of the source makes it easy to fill and place, and the method does not use prior image data or attenuation map assumptions. Our supplemental transmission aided maximum likelihood reconstruction of attenuation and activity (sTX-MLAA) algorithm contains an attenuation map update that maximizes the likelihood of terms representing coincidences originating from tracer in the patient and a weighted expression of counts segmented from the external source alone. Both external source and patient scatter and randoms are fully corrected. We evaluated performance of sTX-MLAA compared to reference standard CT-based AC with FDG PET/CT phantom studies; including modeling a patient with myocardial inflammation. Through an ROI analysis we measured ≤ 5 % bias in activity concentrations for PET images generated with sTX-MLAA and a TX source strength ≥ 12.7 MBq, relative to CT-AC. PET background variability (from noise and sparse sampling) was substantially reduced with sTX-MLAA compared to using counts segmented from the transmission source alone for AC. Results suggest that sTX-MLAA will enable quantitative PET during cardiac PET/CT and PET/MR of human patients.

Image reconstruction of low-count positron emission tomography (PET) data is challenging. Kernel methods address the challenge by incorporating image prior information in the forward model of iterative PET image reconstruction. The kernelized expectation-maximization (KEM) algorithm has been developed and demonstrated to be effective and easy to implement. A common approach for a further improvement of the kernel method would be adding an explicit regularization, which however leads to a complex optimization problem. In this paper, we propose an implicit regularization for the kernel method by using a deep coefficient prior, which represents the kernel coefficient image in the PET forward model using a convolutional neural-network. To solve the maximum-likelihood neural network-based reconstruction problem, we apply the principle of optimization transfer to derive a neural KEM algorithm. Each iteration of the algorithm consists of two separate steps: a KEM step for image update from the projection data and a deep-learning step in the image domain for updating the kernel coefficient image using the neural network. This optimization algorithm is guaranteed to monotonically increase the data likelihood. The results from computer simulations and real patient data have demonstrated that the neural KEM can outperform existing KEM and deep image prior methods.

We present a method that quantifies point-wise changes in surface morphology of the bones of the human wrist. The proposed method, referred to as Registration-based Bone Morphometry (RBM), consists of two steps: an atlas selection step and an atlas warping step. The atlas for individual wrist bones was selected based on the shortest ℓ2 distance to the ensemble of wrist bones from a database of a healthy population of subjects. The selected atlas was then warped to the corresponding bones of individuals in the population using a non-linear registration method based on regularized ℓ2 distance minimization. The displacement field thus calculated showed local differences in bone shape that then were used for the analysis of group differences. Our results indicate that RBM has potential to provide a standardized approach to shape analysis of bones of the human wrist. We demonstrate the performance of RBM for examining group differences in wrist bone shapes based on sex and between those of the right and left wrists in healthy individuals. We also present data to show the application of RBM for tracking bone erosion status in rheumatoid arthritis.