查看表面图{fields}

时间:2014-04-08 19:20:36

标签: r geospatial

我正在尝试在包字段中自定义薄板样条曲线图的外观。例如:

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

如何自定义轴以便获得更多标记和标签? 另外,有没有办法覆盖对应于X轴和Y轴的刻度线的网格? 感谢

1 个答案:

答案 0 :(得分: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")