High Temporal-Resolution Dynamic PET Image Reconstruction Using a New Spatiotemporal Kernel Method

Quantitative PET imaging and modeling of molecular blood-brain barrier permeability

Neuroimaging of blood-brain barrier permeability has been instrumental in identifying its broad involvement in neurological and systemic diseases. However, current methods evaluate the blood-brain barrier mainly as a structural barrier. Here we developed a non-invasive positron emission tomography method in humans to measure the blood-brain barrier permeability of molecular radiotracers that cross the blood-brain barrier through its molecule-specific transport mechanism. Our method uses high-temporal resolution dynamic imaging and kinetic modeling for multiparametric imaging and quantification of the blood-brain barrier permeability-surface area product of molecular radiotracers. We show, in humans, our method can resolve blood-brain barrier permeability across three radiotracers and demonstrate its utility in studying brain aging and brain-body interactions in metabolic dysfunction-associated steatotic liver inflammation. Our method opens new directions to effectively study the molecular permeability of the human blood-brain barrier in vivo using the large catalogue of available molecular positron emission tomography tracers.

Quo Vadis Hyperpolarized 13C MRI?

Over the last two decades, hyperpolarized 13C MRI has gained significance in both preclinical and clinical studies, hereby relying on technologies like PHIP-SAH (ParaHydrogen-Induced Polarization-Side Arm Hydrogenation), SABRE (Signal Amplification by Reversible Exchange), and dDNP (dissolution Dynamic Nuclear Polarization), with dDNP being applied in humans. A clinical dDNP polarizer has enabled studies across 24 sites, despite challenges like high cost and slow polarization. Parahydrogen-based techniques like SABRE and PHIP offer faster, more cost-efficient alternatives but require molecule-specific optimization. The focus has been on imaging metabolism of hyperpolarized probes, which requires long T1, high polarization and rapid contrast generation. Efforts to establish novel probes, improve acquisition techniques and enhance data analysis methods including artificial intelligence are ongoing. Potential clinical value of hyperpolarized 13C MRI was demonstrated primarily for treatment response assessment in oncology, but also in cardiology, nephrology, hepatology and CNS characterization. In this review on biomedical hyperpolarized 13C MRI, we summarize important and recent advances in polarization techniques, probe development, acquisition and analysis methods as well as clinical trials. Starting from those we try to sketch a trajectory where the field of biomedical hyperpolarized 13C MRI might go.

Mucociliary clearance is impaired in small airways of cystic fibrosis pigs

Cystic fibrosis (CF) is a genetic disorder characterized by recurrent airway infections, inflammation, impaired mucociliary clearance, and progressive decline in lung function. The disease may start in the small airways; however, this is difficult to prove due to the limited accessibility of the small airways with the current single-photon mucociliary clearance assay. Here, we developed a dynamic positron emission tomography assay with high spatial and temporal resolution. We tested that mucociliary clearance is abnormal in the small airways of newborn cystic fibrosis pigs. Clearance of [68Ga]-tagged macroaggregated albumin from small airways started immediately after delivery and continued for the duration of the study. Initial clearance was fast but slowed down a few minutes after delivery. Cystic fibrosis pigs' small airways cleared significantly less than non-CF pigs' small airways (non-CF 25.1 ± 3.1% vs. CF 14.6 ± 0.1%). Stimulation of the cystic fibrosis airways with the purinergic secretagogue uridine-5'-triphosphate (UTP) further impaired clearance (non-CF with UTP 20.9 ± 0.3% vs. CF with UTP 13.0 ± 1.8%). None of the cystic fibrosis pigs treated with UTP (n = 6) cleared more than 20% of the delivered dose. These data indicate that mucociliary clearance in the small airways is fast and can easily be missed if the assay is not sensitive enough. The data also indicate that mucociliary clearance is impaired in the small airways of cystic fibrosis pigs. This defect is exacerbated by stimulation of mucus secretions with purinergic agonists.NEW & NOTEWORTHY We developed a novel positron emission tomography scan assay with unprecedented temporal and spatial resolution to measure mucociliary clearance in the small airways. We proved a long-standing but unproven assertion that mucociliary clearance is inherently abnormal in the small airways of newborn cystic fibrosis piglets that are otherwise free of infection or inflammation. This technique can be easily extended to other airway diseases such as asthma, idiopathic pulmonary fibrosis, or chronic obstructive pulmonary disease.

Quantitative PET imaging and modeling of molecular blood-brain barrier permeability

Blood-brain barrier (BBB) disruption is involved in the pathogenesis and progression of many neurological and systemic diseases. Non-invasive assessment of BBB permeability in humans has mainly been performed with dynamic contrast-enhanced magnetic resonance imaging, evaluating the BBB as a structural barrier. Here, we developed a novel non-invasive positron emission tomography (PET) method in humans to measure the BBB permeability of molecular radiotracers that cross the BBB through different transport mechanisms. Our method uses high-temporal resolution dynamic imaging and kinetic modeling to jointly estimate cerebral blood flow and tracer-specific BBB transport rate from a single dynamic PET scan and measure the molecular permeability-surface area (PS) product of the radiotracer. We show our method can resolve BBB PS across three PET radiotracers with greatly differing permeabilities, measure reductions in BBB PS of 18F-fluorodeoxyglucose (FDG) in healthy aging, and demonstrate a possible brain-body association between decreased FDG BBB PS in patients with metabolic dysfunction-associated steatotic liver inflammation. Our method opens new directions to efficiently study the molecular permeability of the human BBB in vivo using the large catalogue of available molecular PET tracers.

