5/3提升小波变换源代码

%-----输入图像-----% clc clear close all fid=fopen('airport_5m_25m.dat','r'); image_i=fread(fid,510*1200*2,'uchar'); image_i=reshape(image_i,510,1200*2); fclose(fid); figure; imagesc(abs(image_i));%display image object colormap(gray); title('原始图像'); disp('压缩前图像的大小：'); whos('image_i'); % % %-----图像扩展-----% % %向四周扩展的尺寸为20个点 % expand_value=10; % image_o=image_expand(image_i,expand_value); % % %-----图像预处理-----% % %无极性处理是对图像数据进行无极性变换和处理，针对无极性分量采 % % 样，去掉信号ＤＣ偏移值。ＳＡＲ图像数据是８位灰度图像，各 % % 像素点采样值I(x,y)（０≤I(x,y)<2B）是Ｂbit正值信 % % 号。经处理后的分量采样变换为以０值为中心的-2B-1≤I'(x,y) % % ≤2B-1范围内的带极性采样值。 image_i=image_i-128; %-----小波变换------% [m,n]=size(image_i); m1=fix(m/2);%截尾取整 n1=fix(n/2);%截尾取整 %------列变换-------% %A.截取 image_Y_O1=image_i((1:2:m-1),:); %奇数矩阵 image_Y_E1=image_i((2:2:m),:); %偶数矩阵 %B.预测 image_Y_H1(1,:)=image_Y_O1(1,:)-image_Y_E1(1,:);%注意，第一行和最后一行的预测公式和其他所有行的不一样 for i=2:m1-1 image_Y_H1(i,:)=image_Y_O1(i,:)-1/2*(image_Y_E1(i-1,:)+image_Y_E1(i,:)); end image_Y_H1(end+1,:)=image_Y_O1(m1,:)-image_Y_E1(m1,:);%注意，第一行和最后一行的预测公式和其他所有行的不一样 %C.更新 image_Y_L1(1,:)=image_Y_E1(1,:)+1/2*image_Y_H1(1,:)+1/2;%注意，第一行和最后一行的预测公式和其他所有行的不一样 for i=2:m1-1 image_Y_L1(i,:)=image_Y_E1(i,:)+1/4*((image_Y_H1(i-1,:)+image_Y_H1(i,:)))+1/2; end image_Y_L1(end+1,:)=image_Y_E1(m1,:)+1/2*image_Y_H1(m1,:)+1/2;%注意，第一行和最后一行的预测公式和其他所有行的不一样 %D.合并 %低频信号在矩阵第1至第m1列，高频矩阵在矩阵第m1+1至m列 image_Y_LH(1:1:m1,:)=image_Y_L1; image_Y_LH(m1+1:1:m,:)=image_Y_H1; %------行变换-------% %A.截取 image_Y_O2=image_Y_LH(:,(1:2:n-1)); %奇数矩阵 image_Y_E2=image_Y_LH(:,(2:2:n)); %偶数矩阵 %B.预测 image_Y_H2(:,1)=image_Y_O2(:,1)-1/2*image_Y_E2(:,1);%注意，第一列和最后一列的预测公式和其他所有列的不一样 for i=2:n1-1 image_Y_H2(:,i)=image_Y_O2(:,i)-1/2*(image_Y_E2(:,i-1)+image_Y_E2(:,i)); end image_Y_H2(:,end+1)=image_Y_O2(:,n1)-1/2*image_Y_E2(:,n1);%注意，第一列和最后一列的预测公式和其他所有行的不一样 %C.更新 image_Y_L2(:,1)=image_Y_E2(:,1)+1/2*image_Y_H2(:,1)+1/2;%注意，第一列和最后一列的预测公式和其他所有行的不一样 for i=2:n1-1 image_Y_L2(:,i)=image_Y_E2(:,i)+1/4*(image_Y_H2(:,i-1)+image_Y_H2(:,i))+1/2; end image_Y_L2(:,end+1)=image_Y_E2(:,n1)+1/2*image_Y_H2(:,n1)+1/2;%注意，第一列和最后一列的预测公式和其他所有行的不一样 %D.合并 %低频信号在矩阵第1至第m1列，高频矩阵在矩阵第m1+1至m列 image_Y_LH(:,1:1:n1)=image_Y_L2; image_Y_LH(:,n1+1:1:n)=image_Y_H2; % image_53_tsxb=image_Y_LH; image_53_tsxb=image_Y_LH(1:m1,1:n1);%压缩后的图像 figure; % imagesc(abs(image_97_tsxb));%display image object imagesc(image_53_tsxb);%display image object colormap(gray); title('一次压缩图像'); disp('一次压缩图像的大小：'); whos('image_97_tsxb');