Denoising of 3D magnetic resonance images using non-local PCA and Transform-Domain Filter

The Magnetic Resonance Imaging (MRI) technologyused in clinical diagnosis demands high Peak Signal-to-Noise ratio(PSNR) and improved resolution for accurate analysis and treatmentmonitoring. However, MRI data is often corrupted by random noisewhich degrades the quality of Magnetic Resonance (MR) images.Denoising is a paramount challenge as removing noise causesreduction in the fine details of MRI images. We have developed anovel algorithm which employs Principal Component Analysis(PCA) decomposition and Wiener filtering. We have proposed a twostage approach. In first stage, non-local PCA thresholding is appliedon noisy image and second stage uses Wiener filter over this filteredimage. Our algorithm is implemented using MATLAB andperformance is measured via PSNR. The proposed approach hasalso been compared with related state-of-art methods. Moreover, wepresent both qualitative and quantitative results which prove thatproposed algorithm gives superior denoising performance

Doctor of Philosophy

dissertationThe gold standard for evaluation of arterial disease using MR continues to be contrast-enhanced MR angiography (MRA) with gadolinium-based contrast agents (Gd-MRA). There has been a recent resurgence in interest in methods that do not rely on gadolinium for enhancement of blood vessels due to associations Gd-MRA has with nephrogenic systemic fibrosis (NSF) in patients with impaired renal function. The risk due to NSF has been shown to be minimized when selecting the appropriate contrast type and dose. Even though the risk of NSF has been shown to be minimized, demand for noncontrast MRA has continued to rise to reduce examination cost, and improve patient comfort and ability to repeat scans. Several methods have been proposed and used to perform angiography of the aorta and peripheral arteries without the use of gadolinium. These techniques have had limitations in transmit radiofrequency field (B1+) inhomogeneities, acquisition time, and specific hardware requirements, which have stunted the utility of noncontrast enhanced MRA. In this work feasibility of noncontrast (NC) MRA at 3T of the femoral arteries using dielectric padding, and using 3D radial stack of stars and compressed sensing to accelerate acquisitions in the abdomen and thorax were tested. Imaging was performed on 13 subjects in the pelvis and thighs using high permittivity padding, and 11 in the abdomen and 19 in the thorax using 3D radial stack of stars with tiny golden angle using gold standards or previously published techniques. Qualitative scores for each study were determined by radiologists who were blinded to acquisition type. Vessel conspicuity in the thigh and pelvis showed significant increase when high permittivity padding was used in the acquisition. No significant difference in image quality was observed in the abdomen and thorax when using undersampling, except for the descending aorta in thoracic imaging. All image quality scores were determined to be of diagnostic quality. In this work it is shown that NC-MRA can be improved through the use of high permittivity dielectric padding and acquisition time can be decreased through the use of 3D radial stack of stars acquisitions

Improving undersampled MRI reconstruction using non-local means

Journal ArticleObtaining high quality images in MR is desirable not only for accurate visual assessment but also for automatic processing to extract clinically relevant parameters. Filtering-based techniques are extremely useful for reducing artifacts caused due to undersampling of k-space (to reduce scan time). The recently proposed Non-Local Means (NLM) filtering method offers a promising means to denoise images. Compared to most previous approaches, NLM is based on a more realistic model of images, which results in little loss of information while removing the noise. Here we extend the NLM method for MR image reconstruction from undersampled k-space data. The method is applied on T1-weighted images of the breast and T2-weighted anatomical brain images. Results show that NLM offers a promising method that can be used for accelerating MR data acquisitions

A CURE for noisy magnetic resonance images: Chi-square unbiased risk estimation

In this article we derive an unbiased expression for the expected mean-squared error associated with continuously differentiable estimators of the noncentrality parameter of a chi-square random variable. We then consider the task of denoising squared-magnitude magnetic resonance image data, which are well modeled as independent noncentral chi-square random variables on two degrees of freedom. We consider two broad classes of linearly parameterized shrinkage estimators that can be optimized using our risk estimate, one in the general context of undecimated filterbank transforms, and another in the specific case of the unnormalized Haar wavelet transform. The resultant algorithms are computationally tractable and improve upon state-of-the-art methods for both simulated and actual magnetic resonance image data.Comment: 30 double-spaced pages, 11 figures; submitted for publicatio

Postreconstruction filtering of 3D PET images by using weighted higher-order singular value decomposition

Additional file 1. Original 3D PET images data used in this work to generate the results

preliminary clinical evaluation of the ASTRA4D algorithm

