我的暑假翻车史一：如何生成年NDVI平均值

我用的是MOD13Q1的16d数据，数据下载参考https://blog.csdn.net/weixin_43519457/article/details/95483905，最后生成的是5.1-9.30的NDVI平均值，年平均值大同小异。
（1）MRT融合

（2）每幅影像都有一个-3000，是无效值，那这个-3000应该如何处理，我做了两种，一种是保留这个-3000，另一种是将其设为无效值并且在计算时忽略无效值，发现最后结果是有差别的。查阅资料过后我认为应该采用第二种的方法，工具用的是spatial analyst工具——局部分析——像元统计——mean，无效值用的是栅格计算器中的setnull。

明显这一切都是在ARCGIS里做的，最后导致了我的翻车，还翻了一上午。
翻车史：envi里利用band math将band<0.1和band>0.9的附成了0，但背景值也变成了0，而且死活去不掉，出现这种问题我认为是ENVI里并没有nodata，只有NAN和INF，所以我又把它弄回到arcgis重新剔除了不需要的区域
公式con(raster,0,raster,“value<0.1|value>0.9”)

(3)解决方案二：首先将-3000设为NAN，计算时忽略无效值，再在ENVI波段计算应该不会出现背景值参与计算的问题，我试了再补充。