代小波harr变换整数小波变换_适合图像处理专业_matlab

%文件名：liftwavedec2.m %编 写：郭林庚 %时 间：2005/11/01 %函数功能：二代小波harr变换，整数小波变换 %输入格式举例：imgwave=liftwavedec2(image,256,3) %参数说明： % image--输入的图像矩阵，要为方阵 % m --输入的图像矩阵大小 % n --小波变换次数 %测试用例： % img=imread('lena.jpg'); % [m,m]=size(img); % imgwave=liftwavedec2(img,m,3); % imshow(imgwave); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function imgwave=liftwavedec2(image,m,n) img=image; M=m; for i=1:n imgwave1=lwavedec2(img,M); imgwave(1:M,1:M)=imgwave1; M=M/2; img=imgwave1(1:M,1:M); end % % 二代小波harr变换，整数小波变换 % function f_row=lwavedec2(image,N) f=image; T=N/2; % 子图像维数 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%% 1.列变换 % A.分裂（奇偶分开） f1=f([1:2:N-1],:); % 奇数 f2=f([2:2:N],:); % 偶数 % f1(:,T+1)=f1(:,1); % 补列 % f2(T+1,:)=f2(1,:); % 补行 % B.预测 for i_hc=1:T; high_frequency_column(i_hc,:)=f1(i_hc,:)-f2(i_hc,:); end; % high_frequency_column(T+1,:)=high_frequency_column(1,:); % 补行 % C.更新 for i_lc=1:T; low_frequency_column(i_lc,:)=f2(i_lc,:)+1/2*high_frequency_column(i_lc,:); end; % D.合并 f_column([1:1:T],:)=low_frequency_column([1:T],:); f_column([T+1:1:N],:)=high_frequency_column([1:T],:); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%% 2.行变换 % A.分裂（奇偶分开） f1=f_column(:,[1:2:N-1]); % 奇数 f2=f_column(:,[2:2:N]); % 偶数 % f2(:,T+1)=f2(:,1); % 补行 % B.预测 for i_hr=1:T; high_frequency_row(:,i_hr)=f1(:,i_hr)-f2(:,i_hr); end; % high_frequency_row(:,T+1)=high_frequency_row(:,1); % 补行 % C.更新 for i_lr=1:T; low_frequency_row(:,i_lr)=f2(:,i_lr)+1/2*high_frequency_row(:,i_lr); end; % D.合并 f_row(:,[1:1:T])=low_frequency_row(:,[1:T]); f_row(:,[T+1:1:N])=high_frequency_row(:,[1:T]); function imgwave=liftwaverec2(image,m,n) M=m/(2^n); imgwave=image; for i=n:-1:1 M=M*2; img=imgwave(1:M,1:M); imgwave1=lwaverec2(img,M); imgwave(1:M,1:M)=imgwave1; end function f_column=lwaverec2(f_row,N); T=N/2; %figure(2); %subplot(221),imshow(f_row),title('变换后图像'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%反变换%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%% 1.行变换 % A.提取（低频高频分开） f1=f_row(:,[T+1:1:N]); % 奇数 f2=f_row(:,[1:1:T]); % 偶数 % f2(:,T+1)=f2(:,1); % 补行 % B.更新 for i_lr=1:T; low_frequency_row(:,i_lr)=f2(:,i_lr)-1/2*f1(:,i_lr); end; % C.预测 for i_hr=1:T; high_frequency_row(:,i_hr)=f1(:,i_hr)+low_frequency_row(:,i_hr); end; % high_frequency_row(:,T+1)=high_frequency_row(:,1); % 补行 % D.合并(奇偶分开合并) f_row(:,[2:2:N])=low_frequency_row(:,[1:T]); f_row(:,[1:2:N-1])=high_frequency_row(:,[1:T]); %subplot(223),imshow(f_row),title('行变换图像'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%% 2.列变换 % A.提取（低频高频分开） f1=f_row([T+1:1:N],:); % 奇数 f2=f_row([1:1:T],:); % 偶数 % f1(:,T+1)=f1(:,1); % 补列 % f2(T+1,:)=f2(1,:); % 补行 % B.更新 for i_lc=1:T; low_frequency_column(i_lc,:)=f2(i_lc,:)-1/2*f1(i_lc,:); end; % C.预测 for i_hc=1:T; high_frequency_column(i_hc,:)=f1(i_hc,:)+low_frequency_column(i_hc,:); end; % high_frequency_column(T+1,:)=high_frequency_column(1,:); % 补行 % D.合并(奇偶分开合并） f_column([2:2:N],:)=low_frequency_column([1:T],:); f_column([1:2:N-1],:)=high_frequency_column([1:T],:); %subplot(224),imshow(f_column),title('行变换图像');