我有一个空间多边形数据框,其坐标已转换为胆立体图类型。但是我现在无法用leaflet()
映射该spdf。我不确定我提供的分辨率是否有问题(因为我不知道这些是什么或做什么),或者与我通过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
的新手。不过这是一个非常酷的程序包。