# Cargar datos desde archivo Excel (ya con comportamiento asignado)
library(readxl)
estaciones_diag <- read_excel("datos_climas_huila.xlsx", sheet = "diagnostico_estaciones")

# Asegurarse de que las coordenadas sean numéricas
estaciones_diag <- estaciones_diag %>%
  mutate(
    latitud_estacion = as.numeric(latitud_estacion),
    longitud_estacion = as.numeric(longitud_estacion)
  )
leaflet(estaciones_diag) %>%
  addTiles() %>%
  addCircleMarkers(
    lng = ~longitud_estacion,
    lat = ~latitud_estacion,
    color = ~case_when(
      comportamiento == "Problema" ~ "red",
      comportamiento == "Posible problema" ~ "orange",
      comportamiento == "Estable" ~ "green",
      TRUE ~ "gray"
    ),
    radius = 7,
    label = ~paste(estacion, "<br>", "Altitud:", altitud_estacion, "m", "<br>", "Estado:", comportamiento),
    fillOpacity = 0.9,
    stroke = FALSE
  ) %>%
  addLegend(
    position = "bottomright",
    colors = c("green", "orange", "red"),
    labels = c("Estable", "Posible problema", "Problema"),
    title = "Comportamiento de Estaciones"
  )