双边滤波C++实现

不同于中值滤波和均值滤波，双边滤波同时使用了空域权重和值域权重，在滤波的同时能够保持边缘信息不被平滑掉。
确定了滤波核大小后空域权重即确定，越靠近核的边缘权重越低。
值域权重取决于周围像素与目标点像素的差的绝对值，此绝对值越大则周围点的权重越低，目标点受它的影响越小，由此保留了图像边缘区域。
设空域权重为ws, 值域权重为wr（代码中为wc）。
计算公式：（图片来自于csdn某博主）

最终权重为ws*wr（wc），所有权重计算完后需要归一化到和为1，即各个权重值除以权重之和。
下图为某次运行的权重矩阵（尚未归一化）：
C++实现双边滤波代码（下图主程序中同样参数使用opencv双边滤波函数实现为3秒，自己搓的实现为8秒）：
滤波前

滤波后

void Bilateral_Filter(Mat src, Mat& dst, int ksize, double sita_r, double sita_s)
{
   
	int h = src.rows;
	int w = src.cols;
	int k = ksize / 2;
	dst = Mat::zeros(h, w, CV_32FC3);
	Mat ws_mat(ksize, ksize, CV_32FC1);//空域权重矩阵
	for (int i = 0; i < ksize; i++)
	{
   
		for (int j = 0; j < ksize; j++)
		{
   
			ws_mat.at<float>(i, j) = exp(-1.0 * ((i - k) * (i - k) + (