图像去噪算法报错求助：沿模糊维度扩展数组失败

错误根源 #

报错“Attempt to grow array along ambiguous dimension.”本质是数组索引越界或维度不明确：

• 当循环中的i-p或j-p小于等于0时，试图给数组A的负索引位置赋值，MATLAB无法确定如何扩展数组维度，直接抛出错误。
• 未提前初始化A导致动态扩展时索引逻辑混乱，或mean(R2)返回向量而非标量，也会间接引发维度问题。

具体修复步骤 #

1. 提前初始化输出数组
   在函数开头根据输入图像尺寸创建同大小的空数组，避免动态扩展的索引歧义：

function A = pmf(I, p)
    [rows, cols] = size(I);
    A = zeros(rows, cols); % 初始化与输入同尺寸的输出数组
    % 后续代码...
end

2. 修正循环边界，避免负索引
   循环的i和j从p+1开始，到图像对应维度的size-p结束，确保i-p和j-p始终在合法索引范围（1到对应维度大小）内：

for i = p+1 : rows - p
    for j = p+1 : cols - p
        R2 = I(i-p:i+p, j-p:j+p); % 提取p大小的邻域
        A(i,j) = mean(R2(:)); % 给当前中心像素赋值，逻辑更直观
    end
end

注：如果你的逻辑是将邻域均值赋值给邻域左上角像素，循环起始值需调整为2p+1，确保i-p不小于1。

3. 处理边缘像素
   上述循环会漏掉图像边缘p范围内的像素，可直接保留原图像边缘，或用镜像/复制边缘填充：

% 保留原图像边缘像素
A(1:p, :) = I(1:p, :);
A(rows-p+1:rows, :) = I(rows-p+1:rows, :);
A(:, 1:p) = I(:, 1:p);
A(:, cols-p+1:cols) = I(:, cols-p+1:cols);

4. 确保均值为标量
   若R2是二维邻域矩阵，默认mean(R2)会按列生成向量，无法赋值给单个元素，需改为计算整个邻域的均值：

A(i,j) = mean(R2(:)); % R2(:)将矩阵转为一维数组，mean后得到标量

验证 #

修改后无论输入图像是512x512还是其他尺寸，只要p的取值不超过图像边长的1/2，就能正常运行，不会再触发维度歧义错误。

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