非局部均值滤波步骤
1. 指定邻域窗口半径d,搜索窗口半径D,高斯函数平滑参数h。
2. 扩展图像边界,保证邻域窗口不越界。
3. 在扩展图像中取邻域窗口W1。
4. 限制搜索窗口的边界,保证可移动的邻域窗口W2不越界。
5. 邻域窗口W2在搜索窗口内滑动,当W2与W1重合时跳过。
6. 在W2滑动过程中使用下述公式计算权值:
其中V(x)和V(y)分别是以x,y为中心的邻域矩阵。
邻域间的距离 与归一化系数Z(x)可用下式计算:
7. 当W1中心像素值遍历到搜索窗口最后一个像素点处时,计算移动的邻域窗口W2内像素的加权和。
8. 将步骤7中的加权和除以Z(x)进行归一化,并替换邻域窗口W1的中心像素值。
9. 移动邻域窗口W1,重复步骤4,直至邻域窗口W1不能再移动。
此外,该程序还要求实现以下功能
1. 显示原图像、移除像素后的图像和恢复出的图像。
2. 求该去噪算法的MSE值。