是否有一种将ShapeFile直接从在线来源中读取为R中的方法？

Question

我正在尝试找到一种从在线存储库/文件夹/url中加载ShapeFiles（.shp）的方法，直接在R中加入我的全局环境，目的是使用geom_sf在GGPLOT2中制作图。首先，我使用Google驱动器来存储这些文件，但理想情况下，我希望找到一个解决方案，该解决方案与具有有效URL且适当访问权限的任何文件夹一起使用。

到目前为止，我已经尝试了一些选项，前两个涉及在Google Drive上汇总源文件夹，然后存储ShapeFiles然后以某种方式下载和解压缩。已包含使用小型测试shapefile的可复制示例：

1. 使用utils :: download.file（）用于检索压缩文件夹并使用base base base :: System（'unzip ..'）<）< /code>或zip :: unzip（）（在此线程中宽松地遵循：从ONS下载County Shapefile ）：

# Create destination data folder (if there isn't one)
if(!dir.exists('data')) dir.create('data')
# Download the zipped file/folder 
download.file("https://drive.google.com/file/d/1BYTCT_VL8EummlAsH1xWCd5rC4bZHDMh/view?usp=sharing", destfile = "data/test_shp.zip")
# Unzip folder using unzip (fails)
unzip(zipfile = "data/test_shp.zip", exdir = "data/test_shp", junkpaths = TRUE)
# Unzip folder using system (also fails)
system("unzip data/test_shp.zip")

如果您无法运行上述代码，则FYI，2个错误消息是：

警告消息： 在unzip（zipfile =“ data/test_shp.zip”，exdir =“ data/test_shp”，： 从zip文件提取

和

中央目录签名中的错误1。这个文件不是 Zipfile或它构成一个多部分档案的一个磁盘。在 后一种情况将在中央目录和Zipfile评论中找到 该存档的最后一个磁盘。 UNZIP：在数据/test_shp.zip或 data/test_shp.zip.zip。

data/test_shp.zip.zip，并且找不到 代码> download.file（）步骤。

2. 使用GoogleRive软件包：

# Create destination data folder (if there isn't one)
if(!dir.exists('data')) dir.create('data')
# Specify googledrive url:
test_shp = drive_get(as_id("https://drive.google.com/file/d/1BYTCT_VL8EummlAsH1xWCd5rC4bZHDMh/view?usp=sharing"))
# Download zipped folder
drive_download(test_shp, path = "data/test_shp.zip")
# Unzip folder
zip::unzip(zipfile = "data/test_shp.zip", exdir = "data/test_shp", junkpaths = TRUE)
# Load test.shp
test_shp <- read_sf("data/test_shp/test.shp")

这可以！

...除了它仍然是一个骇客的解决方法，它要求我zip，下载，解压缩，然后使用单独的功能（例如sf :: read_sf 或st_read）将数据读取到我的全球环境中。而且，由于它使用GoogleRive软件包，因此仅适用于存储在此系统中的文件（不是OneDrive，Dropbox和其他URL）。

3. 我还尝试了sf :: read_sf，st_read和fastshp :: read.shp直接在文件夹URL上一个人可能期待。

因此，我的问题是：是否有一个直接在网上存储在R中的ShapeFile的工作流程，还是我应该停止寻找？如果没有，但是有一种方法可以扩展我的上述解决方案（2）超越Googledrive，我也很感谢它的任何提示！

注意：我还应该补充说，我故意忽略了任何需要包装rgdal的选择在某个时候掉下地图）。提前致谢！

Answer 1

我最近遇到了一个类似的问题，必须直接从Dropbox中读取ShapeFiles到R中。

结果，此解决方案仅适用于Dropbox的情况。

您需要做的第一件事是使用RDROP2为Dropbox创建一个可刷新的令牌，鉴于Dropbox的最新更改将单个令牌使用限制为4小时。您可以关注此 a>。

设置刷新令牌后，请使用以下方式识别Dropbox上空间数据文件夹中的所有文件：

shp_files_on_db<- drop_dir("Dropbox path/to your/spatial data/", dtoken = refreshable_token) %>% 
  filter(str_detect(name, "adm2"))

我的“空间数据”文件夹包含两组ShapeFiles -ADM1和ADM 2。我使用上述代码仅选择与之相关的代码ADM2。

然后创建SHP，CSV，SHX，DBF，CPG文件的名称的向量，如下所示：

shp_filenames<- shp_files_on_db$name

我选择将ShapeFiles读入临时目录中，避免需要存储文件。在我的磁盘上 - 也可用于闪亮的实现。我创建此临时目录如下：

# create a new directory under tempdir
dir.create(dir1 <- file.path(tempdir(), "testdir"))

#If needed later on, you can delete this temporary directory 
unlink(dir1, recursive = T)

#And test that it no longer exists
dir.exists(dir1)

现在将Dropbox文件下载到此临时目录：

for (i in 1: length(shp_filenames)){
  
  drop_download(paste0("Dropbox path/to your/spatial data/",shp_filenames[i]),
                dtoken = refreshable_token,
                local_path = dir1)
  
}

最后，在您的shapefile中阅读如下：

#path to the shapefile in the temporary directory
path1_shp<- paste0(dir1, "/myfile_adm2.shp")

#reading in the shapefile using the sf package - a recommended replacement for rgdal
shp1a <- st_read(path1_shp)