Multidimensional signal processing

In signal processing, multidimensional signal processing covers all signal processing done using multidimensional signals and systems. While multidimensional signal processing is a subset of signal processing, it is unique in the sense that it deals specifically with data that can only be adequately detailed using more than one dimension. In m-D digital signal processing, useful data is sampled in more than one dimension. Examples of this are image processing and multi-sensor radar detection. Both of these examples use multiple sensors to sample signals and form images based on the manipulation of these multiple signals. Processing in multi-dimension (m-D) requires more complex algorithms, compared to the 1-D case, to handle calculations such as the fast Fourier transform due to more degrees of freedom. In some cases, m-D signals and systems can be simplified into single dimension signal processing methods, if the considered systems are separable.

Typically, multidimensional signal processing is directly associated with digital signal processing because its complexity warrants the use of computer modelling and computation. A multidimensional signal is similar to a single dimensional signal as far as manipulations that can be performed, such as sampling, Fourier analysis, and filtering. The actual computations of these manipulations grow with the number of dimensions.

Sampling
Multidimensional sampling requires different analysis than typical 1-D sampling. Single dimension sampling is executed by selecting points along a continuous line and storing the values of this data stream. In the case of multidimensional sampling, the data is selected utilizing a lattice, which is a "pattern" based on the sampling vectors of the m-D data set. These vectors can be single dimensional or multidimensional depending on the data and the application.

Multidimensional sampling is similar to classical sampling as it must adhere to the Nyquist–Shannon sampling theorem. It is affected by aliasing and considerations must be made for eventual Multidimensional Signal Reconstruction.

Fourier Analysis
A multidimensional signal can be represented in terms of sinusoidal components. This is typically done with a type of Fourier transform. The m-D Fourier transform transforms a signal from a signal domain representation to a frequency domain representation of the signal. In the case of digital processing, a discrete Fourier Transform (DFT) is utilized to transform a sampled signal domain representation into a frequency domain representation:
 * $$ X(k_1,k_2,\dots,k_m) = \sum_{n_1=-\infty}^\infty \sum_{n_2=-\infty}^\infty \cdots \sum_{n_m=-\infty}^\infty x(n_1,n_2,\dots,n_m) e^{-j 2 \pi k_1 n_1} e^{-j 2 \pi k_2 n_2} \cdots e^{-j 2 \pi k_m n_m}$$

where X stands for the multidimensional discrete Fourier transform, x stands for the sampled time/space domain signal, m stands for the number of dimensions in the system, n are sample indices and k are frequency samples. Computational complexity is usually the main concern when implementing any Fourier transform. For multidimensional signals, the complexity can be reduced by a number of different methods. The computation may be simplified if there is independence between variables of the multidimensional signal. In general, fast Fourier transforms (FFTs), reduce the number of computations by a substantial factor. While there are a number of different implementations of this algorithm for m-D signals, two often used variations are the vector-radix FFT and the row-column FFT.

Filtering
Filtering is an important part of any signal processing application. Similar to typical single dimension signal processing applications, there are varying degrees of complexity within filter design for a given system. M-D systems utilize digital filters in many different applications. The actual implementation of these m-D filters can pose a design problem depending on whether the multidimensional polynomial is factorable. Typically, a prototype filter is designed in a single dimension and that filter is extrapolated to m-D using a mapping function. One of the original mapping functions from 1-D to 2-D was the McClellan Transform. Both FIR and IIR filters can be transformed to m-D, depending on the application and the mapping function.

Applicable Fields

 * Audio signal processing
 * Image processing
 * Towed array sonar
 * X-ray computed tomography