3

格子パッケージを使用して、R で変数のレベル プロットを作成しました。このグリッドは南アジアに対応します。南アジアの特定の国について、この変数 (エアロゾルの光学的深さ) の値を表示することにのみ関心があります。関心のある国では値 1 を取り、それ以外の場合は 0 を取るダミー変数があります。グリッドのこの部分を黒または他の色にすることはできますか?

stackoverflow の評価が低いため、レベル プロットを表示できません。(rhelp へのクロスポスト メッセージに添付されていた pdf が表示されるはずです:)

ここに私のRコードがあります:

levelplot(aod ~ longitude + latitude | factor(day), data = aod_Jan, 
    aspect="iso", contour = TRUE, layout=c(1,1))

ここに画像の説明を入力

4

2 に答える 2

4

地理データを使用しているため、rasterパッケージが役立つ可能性があります。たとえば、フランスの高度を表示してみましょう (このzipファイルをダウンロードするか、raster::getData関数を使用します)。unzipファイルの後：

library(raster)
fraAlt <- raster('FRA_alt')
plot(fraAlt) ## Not only France is displayed...

フランスの高度のみを表示したい場合は、境界の情報が必要です: このRDataファイルをダウンロード (またはraster::getData関数を使用)。これには、次のように変換できる(という名前の)RDataが含まれています。SpatialPolygonsDataFramegadmRaster

 mk <- rasterize(gadm, fraAlt)

これで、標高ラスターを境界でマスクできます。

 fraAltMask <- mask(fraAlt, x)
 plot(fraAltMask) ##Now only France is displayed

最後に、メソッドを使用する場合は、パッケージlatticeが必要です。rasterVis

 library(rasterVis)
 levelplot(fraAlt)
 levelplot(fraAltMask)

ここで、境界をすべて重ね合わせます。

 s <- stack(fraAlt, fraAltMask)
 layerNames(s) <- c('Alt', 'AltMask')
 boundaries <- as(gadm, 'SpatialLines')
 levelplot(s) + layer(sp.lines(boundaries))
于 2011-09-28T15:10:41.857 に答える
0

levelplot には、subset 引数を使用します。多分：

levelplot(aod ~ longitude + latitude | factor(day), data = aod_Jan, subset = dummy==1,
    aspect="iso", contour = TRUE, layout=c(1,1))
于 2011-09-27T17:46:51.330 に答える