Lendo um shapefile
Shapefiles são uma das formas mais comuns de compartilhar dados espaciais e podem ser facilmente lidos no R usando readOGR() do pacote rgdal. readOGR() tem dois argumentos importantes: dsn e layer. Exatamente o que você passa para esses argumentos depende do tipo de dado que você está lendo. Você viu no vídeo que, para shapefiles, dsn deve ser o caminho para o diretório que contém os arquivos que compõem o shapefile e layer é o nome do arquivo do shapefile específico (sem extensão).
Para o seu mapa, você quer os limites dos bairros. Nós baixamos as Áreas de Tabulação de Bairros (Neighborhood Tabulation Areas), conforme definidas pela cidade de Nova York, da Plataforma de Dados Abertos do Department of City Planning. O download veio em um arquivo zip e colocamos o resultado da extração do arquivo baixado no seu diretório de trabalho.
Você vai usar a função dir() do R base para examinar o conteúdo do seu diretório de trabalho e, em seguida, ler o shapefile no R.
Este exercício faz parte do curso
Visualizando dados geoespaciais em R
Instruções do exercício
- Use
dir()sem argumentos para descobrir o nome do diretório do shapefile. - Use
dir(), passando o caminho para o diretório do shapefile, para ver os arquivos dentro dele. - Agora que você sabe o diretório e o nome do arquivo, use
readOGR()para ler o shapefile dos bairros em um objeto chamadoneighborhoods. - Verifique o conteúdo chamando
summary()emneighborhoods. - Verifique o conteúdo plotando
neighborhoods.
Exercício interativo prático
Experimente este exercício completando este código de exemplo.
library(sp)
library(rgdal)
# Use dir() to find directory name
# Call dir() with directory name
# Read in shapefile with readOGR(): neighborhoods
# summary() of neighborhoods
# Plot neighborhoods