Objectives. To propose and evaluate a four-dimensional (4D) algorithm for joint motion elimination and spatiotemporal noise reduction in low-dose dynamic myocardial computed tomography perfusion (CTP). Methods. Thirty patients with suspected or confirmed coronary artery disease were prospectively included und underwent dynamic contrast-enhanced 320-row CTP. The presented deformable image registration method ASTRA4D identifies a low-dimensional linear model of contrast propagation (by principal component analysis, PCA) of the ex-ante temporally smoothed time-intensity curves (by local polynomial regression). Quantitative (standard deviation, signal-to-noise ratio (SNR), temporal variation, volumetric deformation) and qualitative (motion, contrast, contour sharpness; 1, poor; 5, excellent) measures of CTP quality were assessed for the original and motion-compensated volumes (without and with temporal filtering, PCA/ASTRA4D). Following visual myocardial perfusion deficit detection by two readers, diagnostic accuracy was evaluated using 1.5T magnetic resonance (MR) myocardial perfusion imaging as the reference standard in 15 patients. Results. Registration using ASTRA4D was successful in all 30 patients and resulted in comparison with the benchmark PCA in significantly (p<0.001) reduced noise over time (-83%, 178.5 vs 29.9) and spatially (-34%, 21.4 vs 14.1) as well as improved SNR (+47%, 3.6 vs 5.3) and subjective image quality (motion, contrast, contour sharpness: +1.0, +1.0, +0.5). ASTRA4D resulted in significantly improved per-segment sensitivity of 91% (58/64) and similar specificity of 96% (429/446) compared with PCA (52%, 33/64; 98%, 435/446; p=0.011) and the original sequence (45%, 29/64; 98%, 438/446; p=0.003) in the visual detection of perfusion deficits. Conclusions. The proposed functional approach to temporal denoising and morphologic alignment was shown to improve quality metrics and sensitivity of 4D CTP in the detection of myocardial ischemia.Zielsetzung. Die Entwicklung und Bewertung einer Methode zur simultanen Rauschreduktion und Bewegungskorrektur für niedrig dosierte dynamische CT Myokardperfusion. Methoden. Dreißig prospektiv eingeschlossene Patienten mit vermuteter oder bestätigter koronarer Herzkrankheit wurden einer dynamischen CT Myokardperfusionsuntersuchung unterzogen. Die präsentierte Registrierungsmethode ASTRA4D ermittelt ein niedrigdimensionales Modell des Kontrastmittelflusses (mittels einer Hauptkomponentenanalyse, PCA) der vorab zeitlich geglätteten Intensitätskurven (mittels lokaler polynomialer Regression). Quantitative (Standardabweichung, Signal-Rausch-Verhältnis (SNR), zeitliche Schwankung, räumliche Verformung) und qualitative (Bewegung, Kontrast, Kantenschärfe; 1, schlecht; 5, ausgezeichnet) Kennzahlen der unbearbeiteten und bewegungskorrigierten Perfusionsdatensätze (ohne und mit zeitlicher Glättung PCA/ASTRA4D) wurden ermittelt. Nach visueller Beurteilung von myokardialen Perfusionsdefiziten durch zwei Radiologen wurde die diagnostische Genauigkeit im Verhältnis zu 1.5T Magnetresonanztomographie in 15 Patienten ermittelt. Resultate. Bewegungskorrektur mit ASTRA4D war in allen 30 Patienten erfolgreich und resultierte im Vergleich mit der PCA Methode in signifikant (p<0.001) verringerter zeitlicher Schwankung (-83%, 178.5 gegenüber 29.9) und räumlichem Rauschen (-34%, 21.4 gegenüber 14.1) sowie verbesserter SNR (+47%, 3.6 gegenüber 5.3) und subjektiven Qualitätskriterien (Bewegung, Kontrast, Kantenschärfe: +1.0, +1.0, +0.5). ASTRA4D resultierte in signifikant verbesserter segmentweiser Sensitivität 91% (58/64) und ähnlicher Spezifizität 96% (429/446) verglichen mit der PCA Methode (52%, 33/64; 98%, 435/446; p=0.011) und dem unbearbeiteten Perfusionsdatensatz (45%, 29/64; 98%, 438/446; p=0.003) in der visuellen Beurteilung von myokardialen Perfusionsdefiziten. Schlussfolgerungen. Der vorgeschlagene funktionale Ansatz zur simultanen Rauschreduktion und Bewegungskorrektur verbesserte Qualitätskriterien und Sensitivität von dynamischer CT Perfusion in der visuellen Erkennung von Myokardischämie