在scipy.ndimage中应用旋转后获得不同的范围

Question

我有一个在0,3范围内的随机矩阵A。对于这个矩阵，我将以10的角度旋转。然而，在旋转矩阵中，我得到了不同的范围。我的代码有什么问题？谢谢

from scipy.ndimage import rotate
import numpy as np
A = np.random.randint(0, 4, size=(64,128,128))
print (A.max(), A.shape)
B=rotate(A, angle=10, reshape=False, axes=(0,1))
print (B.max(), B.shape)

输出是

(3, (128, 128, 128))
(4, (128, 128, 128))

Answer 1

scipy.ndimage.rotate (手册)使用样条插值执行矩阵旋转，默认情况下使用样条3阶，这可能导致输出值大于输入数据。

您可以将使用的样条插值的顺序降为1，这将行为更改为线性插值。为此，修改旋转行如下：

B=rotate(A, angle=10, order=1, reshape=False, axes=(0,1))

有关样条插值的更多信息，请参见此介绍性论文。