R语言地图叠加自定义多边形时的CRS坐标系统匹配问题

你创建的自定义多边形属于无CRS标识的sf几何对象，只需要为其赋值和瑞典市镇图层一致的CRS即可正常叠加显示，有两种可行的赋值方式：

• 方式一：直接复用瑞典市镇图层的CRS（最稳妥，避免CRS记混出错）
  完成平滑操作后添加如下代码：

  st_crs(smooth_sweden_polygon) <- st_crs(mun)
  

• 方式二：手动指定WGS 84坐标系（对应EPSG编码为4326）
  完成平滑操作后添加如下代码：

  st_crs(smooth_sweden_polygon) <- 4326
  

赋值完成后即可正常运行叠加绘制的ggplot代码，你也可以给自定义多边形添加样式方便查看效果：

ggplot() +
  geom_sf(data=mun,color="grey20",fill="grey95",size=.3) +
  geom_sf(data=smooth_sweden_polygon, fill="steelblue", alpha=0.4, color="navy") +
  theme_bw()

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