r - 获取光栅 map 中补丁的坐标(R 中的光栅包)

Question

我有一个包含许多补丁(具有相同值的连续单元 block )的栅格 map 。我需要做的是获取每个补丁的中心(或接近中心)的坐标。我对栅格包非常缺乏经验，但似乎只有当我知道 map 中单元格的位置时我才能获得坐标。有什么方法可以让坐标给出单元格的值吗？谢谢

Answer 1

如果您所说的补丁是指团 block ，Raster 包允许您查找和隔离团 block 。以 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)

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

可能没有你想象的那么有趣。

你用 directions = 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, ]

得到

也许还有团 block “质心”

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)

注意事项:

如果你不先尝试使用 clump() 将会出错分离 modeltools 库。 modeltools 被硬币调用，也许其他统计图书馆。