Optimization-derived blood input function using a kernel method and its evaluation with total-body PET for brain parametric imaging

Dynamic PET allows quantification of physiological parameters through tracer kinetic modeling. For dynamic imaging of brain or head and neck cancer on conventional PET scanners with a short axial field of view, the image-derived input function (ID-IF) from intracranial blood vessels such as the carotid artery (CA) suffers from severe partial volume effects. Alternatively, optimization-derived input function (OD-IF) by the simultaneous estimation (SIME) method does not rely on an ID-IF but derives the input function directly from the data. However, the optimization problem is often highly ill-posed. We proposed a new method that combines the ideas of OD-IF and ID-IF together through a kernel framework. While evaluation of such a method is challenging in human subjects, we used the uEXPLORER total-body PET system that covers major blood pools to provide a reference for validation.

METHODS

The conventional SIME approach estimates an input function using a joint estimation together with kinetic parameters by fitting time activity curves from multiple regions of interests (ROIs). The input function is commonly parameterized with a highly nonlinear model which is difficult to estimate. The proposed kernel SIME method exploits the CA ID-IF as a priori information via a kernel representation to stabilize the SIME approach. The unknown parameters are linear and thus easier to estimate. The proposed method was evaluated using 18F-fluorodeoxyglucose studies with both computer simulations and 20 human-subject scans acquired on the uEXPLORER scanner. The effect of the number of ROIs on kernel SIME was also explored.

RESULTS

The estimated OD-IF by kernel SIME showed a good match with the reference input function and provided more accurate estimation of kinetic parameters for both simulation and human-subject data. The kernel SIME led to the highest correlation coefficient (R = 0.97) and the lowest mean absolute error (MAE = 10.5 %) compared to using the CA ID-IF (R = 0.86, MAE = 108.2 %) and conventional SIME (R = 0.57, MAE = 78.7 %) in the human-subject evaluation. Adding more ROIs improved the overall performance of the kernel SIME method.

CONCLUSION

The proposed kernel SIME method shows promise to provide an accurate estimation of the blood input function and kinetic parameters for brain PET parametric imaging.

Mucociliary Clearance is Impaired in Small Airways of Cystic Fibrosis Pigs

Rationale

Cystic fibrosis is a genetic disorder characterized by recurrent airway infections, inflammation, and progressive decline in lung function. Autopsy and spirometry data suggest that cystic fibrosis may start in the small airways which, due to the fractal nature of the airways, account for most of the airway tree surface area. However, they are not easily accessible for testing.

Objectives

Here, we tested the hypothesis that mucociliary clearance is abnormal in the small airways of newborn cystic fibrosis pigs.

Methods

Current mucociliary clearance assays are limited therefore we developed a dynamic positron emission tomography scan assay with high spatial and temporal resolution. Each study was accompanied by a high-resolution computed tomography scan that helped identify the thin outer region of the lung that contained small airways.

Measurements and Main Results

Clearance of aerosolized [ 68 Ga]macro aggregated albumin from distal airways occurred within minutes after delivery and followed a two-phase process. In cystic fibrosis pigs, both early and late clearance rates were slower. Stimulation of the cystic fibrosis airways with the purinergic agonist UTP further impaired late clearance. Only 1 cystic fibrosis pig treated with UTP out of 6 cleared more than 20% of the delivered dose.

Conclusions

These data indicate that mucociliary transport in the small airways is fast and can easily be missed if the acquisition is not fast enough. The data also indicate that mucociliary transport is impaired in small airways of cystic fibrosis pigs. This defect is exacerbated by stimulation of mucus secretions with purinergic agonists.

Using imaging modalities to predict nanoparticle distribution and treatment efficacy in solid tumors: The growing role of ultrasound

Nanomedicine in oncology has not had the success in clinical impact that was anticipated in the early stages of the field's development. Ideally, nanomedicines selectively accumulate in tumor tissue and reduce systemic side effects compared to traditional chemotherapeutics. However, this has been more successful in preclinical animal models than in humans. The causes of this failure to translate may be related to the intra- and inter-patient heterogeneity of the tumor microenvironment. Predicting whether a patient will respond positively to treatment prior to its initiation, through evaluation of characteristics like nanoparticle extravasation and retention potential in the tumor, may be a way to improve nanomedicine success rate. While there are many potential strategies to accomplish this, prediction and patient stratification via noninvasive medical imaging may be the most efficient and specific strategy. There have been some preclinical and clinical advances in this area using MRI, CT, PET, and other modalities. An alternative approach that has not been studied as extensively is biomedical ultrasound, including techniques such as multiparametric contrast-enhanced ultrasound (mpCEUS), doppler, elastography, and super-resolution processing. Ultrasound is safe, inexpensive, noninvasive, and capable of imaging the entire tumor with high temporal and spatial resolution. In this work, we summarize the in vivo imaging tools that have been used to predict nanoparticle distribution and treatment efficacy in oncology. We emphasize ultrasound imaging and the recent developments in the field concerning CEUS. The successful implementation of an imaging strategy for prediction of nanoparticle accumulation in tumors could lead to increased clinical translation of nanomedicines, and subsequently, improved patient outcomes. This article is categorized under: Diagnostic Tools In Vivo Nanodiagnostics and Imaging Therapeutic Approaches and Drug Discovery Nanomedicine for Oncologic Disease Therapeutic Approaches and Drug Discovery Emerging Technologies.

An update on the use of image-derived input functions for human PET studies: new hopes or old illusions?

BACKGROUND

