将图例添加到R中的hist2d图中

Question

如何在通过hist2d生成的二维直方图中添加图例？

我正在绘制气候数据（rdata）的两个重叠参数（温度和降水）的密度图。可以访问here。

这是我的代码：

# load libraries
        library(raster)
        library(gplots)
        library(fields)
# load data
        load("clim_data.rdata")
# assign colors
        col <- (colorRampPalette(c("blue", "green", "yellow", "red")))(8)
# plot
        par(mar=c(5,5,2,15), xpd=T)
        hist2d(clim_data$Prec, clim_data$Temp,
               nbins=250,
               col=col,
               FUN=function(x) log(length(x)),
               main="Title",
               xlab="Precipitation [mm]",
               ylab="Temperature [°C]",
               xlim=c(0, 8000), ylim=c(-40, 40), legend=T)

这将进行绘图（没有图例）并导致此错误：

"legend" is not a graphical parameter

我在hist2d RDocumentation中找不到关于图例的任何信息。但还是要尝试一下，因为根据此post的第二个答案，应该可以将图例生成为hist2d图。 presented result正是我想要的。

如何添加这样的颜色条来指示绘图的密度？我知道使用ggplot2或hexbinplot的选项，但希望使用基本图。我也知道colorbar.plot，add.color.bar，color.bar和legend，但是很难找到hist2d生成的密度值。非常感谢！