ComeçarComece de graça

Unindo dados de diferentes CRS/projeções

Todo objeto espacial tem um sistema de referência de coordenadas (CRS) associado a ele. Geralmente, isso é definido quando os dados são importados e é lido diretamente dos arquivos espaciais. É assim que neighborhoods e nyc_tracts obtiveram suas informações de sistema de coordenadas.

Os pacotes sp e raster têm a função proj4string() que retorna o CRS do objeto em que é chamada.

Tentar trabalhar com dados espaciais usando CRS diferentes é um pouco como trabalhar com um conjunto de dados em milhas e outro em quilômetros. Eles medem a mesma coisa, mas os números não são diretamente comparáveis.

Vamos dar uma olhada nos nossos dois objetos de polígonos.

Este exercício faz parte do curso

Visualizando dados geoespaciais em R

Ver curso

Instruções do exercício

  • Chame proj4string() em neighborhoods e depois em nyc_tracts. Verifique que as duas strings são diferentes.
  • Veja o head() de coordinates() de neighborhoods e repita para nyc_tracts. Consegue ver o problema? nyc_tracts tem coordenadas x por volta de -70, mas neighborhoods está por volta de 1.000.000!
  • Plote neighborhoods, depois plote nyc_tracts com col = "red" e add = TRUE para adicioná-los por cima.

Exercício interativo prático

Experimente este exercício completando este código de exemplo.

library(sp)

# proj4string() on nyc_tracts and neighborhoods



# coordinates() on nyc_tracts and neighborhoods



# plot() neighborhoods and nyc_tracts
Editar e executar o código