Spent a while this morning looking for a generalized question to point duplicates to for questions about as_strided
and/or how to make generalized window functions. There seem to be a lot of questions on how to (safely) create patches, sliding windows, rolling windows, tiles, or views onto an array for machine learning, convolution, image processing and/or numerical integration.
I'm looking for a generalized function that can accept a window
, step
and axis
parameter and return an as_strided
view for over arbitrary dimensions. I will give my answer below, but I'm interested if anyone can make a more efficient method, as I'm not sure using np.squeeze()
is the best method, I'm not sure my assert
statements make the function safe enough to write to the resulting view, and I'm not sure how to handle the edge case of axis
not being in ascending order.
DUE DILIGENCE
The most generalized function I can find is sklearn.feature_extraction.image.extract_patches
written by @eickenberg (as well as the apparently equivalent skimage.util.view_as_windows
), but those are not well documented on the net, and can't do windows over fewer axes than there are in the original array (for example, this question asks for a window of a certain size over just one axis). Also often questions want a numpy
only answer.
@Divakar created a generalized numpy
function for 1-d inputs here, but higher-dimension inputs require a bit more care. I've made a bare bones 2D window over 3d input method, but it's not very extensible.