Rのtmap出力に座標をプロットする方法は？

Question

これまで、tmapパッケージRStudioパッケージの助けを借りて、パキスタンの地区レベルで様々な社会経済的要因のデータを視覚的に表示しました。上記のコードで

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)) 

^、朴は私の空間データフレームの名前であり、条件は、水へのアクセスを示すそのデータフレーム内の変数です。

同時に、パキスタンのさまざまなモバイルネットワークオペレータ（約13000タワー）のさまざまなタワーの地理的座標をプロットするために、ggplot2パッケージを使用しました。出力は以下の通りです。上記のコードで

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マップ関数の出力であり、MNOは、モバイルネットワークオペレータとタワーデータセットが各々のエントリを含む表しますタワーの座標。

次の要件は、2つ目のプロットを最初のプロットにオーバーレイすることです。これは、地区賢明な電力スコアの上にタワーを表示することです。

これを達成するために私が使用できるパッケージ（サンプルコードと同様）に関するアイデアはありますか？

Answer 1

私はあなたのシェープファイルを持っていないので、私はそれを再現することはできませんが、それは間違いなく一人で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：スタックオーバーフローへ

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