2

多くのパッチ（同じ値の連続するセルの塊）を含むラスターマップがあります。私がする必要があるのは、各パッチの中心（または中心に近い）の座標を取得することです。私はラスターパッケージに非常に不慣れですが、マップ内のセルの位置を知っている場合にのみ座標を取得できるようです。代わりにセルの値を与える座標を取得する方法はありますか？ありがとうございました

4

2 に答える 2

5

パッチで凝集塊を意味する場合、ラスターパッケージを使用すると、凝集塊を見つけて分離できます。ラスターパッケージの例を取り上げてclump()、それを拡張します。

library(raster)
library(igraph)
detach("package:coin", unload=TRUE)
r <- raster(ncols=12, nrows=12)
set.seed(0)
r[] <- round(runif(ncell(r))*0.7 )
rc <- clump(r)
clump_id <- getValues(rc)    
xy <- xyFromCell(rc,1:ncell(rc))
df <- data.frame(xy, clump_id, is_clump = rc[] %in% freq(rc, useNA = 'no')[,1])
df[df$is_clump == T, ]

plot(r)

r

plot(rc)
text(df[df$is_clump == T, 1:2], labels = df[df$is_clump == T, 3])

rc

あなたが期待できるほど面白くないかもしれません。

あなたは方向=4でそれをすべてやります

rc <- clump(r, directions = 4)
clump_id <- getValues(rc)    
xy <- xyFromCell(rc,1:ncell(rc))
df <- data.frame(xy, clump_id, is_clump = rc[] %in% freq(rc, useNA = 'no')[,1])
df[df$is_clump == T, ]

取得するため

rc4

そして多分「重心」の塊

dfm <- ddply(df[df$is_clump == T, ], .(clump_id), summarise, xm = mean(x), ym = mean(y))
plot(rc)
text(dfm[, 2:3], labels = dfm$clump_id)

rc4mean

注：

最初にmodeltoolsライブラリclump()をデタッチせずに使用しようとすると、エラーが発生します。modeltoolsは、coinおよびおそらく他の統計ライブラリによって呼び出されます。

于 2013-03-26T10:54:12.343 に答える
2

各パッチの座標の平均を取ることができます：

# some dummy data
m <- matrix(c(
    0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
    0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
    0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,
    0,0,0,0,0,0,0,0,0,1,1,1,0,0,0,0,0,0,0,4,4,0,
    0,0,0,0,0,0,0,0,1,1,1,1,1,0,0,0,0,0,0,4,4,0,
    0,0,0,0,0,0,0,1,1,1,1,1,1,1,0,0,0,0,0,4,4,0,
    0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,
    0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,
    0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,
    0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,
    0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,
    0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,
    0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,
    0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,
    0,0,0,0,0,0,0,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,
    0,0,0,0,0,0,0,0,1,1,1,1,1,0,0,0,0,0,0,0,0,0,
    0,0,0,0,0,0,0,0,0,1,1,1,0,0,0,0,0,0,0,0,0,0,
    0,0,2,3,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,
    0,0,2,3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
    0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0), nrow=20, byrow=T)

# create a raster
r <- raster(m)

# convert raster to points
p <- data.frame(rasterToPoints(r))

# filter out packground
p <- p[p$layer > 0,]

# for each patch calc mean coordinates
sapply(split(p[, c("x", "y")], p$layer), colMeans)
于 2013-03-26T08:55:08.447 に答える