Load Packages

library(sf)
## Linking to GEOS 3.9.1, GDAL 3.2.1, PROJ 7.2.1
library(ggplot2)
library(tigris)
## To enable 
## caching of data, set `options(tigris_use_cache = TRUE)` in your R script or .Rprofile.
library(dplyr)
## 
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union

1. Jabar di Indonesia

Memasukkan Data SHP

Admin2Kabupaten<-"C:/r data spasial/Kabupaten/idn_admbnda_adm2_bps_20200401.shp"

Mengetahui file SHP di R

glimpse(Admin2Kabupaten)
##  chr "C:/r data spasial/Kabupaten/idn_admbnda_adm2_bps_20200401.shp"

Mengubah file

Admin2<-st_read(Admin2Kabupaten)
## Reading layer `idn_admbnda_adm2_bps_20200401' from data source 
##   `C:\r data spasial\Kabupaten\idn_admbnda_adm2_bps_20200401.shp' 
##   using driver `ESRI Shapefile'
## Simple feature collection with 522 features and 14 fields
## Geometry type: MULTIPOLYGON
## Dimension:     XY
## Bounding box:  xmin: 95.01079 ymin: -11.00762 xmax: 141.0194 ymax: 6.07693
## Geodetic CRS:  WGS 84

Melihat file baru

glimpse(Admin2)
## Rows: 522
## Columns: 15
## $ Shape_Leng <dbl> 2.360029, 1.963994, 4.590182, 3.287754, 4.448584, 4.907219,~
## $ Shape_Area <dbl> 0.22896809, 0.15413587, 0.23639581, 0.31616114, 0.34303826,~
## $ ADM2_EN    <chr> "Aceh Barat", "Aceh Barat Daya", "Aceh Besar", "Aceh Jaya",~
## $ ADM2_PCODE <chr> "ID1107", "ID1112", "ID1108", "ID1116", "ID1103", "ID1102",~
## $ ADM2_REF   <chr> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,~
## $ ADM2ALT1EN <chr> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,~
## $ ADM2ALT2EN <chr> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,~
## $ ADM1_EN    <chr> "Aceh", "Aceh", "Aceh", "Aceh", "Aceh", "Aceh", "Aceh", "Ac~
## $ ADM1_PCODE <chr> "ID11", "ID11", "ID11", "ID11", "ID11", "ID11", "ID11", "ID~
## $ ADM0_EN    <chr> "Indonesia", "Indonesia", "Indonesia", "Indonesia", "Indone~
## $ ADM0_PCODE <chr> "ID", "ID", "ID", "ID", "ID", "ID", "ID", "ID", "ID", "ID",~
## $ date       <date> 2019-12-20, 2019-12-20, 2019-12-20, 2019-12-20, 2019-12-20~
## $ validOn    <date> 2020-04-01, 2020-04-01, 2020-04-01, 2020-04-01, 2020-04-01~
## $ validTo    <date> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ geometry   <MULTIPOLYGON [°]> MULTIPOLYGON (((96.26836 4...., MULTIPOLYGON (~

Memasukkan Data via File CSV

JabarIndonesia<-read.csv("https://raw.githubusercontent.com/Nr5D/BatasWilayahID/main/JabarIndonesia.csv",
                         header=TRUE, sep=";")

Merge Data

merged_JabarIndonesia <- geo_join(spatial_data=Admin2, data_frame=JabarIndonesia, by_sp="ADM2_PCODE", by_df="ADM2_PCODE", how = "inner")
## Warning: We recommend using the dplyr::*_join() family of functions instead.

Atur Warna

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

Menampilkan Plot Peta pada Chloropleth

pDATA1<-ggplot()+
  geom_sf(data=merged_JabarIndonesia,aes(fill=DATA))+
  scale_fill_gradientn(colours=mycol)+
  labs(title="Jabar di Indonesia")
pDATA1

2. Jabar di Pulau Jawa

JabarJawa<-read.csv("https://raw.githubusercontent.com/Nr5D/BatasWilayahID/main/JabarJawa.csv", header=TRUE, sep=";")

Merge Data

merged_JabarJawa <- geo_join(spatial_data=Admin2, 
                             data_frame=JabarJawa, by_sp="ADM2_PCODE", by_df="ADM2_PCODE", how = "inner")
## Warning: We recommend using the dplyr::*_join() family of functions instead.

Menampilkan Plot Data

