MATLAB自定义图像模糊函数中累加变量异常归零问题排查

你的问题根源在于uint8类型的数值溢出，并非s被意外重置为0，让我帮你拆解清楚：

为什么会出现这个现象？ #

你的输入图像img是uint8类型，数值范围被限制在0-255之间。当执行s=s+img(r,c);时，虽然s初始是double类型的0，但MATLAB在混合类型运算时，会自动将结果转换为范围更窄的类型（也就是uint8）。一旦累加的和超过255，uint8会自动对数值取模256（简单说就是截断溢出）——比如256会变成0，257变成1，这就会让你误以为s的前值被重置为0了。

如何修复？ #

解决方法非常直接：在累加前把img(r,c)转换为更高精度的数值类型（比如double），这样整个累加过程都会保持高精度，不会触发溢出。只需要修改你标记的问题行：

s = s + double(img(r,c)); % 将uint8转成double后再累加

优化建议（可选） #

如果想要提升代码效率，可以在循环开始前就把整个图像转换成double类型，避免每次循环里重复执行类型转换：

function output=bl(img,w)
row=size(img,1);
col=size(img,2);
g_size=(2*w+1)^2;
output=uint8(zeros(row,col));
img_double = double(img); % 提前转换整个图像为高精度类型
for ii=1:row
    for jj=1:col
        s=0;
        corner=[ii jj]-w;
        for r=corner(1):(corner(1)+(2*w))
            for c=corner(2):(corner(2)+(2*w))
                if((r>0 && r<=row) && (c>0 && c<=col))
                    s=s+img_double(r,c); % 直接使用转换后的矩阵
                    fprintf('The value of s is %d\n',s);
                end
            end
        end
        m=s/g_size;
        output(ii,jj)=fix(m);
        s=0;
    end
end
end

另外，你的边界处理通过判断坐标是否在图像范围内实现，逻辑是通顺的；如果想让代码更简洁，也可以考虑用padarray先给图像填充边界（比如填充0或者复制边缘），这样就不用在循环里做条件判断了，但这属于进阶优化，核心问题还是类型溢出的问题。

内容的提问来源于stack exchange，提问作者Abdullah Shahriar