The need for arterial blood data in quantitative PET research limits the wider usability of this imaging method in clinical research settings. Image-derived input function (IDIF) approaches have been proposed as a cost-effective and non-invasive alternative to gold-standard arterial sampling. However, this approach comes with its own limitations-partial volume effects and radiometabolite correction among the most important-and varying rates of success, and the use of IDIF for brain PET has been particularly troublesome.

MAIN BODY

This paper summarizes the limitations of IDIF methods for quantitative PET imaging and discusses some of the advances that may make IDIF extraction more reliable. The introduction of automated pipelines (both commercial and open-source) for clinical PET scanners is discussed as a way to improve the reliability of IDIF approaches and their utility for quantitative purposes. Survey data gathered from the PET community are then presented to understand whether the field's opinion of the usefulness and validity of IDIF is improving. Finally, as the introduction of next-generation PET scanners with long axial fields of view, ultra-high sensitivity, and improved spatial and temporal resolution, has also brought IDIF methods back into the spotlight, a discussion of the possibilities offered by these state-of-the-art scanners-inclusion of large vessels, less partial volume in small vessels, better description of the full IDIF kinetics, whole-body modeling of radiometabolite production-is included, providing a pathway for future use of IDIF.

CONCLUSION

Improvements in PET scanner technology and software for automated IDIF extraction may allow to solve some of the major limitations associated with IDIF, such as partial volume effects and poor temporal sampling, with the exciting potential for accurate estimation of single kinetic rates. Nevertheless, until individualized radiometabolite correction can be performed effectively, IDIF approaches remain confined at best to a few tracers.

Total-Body Positron Emission Tomography: Adding New Perspectives to Cardiovascular Research

The recent advent of positron emission tomography (PET) scanners that can image the entire human body opens up intriguing possibilities for cardiovascular research and future clinical applications. These new systems permit radiotracer kinetics to be measured in all organs simultaneously. They are particularly well suited to study cardiovascular disease and its effects on the entire body. They could also play a role in quantitatively measuring physiologic, metabolic, and immunologic responses in healthy individuals to a variety of stressors and lifestyle interventions, and may ultimately be instrumental for evaluating novel therapeutic agents and their molecular effects across different tissues. In this review, we summarize recent progress in PET technology and methodology, discuss several emerging cardiovascular applications for total-body PET, and place this in the context of multiorgan and systems medicine. Finally, we discuss opportunities that will be enabled by the technology, while also pointing to some of the challenges that still need to be addressed.

A depth-of-interaction rebinning method based on both geometric and activity weights

BACKGROUND AND OBJECTIVE

For positron emission tomography (PET) scanners with depth-of-interaction (DOI) measurement, the DOI rebinning method that utilizes DOI information to process the projection data is critical to image quality. Current DOI rebinning methods map coincidence events onto the rebinned sinogram based on the correlation of lines of response (LOR). This study aims to incorporate prior radioactivity distribution of the imaging object into DOI rebinning to obtain better image quality.

METHODS

A DOI rebinning method based on both geometric and activity weights was proposed to assign coincidence events to the rebinned sinogram defined by a virtual ring. The geometric weights, representing the correlation between LORs, were calculated based on the areas of intersection. The activity weights, reflecting the activity distribution of the imaging object, were derived from the previous reconstructed image.

RESULTS

Monte Carlo simulation data from four phantoms, including the image quality phantom, Derenzo phantom, and two rat-like ROBY phantoms, was used to evaluate the proposed method. The recovery coefficient (RC), contrast recovery coefficient (CRC), structural similarity index measure (SSIM), and peak signal-to-noise ratio (PSNR) were used as image quality metrics. Compared to other DOI rebinning methods, the proposed method achieved the highest RC (maximum improvement of 32%) and CRC at the same noise level and was also optimal in terms of the SSIM and PSNR. Meanwhile, incorporating the prior activity distribution into DOI rebinning also improved the image reconstruction speed.

CONCLUSIONS

This work developed a new DOI rebinning method combining the correlation of LORs with the prior activity distribution, achieving relatively optimal image quality and reconstruction speed. Furthermore, it still needs to be evaluated on the actual equipment.

Neural KEM: A Kernel Method With Deep Coefficient Prior for PET Image Reconstruction

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.

HYPR4D kernel method on TOF PET data with validations including image-derived input function

Background

Positron emission tomography (PET) images are typically noisy especially in dynamic imaging where the PET data are divided into a number of short temporal frames often with a low number of counts. As a result, image features such as contrast and time–activity curves are highly variable. Noise reduction in PET is thus essential. Typical noise reduction methods tend to not preserve image features/patterns (e.g. contrast and size dependent) accurately. In this work, we report the first application of our HYPR4D kernel method on time-of-flight (TOF) PET data (i.e. PSF-HYPR4D-K-TOFOSEM). The proposed HYPR4D kernel method makes use of the mean 4D high frequency features and inconsistent noise patterns over OSEM subsets as well as the low noise property of the early reconstruction updates to achieve prior-free de-noising. The method was implemented and tested on the GE SIGNA PET/MR and was compared to the TOF reconstructions with PSF resolution modeling available on the system, namely PSF-TOFOSEM with and without standard post filter and PSF-TOFBSREM (TOF Q.Clear) with various beta values (regularization strengths).

Results

Results from experimental contrast phantom and human subject data with various PET tracers showed that the proposed method provides more robust and accurate image features compared to other regularization methods. The preservation of contrast for the PSF-HYPR4D-K-TOFOSEM was observed to be better and less dependent on the contrast and size of the target structures as compared to TOF Q.Clear and PSF-TOFOSEM with filter. At the same contrast level, PSF-HYPR4D-K-TOFOSEM achieved better 4D noise suppression than other methods (e.g. >2 times lower noise than TOF Q.Clear at the highest contrast). We also present a novel voxel search method to obtain an image-derived input function (IDIF) and demonstrate that the obtained IDIF is the most quantitative w.r.t. the measured blood samples when the acquired data are reconstructed with PSF-HYPR4D-K-TOFOSEM.

