使用blockproc实现Laplacian滤波器输出图像过黑问题求助

• 分块逻辑错误：blockproc 默认采用无重叠分块策略，你将块大小设为[3 3]，相当于每3x3个像素为一个独立块计算，相邻块的像素不会参与当前块的滤波计算，完全不符合滑动窗口滤波的要求。而nlfilter是逐像素滑动取3x3邻域计算，逻辑本身是正确的。
• 无效填充引入误差：你在自定义的laplacian_kernel中给3x3的输入块额外填充了一圈0，不仅没有补充到相邻块的真实像素，反而在块边缘引入了0值干扰，导致滤波计算结果整体偏小，这就是拉普拉斯输出图像偏暗的核心原因。
• 多余裁剪逻辑：你后续手动裁剪2列的操作是分块参数设置错误导致输出尺寸不符后的补偿操作，逻辑本身也会引入额外的尺寸误差。

blockproc要实现滑动邻域滤波，需要通过BorderSize参数给每个块配置邻域像素范围，正确的参数设置如下：

1. 块大小设为[1 1]，即每个输出像素对应一个处理块
2. 设置BorderSize = [1 1]，让每个块自动包含上下左右各1行/列的相邻像素，刚好组成3x3的滤波邻域
3. 开启TrimBorder = true，自动裁剪掉每个块的边缘填充部分，输出尺寸和原图保持一致
4. 无需手动填充、手动裁剪，内核函数直接和nlfilter的实现一致即可

修正后的代码如下：

img1 = im2double(imread('Fig0352(a)(blurry_moon).tif'));
kernel = [0 1 0; 1 -4 1; 0 1 0];

% 正确的blockproc参数配置
fun = @(block_struct) sum(block_struct.data .* kernel, 'all');
img2 = blockproc(img1, [1 1], fun, ...
    'BorderSize', [1 1], ...
    'PadPartialBlocks', true, ...
    'TrimBorder', true);

img3 = imsubtract(img1, img2);

% 显示时加[]自动拉伸动态范围，避免负值被截断导致显示偏暗
figure;
subplot(2,2,1); imshow(img1); title('原图');
subplot(2,2,2); imshow(img2, []); title('blockproc得到的拉普拉斯图像');
subplot(2,2,3); imshow(img3); title('增强后图像');

如果需要进一步提升运行速度，可直接使用imfilter实现相同功能，性能远高于blockproc实现，示例：img2 = imfilter(img1, kernel, 'replicate');，输出结果完全一致。

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