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

2013-12-07 Thread Maurizio Marchi
Barcelona data. > Message: 2 > Date: Tue, 3 Dec 2013 11:56:07 +0100 > From: Marc Mar? Dell'Olmo > > To: R-sig-geo mailing list > > Subject: [R-sig-Geo] Problem converting shp ED50 to WGS84 with > spTransform

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 > 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 >>

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 > 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

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 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 wh

[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"