我正在使用新版本的NDVI(物候学中使用的植被指数)数据集来自AVHRR(gimms3g.v1),我需要计算定义为“良好价值”的像素数。
数据集有3个质量标志:
Flag0=good
值(我需要的值)Flag1=spline
插值Flag3=retrieved
图片有9331200,涵盖了整个世界。仅绘制flag0,我显然可以看到更少的像素(因为我没有可视化flag1和flag2),但是当我计算像素数时,R总是给我9331200。
我尝试了length()
,ncell()
,unique()
,但我仍然无法做到。问题是,使用例如str()
,我无法找到如何在数据集中定义标志(格式为Netcdf,version4)。
答案 0 :(得分:0)
我知道有一个更好的解决方案,但是我现在能想到的只是以下几点:
summary(values(yourraster==flag0))
应该为您提供关于flag0标准的“ FALSE”和“ TRUE”像素的计数。