我正在尝试使用readShapeSpatial()
加载Shapefile,但我收到错误
“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
包中的rgdal
加载文件
shp <- readOGR(dsn="./SetorCensitario_2010_Rio/rj", layer="33MUE250GC_SIR")