Cannot load a shapefile in R "Error: no features found"

Cannot load a shapefile in R "Error: no features found"

我正在尝试读取您可以使用此 url 下载的 shapefile。

我有一个类似于下一个自动下载文件的代码：

library("raster")}
url<-"http://www6.gipuzkoa.eus/CATASTRO/Planos/ZIP-A098.zip"
downloader::download(url, dest=paste0(getwd(),"/","my_file.zip"), mode="wb",quiet=T)

zipped_shape_names<-c("098_HELBIDE_SHP/ATRIBUTOAK-A098.cpg","098_HELBIDE_SHP/ATRIBUTOAK-A098.dbf","098_HELBIDE_SHP/ATRIBUTOAK-A098.shp","098_HELBIDE_SHP/ATRIBUTOAK-A098.shx")

unzip("my_file.zip", files=zipped_shape_names)
my_shape<-raster::shapefile("098_HELBIDE_SHP/ATRIBUTOAK-A098.shp")

但是我得到的是以下错误：

Error in rgdal::readOGR(dirname(x), fn, stringsAsFactors = stringsAsFactors,  : 
  no features found
In addition: Warning messages:
1: In .local(x, ...) : .prj file is missing
2: In ogrFIDs(dsn = dsn, layer = layer) : no features found

您可以通过此link并按"Descargar planos"

访问原始网页

我在其他区域没有这个问题，只有这个和另一个，但我不知道这个特定区域发生了什么。

我们将不胜感激。

错误信息很清楚。 no features found 表示您的 shapefile 为空。您可以通过多种方式进行检查。

一种方法是将您的 shapefile 添加到 Q-Gis 或任何其他 GIS 软件工具中。对于 Q-GIS，您的 shapefile 将在图层窗格中弹出，但您看不到任何要素。

也可以查看R中的.dbf文件：

library(foreign)
read.dbf("098_HELBIDE_SHP/ATRIBUTOAK-A098.dbf")

.dbf 应包含与特征一样多的行。在你的情况下 None.

更简单的代码是

library(raster)
url<-"http://www6.gipuzkoa.eus/CATASTRO/Planos/ZIP-A098.zip"
download.file(url, dest="my_file.zip")
unzip("my_file.zip")   
s <- raster::shapefile("098_HELBIDE_SHP/ATRIBUTOAK-A098.shp")

显然该文件是空的。但是，它适用于此文件：

s <-raster::shapefile("098_LANDALUR_SHP/LANDALUR-PARTZELAK-A098.shp")