Peta per Kelurahan

Package

library(tidyverse) #ggplot2
library(sf)
library(viridis)
library(patchwork)

Data kecamatan

shp_smg <- st_read("datasets/KELURAHAN_AR/KELURAHAN_AR.shp")

data_smg_3_kec <- read_csv("datasets/data_smg_3_kec.csv")

shp_smg_kec <- shp_smg |> 
  summarise(geometry = st_union(geometry),
            .by = WADMKC)

data_smg_3_peta_kec <- left_join(shp_smg_kec,
                                 data_smg_3_kec)

Peta kecamatan

data_smg_3_peta_kec |> 
  ggplot(aes(fill = stunting_23)) +
  geom_sf() +
  scale_fill_viridis(direction = -1) +
  geom_sf_text(aes(label = WADMKC),
               size = 2)

ggsave("ggsave/peta_kecamatan_1.png", width = 7, height = 7)

data_smg_3_peta_kec |> 
  ggplot(aes(fill = prev_stu_22)) +
  geom_sf() +
  scale_fill_viridis(direction = -1) +
  geom_sf_text(aes(label = WADMKC),
               size = 2)

ggsave("ggsave/peta_kecamatan_2.png", width = 7, height = 7)

data_smg_3_peta_kec |> 
  ggplot(aes(fill = prev_stu_23)) +
  geom_sf() +
  scale_fill_viridis(direction = -1) +
  geom_sf_text(aes(label = WADMKC),
               size = 2)

ggsave("ggsave/peta_kecamatan_3.png", width = 7, height = 7)

p1 <- data_smg_3_peta_kec |> 
  ggplot(aes(fill = prev_stu_22)) +
  geom_sf() +
  scale_fill_viridis(direction = -1) +
  geom_sf_text(aes(label = WADMKC),
               size = 1) +
  theme_void() +
  theme(legend.position = "top")

p2 <- data_smg_3_peta_kec |> 
  ggplot(aes(fill = prev_stu_23)) +
  geom_sf() +
  scale_fill_viridis(direction = -1) +
  geom_sf_text(aes(label = WADMKC),
               size = 1) +
  theme_void() +
  theme(legend.position = "top")

p1 + p2

ggsave("ggsave/peta_kecamatan_4.png", width = 7, height = 5)