Premium
A simple method for estimating the noise level in a signal region of an MR image
Author(s) -
Steckner Michael C.
Publication year - 2010
Publication title -
medical physics
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.473
H-Index - 180
eISSN - 2473-4209
pISSN - 0094-2405
DOI - 10.1118/1.3480511
Subject(s) - pixel , artifact (error) , diagonal , noise (video) , imaging phantom , computer science , algorithm , signal (programming language) , image (mathematics) , simple (philosophy) , set (abstract data type) , range (aeronautics) , artificial intelligence , mathematics , computer vision , optics , physics , philosophy , geometry , materials science , epistemology , composite material , programming language
Purpose: A simple extension to the NEMA MS‐1 “difference of neighboring pixels” SNR method is presented, which can accurately determine the noise level within a signal region over a wide range of noise levels, image nonuniformities, and artifact levels, as demonstrated by simple simulations and experimental phantom images. Methods: The new method computes difference of neighboring pixels in the read, phase, and diagonal directions. The variance of these three sets of pixel differences appear to contain the simple sum of the underlying variance of noise and any additional component unique to the read and phase directions, respectively, while the diagonal set of pixel differences contains all three components. By solving a set of three equations with three unknowns, it is possible to extract the components and isolate the desired noise variance term. Results: Simulations show that the technique produces good results even when various artifact mechanisms present singly or jointly. Experimental results also demonstrate the technique works well but, depending on the severity of the artifact, cannot be guaranteed to always produce accurate results. Conclusions: Simulations and experimental results show the method to be accurate and robust. This method is applicable to multichannel receiver images, but not parallel reconstructed images.