Description Usage Arguments Details Value See Also Examples
Compute the r-spectrum of a matrix
1 | rspectrum(mat)
|
mat |
A matrix with logical or numeric values |
This functions returns a data.frame with NA
s in the rspec
column if the input matrix has zero variance. Note that if the matrix
is not square, then only the largest square matrix fitting in the upper
right corner is used.
A data.frame with two columns: dist
, the wave number and
rspec
, the normalized value of the r-spectrum
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 | # Spectrum of white noise
rmat <- matrix(runif(100*100) > .5, ncol = 100)
spec <- rspectrum(rmat)
plot(spec, type = "l")
# Add some spatial correlation and compare the two spectra
rmat.cor <- rmat
for (i in seq(1, nrow(rmat)-1)) {
for (j in seq(1, nrow(rmat)-1)) {
rmat.cor[i,j] <- mean(rmat[(i-1):(i+1), (j-1):(j+1)])
}
}
spec.cor <- rspectrum(rmat.cor)
plot(spec.cor, type = "n")
lines(spec, col = "black")
lines(spec.cor, col = "blue")
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.