CommencezCommencez gratuitement

Créer un sous-ensemble des quartiers

Vous n'avez pas besoin de tous ces quartiers supplémentaires à New York; vous allez donc ne conserver que ceux du comté de New York. Vous savez déjà comment faire!

neighborhoods est un SpatialPolygonsDataFrame et, au chapitre 2, vous avez appris à créer un sous-ensemble en fonction d'une colonne dans l'emplacement data. L'idée clé était de créer un vecteur logique, puis de filtrer le SpatialPolygonsDataFrame comme un data frame.

Comment repérer les bons quartiers? Consultez :

head(neighborhoods@data)

CountyFIPS est un code numérique qui identifie le comté. Si vous trouvez le code du comté de New York, vous pourrez garder uniquement les lignes ayant cette valeur.

Cette activité fait partie du cours

Visualiser des données géospatiales avec R

Voir le cours

Instructions de l’exercice

  • L'objet nyc_tracts_merge contient aussi des codes de comté dans la colonne COUNTYFP. Utilisez unique() pour trouver la valeur correspondant au comté de New York.
  • Créez un sous-ensemble de neighborhoods en ajoutant un test logique vérifiant si neighborhoods$CountyFIPS a la bonne valeur.
  • Modifiez votre graphique pour utiliser manhat_hoods au lieu de neighborhoods.
  • Ajoutez une couche tm_text() en associant text à "NTAName".

Exercice interactif pratique

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

library(tmap)

# Find unique() nyc_tracts_merge$COUNTYFP


# Add logical expression to pull out New York County
manhat_hoods <- neighborhoods[___, ]

tm_shape(nyc_tracts_merge) +
  tm_fill(col = "estimate") +
  tm_shape(water) +
  tm_fill(col = "grey90") +
  # Edit to use manhat_hoods instead
  tm_shape(neighborhoods) +
  tm_borders() +
  # Add a tm_text() layer
Modifier et exécuter le code