分形图像压缩代码.docVIP

分形图像压缩（PIFS法）代码 原代码下载自CSDN网站 clear; outfp = fopen(zijifenxingyasuo.txt,w); [DATA,map]=imread(lena.bmp,bmp); DATA=double(DATA); size1=size(DATA); tu=DATA(1:size1(1),1:size1(2),1); xx=size1(2); yy=size1(1); nrx=xx/8; nry=yy/8; ndx=xx/16; ndy=yy/16; DD=zeros(ndy,ndx,8,64); cund1=zeros(16,16); cund2=zeros(8,8); cund3=zeros(1,64); cund33=zeros(1,64); cund4=zeros(8,8); RRR=zeros(nry,nrx,6); % s,o,D(y),D(x),U cunr=zeros(1,64); for i=1:ndy for j=1:ndx % for k=1:16 cund1=tu(1+16*(i-1):16+16*(i-1),1+16*(j-1):16+16*(j-1)); for l=1:8 for m=1:8 cund2(l,m)=(cund1(1+2*(l-1),1+2*(m-1))+cund1(2+2*(l-1),2+2*(m-1))+cund1(2+2*(l-1),1+2*(m-1))+cund1(1+2*(l-1),2+2*(m-1)))/4; end; end; % end; DD(i,j,1,1:64)=reshape(cund2,[1,64]); cund4=fliplr(cund2); DD(i,j,2,1:64)=reshape(cund4,[1,64]); cund4=flipud(cund2); DD(i,j,3,1:64)=reshape(cund4,[1,64]); cund4=flipud(fliplr(cund2)); DD(i,j,4,1:64)=reshape(cund4,[1,64]); cund4=rot90(flipud(cund2)); DD(i,j,5,1:64)=reshape(cund4,[1,64]); cund4=rot90(cund2); DD(i,j,6,1:64)=reshape(cund4,[1,64]); cund4=rot90(rot90(rot90(cund2))); DD(i,j,7,1:64)=reshape(cund4,[1,64]); cund4=cund2; DD(i,j,8,1:64)=reshape(cund4,[1,64]); % DD(i,j,1,1:8,1:8)=cund2; % DD(i,j,2,1:8,1:8)=fliplr(cund2); % DD(i,j,3,1:8,1:8)=flipud(cund2); % DD(i,j,4,1:8,1:8)=flipud(fliplr(cund2)); % DD(i,j,5,1:8,1:8)=Rot90(flipud(cund2)); % DD(i,j,6,1:8,1:8)=Rot90(cund2); % DD(i,j,7,1:8,1:8)=Rot90(Rot90(Rot90(cund2))); % DD(i,j,8,1:8,1:8)=cund2; end; end; for i=1:nry for j=1:nrx cunr=reshape(tu(1+8*(i-1):8+8*(i-1),1+8*(j-1):8+8*(j-1)),[1,64]); sumalpha=sum(cunr); %cunr is ri sumalpha2=norm(cunr)^2; dx=1; dy=1; ut=1;