A Krasnoselskii-Mann Algorithm With an Improved EM Preconditioner for PET Image Reconstruction.
Ontology highlight
ABSTRACT:
SUBMITTER: Lin Y
PROVIDER: S-EPMC7528397 | biostudies-literature | 2019 Sep
REPOSITORIES: biostudies-literature
Publications
A Krasnoselskii-Mann Algorithm With an Improved EM Preconditioner for PET Image Reconstruction.
Lin Yizun Y Schmidtlein C Ross CR Li Qia Q Li Si S Xu Yuesheng Y
IEEE transactions on medical imaging 20190219 9
This paper presents a preconditioned Krasnoselskii-Mann (KM) algorithm with an improved EM preconditioner (IEM-PKMA) for higher-order total variation (HOTV) regularized positron emission tomography (PET) image reconstruction. The PET reconstruction problem can be formulated as a three-term convex optimization model consisting of the Kullback-Leibler (KL) fidelity term, a nonsmooth penalty term, and a nonnegative constraint term which is also nonsmooth. We develop an efficient KM algorithm for so ...[more]
Similar Datasets
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:PurposeThe aim of this study was to evaluate the use of a Bayesian penalized likelihood reconstruction algorithm (Q.Clear) for 89Zr-immunoPET image reconstruction and its potential to improve image quality and reduce the administered activity of 89Zr-immunoPET tracers.MethodsEight 89Zr-immunoPET whole-body PET/CT scans from three 89Zr-immunoPET clinical trials were selected for analysis. On average, patients were imaged 6.3 days (range 5.0-8.0 days) after administration of 69 MBq (range 65-76 MBq) of [89Zr]Zr-DFO-daratumumab, [89Zr]Zr-DFO-pertuzumab, or [89Zr]Zr-DFO-trastuzumab. List-mode PET data was retrospectively reconstructed using Q.Clear with incremental β-values from 150 to 7200, as well as standard ordered-subset expectation maximization (OSEM) reconstruction (2-iterations, 16-subsets, a 6.4-mm Gaussian transaxial filter, "heavy" z-axis filtering and all manufacturers' corrections active). Reduced activities were simulated by discarding 50% and 75% of original counts in each list mode stream. All reconstructed PET images were scored for image quality and lesion detectability using a 5-point scale. SUVmax for normal liver and sites of disease and liver signal-to-noise ratio were measured.ResultsQ.Clear reconstructions with β = 3600 provided the highest scores for image quality. Images reconstructed with β-values of 3600 or 5200 using only 50% or 25% of the original counts provided comparable or better image quality scores than standard OSEM reconstruction images using 100% of counts.ConclusionThe Bayesian penalized likelihood reconstruction algorithm Q.Clear improved the quality of 89Zr-immunoPET images. This could be used in future studies to improve image quality and/or decrease the administered activity of 89Zr-immunoPET tracers.
| S-EPMC7815860 | 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:INTRODUCTION:11C-acetate (ACE)-PET/CT is used for staging of high-risk prostate cancer. PET data is reconstructed with iterative algorithms, such as VUEPointHD ViP (VPHD) and VUEPoint HD Sharp IR (SharpIR), the latter with additional resolution recovery. It is expected that the resolution recovery algorithm should render more accurate maximum and mean standardized uptake values (SUVmax and SUVmean) and functional tumor volumes (FTV) than the ordinary OSEM. Performing quantitative analysis, choice of volume-of-interest delineation algorithm (SUV threshold) may influence FTV. Optimizing PET acquisition time is justified if image quality and quantitation do not deteriorate. The aim of this study is to identify the optimal reconstruction algorithm, SUV threshold and acquisition time for ACE-PET/CT. METHODS:ACE-PET/CT data acquired with a General Electric Discovery 690 PET/CT from 16 consecutive high-risk prostate cancer patients was reconstructed with VPHD and SharpIR. Forty pelvic lymph nodes (LNs) and 14 prostate glands were delineated with 42% and estimated threshold. SUVmax, SUVmean, FTV and total lesion uptake were measured. Default acquisition time was four minutes per bed position. In a subset of lesions, acquisition times of one, two and four minutes were evaluated. Structural tumor volumes (STV) of the LNs were measured with CT for correlation with functional volumetric parameters. To validate SUV quantification under different conditions with SharpIR 42%, recovery coefficients (RCs) of SUVmean and FTV were calculated from a phantom with 18F-fluoro-deoxy-glucose (FDG)-filled volumes 0.1-9.2cm3 and signal-to-background (S/B) ratios 4.3-15.9. RESULTS:With SharpIR, SUVmax and SUVmean were higher and FTV lower compared with VPHD, regardless of threshold method, in both prostates and LNs. Total lesion uptake determined with both threshold methods was lower with SharpIR compared with VPHD with both threshold methods, except in subgroup analysis of prostate targets where estimated threshold returned higher values. Longer acquisition times returned higher FTV for both threshold methods, regardless of reconstruction algorithm. The FTV difference was most pronounced with one minute's acquisition per bed position, which also produced visually the highest noise. SUV parameters were unaffected by varying acquisition times. FTV with SharpIR 42% showed the best correspondence with STV. SharpIR 42% gave higher RCs of SUVmean and FTV with increasing phantom size and S/B-ratio, as expected. CONCLUSIONS:Delineation with SharpIR 42% seems to provide the most accurate combined information from SUVmax, SUVmean, FTV and total lesion uptake. Acquisition time may be shortened to two minutes per bed position with preserved image quality.
| S-EPMC6292629 | biostudies-literature
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
Project description:The data presented here provide information about the role of reconstruction parameters on Positron Emission Tomography (PET) image quantification. Multiple phantom measurements in four different Spheres to Background Ratio (SBR) were performed on Biograph 6 TruePoint TrueV PET/CT scanner. PET raw data were reconstructed with/without resolution recovery algorithm using six various iteration x subsets with five different Full-Width Half-Maximum (FWHM) values of Gaussian post-smoothing filter. The Recovery Coefficient (RC) of six spheres using three common Volume of Interest (VOI) methods (max, 3D-50% Isocontour, and peak) were calculated. Moreover, SUVmax, SUVmean, and SUVpeak and volumetric indices, such as metabolic tumor volume (MTV), volume recovery coefficient (VRC), and total lesion glycolysis (TLG) were measured. RCmax, RC50%, and RCpeak as a function of sphere size were plotted in all reconstruction methods considering different SBRs. The data could be noticeable for standardization and optimization of quantitative metrics in PET imaging.
| S-EPMC8237526 | biostudies-literature
Project description:To address the issue of image denoising algorithms blurring image details during the denoising process, we propose an adaptive denoising algorithm for the 3D reconstruction of 2D images. This algorithm takes into account the inherent visual characteristics of human eyes and divides the image into regions based on the entropy value of each region. The background region is subject to threshold denoising, while the target region undergoes processing using an adversarial generative network. This network effectively handles 2D target images with noise and generates a 3D model of the target. The proposed algorithm aims to enhance the noise immunity of 2D images during the 3D reconstruction process and ensure that the constructed 3D target model better preserves the original image's detailed information. Through experimental testing on 2D images and real pedestrian videos contaminated with noise, our algorithm demonstrates stable preservation of image details. The reconstruction effect is evaluated in terms of noise reduction and the fidelity of the 3D model to the original target. The results show an average noise reduction exceeding 95% while effectively retaining most of the target's feature information in the original image. In summary, our proposed adaptive denoising algorithm improves the 3D reconstruction process by preserving image details that are often compromised by conventional denoising techniques. This has significant implications for enhancing image quality and maintaining target information fidelity in 3D models, providing a promising approach for addressing the challenges associated with noise reduction in 2D images during 3D reconstruction.
| S-EPMC10557518 | biostudies-literature
Project description:We investigated the effect of deep learning-based image reconstruction (DLIR) compared to iterative reconstruction on image quality in CT pulmonary angiography (CTPA) for suspected pulmonary embolism (PE). For 220 patients with suspected PE, CTPA studies were reconstructed using filtered back projection (FBP), adaptive statistical iterative reconstruction (ASiR-V 30%, 60% and 90%) and DLIR (low, medium and high strength). Contrast-to-noise ratio (CNR) served as the primary parameter of objective image quality. Subgroup analyses were performed for normal weight, overweight and obese individuals. For patients with confirmed PE (n = 40), we further measured PE-specific CNR. Subjective image quality was assessed independently by two experienced radiologists. CNR was lowest for FBP and enhanced with increasing levels of ASiR-V and, even more with increasing strength of DLIR. High strength DLIR resulted in an additional improvement in CNR by 29-67% compared to ASiR-V 90% (p < 0.05). PE-specific CNR increased by 75% compared to ASiR-V 90% (p < 0.05). Subjective image quality was significantly higher for medium and high strength DLIR compared to all other image reconstructions (p < 0.05). In CT pulmonary angiography, DLIR significantly outperforms iterative reconstruction for increasing objective and subjective image quality. This may allow for further reductions in radiation exposure in suspected PE.
| S-EPMC10827738 | biostudies-literature