Re: [R-sig-Geo] Conversion geographical coodinate in UTM

2014-07-01 Thread Gilles Benjamin Leduc
Dear R FAQ Team :p So I have been able to use the readOGR function, could anyone transfer to its maintainer a request to add writeOGR in the "see also" section of help? Anyway, after having done that for my shapefile, I tried to deal with the points I had to plot on, here is what I tried:

Re: [R-sig-Geo] Conversion geographical coodinate in UTM

2014-06-30 Thread Roger Bivand
On Mon, 30 Jun 2014, Barry Rowlingson wrote: Not your fault that readShapeSpatial ignores the projection! I think Roger recently spoke about fixing these things... The underlying shapelib code will not be maintained, anmd maptools::readShape* will be deprecated very soon. Always use rgdal::

Re: [R-sig-Geo] Conversion geographical coodinate in UTM

2014-06-30 Thread Barry Rowlingson
Not your fault that readShapeSpatial ignores the projection! I think Roger recently spoke about fixing these things... Barry On Mon, Jun 30, 2014 at 4:33 PM, Gilles Benjamin Leduc wrote: > Thanks Barry, > > That's it… Sorry if my question was too simple… I'm not a geograph, I'm a > molecular

Re: [R-sig-Geo] Conversion geographical coodinate in UTM

2014-06-30 Thread Gilles Benjamin Leduc
Thanks Barry, That's it… Sorry if my question was too simple… I'm not a geograph, I'm a molecular biologist… so I'm not used to these packages … I started from the simple script I had for normal maps… (and it is monday, brain isn't fully operational!) Benjamin Le Lundi 30 Juin 2014 15

Re: [R-sig-Geo] Conversion geographical coodinate in UTM

2014-06-30 Thread Barry Rowlingson
On Mon, Jun 30, 2014 at 3:18 PM, Sarah Goslee wrote: > Hi, > > At a guess, your Front dataset doesn't have a CRS set, so Probably because you've used readShapeSpatial which ignores projection info from shapefiles... >> Front<-readShapeSpatial("ISL_adm0.shp") >> class(Front) >> [1] "SpatialPolygo

Re: [R-sig-Geo] Conversion geographical coodinate in UTM

2014-06-30 Thread Sarah Goslee
Hi, At a guess, your Front dataset doesn't have a CRS set, so spTransform() doesn't know what to transform *from*. See the examples in ?spTransform for the whole process, including setting the CRS for your data. Sarah On Mon, Jun 30, 2014 at 10:10 AM, Gilles Benjamin Leduc wrote: > Hi, > > I n

[R-sig-Geo] Conversion geographical coodinate in UTM

2014-06-30 Thread Gilles Benjamin Leduc
Hi, I need to convert some geographical coordinate it UTM. First I have a shapefile: Front<-readShapeSpatial("ISL_adm0.shp") class(Front) [1] "SpatialPolygonsDataFrame" then… I tried this after googling: DataUTM <- spTransform(Front, CRS("+proj=utm +zone=27 +north +datum=WGS84")) Erreur dans