使用numpy进行数据比较

Question

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

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

数据数组包含介于0和100之间的值，并表示浓度级别(100表示深蓝色)。数据如下所示：

​

​

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

Answer 1

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

diff = img1 - img2

Numpy基本上是底层C代码库的Python包装器，专为这些类型的操作而设计。尽管在下面它是一个元素一个元素地比较(正如你上面所说的)；但在这些类型的操作中它的速度要快得多。