Pre-Informe Asesoría 1

Como parte de la caracterización descriptiva, se generaron visualizaciones espaciales de los sismos mediante mapas interactivos para la zona transformante de California y para el margen de subducción de Chile, en función de sus magnitudes, profundidades y otras características. Este enfoque multivariante permite identificar visualmente las diferencias estructurales subyacentes entre ambos regímenes tectónicos.

# LIBRERIAS ----
library(readr)
library(leaflet)
library(dplyr)
library(maps)
library(htmltools)

Para California se tiene:

mapa_ca <- leaflet(california) %>%
  addProviderTiles(providers$CartoDB.Positron) %>%
  addCircleMarkers(
    lng = ~longitude,
    lat = ~latitude,
    radius = ~pmax((mag - 4) * 4, 3),
    color = ~pal_depth_ca(depth),
    stroke = TRUE,
    weight = 1,
    opacity = 1,
    fillOpacity = 0.75,
    popup = popup_ca
  ) %>%
  addLegend(
    position = "topright",
    pal = pal_depth_ca,
    values = ~depth,
    title = "Profundidad (km)",
    opacity = 1
  ) %>%
  addControl(
    leyenda_mag_ca,
    position = "topright"
  ) %>%
  setView(lng = -119.5, lat = 37.2, zoom = 6)

mapa_ca

Y para Chile:

mapa_chile_leaflet <- leaflet(chile) %>%
  addProviderTiles(providers$CartoDB.Positron) %>%
  addCircleMarkers(
    lng = ~longitude,
    lat = ~latitude,
    radius = ~pmax((mag - 4) * 4, 3),
    color = ~pal_depth_chile(depth),
    stroke = TRUE,
    weight = 1,
    opacity = 1,
    fillOpacity = 0.75,
    popup = popup_chile
  ) %>%
  addLegend(
    position = "topright",
    pal = pal_depth_chile,
    values = ~depth,
    title = "Profundidad (km)",
    opacity = 1
  ) %>%
  addControl(
    leyenda_mag_chile,
    position = "topright"
  ) %>%
  setView(lng = -71, lat = -35, zoom = 4)

mapa_chile_leaflet