np.diff 使用陷阱 数据溢出

I_diff = np.diff(warped.astype(np.int16), axis=1)  # 计算列方向的梯度差分j
计算差分的对象的数据类型时int16.

astype(np.int16)

将图像（或数组）warped的元素转换为int16类型。在差分计算之前进行这个转换，意味着差分结果不会受到溢出问题的影响，因为int16具有比默认uint8更大的数值范围（-32768到32767）。 这对于图像中可能有负值的情况（比如计算梯度或差分时）非常重要，因为uint8只能表示0到255的值，进行差分计算时如果存在溢出，结果可能会不正确。