# # ================================ # Load required packages # ================================ # install.packages(c("dplyr", "purrr", "readr", "ggmap", "tidygeocoder")) library(dplyr) library(purrr) library(readr) library(ggmap) library(tidygeocoder) # ================================ # CONFIGURATION # ================================ use_google <- FALSE # Set to TRUE to use Google API; FALSE to use OpenStreetMap (free) google_api_key <- "YOUR_GOOGLE_API_KEY" # Replace with your actual Google API key # ================================ # Load earthquake dataset # ================================ df <- read_csv("F:/RR/earthquakeM.csv") # Change to your file path if needed # Filter earthquakes with magnitude ≥ 5.5 df_filtered <- df %>% filter(magnitude >= 5.5) # Sample 500 rows (you can adjust this) df_sample <- df_filtered[1:500, ] # ================================ # Geocoding section # ================================ if (use_google) { # Register Google API key register_google(key = google_api_key) # Reverse geocode using Google Maps API df_sample$full_address <- map_chr(1:nrow(df_sample), function(i) { tryCatch({ Sys.sleep(0.2) # Rate limit delay revgeocode( location = c(df_sample$longitude[i], df_sample$latitude[i]), output = "address", override_limit = TRUE ) }, error = function(e) { message(sprintf("Error at row %d: %s", i, e$message)) return(NA) }) }) # Extract country from address df_sample$country <- sapply(strsplit(df_sample$full_address, ", "), function(x) { if (is.na(x[1])) return(NA) tail(x, 1) }) } else { # Reverse geocode using OpenStreetMap via tidygeocoder df_sample <- df_sample %>% reverse_geocode(lat = latitude, long = longitude, method = 'osm', full_results = TRUE) # Extract country directly df_sample$country <- df_sample$country } # ================================ # Analyze results # ================================ # Count earthquakes by country country_counts <- df_sample %>% filter(!is.na(country)) %>% group_by(country) %>% summarise(earthquake_count = n(), .groups = "drop") %>% arrange(desc(earthquake_count)) # View top 10 countries print(head(country_counts, 20)) # Check Japan's rank japan_rank <- which(country_counts$name == "Japan") cat("Japan's rank:", japan_rank, "\n")