python - 如何在屏蔽数组上运行 numpy 函数 percentile()？

Question

我尝试从具有 NoData 值的数组中检索百分位数。在我的例子中，Nodata 值由 -3.40282347e+38 表示。我认为屏蔽数组会将此值排除在进一步计算之外。我成功创建了掩码数组，但对于 np.percentile() 函数，掩码无效。

>>> DataArray = np.array(data)
>>> DataArray

([[ value, value...]], dtype=float32)

>>> masked_data = ma.masked_where(DataArray < 0, DataArray)
>>> p5 = np.percentile(masked_data, 5)
>>> print p5

 -3.40282347e+38

Answer 1

如果将掩码值填充为 np.nan，则可以使用 np.nanpercentile

import numpy as np
data = np.arange(-5.5,10.5) # Note that you need a non-integer array to store NaN
mdata = np.ma.masked_where(data < 0, data)
mdata = np.ma.filled(mdata, np.nan)
np.nanpercentile(mdata, 50) # 50th percentile