Compute spatial Fourier power spectrum
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 |
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
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.