Используя qtm из пакета tmap, я могу создать базовый c сюжет боевых мест:
library(tidyverse)
library(tmap)
library(leaflet)
library(leaflet.extras)
library(sp)
GeoBattleData_YZ <- read.csv("GeoBattleData_YZ.csv") # from https://doi.org/10.7910/DVN/1KCCX2
battle_dat<-GeoBattleData_YZ %>%
select(war,atkr,defr,date,long,lat,dist_att,dist_def,duration) %>%
mutate(year=substr(date,1,4))
battle_dat<-battle_dat[!is.na(battle_dat$long),]
battle_dat<-SpatialPointsDataFrame(coords = battle_dat[,c(5:6)],
data = battle_dat,
proj4string = CRS("+init=EPSG:4326"))
tm <-qtm(battle_dat,
symbols.col = 'duration')
tm
map<-tmap_leaflet(tm) %>%
addDrawToolbar(
targetGroup='draw',
polylineOptions=FALSE,
markerOptions = FALSE,
circleOptions = T,
polygonOptions=F,
singleFeature=T,
circleMarkerOptions=F)
map
Ваш Например, я думаю, что в исходной строке tm_shape (dat) следовало вызвать battle_dat, а не dat.