1.Load Package

library(sp)

## Warning: package 'sp' was built under R version 4.1.3

library(magrittr)

## Warning: package 'magrittr' was built under R version 4.1.3

library(leaflet)

## Warning: package 'leaflet' was built under R version 4.1.3

Membaca data

library(readxl)

## Warning: package 'readxl' was built under R version 4.1.3

quakes <- read_excel("datagempa.xlsx")

head(quakes)

Melihat Plot Data

plot(quakes,     pch=20, col= rgb(0,0,0,.1) , lwd=.6)

2.Membuat plot kejadian gempa bumi

LAT <-as.numeric(quakes$lat)
require(ggplot2)

## Loading required package: ggplot2

## Warning: package 'ggplot2' was built under R version 4.1.3

qplot(data = quakes, x = long, y = LAT, size = exp(mag), color = mag, alpha = I(.8))

str(quakes)

## tibble [707 x 4] (S3: tbl_df/tbl/data.frame)
##  $ long : num [1:707] 111 113 113 110 111 ...
##  $ lat  : chr [1:707] "-10.6" "-10.09" "-10.039999999999999" "-9.5399999999999991" ...
##  $ depth: num [1:707] 40 71 14 10 9 13 40 121 166 563 ...
##  $ mag  : num [1:707] 4 3.2 3.2 3 2.75 2.8 3.8 3.6 4.2 4.9 ...

3.Membuat Plot Gempa Bumi dengan Basemap leaflet

leaflet(data = quakes[1:707,]) %>% addTiles() %>%
  addMarkers(~long, ~LAT, popup = ~as.character(mag), label = ~as.character(mag))

4.Referensi:

https://api.rpubs.com/ellysalutfi/525171 Data: BMKG Stasiun Geofisika Nganjuk