0

我正在尝试使用加载 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"))

奇怪的是我可以使用readOGRfrom package加载文件rgdal

shp <- readOGR(dsn="./SetorCensitario_2010_Rio/rj", layer="33MUE250GC_SIR")
4

0 回答 0