如何移除R语言自定义绘图函数输出图像的右侧边距？

我来帮你搞定这个右侧边距的问题！你当前的代码已经设置了mar和oma为0，但还是有边距，大概率是因为默认的绘图轴边距设置和PNG设备的尺寸匹配问题。下面是修改后的完整函数和关键细节解释：

修改后的函数代码 #

NDVI_Result = function(x, y, output_path) {
  # 读取输入栅格
  r_x = raster::raster(base::normalizePath(x))
  r_y = raster::raster(base::normalizePath(y))
  
  # 计算NDVI
  output = raster::overlay(r_x, r_y, fun = NDVI)
  
  # 可选：裁剪掉栅格边缘的无数据区域（如果有的话）
  # output = raster::trim(output)
  
  # 获取栅格的实际像素尺寸，保证导出图像比例完全匹配
  n_cols = ncol(output)
  n_rows = nrow(output)
  
  # 打开PNG设备，设置与栅格像素一致的尺寸，透明背景
  png(
    filename = output_path,
    width = n_cols,
    height = n_rows,
    bg = NA,
    res = 96  # 分辨率根据需求调整，保持和屏幕一致用96即可
  )
  
  # 配置绘图参数，彻底移除所有额外边距
  graphics::par(
    bg = NA,
    mar = c(0, 0, 0, 0),
    oma = c(0, 0, 0, 0),
    xaxs = "i",  # 强制x轴严格贴合数据范围，无默认的5%边距
    yaxs = "i"   # 强制y轴严格贴合数据范围，无默认的5%边距
  )
  
  # 绘制NDVI栅格，仅保留边框
  raster::plot(
    output,
    axes = FALSE,
    box = TRUE,  # 只绘制栅格的边框
    legend = FALSE,
    frame = FALSE,
    xlim = raster::extent(output)[c(1, 2)],  # 明确限定x范围为栅格的实际范围
    ylim = raster::extent(output)[c(3, 4)],  # 明确限定y范围为栅格的实际范围
    useRaster = TRUE  # 使用光栅化绘制，避免矢量渲染的微小边距误差
  )
  
  # 关闭绘图设备，完成导出
  dev.off()
}

关键修改点说明 #

1. PNG设备尺寸匹配：直接用栅格的列数/行数作为导出图像的宽/高，确保图像和栅格像素1:1对应，不会因为比例不匹配产生额外边距。
2. xaxs="i"和yaxs="i"：这是解决边距的核心！R默认的xaxs="r"会在数据范围外自动添加5%的空白边距，这就是你看到右侧边距的主要原因。设置为"i"会强制绘图区域完全贴合栅格的实际范围。
3. 明确限定xlim和ylim：直接绑定到栅格的extent，进一步确保绘图范围没有多余空间。
4. 可选的trim()操作：如果你的原始栅格边缘存在无数据区域，用trim()裁剪后可以让导出的图像只保留有效数据区域。

如果按上面的修改后还有问题，可以检查下你的NDVI函数是否返回了带有额外无数据边缘的栅格，或者尝试调整res参数（分辨率）来匹配你的需求。

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