R语言批量生成地图投影失败问题求助

1. CRS未找到错误 #

• 原因：proj.org上的部分投影仅提供方法名（如merc），但sf需要完整的PROJ参数字符串或标准EPSG代码才能识别，单独的方法名无法被解析为有效坐标系。
• 解决：
  • 不要直接用投影方法名，构造完整的PROJ4字符串。比如墨卡托投影的完整参数是："+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs"
  • 批量处理前先验证CRS有效性，无效的直接跳过：

    check_crs <- function(proj_str) {
      !is.na(sf::st_crs(proj_str))
    }
    valid_projs <- Filter(check_crs, your_proj_list)
    

2. LinearRing点数量无效错误 #

• 原因：部分投影会让全球多边形（如国家边界）出现拓扑错误，比如自相交、顶点顺序颠倒，或者跨日界线的几何体投影后断裂，导致LinearRing不符合“至少4个点且首尾重合”的要求。
• 解决：
  • 投影前修复几何体拓扑：

    world_data <- sf::st_make_valid(world_data)
    

  • 处理跨日界线的几何体，避免投影后异常：

    world_data <- sf::st_wrap_dateline(world_data, options = c("WRAPDATELINE=YES", "DATELINEOFFSET=180"))
    

  • 筛选有效几何体，跳过无法修复的部分：

    world_valid <- world_data[sf::st_is_valid(world_data), ]
    

3. 地图出现异常线条 #

• 原因：很多投影（如方位、圆锥类）有特定的有效范围，全球数据强行投影会导致超出范围的部分被拉伸折叠；另外未修复的拓扑错误也会产生乱线。
• 解决：
  • 绘图时开启裁剪，去掉超出可视区域的内容：

    ggplot(world_valid) +
      geom_sf() +
      coord_sf(crs = proj_str, clip = "on")
    

  • 给需要特定参数的投影补全参数，比如阿尔伯斯投影必须指定标准纬线：

    albers_proj <- "+proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=37.5 +lon_0=-96 +datum=NAD83 +units=m +no_defs"
    

  • 跳过局部投影（如UTM分带），这类投影本来就不适合绘制全球地图。

4. 批量处理优化建议 #

• 用安全函数包裹绘图逻辑，避免单个报错中断整个批量任务：

  safe_map_plot <- purrr::safely(function(proj_str, data, save_path) {
    plot <- ggplot(data) +
      geom_sf() +
      coord_sf(crs = proj_str, clip = "on") +
      labs(title = proj_str)
    ggsave(save_path, plot, width = 10, height = 8)
    return(plot)
  })
  
  # 遍历有效投影生成地图
  map_results <- purrr::map(valid_projs, function(p) {
    filename <- paste0("map_", gsub("[^a-z0-9]", "_", tolower(p)), ".png")
    safe_map_plot(p, world_valid, filename)
  })
  
  # 查看报错的投影详情
  error_details <- purrr::keep(map_results, ~!is.null(.x$error))
  

内容的提问来源于stack exchange，提问作者Jean-Paul Bardou