手动实现Matlab图像LoG卷积（弃用conv2）：代码报错求助

我来帮你排查代码里的问题，你的手动卷积实现有几个关键错误，我会一步步帮你修正：

错误点分析 #

• 索引范围逻辑错误：你在提取图像子块时，行和列的索引完全写反了！5×5滤波器需要取以(r,c)为中心的区域，正确的行范围是r-2 : r+2，列范围是c-2 : c+2，而不是代码里的(r-2):(c-2)和(r+2):(c+2)，这会导致子块维度不匹配，直接触发报错。
• 灰度转换公式错误：RGB转灰度时，你误将蓝色通道（第三个分量）写成了红色通道（第一个分量），标准公式是G(x,y) = 0.3*R + 0.59*G + 0.11*B，你的代码里第三个项用了Img(x,y,1)，应该改成Img(x,y,3)。
• 原图像被覆盖：直接在原始图像矩阵I上修改像素值，后续循环会用到已修改的结果，导致卷积计算完全错误。需要创建一个新的输出矩阵来存储结果。
• 图像显示异常：LoG卷积结果会包含正负值，直接转uint8会截断有效信息，需要先将结果归一化到0-255范围再显示。

修正后的完整代码 #

clc 
clear; 
Img = imread('11.jpg'); % 读取图像

% 转换为灰度图（修正公式错误）
G = zeros(size(Img,1), size(Img,2));
for x = 1:size(Img,1)
    for y = 1:size(Img,2)
        G(x,y) = Img(x,y,1)*0.3 + Img(x,y,2)*0.59 + Img(x,y,3)*0.11;
    end
end

Mask = [0,0,-1,0,0;0,-1,-2,-1,0;-1,-2,16,-2,-1;0,-1,-2,-1,0;0,0,-1,0,0]; % LoG滤波器模板
I = double(G);
Mask = double(Mask); 

% 创建输出矩阵，避免覆盖原图像
output = zeros(size(I));

% 手动实现卷积（修正索引错误）
for r = 3: size(I,1)-2
    for c = 3: size(I,2)-2
        % 提取以(r,c)为中心的5×5子块
        sub_block = I(r-2:r+2, c-2:c+2);
        % 计算模板与子块的点积和
        output(r,c) = sum(sum(Mask .* sub_block));
    end
end

% 归一化结果到0-255范围，避免显示异常
output = mat2gray(output);
output = uint8(output * 255);

% 显示图像
imshow(Img);title('Original Image'); 
figure, imshow(G);title('Gray Image'); 
figure, imshow(output);title('Processed Image');

修正说明 #

1. 灰度转换：修正了蓝色通道的索引，确保灰度图的计算符合标准公式。
2. 输出矩阵：新增output矩阵存储卷积结果，彻底避免修改原图像I导致的连锁计算错误。
3. 索引修正：将子块提取改为r-2:r+2, c-2:c+2，确保每次都能取到以当前像素为中心的5×5区域。
4. 结果归一化：使用mat2gray将卷积结果归一化到0-1范围，再缩放为0-255的uint8格式，保证图像显示的细节完整。

内容的提问来源于stack exchange，提问作者محمد راضي عزيز البديري