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
