how can a wiener filter be designed for reduction of noise on images; the output should be the one with reduced noise on image and the input being a noisy image with gaussian noise

What I have done is......I took 2 images....one image corrupted by Gaussian noise and the other is the original image.....and then I calculated the mean square error and peak signal to noise ratio.....these values are produced as output.... But I require the output as an image with reduced noise.....the solution what I found is taking some pixels from the noisy image and calculating the MSE and PSNR.....and replacing these pixels with their values so that the particular pixel values make some difference in the final image.....so can anyone of you help me with the code in MATLAB

Answer 1

[Image Analyst]

But the code is so, so very simple, even for complete beginners - it's like 3 lines of code: a call to each function. Did you look up the functions in the help and try any of the example code? Look up wiener2(), immse(), and psnr() and try the example code. If you really can't figure it out, come back here with your attempt at code.

[Image Analyst]

Pyla's "Answer" moved here since it's a reply to me rather than an answer to the original question at the top of the page:

can you give me the code without using any of the functions directly available in MATLAB.......that is according to its algorithm

[Image Analyst]

Attached is a function for psnr. I think it uses MSE inside. Why do you not want to use any of the built-in functions that have been tested and validated?

[Pyla Sahitya]

it is because my project is designed in such a way that I don't want to use the commands directly available in MATLAB... so can you please help me out....

[Image Analyst]

I had that psnr function before they included it in MATLAB, so that's help. Did you not notice the psnr.m file that was attached? I don't really want to code up my own wiener2 function for you so I won't. You can do that yourself with fft2() if you want, though I warn you, fft2() is a built in function "directly available in MATLAB" so that may mean you don't want to use it. You can program up your own fft2() function if you want. Please read this link: http://www.mathworks.com/matlabcentral/answers/38787-what-can-be-programmed-without-any-built-in-functions

[Jasper Schilling]

Pyla, were you able to ever figure out how to code the Wiener2 function without built in functions?