R Leaflet中SP对象转SF格式实现addPolygons正常绘图的方法

• 方法调度报错根源：geometry()为sf包定义的通用几何提取方法，sp包定义的SP类空间对象无对应方法绑定，因此触发invalid call in method dispatch to 'geometry' (no default method)调度报错
• 视野过度缩放、图层不显示的核心诱因：SP对象与Leaflet交互时存在坐标参考系不匹配、几何结构不合规问题，导致Leaflet无法正确识别空间范围

• 替换报错的属性提取逻辑
  查看SP对象几何属性时，直接调用sp包专属方法，不使用sf体系的通用geometry()函数：
  • 提取空间范围：sp::bbox(目标SP对象)
  • 提取CRS定义：sp::proj4string(目标SP对象)
• 校验并修正坐标参考系
  绝大多数视野过度缩放问题来自SP对象使用投影坐标系（单位为米的平面坐标）而非WGS84经纬度坐标：
  • 坐标值绝对值超过经纬度取值范围（经度绝对值>180、纬度绝对值>90）是坐标系不匹配的核心判断标志
  • 执行格式转换：wgs84_sp_obj <- sp::spTransform(目标SP对象, CRSobj = sp::CRS("+proj=longlat +datum=WGS84 +no_defs"))
• 修正不合规的SP几何结构
  若SP对象的几何字段被手动修改为普通列表结构，Leaflet将无法识别其空间属性，需通过sp::SpatialPolygons()重新构造合规的SP对象后再传入
• 按规则调用addPolygons完成绘制
  无需将SP对象强制转换为sf格式，直接传入合规的WGS84 SP多边形对象即可完成渲染，示例代码：

  leaflet() %>%
    addTiles() %>%
    addPolygons(
      data = wgs84_sp_obj,
      fillColor = "#3182bd",
      color = "#ffffff",
      weight = 1
    )
  

• 手动修正异常视野
  完成坐标转换后若仍出现视野过度缩放，可在图层加载后手动指定地图视野，替代错误的自动边界拟合：

  # 替换为SP对象的实际经纬度范围：最小经度、最小纬度、最大经度、最大纬度
  leaflet() %>%
    addTiles() %>%
    addPolygons(data = wgs84_sp_obj) %>%
    fitBounds(lng1 = 最小经度值, lat1 = 最小纬度值,
              lng2 = 最大经度值, lat2 = 最大纬度值)
  

• 不要直接将SP对象的@polygons插槽内容作为data参数传入addPolygons，必须传入完整的SP空间对象本体
• 加载shapefile时可统一使用sf::st_read()读入，通过sf::as_Spatial()转换为SP对象，可规避rgdal老版本读取时的CRS识别错误

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