CommencerCommencer gratuitement

Fusionner des données avec des SCR/projections différents

Chaque objet spatial est associé à un système de coordonnées de référence (SCR). En général, il est défini lors de l’import des données et lu directement à partir des fichiers spatiaux. C’est ainsi que neighborhoods et nyc_tracts ont obtenu leurs informations de système de coordonnées.

Les packages sp et raster proposent une fonction proj4string() qui renvoie le SCR de l’objet sur lequel elle est appelée.

Essayer de travailler avec des données spatiales ayant des SCR différents, c’est un peu comme manipuler un jeu de données en miles et un autre en kilomètres. Ils mesurent la même chose, mais les valeurs ne sont pas directement comparables.

Jetons un œil à nos deux objets polygones.

Cet exercice fait partie du cours

Visualiser des données géospatiales avec R

Afficher le cours

Instructions

  • Appelez proj4string() sur neighborhoods, puis à 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, tandis 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 cet exemple de code.

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