在栅格tiff文件中查找次最小值

Question

如何在光栅中找到第二个最小值？我有一个光栅，我想计算值的范围。最小值是0，所以我想取第二个最小值。

Example data
set.seed(42)
require(raster)
r1 <- raster(nrows=10, ncols=10)

Answer 1

你应该能够对这些值进行排序，然后选择你喜欢的索引。例如

sort(values(r1))[2]

Answer 2

如果您不能将所有值加载到RAM中(正如MrFlick的解决方案所要求的那样)，这可能会起作用：

keep <- c(Inf, Inf)
f <- function(x) {
   sort(c(keep, na.omit(x)))[1:2]
}

library(terra)
r <- rast()
values(r) <- 1:ncell(r)
g <- global(r, f)
g
#      global_1 global_2
#lyr.1        1        2

然后使用第二个值