Conclusions

The overall results support superior performance of the PSF-HYPR4D-K-TOFOSEM for TOF PET data and demonstrate that the proposed method is likely suitable for all imaging tasks including the generation of IDIF without requiring any prior information as well as further improving the effective sensitivity of the imaging system.

Deep Kernel Representation for Image Reconstruction in PET

Image reconstruction for positron emission tomography (PET) is challenging because of the ill-conditioned tomographic problem and low counting statistics. Kernel methods address this challenge by using kernel representation to incorporate image prior information in the forward model of iterative PET image reconstruction. Existing kernel methods construct the kernels commonly using an empirical process, which may lead to unsatisfactory performance. In this paper, we describe the equivalence between the kernel representation and a trainable neural network model. A deep kernel method is then proposed by exploiting a deep neural network to enable automated learning of an improved kernel model and is directly applicable to single subjects in dynamic PET. The training process utilizes available image prior data to form a set of robust kernels in an optimized way rather than empirically. The results from computer simulations and a real patient dataset demonstrate that the proposed deep kernel method can outperform the existing kernel method and neural network method for dynamic PET image reconstruction.

Spatiotemporal Kernel Reconstruction for Linear Parametric Neurotransmitter PET Kinetic Modeling in Motion Correction Brain PET of Awake Rats

The linear parametric neurotransmitter positron emission tomography (lp-ntPET) kinetic model can be used to detect transient changes (activation) in endogenous neurotransmitter levels. Preclinical PET scans in awake animals can be performed to investigate neurotransmitter transient changes. Here we use the spatiotemporal kernel reconstruction (Kernel) for noise reduction in dynamic PET, and lp-ntPET kinetic modeling. Kernel is adapted for motion correction reconstruction, applied in awake rat PET scans. We performed 2D rat brain phantom simulation using the ntPET model at 3 different noise levels. Data was reconstructed with independent frame reconstruction (IFR), IFR with HYPR denoising, and Kernel, and lp-ntPET kinetic parameters (k 2a : efflux rate, γ: activation magnitude, t d : activation onset time, and t p : activation peak time) were calculated. Additionally, significant activation magnitude (γ) difference with respect to a region with no activation (rest) was calculated. Finally, [11C]raclopride experiments were performed in anesthetized and awake rats, injecting cold raclopride at 20 min after scan start to simulate endogenous neurotransmitter release. For simulated data at the regional level, IFR coefficient of variation (COV) of k 2a , γ, t d and t p was reduced with HYPR denoising, but Kernel showed the lowest COV (2 fold reduction compared with IFR). At the pixel level the same trend is observed for k 2a , γ, t d and t p COV, but reduction is larger with Kernel compared with IFR (10-14 fold). Bias in γ with respect with noise-free values was additionally reduced using Kernel (difference of 292, 72.4, and -6.92% for IFR, IFR+KYPR, and Kernel, respectively). Significant difference in activation between the rest and active region could be detected at a simulated activation of 160% for IFR and IFR+HYPR, and of 120% for Kernel. In rat experiments, lp-ntPET parameters have better confidence intervals using Kernel. In the γ, and t d parametric maps, the striatum structure can be identified with Kernel but not with IFR. Striatum voxel-wise γ, t d and t p values have lower variability using Kernel compared with IFR and IFR+HYPR. The spatiotemporal kernel reconstruction adapted for motion correction reconstruction allows to improve lp-ntPET kinetic modeling noise in awake rat studies, as well as detection of subtle neurotransmitter activations.

Hybrid kernelised expectation maximisation for Bremsstrahlung SPECT reconstruction in SIRT with 90Y micro-spheres

Background

Selective internal radiation therapy with Yttrium-90 microspheres is an effective therapy for liver cancer and liver metastases. Yttrium-90 is mainly a high-energy beta particle emitter. These beta particles emit Bremsstrahlung radiation during their interaction with tissue making post-therapy imaging of the radioactivity distribution feasible. Nevertheless, image quality and quantification is difficult due to the continuous energy spectrum which makes resolution modelling, attenuation and scatter estimation challenging and therefore the dosimetry quantification is inaccurate. As a consequence a reconstruction algorithm able to improve resolution could be beneficial.

Methods

In this study, the hybrid kernelised expectation maximisation (HKEM) is used to improve resolution and contrast and reduce noise, in addition a modified HKEM called frozen HKEM (FHKEM) is investigated to further reduce noise. The iterative part of the FHKEM kernel was frozen at the 72nd sub-iteration. When using ordered subsets algorithms the data is divided in smaller subsets and the smallest algorithm iterative step is called sub-iteration. A NEMA phantom with spherical inserts was used for the optimisation and validation of the algorithm, and data from 5 patients treated with Selective internal radiation therapy were used as proof of clinical relevance of the method.

Results

The results suggest a maximum improvement of 56% for region of interest mean recovery coefficient at fixed coefficient of variation and better identification of the hot volumes in the NEMA phantom. Similar improvements were achieved with patient data, showing 47% mean value improvement over the gold standard used in hospitals.

Conclusions

Such quantitative improvements could facilitate improved dosimetry calculations with SPECT when treating patients with Selective internal radiation therapy, as well as provide a more visible position of the cancerous lesions in the liver.

Dynamic PET Imaging Using Dual Texture Features

