Description Usage Arguments Value Author(s) Examples
The input fields are transformed to Fourier space,
their Fourier spectra are merged at wavenumber kmax
after adapting the
variance of the fine field and an inverse Fourier transform is performed.
1 | mergespec(ra, r, kmax = 0)
|
ra |
matrix containing a coarse field of size |
r |
matrix containing a field at higher resolution. |
kmax |
wavenumber to use for merging (default |
The merged field, in physical space.
Jost von Hardenberg, j.vonhardenberg@isac.cnr.it
1 2 3 4 5 6 7 8 9 10 11 12 | # Make a coarse field with power-law Fourier spectrum
fa <- initmetagauss(1.7, 8)
ra <- metagauss(fa)
# Make a fine power-law Fourier spectrum
f <- initmetagauss(1.7, 32)
r <- metagauss(f)
# Merge the two fields in spectral space
rm <- mergespec(ra, r, kmax = 4)
# Check spectral slope of the resulting field
fx <- fft2d(rm)
fitslope(fx)
# 1.678554
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.