[R-sig-Geo] "Geographical CRS given to non-conformant data" for CRS of epsg:4326
Agustin Lobo
alobolistas at gmail.com
Fri Dec 16 13:24:08 CET 2011
Hi!
I' trying to read this shapefile that is just in lon,lat WGS84
> ogrInfo(dsn=datad,layer="world_adm0v2")
Source: "/media/Iomega_HDD/JACOB/ExIE2011DATA/ExFireBorneoDATA/world_adm0",
layer: "world_adm0v2"
Driver: ESRI Shapefile number of rows 209
Feature type: wkbPolygon with 2 dimensions
+proj=longlat +datum=WGS84 +no_defs
Number of fields: 3
name type length typeName
1 NAME 4 40 String
2 GMI_CNTRY 4 3 String
3 REGION 4 25 String
but I get:
> world <-readOGR(dsn=datad,layer="world_adm0v2")
OGR data source with driver: ESRI Shapefile
Source: "/media/Iomega_HDD/JACOB/ExIE2011DATA/ExFireBorneoDATA/world_adm0",
layer: "world_adm0v2"
with 209 features and 3 fields
Feature type: wkbPolygon with 2 dimensions
Error in validityMethod(as(object, superClass)) :
Geographical CRS given to non-conformant data
I've tried
> world <-readOGR(dsn=datad,layer="world_adm0v2",p4s=CRSargs(CRS("+init=epsg:4326")))
but get the same result
can't this CRS (epsg:4326)
+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +towgs84=0,0,0
be read by readOGR() ?
Thanks
Agus
More information about the R-sig-Geo
mailing list