How to apply noise reduction filter to my DICOM image series

HI all

I am trying to apply noise reduction filter to my DICOM image series (n = 16). let's suppose my image name is (MM.dcm). Can any one help me on how to do this ?

Note that i was recieving this massege every time i want to apply Wiener2 filter:

WIENER2 does not support 3D truecolor images as an input.

The code i used is below:

imvol = dicomread(MM.dcm);
imvol = squeeze(imvol);
imvol = wiener2(imvol,[5 5]);

Answer 1

[Kritika Bansal]

You can possibly apply the 2D Wiener filter on each channel of your 3D DICOM image and see how it works. Probably, something like this may work:

x = wiener2(imvol(:,:,1), [5 5]);
y = wiener2(imvol(:,:,2), [5 5]);
z = wiener2(imvol(:,:,3), [5 5]);
Inew(:,:,1) = x;
Inew(:,:,2) = y;
Inew(:,:,3) = z;

[Muhammad Kamran]

Thanks for this suggestion, it worked for me.