曲面图外观{fields}

Question

我正尝试在封装字段中自定义薄板样条曲线的表面外观。例如：

require(fields)
fit<- Tps(ozone$x, ozone$y)
set.panel()
surface(fit, type="C", xlab=("East"), ylab=("North"))

如何自定义轴，以便获得更多刻度线和标签？另外，有没有办法覆盖X轴和Y轴上与刻度线对应的网格？谢谢

Answer 1

使用axes=F绘制不带轴的绘图，然后使用axis手动添加它们。

例如

surface(fit, type="C", xlab=("East"), ylab=("North"), axes=F)
# X axis every 5
axis(1, at=seq(-20, 30, 5))
# Y axis, at various points with custom labels. las=1 puts labels horizontally
axis(2, at=c(-20, -10, -5, 0, 20), labels=letters[1:5], las=1)

对于栅格，对于规则栅格使用grid命令，对于任意间距的栅格使用abline。

# 5x5 grid
grid(nx=5, ny=5, col="black", lty="dotted")
# Arbitrary gridlines
abline(v=seq(-20, 10, 5), h=seq(-10, 10, 2), col="gray")