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

私の次の要件は、最初のプロットの上に 2 番目のプロットをオーバーレイすることです。基本的に、地区ごとの電力スコアの上に塔を表示します。

これを達成するために、どのパッケージ (およびサンプル コード) を使用できるかについてのアイデアはありますか?

4

1 に答える 1

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: バージョン 1.6 が必要で、今日 CRAN に提出しました:) 開発バージョンをインストールすることもできます:

library(devtools)
install_github("mtennekes/tmap", subdir = "pkg")
于 2016-10-21T19:52:56.037 に答える