使用 numpy 进行数据比较

Question

我正在尝试仅使用 numpy 制作一个算法（我看到其他人使用 PIL，但它有一些缺点）可以比较和绘制两张显示不同年份冰层的地图之间的差异。我加载图像并将 NaN 设置为零，因为我有一些。

data = np.load(filename)
data[np.isnan(data)]=0

数据数组包含 0 到 100 之间的值，代表浓度水平（100 是深蓝色）。 数据如下所示：

我正在尝试计算差异，以便随着时间的推移，冰的损失将对应于负值，而冰的增加将对应于正值。冰在上图中用蓝色表示。 有什么提示吗？逐个元素比较似乎不是最好的主意...

Answer 1

要获得 2 个相同大小的 numpy 数据数组之间的差异，只需从另一个中取一个：

diff = img1 - img2

Numpy 基本上是底层 C 代码库的 Python 包装器，专为此类操作而设计。虽然在它下面是元素之间的比较（如你上面所说）；在这些类型的操作中速度明显更快。