This is a matrix-vector product. Alternatively, if we are viewing the rgb vector as a row vector then we would compute: img[i,j,:] @ A.T. or using dot we could do. img[i,j,:].dot(A.T) The latter is ...