So beschriften Sie einen individuellen Status auf der Karte, während die anderen auf Unterteilung

Question

Ich habe es geschafft, eine Karte zu erstellen, aber ich muss ein Label für einen Staat (Ebene 2), die Unterteilungen (Ebene 3) enthält hinzufügen , anstatt jede Unterteilung zu kennzeichnen (nur für diesen Zustand). In den Daten "newpak" entsprechen Zeilen 641-664 diesem Zustand, gibt es eine Möglichkeit, nur einen Namen über diesen Zustand zu setzen.

library(dplyr) 
library(raster) 
library(sf) 
library(tidyverse) 
library(ggrepel) 
devtools::install_github("tidyverse/ggplot2", force = TRUE) 
library(ggplot2) 

pak <- getData("GADM",country="PAK",level=3) 

pak <- st_as_sf(pak) %>% 
    mutate(
    lon = map_dbl(geometry, ~st_centroid(.x)[[1]]), 
    lat = map_dbl(geometry, ~st_centroid(.x)[[2]])) 

ggplot(pak) + geom_sf() + geom_text(aes(label = NAME_3, x = lon, y = lat), size = 2) 

ind <- getData("GADM",country="IND",level=3) 

ind <- st_as_sf(ind) %>% 
    mutate(
    lon = map_dbl(geometry, ~st_centroid(.x)[[1]]), 
    lat = map_dbl(geometry, ~st_centroid(.x)[[2]])) 

jnk <- subset(ind, OBJECTID >= 641 & OBJECTID <= 664) 

newpak <- rbind(pak, jnk) 


regionalValues <- runif(165) # Simulate a value for each region between 0 and 1 

ggplot(newpak) + geom_sf(aes(fill = regionalValues)) + geom_text(aes(label = NAME_3, x = lon, y = lat), size = 2) 

Answer 1

Hier ist eine komplette Lösung, die sf Paket.

library(raster) 
library(sf) 
library(tidyverse) 

# downlaod PAK data and convert to sf 
pak <- getData("GADM",country="PAK",level=3) %>% 
    st_as_sf() 

# download IND data, convert to sf, filter out 
# desired area, and add NAME_3 label 
jnk <- getData("GADM",country="IND",level=3) %>% 
    st_as_sf() %>% 
    filter(OBJECTID %>% between(641, 664)) %>% 
    group_by(NAME_0) %>% 
    summarize() %>% 
    mutate(NAME_3 = "Put desired region name here") 


regionalValues <- runif(142) # Simulate a value for each region between 0 and 1 

# combine the two dataframes, find the center for each 
# region, and the plot with ggplot 
pak %>% 
    select(NAME_0, NAME_3, geometry) %>% 
    rbind(jnk) %>% 
    mutate(
    lon = map_dbl(geometry, ~st_centroid(.x)[[1]]), 
    lat = map_dbl(geometry, ~st_centroid(.x)[[2]]) 
    ) %>% 
    ggplot() + 
    geom_sf(aes(fill = regionalValues)) + 
    geom_text(aes(label = NAME_3, x = lon, y = lat), size = 2) + 
    scale_fill_distiller(palette = "Spectral") 

Einige Anmerkungen:

• benutzte ich sf::filter statt raster::subset die gewünschte Teilmenge der IND-Daten zu bekommen, weil ich es mehr idiomatische tidyverse Code fühlen.

• Um Bereiche kombinieren mit sf können Sie gruppieren die verschiedenen Regionen durch eine gemeinsame Gruppe mit group_by und dann einfach summarize nennen. Dies ist die Methode, die ich in meiner obigen Lösung verwendet habe. Es gibt andere Funktionen in dem Paket sf, die ähnliche Ergebnisse erzielen, die es wert sind, betrachtet zu werden. Sie sind st_combine und st_union.

• Mit st_centroid zum Zweck der Region Etiketten von der grafischen Darstellung ist nicht unbedingt die beste Methode für eine gute Lage für den Bereich Etiketten zu finden. Ich habe es benutzt, weil es am bequemsten ist. Sie können andere Methoden ausprobieren, einschließlich der manuellen Platzierung von Etiketten.

• Ich habe die Füllpalette in eine abweichende Farbpalette geändert, weil ich denke, dass sie den Unterschied zwischen einer Region und der nächsten deutlicher zeigt. Sie können einige der verfügbaren Farbpaletten mit RColorBrewer::display.brewer.all()