0

levelplotパッケージの関数を使用して、RasterVis1 つのラスターをプロットしています。GeoTiff で次のコードを使用します。

require(raster)
require(rasterVis)
data <- raster("mytiff.tif")
levelplot(data, layers=1, par.settings=RdBuTheme)

デフォルトでは、余白に表示されるグラフは何を示していますか? それは各列/行の平均ですか、中央値ですか、累積数ですか?

ヘルプ情報でこれを見つけることができないので、どんな啓発も大歓迎です。

ありがとう!

4

1 に答える 1

3

これらの図は、 の行と列の要約です RasterLayerFUN.margin要約は、 (meanデフォルト値として使用する)で定義された関数で計算さ れます。

例を挙げて説明しましょう:

library(raster)
library(rasterVis)

f <- system.file("external/test.grd", package="raster")
r <- raster(f)

levelplot(r)

レベルプロット

余白に表示されているグラフィックは、zonal 関数で作成できます。行と列の番号を使用して 2 つ作成し、要約するゾーンを定義しますinitRasterLayer

rows <- init(r, v='row')
cols <- init(r, v='col')

rAvg <- zonal(r, rows, fun='mean')
cAvg <- zonal(r, cols, fun='mean')

結果は同じです。

plot(rAvg, type='l')

行平均

plot(cAvg, type='l')

列平均

于 2013-09-03T14:19:50.273 に答える