Conversion geographical coodinate in UTM
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.rowlingson at 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 <gbl1 at 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.rowlingson at lancaster.ac.uk> a ?crit:
On Mon, Jun 30, 2014 at 3:18 PM, Sarah Goslee <sarah.goslee at 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