Создание 4D-представления на 2D-массиве для разделения его на ячейки фиксированного размера

У меня есть 2D-массив t в numpy:

>>> t = numpy.array(range(81)).reshape((9,9))
>>> t
array([[ 0,  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, 28, 29, 30, 31, 32, 33, 34, 35],
       [36, 37, 38, 39, 40, 41, 42, 43, 44],
       [45, 46, 47, 48, 49, 50, 51, 52, 53],
       [54, 55, 56, 57, 58, 59, 60, 61, 62],
       [63, 64, 65, 66, 67, 68, 69, 70, 71],
       [72, 73, 74, 75, 76, 77, 78, 79, 80]])

Он индексируется двумя числами: индексом строки и столбца.

>>> t[2,3]
21
>>> t.shape
(9, 9)
>>> t.strides
(72, 8)

Я хочу разделить массив на прямоугольные ячейки фиксированный размер , например 3 × 3. Я бы хотел избежать копирования памяти. Я пытаюсь добиться этого, создавая представление на t с соответствующей формой и шагов ( (3,3,3,3) и (216,24,72,8) соответственно). Таким образом, первые два индекса представления будут означать положение Ячейка 3 × 3 в большей сетке, а последние две будут означать положение элемента внутри ячейки. Например, t [0,1,:,:] вернет

array([[ 3,  4,  5],
       [12, 13, 14],
       [21, 22, 23]])

Итак, мой вопрос: - как создать описанный вид? Мне не хватает более простого метода? Можно ли это сделать элегантно с синтаксисом нарезки?

8
задан ulidtko 31 March 2016 в 08:14
поделиться