如何用rerddap的griddap处理无PM180版本的0-360°经度数据集

当遇到只有0°-360°经度格式的ERDDAP数据集，又不想手动转换经度范围的话，可以用下面这个封装函数自动处理——核心是利用rerddap::info()拿到数据集的经度元信息，自动适配请求范围，最后再把结果转回到-180°-180°格式。

通用实现代码 #

library(rerddap)
library(dplyr)

griddap_lon_convert <- function(dataset_id, ..., lon = c(-180, 180), convert_back = TRUE) {
  # 拉取数据集的元信息，获取经度范围
  ds_info <- info(dataset_id)
  lon_range <- ds_info$alldata$longitude$range
  
  # 判断数据集是不是0-360的经度格式
  is_360 <- lon_range[2] > 180
  
  if (is_360) {
    user_lon <- lon
    # 处理跨0度经线的请求（比如用户要170°到-170°的范围）
    if (user_lon[1] > user_lon[2]) {
      # 拆成两段范围分别请求
      lon_part1 <- c(user_lon[1], 360)
      lon_part2 <- c(0, user_lon[2] + 360)
      
      data1 <- griddap(dataset_id, ..., lon = lon_part1)
      data2 <- griddap(dataset_id, ..., lon = lon_part2)
      
      # 合并两段数据
      combined_data <- bind_rows(data1$data, data2$data)
    } else {
      # 常规范围：把负数经度转成0-360格式
      converted_lon <- user_lon + ifelse(user_lon < 0, 360, 0)
      combined_data <- griddap(dataset_id, ..., lon = converted_lon)$data
    }
    
    # 把结果的经度转回-180到180格式（默认开启）
    if (convert_back) {
      combined_data <- combined_data %>%
        mutate(longitude = ifelse(longitude > 180, longitude - 360, longitude))
    }
    
    # 保留原griddap返回对象的结构和属性
    result <- list(data = combined_data, summary = ds_info, url = attr(data1, "url"))
    class(result) <- class(data1)
    return(result)
  } else {
    # 数据集本身就是-180到180，直接调用原函数
    return(griddap(dataset_id, ..., lon = lon))
  }
}

使用示例 #

比如要访问某个只有0-360版本的数据集，直接传你习惯的-180°-180°范围就行：

# 替换成你要访问的数据集ID
my_data <- griddap_lon_convert(
  dataset_id = "your_target_dataset",
  time = c("2020-01-01", "2020-01-02"),
  latitude = c(20, 40),
  lon = c(-120, -80)  # 不用手动转成240-280，函数自动处理
)

# 检查处理后的经度范围，应该是-120到-80
range(my_data$data$longitude)

关键逻辑说明 #

• 自动识别数据集经度格式：通过info()返回的经度最大值判断是不是0-360格式
• 兼容两种请求场景：不管是常规西经范围，还是跨0°经线的范围（比如太平洋跨180°的区域），都能自动拆分/转换请求
• 保留原API的使用体验：返回的对象和griddap()原生返回的结构一致，后续处理不用改代码

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