借助于opencv的数据类型Mat简单实现均值滤波,该代码并未考虑边缘影响
#include <iostream>
#include "vector"
#include "algorithm"
#include "opencv2/highgui.hpp"
#include "opencv2/core/core.hpp"
#include "opencv2/opencv.hpp"
int mymean(int n1,int n2,int n3,int n4,int n5,int n6,int n7,int n8,int n9)
{
int sum=0;
int arr[9]={n1,n2,n3,n4,n5,n6,n7,n8,n9};
for(auto elem:arr)
{
sum+=elem;
}
return (int)sum/9;
}
void mymean_Blur(cv::Mat& src,cv::Mat& dst)
{
cv::Mat img;
src.copyTo(img);
cv::Mat dst_(img.size(),img.type());
int rows=img.rows,cols=img.cols;
for(int r=1;r<rows-1;++r)
{
for(int c=1;c<cols-1;++c)
{
dst_.at<uchar>(r,c)= mymean(img.at<uchar>(r-1,c-1),img.at<uchar>(r-1,c),img.at<uchar>(r-1,c+1),img.at<uchar>(r,c-1),img.at<uchar>(r,c),img.at<uchar>(r,c+1),img.at<uchar>(r+1,c-1), img.at<uchar>(r+1,c),img.at<uchar>(r+1,c+1));
}
}
dst_.copyTo(dst);
}
void salt(cv::Mat& src,int num)
{
cv::RNG rng(12345);//使用种子seed产生一个64位随机整数
//rng.uniform(a,b);//uniform产生一个[a,b)范围内均匀分布的随机数
//rng.gaussian(a);//gaussian() 产生一个均值为0,标准差为a的随机数(高斯分布)
if(!src.data)
{
return;
}
int rows=src.rows,cols=src.cols;
for(int i=0;i<num;++i)
{
int c=rng.uniform(0,cols);//uniform产生一个均匀分布的随机数
int r=rng.uniform(0,rows);
src.at<cv::Vec3b>(r,c)=cv::Vec3b(255,255,255);
src.at<uchar>(r,c);
}
// cv::imshow("salt",src);
}
int main()
{
std::string imgDir="/home/.../images/anni.jpg";
cv::Mat img=cv::imread(imgDir);
if(img.empty())
{
std::cout<<"file fail to load!"<<std::endl;
return -1;
}
salt(img,10000);
cv::resize(img,img,cv::Size(640,480),0,0,cv::INTER_LINEAR);
cv::imshow("Raw",img);
cv::cvtColor(img,img,cv::COLOR_BGR2GRAY);
cv::Mat dst;
mymean_Blur(img,dst);
cv::imshow("Test",dst);
cv::waitKey(0);
std::cout << "Hello, World!" << std::endl;
return 0;
}
运行结果:
参考:
https://blog.csdn.net/weixin_46196863/article/details/112255510 椒盐噪声
https://blog.csdn.net/qq_33485434/article/details/78980587 随机数生成
https://blog.csdn.net/weixin_37720172/article/details/72627543 均值滤波