R: map 上的颜色点基于它们的值(value)并添加图例

Question

我有一个包含经度、纬度和显示湿度的值列的数据集，长度大约为 300 行。每个点显示不同位置的湿度。我想将它们全部绘制在 map 上并根据它们的值(例如渐变颜色)为它们着色并添加图例。它有点类似于 question here ，但我无法让它工作。代码基本上就在那里，但只有着色和在图例中正确显示它并没有真正起作用。这些点代表非洲的一条线，湿度值最初是从栅格数据集中提取的，它们包含几位数字。我创建了一些样本数据来说明我被困在哪里。

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)

结果图如下所示:图例显然有问题，我想在图例中显示一些具有相应颜色和值的点，甚至使用漂亮的颜色条。我不确定为什么图例中的颜色只是蓝色。我还怀疑点线没有根据它们的实际值着色，只是显示整个颜色渐变。感谢您的任何建议!更新来自 Alex 的代码:

n <- 10

colfunc<-colorRampPalette(c("dodgerblue2","khaki","orangered")) # create colours
mycol <- function(x, myrange, n=10) round( 1+(x-myrange[1])/diff(myrange) * (n-1))
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(n)[mycol(humidity, range(humidity), n)])
mylist <- c(10,11,11.5,12,12.5,13)
legend("topleft",title="q (g/kg)",legend=mylist,col = colfunc(n)[mycol(mylist,range(humidity), n)], pch=20)

生成此图:点是重叠的，很难看到点的整体值，有没有办法根据定义的范围使用色带给点上色？例如“红色”代表值 10 到 11，“绿色”代表 11 到 12 等等？

Answer 1

您可能混淆了一些东西。在您的代码中，您绘制的点的颜色只是点顺序的结果(第一个点在列表中获得第一个颜色等)。颜色不取决于值。

现在，在整个范围 0:100 的湿度值的颜色渐变中，坦率地说，您将看不到值 11 和 13 之间的任何差异。您需要更多的对比度。所以你应该先做

mycol <- function(x, myrange, n=100) round( 1+(x-myrange[1])/diff(myrange) * (n-1))

现在 mycol(x, range(humidity), n) 将返回一个整数，最小值是 1，最大值是 n。

n=100
points(data$lon, data$lat, cex=.5, pch=19, col=colfunc(n)[mycol(humidity, range(humidity), n)])

mylist <- c(11,11.5,12,12.5,13)
legend("topleft",title="q (g/kg)",legend=mylist,col = colfunc(n)[mycol(mylist,range(humidity), n)], pch=20)