[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