Инкремент массива Numpy multi-d с повторяющимися индексами - PullRequest
3 голосов
/ 15 сентября 2011

Меня интересует многомерный случай Увеличение массива Numpy с повторяющимися индексами .

У меня есть N-мерный массив и множество N индексных массивов, чьи значения яхочу увеличить.В индексных массивах могут быть повторяющиеся записи.

Без повторов решение:

a = arange(24).reshape(2,3,4)
i = array([0,0,1])
j = array([0,1,1])
k = array([0,0,3])
a[i,j,k] += 1

С повторениями (например, j=array([0,0,2])) я не могу сделать пустое приращениеразмножается.

Ответы [ 2 ]

3 голосов
/ 15 сентября 2011

Как насчет этого:

import numpy as np
a = np.zeros((2,3,4))
i = np.array([0,0,1])
j = np.array([0,0,1])
k = np.array([0,0,3])

ijk = np.vstack((i,j,k)).T
H,edge = np.histogramdd(ijk,bins=a.shape)
a += H  
0 голосов
/ 15 сентября 2011

Я не знаю, есть ли более простое решение с прямой индексацией массива, но это работает:

for x,y,z in zip(i,j,k):
    a[x,y,z] +=1
Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...