Become an expert in R — Interactive courses, Cheat Sheets, certificates and more!
Get Started for Free

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

Author(s)

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

Examples

# 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

rainfarmr

Stochastic Precipitation Downscaling with the RainFARM Method

v0.1
Apache License 2.0
Authors
Jost von Hardenberg [aut, cre, cph] (<https://orcid.org/0000-0002-5312-8070>)
Initial release

We don't support your browser anymore

Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.