Matlab：imread读取PGM图调用blockproc报错“Dot indexing...”求助

核心问题 #

imread读取JPG图片时，返回的是包含coef_arrays、quant_tables等字段的结构体；但读取PGM图片时，返回的只是普通的数值矩阵，根本没有这些属性。你代码里用stego.coef_arrays、stego.quant_tables这种点索引操作，对数值矩阵完全不支持，这就是报错的直接原因。

解决步骤 #

PGM是无压缩灰度格式，本身没有JPEG专属的量化表、DCT系数数组，得手动完成这些流程：

1. 把图像转成适合处理的数值格式
2. 手动定义量化表（可复用JPEG标准灰度量化表，或自定义）
3. 分块执行DCT变换、量化操作

修改后的代码示例：

% 读取PGM图像，得到数值矩阵
stego = imread('1.pgm');
% 转换为double类型，方便后续变换计算
stego_double = im2double(stego);
quality_factor = 75;

% 定义JPEG标准灰度量化表，根据质量因子调整
% 若MATLAB版本支持，可直接用内置函数生成
quant_table = jpeg_read_quant_table(quality_factor, 'grayscale');
% 若没有内置函数，手动定义标准表后缩放
% quant_table = [16 11 10 16 24 40 51 61; 12 12 14 19 26 58 60 55; ...
%                14 13 16 24 40 57 69 56; 14 17 22 29 51 87 80 62; ...
%                18 22 37 56 68 109 103 77; 24 35 55 64 81 104 113 92; ...
%                49 64 78 87 103 121 120 101; 72 92 95 98 112 100 103 99];
% scale = quality_factor < 50 ? 5000/quality_factor : 200 - 2*quality_factor;
% quant_table = round((quant_table * scale + 50)/100);
% quant_table(quant_table < 1) = 1;
% quant_table(quant_table > 255) = 255;

% 定义块处理函数：DCT变换→量化→逆DCT变换
fun = @(x) idct2(dct2(x.data) .* quant_table);
% 分8x8块处理图像
I_spatial = blockproc(stego_double, [8 8], fun);
% 转换回uint8格式（按需选择）
I_spatial = im2uint8(I_spatial);

补充说明 #

如果你的需求是对PGM做类JPEG的量化压缩，必须手动完成DCT、量化、逆变换的完整流程——毕竟PGM本身不携带任何JPEG特有的元数据。

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