如何为MODIS LST栅格创建多值筛选的质量掩码

问题修正 #

你当前代码的最后一行maskqa1[maskqa1 > 17]<-NA会误将65设为NA，因为65满足>17的条件。以下是两种修正方案：

方案一：直接匹配目标值集合（推荐）

这种方式逻辑清晰，避免范围判断错误：

# 复制原QA栅格并处理NA值
maskqa1 <- qa1
maskqa1[is.na(maskqa1)] <- 0

# 仅保留0、17、65，其余值设为NA
maskqa1[!(maskqa1 %in% c(0, 17, 65))] <- NA

方案二：精准范围条件判断

如果需要用范围逻辑实现，可调整为：

maskqa1 <- qa1
maskqa1[is.na(maskqa1)] <- 0

# 屏蔽2、3，以及17到65之间（不含两端）、大于65的值
maskqa1[maskqa1 %in% c(2,3) | (maskqa1 > 17 & maskqa1 < 65) | maskqa1 > 65] <- NA

通用优化方案 #

针对后续处理不同标志栅格的需求，推荐以下两种灵活方案：

• 目标值集合匹配
  把需要保留的标志值存入向量，用%in%做判断，后续修改仅需调整向量内容，扩展性极强：

  # 定义需保留的质量标志值
  keep_values <- c(0, 17, 65)
  maskqa1 <- qa1
  maskqa1[is.na(maskqa1)] <- 0
  maskqa1[!(maskqa1 %in% keep_values)] <- NA
  

• 位掩码解析（适配MODIS QA编码特性）
  MODIS QA标志多为二进制位编码，直接解析位掩码更贴合设计逻辑，避免范围判断误差：

  # 示例：提取QA值的第0位和第4位（对应17的二进制00010001）
  bit0 <- bitwAnd(qa1, 1)    # 第0位对应值1
  bit4 <- bitwAnd(qa1, 16)   # 第4位对应值16
  
  maskqa1 <- qa1
  maskqa1[is.na(maskqa1)] <- 0
  # 保留符合位条件的17，以及0、65
  maskqa1[!( (bit0 == 1 & bit4 == 16) | maskqa1 %in% c(0,65) )] <- NA
  

内容的提问来源于stack exchange，提问作者Martín Francia