[R] conversion of LL coordenates to UTM problems (ED50-WGS84 format)

From: Berta <ibanez_at_bioef.org>
Date: Wed 18 Oct 2006 - 10:27:21 GMT


Hi R-Users,
I have plotted a region whose polygon coordinates are given in shp format ED50 UTM (zone=30) ) using "readShapePoly" in library(maptools).

Now I need to plot a set of points in that region (my.dataframe, with X and Y geographic coordinates), which have been read using GPS in Longitud-Latitud form (using WGS84 system), so I first need to convert these Longitud-Latitud data into UTM data. To do so, I use package PBSmapping:

attr(my.dataframe, "projection") <-"LL" attr(my.dataframe, "zone") <- 30
my.dataframe<-as.PolySet(my.dataframe)
my.dataframeUL<-convUL(my.dataframe)

obtaining the UTM coordinates. But when I plot them in the polygon region, points seem to be badly located. I suspect it is due to the ED50 UTM of the polygon versus WGS84 coordenates of the points. If this is the problem, is there a way to convert them into the same format? It it is not the problem, what am I doing wrong?

Thanks a lot in advance.
Berta.

        [[alternative HTML version deleted]]



R-help@stat.math.ethz.ch 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. Received on Wed Oct 18 20:32:23 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Wed 18 Oct 2006 - 11:30:11 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.