R语言sf包如何将非度计量地理坐标转换为度单位坐标

你这是典型的CRS元数据标记错误：现在数据里挂的WGS84（EPSG:4326，单位度）是错的。
经纬度的合理取值范围是经度-180_{180、纬度-90}90，你现在的边界框x从-999到9851、y从4660到9851，数值完全超出度的范围，本质是数据本身是投影坐标（单位大概率是米/公里或者做过缩放），但被错误贴了WGS84的标签，mapview按经纬度去解析当然会画错。
注意：这时候绝对不能直接用st_transform做转换，st_transform的前提是现有CRS标记是正确的，标签错了直接转只会得到完全无效的结果。

1. 先清掉错误的CRS标签 #

先把错挂的4326去掉，避免后续操作被误导：

map <- st_set_crs(map, NA)

2. 给数据匹配真实的CRS #

从你给的坐标范围和国家位置对应关系看，最可能的情况是数据原本是Web墨卡托（EPSG:3857，单位米），被做了单位缩放（转成公里）或者线性归一化，先从最常见的情况开始试：

# 先试标准Web墨卡托（单位米）
map_test <- st_set_crs(map, 3857)
plot(st_geometry(map_test), max.plot = 1)

画出来如果大陆轮廓歪、国家位置不对，就试自定义单位为公里的Web墨卡托：

# 自定义单位为公里的Web墨卡托CRS
crs_3857_km <- st_crs("+proj=merc +a=6378137 +b=6378137 +lat_ts=0 +lon_0=0 +x_0=0 +y_0=0 +k=1 +units=km +nadgrids=@null +wktext +no_defs")
map_test <- st_set_crs(map, crs_3857_km)
plot(st_geometry(map_test), max.plot = 1)

如果试完还是不对，就换其他常用全球投影（比如EPSG:3395世界墨卡托、EPSG:4087世界等距圆柱）挨个试，每次试完画个图看国家轮廓和位置对不对——比如中国应该在东亚、美国在北美、欧洲在0度经线附近，位置对了就说明找对了CRS。

不建议盲目猜CRS，最稳妥的方式是找这份数据的原始来源说明，正规发布的空间数据都会明确标注坐标系参数，你大概率是读数据时丢了投影文件，或者之前导出数据的人误改了CRS标签。

3. 转成度为单位的WGS84 #

确认CRS匹配正确后，再做坐标转换到EPSG:4326，得到单位为度的经纬度：

map_deg <- st_transform(map_test, 4326)
# 校验边界框，数值应该落在经度-180~180、纬度-90~90区间
st_bbox(map_deg)

校验通过后再用mapview绘图就不会出错了。

• 导出数据时只改了CRS元数据标签，没有对坐标值做对应的转换，相当于给苹果贴了个梨的标签
• 旧工具导出数据时没有生成.prj投影文件，读取时被自动默认赋值为WGS84
• 坐标做过缩放（比如米转公里、自定义范围归一化），但没有在CRS参数里标注单位换算比例

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