A simple noniterative principal component technique for rapid noise reduction in parallel MR images.
Ontology highlight
ABSTRACT:
SUBMITTER: Patel AS
PROVIDER: S-EPMC3170692 | biostudies-literature | 2012 Jan
REPOSITORIES: biostudies-literature
Publications
A simple noniterative principal component technique for rapid noise reduction in parallel MR images.
Patel Anand S AS Duan Qi Q Robson Philip M PM McKenzie Charles A CA Sodickson Daniel K DK
NMR in biomedicine 20110504 1
The utilization of parallel imaging permits increased MR acquisition speed and efficiency; however, parallel MRI usually leads to a deterioration in the signal-to-noise ratio when compared with otherwise equivalent unaccelerated acquisitions. At high accelerations, the parallel image reconstruction matrix tends to become dominated by one principal component. This has been utilized to enable substantial reductions in g-factor-related noise. A previously published technique achieved noise reductio ...[more]
Similar Datasets
Project description:Non-destructive measurement of acceleration-induced displacement fields within a closed object is a fundamental challenge. Inferences of how the brain deforms following skull impact have thus relied largely on indirect estimates and course-resolution cadaver studies. We developed a magnetic resonance technique to quantitatively identify the modes of displacement of an accelerating soft object relative to an object enclosing it, and applied it to study acceleration-induced brain deformation in human volunteers. We show that, contrary to the prevailing hypotheses of the field, the dominant mode of interaction between the brain and skull in mild head acceleration is one of sliding arrested by meninges.
| S-EPMC3136495 | biostudies-literature
Project description:PurposeBrief bursts of RF noise during MR data acquisition ("k-space spikes") cause disruptive image artifacts, manifesting as stripes overlaid on the image. RF noise is often related to hardware problems, including vibrations during gradient-heavy sequences, such as diffusion-weighted imaging. In this study, we present an application of the Robust Principal Component Analysis (RPCA) algorithm to remove spike noise from k-space.MethodsCorrupted k-space matrices were decomposed into their low-rank and sparse components using the RPCA algorithm, such that spikes were contained within the sparse component and artifact-free k-space data remained in the low-rank component. Automated center refilling was applied to keep the peaked central cluster of k-space from misclassification in the sparse component.ResultsThis algorithm was demonstrated to effectively remove k-space spikes from four data types under conditions generating spikes: (i) mouse heart T1 mapping, (ii) mouse heart cine imaging, (iii) human kidney diffusion tensor imaging (DTI) data, and (iv) human brain DTI data. Myocardial T1 values changed by 86.1 ± 171 ms following despiking, and fractional anisotropy values were recovered following despiking of DTI data.ConclusionThe RPCA despiking algorithm will be a valuable postprocessing method for retrospectively removing stripe artifacts without affecting the underlying signal of interest. Magn Reson Med 75:2517-2525, 2016. © 2015 The Authors. Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.
| S-EPMC4720596 | biostudies-literature
Project description:Despite the utility of tumour characterisation using quantitative parameter maps from multi-b-value diffusion-weighted MRI (DWI), clinicians often prefer the use of the image with highest diffusion-weighting (b-value), for instance for defining regions of interest (ROIs). However, these images are typically degraded by noise, as they do not utilize the information from the full acquisition. We present a principal component analysis (PCA) approach for model-free denoising of DWI data. PCA-denoising was compared to synthetic MRI, where a diffusion model is fitted for each voxel and a denoised image at a given b-value is generated from the model fit. A quantitative comparison of systematic and random errors was performed on data simulated using several diffusion models (mono-exponential, bi-exponential, stretched-exponential and kurtosis). A qualitative visual comparison was also performed for in vivo images in six healthy volunteers and three pancreatic cancer patients. In simulations, the reduction in random errors from PCA-denoising was substantial (up to 55%) and similar to synthetic MRI (up to 53%). Model-based synthetic MRI denoising resulted in substantial (up to 29% of signal) systematic errors, whereas PCA-denoising was able to denoise without introducing systematic errors (less than 2%). In vivo, the signal-to-noise ratio (SNR) and sharpness of PCA-denoised images were superior to synthetic MRI, resulting in clearer tumour boundaries. In the presence of motion, PCA-denoising did not cause image blurring, unlike image averaging or synthetic MRI. Multi-b-value MRI can be denoised model-free with our PCA-denoising strategy that reduces noise to a level similar to synthetic MRI, but without introducing systematic errors associated with the synthetic MRI method.
| S-EPMC7655121 | biostudies-literature
Project description:We have previously developed a retrospective 4D-MRI technique using body area as the respiratory surrogate, but generally, the reconstructed 4D MR images suffer from severe or mild artifacts mainly caused by irregular motion during image acquisition. Those image artifacts may potentially affect the accuracy of tumor target delineation or the shape representation of surrounding nontarget tissues and organs. So the purpose of this study is to propose an approach employing principal component analysis (PCA), combined with a linear polynomial fitting model, to remodel the displacement vector fields (DVFs) obtained from deformable image registration (DIR), with the main goal of reducing the motion artifacts in 4D MR images. Seven patients with hepatocellular carcinoma (2/7) or liver metastases (5/7) in the liver, as well as a patient with non-small cell lung cancer (NSCLC), were enrolled in an IRB-approved prospective study. Both CT and MR simulations were performed for each patient for treatment planning. Multiple-slice, multiple-phase, cine-MRI images were acquired in the axial plane for 4D-MRI reconstruction. Single-slice 2D cine-MR images were acquired across the center of the tumor in axial, coronal, and sagittal planes. For a 4D MR image dataset, the DVFs in three orthogonal direction (inferior–superior (SI), anterior–posterior (AP), and medial–lateral (ML)) relative to a specific reference phase were calculated using an in-house DIR algorithm. The DVFs were preprocessed in three temporal and spatial dimensions using a polynomial fitting model, with the goal of correcting the potential registration errors introduced by three-dimensional DIR. Then PCA was used to decompose each fitted DVF into a linear combination of three principal motion bases whose spanned subspaces combined with their projections had been validated to be sufficient to represent the regular respiratory motion. By wrapping the reference MR image using the remodeled DVFs, 'synthetic' MR images with reduced motion artifacts were generated at selected phase. Tumor motion trajectories derived from cine-MRI, 4D CT, original 4D MRI, and 'synthetic' 4D MRI were analyzed in the SI, AP, and ML directions, respectively. Their correlation coefficient (CC) and difference (D) in motion amplitude were calculated for comparison. Of all the patients, the means and standard deviations (SDs) of CC comparing 'synthetic' 4D MRI and cine-MRI were 0.98 ± 0.01, 0.98 ± 0.01, and 0.99 ± 0.01 in SI, AP, and ML directions, respectively. The mean ± SD Ds were 0.59 ± 0.09 mm, 0.29± 0.10 mm, and 0.15 ± 0.05 mm in SI, AP and ML directions, respectively. The means and SDs of CC comparing 'synthetic' 4D MRI and 4D CT were 0.96 ± 0.01, 0.95± 0.01, and 0.95 ± 0.01 in SI, AP, and ML directions, respectively. The mean ± SD Ds were 0.76 ± 0.20 mm, 0.33 ± 0.14 mm, and 0.19± 0.07 mm in SI, AP, and ML directions, respectively. The means and SDs of CC comparing 'synthetic' 4D MRI and original 4D MRI were 0.98 ± 0.01, 0.98± 0.01, and 0.97± 0.01 in SI, AP, and ML directions, respectively. The mean ± SD Ds were 0.58 ± 0.10 mm, 0.30 ± 0.09mm, and 0.17 ± 0.04 mm in SI, AP, and ML directions, respectively. In this study we have proposed an approach employing PCA combined with a linear polynomial fitting model to capture the regular respiratory motion from a 4D MR image dataset. And its potential usefulness in reducing motion artifacts and improving image quality has been demonstrated by the preliminary results in oncological patients.
| S-EPMC5690092 | biostudies-literature
Project description:STEM XEDS spectrum images can be drastically denoised by application of the principal component analysis (PCA). This paper looks inside the PCA workflow step by step on an example of a complex semiconductor structure consisting of a number of different phases. Typical problems distorting the principal components decomposition are highlighted and solutions for the successful PCA are described. Particular attention is paid to the optimal truncation of principal components in the course of reconstructing denoised data. A novel accurate and robust method, which overperforms the existing truncation methods is suggested for the first time and described in details.
| S-EPMC6456488 | biostudies-literature
Project description:BACKGROUND:In this paper we apply the principal-component analysis filter (Hotelling filter) to reduce noise from dynamic positron-emission tomography (PET) patient data, for a number of different radio-tracer molecules. We furthermore show how preprocessing images with this filter improves parametric images created from such dynamic sequence.We use zero-mean unit variance normalization, prior to performing a Hotelling filter on the slices of a dynamic time-series. The Scree-plot technique was used to determine which principal components to be rejected in the filter process. This filter was applied to [11C]-acetate on heart and head-neck tumors, [18F]-FDG on liver tumors and brain, and [11C]-Raclopride on brain. Simulations of blood and tissue regions with noise properties matched to real PET data, was used to analyze how quantitation and resolution is affected by the Hotelling filter. Summing varying parts of a 90-frame [18F]-FDG brain scan, we created 9-frame dynamic scans with image statistics comparable to 20 MBq, 60 MBq and 200 MBq injected activity. Hotelling filter performed on slices (2D) and on volumes (3D) were compared. RESULTS:The 2D Hotelling filter reduces noise in the tissue uptake drastically, so that it becomes simple to manually pick out regions-of-interest from noisy data. 2D Hotelling filter introduces less bias than 3D Hotelling filter in focal Raclopride uptake. Simulations show that the Hotelling filter is sensitive to typical blood peak in PET prior to tissue uptake have commenced, introducing a negative bias in early tissue uptake. Quantitation on real dynamic data is reliable. Two examples clearly show that pre-filtering the dynamic sequence with the Hotelling filter prior to Patlak-slope calculations gives clearly improved parametric image quality. We also show that a dramatic dose reduction can be achieved for Patlak slope images without changing image quality or quantitation. CONCLUSIONS:The 2D Hotelling-filtering of dynamic PET data is a computer-efficient method that gives visually improved differentiation of different tissues, which we have observed improve manual or automated region-of-interest delineation of dynamic data. Parametric Patlak images on Hotelling-filtered data display improved clarity, compared to non-filtered Patlak slope images without measurable loss of quantitation, and allow a dramatic decrease in patient injected dose.
| S-EPMC3636030 | biostudies-literature
Project description:Quantification of the angular orientation distribution of fibrous tissue structures in scientific images benefits from the Fourier image analysis to obtain quantitative information. Measurement uncertainties represent a major challenge and need to be considered by propagating them in order to determine an adaptive anisotropic Fourier filter. Our adaptive filter method (AF) is based on the maximum relative uncertainty ?cut of the power spectrum as well as a weighted radial sum with weighting factor ?. We use a Monte-Carlo simulation to obtain realistic greyscale images that include defined variations in fiber thickness, length, and angular dispersion as well as variations in noise. From this simulation the best agreement between predefined and derived angular orientation distribution is found for evaluation parameters ?cut = 2.1% and ? = 1.5. The resulting cumulative orientation distribution was modeled by a sigmoid function to obtain the mean angle and the fiber dispersion. A comparison to a state-of-the-art band-pass method revealed that the AF method is more suitable for the application on greyscale fiber images, since the error of the fiber dispersion significantly decreased from (33.9 ± 26.5)% to (13.2 ± 12.7)%. Both methods were found to accurately quantify the mean fiber orientation with an error of (1.9 ± 1.5)° and (2.3 ± 2.1)° in case of the AF and the band-pass method, respectively. We demonstrate that the AF method is able to accurately quantify the fiber orientation distribution in in vivo second-harmonic generation images of dermal collagen with a mean fiber orientation error of (6.0 ± 4.0)° and a dispersion error of (9.3 ± 12.1)%.
| S-EPMC6964846 | biostudies-literature
Project description:"Principal Component Analysis" (PCA) is an established linear technique for dimensionality reduction. It performs an orthonormal transformation to replace possibly correlated variables with a smaller set of linearly independent variables, the so-called principal components, which capture a large portion of the data variance. The problem of finding the optimal number of principal components has been widely studied for offline PCA. However, when working with streaming data, the optimal number changes continuously. This requires to update both the principal components and the dimensionality in every timestep. While the continuous update of the principal components is widely studied, the available algorithms for dimensionality adjustment are limited to an increment of one in neural network-based and incremental PCA. Therefore, existing approaches cannot account for abrupt changes in the presented data. The contribution of this work is to enable in neural network-based PCA the continuous dimensionality adjustment by an arbitrary number without the necessity to learn all principal components. A novel algorithm is presented that utilizes several PCA characteristics to adaptivly update the optimal number of principal components for neural network-based PCA. A precise estimation of the required dimensionality reduces the computational effort while ensuring that the desired amount of variance is kept. The computational complexity of the proposed algorithm is investigated and it is benchmarked in an experimental study against other neural network-based and incremental PCA approaches where it produces highly competitive results.
| S-EPMC8009402 | biostudies-literature
Project description:Speckle noise limits the diagnostic capabilities of optical coherence tomography (OCT) images, causing both a reduction in contrast and a less accurate assessment of the microstructural morphology of the tissue. We present a speckle-noise reduction method for OCT volumes that exploits the advantages of adaptive-noise wavelet thresholding with a wavelet compounding method applied to several frames acquired from consecutive positions. The method takes advantage of the wavelet representation of the speckle statistics, calculated properly from a homogeneous sample or a region of the noisy volume. The proposed method was first compared quantitatively with different state-of-the-art approaches by being applied to three different clinical dermatological OCT volumes with three different OCT settings. The method was also applied to a public retinal spectral-domain OCT dataset to demonstrate its applicability to different imaging modalities. The results based on four different metrics demonstrate that the proposed method achieved the best performance among the tested techniques in suppressing noise and preserving structural information. The proposed OCT denoising technique has the potential to adapt to different image OCT settings and noise environments and to improve image quality prior to clinical diagnosis based on visual assessment.
| S-EPMC8211087 | biostudies-literature
Project description:Although principal component analysis is frequently used in multivariate/ analysis, it has disadvantages when applied to experimental or diagnostic data. First, the identified principal components have poor generality; since the size and directions of the components are dependent on the particular data set, the components are valid only within the set. Second, the method is sensitive to experimental noise and bias between sample groups, since it cannot reflect the design of experiments; rather, it estimates the same weight and independence of all the samples in the matrix. Third, the resulting components are often difficult to interpret. To address these issues, several options were introduced to the methodology. The resulting components were scaled to unify their size unit. Also, the principal axes were identified using training data sets and shared among experiments. This training data reflects the design of experiments, and its preparation allows noise to be reduced and group bias to be removed. The effects of these options were observed in microarray experiments, and showed an improvement in the separation of groups and robustness to noise. Additionally, unknown samples were appropriately classified using pre-arranged axes, and principal axes well reflected the characteristics of groups in the experiments. This SuperSeries is composed of the SubSeries listed below.
2011-08-15 | GSE31375 | GEO