为什么当我通过 rasterio 或 gdal 加载 NDVI 图像时所有数据都是 Nan

Question

我在此页面之后计算了 NDVI 值并保存了单波段图像； https://developers.planet.com/docs/planetschool/calculate-an-ndvi-in-python/ 。 之后我想通过rasterio或gdal将这些图像加载到python上，加载图像成功，但不幸的是所有值都是Nan。 为什么所有的值都加载为 Nan？ 以及如何解决这种情况？

我用谷歌搜索但找不到解决方案。 我需要你的帮助，请...

Answer 1

您可以按照此答案的建议使用掩码数组：

mar = np.ma.masked_array(ar, np.isnan(ar))
print(mar.min(), np.median(mar), mar.mean(), mar.std(), mar.max())

这个掩码数组允许在不考虑无效值的情况下对数组进行计算。