在r中计数相同值的像素

Question

我需要在R中计算相同值的像素的帮助。我正在做一个主论文，其中一部分是在洪水之前和之后计算一个区域的NDWI，因此我可以计算出水下的表面。从高于0之后的洪水值之后的NDWI的结果中，在水下，低于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

在这种情况下总和这些区域，这类似于（错误）假设 naminal naminal 空间分辨率对每个网格单元都是有效的。

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