dwt_2d.m

来自「小波压缩编码平台软件包」· M 代码 · 共 36 行

function [A,H,V,D] = dwt_2D(X,wavelet)
%[A,H,V,D] = dwt_2D(X,wavelet)
%Version: 1.00, Date: 2006/01/31, author: Nikola Sprljan
%Two-dimensional separable DWT
%
%Input: 
% X - matrix to be transformed containing the input signal
% wavelet - wavelet identification string, or wavelet data structure
%     
%Output: 
% A,H,V,D - approximation signal, horizontal, vertical and diagonal details
%           signal
%
%Uses:
% load_wavelet.m (Wavelet Toolbox)
% dwt_dim.m (Wavelet Toolbox)
% subband.m (Wavelet Toolbox)
%
%Example:
% [A,H,V,D] = dwt_2D(X,'CDF_9x7');

%load the wavelet here
if (isstr(wavelet))
    wvf = load_wavelet(wavelet);
else
    wvf = wavelet;
end;

X = dwt_dim(X,1,wvf); %columns
X = dwt_dim(X,2,wvf); %rows

A = subband(X,1,'ll');
H = subband(X,1,'hl');
V = subband(X,1,'lh');
D = subband(X,1,'hh');