[英]plot ggmap and sf points
I'm struggling to plot a basemap from ggmap with points (sf object).我正在努力 plot 来自 ggmap 的带有点(sf 对象)的底图。 The most promising solution I came across is this SO answer from @andyteucher.
我遇到的最有希望的解决方案是来自@andyteucher 的SO 答案。 I tried to reproduce it below, but I'm not having much luck.
我试图在下面重现它,但我运气不佳。 I get the basemap to print but not the points.
我得到要打印的底图,但不是点。
library(tidyverse)
library(sf)
library(ggmap)
# Define a function to fix the bbox to be in EPSG:3857
# https://stackoverflow.com/a/50844502/841405
ggmap_bbox <- function(map) {
if (!inherits(map, "ggmap")) stop("map must be a ggmap object")
# Extract the bounding box (in lat/lon) from the ggmap to a numeric vector,
# and set the names to what sf::st_bbox expects:
map_bbox <- setNames(unlist(attr(map, "bb")),
c("ymin", "xmin", "ymax", "xmax"))
# Coonvert the bbox to an sf polygon, transform it to 3857,
# and convert back to a bbox (convoluted, but it works)
bbox_3857 <- st_bbox(st_transform(st_as_sfc(st_bbox(map_bbox, crs = 4326)), 3857))
# Overwrite the bbox of the ggmap object with the transformed coordinates
attr(map, "bb")$ll.lat <- bbox_3857["ymin"]
attr(map, "bb")$ll.lon <- bbox_3857["xmin"]
attr(map, "bb")$ur.lat <- bbox_3857["ymax"]
attr(map, "bb")$ur.lon <- bbox_3857["xmax"]
map
}
# requires API key
basemap <- get_map(location=c(lon = 75.85199398072335,
lat = 22.7176905515565),
zoom=9, maptype = 'toner-hybrid',
source = 'google')
basemap_3857 <- ggmap_bbox(basemap)
points <- tribble(
~name, ~lat, ~lon,
"test1", 22.7176905515565, 75.85199398072335,
"test2", 22.71802612842761, 75.84848927237663,
) %>%
st_as_sf(coords = c("lat", "lon"),
crs = 3857)
ggmap(basemap_3857) +
coord_sf(crs = st_crs(3857)) +
geom_sf(data = points,
inherit.aes = FALSE)
I believe you had an issue with your coordinate reference systems - you seem to have used degrees in context of CRS 3857, which is defined in meters (and so several degrees of magnitude off...)我相信您的坐标参考系统存在问题-您似乎在 CRS 3857 的上下文中使用了度数,该度数以米为单位定义(因此相差了几个度数……)
If my hunch is correct you need to first declare your sf object to be in 4326 (WGS84 = the CRS of GPS coordinates) and then - and only then - apply transformation to 3857 (from a known start).如果我的预感是正确的,您需要首先将您的 sf object 声明为 4326(WGS84 = GPS 坐标的 CRS),然后 - 并且只有这样 - 应用转换到 3857(从已知的开始)。
Does this code and map fit your expectations?此代码和 map 是否符合您的期望? (I also cleaned up the get_map call somewhat, as it had mixed Google and Stamen terms; no big deal there)
(我还稍微清理了 get_map 调用,因为它混合了 Google 和 Stamen 术语;没什么大不了的)
# requires API key
basemap <- get_map(location=c(lon = 75.85199398072335,
lat = 22.7176905515565),
zoom=9,
source = 'google')
basemap_3857 <- ggmap_bbox(basemap)
points <- tribble(
~name, ~lat, ~lon,
"test1", 22.7176905515565, 75.85199398072335,
"test2", 22.71802612842761, 75.84848927237663,
) %>%
st_as_sf(coords = c("lon", "lat"),
crs = 4326) %>% # this is important! first declare wgs84
st_transform(3857) # and then transform to web mercator
ggmap(basemap_3857) +
geom_sf(data = points,
color = "red",
inherit.aes = F)
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.