Question

我正在尝试读取可通过此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”来访问原始网页

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

任何帮助将不胜感激。

Answer 1

错误消息很清楚。 DataTabe dt = objProfitLossDT.Select("AppBalance <= 0");表示您的shapefile为空。您可以通过多种方式进行检查。

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

还可以检查R中的no features found文件：

.dbf

library(foreign) read.dbf("098_HELBIDE_SHP/ATRIBUTOAK-A098.dbf")应该包含的行数与要素的数量一样多。就您而言，.dbf。

Answer 2

更简单的代码是

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")

无法在R中加载shapefile“错误：找不到特征”

2 个答案: