[R-sig-Geo] Problem converting shp ED50 to WGS84 with spTransform and writeOGR

2013-12-03 Thread Marc Marí Dell'Olmo
Dear all, I'm trying to convert a map of Barcelona (in shp format) with reference ED50 to a kml with Google Maps reference. I have used the following syntax, but it creates a slightly displaced map (maybe just a few meters) and do not understand why! setwd(C:/dir) barcelona - readOGR(C:/dir,

Re: [R-sig-Geo] Problem converting shp ED50 to WGS84 with spTransform and writeOGR

2013-12-03 Thread Roger Bivand
On Tue, 3 Dec 2013, Marc Marí Dell'Olmo wrote: Dear all, I'm trying to convert a map of Barcelona (in shp format) with reference ED50 to a kml with Google Maps reference. I have used the following syntax, but it creates a slightly displaced map (maybe just a few meters) and do not understand

Re: [R-sig-Geo] Problem converting shp ED50 to WGS84 with spTransform and writeOGR

2013-12-03 Thread Marc Marí Dell'Olmo
Dear Roger, thank you very much for your answer, but the problem persists... Now I'm using R3.0.1 and I have updated rgdal. Therefore now I see: CRS(+init=epsg:23031) CRS arguments: +init=epsg:23031 +proj=utm +zone=31 +ellps=intl +towgs84=-87,-98,-121,0,0,0,0 +units=m +no_defs I have tried

Re: [R-sig-Geo] Problem converting shp ED50 to WGS84 with spTransform and writeOGR

2013-12-03 Thread Marc Marí Dell'Olmo
In my previous email I was wrong, the proj4string is working ok! sorry... 2013/12/3 Roger Bivand roger.biv...@nhh.no On Tue, 3 Dec 2013, Marc Marí Dell'Olmo wrote: Dear all, I'm trying to convert a map of Barcelona (in shp format) with reference ED50 to a kml with Google Maps reference.

Re: [R-sig-Geo] Problem converting shp ED50 to WGS84 with spTransform and writeOGR

2013-12-03 Thread Marc Marí Dell'Olmo
Thank you very much for your answer. Now it works! Marc 2013/12/3 Roger Bivand roger.biv...@nhh.no On Tue, 3 Dec 2013, Marc Marí Dell'Olmo wrote: Dear all, I'm trying to convert a map of Barcelona (in shp format) with reference ED50 to a kml with Google Maps reference. I have used the

[R-sig-Geo] 3rd+ rasters added using plot(r, add=TRUE) are misregistered

2013-12-03 Thread Josh O'Brien
Hi all, Briefly, running the following code will not produce a single grey ring, as would be expected. The first three raster layers all stack up neatly, but the 4th and subsequent raster layers are all shifted/skewed to the right. It appears that the plotting engine loses track of the fact that

[R-sig-Geo] rasterEngine at AGU

2013-12-03 Thread Jonathan Greenberg
If any of you are headed to the American Geophysical Union in San Francisco next week, please drop by my poster and say hi! I'll be manning a poster about rasterEngine on Monday from 1:40pm to 6pm: B13D-0525. rasterEngine: an easy-to-use R function for applying complex geostatistical models to