r - 如何在R中绘制双曲向量场？

Question

我正在尝试绘制这样的矢量场：

http://www.aanda.org/articles/aa/full/2006/11/aa4228-05/img58.gif

这将有一个形式为 (y2 /b)-(x2 /k)=1 的方程，

在这样的情节中（http://math.haifa.ac.il/ROVENSKI/rovenski/images/ROV05_628.gif）

有漂亮的小箭头之类的。

有谁知道如何做到这一点？谢谢！

Answer 1

df <- expand.grid(x=seq(-2, 2, .1), y=seq(-2, 2, .1)); k=2; b=4
df$z <- with(df, (y^2)-(x^2) )
library(raster)
library(rasterVis)
r <- rasterFromXYZ(df)
projection(r) <- CRS("+proj=longlat +datum=WGS84")
vectorplot(r, par.settings=RdBuTheme)

Ypu 可以玩弄箭头大小。这个设置在我的屏幕上看起来更令人愉悦：lwd.arrows=0.3,length=unit(.04, 'npc'))