Positron Emission Tomography Data Sets Research Articles

TAI-GAN: A Temporally and Anatomically Informed Generative Adversarial Network for early-to-late frame conversion in dynamic cardiac PET inter-frame motion correction

Inter-frame motion in dynamic cardiac positron emission tomography (PET) using rubidium-82 (82Rb) myocardial perfusion imaging impacts myocardial blood flow (MBF) quantification and the diagnosis accuracy of coronary artery diseases. However, the high cross-frame distribution variation due to rapid tracer kinetics poses a considerable challenge for inter-frame motion correction, especially for early frames where intensity-based image registration techniques often fail. To address this issue, we propose a novel method called Temporally and Anatomically Informed Generative Adversarial Network (TAI-GAN) that utilizes an all-to-one mapping to convert early frames into those with tracer distribution similar to the last reference frame. The TAI-GAN consists of a feature-wise linear modulation layer that encodes channel-wise parameters generated from temporal information and rough cardiac segmentation masks with local shifts that serve as anatomical information. Our proposed method was evaluated on a clinical 82Rb PET dataset, and the results show that our TAI-GAN can produce converted early frames with high image quality, comparable to the real reference frames. After TAI-GAN conversion, the motion estimation accuracy and subsequent myocardial blood flow (MBF) quantification with both conventional and deep learning-based motion correction methods were improved compared to using the original frames. The code is available at https://github.com/gxq1998/TAI-GAN.

Population-based deep image prior for dynamic PET denoising: A data-driven approach to improve parametric quantification

The high noise level of dynamic Positron Emission Tomography (PET) images degrades the quality of parametric images. In this study, we aim to improve the quality and quantitative accuracy of Ki images by utilizing deep learning techniques to reduce the noise in dynamic PET images. We propose a novel denoising technique, Population-based Deep Image Prior (PDIP), which integrates population-based prior information into the optimization process of Deep Image Prior (DIP). Specifically, the population-based prior image is generated from a supervised denoising model that is trained on a prompts-matched static PET dataset comprising 100 clinical studies. The 3D U-Net architecture is employed for both the supervised model and the following DIP optimization process. We evaluated the efficacy of PDIP for noise reduction in 25%-count and 100%-count dynamic PET images from 23 patients by comparing with two other baseline techniques: the Prompts-matched Supervised model (PS) and a conditional DIP (CDIP) model that employs the mean static PET image as the prior. Both the PS and CDIP models show effective noise reduction but result in smoothing and removal of small lesions. In addition, the utilization of a single static image as the prior in the CDIP model also introduces a similar tracer distribution to the denoised dynamic frames, leading to lower Ki in general as well as incorrect Ki in the descending aorta. By contrast, as the proposed PDIP model utilizes intrinsic image features from the dynamic dataset and a large clinical static dataset, it not only achieves comparable noise reduction as the supervised and CDIP models but also improves lesion Ki predictions.

A multi-instance tumor subtype classification method for small PET datasets using RA-DL attention module guided deep feature extraction with radiomics features

Background:Positron emission tomography (PET) is extensively employed for diagnosing and staging various tumors, including liver cancer, lung cancer, and lymphoma. Accurate subtype classification of tumors plays a crucial role in formulating effective treatment plans for patients. Notably, lymphoma comprises subtypes like diffuse large B-cell lymphoma and Hodgkin’s lymphoma, while lung cancer encompasses adenocarcinoma, small cell carcinoma, and squamous cell carcinoma. Similarly, liver cancer consists of subtypes such as cholangiocarcinoma and hepatocellular carcinoma. Consequently, the subtype classification of tumors based on PET images holds immense clinical significance. However, in clinical practice, the number of cases available for each subtype is often limited and imbalanced. Therefore, the primary challenge lies in achieving precise subtype classification using a small dataset. Method:This paper presents a novel approach for tumor subtype classification in small datasets using RA-DL (Radiomics-DeepLearning) attention. To address the limited sample size, Support Vector Machines (SVM) is employed as the classifier for tumor subtypes instead of deep learning methods. Emphasizing the importance of texture information in tumor subtype recognition, radiomics features are extracted from the tumor regions during the feature extraction stage. These features are compressed using an autoencoder to reduce redundancy. In addition to radiomics features, deep features are also extracted from the tumors to leverage the feature extraction capabilities of deep learning. In contrast to existing methods, our proposed approach utilizes the RA-DL-Attention mechanism to guide the deep network in extracting complementary deep features that enhance the expressive capacity of the final features while minimizing redundancy. To address the challenges of limited and imbalanced data, our method avoids using classification labels during deep feature extraction and instead incorporates 2D Region of Interest (ROI) segmentation and image reconstruction as auxiliary tasks. Subsequently, all lesion features of a single patient are aggregated into a feature vector using a multi-instance aggregation layer. Result:Validation experiments were conducted on three PET datasets, specifically the liver cancer dataset, lung cancer dataset, and lymphoma dataset. In the context of lung cancer, our proposed method achieved impressive performance with Area Under Curve (AUC) values of 0.82, 0.84, and 0.83 for the three-classification task. For the binary classification task of lymphoma, our method demonstrated notable results with AUC values of 0.95 and 0.75. Moreover, in the binary classification task of liver tumor, our method exhibited promising performance with AUC values of 0.84 and 0.86. Conclusion:The experimental results clearly indicate that our proposed method outperforms alternative approaches significantly. Through the extraction of complementary radiomics features and deep features, our method achieves a substantial improvement in tumor subtype classification performance using small PET datasets.

