Mappa coropletica

Ora che sai come disegnare i poligoni, mettiamo i tuoi poligoni su una mappa. Ricorda: sostituisci la chiamata a ggplot() con una a ggmap() e la chiamata originale a ggplot() va nell’argomento base_layer(), poi aggiungi il layer dei poligoni come al solito:

ggmap(corvallis_map_bw,
      base_layer = ggplot(ward_sales,
                          aes(lon, lat))) +
  geom_polygon(aes(group = group, fill = ward))

Provalo ora nella console!

Ops, qualcosa non torna. I distretti 1, 3 e 8 sembrano frastagliati e sbagliati. Che è successo? Parte dei confini dei distretti va oltre il bordo della mappa. A causa delle impostazioni predefinite di ggmap(), tutti i dati fuori mappa vengono scartati prima del plot, quindi alcuni bordi dei poligoni vengono eliminati e, quando i punti rimanenti vengono uniti, ottieni forme errate.

Nessun problema, c’è una soluzione: ggmap() offre alcuni argomenti per controllare questo comportamento. Gli argomenti extent = "normal" insieme a maprange = FALSE obbligano il grafico a usare l’intervallo dei dati invece dell’intervallo della mappa per definire i limiti del plotting.

Questo esercizio fa parte del corso

Visualizzare dati geospaziali in R

Visualizza il corso

Esercizio pratico interattivo

Prova a risolvere questo esercizio completando il codice di esempio.

# Fix the polygon cropping
ggmap(corvallis_map_bw, 
      base_layer = ggplot(ward_sales, aes(lon, lat))) +
  geom_polygon(aes(group = group, fill = ward))

Modifica ed esegui il codice

Questo esercizio fa parte del corso

Visualizzare dati geospaziali in R

IntermediárioNível de habilidade

4.6+

Inizia il corso gratis

We'll dive in by displaying some spatial data -- property sales in a small US town -- using ggplot2 and we'll introduce you to the ggmap package as a quick way to add spatial context to your plots. We'll talk about what makes spatial data special and introduce you to the common types of spatial data we'll be working with throughout the course.

Exercise 1: Introduzione ai dati spaziali Exercise 2: Recuperare una mappa di sfondo Exercise 3: Mettere tutto insieme Exercise 4: Intuizioni attraverso le estetiche Exercise 5: Opzioni utili di get_map() e ggmap()Exercise 6: Mappe diverse Exercise 7: Sfruttare i punti di forza di ggplot2 Exercise 8: Un'alternativa rapida Exercise 9: Tipi comuni di dati spaziali Exercise 10: Disegnare poligoni Exercise 11: Mappa coropletica

Esercizio in corso

Exercise 12: Dati raster come heatmap

You can get a long way with spatial data stored in data frames, but it makes life easier if they are stored in special spatial objects. In this chapter we'll introduce you to the spatial object classes provided by the sp package, particularly for point and polygon data. You'll learn how to explore and subset these objects by exploring a world map. The reward for learning about these object classes: we'll show you the package tmap which requires spatial objects as input, but makes creating maps really easy! You'll finish up by making a map of the world's population.

Exercise 1: Introducing sp objects Exercise 2: Let's take a look at a spatial object Exercise 3: What's inside a spatial object?Exercise 4: A more complicated spatial object Exercise 5: sp and S4 Exercise 6: Walking the hierarchy Exercise 7: Further down the rabbit hole Exercise 8: More sp classes and methods Exercise 9: Subsetting by index Exercise 10: Accessing data in sp objects Exercise 11: Subsetting based on data attributes Exercise 12: tmap, a package that works with sp objects Exercise 13: Introduction to tmap Exercise 14: Building a plot in layers Exercise 15: Why is Greenland so big?Exercise 16: Saving a tmap plot

While the sp package provides some classes for raster data, the raster package provides more useful classes. You'll be introduced to these classes and their advantages and then learn to display them. The examples continue with the theme of population from Chapter 2, but you'll look at some much finer detail datasets, both spatially and demographically. In the second half of the chapter you'll learn about color -- an essential part of any visual display, but especially important for maps.

Exercise 1: The raster package Exercise 2: What's a raster object?Exercise 3: Some useful methods Exercise 4: A more complicated object Exercise 5: A package that uses Raster objects Exercise 6: Color scales Exercise 7: Pick the right palette Exercise 8: Adding a custom continuous color palette to ggplot2 plots Exercise 9: Custom palette in tmap Exercise 10: More about color scales Exercise 11: An interval scale example Exercise 12: A diverging scale example Exercise 13: A qualitative example

In this chapter you'll follow the creation of a visualization from raw spatial data files to adding a credit to a map. Along the way, you'll learn how to read spatial data into R, more about projections and coordinate reference systems, how to add additional data to a spatial object, and some tips for polishing your maps.

Exercise 1: Reading in spatial data Exercise 2: Reading in a shapefile Exercise 3: Reading in a raster file Exercise 4: Getting data using a package Exercise 5: Coordinate reference systems Exercise 6: Merging data from different CRS/projections Exercise 7: Converting from one CRS/projection to another Exercise 8: Adding data to spatial objects Exercise 9: The wrong way Exercise 10: Checking data will match Exercise 11: Merging data attributes Exercise 12: A first plot Exercise 13: Polishing a map Exercise 14: Subsetting the neighborhoods Exercise 15: Adding neighborhood labels Exercise 16: Tidying up the legend and some final tweaks Exercise 17: Wrap up