waveup3.m

来自「matlab实现s小波变换代码,m文件。输入图像矩阵」· M 代码 · 共 43 行

function wvlet=waveup3(x)
%function [ll,hl,lh,hh]=waveup3(x)
x=double(x);
wvlet=x;
[hp,lp]=size(x);
hc=hp/2;
lc=lp/2;
%奇偶列重排
for n=1:lc
	j=n*2-1;
	wvlet(:,n)=x(:,j);
	wvlet(:,lc+n)=x(:,j+1);
end
%行变换开始
for n=1:lc	
	wvlet(:,lc+n)=wvlet(:,lc+n)-wvlet(:,n);
end
for n=1:lc
    k=wvlet(:,lc+n)/2-mod(wvlet(:,lc+n)/2,1);
	wvlet(:,n)=wvlet(:,n)+k;
end
%行变换结束
x=wvlet;
%奇偶行重排
for n=1:hc
	j=n*2-1;
	wvlet(n,:)=x(j,:);
	wvlet(hc+n,:)=x(j+1,:);
end
%列重排开始	
for n=1:hc	
	wvlet(hc+n,:)=wvlet(hc+n,:)-wvlet(n,:);
end
for n=1:hc
	k=wvlet(hc+n,:)/2-mod(wvlet(hc+n,:)/2,1);
	wvlet(n,:)=wvlet(n,:)+k;
end
%列重排结束
%wvlet=uint8(wvlet);
ll=wvlet(1:hc,1:lc);
hl=wvlet(1:hc,lc+1:lp);
lh=wvlet(hc+1:hp,1:lc);
hh=wvlet(hc+1:hp,lc+1:lp);