我正在尝试使用加载 Shapefile,readShapeSpatial()
但出现错误
“readShapeSpatial(shp) 中的错误:无法读取文件类型”。
关于这里出了什么问题的任何想法?这是我正在使用的 shapefile 的代码,它来自 2010 年巴西人口普查。
# Load Libraries
library(maptools)
library(rgdal)
# Set Working Directory
setwd("C:/Users/rafa/Desktop/test_shape")
# Download official shape files
download.file("ftp://geoftp.ibge.gov.br/malhas_digitais/censo_2010/setores_censitarios/rj.zip", "SetorCensitario_2010_Rio.zip", quiet = FALSE)
# Unzip files
unzip("SetorCensitario_2010_Rio.zip", exdir="SetorCensitario_2010_Rio")
# Load shapefile using `readShapeSpatial`
shp<- "C:/Users/rafa/Desktop/test_shape/SetorCensitario_2010_Rio/rj/33MUE250GC_SIR.shp"
shp <- readShapeSpatial(rioshape, proj4string = CRS("+proj=longlat +datum=WGS84"))
奇怪的是我可以使用readOGR
from package加载文件rgdal
shp <- readOGR(dsn="./SetorCensitario_2010_Rio/rj", layer="33MUE250GC_SIR")