试图实现imfilter，但图像变红

Question

所以我想在matlab中实现imfilter，而不使用函数本身。这是我到现在为止所做的：

   % Menampilkan gambar sebelum di smoothing
    gambar = imread('bird.jpg');
    ukuran = size(gambar);
    figure
    imshow(gambar)

    % Menampilkan gambar setelah di smoothing
    gambardouble = im2double(gambar);
    filter = ones(3) * 1/9;
    gambarfinal = zeros(ukuran);

for i = 2:ukuran(1)-1
    for j = 2:ukuran(2)-1
        temp= gambardouble(i-1:i+1,j-1:j+1) .* filter;
        gambarfinal(i,j) = sum(temp(:));
    end
end

disp(gambarfinal);
figure
imshow(gambarfinal);

现在，正常的图像很好。但输出变成这样的红色：

任何想法出了什么问题？

Answer 1

gambar似乎是RGB图像，大小为NxMx3。在gambarfinal = zeros(ukuran);中，您生成一个大小为NxMx3的零矩阵。但是，您只将滤镜应用于红色波段，而将其他两个波段设为零=黑色，从而产生红色图像。

要将滤镜应用于其他乐队，您可以立即对所有3个波段应用滤波器：

...

% generate 3-d filter
filter = ones(3,3,3) * 1/9;

...

% extract a 3x3x3 cube instad of a 3x3 patch
temp= gambardouble(i-1:i+1,j-1:j+1,:) .* filter;

% sum first along 1. dimension, then along the 2.dimension. This gives a 1x1x3 column.
gambarfinal(i,j,:) = sum(sum(temp,1),2); 

如果你对它感觉更舒服，你可以在其他两个乐队上产生一个循环 - 但这会慢得多。

修改

冒号运算符:返回该维度中矩阵的所有元素。例如。对于二维矩阵A，A(:,2)将返回A的第二列。同样，对于3维矩阵A(2,5,:)将返回沿3. Dimension（向量）的所有元素。另请注意，A(2,5)与A(2,5,1)隐含相同。

在您的情况下，gambardouble(i-1:i+1,j-1:j+1)将返回3维矩阵的1.层（红色）的3x3窗口。通过将其更改为gambardouble(i-1:i+1,j-1:j+1,:)，您将希望上面的示例在其他图层中获得相同的像素。因此，在您的情况下，temp将是一个3x3x3多维数据集。

sum(temp,1)沿着1.维度对立方体求和，返回1x3x3矩阵。然后通过使用sum(... , 2)沿着2.维度获取总和，将其转换为1x1x3向量，其中每个元素属于一个波段。

gambarfinal(i,j,:)是位置(i,j)处的RGB值，我们将其替换为新值。

有关详情，请参阅colon documentation