This:

     np.zeros((image.shape[0] + 1, image.shape[1]), order='C')

could be integrated in the _slanted_integral_image function as it is the same for all calls.