Haar2MA.diag: Generate 2-D Haar MA process (diagonal/vertical/horizontal... In LS2W: Locally Stationary Two-Dimensional Wavelet Process Estimation Scheme

 Haar2MA.diag R Documentation

Generate 2-D Haar MA process (diagonal/vertical/horizontal direction).

Description

These functions generate an arbitrary number of observations from a Haar MA process of any order with a particular variance. We will focus here on `Haar2MA.diag` — the routine which generates processes having spectral structure solely in the diagonal decomposition direction.

Usage

```Haar2MA.diag(n, sd = 1, order = 5)
```

Arguments

 `n` The dimension of the realisation that you want to create. Note that `n` does NOT have to be a power of two, though it is square (nxn). `sd` The standard deviation of the innovations. `order` The order of the MA process.

Details

A two-dimensional Haar MA process is a special kind of moving-average (MA) field. A diagonal Haar MA process of order \$k\$ is a MA field of order \$2^k-1\$, the coefficients of the process being given by the filter coefficients of the two-dimensional, discrete Haar wavelet at various scales within the diagonal direction. For example: the diagonal Haar MA field of order 1 is an MA process of order 1. It is possible to define such processes for other wavelets as well.

Value

A matrix containing a realisation of the specified dimension, order and standard deviation.

Idris Eckley

References

Eckley, I.A., Nason, G.P. and Treloar, R.L. (2010) Locally stationary wavelet fields with application to the modelling and analysis of image texture. Journal of the Royal Statistical Society (Series C), 59, 595 - 616.

Eckley, I.A. and Nason, G.P. (2011). LS2W: Implementing the Locally Stationary 2D Wavelet Process Approach in R, Journal of Statistical Software, 43(3), 1-23. URL http://www.jstatsoft.org/v43/i03/.

`HaarMontage`

Examples

```#
# Generate a realisation of a diagonal component 2-D MA field
# of order 4.
#
image1 <- Haar2MA.diag(n=128, sd=3, order=4)
#
#
#
```

LS2W documentation built on Nov. 2, 2022, 1:06 a.m.