2016-10-20 77 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")