R中同值像素的计数

Question

我需要帮助计算R中相同值的像素，我正在做一篇硕士论文，其中一部分是计算洪水前后某一区域的NDWI，这样我就可以计算水下地表。从NDWI的结果来看，洪水值在0以上的是水下的，0以下的不是。所以我被困在bcs里，我不知道如何计算0以上的所有像素。

Answer 1

使用dropbox数据：

library(terra)
ndwi <- rast('Clip_ndwi_poslije1.tif')
ndwi_vec <- values(ndwi, mat = FALSE)
> length(which(ndwi_vec > 0))
[1] 21157
> length(which(ndwi_vec <= 0))
[1] 13217
> dim(ndwi)[1]
[1] 220
> dim(ndwi)[1]*dim(ndwi)[2]
[1] 58300
> sum(length(which(ndwi_vec > 0)), length(which(ndwi_vec <= 0)))
[1] 34374
> table(is.na(ndwi[])) #  see below

FALSE  TRUE 
34374 23926 

表方法来自SOF-GIS ?s，并让您开始。

Answer 2

你的数据

library(terra)
r <- rast("Clip_ndwi_poslije1.tif")

您可以使用大于零(或不大于零)的值来计数单元格，如下所示：

freq(r > 0)
#     layer value count
#[1,]     1     0 13217
#[2,]     1     1 21157

TRUE == 1，所以您有21,157个单元，其值大于零。

但是，要精确计算这些单元格覆盖的面积，首先需要得到每个单元格的面积：

a <- cellSize(r)

然后把这些区域加起来

b <- ifel(r > 0, a, NA)
# or the equivalent
# b <- mask(a, r>0, maskvalue=0)

global(b, "sum", na.rm=TRUE)
#            sum
#ndwiii 19013036

在这种情况下，这类似于(错误地)假设名义空间分辨率对每个网格单元有效。

global(r > 0, "sum", na.rm=TRUE) * prod(res(r))
#             sum
#ndwiii 19041300