WMDCT2 - 2D Discrete windowed MDCT transform

Usage

c=wmdct2(f,g,M);
c=wmdct2(f,g1,g2,[M1,M2]);
c=wmdct2(f,g1,g2,[M1,M2],[L1,L2]);
[c,L]=wmdct2(f,g1,g2,[M1,M2],[L1,L2]);

Input parameters

`f`	Input data, matrix.
`g, g1, g2`
	Window functions.
`M, M1, M2`
	Number of bands.
`L1, L2`	Length of transform to do.

Output parameters

`c`	array of coefficients.
`Ls`	Original size of input matrix.

Description

wmdct2(f,g,M) calculates a two dimensional Modified Discrete Cosine transform of the input signal f using the window g and parameter M along each dimension.

For each dimension, the length of the transform will be the smallest possible that is larger than the length of the signal along that dimension. f will be appropriately zero-extended.

All windows must be whole-point even.

wmdct2(f,g,M,L) computes a 2D windowed MDCT as above, but does a transform of length L along each dimension. f will be cut or zero-extended to length L before the transform is done.

[c,Ls]=wmdct(f,g,M) or [c,Ls]=wmdct(f,g,M,L) additionally return the length of the input signal f. This is handy for reconstruction.

c=wmdct2(f,g1,g2,M) makes it possible to use different windows along the two dimensions.

The parameters L, M and Ls can also be vectors of length 2. In this case the first element will be used for the first dimension and the second element will be used for the second dimension.

The output c has 4 or 5 dimensions. The dimensions index the following properties:

Number of translation along 1st dimension of input.

Number of channel along 1st dimension of input

Number of translation along 2nd dimension of input.

Number of channel along 2nd dimension of input

Plane number, corresponds to 3rd dimension of input.