imicola schreef:
Hi,
I read somewhere that when carrying out geostatistical analysis in R you
should not use latitude and longitude...can anyone expand on this a little
for me, and what would be the best coordinate system to use?
I have my data in a geographic coordinate system, WGS84, decimal
degrees....is this the wrong format for such analyses?
I have also converted my data in the UTM projection and so have it in
metres....(ranging from 480,000 to 550,000 E and 170,000 to 230,000 N).
If I was to use the UTM coordinates, should I be using the actual
coordinates in metres, or should I convert this into an arbitrary coordinate
system (i.e. from 0 - 1) somehow?
I have noticed that running an analysis on the data gives different results
depending on which type of system you use, so I want to make sure I have the
correct system. I should also probably note that I am a geostatistical
novice!
Thanks,
Hi,
I use the gstat package for geostatistics. For doing the analysis I
don't think it is necessary to convert to UTM. But maybe just do it to
be on the safe side. If you use the spatial objects provided by the
sp-package (http://cran.r-project.org/web/packages/sp/vignettes/sp.pdf)
you transform your data to other projections using the spTransform package.
Questions regarding geostatistics and spatial data will result in more
answers on the r-sig-geo list.
cheers,
Paul
--
Drs. Paul Hiemstra
Department of Physical Geography
Faculty of Geosciences
University of Utrecht
Heidelberglaan 2
P.O. Box 80.115
3508 TC Utrecht
Phone: +31302535773
Fax: +31302531145
http://intamap.geo.uu.nl/~paul
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