Purpose: This study aims to explore the impact of adding texture features in dynamic positron emission tomography (PET) reconstruction of imaging results. Methods: We have improved a reconstruction method that combines radiological dual texture features. In this method, multiple short time frames are added to obtain composite frames, and the image reconstructed by composite frames is used as the prior image. We extract texture features from prior images by using the gray level-gradient cooccurrence matrix (GGCM) and gray-level run length matrix (GLRLM). The prior information contains the intensity of the prior image, the inverse difference moment of the GGCM and the long-run low gray-level emphasis of the GLRLM. Results: The computer simulation results show that, compared with the traditional maximum likelihood, the proposed method obtains a higher signal-to-noise ratio (SNR) in the image obtained by dynamic PET reconstruction. Compared with similar methods, the proposed algorithm has a better normalized mean squared error (NMSE) and contrast recovery coefficient (CRC) at the tumor in the reconstructed image. Simulation studies on clinical patient images show that this method is also more accurate for reconstructing high-uptake lesions. Conclusion: By adding texture features to dynamic PET reconstruction, the reconstructed images are more accurate at the tumor.

Dynamic PET image reconstruction incorporating a median nonlocal means kernel method

In dynamic positron emission tomography (PET) imaging, the reconstructed image of a single frame often exhibits high noise due to limited counting statistics of projection data. This study proposed a median nonlocal means (MNLM)-based kernel method for dynamic PET image reconstruction. The kernel matrix is derived from median nonlocal means of pre-reconstructed composite images. Then the PET image intensities in all voxels were modeled as a kernel matrix multiplied by coefficients and incorporated into the forward model of PET projection data. Then, the coefficients of each feature were estimated by the maximum likelihood method. Using simulated low-count dynamic data of Zubal head phantom, the quantitative performance of the proposed MNLM kernel method was investigated and compared with the maximum-likelihood method, conventional kernel method with and without median filter, and nonlocal means (NLM) kernel method. Simulation results showed that the MNLM kernel method achieved visual and quantitative accuracy improvements (in terms of the ensemble mean squared error, bias versus variance, and contrast versus noise performances). Especially for frame 2 with the lowest count level of a single frame, the MNLM kernel method achieves lower ensemble mean squared error (10.43%) than the NLM kernel method (13.68%), conventional kernel method with and without median filter (11.88% and 23.50%), and MLEM algorithm (24.77%). The study on real low-dose ¹⁸F-FDG rat data also showed that the MNLM kernel method outperformed other methods in visual and quantitative accuracy improvements (in terms of regional noise versus intensity mean performance).

Multiparametric cardiac18F-FDG PET in humans: pilot comparison of FDG delivery rate with82Rb myocardial blood flow

Myocardial blood flow (MBF) and flow reserve are usually quantified in the clinic with positron emission tomography (PET) using a perfusion-specific radiotracer (e.g.⁸²Rb-chloride). However, the clinical accessibility of existing perfusion tracers remains limited. Meanwhile,¹⁸F-fluorodeoxyglucose (FDG) is a commonly used radiotracer for PET metabolic imaging without similar limitations. In this paper, we explore the potential of¹⁸F-FDG for myocardial perfusion imaging by comparing the myocardial FDG delivery rateK₁with MBF as determined by dynamic⁸²Rb PET in fourteen human subjects with heart disease. Two sets of FDGK₁were derived from one-hour dynamic FDG scans. One was the original FDGK₁estimates and the other was the correspondingK₁values that were linearly normalized for blood glucose levels. A generalized Renkin-Crone model was used to fit FDGK₁with Rb MBF, which then allowed for a nonlinear extraction fraction correction for converting FDGK₁to MBF. The linear correlation between FDG-derived MBF and Rb MBF was moderate (r= 0.79) before the glucose normalization and became much improved (r> 0.9) after glucose normalization. The extraction fraction of FDG was also similar to that of Rb-chloride in the myocardium. The results from this pilot study suggest that dynamic cardiac FDG-PET with tracer kinetic modeling has the potential to provide MBF in addition to its conventional use for metabolic imaging.

MRI-aided kernel PET image reconstruction method based on texture features

We investigate the reconstruction of low-count positron emission tomography (PET) projection, which is an important, but challenging, task. Using the texture feature extraction method of radiomics, i.e. the gray-level co-occurrence matrix (GLCM), texture features can be extracted from magnetic resonance imaging images with high-spatial resolution. In this work, we propose a kernel reconstruction method combining autocorrelation texture features derived from the GLCM. The new kernel function includes the correlations of both the intensity and texture features from the prior image. By regarding the GLCM as a discrete approximation of a probability density function, the asymptotically gray-level-invariant autocorrelation texture feature is generated, which can maintain the accuracy of texture features extracted from small image regions by reducing the number of quantized image gray levels. A computer simulation shows that the proposed method can effectively reduce the noise in the reconstructed image compared to the maximum likelihood expectation maximum method and improve the image quality and tumor region accuracy compared to the original kernel method for low-count PET reconstruction. A simulation study on clinical patient images also shows that the proposed method can improve the whole image quality and that the reconstruction of a high-uptake lesion is more accurate than that achieved by the original kernel method.

Technical note: A preliminary study of dual-tracer PET image reconstruction guided by FDG and/or MR kernels

PURPOSE

Clinically, single radiotracer positron emission tomography (PET) imaging is a commonly used examination method; however, since each radioactive tracer reflects the information of only one kind of cell, it easily causes false negatives or false positives in disease diagnosis. Therefore, reasonably combining two or more radiotracers is recommended to improve the accuracy of diagnosis and the sensitivity and specificity of the disease when conditions permit.

