栅格计算器求平均

我的多幅栅格图像存在大量的无效值，所以想做气溶胶月平均数据。然后我打算先将栅格数据转为浮点型并且把无效值-9999赋值为0，用公式（b1+b2+b3…）/((b1>0)+(b2>0)+(b3>0)…)求平均，例如：图像中某个点在第一张图b1=0,在第二张图b2=3,在第三张图b3=5，这样平均值为4。但是这样需要自己手动输入每个图层，不知道有没有更好的办法呢？要怎么样求有效的图层的平均？