2016-10-20 111 views
0

到目前为止,我已经设法在巴基斯坦的一个地区级视觉呈现不同社会经济因素的数据,在RStudio包的帮助下创建以下输出。如何在R中绘制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)) 

^在上面的代码,是我的空间数据帧的名称和条件是数据帧表示获得水内的变量。同时,我用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是和输出的得到映射功能,MNO表示移动网络运营商和数据集包括针对每个条目塔的坐标。

我的下一个要求是将第二个地块叠加在第一个地块上,基本上显示了智能电力分数之上的塔楼。

为了达到这个目的,我可以使用哪些软件包(以及示例代码)?

+0

欢迎堆栈溢出。这是一个有趣的问题,但没有代码和数据很难回答。请编辑您的问题并添加您的代码和示例数据以创建一个最小可重现的示例。谢谢。 – Uwe

+0

嘿,我已经添加了更多的代码数据。我还可以如何让我的问题更全面? –

回答

1

我没有你的shape文件,所以我不能复制它,但它绝对有可能需要单独的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")