我想请教如何在
我需要在leaflet中重现第15行的图表,我在这里找到了一些想法:https://gis.stackexchange.com/questions/239118/r-convert-sf-object-back-to-spatialpolygonsdataframe。然而,使用这种方法并不起作用。
leaflet包
中绘制sf对象
。我知道mapview包
可以绘制它,但是我更喜欢使用leaflet包。下面是我的示例:library(leaflet)
library(eurostat)
library(dplyr)
library(ggplot2)
options(readr.default_locale=readr::locale(tz="Europe/Berlin"))
df60 <- get_eurostat_geospatial(resolution = 60)
CE.sf <- df60 %>%
dplyr::filter(LEVL_CODE == 2 & CNTR_CODE %in% c("AT","CZ","DE","HU","PL","SK")) %>%
dplyr::select(NUTS_ID)
plot(CE.sf)
CE.sf %>%
ggplot() +
geom_sf(color = "black", size = 0.4)
CE = sf::as_Spatial(CE.sf)
leaflet() %>%
addProviderTiles("CartoDB.Positron") %>%
addPolygons(data= CE, color = "green")
我需要在leaflet中重现第15行的图表,我在这里找到了一些想法:https://gis.stackexchange.com/questions/239118/r-convert-sf-object-back-to-spatialpolygonsdataframe。然而,使用这种方法并不起作用。
leaflet
函数之前从sf
对象中删除 CRS。例如,运行:CE.sf = st_set_crs(CE.sf, NA)
然后再运行之前的代码。 - agila