R语言中ggswissmaps包输出的未知坐标格式识别问题
ggswissmaps包中g1b15图层的坐标系说明
你看到的坐标是瑞士旧平面坐标系(LV03,EPSG代码21781),单位为米,是瑞士传统的测绘专用坐标系,并非通用的UTM或WGS84经纬度格式。
补充说明:
- LV03是瑞士之前的官方平面坐标系,坐标范围完全匹配你输出的数值:x轴(对应你的
long字段)范围约600000-1000000米,y轴(对应lat字段)范围约0-350000米。 - 如果需要转换为Google Maps常用的WGS84经纬度(EPSG:4326),可以用
sf包实现,示例代码如下:
library(sf) library(ggswissmaps) data("shp_df") # 将g1b15转换为sf对象,指定原坐标系为EPSG:21781 g1b15_sf <- st_as_sf(shp_df[["g1b15"]], coords = c("long", "lat"), crs = 21781) # 转换为WGS84经纬度坐标系 g1b15_wgs84 <- st_transform(g1b15_sf, crs = 4326) # 查看转换后的坐标 head(st_coordinates(g1b15_wgs84))
内容的提问来源于stack exchange,提问作者jkortner
相关产品推荐
相关产品推荐

