0

我正在尝试在 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 过滤器对每个像素进行平均......对你们来说有什么明显的为什么会这样？

4

2 回答 2

2

您需要查看 MATLAB 的 conv2 函数。以下用于 2D 卷积的函数已从 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)
于 2013-02-05T19:10:39.080 回答
0

简单的算术错字。

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 正确执行了该函数，但是即使过滤器使用不同的算术运算符正确运行，也会引入异常（或意外结果）。

问题解决了。

于 2013-02-05T19:58:00.623 回答