如何在python中找到3D数组的局部最大值？

import numpy as np

arr = np.array([[[1, 2, 3], [4, 5, 6], [7, 8, 9]],
               [[10, 11, 12], [13, 14, 15], [16, 17, 18]],
               [[19, 20, 21], [22, 23, 24], [25, 26, 27]]])

local_max = np.maximum(np.maximum(arr[:-1, :-1, :-1], arr[:-1, :-1, 1:]), np.maximum(arr[:-1, 1:, :-1], arr[:-1, 1:, 1:]))
local_max = np.maximum(local_max, np.maximum(np.maximum(arr[1:, :-1, :-1], arr[1:, :-1, 1:]), np.maximum(arr[1:, 1:, :-1], arr[1:, 1:, 1:])))
max_value = np.max(local_max)
max_index = np.unravel_index(np.argmax(local_max), local_max.shape)

print("局部最大值:", max_value)
print("最大值索引:", max_index)