我有一个lat/lon数据集和一个时间戳.我希望标记的颜色用连续的调色板显示时间.我正在使用colorNumeric()之前创建的朱利安日期julian(x, "2015-01-01").
data = structure(list(timestamp = structure(c(1434056453, 1434148216, 1434153635, 1434245436, 1434358840,
1434364288, 1434369611, 1434461435, 1434466830, 1434558725), class = c("POSIXct", "POSIXt"), tzone = ""),
lon = c(-119.8777, -119.9614, -119.8769, -119.8775, -120.2283,
-120.2285, -119.8429, -120.0954, -120.3957, -120.4421),
lat = c(34.4041,34.376, 34.4061, 34.4021, 34.4696,
34.4697, 34.1909, 34.4328, 34.4554, 34.4456),
ID = as.factor(c("Z11","Z05","Z01", "Z04", "Z11", "Z04","Z01","Z05","Z05","Z11"))),
.Names = c("timestamp", "lon", "lat", "ID"),
row.names = c(1:10),
class = "data.frame")
data$julian = as.numeric(julian(data$timestamp, origin = "2015-01-01"))
pal = colorNumeric( palette = rainbow(7), …Run Code Online (Sandbox Code Playgroud) R的传单现在可以选择用于addMapPane按确定的顺序堆叠传单对象。这适用于矢量数据,但使用它会addRasterImage返回错误消息。options使用pane带有栅格的a 是否有特殊之处?
library(leaflet)
library(raster)
r <- raster(ncol=1000, nrow=1000)
r[] <- runif(ncell(r),0,1)
extent(r) <- matrix(c(172, -37, 175, -38), nrow=2)
crs(r)<-sp::CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")
pal <- colorNumeric(c("#0C2C84", "#41B6C4", "#FFFFCC"), values(r),
na.color = "transparent")
m <- leaflet(m) %>% addTiles() %>%
addMapPane("baseMap", zIndex = 410) %>%
addMapPane("baseSat", zIndex = 420) %>%
addMapPane("data", zIndex = 425) %>%
addMapPane("r", zIndex = 430) %>%
addProviderTiles(providers$Esri.WorldImagery, group = "Satellite imagery",
options = pathOptions(pane = "baseSat")) %>%
addProviderTiles(providers$Esri.WorldTerrain,
options = pathOptions(pane …Run Code Online (Sandbox Code Playgroud)