web-dev-qa-db-fra.com

Lire le fichier de forme avec les versets readOGR readShapePoly

J'ai lu un fichier de formes utilisant readShapePoly dans le package maptools, mais je ne peux pas lire le même fichier avec readOGR. J'espère que quelqu'un pourra m'aider à lire le fichier de formes avec readOGR.

J'ai téléchargé le fichier orcounty.shp à partir d'ici: http://geography.uoregon.edu/geogr/topics/maps.htm

J'ai également téléchargé les fichiers associés: orcounty.shx, orcounty.sbx, orcounty.sbn et orcounty.dbf et ai placé les cinq fichiers dans le dossier: c:/users/mark w miller/gis_in_R/shapefile_example/.

Le code suivant lit le fichier de formes et affiche certains attributs:

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:

Cependant, lorsque j'essaie de lire le même fichier de formes à l'aide du code suivant, un message d'erreur s'affiche:

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)

Le message d'erreur dit:

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

Je peux lire Natural Earth http://www.naturalearthdata.com/ shapefiles en utilisant:

library(rgdal)

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

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

Donc, apparemment, il existe une différence entre les fichiers de formes Natural Earth et les fichiers de formes Oregon orcounty.shp.

Merci pour tout conseil sur la façon de lire orcounty.shp avec readOGR. Ma question est semblable à la question ici: rgdal/readOGR - impossible de lire le fichier de formes à partir de .Zip

rshapefilergdal
12
Mark Miller

Essayez de supprimer votre dernier '/' du chemin du fichier.

readOGR(dsn = 'c:/users/mark w miller/gis_in_R/shapefile_example',
        layer = 'orcounty')
13
Paulo E. Cardoso

Pour tous ceux qui se retrouvent ici avec cette erreur sur une machine Linux, j’ai trouvé que le problème était l’utilisation d’un raccourci vers le chemin principal. c'est à dire.

# Works
readOGR(dsn="/home/user/dir", layer="file")

# Doesn't work
readOGR(dsn="~/dir", layer="file")

Je ne sais pas pourquoi.

3
MikeRSpencer

J'ai utilisé le fichier ne_110m_land

Essayez avec ceci:

setwd('D:/JMSR/codes.R/mapas')
unzip("ne_110m_land.Zip")
ogrInfo(".", "ne_110m_land")
wmap <- readOGR(".", "ne_110m_land")
1
EdJo1924

raster::shapefile encapsule readOGR pour prendre en charge les chemins et les tildes; il suffit de passer le nom de fichier complet.

 library(raster)
 x <- shapefile("c:/users/orcounty.shp')

ou 

 y <- shapefile("~/users/orcounty.shp")
0
Robert Hijmans