pDATA2<-ggplot()+
  geom_sf(data=merged_JabarJawa,aes(fill=DATA))+
  scale_fill_gradientn(colours=mycol)+
  labs(title="Jabar di Pulau Jawa")
pDATA2

3. Bogor di Jawa Barat

Memasukkan data penduduk Jawa Barat

BogorJabar<-read.csv("C:/r data spasial/Book3.csv", header=TRUE, sep=";")

Merge Data

merged_BogorJabar <- geo_join(spatial_data=Admin2, 
                             data_frame=BogorJabar, by_sp="ADM2_PCODE", by_df="ADM2_PCODE", how = "inner")
## Warning: We recommend using the dplyr::*_join() family of functions instead.

Menampilkan Plot Data

pDATA3<-ggplot()+
  geom_sf(data=merged_BogorJabar,aes(fill=DATA))+
  scale_fill_gradientn(colours=mycol)+
  labs(title="Bogor di Jawa Barat")
pDATA3

4. Kecamatan di Bogor

Memasukkan Data SHP

Admin3Kecamatan<-"C:/r data spasial/Kecamatan/idn_admbnda_adm3_bps_20200401.shp"

Mengetahui file SHP di R

glimpse(Admin3Kecamatan)
##  chr "C:/r data spasial/Kecamatan/idn_admbnda_adm3_bps_20200401.shp"

Mengubah file

Admin3<-st_read(Admin3Kecamatan)
## Reading layer `idn_admbnda_adm3_bps_20200401' from data source 
##   `C:\r data spasial\Kecamatan\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

Melihat file baru

glimpse(Admin3)
## Rows: 7,069
## Columns: 17
## $ Shape_Leng <dbl> 0.2798656, 0.7514001, 0.6900061, 0.6483629, 0.2437073, 1.35~
## $ Shape_Area <dbl> 0.003107633, 0.016925540, 0.024636382, 0.010761277, 0.00116~
## $ ADM3_EN    <chr> "2 X 11 Enam Lingkung", "2 X 11 Kayu Tanam", "Abab", "Abang~
## $ ADM3_PCODE <chr> "ID1306050", "ID1306052", "ID1612030", "ID5107050", "ID7471~
## $ ADM3_REF   <chr> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,~
## $ ADM3ALT1EN <chr> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,~
## $ ADM3ALT2EN <chr> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA,~
## $ ADM2_EN    <chr> "Padang Pariaman", "Padang Pariaman", "Penukal Abab Lematan~
## $ ADM2_PCODE <chr> "ID1306", "ID1306", "ID1612", "ID5107", "ID7471", "ID9432",~
## $ ADM1_EN    <chr> "Sumatera Barat", "Sumatera Barat", "Sumatera Selatan", "Ba~
## $ ADM1_PCODE <chr> "ID13", "ID13", "ID16", "ID51", "ID74", "ID94", "ID94", "ID~
## $ ADM0_EN    <chr> "Indonesia", "Indonesia", "Indonesia", "Indonesia", "Indone~
## $ ADM0_PCODE <chr> "ID", "ID", "ID", "ID", "ID", "ID", "ID", "ID", "ID", "ID",~
## $ date       <date> 2019-12-20, 2019-12-20, 2019-12-20, 2019-12-20, 2019-12-20~
## $ validOn    <date> 2020-04-01, 2020-04-01, 2020-04-01, 2020-04-01, 2020-04-01~
## $ validTo    <date> NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA~
## $ geometry   <MULTIPOLYGON [°]> MULTIPOLYGON (((100.2811 -0..., MULTIPOLYGON (~

Memasukkan Data Penduduk per Kecamatan di Bogor via File CSV

KecamatanBogor<-read.csv("C:/r data spasial/Book2.csv", header=TRUE, sep=";")

Merge Data

merged_KecamatanBogor <- geo_join(spatial_data=Admin3, data_frame=KecamatanBogor, by_sp="ADM3_PCODE", by_df="ADM3_PCODE", how = "inner")
## Warning: We recommend using the dplyr::*_join() family of functions instead.

Atur Warna

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

Menampilkan Plot Peta pada Chloropleth

pDATA4<-ggplot()+
  geom_sf(data=merged_KecamatanBogor,aes(fill=DATA))+
  scale_fill_gradientn(colours=mycol)+
  labs(title="Kecamatan di Bogor")
pDATA4