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: > uniqueGPS(datgps) lat lon B-40 64.76 -21.59 D-01 66.01 -16.50 E-01 65.32 -14.36 G-04 64.43 -14.90 H-03 64.03 -16.98 M-02 66.10 -22.39 > SpatialPoints(uniqueGPS(datgps)) SpatialPoints: lat lon [1,] 64.76 -21.59 [2,] 66.01 -16.50 [3,] 65.32 -14.36 [4,] 64.43 -14.90 [5,] 64.03 -16.98 [6,] 66.10 -22.39 Coordinate Reference System (CRS) arguments: NA > SpatialPoints(uniqueGPS(datgps))->position > positionUTM <- spTransform(position, CRS("+proj=utm +zone=28 +north > +datum=WGS84")) Erreur dans spTransform(position, CRS("+proj=utm +zone=28 +north +datum=WGS84")) : No transformation possible from NA reference system >
What have I missed? Benjamin Le Lundi 30 Juin 2014 16:19 GMT, Barry Rowlingson <b.rowling...@lancaster.ac.uk> a écrit: > 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 <g...@hi.is> wrote: > > 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:12 GMT, Barry Rowlingson > > <b.rowling...@lancaster.ac.uk> a écrit: > > > >> On Mon, Jun 30, 2014 at 3:18 PM, Sarah Goslee <sarah.gos...@gmail.com> > >> 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] "SpatialPolygonsDataFrame" > >> > >> Use rgdal: > >> > >> require(rgdal) > >> Front = readOGR(".","ISL_adm0") > >> > >> then check Front has a proj4string. summary(Front) should mention it. > >> > >> This is getting to be the R-sig-geo equivalent of R FAQ 7.31.. > >> > >> Barry > > > > > > _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo