Question

到目前为止，我已经设法在 RStudio tmap 包的帮助下，在巴基斯坦的地区层面直观地呈现不同社会经济因素的数据。 >创建以下输出。

Access to Water

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 是我的空间数据框的名称，而条件是该数据框中显示对水的访问权限的变量。

与此同时，我使用 ggplot2 包绘制了巴基斯坦不同移动网络运营商（大约13000个塔楼）的各种塔楼的地理坐标。输出如下。

Tower Presence

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

我没有你的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")

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

1 个答案: