for i in range(1, np.shape(data)[0] - 1):
    for j in range(1, np.shape(data)[1] - 1):
        if np.all(data[i - 1:i + 2, j - 1:j + 2] == data[i, j]):
            pass
Retrieve position of elements with setting some criteria in numpy