%預(yù)測編碼的解壓縮過程
fid=fopen('mydata.dat','r');
I1=fread(fid,inf,'int8');
tt=1;
for l=1:n
    for k=1:m
        I(k,l) = I1(tt);
        tt=tt+1;
    end
end
I=double(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)=I(1:m,n);
J(m,1:n)=I(m,1:n);
for k=2:m-1
    for l=2:n-1
    J(k,l)=I(k,l)+(J(k,l-1)/2+J(k-1,l)/4+J(k-1,l-1)/8+J(k-1,l+1)/8);
    end
end
cc=fclose(fid);
J=uint8(J);
subplot(1,2,1),imshow(I2);
subplot(1,2,2),imshow(J);