An iterative reconstruction algorithm for unsupervised PET image

Objective. In recent years, convolutional neural networks (CNNs) have shown great potential in positron emission tomography (PET) image reconstruction. However, most of them rely on many low-quality and high-quality reference PET image pairs for training, which are not always feasible in clinical practice. On the other hand, many works improve the quality of PET image reconstruction by adding explicit regularization or optimizing the network structure, which may lead to complex optimization problems. Approach. In this paper, we develop a novel iterative reconstruction algorithm by integrating the deep image prior (DIP) framework, which only needs the prior information (e.g. MRI) and sinogram data of patients. To be specific, we construct the objective function as a constrained optimization problem and utilize the existing PET image reconstruction packages to streamline calculations. Moreover, to further improve both the reconstruction quality and speed, we introduce the Nesterov’s acceleration part and the restart mechanism in each iteration. Main results. 2D experiments on PET data sets based on computer simulations and real patients demonstrate that our proposed algorithm can outperform existing MLEM-GF, KEM and DIPRecon methods. Significance. Unlike traditional CNN methods, the proposed algorithm does not rely on large data sets, but only leverages inter-patient information. Furthermore, we enhance reconstruction performance by optimizing the iterative algorithm. Notably, the proposed method does not require much modification of the basic algorithm, allowing for easy integration into standard implementations.

Partial volume correction of PET image data using geometric transfer matrices based on uniform B-splines

Objective. Most methods for partial volume correction (PVC) of positron emission tomography (PET) data employ anatomical segmentation of images into regions of interest. This approach is not optimal for exploratory functional imaging beyond regional hypotheses. Here, we describe a novel method for unbiased voxel-wise PVC. Approach. B-spline basis functions were combined with geometric transfer matrices to enable a method (bsGTM) that provides PVC or alternatively provides smoothing with minimal regional crosstalk. The efficacy of the proposed method was evaluated using Monte Carlo simulations, human PET data, and murine functional PET data. Main results. In simulations, bsGTM provided recovery of partial volume signal loss comparable to iterative deconvolution, while demonstrating superior resilience to noise. In a real murine PET dataset, bsGTM yielded much higher sensitivity for detecting amphetamine-induced reduction of [11C]raclopride binding potential. In human PET data, bsGTM smoothing enabled increased signal-to-noise ratios with less degradation of binding potentials relative to Gaussian convolution or non-local means. Significance. bsGTM offers improved performance for PVC relative to iterative deconvolution, the current method of choice for voxel-wise PVC, especially in the common PET regime of low signal-to-noise ratio. The new method provides an anatomically unbiased way to compensate partial volume errors in cases where anatomical segmentation is unavailable or of questionable relevance or accuracy.

Comparison of simple augmentation transformations for a convolutional neural network classifying medical images

Simple image augmentation techniques, such as reflection, rotation, or translation, might work differently for medical images than they do for regular photographs due to the fundamental properties of medical imaging techniques and the bilateral symmetry of the human body. Here, we compare the predictions of a convolutional neural network (CNN) trained for binary classification by using either no augmentation or one of seven usual types augmentation. We have 11 different medical data sets, mostly related to lung infections or cancer, with X-rays, ultrasound (US) images, and images from positron emission tomography (PET) and magnetic resonance imaging (MRI). According to our results, the augmentation types do not produce statistically significant differences for US and PET data sets, but, for X-rays and MRI images, the best augmentation technique is adding Gaussian blur to images.

