R语言处理RasterLayer栅格时如何筛选指定范围像素值

你之前的代码失效核心是条件判断的逻辑完全写错了：

• 第一版代码里的判断条件(R<=0.1) & (R>=0.2)，要求像素值同时满足「小于等于0.1」和「大于等于0.2」，不存在任何实数能符合这个互斥的要求，因此要么匹配不到任何像素，赋值NA的操作完全不生效；如果不小心把符号写反，又会把所有像素都匹配到，导致整个栅格变成NA。
• 第二版单值筛选的代码(R<=0.1) & (R>=0.1)逻辑上等价于R == 0.1，但栅格值默认是浮点型存储，直接用等号判断很容易因为存储精度误差，匹配不到你想要的目标值。

保留0.1~0.2区间像素值 #

核心逻辑是：把不在目标区间内的像素赋值为NA，也就是把「值小于0.1」或者「值大于0.2」的部分设为NA，两个排除条件是互斥关系，必须用或运算符|连接，不能用与运算符&。
如果你用的是经典raster包处理RasterLayer对象，直接按索引赋值即可：

R <- myraster
R[R < 0.1 | R > 0.2] <- NA

如果你用的是新版terra包处理SpatRaster对象，可以直接用内置的clamp函数一步完成，不需要手动写索引判断：

library(terra)
# values=FALSE表示区间外的值直接设为NA，而非截断到0.1或0.2的边界值
R <- clamp(myraster, lower = 0.1, upper = 0.2, values = FALSE)

单值筛选的正确写法 #

如果需要筛选等于某个特定值的像素，不要直接用==做浮点判断，要预留浮点精度误差阈值，比如要保留值为0.1的像素：

R <- myraster
# 允许1e-7的精度误差，匹配存储层面近似等于0.1的像素
R[abs(R - 0.1) > 1e-7] <- NA

注意：做区间筛选时很容易混淆逻辑运算符：要保留[a,b]区间的值，排除的是<a和>b两个独立范围的数值，两个排除条件必须用|连接，用&连接永远无法匹配到有效像素。

内容的提问来源于stack exchange，提问作者Maik