function s = WienerFilterImage (x, h, k)

H=fft2(h);
H2=conj(H).*H;
s = real(ifft2((H2./((H2+k).*(H+eps))).*fft2(x)));