我正在尝试在 mapview 中使用不同的比例主题来帮助可视化收益与损失,包括:
有任何想法吗?
library(tidyverse)
library(mapview)
library(sf)
lat <- rep(34,16)
lon <- seq(-128, -126, length = 16)
value <- c(-1000, -800, -600, -400, -200, -100, -50,
-25, 25, 50, 100, 200, 400, 600, 800, 1000)
#make data.frame
df <- data.frame(lat, lon, value)
#make spatial object for mapview
df <- st_as_sf(df, coords = c("lon", "lat"), crs = 4326) %>%
mutate(value_abs = abs(value)) #value_abs intended for `cex` argument
pal <- mapviewPalette("mapviewSpectralColors") #from mapview doc. example
m <- mapview(df["value"], #sets …Run Code Online (Sandbox Code Playgroud) 有没有办法在连续渐变图例中添加一条线?我无法找到任何使用 ggplot2 的示例。
例如,如何仅向 1.7 处的图例添加红色水平线?
library(ggplot2)
x <- seq(1:1000)
y <-rnorm(1000,0,1)
df <- data.frame(x,y)
ggplot(df, aes(x, y, color = y)) + geom_point()
Run Code Online (Sandbox Code Playgroud)
像这样:
我为什么会喜欢这个可能没有多大意义。为了获得更多背景信息,我希望在这些图的图例中添加一条红线(在 4552 处),以在不同气候变化情景下的不同年度水库流入预测的背景下显示水库的容量,在幻灯片上讨论目的(我不需要注释红线)。有相当多的水库,所以如果可能的话我想用 R 来完成这一切。
谢谢你的任何想法。
这个块在下面运行,并按df_all预期生成,但是当我取消注释顶部的单个函数时(甚至不在这里应用它,但我确实需要其他东西)并重新运行相同的块,我得到:Error in bind_rows_(x, .id): Argument 1 must be a data frame or a named atomic vector, not a function
library(data.table)
# addxtoy_newy_csv <- function(df) {
# zdf1 <- df %>% filter(Variable == "s44")
# setDT(df)
# setDT(zdf1)
# df[zdf1, Value := Value + i.Value, on=.(tstep, variable, Scenario)]
# setDF(df)
#}
tstep <- rep(c("a", "b", "c", "d", "e"), 5)
Variable <- c(rep(c("v"), 5), rep(c("w"), 5), rep(c("x"), 5), rep(c("y"), 5), rep(c("x"), 5))
Value <- c(1,2,3,4,5,10,11,12,13,14,33,22,44,57,5,3,2,1,2,3,34,24,11,11,7)
Scenario <- c(rep(c("i"), 20), …Run Code Online (Sandbox Code Playgroud)