lattice.layers

install.packages("lattice.layers")

library(lattice.layers)
# Registered S3 method overwritten by 'lattice.layers':
#   method    from        
#   +.trellis latticeExtra
# 
# 载入程辑包：'lattice.layers'
# The following object is masked from 'package:Ipaper':
# 
#     dev_off
library(Ipaper) # install_githubrpkgs/Ipaper
library(sp)

sp_plot两种绘图语法：

• 如果不提供formula, 采用的是sp::spplot的绘图方法；
• 如果提供formula，采用的是lattice::levelplot的绘图方法。

formula的书写方式与lattice::levelplot的语法格式一致，只是这里的横轴和纵轴坐标是固定lon+lat。以formula = slope ~ lon + lat | type为例：~前的变量slope是绘图变量z值， |之后的变量type表示group。

sp_plot提供了给显著性的格点打patch的功能。需要df中含有mask变量，TRUE表示显著。

spplot style

grid2 = grid_slope
df = grid2@data %>% mutate(mask = pvalue <= 0.05)
df = list(a = df, b = df %>% mutate(slope = slope + 0.1)) %>%
    melt_list("type")

p <- sp_plot(grid2, 
             df, df.mask = df[, .(type, mask)],
             formula = slope ~ lon + lat | type,
             ylim = c(-60, 95), xlim = c(-180, 180) + c(-1, 1)*10,
             layout = c(2, 1),
             colorkey = list(space = "bottom"),
             density = 0.1) + 
    theme_lattice(
        key.margin = c(1, 1, 1, 1),
        plot.margin = c(0, 1, 2, 1))
p

# write_fig(p, "a.pdf", 7, 6)

# add layers
lyrs = layers(
    layer_signDist(density  = 0.1, lwd = 0.2), 
    layer_signPerc(x = 0.01, y = 0.72, fill = "white"),
    layer_barchart(y = 0.12))
p + lyrs