Я пытаюсь прочитать шейп-файл, который вы можете скачать с помощью этого url.
У меня есть код, похожий на следующий, для автоматической загрузки файлов:
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
Вы можете получить доступ к исходной веб-странице, нажав ссылка и нажав «Descargar planos»
У меня нет этой проблемы с другими областями, только с этой и еще одной, но я не знаю, что происходит с этой конкретной областью.
Любая помощь будет оценена по достоинству.
Я загрузил форму в Q-gis. В шейп-файле просто нет функций (он пуст, нет данных), сообщение об ошибке gdal довольно четкое.
Вы также получаете предупреждение о том, что ваш файл проекции отсутствует - его следует добавить, чтобы получить правильную географическую привязку (например, WGS84 или метрическую проекцию?) Шейп-файла.
Сообщение об ошибке довольно ясное. no features found
означает, что ваш шейп-файл пуст. Вы можете проверить это несколькими способами.
Один из них - добавить свой шейп-файл в Q-Gis или любой другой программный инструмент ГИС. В случае Q-GIS ваш шейп-файл появится на панели слоев, но вы не увидите никаких функций.
Также можно проверить файл .dbf
в R:
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")
Я не хочу распаковывать всю папку, поэтому я создаю переменную под названием zipped_shape_names
.
очевидно, это неполный шейп-файл. Можете ли вы загрузить его в другой программный инструмент ГИС?