1

我有一个空间多边形数据框,其坐标已转换为胆立体类型。但是我现在无法映射该 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一般来说。这是一个非常酷的包。

4

0 回答 0