R/传单-绘制多个多边形
我试图用传单包装画出许多多边形,但我不明白出了什么问题 我使用的shapefile可以在这里找到:R/传单-绘制多个多边形,r,leaflet,projection,R,Leaflet,Projection,我试图用传单包装画出许多多边形,但我不明白出了什么问题 我使用的shapefile可以在这里找到: 图书馆(传单) 图书馆(rgdal) df查看df@proj4string和绘图(df)的输出;轴(1);轴(2)。您的形状文件使用特定的CRS。您需要使用一个通用的CRSobj(我从这里获得了CRS代码:)来转换您的空间多边形数据帧 库(sp) pj% 添加多边形(fillOpacity=0.8,color=“Blue”,weight=1) 微管相关蛋白2 它对这些数据有效吗:df谢谢,这很有
图书馆(传单)
图书馆(rgdal)
df查看df@proj4string
和绘图(df)的输出;轴(1);轴(2)
。您的形状文件使用特定的CRS。您需要使用一个通用的CRSobj
(我从这里获得了CRS代码:)来转换您的空间多边形数据帧
库(sp)
pj%
添加多边形(fillOpacity=0.8,color=“Blue”,weight=1)
微管相关蛋白2
它对这些数据有效吗:df谢谢,这很有效。请指定如何在CRS()中找到投影参数?@pierreduk;我知道这是最常见的一个,所以我使用了它(我回答时没有搜索它)。但现在我在R的传单描述中找到了CRS,我将CRS重写到其中(它们的区别只是一个表达)(见:)
library(leaflet)
library(rgdal)
df <- readOGR("C:/Users/me/codes_postaux","codes_postaux_region")
plot(df)
map <- leaflet(df) %>%
addProviderTiles("CartoDB.Positron")%>%
fitBounds(10,38,10,55) %>%
addPolygons(fillOpacity = 0.8, color = "Blue", weight = 1)
map
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