如何在水平图（晶格）中使用非常规数字制作常规标签？

Question

我需要制作带有包装格子的气象图。但是我在comand levelplot（）中发现了一个问题。我可以使用常规标签使用comand制作地图。例如：标签为0、0.1、0.15、0.2、0.25、0.3、0.35 0.4 ...的相关图...在此示例中，标签跳过了五乘五（图1）。

h1<- levelplot(var~x*y,data = idw.msk.dfr,contour=F,at=seq(0,0.5,0.05),
           par.settings = paleta1,main = "correlation map",
           xlab = NULL, ylab = NULL, ylim = c(-60,15), xlim = c(-90,-30))

图1：

但是，我需要使用非常规值制作地图。例如：0、0.1、0.15、0.2、0.22、0.25、0.40 ...将这些值放入代码中后，我得到此结果（图2）：

h1<- levelplot(var~x*y,data = idw.msk.dfr,contour=F,at=c(0,0.1,0.15,0.2,0.22,0.25,0.4,0.5),
           par.settings = paleta1,main = "correlation map",
           xlab = NULL, ylab = NULL, ylim = c(-60,15), xlim = c(-90,-30))

图2

请注意，地图的标签非常奇怪且不规则。

所以。我该如何解决这个问题？我会感谢您的帮助。

Answer 1

您需要指定自定义色键。将色键添加到levelplot函数中。

x <- seq(pi/4, 5 * pi, length = 100)
y <- seq(pi/4, 5 * pi, length = 100)
r <- as.vector(sqrt(outer(x^2, y^2, "+")))
grid <- expand.grid(x=x, y=y)
grid$z <- cos(r^2) * exp(-r/(pi^3))

breaks <- c(0, 0.1, 0.15, 0.2, 0.22, 0.25, 0.4, 0.5)

levelplot(z~x*y, grid, at=breaks)

myColorkey=list(at=breaks, labels=list(at=breaks, labels=breaks))

levelplot(z~x*y, grid, colorkey=myColorkey)

要使间隔具有相同的大小，请更改at参数：

ats=seq(0, 0.5, by=0.07)
myColorkey=list(at=ats, labels=list(at=ats, labels=breaks))
levelplot(z~x*y, grid, colorkey=myColorkey)