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包装器，专为这些类型的操作而设计。虽然它下面是元素与元素的比较（如上所述）;在这些类型的操作中它明显更快。