数据图像实验四：图像滤波c++实现

实验4.1：高斯滤波

实现图像的高斯滤波：
通过调整高斯函数的标准差(sigma)来控制平滑程度；
滤波窗口大小取为[6*sigma-1]/2*2+1，[.]表示取整；
利用二维高斯函数的行列可分离性进行加速；
先对每行进行一维高斯滤波，再对结果的每列进行同样的一维高斯滤波；

实验4.2 快速均值滤波

实现图像的均值滤波
滤波窗口大小通过参数来指定：
采用积分图进行加速，实现与滤波窗口大小无关的效率；

4.1代码：

#include<opencv2/core/core.hpp>
#include<opencv2/highgui/highgui.hpp>
#include "math.h"
using namespace cv;
using namespace std;


void Gaussian(Mat img, Mat output,double sigma)
{
    img.copyTo(output);
    int size = sigma * 6 - 1;
    double *weight = new double[size];
    double sum = 0;
    for (int i = 0; i < size; i++)
    {
        weight[i] = exp(((i - size / 2) ^ 2) / -2 / pow(sigma, 2));
        sum = sum + weight[i];
    }
    for (int i = 0; i < size; i++)
    {
        weight[i] = weight[i]/sum;
    }
    for (int num = 0; num < img.channels(); num++)
    {
        for (int x = 0; x < img.rows; x++)
        {
            for (int y = 0; y < img.cols; y++)
            {
                double temp = 0;
                for (int i = 0; i <size; i++)
                {
                    int j = y + i - size / 2;
                    if (j < 0 )
                    {
                        j = -j;
                    }
                    if (j >= img.cols)
                    {
                        j = 2 * img.cols - 1 - j;
                    }
                    temp = temp + img.ptr(x, j)[num] * weight[i];
                }
                output.ptr(x, y)[num] = temp;
            }
        }
    }
    for (int num = 0; num < img.channels(); num++)
    {
        for (int x = 0; x < img.rows; x++)
        {
            for (int y = 0; y < img.cols; y++)
            {
                double temp = 0;
                for (int i = size / 2; i >= -size / 2; i--)
                {
                    int j = x + i - size / 2;
                    if (j < 0)
                    {
                        j = -j;
                    }
                    if (j >= img.rows)
                    {
                        j = 2 * img.rows - 1 - j;
                    }
                    temp = temp + output.ptr(j, y)[num] * weight[size / 2 - i];
                }
                output.ptr(x, y)[num] = (int)temp;
            }
        }
    }
}

int main()
{
    Mat img = imread("god.jpg");
    Mat output;
    Gaussian(img, output, 2);
    imshow("原图", img);
    imshow("滤波后",output);
    waitKey();
    return 0;
}

4.2代码：

#include<opencv2/core/core.hpp>
#include<opencv2/highgui/highgui.hpp>
#include<iostream>
using namespace cv;
using namespace std;

void MeanFilter(Mat simg, Mat oimg, int size)
{
    simg.copyTo(oimg);
    vector<vector<int> > s(simg.rows, vector<int>(simg.cols));
    for (int num = 0; num < simg.channels(); num++)
    {
        for (int i = 0; i < simg.rows; i++)
        {
            for (int j = 0; j < simg.cols; j++)
            {
                s[i][j] = 0;
            }
        }
        for (int x = 0; x < simg.rows; x++)
        {
            int sum = 0;
            for (int y = 0; y < simg.cols; y++)
            {
                sum += (int)simg.ptr(x, y)[num];
                if (x - 1 < 0)
                    s[x][y] = sum;
                else
                    s[x][y] = sum + s[x - 1][y];
            }
        }
        for (int i = 0; i < oimg.rows; i++)
        {
            for (int j = 0; j < oimg.cols; j++)
            {
                int x = i;
                int y = j;
                int x1 = i + size;
                int x2 = i - size - 1;
                int y1 = j + size;
                int y2 = j - size - 1;
                if (x1 >= simg.rows)
                    x1=x1-size;
                if (y1 >= simg.cols)
                    y1=y1-size;
                if (x2 < 0)
                    x2 = x2+size+1;
                if (y2 < 0)
                    y2 = y2+size+1;
                oimg.ptr(x, y)[num] = (s[x1][y1] + s[x2][y2] - s[x2][y1] - s[x1][y2]) / pow(2 * size + 1, 2);
            }
        }
    }
}

int main()
{
    Mat img = imread("god.jpg");
    Mat oimg;
    img.copyTo(oimg);
    MeanFilter(img, oimg, 1);
    imshow("原图", img);
    imshow("滤波后", oimg);
    waitKey();
    return 0;
}

注意：运行前要配置好opencv。