opencv 高斯滤波自己实现

最新推荐文章于 2024-08-03 14:54:02 发布

zhzcc

最新推荐文章于 2024-08-03 14:54:02 发布

阅读量1.2k

点赞数 1

分类专栏： opencv 文章标签： opencv

本文链接：https://blog.csdn.net/zhzcc/article/details/70146709

版权

opencv 专栏收录该内容

4 篇文章 0 订阅

订阅专栏

本文展示了如何在C++中利用OpenCV库自定义实现高斯滤波器。通过`CreateMuban`函数创建高斯核，然后在`operateimage`函数中应用该核对图像进行滤波操作。实验中分别使用了不同标准差的高斯核对同一张图片进行处理，以观察不同参数下滤波效果的差异。

摘要由CSDN通过智能技术生成

#include <iostream>
#include <opencv2/opencv.hpp>
using namespace std;
using namespace cv;
double* CreateMuban(int iSize,double sigma)
{
double *gauss=new double[iSize];
int radius=(iSize-1)/2;
double MySigma=2*sigma*sigma;
double value=0;
for (int i = 0; i <iSize; i++)
{
gauss[i]=exp(-(i-radius)*(i-radius)/MySigma);
value=value+gauss[i];
}
for (int i = 0; i < iSize; i++)
{
gauss[i]=gauss[i]/value;
}
return gauss;
}
Mat operateimage(Mat _img,double*Muban,int iSize)
{
Mat img=_img;
Mat img1=_img;
Mat img2=_img;
int radius=(iSize-1)/2;
int r=0;
int g=0;
int b=0;
for (int i = (iSize+1)/2; i < img.rows-(iSize+1)/2; i++)
{
for (int j = (iSize+1)/2; j <img.cols-(iSize+1)/2; j++)
{
for (int k = 0; k < iSize; k++)
{
b=b+(img.at<Vec3b>(i,j-radius+k)[0])*Muban[k];
g=g+(img.at<Vec3b>(i,j-radius+k)[1])*Muban[k];
r=r+(img.at<Vec3b>(i,j-radius+k)[2])*Muban[k];
}
img1.at<Vec3b>(i,j)=Vec3b(b,g,r);
b=0;
g=0;
r=0;
}
}
for (int i = (iSize+1)/2; i < img1.cols-(iSize+1)/2; i++)
{
for (int j = (iSize+1)/2; j <img1.rows-(iSize+1)/2; j++)
{
for (int k = 0; k < iSize; k++)
{
b=b+(img1.at<Vec3b>(j-radius+k,i)[0])*Muban[k];
g=g+(img1.at<Vec3b>(j-radius+k,i)[1])*Muban[k];
r=r+(img1.at<Vec3b>(j-radius+k,i)[2])*Muban[k];
}
img2.at<Vec3b>(j,i)=Vec3b(b,g,r);
b=0;
g=0;
r=0;
}
}
return img2;
}
int main()
{
double* Muban=CreateMuban(3,0.01);
double* Muban1=CreateMuban(3,1);
Mat image=imread("Rpic0.jpg",1);
if (!image.data)
{
return -1;
}
Mat a;
image.copyTo(a);
Mat b;
image.copyTo(b);
Mat image1=operateimage(a,Muban,3);
Mat image2=operateimage(b,Muban1,3);
imshow("Guass",image1);
imshow("Guass1",image2);
imshow("原图",image);
waitKey(0);
return 0;
}