使用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的基本上是一个Python用于底层C代码库的封装器，专为这些操作而设计。虽然它下面是比较元素元素（如你所说的）;这些操作的速度要快得多。