1

到目前为止，我已经成功地在包的帮助下直观地呈现了巴基斯坦地区不同社会经济因素的数据，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数据集包括每个塔坐标的条目。

我的下一个要求是将第二个地块覆盖在第一个地块上，本质上是在地区明智的电力分数之上显示塔。

关于我可以使用哪些包（以及示例代码）来实现这一点的任何想法？

4

1 回答 1

1

我没有你的 shapefile，所以我无法复制它，但绝对可以单独使用 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")
于 2016-10-21T19:52:56.037 回答