无法弄清楚如何用胆立体投影在传单r中绘制数据

时间:2019-04-15 03:00:24

标签: r map-projections r-leaflet proj4 spdf

我有一个空间多边形数据框,其坐标已转换为胆立体图类型。但是我现在无法用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的新手。不过这是一个非常酷的程序包。

0 个答案:

没有答案