File Exchange

image thumbnail

LMMSE filter for Rician MRI data

version 1.2.0.0 (3 KB) by SANTIAGO AJA-FERNANDEZ
A filtering scheme for denoising MR data assuming an underlying Rician model.
4.6
5 Ratings

4 Downloads

Updated 12 Jun 2012

View Version History

View License

The linear minimum mean square error (LMMSE) estimator estimates the signal out of a noisy MR adquisition. The method automatically estimates the noise level from background. It is based on a Rician noise assumption, which is valid for single coil adquisitions (and sometimes in SENSE).

The whole method is described in

S Aja-Fernández, M. Niethammer, M. Kubicki, M. E. Shenton, C.F. Westin, "Restoration of DWI data using a Rician LMMSE estimator". IEEE Tr. on Medical Imag. Vol. 27, No. 10, Oct. 2008, pp. 1389-1403.

Cite As

SANTIAGO AJA-FERNANDEZ (2021). LMMSE filter for Rician MRI data (https://www.mathworks.com/matlabcentral/fileexchange/36741-lmmse-filter-for-rician-mri-data), MATLAB Central File Exchange. Retrieved .

Comments and Ratings (10)

Roger Mullins

MATLAB Version: 9.2.0.556344 (R2017a)
Array dimensions must match for binary array op.

Error in MRI_lmmse (line 149)
I_est=sqrt(Squ-2.*sigma2+K1.*(Im.^2-Squ));

NUR LYANA SHAHFIQA ALBASHAH

Not enough input arguments.

Error in mean (line 12)
a=find(image>thr);%m0
How to fix this error?

Error in MRI_lmmse>moda (line 238)
M1=mean(u(:));

Error in MRI_lmmse (line 103)
sigma2=(prod(Ws)/(prod(Ws)-1)).*(moda(En,1000)./2);

max Liu

I would like to know whether this method can be successfully deal with the T1 weighted MRI data

afshin shoeibi

thanks

majid sadeghi

Thanks

Ram Singh

very good stuff

Jorge Rubio Álvarez

Christian

very useful! thank you!

SANTIAGO AJA-FERNANDEZ

The method now requires to choose one of the existing noise estimation methods as an input: (See help)

I_est=MRI_lmmse(I,[7,7],'mode2N');

However, I have just realized that the function moda.m is not in the toolbox. I upload now a new version with this fixed, and a default option

Alex

I get the error:

I_est = MRI_lmmse(img,[10, 10]);
Undefined function or variable "sigma2".

Error in MRI_lmmse (line 137)
K1=1+(4.*sigma2^2-4.*sigma2.*Squ)./(Qua-Squ.^2);

Did I input something incorrectly?

MATLAB Release Compatibility
Created with R11
Compatible with any release
Platform Compatibility
Windows macOS Linux

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!