Anatomically Guided PET Image Reconstruction Using Conditional Weakly-Supervised Multi-Task Learning Integrating Self-Attention.

To address the lack of high-quality training labels in positron emission tomography (PET) imaging, weakly-supervised reconstruction methods that generate network-based mappings between prior images and noisy targets have been developed. However, the learned model has an intrinsic variance proportional to the average variance of the target image. To suppress noise and improve the accuracy and generalizability of the learned model, we propose a conditional weakly-supervised multi-task learning (MTL) strategy, in which an auxiliary task is introduced serving as an anatomical regularizer for the PET reconstruction main task. In the proposed MTL approach, we devise a novel multi-channel self-attention (MCSA) module that helps learn an optimal combination of shared and task-specific features by capturing both local and global channel-spatial dependencies. The proposed reconstruction method was evaluated on NEMA phantom PET datasets acquired at different positions in a PET/CT scanner and 26 clinical whole-body PET datasets. The phantom results demonstrate that our method outperforms state-of-the-art learning-free and weakly-supervised approaches obtaining the best noise/contrast tradeoff with a significant noise reduction of approximately 50.0% relative to the maximum likelihood (ML) reconstruction. The patient study results demonstrate that our method achieves the largest noise reductions of 67.3% and 35.5% in the liver and lung, respectively, as well as consistently small biases in 8 tumors with various volumes and intensities. In addition, network visualization reveals that adding the auxiliary task introduces more anatomical information into PET reconstruction than adding only the anatomical loss, and the developed MCSA can abstract features and retain PET image details.

Deep learning approach using SPECT-to-PET translation for attenuation correction in CT-less myocardial perfusion SPECT imaging

