如何在tmap输出上绘制坐标？

Question

到目前为止，我已经在tmap包的帮助下，成功地在巴基斯坦的地区一级直观地展示了不同社会经济因素的数据，以创建以下输出：

​

​

accesstowater <- tm_shape(pak) + 
    tm_fill("Condition", n = 5, palette = "-Blues", colorNA = "grey",
             textNA = "Missing Values", title = "Access to Water") + 
    tm_borders(alpha = 0.5) +
    tm_credits("Source: PBS - PSLM 2014 - 15\nDH Corp Ltd.", size = .64,
                position = c(.62,0.02))

在上面的代码中，pak是我的空间数据的名称，而Condition是该数据帧中的一个变量，显示了对水的访问。

同时，我使用ggplot2软件包绘制了巴基斯坦不同移动网络运营商(大约13000个塔)的各个塔的地理坐标。产出如下：

​

​

pakTowers <- ggmap(pakmap6) + 
               geom_point(aes(Longitude, Latitude, colour = MNO), 
                          data = Towers, size = 0.7, alpha = 0.5,
                          show.legend = FALSE)

在上面的代码中，pakmap6是get-map()函数的输出，MNO表示移动网络运营商，Towers数据集包括每个塔的坐标条目。

我的下一个要求是将第二个地块覆盖在第一个地块上，实质上是显示高楼在地区智能电力得分的顶部。

为了实现这一点，我可以使用哪些包(以及示例代码)，有什么想法吗？

Answer 1

我没有您的shapefiles，所以我无法复制它，但是只使用tmap绝对是可能的：

pak_osm <- read_osm(pak, ext=1.1) # reads OSM data based on the bounding box of pak with 10 percent extension

qtm(pak_osm) +
tm_shape(pak) + 
tm_fill("Condition", n = 5, palette = "-Blues", colorNA = "grey", 
    textNA = "Missing Values", title = "Access to Water") + 
tm_borders(alpha = 0.5) + 
tm_credits("Source: PBS - PSLM 2014 - 15\nDH Corp Ltd.", size = .64, position = c(.62,0.02)) +
tm_shape(Towers) +
tm_symbols(color="MNO", shape="MNO", alpha=0.5)

如果这是你想要的，请告诉我。

PS:它需要1.6版本，我今天向CRAN提交了该版本：)您还可以安装开发版本：

library(devtools)
install_github("mtennekes/tmap", subdir = "pkg")