METHODS

This paper proposes incorporating ¹⁸ F-fluorodeoxyglucose (FDG) as a higher-quality PET image to guide the reconstruction of other lower-count ¹¹ C-methionine (MET) PET datasets to compensate for the lower image quality by a popular kernel algorithm. Specifically, the FDG prior is needed to extract kernel features, and these features were used to build a kernel matrix using a k-nearest-neighbor (kNN) search for MET image reconstruction. We created a 2-D brain phantom to validate the proposed method by simulating sinogram data containing Poisson random noise and quantitatively compared the performance of the proposed FDG-guided kernelized expectation maximization (KEM) method with the performance of Gaussian and non-local means (NLM) smoothed maximum likelihood expectation maximization (MLEM), MR-guided KEM, and multi-guided-S KEM algorithms. Mismatch experiments between FDG/MR and MET data were also carried out to investigate the outcomes of possible clinical situations.

RESULTS

In the simulation study, the proposed method outperformed the other algorithms by at least 3.11% in the signal-to-noise ratio (SNR) and 0.68% in the contrast recovery coefficient (CRC), and it reduced the mean absolute error (MAE) by 8.07%. Regarding the tumor in the reconstructed image, the proposed method contained more pathological information. Furthermore, the proposed method was still superior to the MR-guided KEM method in the mismatch experiments.

CONCLUSIONS

The proposed FDG-guided KEM algorithm can effectively utilize and compensate for the tissue metabolism information obtained from dual-tracer PET to maximize the advantages of PET imaging.

(An overview of) Synergistic reconstruction for multimodality/multichannel imaging methods

Imaging is omnipresent in modern society with imaging devices based on a zoo of physical principles, probing a specimen across different wavelengths, energies and time. Recent years have seen a change in the imaging landscape with more and more imaging devices combining that which previously was used separately. Motivated by these hardware developments, an ever increasing set of mathematical ideas is appearing regarding how data from different imaging modalities or channels can be synergistically combined in the image reconstruction process, exploiting structural and/or functional correlations between the multiple images. Here we review these developments, give pointers to important challenges and provide an outlook as to how the field may develop in the forthcoming years. This article is part of the theme issue 'Synergistic tomographic image reconstruction: part 1'.

Applications of artificial intelligence in nuclear medicine image generation

Recently, the application of artificial intelligence (AI) in medical imaging (including nuclear medicine imaging) has rapidly developed. Most AI applications in nuclear medicine imaging have focused on the diagnosis, treatment monitoring, and correlation analyses with pathology or specific gene mutation. It can also be used for image generation to shorten the time of image acquisition, reduce the dose of injected tracer, and enhance image quality. This work provides an overview of the application of AI in image generation for single-photon emission computed tomography (SPECT) and positron emission tomography (PET) either without or with anatomical information [CT or magnetic resonance imaging (MRI)]. This review focused on four aspects, including imaging physics, image reconstruction, image postprocessing, and internal dosimetry. AI application in generating attenuation map, estimating scatter events, boosting image quality, and predicting internal dose map is summarized and discussed.

Low activity [11C]raclopride kinetic modeling in the mouse brain using the spatiotemporal kernel method

Depending on the molar activity of the tracer, the maximal allowable injected activity in mouse brain PET studies can be extremely low in order to avoid receptor saturation. Therefore, a high level of noise can be present in the image. We investigate several dynamic PET reconstruction methods in reduced counts, or equivalently in reduced injected activity, data exemplified in [¹¹C]racloprideBP_NDandR1quantification using the simplified reference tissue model (SRTM). We compared independent frame reconstruction (IFR), post-reconstruction HYPR denoising (IFR + HYPR), direct reconstruction using the SRTM model (DIR-SRTM), and the spatial (KERS) and spatiotemporal kernel reconstruction (KERST). Additionally, HYPR denoising of the frames used as features for the calculation of the spatial kernel matrix, was investigated (KERS-HYPR and KERST-HYPR).In vivodata of 11 mice, was used to generate list-mode data for five reduced count levels corresponding to reductions by a factor 4, 8, 12, 16 and 32 (equivalently 2.07, 1.04, 0.691, 0.518, and 0.260 MBq). Correlation of regionalBP_NDandR1values (reduced versus full counts reconstructions) was high (r > 0.94) for all methods, with KERS-HYPR and KERST-HYPR reaching the highest correlation (r > 0.96). Among methods with regularization, DIR-SRTM showed the largest variability inBP_ND(Bland-Altman SD from 3.0% to 12%), while IFR showed it forR1(5.1%-14.6%). KERST and KERST-HYPR were the only methods with Bland-Altman bias and SD below 5% for noise level up to a reduction factor of 16. At the voxel level,BP_NDandR1correlation was gradually decreased with increasing noise, with the largest correlation (BP_NDr > 0.88,R1r > 0.62) for KERS-HYPR and KERST-HYPR. The spatial and the spatiotemporal kernel methods performed similarly, while using only temporal regularization with direct reconstruction showed more variability. AlthoughR1 values present noise, using the spatiotemporal kernel reconstruction, accurate estimates of binding potential could be obtained with mouse injected activities as low as 0.26-0.518 MBq. This is desirable in order to maintain the tracer kinetics principle in mouse studies.

3D/4D Reconstruction and Quantitative Total Body Imaging

Total-body PET image reconstruction follows a similar procedure to the image reconstruction process for standard whole-body PET scanners. One unique aspect of total-body imaging is simultaneous coverage of the entire human body, which makes it convenient to perform total-body dynamic PET scans. Therefore, four-dimensional dynamic PET reconstruction and parametric imaging are of great interest in total-body imaging. This article covers some basics of PET image reconstruction and then focuses on three- and four-dimensional PET reconstruction for total-body imaging. Methods for image formation from raw measurements in total-body PET are described. Challenges and opportunities in total-body PET image reconstruction are discussed.

