[R-sig-Geo] finding CRS from latitude/longitude

Tom Gottfried tom.gottfried at o2mail.de
Sat Oct 6 15:51:27 CEST 2012


"+proj=longlat" can not be a valid CRS, as it only specifies that your 
coordinates are not projected. Probably you want to use the WGS84 
ellipsoid and datum:
http://www.spatialreference.org/ref/epsg/4326/

HTH,
Tom

Am 06.10.2012 14:48, schrieb alemante zd:
> Dear all,
>
> I would like to sp transform my latitude/longitude vector to by first
> telling it what CRS it is. This is what I did.
>
> F1<-data.frame(x, y)
> coordinates(F1)=~x+y
> SP<-SpatialPoints(F1,proj4string=CRS("+proj=longlat"))
> pts = SpatialPixels(SpatialPoints(F1))
> gridded(F1) = TRUE
>
> but I keep on getting the message that CRS is wrong.
> Coordinate Reference System (CRS) arguments: NA
>
> Could you please advice me on how I can get it right?
>
> cheers
> minti
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo



More information about the R-sig-Geo mailing list