fd.get: Access Method for Objects of Class FractalDim

In fractaldim: Estimation of Fractal Dimensions

Description

For given method it returns the corresponding estimates.

Usage

fd.get(fractaldim, method)

Arguments

fractaldim	object of class FractalDim as defined in fd.estimate.
method	character string specifying the method. For 1-d estimators, possible values are “variogram”, “madogram”, “rodogram”, “variation”, “incr1”, “boxcount”, “hallwood”, “periodogram”, “wavelet”, “dctII”, and “genton”. For 2-d estimators, possible values are “transect.var”, “transect.incr1”, “isotropic”, “squareincr”, and “filter1” (see fd.get.available.methods).

Value

A FractalDim object. The original fractaldim arrays fd and scale are reduced in the last dimension into only one method, namely the given method.

See Also

fd.estimate, fd.get.available.methods

Examples

## Not run: 
library(RandomFields)
x <- seq(0, 10000)
# generate a random field
truealpha <- 1.5
rf <- GaussRF(x = x, model = "stable", grid = TRUE,
    param = c(mean=0, variance=1, nugget=0, scale=100, 
        alpha=truealpha))

#compute fractal dimension using various methods
methods <- c("madogram", "variogram", "hallwood", "boxcount",
                       "periodogram","dctII", "wavelet")
fdts <- fd.estimate (rf, methods = methods, window.size = 500, 
    step.size = 100, nlags = 10, trim = FALSE, debuglevel = 3)

# plot the variation
cols <- rainbow(length(methods))
plot(ts(fd.get (fdts, methods[1])$fd),ylim=c(min(fdts$fd), max(fdts$fd)),
             ylab="fd", col=cols[1])
for (imeth in 2:length(methods)) 
    lines(ts(fd.get (fdts, methods[imeth])$fd), col=cols[imeth])
legend('topleft', legend=methods, col=cols, lwd=1)
abline(h=2-truealpha/2)

## End(Not run)

fractaldim documentation built on Oct. 7, 2021, 9:10 a.m.