fft2d: Compute spatial Fourier power spectrum
In rainfarmr: Stochastic Precipitation Downscaling with the RainFARM Method

Description Usage Arguments Value Author(s) Examples

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

fft2d(z)

`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).

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

Jost von Hardenberg, j.vonhardenberg@isac.cnr.it

# 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