PET-enabled dual-energy CT: image reconstruction and a proof-of-concept computer simulation study

Standard dual-energy computed tomography (CT) uses two different x-ray energies to obtain energy-dependent tissue attenuation information to allow quantitative material decomposition. The combined use of dual-energy CT and positron emission tomography (PET) may provide a more comprehensive characterization of disease states in cancer and other diseases. However, the integration of dual-energy CT with PET is not trivial, either requiring costly hardware upgrades or increasing radiation exposure. This paper proposes a different dual-energy CT imaging method that is enabled by PET. Instead of using a second x-ray CT scan with a different energy, this method exploits time-of-flight PET image reconstruction via the maximum likelihood attenuation and activity (MLAA) algorithm to obtain a 511 keV gamma-ray attenuation image from PET emission data. The high-energy gamma-ray attenuation image is then combined with the low-energy x-ray CT of PET/CT to provide a pair of dual-energy CT images. A major challenge with the standard MLAA reconstruction is the high noise present in the reconstructed 511 keV attenuation map, which would not compromise the PET activity reconstruction too much but may significantly affect the performance of the gamma-ray attenuation image for material decomposition. To overcome the problem, we further propose a kernel MLAA algorithm to exploit the prior information from the available x-ray CT image. We conducted a computer simulation to test the concept and algorithm for the task of material decomposition. The simulation results demonstrate that this PET-enabled dual-energy CT method is promising for quantitative material decomposition. The proposed method can be readily implemented on time-of-flight PET/CT scanners to enable simultaneous PET and dual-energy CT imaging.

PET Parametric Imaging: Past, Present, and Future

Positron emission tomography (PET) is actively used in a diverse range of applications in oncology, cardiology, and neurology. The use of PET in the clinical setting focuses on static (single time frame) imaging at a specific time-point post radiotracer injection and is typically considered as semi-quantitative; e.g. standardized uptake value (SUV) measures. In contrast, dynamic PET imaging requires increased acquisition times but has the advantage that it measures the full spatiotemporal distribution of a radiotracer and, in combination with tracer kinetic modeling, enables the generation of multiparametric images that more directly quantify underlying biological parameters of interest, such as blood flow, glucose metabolism, and receptor binding. Parametric images have the potential for improved detection and for more accurate and earlier therapeutic response assessment. Parametric imaging with dynamic PET has witnessed extensive research in the past four decades. In this paper, we provide an overview of past and present activities and discuss emerging opportunities in the field of parametric imaging for the future.

Multiparametric Cardiac 18F-FDG PET in Humans: Kinetic Model Selection and Identifiability Analysis

Cardiac ¹⁸F-FDG PET has been used in clinics to assess myocardial glucose metabolism. Its ability for imaging myocardial glucose transport, however, has rarely been exploited in clinics. Using the dynamic FDG-PET scans of ten patients with coronary artery disease, we investigate in this paper appropriate dynamic scan and kinetic modeling protocols for efficient quantification of myocardial glucose transport. Three kinetic models and the effect of scan duration were evaluated by using statistical fit quality, assessing the impact on kinetic quantification, and analyzing the practical identifiability. The results show that the kinetic model selection depends on the scan duration. The reversible two-tissue model was needed for a one-hour dynamic scan. The irreversible two-tissue model was optimal for a scan duration of around 10-15 minutes. If the scan duration was shortened to 2-3 minutes, a one-tissue model was the most appropriate. For global quantification of myocardial glucose transport, we demonstrated that an early dynamic scan with a duration of 10-15 minutes and irreversible kinetic modeling was comparable to the full one-hour scan with reversible kinetic modeling. Myocardial glucose transport quantification provides an additional physiological parameter on top of the existing assessment of glucose metabolism and has the potential to enable single tracer multiparametric imaging in the myocardium.

Spatial and spatio-temporal statistical analyses of retinal images: a review of methods and applications

Background

Clinical research and management of retinal diseases greatly depend on the interpretation of retinal images and often longitudinally collected images. Retinal images provide context for spatial data, namely the location of specific pathologies within the retina. Longitudinally collected images can show how clinical events at one point can affect the retina over time. In this review, we aimed to assess statistical approaches to spatial and spatio-temporal data in retinal images. We also review the spatio-temporal modelling approaches used in other medical image types.

Methods

We conducted a comprehensive literature review of both spatial or spatio-temporal approaches and non-spatial approaches to the statistical analysis of retinal images. The key methodological and clinical characteristics of published papers were extracted. We also investigated whether clinical variables and spatial correlation were accounted for in the analysis.

Results

Thirty-four papers that included retinal imaging data were identified for full-text information extraction. Only 11 (32.4%) papers used spatial or spatio-temporal statistical methods to analyse images, others (23 papers, 67.6%) used non-spatial methods. Twenty-eight (82.4%) papers reported images collected cross-sectionally, while 6 (17.6%) papers reported analyses on images collected longitudinally. In imaging areas outside of ophthalmology, 19 papers were identified with spatio-temporal analysis, and multiple statistical methods were recorded.

Conclusions

In future statistical analyses of retinal images, it will be beneficial to clearly define and report the spatial distributions studied, report the spatial correlations, combine imaging data with clinical variables into analysis if available, and clearly state the software or packages used.

Dynamic Cell Imaging in PET With Optimal Transport Regularization

