View source: R/estimateAnisotropy.R
rotateAnisotropicData | R Documentation |
This function applies an isotropic transformation of
the coordinates specified in object
.
rotateAnisotropicData(object,anisPar)
object |
(i) An Intamap type object (see |
anisPar |
An array containing the anisotropy parameters (anisotropy ratio and axes orientation)
(see |
This function performs a rotation and rescaling of the
coordinate axes in order to obtain a new coordinate system, in which
the observations become statistically isotropic. This assumes that
the estimates of the anisotropy ratio and the orientation angle
provided in anisPar
are accurate.
(i) A modified object with transformed coordinates if
rotateAnisotropicData is called with an Intamap object as input (see
intamap-package
) or (ii) the transformed coordinates if a
SpatialPointsDataFrame
is used as input or (iii)
the transformed coordinates if a SpatialPoints
object is the input.
Hristopulos Dionisis, Spiliopoulos Giannis
[1] Pebesma, E., Cornford, D., Dubois, G., Heuvelink, G.B.M., Hristopulos, D., Pilz, J., Stohlker, U., Morin, G., Skoien, J.O. INTAMAP: The design and implementation f an interoperable automated interpolation Web Service. Computers and Geosciences 37 (3), 2011.
[2] A. Chorti and D. T. Hristopulos (2008). Non-parametric Identification of Anisotropic (Elliptic) Correlations in Spatially Distributed Data Sets, IEEE Transactions on Signal Processing, 56(10), 4738-4751 (2008).
estimateAnisotropy
library(gstat)
data(sic2004)
coordinates(sic.val)=~x+y
sic.val$value=sic.val$dayx
anisPar <- estimateAnisotropy(sic.val)
print(anisPar)
rotatedObs <- rotateAnisotropicData(sic.val,anisPar)
newAnisPar <- estimateAnisotropy(rotatedObs)
print(newAnisPar)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.