1

ラスターファイルがあります

r <- raster('ras')

とシェープファイル

abys <- readShapeSpatial('abys')

次の方法で、シェープファイルによって定義された平均値を計算しました。

r.vals<- extract(r,abys)
r.mean <- lapply(r.vals,FUN=mean)

ただし、出力を返すときにいくつかのシェープファイルを使用すると、複数の結果が得られます。例：

[[1]]
[1] 9321

[[2]]
[1] 6616

[[3]]
[1] 8348

それは私が通常得るものであるものを返す必要があります。これは、シェープファイルの特性によるものですか、それとも方法論の問題によるものですか？

ご入力いただきありがとうございます

4

1 に答える 1

2

問題は、に3つのポリゴンがあることabysです。

最善の解決策は、結果を平均化するのではなく、最初にポリゴンを結合することです。

library(rgeos)
abys.single <- gUnaryUnion(abys)
r.vals<- extract(r,abys.single)
r.mean <- lapply(r.vals,FUN=mean)
于 2012-10-19T22:03:45.707 に答える