R /传单 - 绘制多个多边形
R / leaflet - drawing numerous polygons
我正在尝试使用传单包绘制大量多边形,但我不明白出了什么问题。
我使用的 shapefile 可以在这里找到:https://www.data.gouv.fr/en/datasets/fond-de-carte-des-codes-postaux/
library(leaflet)
library(rgdal)
df <- readOGR("C:/Users/me/codes_postaux","codes_postaux_region")
plot(df)
shapefile 对我来说似乎没问题,我使用的代码也很简单。但是我只得到地图作为输出而没有多边形。我已经为这个问题苦苦挣扎了很长时间,如果有人能在这里帮助我,我将不胜感激。
map <- leaflet(df) %>%
addProviderTiles("CartoDB.Positron")%>%
fitBounds(10,38,10,55) %>%
addPolygons(fillOpacity = 0.8, color = "Blue", weight = 1)
map
查看 df@proj4string
和 plot(df); axis(1); axis(2)
的输出。您的 shapefile 使用特定的 CRS。你需要用一个普通的 CRSobj
来转换你的 SpatialPolygonsDataFrame
(我从这里得到了 CRS 代码:Leaflet for R: Raster Images)。
library(sp)
pj <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")
df2 <- spTransform(df, pj)
map2 <- leaflet(df2) %>%
addProviderTiles("CartoDB.Positron")%>%
fitBounds(10,38,10,55) %>%
addPolygons(fillOpacity = 0.8, color = "Blue", weight = 1)
map2
我正在尝试使用传单包绘制大量多边形,但我不明白出了什么问题。
我使用的 shapefile 可以在这里找到:https://www.data.gouv.fr/en/datasets/fond-de-carte-des-codes-postaux/
library(leaflet)
library(rgdal)
df <- readOGR("C:/Users/me/codes_postaux","codes_postaux_region")
plot(df)
shapefile 对我来说似乎没问题,我使用的代码也很简单。但是我只得到地图作为输出而没有多边形。我已经为这个问题苦苦挣扎了很长时间,如果有人能在这里帮助我,我将不胜感激。
map <- leaflet(df) %>%
addProviderTiles("CartoDB.Positron")%>%
fitBounds(10,38,10,55) %>%
addPolygons(fillOpacity = 0.8, color = "Blue", weight = 1)
map
查看 df@proj4string
和 plot(df); axis(1); axis(2)
的输出。您的 shapefile 使用特定的 CRS。你需要用一个普通的 CRSobj
来转换你的 SpatialPolygonsDataFrame
(我从这里得到了 CRS 代码:Leaflet for R: Raster Images)。
library(sp)
pj <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")
df2 <- spTransform(df, pj)
map2 <- leaflet(df2) %>%
addProviderTiles("CartoDB.Positron")%>%
fitBounds(10,38,10,55) %>%
addPolygons(fillOpacity = 0.8, color = "Blue", weight = 1)
map2