mapping.m

采用小波零树的EZW压缩编码算法

function A = mapping(n)

% To create a Morton Scan order matrix

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