求助：如何在R中基于UTM坐标绘制带点位的地图？

解决UTM坐标转经纬度+地图绘制的完整方案 #

1. 核心问题排查

你当前代码的两个关键错误：

• 没有将UTM坐标转换为WGS84经纬度（leaflet/大多数地图服务要求的坐标系），直接把UTM的x/y当作经纬度传入addMarkers，导致点位完全错位
• 若要加载Google地图，需要配置API密钥（当前Google Maps API需付费且需密钥），OpenStreetMap可直接通过leaflet的provider tiles调用

下面用sf包（替代老旧的sp包，更简洁）配合leaflet、ggplot2实现需求：

2. 准备工作：加载包与数据

# 安装所需包（首次运行）
install.packages(c("sf", "leaflet", "ggplot2", "ggmap"))

# 加载包
library(sf)
library(leaflet)
library(ggplot2)
library(ggmap)

# 你的示例数据
coordenadas <- data.frame(
  x = c(474731, 518127, 488067, 459943, 497732, 502623),
  y = c(4539196, 4539351, 4521530, 4528829, 4539816, 4531443),
  huso = rep(30, 6)
)

3. UTM转经纬度（关键步骤）

用sf包创建空间对象并转换坐标系：

# 创建UTM空间对象（EPSG:32630对应UTM 30N WGS84）
coords_sf <- st_as_sf(coordenadas, 
                      coords = c("x", "y"),
                      crs = st_crs("+proj=utm +zone=30 +north +ellps=WGS84 +units=m +no_defs"))

# 转换为WGS84经纬度（EPSG:4326，地图服务通用）
coords_wgs84 <- st_transform(coords_sf, crs = 4326)

# 提取经纬度到数据框（方便调用）
coords_wgs84_df <- coords_wgs84 %>% 
  mutate(lng = st_coordinates(.)[,1],
         lat = st_coordinates(.)[,2])

4. 用leaflet绘制地图（解决地图加载问题）

# 初始化leaflet，用OpenStreetMap底图
m <- leaflet(coords_wgs84_df) %>%
  setView(lng = mean(coords_wgs84_df$lng), 
          lat = mean(coords_wgs84_df$lat), 
          zoom = 10) %>%  # 自动聚焦到点位中心
  addProviderTiles(providers$OpenStreetMap)  # 加载OpenStreetMap底图

# 添加点位标记
m %>% addMarkers(lng = ~lng, lat = ~lat)

• 若要加载Google地图：需先获取Google Maps API密钥，然后使用addProviderTiles(providers$Google.Maps, options = providerTileOptions(apiKey = "你的密钥"))，注意Google API当前需付费

5. 用ggplot2+ggmap绘制地图

# 获取点位所在区域的底图（用Stamen免费底图替代Google）
map_base <- get_map(location = c(lon = mean(coords_wgs84_df$lng), 
                                 lat = mean(coords_wgs84_df$lat)),
                    zoom = 10,
                    maptype = "terrain",
                    source = "stamen")

# 绘制地图+点位
ggmap(map_base) +
  geom_point(data = coords_wgs84_df, 
             aes(x = lng, y = lat), 
             color = "red", size = 3) +
  labs(title = "UTM点位地图")

关键提示

• 若你的数据集包含多个UTM带号（huso不同），需按带号分组转换坐标系，再合并为WGS84的数据集
• sf包是当前R空间数据处理的标准工具，比sp包更直观且维护更积极

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