[kiko2009]

[c,cs]=wavedec2(orig,3,’haar’);
s=c(32*32+1:64*64);m=size(s,2);
[s1,ind]=sort(s);