9

R のパッケージを使用してraster生成されたマップがあります。これらのラスター レイヤーは、パッケージの関数rasterを使用して視覚化できます。rasterVislevelplot

levelplot(rasterstack, layout=c(1, 2), 
          col.regions=colorRampPalette(c('darkred', 'red3', 'orange2', 'orange', 
                                         'yellow', 'lightskyblue', 'steelblue3', 
                                         'royalblue3', 'darkblue')))

ここで、xy 座標で定義された z 値を levelplot マップに追加したいと思います。z 値を含むデータフレームには 4 つの列があります。列 1 & 2 には x & y 座標が含まれ、列 3 にはレイアウト (1, 1) のマップ 1 の z 値が含まれ、列 4 にはレイアウト (1, 2) の z 値が含まれます。

マップごとのポイントは、z < 0.05 のpch=2場合、z > 0.05 の場合、pch=3.

Web を検索したところ、Ripley による解決策が見つかりましたが、私の場合は機能しません。

levelplot(rcp852, xlab = "", ylab = "",
          panel = function(x, y, subscripts, ...) {
            panel.levelplot(x, y, subscripts, ...)
            panel.xyplot(topo$x,topo$y, cex = 0.5, col = 1)
          }
)

他の多くのオプションを試しましたが、ポイントが で生成されたマップと一致しませんlevelplot

4

2 に答える 2