El objetivo es crear una "matriz mascara" a partir de una matriz conocida de números reales con dimensiones [H x W].
La matriz máscara estará compuesta de BLOCK_SIZES
, donde cada BLOCK_SIZE será de dimensiones BSxBS
. Las dimensiones de la matriz máscara serán de [BS x n·BS]
donde n = Nº de BLOCK_SIZES.
Supongamos que se tiene la siguiente matriz:
|5 0 0 0 2 0 0 0|
|0 0 0 0 0 0 0 0|
|0 0 3 0 0 0 0 0|
A = |0 4 0 0 0 0 0 0|
|0 0 0 0 0 0 0 0|
|0 0 0 6 0 0 0 0|
|0 0 0 0 0 0 0 0|
|0 0 0 0 0 0 0 8|
La nueva matriz que se desea es:
B = |5 0 2 0 0 0 3 0 0 0 0 0|
|0 0 0 0 0 4 0 0 0 6 0 8|
Donde se aprecia que se han cogido bloques con N = 2
y al menos uno de sus elementos es no nulo (condición de bloque no nulo), entendiendo como nulo el 0 a la hora de realizar el producto matricial.
El código que yo tengo recorre la matriz original y crea una máscara con ceros y unos en función de si el bloque obtenido cumple la condición de bloque no nulo:
def split(array, nrows, ncols):
r, h = array.shape
r_, h_ = int(r / nrows), int(h / ncols)
mask = np.zeros((r_, h_))
for i in range(r_):
for j in range(h_):
if np.sum(array[nrows*i:nrows*(i+1),
ncols*j:ncols*(j+1)]) != 0:
mask[i, j] = 1
return mask
La matriz máscara sería (Con el mismo algoritmo puedo generar facilmente la matriz B
):
|1 0 1 0|
mask = |1 1 0 0|
|0 1 0 0|
|0 0 0 1|
El problema de mi método es que es tremendamente ineficaz, y como entrada (sin exagerar) tengo matrices de 16000x17000 y dimensiones similares o mayores.
Es posible realizar esto mediante matrices dispersas, como podría hacer lo mas eficiente posible este algoritmo?