r - 计算值 <0 的最小值？

问题描述

library(raster)
r1 <- r2 <- r3 <- raster(ncol=10, nrow=10)
r1[] <- runif(ncell(r1))
r2[] <- runif(ncell(r2)) / 2
r3[] <- runif(ncell(r3)) * 1.5
s <- stack(r1, r2, r3)

我可以使用minandmax来获得 s 的最小值和最大值

 res=max(s); res=min(s)

我需要的是仅对大于 0 的值取最大值，对小于 0 的值取最小值。

如果一个像素的所有值 < 0，则返回 NA

标签： rraster

我误解了你的问题，这是修改后的答案

您可以首先编写一个函数来为向量执行此操作

f <- function(x) {
    s <- sign(x)
    if (all(is.na(s))) return( c(NA, NA))
    if (all(s == -1)) return( c(NA, NA))
    mn <- min(x[which(s==-1)])
    mx <- max(x[which(s==1)])
    r <- c(mn, mx)
    r[!is.finite(r)] <- NA
    r
}

并测试它

a <- c(5,-1,2)
b <- c(-5,-1,-2)
f(a)
#[1] -1  5
f(b)
#[1] NA NA
f(c(NA, NA, NA))
#[1] NA NA
f(c(0, 0, 0))
#[1] NA NA

有可以忽略（或关闭）的警告现在，使用栅格数据：

library(raster)
r1 <- r2 <- r3 <- raster(ncol=10, nrow=10)
set.seed(1)    
values(r1) <- runif(ncell(r1), -1, 1)
values(r2) <- runif(ncell(r2), -1, 1) / 2
values(r3) <- runif(ncell(r3), -1, 1) * 1.5
s <- stack(r1, r2, r3)

x <- calc(s, f)
x

#class       : RasterBrick 
#dimensions  : 10, 10, 100, 2  (nrow, ncol, ncell, nlayers)
#resolution  : 36, 18  (x, y)
#extent      : -180, 180, -90, 90  (xmin, xmax, ymin, ymax)
#coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
#data source : in memory
#names       :       layer.1,       layer.2 
#min values  : -1.3561226050,  0.0004409726 
#max values  :   -0.04476076,    1.44469038

或者使用光栅代数，像这样：

pos <- reclassify(s, cbind(-Inf, 0, NA))
neg <- reclassify(s, cbind(0, Inf, NA))
xp <- max(pos, na.rm=TRUE)
nn <- min(neg, na.rm=TRUE)
nn <- mask(nn, xp)
r <- stack(nn, xp)
r

（我不完全确定你想如何处理只有零值的单元格，所以我不担心）

r - 计算值 <0 的最小值？

问题描述

解决方案

推荐阅读