compress.m

用matlab写的图像预测编码

%用 J(k,l)=I(k,l)-(I(k,l-1)/2+I(k-1,l)/4+I(k-1,l-1)/8+I(k-1,l+1)/8)进行预测编码
I = imread('F:\标准图片\mri.pgm');
I2=I;
I=double(I);
fid = fopen('mydata.dat','w');
[m n] = size(I);
J=ones(m,n);
J(1:m,1)=I(1:m,1);
J(1,1:n)=I(1,1:n);
J(1:m,n)=J(1:m,n);
J(m,1:n)=J(m,1:n);
for k=2:m-1
    for l=2:n-1
        J(k,l)=I(k,l)-(I(k,l-1)/2+I(k-1,l)/4+I(k-1,l-1)/8+I(k-1,l+1)/8);
    end
end
J=round(J);
cout=fwrite(fid,J,'int8');
cont=fwrite(fid,J,'int8');
cc=fclose(fid);