Векторизованное движущееся окно в 2D-массиве в numpy

Я применяю операцию к движущемуся окну постоянного размера в 2D-массиве. Есть ли эффективная операция, подобная векторизации, которую я могу реализовать, чтобы сделать это без цикла в Python? Моя текущая структура выглядит примерно так

 for i in range(1,xmax-1):
     for j in range(1,ymax-1):
        out[i][j] = f(in[i][j],in[i+1][j],in[i-1][j],in[i][j+1],in[i][j-1],...)

Комментарии, которые eat оставили в этом вопросе, намекают на возможность векторизации этой операции this, но без дополнительных деталей векторизованное индексирование / нарезка в numpy / scipy?