我有一个包含经度、纬度和一个显示湿度的值列的数据集,长度约为 300 行。每个点显示不同位置的湿度。我想将它们全部绘制在地图上,并根据它们的值(例如渐变颜色)为它们着色并添加图例。它有点类似于这里的问题,但我无法让它工作。代码基本上就在那里,但只有在图例中正确地着色和显示它并不能真正起作用。这些点代表非洲的一条线,湿度值最初是从栅格数据集中提取的,它们包含几个数字。我创建了一些示例数据来说明我卡在哪里。
library("maps")
library("raster")
# create sample data
lon <- seq(from=35.6, to=43.2, by=0.2)
lat <- seq(from=10.5, to=22.2, by=0.2)
humidity <- runif(59, min=9.6, max=13.5)
data <- data.frame(lon,lat, humidity)
colfunc<-colorRampPalette(c("dodgerblue2","khaki","orangered")) # create colours
map('world', xlim = c(20, 80), ylim = c(5, 30), lwd=0.5, col = "grey95", fill = T, interior = FALSE)
title("specific humidity along line")
map.axes()
points(data$lon, data$lat, cex=.5, pch=19, col=colfunc(100))
legend("topleft",title="q (g/kg)",legend=c(11,11.5,12,12.5,13),col =colfunc(100), pch=20)
Run Code Online (Sandbox Code Playgroud)
结果图如下所示:
图例显然有问题,我想在图例中显示一些具有相应颜色和值的点,甚至使用漂亮的颜色条。我不确定为什么图例中的颜色只是蓝色。我还怀疑点线没有根据它们的实际值着色,只是显示整个颜色渐变。感谢您的任何建议!
更新来自 Alex 的代码:
n <- 10
colfunc<-colorRampPalette(c("dodgerblue2","khaki","orangered")) # create …Run Code Online (Sandbox Code Playgroud)