www.pudn.com > SPIHT_bandelet.rar > btc_image_encode.asv

function out_put=btc_image_encode(in_put,block_size) 
X=in_put; 
K=block_size; 
X1=double(X); 
y1=size(X); 
n=y1(1); 
m=y1(2); 
k=1;l=1; 
 
for i=1:K:n 
    for j=1:K:m 
        tmp([1:K],[1:K])=X1([i:i+(K-1)],[j:j+(K-1)]); 
        mn=mean(mean(tmp));%mn为块内所有像素的均值 
        tmp1([i:i+(K-1)],[j:j+(K-1)])=tmp>mn;%位元图tmp1 
        Lsmat=(tmp=mn); 
        Lsmn=sum(sum(Lsmat));%相当于算法中所说的w 
        Mrmn=sum(sum(Mrmat));%相当于算法中所说的q 
        Mu(k)=sum(sum(Lsmat.*tmp))/(Lsmn+.5);k=k+1; 
        Mi(l)=sum(sum(Mrmat.*tmp))/Mrmn;l=l+1; 
    end 
end 
out_put=tmp1;