如何使geom_hex在地图上覆盖 我试图在地图上覆盖跨六角形的lunge_count平均值；但是，我只能产生两个单独的图（一个GEOM_HEX和一个地图之一）。 #Libraries 图书馆（dplyr）＃for

Question

      ptt depth_m duration_s duration_m lunge_count year               id       lon      lat
9946 5883      51        778   12.96667           0 2023 2023MX-Mno-05883 -107.2284 22.81825
9947 5883      53        714   11.90000           0 2023 2023MX-Mno-05883 -107.2284 22.81825
9948 5883      49        696   11.60000           1 2023 2023MX-Mno-05883 -107.2284 22.81825
9949 5883      45        756   12.60000           0 2023 2023MX-Mno-05883 -107.2284 22.81825
9950 5883      45        780   13.00000           0 2023 2023MX-Mno-05883 -107.2284 22.81825
9951 5883      45        752   12.53333           0 2023 2023MX-Mno-05883 -107.2284 22.81825

我正在使用的代码：＃海洋地图信息

esri_ocean <- paste0('https://services.arcgisonline.com/arcgis/rest/services/','Ocean/World_Ocean_Base/MapServer/tile/${z}/${y}/${x}.jpeg') ＃创建海洋基地

basemap36 <- ggplot() + annotation_map_tile(type = esri_ocean, zoomin = 1, progress = "none") + coord_sf(crs = 4326,xlim = c(-100, -125),ylim = c(17, 36))

＃绘图地图和六边形弓步计数平均

basemap36 & ggplot(breeding_redo_5801, aes(lon, lat, z = lunge_count)) & stat_summary_hex(fun = mean) ＃结果

（

Https：//i.sstatic.net/old8ben1.png

）

＃其他尝试

i我能够将图与下面的代码覆盖地图覆盖地图；但是，它产生的总和而不是六角形的平均值。

basemap36 & geom_hex(data = breeding_redo_5801, aes(x = lon, y = lat, z = lunge_count)) & stat_summary_hex(fun = mean)

＃这次尝试的结果

输入图像描述在这里

我非常感谢您的帮助！

我相信这会起作用，但是我认为我们没有足够的示例数据来尝试它：

basemap36 & stat_summary_hex(
  data = breeding_redo_5801, 
  mapping = aes(lon, lat, z = lunge_count), 
  fun = mean
)