CommencezCommencez gratuitement

Fusionner des données provenant de différents SCR/projections

Chaque objet spatial est associé à un système de référence de coordonnées (SRC/CRS). En général, il est défini à l'importation des données et se lit directement dans les fichiers spatiaux. C'est ainsi que neighborhoods et nyc_tracts ont obtenu l'information sur leur système de coordonnées.

Les forfaits sp et raster offrent tous deux une fonction proj4string() qui renvoie le CRS de l'objet sur lequel elle est appelée.

Tenter de travailler avec des données spatiales utilisant des CRS différents, c'est un peu comme travailler avec un jeu de données en milles et un autre en kilomètres. Ils mesurent la même chose, mais les nombres ne sont pas directement comparables.

Jetons un coup d'œil à nos deux objets de polygones.

Cette activité fait partie du cours

Visualiser des données géospatiales avec R

Voir le cours

Instructions de l’exercice

  • Appelez proj4string() sur neighborhoods, puis de nouveau sur nyc_tracts. Vérifiez que les deux chaînes sont différentes.
  • Examinez le head() des coordinates() de neighborhoods et répétez pour nyc_tracts. Voyez-vous le problème ? nyc_tracts a des coordonnées x autour de -70, alors que neighborhoods est autour de 1 000 000 !
  • Tracez neighborhoods, puis tracez nyc_tracts avec col = "red" et add = TRUE pour les superposer.

Exercice interactif pratique

Essayez cet exercice en complétant ce code d’exemple.

library(sp)

# proj4string() on nyc_tracts and neighborhoods



# coordinates() on nyc_tracts and neighborhoods



# plot() neighborhoods and nyc_tracts
Modifier et exécuter le code