我有一个空间多边形数据框,其坐标已转换为胆立体类型。但是我现在无法映射该 spdf leaflet()
。我不确定这是否是我提供的分辨率的问题(因为我不知道这些是什么或做什么),或者它是否与我提供的平铺地图有关addTiles()
,或者它是否完全是其他东西。
因此,在将数据转换为胆立体投影之前,一切都很好,我可以轻松地绘制它。输出很漂亮。之后……没有那么多。
我手头没有数据可以准确地向您展示我在做什么,但是以下代码行应该可以使用任何 spdf 重现。我很确定我在做一些非常明显的错误 - 对任何熟悉的人来说都很明显leaflet()
- 但是关于投影的在线资源是有限的,而且我对地图投影和包的熟悉程度 - 更是如此。无论如何,这就是我正在做的事情:
spdf2 <- spTransform(spdf, CRS = CRS("+proj=gall"))
sp::proj4string(spdf2)
leaflet(spdf2, options =
leafletOptions(
crs=leafletCRS(crsClass="L.Proj.CRS",
code="ESRI:54016",
proj4def= "+proj=gall +ellps=WGS84",
resolutions = c(65536, 32768, 16384, 8192, 4096, 2048)
))) %>% addTiles()
也试过
addGraticule(sphere = TRUE, style= list(color= '#777', weight= 1, opacity= 0.25))
而不是addTiles
,但我不确定它做了什么,也没有帮助。显然,我想要addPolygons
之后,但现在我只想要地图;我可以自己弄清楚多边形。
预期的:
胆立体格式的漂亮地图
结果:
空空如也。
我希望这是由于该addTiles
功能,它不应该为您提供带有胆投影的地图。如果是这样,任何人都可以指导我找到解决方案,例如在哪里可以找到使用胆投影的地图?
谢谢,伙计们,我真的很感激,我一直在努力让它工作几个小时,似乎除了把我的狗献给撒旦之外,我什么都做。对不起,如果我有点浪费你的时间,我只是新手leaflet()
,crs
一般来说。这是一个非常酷的包。