Buen día,
Podrías utilizar una combinación de numpy.reshape
con numpy.swapaxes
para crear los bloques, pero hay otro problema, como puse en el comentario, hay casos que la matriz (M, N) no sea divisible en el tamaño de bloques (m, n) configurado.
Un ejemplo sería, teniendo una matriz de 4 por 4, al intentar separarla en bloques de 2 por 2 no habría problemas, pero al intentar separarla en bloques de 3 por 3 únicamente podríamos hacer el primer bloque porque no hay elementos suficientes. Otra opción que creo sería mejor, es agregar filas y columnas con un valor determinado que pueda ser elegido por el usuario.
Hice una función que realiza lo comentado anteriormente
Los métodos numpy.r_
y numpy.c_
se utilizan para concatenar filas (Rows) y columnas respectivamente.
Los bloques se crean de izquierda a derecha y de arriba hacia abajo
import numpy as np
def bloques(matriz, nrows, ncols, filler):
#Dividir una matriz en sub matrizes
r, h = matriz.shape
#Validar si la matriz es divisible por el número de filas y columnas
#En caso contrario agregar filas y/o columnas según sea necesario
while h % nrows != 0:
# Agregar columnas utilizando "filler" para cada nuevo elemento
matriz = np.c_[matriz, np.full(r, filler)]
r, h = matriz.shape
while r % ncols != 0:
#Agregar filas utilizando "filler" para cada nuevo elemento
matriz = np.r_[matriz, [np.full(h, filler)]]
r, h = matriz.shape
# Regresar los bloques de tamaño [nrows, ncols]
return (matriz.reshape(h//nrows, nrows, -1, ncols)
.swapaxes(1, 2)
.reshape(-1, nrows, ncols))
matriz = np.random.rand(4,4)
print(matriz)
b = bloques(matriz, 3, 3, 0)
for i in b:
print(i)
Utilizando una matriz de ejemplo de 4 por 4
[[0.52415764 0.40264294 0.05015354 0.18056611]
[0.22324492 0.76061367 0.83139778 0.71840166]
[0.64091326 0.89751854 0.4006832 0.22911819]
[0.38059384 0.12351579 0.72566029 0.34900971]]
Obtenemos los siguientes 4 bloques de 3 por 3
#Bloque 1
[[0.52415764 0.40264294 0.05015354]
[0.22324492 0.76061367 0.83139778]
[0.64091326 0.89751854 0.4006832 ]]
#Bloque 2
[[0.18056611 0. 0. ]
[0.71840166 0. 0. ]
[0.22911819 0. 0. ]]
#Bloque 3
[[0.38059384 0.12351579 0.72566029]
[0. 0. 0. ]
[0. 0. 0. ]]
#Bloque 4
[[0.34900971 0. 0. ]
[0. 0. 0. ]
[0. 0. 0. ]]