【C#】图像滤镜(一):图像平滑

一、高斯平滑

1. 高斯核计算

    private double[,] GaussianKernel(int N)
        {
            int i, j;
            double sigma = 1;
            double[,] gussian = new double[N, N];
            double sum = 0.0;
            for (i = 0; i < N; i++)
            {
                for (j = 0; j < N; j++)
                {
                    gussian[i, j] = Math.Exp(-((i - N / 2) * (i - N / 2) + (j - N / 2) * (j - N / 2)) / (2.0 * sigma * sigma));
                    sum += gussian[i, j];
                }
            }
            for (i = 0; i < N; i++)
            {
                for (j = 0; j < N; j++)
                {
                    gussian[i, j] = gussian[i, j] / sum;

                }
            }
            return gussian;
        }

2. 图像扩展

    private Image<Bgr, byte> ExtendImg(Image<Bgr, byte> img, int size)
        {
            int w = img.Width;
            int h = img.Height;
            int nw = w + size * 2;
            int nh = h + size * 2;

            Image<Bgr, byte> nimg = new Image<Bgr, byte>(nw, nh);

            for (int i = 0; i < size; i++)
                for (int j = size; j < nh - size; j++)
                    nimg[j, i] = img[j - size, i];
            for (int i = nw - size; i < nw; i++)
                for (int j = 0; j < h; j++)
                    nimg[j, i] = img[j, i - size * 2];
            for (int i = size; i < nw - size; i++)
                for (int j = 0; j < size; j++)
                    nimg[j, i] = img[j, i - size];
            for (int i = size; i < nw - size; i++)
                for (int j = nh - size; j < nh; j++)
                    nimg[j, i] = img[j - size * 2, i - size];

            for (int i = 0; i < size; i++)
                for (int j = 0; j < size; j++)
                    nimg[j, i] = img[j, i];
            for (int i = nw - size; i < nw; i++)
                for (int j = 0; j < size; j++)
                    nimg[j, i] = img[j, i - size * 2];
            for (int i = nw - size; i < nw; i++)
                for (int j = nh - size; j < nh; j++)
                    nimg[j, i] = img[j - size * 2, i - size * 2];
            for (int i = 0; i < size; i++)
                for (int j = nh - size; j < nh; j++)
                    nimg[j, i] = img[j - size * 2, i];
            for (int i = 0; i < w; i++)
                for (int j = 0; j < h; j++)
                    nimg[j + size, i + size] = img[j, i];
            return nimg;
        }

3. 卷积计算

 private Image<Bgr, byte> Convolution(Image<Bgr, byte> pic, double[,] k)
        {
            int d = k.GetLength(0) / 2;
            Image<Bgr, byte> temp = ExtendImg(pic, d);
            Image<Bgr, byte> output = new Image<Bgr, byte>(pic.Size);
            double b = 0, g = 0, r = 0;
            for (int i = d; i < temp.Width - d; i++)
            {
                for (int j = d; j < temp.Height - d; j++)
                {
                    b = 0; g = 0; r = 0;
                    for (int m = -d; m < d; m++)
                    {
                        for (int n = -d; n < d; n++)
                        {
                            b += temp[j + m, i + n].Blue * k[m + d, n + d];
                            g += temp[j + m, i + n].Green * k[m + d, n + d];
                            r += temp[j + m, i + n].Red * k[m + d, n + d];
                        }
                    }
                    output[j - d, i - d] = new Bgr(b, g, r);
                }
            }
            return output;
        }

3. 结果

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

勇敢牛牛_

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值