stkde.R

data("storms", package = "dplyr")
## insert data here with approximately similar scales for x, y, t
xyt <- data.frame (x = storms$long, y = storms$lat, 
                   t = unclass(with(storms, ISOdatetime(year, month, day, hour, 0, 0)))/1e6)
library(stkde)  ##  (works on MarinePredators) https://github.com/karpfen/stkde
library(MASS)
stkde_val <- stkde (xyt_in = xyt, hs = mean(c(bandwidth.nrd(xyt$x), bandwidth.nrd(xyt$x)))/3, 
                    ht = bandwidth.nrd(xyt$t)/3, x_size = 30, y_size = 30,
                    t_size = 30)
library(mapview)

## Loading required package: leaflet

## use PageUp/PageDown, arrows up/down, left/right
cubeView(raster::brick(stkde_val))