ObjectiveDeep learning approaches have attracted attention for improving the scoring accuracy in computed tomography-less single photon emission computed tomography (SPECT). In this study, we proposed a novel deep learning approach referring to positron emission tomography (PET). The aims of this study were to analyze the agreement of representative voxel values and perfusion scores of SPECT-to-PET translation model-generated SPECT (SPECTSPT) against PET in 17 segments according to the American Heart Association (AHA).MethodsThis retrospective study evaluated the patient-to-patient stress, resting SPECT, and PET datasets of 71 patients. The SPECTSPT generation model was trained (stress: 979 image pairs, rest: 987 image pairs) and validated (stress: 421 image pairs, rest: 425 image pairs) using 31 cases of SPECT and PET image pairs using an image-to-image translation network. Forty of 71 cases of left ventricular base-to-apex short-axis images were translated to SPECTSPT in the stress and resting state (stress: 1830 images, rest: 1856 images). Representative voxel values of SPECT and SPECTSPT in the 17 AHA segments against PET were compared. The stress, resting, and difference scores of 40 cases of SPECT and SPECTSPT were also compared in each of the 17 segments.ResultsFor AHA 17-segment-wise analysis, stressed SPECT but not SPECTSPT voxel values showed significant error from PET at basal anterior regions (segments #1, #6), and at mid inferoseptal regions (segments #8, #9, and #10). SPECT, but not SPECTSPT, voxel values at resting state showed significant error at basal anterior regions (segments #1, #2, and #6), and at mid inferior regions (segments #8, #9, and #11). Significant SPECT overscoring was observed against PET in basal-to-apical inferior regions (segments #4, #10, and #15) during stress. No significant overscoring was observed in SPECTSPT at stress, and only moderate over and underscoring in the basal inferior region (segment #4) was found in the resting and difference states.ConclusionsOur PET-supervised deep learning model is a new approach to correct well-known inferior wall attenuation in SPECT myocardial perfusion imaging. As standalone SPECT systems are used worldwide, the SPECTSPT generation model may be applied as a low-cost and practical clinical tool that provides powerful auxiliary information for the diagnosis of myocardial blood flow.

Positron emission tomography dataset of [11C]carbon dioxide storage in coal for geo-sequestration application

Positron Emission Tomography (PET) imaging has demonstrated its capability in providing time-lapse fluid flow visualisation for improving the understanding of flow properties of geologic media. To investigate the process of CO2 geo-sequestration using PET imaging technology, [11C]CO2 is the most optimal and direct radiotracer. However, it has not been extensively used due to the short half-life of Carbon-11 (20.4 minutes). In this work, a novel laboratory protocol is developed to use [11C]CO2 as radiolabelled tracer to visualise and quantify in-situ CO2 adsorption, spreading, diffusion, and advection flow in coal. This protocol consists of generation and delivering of [11C]CO2, lab-based PET scanning, subsequent micro-CT scanning, and data processing. The lab-based PET scanning setup integrates in-situ core flooding tests with PET scanning. The real-time PET images are acquired under different storage conditions, including early gas production stage, depleted stage, and late storage stage. These datasets can be used to study across-scale theoretical and experimental study of CO2 flow behaviour in coal with the application to CO2 geo-sequestration.

PET/CT Fluoroscopy during PET/CT-Guided Interventions: Initial Experience

This study assessed the feasibility and functionality of the use of a high-speed image fusion technology to generate and display positron emission tomography (PET)/computed tomography (CT) fluoroscopic images during PET/CT-guided tumor ablation procedures. Thirteen patients underwent 14 PET/CT-guided ablations for the treatment of 20 tumors. A Food and Drug Administration–cleared multimodal image fusion platform received images pushed from a scanner, followed by near–real-time, nonrigid image registration. The most recent intraprocedural PET dataset was fused to each single-rotation CT fluoroscopy dataset as it arrived, and the fused images were displayed on an in-room monitor. PET/CT fluoroscopic images were generated and displayed in all procedures and enabled more confident targeting in 3 procedures. The mean lag time from CT fluoroscopic image acquisition to the in-room display of the fused PET/CT fluoroscopic image was 21 seconds ± 8. The registration accuracy was visually satisfactory in 13 of 14 procedures. In conclusion, PET/CT fluoroscopy was feasible and may have the potential to facilitate PET/CT-guided procedures.

PET scatter estimation using deep learning U-Net architecture

Objective. Positron emission tomography (PET) image reconstruction needs to be corrected for scatter in order to produce quantitatively accurate images. Scatter correction is traditionally achieved by incorporating an estimated scatter sinogram into the forward model during image reconstruction. Existing scatter estimated methods compromise between accuracy and computing time. Nowadays scatter estimation is routinely performed using single scatter simulation (SSS), which does not accurately model multiple scatter and scatter from outside the field-of-view, leading to reduced qualitative and quantitative PET reconstructed image accuracy. On the other side, Monte-Carlo (MC) methods provide a high precision, but are computationally expensive and time-consuming, even with recent progress in MC acceleration. Approach. In this work we explore the potential of deep learning (DL) for accurate scatter correction in PET imaging, accounting for all scatter coincidences. We propose a network based on a U-Net convolutional neural network architecture with 5 convolutional layers. The network takes as input the emission and computed tomography (CT)-derived attenuation factor (AF) sinograms and returns the estimated scatter sinogram. The network training was performed using MC simulated PET datasets. Multiple anthropomorphic extended cardiac-torso phantoms of two different regions (lung and pelvis) were created, considering three different body sizes and different levels of statistics. In addition, two patient datasets were used to assess the performance of the method in clinical practice. Main results. Our experiments showed that the accuracy of our method, namely DL-based scatter estimation (DLSE), was independent of the anatomical region (lungs or pelvis). They also showed that the DLSE-corrected images were similar to that reconstructed from scatter-free data and more accurate than SSS-corrected images. Significance. The proposed method is able to estimate scatter sinograms from emission and attenuation data. It has shown a better accuracy than the SSS, while being faster than MC scatter estimation methods.

TAI-GAN: Temporally and Anatomically Informed GAN for Early-to-Late Frame Conversion in Dynamic Cardiac PET Motion Correction.

The rapid tracer kinetics of rubidium-82 (82Rb) and high variation of cross-frame distribution in dynamic cardiac positron emission tomography (PET) raise significant challenges for inter-frame motion correction, particularly for the early frames where conventional intensity-based image registration techniques are not applicable. Alternatively, a promising approach utilizes generative methods to handle the tracer distribution changes to assist existing registration methods. To improve frame-wise registration and parametric quantification, we propose a Temporally and Anatomically Informed Generative Adversarial Network (TAI-GAN) to transform the early frames into the late reference frame using an all-to-one mapping. Specifically, a feature-wise linear modulation layer encodes channel-wise parameters generated from temporal tracer kinetics information, and rough cardiac segmentations with local shifts serve as the anatomical information. We validated our proposed method on a clinical 82Rb PET dataset and found that our TAI-GAN can produce converted early frames with high image quality, comparable to the real reference frames. After TAI-GAN conversion, motion estimation accuracy and clinical myocardial blood flow (MBF) quantification were improved compared to using the original frames. Our code is published at https://github.com/gxq1998/TAI-GAN.

Virtual high-count PET image generation using a deep learning method.

Recently, deep learning-based methods have been established to denoise the low-count positron emission tomography (PET) images and predict their standard-count image counterparts, which could achieve reduction of injected dosage and scan time, and improve image quality for equivalent lesion detectability and clinical diagnosis. In clinical settings, the majority scans are still acquired using standard injection dose with standard scan time. In this work, we applied a 3D U-Net network to reduce the noise of standard-count PET images to obtain the virtual-high-count (VHC) PET images for identifying the potential benefits of the obtained VHC PETimages. The training datasets, including down-sampled standard-count PET images as the network input and high-count images as the desired network output, were derived from 27 whole-body PET datasets, which were acquired using 90-min dynamic scan. The down-sampled standard-count PET images were rebinned with matched noise level of 195 clinical static PET datasets, by matching the normalized standard derivation (NSTD) inside 3D liver region of interests (ROIs). Cross-validation was performed on 27 PET datasets. Normalized mean square error (NMSE), peak signal to noise ratio (PSNR), structural similarity index (SSIM), and standard uptake value (SUV) bias of lesions were used for evaluation on standard-count and VHC PET images, with real-high-count PET image of 90 min as the gold standard. In addition, the network trained with 27 dynamic PET datasets was applied to 195 clinical static datasets to obtain VHC PET images. The NSTD and mean/max SUV of hypermetabolic lesions in standard-count and VHC PET images were evaluated. Three experienced nuclear medicine physicians evaluated the overall image quality of randomly selected 50 out of 195 patients' standard-count and VHC images and conducted 5-score ranking. A Wilcoxon signed-rank test was used to compare differences in the grading of standard-count and VHCimages. The cross-validation results showed that VHC PET images had improved quantitative metrics scores than the standard-count PET images. The mean/max SUVs of 35 lesions in the standard-count and true-high-count PET images did not show significantly statistical difference. Similarly, the mean/max SUVs of VHC and true-high-count PET images did not show significantly statistical difference. For the 195 clinical data, the VHC PET images had a significantly lower NSTD than the standard-count images. The mean/max SUVs of 215 hypermetabolic lesions in the VHC and standard-count images showed no statistically significant difference. In the image quality evaluation by three experienced nuclear medicine physicians, standard-count images and VHC images received scores with mean and standard deviation of 3.34±0.80 and 4.26 ± 0.72 from Physician 1, 3.02 ± 0.87 and 3.96 ± 0.73 from Physician 2, and 3.74 ± 1.10 and 4.58 ± 0.57 from Physician 3, respectively. The VHC images were consistently ranked higher than the standard-count images. The Wilcoxon signed-rank test also indicated that the image quality evaluation between standard-count and VHC images had significantdifference. A DL method was proposed to convert the standard-count images to the VHC images. The VHC images had reduced noise level. No significant difference in mean/max SUV to the standard-count images was observed. VHC images improved image quality for better lesion detectability and clinicaldiagnosis.

Kernel graph filtering-A new method for dynamic sinogram denoising.

Low count PET (positron emission tomography) imaging is often desirable in clinical diagnosis and biomedical research, but its images are generally very noisy, due to the very weak signals in the sinograms used in image reconstruction. To address this issue, this paper presents a novel kernel graph filtering method for dynamic PET sinogram denoising. This method is derived from treating the dynamic sinograms as the signals on a graph, and learning the graph adaptively from the kernel principal components of the sinograms to construct a lowpass kernel graph spectrum filter. The kernel graph filter thus obtained is then used to filter the original sinogram time frames to obtain the denoised sinograms for PET image reconstruction. Extensive tests and comparisons on the simulated and real life in-vivo dynamic PET datasets show that the proposed method outperforms the existing methods in sinogram denoising and image enhancement of dynamic PET at all count levels, especially at low count, with a great potential in real life applications of dynamic PET imaging.

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

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. This paper proposes incorporating 18 F-fluorodeoxyglucose (FDG) as a higher-quality PET image to guide the reconstruction of other lower-count 11 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. 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. 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.

Convolutional neural networks for PET functional volume fully automatic segmentation: development and validation in a multi-center setting

PurposeIn this work, we addressed fully automatic determination of tumor functional uptake from positron emission tomography (PET) images without relying on other image modalities or additional prior constraints, in the context of multicenter images with heterogeneous characteristics.MethodsIn cervical cancer, an additional challenge is the location of the tumor uptake near or even stuck to the bladder. PET datasets of 232 patients from five institutions were exploited. To avoid unreliable manual delineations, the ground truth was generated with a semi-automated approach: a volume containing the tumor and excluding the bladder was first manually determined, then a well-validated, semi-automated approach relying on the Fuzzy locally Adaptive Bayesian (FLAB) algorithm was applied to generate the ground truth. Our model built on the U-Net architecture incorporates residual blocks with concurrent spatial squeeze and excitation modules, as well as learnable non-linear downsampling and upsampling blocks. Experiments relied on cross-validation (four institutions for training and validation, and the fifth for testing).ResultsThe model achieved good Dice similarity coefficient (DSC) with little variability across institutions (0.80 ± 0.03), with higher recall (0.90 ± 0.05) than precision (0.75 ± 0.05) and improved results over the standard U-Net (DSC 0.77 ± 0.05, recall 0.87 ± 0.02, precision 0.74 ± 0.08). Both vastly outperformed a fixed threshold at 40% of SUVmax (DSC 0.33 ± 0.15, recall 0.52 ± 0.17, precision 0.30 ± 0.16). In all cases, the model could determine the tumor uptake without including the bladder. Neither shape priors nor anatomical information was required to achieve efficient training.ConclusionThe proposed method could facilitate the deployment of a fully automated radiomics pipeline in such a challenging multicenter context.

Use of radiomics based on 18F-FDG PET/CT and machine learning methods to aid clinical decision-making in the classification of solitary pulmonary lesions: an innovative approach.

This study was designed and performed to assess the ability of 18F-fluorodeoxyglucose (FDG) positron emission tomography (PET) and computed tomography (CT) radiomics features combined with machine learning methods to differentiate between primary and metastatic lung lesions and to classify histological subtypes. Moreover, we identified the optimal machine learning method. A total of 769 patients pathologically diagnosed with primary or metastatic lung cancers were enrolled. We used the LIFEx package to extract radiological features from semiautomatically segmented PET and CT images within the same volume of interest. Patients were randomly distributed in training and validation sets. Through the evaluation of five feature selection methods and nine classification methods, discriminant models were established. The robustness of the procedure was controlled by tenfold cross-validation. The model's performance was evaluated using the area under the receiver operating characteristic curve (AUC). Based on the radiomics features extracted from PET and CT images, forty-five discriminative models were established. Combined with appropriate feature selection methods, most classifiers showed excellent discriminative ability with AUCs greater than 0.75. In the differentiation between primary and metastatic lung lesions, the feature selection method gradient boosting decision tree (GBDT) combined with the classifier GBDT achieved the highest classification AUC of 0.983 in the PET dataset. In contrast, the feature selection method eXtreme gradient boosting combined with the classifier random forest (RF) achieved the highest AUC of 0.828 in the CT dataset. In the discrimination between squamous cell carcinoma and adenocarcinoma, the combination of GBDT feature selection method with GBDT classification had the highest AUC of 0.897 in the PET dataset. In contrast, the combination of the GBDT feature selection method with the RF classification had the highest AUC of 0.839 in the CT dataset. Most of the decision tree (DT)-based models were overfitted, suggesting that the classification method was not appropriate for practical application. 18F-FDG PET/CT radiomics features combined with machine learning methods can distinguish between primary and metastatic lung lesions and identify histological subtypes in lung cancer. GBDT and RF were considered optimal classification methods for the PET and CT datasets, respectively, and GBDT was considered the optimal feature selection method in our analysis.

An efficient method for PET image denoising by combining multi-scale transform and non-local means

The diagnosis of dementia, particularly in the early stages is very much helpful with Positron emission tomography (PET) image processing. The most important challenges in PET image processing are noise removal and region of interests (ROIs) segmentation. Although denoising and segmentation are performed independently, but the performance of the denoising process significantly affects the performance of the segmentation process. Due to the low signals to noise ratio and low contrast, PET image denoising is a challenging task. Individual wavelet, curvelet and non-local means (NLM) based methods are not well suited to handle both isotropic (smooth details) and anisotropic (edges and curves) features due to its restricted abilities. To address these issues, the present work proposes an efficient denoising framework for reducing the noise level of brain PET images based on the combination of multi-scale transform (wavelet and curvelet) and tree clustering non-local means (TNLM). The main objective of the proposed method is to extract the isotropic features from a noisy smooth PET image using tree clustering based non-local means (TNLM). Then curvelet-based denoising is applied to the residual image to extract the anisotropic features such as edges and curves. Finally, the extracted anisotropic features are inserted back into the isotropic features to obtain an estimated denoised image. Simulated phantom and clinical PET datasets have been used in this proposed work for testing and measuring the performance in the medical applications, such as gray matter segmentation and precise tumor region identification without any interaction with other structural images like MRI or CT. The results in the experimental section show that the proposed denoising method has obtained better performance than existing wavelet, curvelet, wavelet-curvelet, non-local means (NLM) and deep learning methods based on the preservation of the edges. Qualitatively, a notable gain is achieved in the proposed denoised PET images in terms of contrast enhancement than other existing denoising methods.

Myocardial creep-induced misalignment artifacts in PET/MR myocardial perfusion imaging

PurposeMisalignment between positron emission tomography (PET) datasets and attenuation correction (AC) maps is a potential source of artifacts in myocardial perfusion imaging (MPI). We assessed the impact of adenosine on the alignment of AC maps derived from magnetic resonance (MR) and PET datasets during MPI on a hybrid PET/MR scanner.MethodsTwenty-eight volunteers underwent adenosine stress and rest 13N-ammonia MPI on a PET/MR. We acquired Dixon sequences for the creation of MRAC maps. After reconstruction of the original non-shifted PET images, we examined MRAC and PET datasets for cardiac spatial misalignment and, if necessary, reconstructed a second set of shifted PET images after manually adjusting co-registration. Summed rest, stress, and difference scores (SRS, SSS, and SDS) were compared between shifted and non-shifted PET images. Additionally, we measured the amount of cranial movement of the heart (i.e., myocardial creep) after termination of adenosine infusion.ResultsRealignment was necessary for 25 (89.3%) stress and 12 (42.9%) rest PET datasets. Median SRS, SSS, and SDS of the non-shifted images were 6 (IQR = 4–7), 12 (IQR = 7–18), and 8 (IQR = 2–11), respectively, and of the shifted images 2 (IQR = 1–6), 4 (IQR = 7–18), and 1 (IQR = 0–2), respectively. All three scores were significantly higher in non-shifted versus shifted images (all p < 0.05). The difference in SDS correlated moderately but significantly with the amount of myocardial creep (r = 0.541, p = 0.005).ConclusionMisalignment of MRAC and PET datasets commonly occurs during adenosine stress MPI on a hybrid PET/MR device, potentially leading to an increase in false-positive findings. Our results suggest that myocardial creep may substantially account for this and prompt for a careful review and correction of PET/MRAC data.

Brain activity detection in single- and multi-subject PET data by Bayesian analysis

Positron emission tomography (PET) is a functional neuroimaging method that maps brain activity non-invasively. Statistical methods play an essential part in understanding and analysing functional PET data. Several Bayesian approaches have been proposed for neuroimaging techniques that arrange information on the brain structure or activity function. In this paper, Bayesian analysis was used to detect functional brain activity in single- and multiple-subject PET data. Free PET dataset analyses for a single subject and five multiple subjects were conducted. A total of 72 functional PET images were processed, 12 for each one of the five multiple subjects and for the single subject. Several ways to design multiple-subject PET data were introduced in this work. Bayesian analysis was performed on the five multiple-subject and the single-subject PET data. A comparison was presented to determine which statistical matrix design is applicable for brain detection activity in PET data. The results of the design matrix and brain activity detection were presented for each selected design matrix. The Bayesian estimation of each case of the PET dataset for all the subjects was plotted. The brain activity was plotted as voxels on a transparent brain image in three orthogonal planes. The voxels were visualised using the maximum intensity projection method. Results showed that brain activity could not be detected easily in single-subject PET data. Finding the activity in multiple subjects depended on the design matrix used for PET data analysis.