mapping.m

ezw编码

%生成一个与图像同样大小、存储着扫描序号的矩阵
function A = mapping(n)
if n == 2
    A = [1 2; 3 4];
else
    B = mapping(n/2);
    A = [B B+(n/2)^2; B+(n/2)^2*2 B+(n/2)^2*3];
end