rainfarmr (version 0.1)

fft2d: Compute spatial Fourier power spectrum

Description

The Fourier power spectrum of the input field is computed and averaged over shells (bins with width 1) of the modulus of the wavenumber.

Usage

fft2d(z)

Arguments

z

matrix or array of input data with spatial dimensions c(N, N). If z is an array the power spectra are averaged over the third dimension. The first two dimensions must be the same (the input fields must be square).

Value

Spectral power (average square of absolute value of spectral amplitudes) for wavenumbers k=1:(N/2).

Examples

# NOT RUN {
# Make a synthetic rainfall field with prescribed logarithmic spectral slope
f = initmetagauss(1.7, 64)
r = metagauss(f)
# Check slope of the resulting field
fx <- fft2d(r)
fitslope(fx)
# 1.640373
# }