We propose a novel dynamic image reconstruction method from PET listmode data that could be particularly suited to tracking single or small numbers of cells. In contrast to conventional PET reconstruction our method combines the information from all detected events not only to reconstruct the dynamic evolution of the radionuclide distribution, but also to improve the reconstruction at each single time point by enforcing temporal consistency. This is achieved via optimal transport regularization where in principle, among all possible temporally evolving radionuclide distributions consistent with the PET measurement, the one is chosen with least kinetic motion energy. The reconstruction is found by convex optimization so that there is no dependence on the initialization of the method. We study its behaviour on simulated data of a human PET system and demonstrate its robustness even in settings with very low radioactivity. In contrast to previously reported cell tracking algorithms, our technique is oblivious to the number of tracked cells. Without any additional complexity one or multiple cells can be reconstructed, and the model automatically determines the number of particles. For instance, four radiolabelled cells moving at a velocity of 3.1 mm/s and a PET recorded count rate of 1.1 cps (for each cell) could be simultaneously tracked with a tracking accuracy of 5.3 mm inside a simulated human body.

VVBP-Tensor in the FBP Algorithm: Its Properties and Application in Low-Dose CT Reconstruction

For decades, commercial X-ray computed tomography (CT) scanners have been using the filtered backprojection (FBP) algorithm for image reconstruction. However, the desire for lower radiation doses has pushed the FBP algorithm to its limit. Previous studies have made significant efforts to improve the results of FBP through preprocessing the sinogram, modifying the ramp filter, or postprocessing the reconstructed images. In this paper, we focus on analyzing and processing the stacked view-by-view backprojections (named VVBP-Tensor) in the FBP algorithm. A key challenge for our analysis lies in the radial structures in each backprojection slice. To overcome this difficulty, a sorting operation was introduced to the VVBP-Tensor in its z direction (the direction of the projection views). The results show that, after sorting, the tensor contains structures that are similar to those of the object, and structures in different slices of the tensor are correlated. We then analyzed the properties of the VVBP-Tensor, including structural self-similarity, tensor sparsity, and noise statistics. Considering these properties, we have developed an algorithm using the tensor singular value decomposition (named VVBP-tSVD) to denoise the VVBP-Tensor for low-mAs CT imaging. Experiments were conducted using a physical phantom and clinical patient data with different mAs levels. The results demonstrate that the VVBP-tSVD is superior to all competing methods under different reconstruction schemes, including sinogram preprocessing, image postprocessing, and iterative reconstruction. We conclude that the VVBP-Tensor is a suitable processing target for improving the quality of FBP reconstruction, and the proposed VVBP-tSVD is an effective algorithm for noise reduction in low-mAs CT imaging. This preliminary work might provide a heuristic perspective for reviewing and rethinking the FBP algorithm.

Intercomparison of MR-informed PET image reconstruction methods

PURPOSE

Numerous image reconstruction methodologies for positron emission tomography (PET) have been developed that incorporate magnetic resonance (MR) imaging structural information, producing reconstructed images with improved suppression of noise and reduced partial volume effects. However, the influence of MR structural information also increases the possibility of suppression or bias of structures present only in the PET data (PET-unique regions). To address this, further developments for MR-informed methods have been proposed, for example, through inclusion of the current reconstructed PET image, alongside the MR image, in the iterative reconstruction process. In this present work, a number of kernel and maximum a posteriori (MAP) methodologies are compared, with the aim of identifying methods that enable a favorable trade-off between the suppression of noise and the retention of unique features present in the PET data.

METHODS

The reconstruction methods investigated were: the MR-informed conventional and spatially compact kernel methods, referred to as KEM and KEM largest value sparsification (LVS) respectively; the MR-informed Bowsher and Gaussian MR-guided MAP methods; and the PET-MR-informed hybrid kernel and anato-functional MAP methods. The trade-off between improving the reconstruction of the whole brain region and the PET-unique regions was investigated for all methods in comparison with postsmoothed maximum likelihood expectation maximization (MLEM), evaluated in terms of structural similarity index (SSIM), normalized root mean square error (NRMSE), bias, and standard deviation. Both simulated BrainWeb (10 noise realizations) and real [18 F] fluorodeoxyglucose (FDG) three-dimensional datasets were used. The real [18 F]FDG dataset was augmented with simulated tumors to allow comparison of the reconstruction methodologies for the case of known regions of PET-MR discrepancy and evaluated at full counts (100%) and at a reduced (10%) count level.

RESULTS

For the high-count simulated and real data studies, the anato-functional MAP method performed better than the other methods under investigation (MR-informed, PET-MR-informed and postsmoothed MLEM), in terms of achieving the best trade-off for the reconstruction of the whole brain and PET-unique regions, assessed in terms of the SSIM, NRMSE, and bias vs standard deviation. The inclusion of PET information in the anato-functional MAP method enables the reconstruction of PET-unique regions to attain similarly low levels of bias as unsmoothed MLEM, while moderately improving the whole brain image quality for low levels of regularization. However, for low count simulated datasets the anato-functional MAP method performs poorly, due to the inclusion of noisy PET information in the regularization term. For the low counts simulated dataset, KEM LVS and to a lesser extent, HKEM performed better than the other methods under investigation in terms of achieving the best trade-off for the reconstruction of the whole brain and PET-unique regions, assessed in terms of the SSIM, NRMSE, and bias vs standard deviation.

CONCLUSION

For the reconstruction of noisy data, multiple MR-informed methods produce favorable whole brain vs PET-unique region trade-off in terms of the image quality metrics of SSIM and NRMSE, comfortably outperforming the whole image denoising of postsmoothed MLEM.