He leído un shapefile usandoreadShapePoly en elmaptools paquete, pero no puede leer el mismo archivo conreadOGR. Espero que alguien pueda ayudarme a leer el archivo de forma conreadOGR.

Descargué el archivoorcounty.shp de aquí:http://geography.uoregon.edu/geogr/topics/maps.htm

También descargué los archivos asociados:orcounty.shx, orcounty.sbx, orcounty.sbnyorcounty.dbf y poner los cinco archivos en la carpeta:c:/users/mark w miller/gis_in_R/shapefile_example/

El siguiente código lee el archivo de forma y muestra algunos atributos:

library(maptools)

setwd('c:/users/mark w miller/gis_in_R/shapefile_example/')

# Oregon county census data (polygons)
orcounty.poly <- readShapePoly('orcounty.shp', proj4string=CRS("+proj=longlat"))
orcounty.line <- readShapeLines('orcounty.shp', proj4string=CRS("+proj=longlat"))

# see projection
summary(orcounty.poly)

Object of class SpatialPolygonsDataFrame
Coordinates:
         min        max
x -124.55840 -116.46944
y   41.98779   46.23626
Is projected: FALSE 
proj4string : [+proj=longlat]
Data attributes:

Sin embargo, cuando intento leer el mismo archivo de forma usando el siguiente código, recibo un error:

library(rgdal)

# read shapefile
oregon.map <- readOGR(dsn="c:/users/mark w miller/gis_in_R/shapefile_example/", layer="orcounty")

# convert to dataframe
oregon.map_df <- fortify(oregon.map)

El mensaje de error dice:

Error in ogrInfo(dsn = dsn, layer = layer, encoding = encoding, use_iconv = use_iconv) : 
  Cannot open file

Puedo leer la tierra naturalhttp://www.naturalearthdata.com/ archivos de forma usando:

library(rgdal)

setwd("c:/users/mark w miller/gis_in_R/")

# read shapefile
wmap <- readOGR(dsn="ne_110m_physical", layer="ne_110m_land")

Entonces, aparentemente hay una diferencia entre los archivos de forma Natural Earth y el archivo de forma Oregonorcounty.shp.

Gracias por cualquier consejo sobre cómo leerorcounty.shp conreadOGR. Mi pregunta es similar a la pregunta aquí:rgdal / readOGR: no se puede leer el archivo de forma de .zip