Praktikum 3b

Library

#loading package
library(sf)
library(ggplot2)
library(tigris)
library(tidyverse)

Membaca Data Spasial

#memasukkan data shp
Admin3Kecamatan <- "C:/datashp/idn_admbnda_adm3_bps_20200401.shp"

#ubah data shp
Admin3 <- st_read(Admin3Kecamatan)
Reading layer `idn_admbnda_adm3_bps_20200401' from data source 
  `C:\datashp\idn_admbnda_adm3_bps_20200401.shp' using driver `ESRI Shapefile'
Simple feature collection with 7069 features and 16 fields
Geometry type: MULTIPOLYGON
Dimension:     XY
Bounding box:  xmin: 95.01079 ymin: -11.00762 xmax: 141.0194 ymax: 6.07693
Geodetic CRS:  WGS 84

Memasukkan Data Jumlah Bencana Alam Banjir di Kota Bekasi

#memasukkan data csv
titikBanjirBekasi<-read.csv("C:/datashp/data_titik_banjir_bekasi.csv", header=TRUE, sep=";")

titikBanjirBekasi
   Shape_Leng  Shape_Area        ADM3_EN ADM3_PCODE ADM3_REF ADM3ALT1EN
1 0,282645428 0,001551537   Bantargebang  ID3275030       NA         NA
2  0,17084769 0,001206115   Bekasi Barat  ID3275060       NA         NA
3 0,197664394 0,001289219 Bekasi Selatan  ID3275050       NA         NA
4 0,162630726 0,001202486   Bekasi Timur  ID3275040       NA         NA
  ADM3ALT2EN     ADM2_EN ADM2_PCODE    ADM1_EN ADM1_PCODE   ADM0_EN ADM0_PCODE
1         NA Kota Bekasi     ID3275 Jawa Barat       ID32 Indonesia         ID
2         NA Kota Bekasi     ID3275 Jawa Barat       ID32 Indonesia         ID
3         NA Kota Bekasi     ID3275 Jawa Barat       ID32 Indonesia         ID
4         NA Kota Bekasi     ID3275 Jawa Barat       ID32 Indonesia         ID
        date  validOn validTo TitikBanjir
1 12/20/2019 4/1/2020      NA          13
2 12/20/2019 4/1/2020      NA           6
3 12/20/2019 4/1/2020      NA          15
4 12/20/2019 4/1/2020      NA           6
 [ reached 'max' / getOption("max.print") -- omitted 8 rows ]

Menggabungkan Data Spasial dan Data

#merge data
mergedata <- geo_join(spatial_data = Admin3, data_frame = titikBanjirBekasi, by_sp="ADM3_PCODE", by_df="ADM3_PCODE", how="inner")

#atur warna
mycol <- c("green","yellow","red","red4")

Menampilkan Plot Peta Kota Bekasi

#menampilkan plot peta
peta <- ggplot() +
        geom_sf(data=mergedata,aes(fill=TitikBanjir)) +
        scale_fill_gradientn(colours = mycol) +
        labs(title = "Jumlah Titik Banjir di Kota Bekasi")
peta

Plot Peta Kota Bekasi dan Sekitarnya

#memasukkan data csv
BekasiJakarta<-read.csv("C:/datashp/bekasi_jakarta.csv", header=TRUE, sep=";")

#merge data
mergedata <- geo_join(spatial_data = Admin3, data_frame = BekasiJakarta, by_sp="ADM3_PCODE", by_df="ADM3_PCODE", how="inner")

#atur warna
mycol <- c("green","yellow","red","red4")

#menampilkan plot peta
peta1 <- ggplot() +
        geom_sf(data=mergedata,aes(fill=TitikBanjir)) +
        scale_fill_gradientn(colours = mycol) +
        labs(title = "Jumlah Titik Banjir di Kota Bekasi")
peta1