维纳滤波实现

参考链接：Matlab Wiener2函数

一、算法原理及公式：

二、算法实现：

步骤一：计算局部均值图localMean与局部方差图localVar，可采用积分图加速;

步骤二：设noiseVal为噪声方差，若未给出，则参考Matlab wiener2实现，以局部方差图localVar的均值代替；

步骤三：根据公式计算结果:

令：I为源灰度图

f = I - localMean;

err = localVar - noiseVal;

err = max(err，0);

localVar = max(localVar，noise);

f = f ./ localVar;

f = f .* g;

f = f + localMean;

[1]Lim, Jae S., Two-Dimensional Signal and Image Processing, Englewood Cliffs, NJ, Prentice Hall, 1990, p. 548, equations 9.26, 9.27, and 9.29.