如何导出无白边且保留NA值透明度的rasterVis绘图？

我来帮你搞定这两个棘手的问题！针对你用rasterVis导出栅格时遇到的白边框无法去除、NA值透明度无法保留的问题，下面是具体的解决方案和完整代码示例：

关键注意点先提一句 #

首先要明确：JPG格式不支持透明效果，如果需要保留NA值的透明度，必须改用PNG格式导出——这是格式本身的限制，没法绕开。

1. 彻底去除白色边框的设置 #

要完全去掉白边框，除了你已经用到的margin=F、scales=list(draw=FALSE)，还需要调整levelplot的布局边距和轴线样式，具体如下：

• 在par.settings里设置layout.heights和layout.widths的所有padding为0，彻底消除图表内部的空白边距
• 将axis.line的颜色设为透明，隐藏默认的轴线
• 必要时用grid包的工具移除残留的边框元素

2. 保留NA值透明度的设置 #

• 改用png()设备初始化导出，指定bg="transparent"确保背景透明
• 在levelplot中添加na.color="transparent"参数，直接将NA对应的区域设为透明

完整可运行代码 #

library(raster)
library(rasterVis)
library(grid)

# 构造带NA值的测试数据（方便验证透明度效果）
m <- matrix(seq(0,100,length.out=15000),150,100)
m[sample(1:15000, 2000)] <- NA # 随机生成一些NA值
r <- raster(m, xmn=0, ymn=0, xmx=ncol(m), ymx=nrow(m))

# 初始化PNG导出设备（注意：JPG不支持透明，必须用PNG）
png(file = "test_transparent.png", 
    bg = "transparent", 
    height = nrow(r), 
    width = ncol(r), 
    res = 96) # res设置为96和屏幕分辨率一致，避免导出后缩放变形

# 绘制栅格图，去除所有边框并保留NA透明度
levelplot(t(r), 
          contour = TRUE, 
          margin = FALSE, 
          scales = list(draw = FALSE), 
          colorkey = NULL,
          par.settings = list(
            axis.line = list(col = "transparent"), # 隐藏轴线
            layout.heights = list(top.padding = 0, bottom.padding = 0,
                                 left.padding = 0, right.padding = 0), # 消除内部边距
            layout.widths = list(left.padding = 0, right.padding = 0)
          ),
          col.regions = colorRampPalette(c("blue", "white", "red")), # 自定义颜色渐变
          na.color = "transparent") # 设置NA值为透明

# 移除可能残留的边框元素（如果有的话）
grid.remove("border", grep = TRUE)

# 关闭设备
dev.off()

运行这段代码后，你会得到一个没有白边框、NA区域完全透明的PNG图片。如果确实需要JPG格式，只需要把png()换成jpeg()，但NA区域会变成白色（因为JPG不支持透明），同时边框仍然可以通过上述方法去除。

内容的提问来源于stack exchange，提问作者freds