Question

我正在尝试在Matlab中创建一个模糊我图像的函数。我正在使用Matlabs演示图像peppers.png。

这是我的功能：

    function g = myfilter(f, h)

    f = double(f); %convert to double
    g = zeros(size(f)); %new array (size of f)
    a = (size(h, 1) - 1) / 2; %padding on edges

    for row = (a + 1) : (size(f,1) - a)
        for col = (a + 1) : (size(f,2) - a)

            gxy = 0; %running sum

            for m = -a:a
                for n = -a:a

                    gxy = gxy + f(row - m, col - n) + h(m + a+1, n + a+1);
                end
            end

            g(row, col) = gxy;
        end
    end

    g = uint8(g); %convert back to int

以下是我的命令：

    >> img = imread('peppers.png');
    >> imshow(img)
    >> imgGray = rgb2gray(img);
    >> imshow(imgGray)
    >> 
    >> filt1 = (1/9)*ones(3)

       filt1 =

        0.1111    0.1111    0.1111
        0.1111    0.1111    0.1111
        0.1111    0.1111    0.1111

   >> test = myfilter(imgGray, filt1);
   >> imshow(test)

成功将彩色图像转换为灰色并应用滤镜。

不幸的是，滤镜只会创建一个几乎完整的白色图像（太亮）...我根本看不出原因......应该使用3x3滤镜对每个像素进行平均... 对你们来说，为什么会发生这种情况有什么明显的？

Answer 1

您需要查看MATLAB的conv2函数。 2D convolution的以下函数已经从conv2中提取出来，对你给定的滤波器非常有用。

function c = myfilter(a, b)
 [ma, na] = size(a);
 [mb, nb] = size(b);
c = zeros( ma+mb-1, na+nb-1 );
    for i = 1:mb
        for j = 1:nb
            r1 = i;
            r2 = r1 + ma - 1;
            c1 = j;
            c2 = c1 + na - 1;
            c(r1:r2,c1:c2) = c(r1:r2,c1:c2) + b(i,j) * a;
        end
    end
c = uint8(c)

Answer 2

简单的算术错字。

gxy = gxy + f(row-m, col-n) + h(m+a+1, n+a+1);

应该是：gxy = gxy + f(row-m, col-n) * h(m+a+1, n+a+1);。

它工作正常，现在会创建一个模糊的图像。

它们不是乘以f和h，而是在上面的代码中求和，它们不符合2D卷积定义的空间域图像滤镜。 Matlab正在正确执行该函数，但是引入了异常（或意外结果），即使过滤器使用不同的算术运算符正确运行。

问题解决了。

使用函数的Matlab模糊滤镜

2 个答案: