如何在R语言中批量转换Sentinel 5p的methane_mixing_ratio_bias_corrected为GeoTIFF？

1. 安装并加载所需R包 #

处理Sentinel-5P的NetCDF格式数据，需要用到读取NetCDF的ncdf4包，以及栅格数据处理的terra包（替代旧版raster包，性能更优）：

install.packages(c("ncdf4", "terra"))
library(ncdf4)
library(terra)

2. 准备数据 #

将获取到的两个Sentinel-5P NetCDF文件放到同一个本地文件夹中，记录该文件夹的绝对路径。

3. 批量转换脚本 #

以下脚本会自动遍历目标文件夹中的所有NetCDF文件，提取PRODUCT/methane_mixing_ratio_bias_corrected变量，转换为WGS84坐标系的GeoTIFF并保存：

# 替换为你的数据文件夹路径
data_dir <- "/path/to/your/nc/files"

# 获取文件夹内所有NetCDF文件（匹配.nc后缀）
nc_files <- list.files(data_dir, pattern = "\\.nc$", full.names = TRUE)

# 循环处理每个文件
for (file in nc_files) {
  # 打开NetCDF文件连接
  nc_conn <- nc_open(file)
  
  # 提取目标变量（注意变量的完整路径）
  methane <- ncvar_get(nc_conn, varid = "PRODUCT/methane_mixing_ratio_bias_corrected")
  # 提取经纬度数据
  lon <- ncvar_get(nc_conn, varid = "PRODUCT/longitude")
  lat <- ncvar_get(nc_conn, varid = "PRODUCT/latitude")
  
  # 创建栅格对象，设置坐标系为WGS84
  raster_obj <- rast(methane,
                     xmin = min(lon), xmax = max(lon),
                     ymin = min(lat), ymax = max(lat),
                     crs = "EPSG:4326")
  
  # 翻转栅格（Sentinel-5P原始数据的纬度是倒序的）
  raster_obj <- flip(raster_obj, "vertical")
  
  # 生成输出文件名（替换.nc为.tif）
  output_path <- gsub("\\.nc$", ".tif", file)
  
  # 保存为GeoTIFF
  writeRaster(raster_obj, output_path, format = "GTiff", overwrite = TRUE)
  
  # 关闭NetCDF连接
  nc_close(nc_conn)
  
  cat("转换完成：", basename(file), "→", basename(output_path), "\n")
}

关键注意事项 #

• 若变量路径存在差异（比如部分文件的变量未嵌套在PRODUCT组下），需调整varid参数的路径。
• 如果经纬度是二维数组，min(lon)报错的话，可以改用range(lon, na.rm = TRUE)来获取范围值。
• 大文件处理时可添加memfrac参数控制内存占用，比如writeRaster(..., memfrac = 0.5)。

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