R：绘制垂直于等高线的箭头

Question

我有两个向量表示散点图中的 x 和 y 坐标，每个 (x,y) 坐标都有一个第三变量 (z)，表示从中绘制等高线的变量。示例数据如下：

df<-data.frame(x=runif(n=30,min=-6,max=6),
               y=runif(n=30,min=-6,max=10),
               z=seq(1,100,length.out=30))

我使用 R 包 akima 为等高线图生成 z 矩阵

library(akima)
M1 <- interp(x=df$x,y=df$y,z=df$z)
contour(x=M1$x,y=M1$y,z=M1$z)

我现在想绘制垂直于轮廓线的箭头，最好使用 R 包 pracma 中的函数“quiver”之类的东西，每个 (x,y) 处都有一个箭头的原点- 坐标，箭头指向轮廓线的渐变方向。有没有办法做到这一点？

到目前为止，我最好的想法是以某种方式提取轮廓线的 (x,y)-梯度，并将其用作 quiver 函数中的速度。

感谢您的帮助。

Answer 1

pracma 包有一个gradient 函数，可以使用原始的M1$z 值为您执行此操作。例如，使用你的代码在set.seed(123)之后得到M1：

contour(x=M1$x,y=M1$y,z=M1$z, asp = 1) # asp = 1 needed so things look perpendicular

library(pracma)
g <- gradient(M1$z, M1$x, M1$y)

x <- outer(M1$x, M1$y, function(x, y) x)
y <- outer(M1$x, M1$y, function(x, y) y)

quiver(x, y, g$Y, g$X, scale = 0.02, col = "blue")

请注意，quiver 图中的渐变标签已被交换。也许我设置了 x 和 y 值，这些值是从包的预期方式转换而来的。这是你得到的：