[R-sig-Geo] Distances in geoR package

Ernesto Jardim ernesto at ipimar.pt
Thu Jun 16 17:23:39 CEST 2005


Ruben Roa wrote:
>>-----Original Message-----
>>From: r-sig-geo-bounces at stat.math.ethz.ch
>>[mailto:r-sig-geo-bounces at stat.math.ethz.ch]On Behalf Of giovanna jona
>>lasinio
>>Sent: 16 June 2005 12:59
>>To: r-sig-geo at stat.math.ethz.ch
>>Subject: [R-sig-Geo] Distances in geoR package
>>
>>
>>Dear all, this may be a stupid question, but I wasn't able to find an
>>answer in the FAQ and in the mail archive.
>>I have a huge amount of  data, 644 sites for 30 different  
>>times, which coordinates are expressed in decimal longitude and latitude, 
>>I'd like to simply compute the empirical variogram and use "variofit" to estimate
>>variogram's parameters for each time. The "variog" function uses
>>euclidean distances and then it is not appropriate for lon-lat
>>coordinate. Is there any way to change the distance type in 
>>"variog"? Or which information are necessary to "variofit" in order to estimate
>>variogram's parameters?
>>
>>Thanks for any help
>>Giovanna Jona Lasinio
> 
> 
> I don't think you could do coordinate transformation in geostatistical packages.
> Ernesto Jardim (on a message sent to this list the 9th of June) mentioned a package 
> called spproj but today i could not find it in CRAN/bin/windows/contrib/2.0,2.1,2.2.
> For coordinate transformation i use GeoConv. It is free and its author,
> Eino Uikkanen, will be very glad to help you (me too). It works in batch mode 
> so you can send all your 30 data files and get the output in a few seconds.
> http://www.kolumbus.fi/eino.uikkanen/geoconvgb/
> If all your data are all in the same UTM zone, then GeoConv will automatically identify
> that zone.
> Ruben
> 
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo


Hi,

You can do transformations into UTM with PBSmapping. I never found out 
how to do it with spproj mainly because I was able to do it with 
PBSmapping :-)

Regards

EJ

ps: spproj is not in CRAN, look for it here 
http://r-spatial.sourceforge.net/




More information about the R-sig-Geo mailing list