R:在filled.contour图上绘制散点图

Question

我对R非常陌生，我用插值数据制作了一个filled.contour图，就像在不规则网格上绘制等高线中的数据一样。使用来自在不规则网格上绘制等高线的一些示例数据，我使用以下代码制作了一个filled.contour和简单的散射图

x <- datr$Lat
y <- datr$Lon
z <- datr$Rain

require(akima)
fld <- interp(x,y,z)
filled.contour(fld)
plot(x,y)

是否有办法在同一地块(覆盖)上制作地块(x，y)和filled.contour(fld)？我试过点(x，y)，但这与x和y轴不匹配。在Matlab中，我相信我会用暂停来完成这个任务，但我不知道如何在R上完成。

谢谢!

Answer 1

为此，可以使用参数plot.title或plot.axes：

x <- 10*1:nrow(volcano)
y <- 10*1:ncol(volcano)
filled.contour(x, y, volcano, plot.title = { 
  points(x = 200, y = 200)
})

(通过)

Answer 2

一种方法是读取filled.contour的代码，并进行如下所示的一些黑客操作：

做好你的身材：

filled.contour(fld)

通过从参数列表复制这些常量来定义它们。

nlevels = 20
zlim = range(z, finite = TRUE)
las = 1 
levels = pretty(zlim, nlevels)
xlim = range(x, finite = TRUE)
ylim = range(y, finite = TRUE)
xaxs = "i"
yaxs = "i"
asp = NA

通过从函数体复制代码来计算这些值。

mar.orig <- (par.orig <- par(c("mar", "las", "mfrow")))$mar
w <- (3 + mar.orig[2L]) * par("csi") * 2.54

通过从函数体复制代码来设置布局

layout(matrix(c(2, 1), ncol = 2L), widths = c(1, lcm(w)))

注意到该图形实际上是在颜色标度之后绘制的，但是我们不需要反转布局的顺序，因为layout实际上将“当前”区域设置为最后一个区域，因为第一次调用plot.new将导致当前区域绕到第一个区域。因此，当您设置绘图窗口并通过以下方法绘制点时：

plot.window(ylim=ylim,xlim=xlim)
points(x,y)
title(main='title',
      sub='Sub-Title',
      xlab='This is the x axis',
      ylab='This is the y axis')

它们按需要叠加图形。