非线性滤波:中值滤波,双边滤波
1.1 非线性滤波
滤波器是线性的,即两个信号之和的响应是它们各自响应的和(每个像素的输出值是一些输入像素的加权和。线性滤波器易于构造,但是很多情况下,使用邻域像素的非线性滤波会达到更好的效果。比如:在噪声是散粒噪声而不是高斯噪声时,图像偶尔出现很大的值的时候,用高斯滤波器对图像进行模糊的话,噪声像素是不会被去除的,他们只会转换为更为柔和的但仍可见的散粒,此时可以使用中值滤波;
1.2 中值滤波
中值滤波是一种典型的非滤波技术,基本方法就是:用像素点领域灰度值的中值来代替该像素点的灰度值。它在去除脉冲噪声,椒盐噪声的同时又能保留图像的边缘细节。
中值滤波是一种基于排序统计理论的一种能有效抑制噪声的非线性信号处理技术,原理:把数字图像或数字序列中的一点的值用该点的一个邻域中个各值的中值代替,让周围的像素值接近真实值,从而消除孤立的噪声点。这对于斑点噪声和椒盐噪声来说很重要,因为它们不依赖与邻域内那些与典型值相差很大的值。中值滤波器在处理连续图像窗口函数时与线性滤波器的工作方式一样,但滤波过程不再加权。
中值滤波在一定的条件下可以克服常见线性滤波器(最小均方滤波,方框滤波,均值滤波)等带来的图像细节模糊,对滤除脉冲干扰,图像扫描噪声有效,也用于保护边缘信息。
1.3 中值滤波与均值滤波的比较
1.4 双边滤波
是结合图像的空间邻近度和像素值相似的一种折中处理,同时考虑空域信息和灰度相似性,达到保边取噪声的目的,具有简单,非迭代,局部的特点。
1.5 非线性滤波相关函数
(1)中值滤波:medianBlur函数
作用:使用中值滤波器来平滑处理一张图像。对于多通道图像,它对每一个通道都进行单独的处理,并且支持就地操作。
原型:void medianBlur(InputArray src, OutputArray dst, int ksize)
#include "opencv2/opencv.hpp"
#include "opencv2/highgui/highgui.hpp"
#include "opencv2/imgproc/imgproc.hpp"
#include
using namespace std;
using namespace cv;
int main()
{
Mat image = imread("../1.jpg");
namedWindow("原图");
namedWindow("中值滤波图");
imshow("原图", image);
Mat out;
medianBlur(image, out, 7);
imshow("中值滤波图",out );
while (char(waitKey(1)) != 'q') {}
}
(2)双边滤波:bilateraFilter函数
作用:用双边滤波器来模糊一张图像;
原型:void bilateralFilter(InputArray src, OutputArray dst, int d, double sigmaColor, double sigmaSpace , int borderType=BORDER_BEFAULT)
#include "opencv2/opencv.hpp"
#include "opencv2/highgui/highgui.hpp"
#include "opencv2/imgproc/imgproc.hpp"
#include
using namespace std;
using namespace cv;
int main()
{
Mat image = imread("../1.jpg");
namedWindow("原图");
namedWindow("双边滤波图");
imshow("原图", image);
Mat out;
bilateralFilter(image, out, 25, 25*2, 25/2);
imshow("双边滤波图",out );
while (char(waitKey(1)) != 'q') {}
}
#include "opencv2/opencv.hpp"
#include "opencv2/highgui/highgui.hpp"
#include "opencv2/imgproc/imgproc.hpp"
#include
using namespace std;
using namespace cv;
Mat image, dst_medianImage, dst_bilateraImage;
int g_nMedianBlurValue = 10;
int g_nBilaterValue = 10;
static void On_MedianBlur(int, void*);
static void On_BilaterBlur(int , void*);
int main()
{
image = imread("../1.jpg");
if(!image.data) {
printf("获取图像路径错误");
return false;
}
dst_medianImage = image.clone();
dst_bilateraImage = image.clone();
namedWindow("原图");
imshow("原图", image);
namedWindow("中值滤波");
namedWindow("双边滤波");
createTrackbar("参数值:", "中值滤波", &g_nMedianBlurValue, 50, On_MedianBlur);
createTrackbar("参数值:", "双边滤波", &g_nBilaterValue, 50, On_BilaterBlur);
On_MedianBlur(g_nMedianBlurValue, 0);
On_BilaterBlur(g_nBilaterValue, 0);
while (char (waitKey(1)) != 'q') {}
return 0;
}
static void On_MedianBlur(int , void*)
{
medianBlur(image, dst_medianImage, g_nMedianBlurValue*2 + 1);
imshow("中值滤波", dst_medianImage);
}
static void On_BilaterBlur(int , void*)
{
bilateralFilter(image, dst_bilateraImage, g_nBilaterValue, g_nBilaterValue*2, g_nBilaterValue/2);
imshow("双边滤波", dst_bilateraImage);
}