In Python voglio trovare un modo veloce (senza ciclo for) per assegnare indici di array di riferimento. Questo è il risultato desiderato utilizzando un ciclo for:Assegnazione contemporanea di indici di array identici in Python/Numpy
import numpy as np
a=np.arange(9, dtype=np.float64).reshape((3,3))
# The array indices: [2,3,4] are identical.
Px = np.uint64(np.array([0,1,1,1,2]))
Py = np.uint64(np.array([0,0,0,0,0]))
# The array to be added at the array indices (may also contain random numbers).
x = np.array([.1,.1,.1,.1,.1])
for m in np.arange(len(x)):
a[Px[m]][Py[m]] += x
print a
%[[ 0.1 1. 2.]
%[ 3.3 4. 5.]
%[ 6.1 7. 8.]]
Quando provo ad aggiungere x
-a
gli indici Px,Py
io, ovviamente, non si ottiene lo stesso risultato (3.3 vs 3.1):
a[Px,Py] += x
print a
%[[ 0.1 1. 2.]
%[ 3.1 4. 5.]
%[ 6.1 7. 8.]]
C'è un modo per farlo con Numpy? Grazie.
Prima di tutto, vorrei raggruppare i valori insieme, in modo da avere una lista di tuple (Px, Py). Di ordinare questo elenco, contare le occorrenze, moltiplicare x con quel numero e aggiungerlo all'array. Ma in qualche modo intorpidito sembra saltare i doppi ingressi .... strano. – Dschoni