Не удается открыть шейп-файл в R - PullRequest
Новая
       6

Не удается открыть шейп-файл в R

0 голосов
/ 26 ноября 2018

Я не уверен, почему, но я вообще не могу читать шейп-файлы в r!

На MacOS Mojave,

Я пробовал: 

afg <- readOGR("afg.shp")
afg <- readOGR(dsn = ".", layer = "afg")

Шейп-файлыс этого сайта: http://www.diva -gis.org / gdata , дороги Афганистана.

Эта ошибка повторяется: 

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

Однако этот код работает: 

countries <- readOGR("countries.geojson")

Может кто-нибудь выяснить, почему?

1 Ответ

0 голосов
/ 26 ноября 2018

Я нахожусь на Мохаве, и это прекрасно работает: 

httr::GET(
  url = "http://biogeo.ucdavis.edu/data/diva/rds/AFG_rds.zip",
  httr::write_disk("AFG_rds.zip"),
  httr::progress()
) -> res

fils <- unzip("AFG_rds.zip", exdir = "afg-roads")

fils
## [1] "afg-roads/AFG_roads.dbf" "afg-roads/AFG_roads.prj" 
## [3] "afg-roads/AFG_roads.shp" "afg-roads/AFG_roads.shx"

rgdal::readOGR("afg-roads/AFG_roads.shp")
## OGR data source with driver: ESRI Shapefile 
## Source: "/Users/bob/Development/afg-roads/AFG_roads.shp", layer: "AFG_roads"
## with 1409 features
## It has 5 fields
## class       : SpatialLinesDataFrame 
## features    : 1409 
## extent      : 60.6075, 74.67168, 29.38582, 38.48983  (xmin, xmax, ymin, ymax)
## coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 
## variables   : 5
## names       :     MED_DESCRI,    RTT_DESCRI, F_CODE_DES, ISO,  ISOCOUNTRY 
## min values  :        Unknown, Primary Route,       Road, AFG, AFGHANISTAN 
## max values  : Without Median,       Unknown,      Trail, AFG, AFGHANISTAN
...