library(tidyverse)
library(readxl)
library(mapview)
library(sf)# Shared Table from: https://docs.google.com/spreadsheets/d/1rMqHBqDd8S8Gi969yqZK6NH8-kCXZnJS/edit?usp=sharing&ouid=113509074802217952757&rtpof=true&sd=true
pasivas <- read_excel("D:/BoxFiles/Box Sync/CodigoR/Bats/data/grabaciones de Murcielagos.xlsx",
sheet = "Grabaciones Pasivas")
pasivas$Lat <-as.numeric(pasivas$Lat)
pasivas$Lon <-as.numeric(pasivas$Lon)
referencia <- read_excel("D:/BoxFiles/Box Sync/CodigoR/Bats/data/grabaciones de Murcielagos.xlsx",
sheet = "Grabaciones Referencia")
pasivas_sf <- st_as_sf(pasivas, coords = c("Lon", "Lat"), crs = "+proj=longlat +datum=WGS84 +no_defs")
referencia_sf <- st_as_sf(referencia, coords = c("Lon", "Lat"), crs = "+proj=longlat +datum=WGS84 +no_defs")El tamaño del círculo representa el número de noches.
pasivas_sf %>%
# filter (Tipo == "Grabación Pasiva") %>%
select(Nombre, Ecosistema, Equipo_marca, Noches) %>%
mapview(zcol = "Equipo_marca", cex = "Noches", #at = seq(1, 5, 10, 20, 300),
alpha = 0.1,
map.types = c("CartoDB.DarkMatter","Esri.WorldImagery", "OpenStreetMap.DE")
)El color representa la familia. La localización no es exacta, es aproximada con el centroide del municipio.
referencia_sf %>%
# filter (Tipo == "Grabación Pasiva") %>%
select(Nombre, Ecosistema, Familia, Especimen_museo, Equipo_marca_ref) %>%
st_jitter(factor = 0.001) %>%
mapview(zcol = "Familia",# cex = "Equipo_marca_ref", #at = seq(1, 5, 10, 20, 300),
alpha = 0.2,
map.types = c("CartoDB.DarkMatter","Esri.WorldImagery", "OpenStreetMap.DE")
)