使用 sf 列调整 Leaflet (R) 中的边界
Adjusting bounds in Leaflet (R) with sf column
我正在用 Shiny 制作一个小应用程序,它保存国家和地区的数据,用户可以在其中选择一个地区。然后我的想法是我在应用程序中拥有的传单地图将放大并专注于所选区域(即用户单击 "Europe" 并且地图放大到欧洲)。
我不知道应该如何使用简单的特征geometry
列作为传单地图的过滤器。这是一个简单的例子(不是在 Shiny 中,但我想问题与 Shiny 无关)。
library(rnaturalearth)
library(dplyr)
library(leaflet)
# sf data:
earth <- ne_countries(returnclass = "sf") %>%
select(Region = region_un, geometry)
# little dataset:
df <- data_frame(
Region = c("Europe", "Africa", "Asia", "Oceania", "Americas"),
Score = c(0.85, 0.77, 0.81, 0.93, 0.79)
)
# join:
df <- full_join(df, earth)
# simulate what I'm doing in Shiny:
input <- list()
input$region <- "Europe"
df2 <- filter(df, Region == input$region)
leaflet(df2) %>% addTiles()
这会产生:
这与我使用 df
(未过滤的数据帧)相同。关于如何解决这个问题的任何想法?我在 Shiny/leaflet 文档中找不到它。
我们可以使用 sf::st_coordinates
从 df2$geometry
中提取值,获取平均纬度和经度,然后使用 leaflet::setView()
关注我们感兴趣的区域:
library(sf)
coords <- st_coordinates(df2$geometry)
lng <- mean(coords[,1])
lat <- mean(coords[,2])
leaflet(df2) %>% addTiles() %>%
setView(lng, lat, zoom = 3) # 3 for "continent" level
有一些特殊的方法可以在 Shiny 设置中更新 leaflet
地图,主要是 leafletProxy()
,这些在 docs.
中有很好的描述
我正在用 Shiny 制作一个小应用程序,它保存国家和地区的数据,用户可以在其中选择一个地区。然后我的想法是我在应用程序中拥有的传单地图将放大并专注于所选区域(即用户单击 "Europe" 并且地图放大到欧洲)。
我不知道应该如何使用简单的特征geometry
列作为传单地图的过滤器。这是一个简单的例子(不是在 Shiny 中,但我想问题与 Shiny 无关)。
library(rnaturalearth)
library(dplyr)
library(leaflet)
# sf data:
earth <- ne_countries(returnclass = "sf") %>%
select(Region = region_un, geometry)
# little dataset:
df <- data_frame(
Region = c("Europe", "Africa", "Asia", "Oceania", "Americas"),
Score = c(0.85, 0.77, 0.81, 0.93, 0.79)
)
# join:
df <- full_join(df, earth)
# simulate what I'm doing in Shiny:
input <- list()
input$region <- "Europe"
df2 <- filter(df, Region == input$region)
leaflet(df2) %>% addTiles()
这会产生:
这与我使用 df
(未过滤的数据帧)相同。关于如何解决这个问题的任何想法?我在 Shiny/leaflet 文档中找不到它。
我们可以使用 sf::st_coordinates
从 df2$geometry
中提取值,获取平均纬度和经度,然后使用 leaflet::setView()
关注我们感兴趣的区域:
library(sf)
coords <- st_coordinates(df2$geometry)
lng <- mean(coords[,1])
lat <- mean(coords[,2])
leaflet(df2) %>% addTiles() %>%
setView(lng, lat, zoom = 3) # 3 for "continent" level
有一些特殊的方法可以在 Shiny 设置中更新 leaflet
地图,主要是 leafletProxy()
,这些在 docs.