PET NEMA IQ Phantom dataset: image reconstruction settings for quantitative PET imaging.
Ontology highlight
ABSTRACT:
SUBMITTER: Vosoughi H
PROVIDER: S-EPMC8237526 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:PurposeMany neurological diseases affect small structures in the brain and, as such, reliable visual evaluation and accurate quantification are required. Recent technological developments made the clinical use of hybrid positron emission tomography/magnetic resonance (PET/MR) systems possible, providing both functional and anatomical information in a single imaging session. Nevertheless, there is a trade-off between spatial resolution and image quality (contrast and noise), which is dictated mainly by the chosen acquisition and reconstruction protocols. Image reconstruction algorithms using spherical symmetric basis functions (blobs) for image representation have a number of additional parameters that impact both the qualitative and quantitative image characteristics. Hence, a detailed investigation of the blob-based reconstruction characteristics using different parameters is needed to achieve optimal reconstruction results.ProceduresThis work evaluated the impact of a range of blob parameters on image quality and quantitative accuracy of brain PET images acquired on the Ingenuity Time-of-Flight (TOF) PET/MR system. Two different phantoms were used to simulate brain imaging applications. Image contrast and noise characteristics were assessed using an image quality phantom. Quantitative performance in a clinical setting was investigated using the Hoffman 3D brain phantom at various count levels. Furthermore, the visual quality of four clinical studies was scored blindly by two experienced physicians to qualitatively evaluate the influence of different reconstruction protocols, hereby providing indications on parameters producing the best image quality.ResultsQuantitative evaluation using the image quality phantom showed that larger basis function radii result in lower contrast recovery (∼2%) and lower variance levels (∼15%). The brain phantom and clinical studies confirmed these observations since lower contrast was seen between anatomical structures. High and low count statistics gave comparable values. The qualitative evaluation of the clinical studies, based on the assessment performed by the physicians, showed a preference towards lower image variance levels with a slightly lower contrast, favoring higher radii and four iterations.ConclusionThis study systematically evaluated a number of basis function parameters and their quantitative and qualitative effect within PET image reconstruction in the context of brain imaging. A range of blob parameters can minimize error and provided optimal image quality, where the anatomical structures could be recognized but the exact delineation of these structures is simplified in scans with lower image variance levels and thus, higher radii should be preferred. With the optimization of blob parameters, the reconstructed images were found to be qualitatively improved (optimum parameters {d = 2.0375, alpha = 10.4101, radius = 3.9451}) as assessed by the physicians compared to the current clinical protocol. However, this qualitative improvement is task specific, depending on the desired image characteristics to be extracted. Finally, this work could be used as a guide for application-specific optimal parameter selection.
| S-EPMC4768229 | biostudies-literature
Project description:Recently, deep neural networks have been widely and successfully applied in computer vision tasks and have attracted growing interest in medical imaging. One barrier for the application of deep neural networks to medical imaging is the need for large amounts of prior training pairs, which is not always feasible in clinical practice. This is especially true for medical image reconstruction problems, where raw data are needed. Inspired by the deep image prior framework, in this paper, we proposed a personalized network training method where no prior training pairs are needed, but only the patient's own prior information. The network is updated during the iterative reconstruction process using the patient-specific prior information and measured data. We formulated the maximum-likelihood estimation as a constrained optimization problem and solved it using the alternating direction method of multipliers algorithm. Magnetic resonance imaging guided positron emission tomography reconstruction was employed as an example to demonstrate the effectiveness of the proposed framework. Quantification results based on simulation and real data show that the proposed reconstruction framework can outperform Gaussian post-smoothing and anatomically guided reconstructions using the kernel method or the neural-network penalty.
| S-EPMC6584077 | biostudies-literature
Project description:Positron emission tomography (PET) with 15O-tracers is commonly used to measure brain hemodynamic parameters such as cerebral blood flow, cerebral blood volume, and cerebral metabolic rate of oxygen. Conventionally, the absolute quantification of these parameters requires an arterial input function that is obtained invasively by sampling blood from an artery. In this work, we developed and validated an image-derived arterial input function technique that avoids the unreliable and burdensome arterial sampling procedure for full quantitative 15O-PET imaging. We then compared hemodynamic PET imaging performed on a PET/MR hybrid scanner against a conventional PET only scanner. We demonstrated the proposed imaging-based technique was able to generate brain hemodynamic parameter measurements in strong agreement with the traditional arterial sampling based approach. We also demonstrated that quantitative 15O-PET imaging can be successfully implemented on a PET/MR hybrid scanner.
| S-EPMC5453463 | biostudies-literature
Project description:ObjectiveSimultaneous PET/MRIs vary in their quantitative PET performance due to inherent differences in the physical systems and differences in the image reconstruction implementation. This variability in quantitative accuracy confounds the ability to meaningfully combine and compare data across scanners. In this work, we define image reconstruction parameters that lead to comparable contrast recovery curves across simultaneous PET/MRI systems.MethodThe NEMA NU-2 image quality phantom was imaged on one GE Signa and on one Siemens mMR PET/MRI scanner. The phantom was imaged at 9.7:1 contrast with standard spheres (diameter 10, 13, 17, 22, 28, 37 mm) and with custom spheres (diameter: 8.5, 11.5, 15, 25, 32.5, 44 mm) using a standardized methodology. Analysis was performed on a 30 min listmode data acquisition and on 6 realizations of 5 min from the listmode data. Images were reconstructed with the manufacturer provided iterative image reconstruction algorithms with and without point spread function (PSF) modeling. For both scanners, a post-reconstruction Gaussian filter of 3-7 mm in steps of 1 mm was applied. Attenuation correction was provided from a scaled computed tomography (CT) image of the phantom registered to the MR-based attenuation images and verified to align on the non-attenuation corrected PET images. For each of these image reconstruction parameter sets, contrast recovery coefficients (CRCs) were determined for the SUVmean, SUVmax and SUVpeak for each sphere. A hybrid metric combining the root-mean-squared discrepancy (RMSD) and the absolute CRC values was used to simultaneously optimize for best match in CRC between the two scanners while simultaneously weighting toward higher resolution reconstructions. The image reconstruction parameter set was identified as the best candidate reconstruction for each vendor for harmonized PET image reconstruction.ResultsThe range of clinically relevant image reconstruction parameters demonstrated widely different quantitative performance across cameras. The best match of CRC curves was obtained at the lowest RMSD values with: for CRCmean, 2 iterations-7 mm filter on the GE Signa and 4 iterations-6 mm filter on the Siemens mMR, for CRCmax, 4 iterations-6 mm filter on the GE Signa, 4 iterations-5 mm filter on the Siemens mMR and for CRCpeak, 4 iterations-7 mm filter with PSF on the GE Signa and 4 iterations-7 mm filter on the Siemens mMR. Over all reconstructions, the RMSD between CRCs was 1.8%, 3.6% and 2.9% for CRC mean, max and peak, respectively. The solution of 2 iterations-3 mm on the GE Signa and 4 iterations-3 mm on Siemens mMR, both with PSF, led to simultaneous harmonization and with high CRC and low RMSD for CRC mean, max and peak with RMSD values of 2.8%, 5.8% and 3.2%, respectively.ConclusionsFor two commercially available PET/MRI scanners, user-selectable parameters that control iterative updates, image smoothing and PSF modeling provide a range of contrast recovery curves that allow harmonization in harmonization strategies of optimal match in CRC or high CRC values. This work demonstrates that nearly identical CRC curves can be obtained on different commercially available scanners by selecting appropriate image reconstruction parameters.
| S-EPMC8571452 | biostudies-literature
Project description:Multi-tracer positron emission tomography (PET) has the potential to enhance PET imaging by providing complementary information from different physiological processes. However, one or more of the images may present high levels of noise. Guided image reconstruction methods transfer information from a guide image into the PET image reconstruction to encourage edge-preserving noise reduction. In this work we aim to reduce noise in poorer quality PET datasets via guidance from higher quality ones by using a weighted quadratic penalty approach. In particular, we applied this methodology to [18F]fluorodeoxyglucose (FDG) and [11C]methionine imaging of gliomas. 3D simulation studies showed that guiding the reconstruction of methionine datasets using pre-existing FDG images reduced reconstruction errors across the whole-brain (-8%) and within a tumour (-36%) compared to maximum likelihood expectation-maximisation (MLEM). Furthermore, guided reconstruction outperformed a comparable non-local means filter, indicating that regularising during reconstruction is preferable to post-reconstruction approaches. Hyperparameters selected from the 3D simulation study were applied to real data, where it was observed that the proposed FDG-guided methionine reconstruction allows for better edge preservation and noise reduction than standard MLEM. Overall, the results in this work demonstrate that transferring information between datasets in multi-tracer PET studies improves image quality and quantification performance.
| S-EPMC6130802 | biostudies-other
Project description:Unmanned aerial vehicles (UAVs) have become popular platforms for collecting various types of geospatial data for various mapping, monitoring and modelling applications. With the advancement of imaging and computing technologies, a vast variety of photogrammetric, computer-vision and, nowadays, end-to-end learning workflows are introduced to produce three-dimensional (3D) information in form of digital surface and terrain models, textured meshes, rectified mosaics, CAD models, etc. These 3D products might be used in applications where accuracy and precision play a vital role, e.g. structural health monitoring. Therefore, extensive tests against data with relevant characteristics and reliable ground-truth are required to assess and ensure the performance of 3D modelling workflows. This article describes the images collected by a customized unmanned aerial vehicle (UAV) system from an open-pit gravel mine accompanied with additional data that will allow implementing and evaluating any structure-from-motion or photogrammetric approach for sparse or dense 3D reconstruction. This dataset includes total of 158 high-quality images captured with more than 80% endlap and spatial resolution higher than 1.5 cm, the 3D coordinates of 109 ground control points and checkpoints, 2D coordinates of more than 40K corresponding points among the images, a subset of 25 multi-view stereo images selected from an area of approximately 30 m × 40 m within the scene accompanied with a dense point cloud measured by a terrestrial laser scanner.
| S-EPMC6554229 | biostudies-literature
Project description:BackgroundPositron emission tomography (PET) with prostate specific membrane antigen (PSMA) have shown superior performance in detecting metastatic prostate cancers. Relative to [18F]fluorodeoxyglucose ([18F]FDG) PET images, PSMA PET images tend to visualize significantly higher-contrast focal lesions. We aim to evaluate segmentation and reconstruction algorithms in this emerging context. Specifically, Bayesian or maximum a posteriori (MAP) image reconstruction, compared to standard ordered subsets expectation maximization (OSEM) reconstruction, has received significant interest for its potential to reach convergence with minimal noise amplifications. However, few phantom studies have evaluated the quantitative accuracy of such reconstructions for high contrast, small lesions (sub-10 mm) that are typically observed in PSMA images. In this study, we cast 3 mm-16-mm spheres using epoxy resin infused with a long half-life positron emitter (sodium-22; 22Na) to simulate prostate cancer metastasis. The anthropomorphic Probe-IQ phantom, which features a liver, bladder, lungs, and ureters, was used to model relevant anatomy. Dynamic PET acquisitions were acquired and images were reconstructed with OSEM (varying subsets and iterations) and BSREM (varying β parameters), and the effects on lesion quantitation were evaluated.ResultsThe 22Na lesions were scanned against an aqueous solution containing fluorine-18 (18F) as the background. Regions-of-interest were drawn with MIM Software using 40% fixed threshold (40% FT) and a gradient segmentation algorithm (MIM's PET Edge+). Recovery coefficients (RCs) (max, mean, peak, and newly defined "apex"), metabolic tumour volume (MTV), and total tumour uptake (TTU) were calculated for each sphere. SUVpeak and SUVapex had the most consistent RCs for different lesion-to-background ratios and reconstruction parameters. The gradient-based segmentation algorithm was more accurate than 40% FT for determining MTV and TTU, particularly for lesions [Formula: see text] 6 mm in diameter (R2 = 0.979-0.996 vs. R2 = 0.115-0.527, respectively).ConclusionAn anthropomorphic phantom was used to evaluate quantitation for PSMA PET imaging of metastatic prostate cancer lesions. BSREM with β = 200-400 and OSEM with 2-5 iterations resulted in the most accurate and robust measurements of SUVmean, MTV, and TTU for imaging conditions in 18F-PSMA PET/CT images. SUVapex, a hybrid metric of SUVmax and SUVpeak, was proposed for robust, accurate, and segmentation-free quantitation of lesions for PSMA PET.
| S-EPMC8761183 | biostudies-literature
Project description:Iterative reconstruction algorithms are becoming more commonly employed in positron emission tomography (PET) imaging; however, the quantitative accuracy of the reconstructed images still requires validation for various levels of contrast and counting statistics.The authors present an evaluation of the quantitative accuracy of the 3D maximum a posteriori (3D-MAP) image reconstruction algorithm for dynamic PET imaging with comparisons to two of the most widely used reconstruction algorithms: the 2D filtered-backprojection (2D-FBP) and 2D-ordered subsets expectation maximization (2D-OSEM) on the Siemens microPET scanners. The study was performed for various levels of count density encountered in typical dynamic scanning as well as the imaging of cardiac activity concentration in small animal studies on the Focus 120. Specially designed phantoms were used for evaluation of the spatial resolution, image quality, and quantitative accuracy. A normal mouse was employed to evaluate the accuracy of the blood time activity concentration extracted from left ventricle regions of interest (ROIs) within the images as compared to the actual blood activity concentration measured from arterial blood sampling.For MAP reconstructions, the spatial resolution and contrast have been found to reach a stable value after 20 iterations independent of the β values (i.e., hyper parameter which controls the weight of the penalty term) and count density within the frame. The spatial resolution obtained with 3D-MAP reaches values of ∼1.0 mm with a β of 0.01 while the 2D-FBP has value of 1.8 mm and 2D-OSEM has a value of 1.6 mm. It has been observed that the lower the hyper parameter β used in MAP, more iterations are needed to reach the stable noise level (i.e., image roughness). The spatial resolution is improved by using a lower β value at the expense of higher image noise. However, with similar noise level the spatial resolution achieved by 3D-MAP was observed to be better than that by 2D-FBP or 2D-OSEM. Using an image quality phantom containing hot spheres, the estimated activity concentration in the largest sphere has the expected concentration relative to the background area for all the MAP images. The obtained recovery coefficients have been also shown to be almost independent of the count density. 2D-FBP and 2D-OSEM do not perform as well, yielding recovery coefficients lower than those observed with 3D-MAP (approximately 33% lower for the smallest sphere). However, a small positive bias was observed in MAP reconstructed images for frames of very low count density. This bias is present in the uniform area for count density of less than 0.05 × 10(6) counts/ml. For the dynamic mouse study, it was observed that 3D-MAP (even gated at diastole) cannot predict accurately the blood activity concentration due to residual spill-over activity from the myocardium into the left ventricle (approximately 15%). However, 3D-MAP predicts blood activity concentration closer to blood sampling than 2D-FBP.The authors observed that 3D-MAP produces more accurate activity concentration estimates than 2D-FBP or 2D-OSEM at all practical levels of statistics and contrasts due to improved spatial resolution leading to lesser partial volume effect.
| S-EPMC3281972 | biostudies-other
Project description: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.
| S-EPMC6899618 | biostudies-literature
Project description:In positron emission tomography (PET) image reconstruction, the Bayesian framework with various regularization terms has been implemented to constrain the radio tracer distribution. Varying the regularizing weight of a maximum a posteriori (MAP) algorithm specifies a lower bound of the tradeoff between variance and spatial resolution measured from the reconstructed images. The purpose of this paper is to build a patch-based image enhancement scheme to reduce the size of the unachievable region below the bound and thus to quantitatively improve the Bayesian PET imaging. We cast the proposed enhancement as a regression problem which models a highly nonlinear and spatial-varying mapping between the reconstructed image patches and an enhanced image patch. An artificial neural network model named multilayer perceptron (MLP) with backpropagation was used to solve this regression problem through learning from examples. Using the BrainWeb phantoms, we simulated brain PET data at different count levels of different subjects with and without lesions. The MLP was trained using the image patches reconstructed with a MAP algorithm of different regularization parameters for one normal subject at a certain count level. To evaluate the performance of the trained MLP, reconstructed images from other simulations and two patient brain PET imaging data sets were processed. In every testing cases, we demonstrate that the MLP enhancement technique improves the noise and bias tradeoff compared with the MAP reconstruction using different regularizing weights thus decreasing the size of the unachievable region defined by the MAP algorithm in the variance/resolution plane.
| S-EPMC6132251 | biostudies-literature