今天在做多媒体作业的时候需要实现JPEG压缩，其中一步就是将量化后的块做zigzag扫描。我在网上找了很久，发现反复流传的就一份MATLAB代码，无奈这份代码跑起来结果和我预料的不一致，于是进行了修改，下面是修改后的代码。
function b=zigzag(a)
% 输入是一个8*8的矩阵，输出一个1*64的矩阵(行向量)
% 这是参照 University of California 提供的 MPEG 源代码的基础上编制的。
% Copyright (c) 1995 The Regents of the University of California.
[n,m]=size(a);
if(n~=8 && m~=8)
error(‘Input array is NOT 8-by-8’);
end
% Set up array for fast conversion from row/column coordinates to
zigzag=[1 2 9 17 10 3 4 11
18 25 33 26 19 12 5 6
13 20 27 34 41 49 42 35
28 21 14 7 8 15 22 29
36 43 50 57 58 51 44 37
30 23 16 24 31 38 45 52
59 60 53 46 39 32 40 47
54 61 62 55 48 56 63 64];
aa = reshape(a’,1,64); % 将输入块变成1×64的向量
zigzagR = reshape(zigzag’,1,64);
b = aa(zigzagR); % 对 aa 按照查表方式取元素，得到 zig-zag 扫描结果
end
验证结果如图：
然后下图是手写推理验证：
希望能对大家有所帮助
发布者：全栈程序员-站长，转载请注明出处：https://javaforall.net/213710.html原文链接：https://javaforall.net
