¿Cómo hago esta indexación de matrices en numpy

Question

Dada una matriz de índices I, cómo puedo establecer los valores de una matriz de datos D cuyos índices no existen en I?

Ejemplo: ¿Cómo obtengo A de I y D?

I = array([[1,1], [2,2], [3,3]]) 

D = array([[ 1, 2, 3, 4, 5, 6], 
      [ 7, 8, 9, 1, 2, 3], 
      [ 4, 5, 6, 7, 8, 9], 
      [ 1, 2, 3, 4, 5, 6], 
      [ 7, 8, 9, 1, 2, 3]]) 

A = array([[ 0, 0, 0, 0, 0, 0], 
      [ 0, 8, 0, 0, 0, 0], 
      [ 0, 0, 6, 0, 0, 0], 
      [ 0, 0, 0, 4, 0, 0], 
      [ 0, 0, 0, 0, 0, 0]]) 

Editar: Estoy buscando la manera de hacerlo de una sola vez para los casos en Id y son grandes.

Answer 1

solución simple:

A = zeros(D.shape) 
for i, j in I: 
    A[i, j] = D[i, j] 

vectorizado:

A = zeros(D.shape) 
i, j = I.T 
A[i, j] = D[i, j]