R - NDVI3g.v1标记像素数

Question

我正在使用新版本的NDVI（物候学中使用的植被指数）数据集来自AVHRR（gimms3g.v1），我需要计算定义为“良好价值”的像素数量。

数据集有3个质量标志：

• Flag0=good价值观（我需要的）
• Flag1=spline插值
• 来自季节性情况的Flag3=retrieved

图像有9331200，覆盖了整个世界。仅绘制flag0，我可以明显看到更少的像素（因为我没有可视化flag1和flag2），但是当我计算像素数时，R总是给我9331200。

我试过length()，ncell()，unique()，但我仍然不能这样做。问题是，使用例如str()，我无法找到如何在数据集中定义标志（格式为Netcdf，version4）。

Answer 1

我知道有一个更好的解决方案，但我现在能想到的就是：

summary(values(yourraster==flag0))

哪个应该给你与flag0标准相对应的“FALSE”和“TRUE”像素的计数。