2016-03-29 26 views
0

I hist.max() (x, y) bilmek isteyen, bu bir 2d histogramNumpy: Max bin konumu nasıl alınır?

from mpl_toolkits.mplot3d import Axes3D 
import matplotlib.pyplot as plt 
import numpy as np 

fig = plt.figure() 
ax = fig.add_subplot(111, projection='3d') 
x, y = np.random.rand(2, 100) * 10 
hist, xedges, yedges = np.histogram2d(x, y, bins=20) 

elements = (len(xedges) - 1) * (len(yedges) - 1) 
xpos, ypos = np.meshgrid(xedges[:-1]+0.25, yedges[:-1]+0.25) 

xpos = xpos.flatten() 
ypos = ypos.flatten() 
zpos = np.zeros(elements) 
dx = 0.5 * np.ones_like(zpos) 
dy = dx.copy() 
dz = hist.flatten() 

ax.bar3d(xpos, ypos, zpos, dx, dy, dz, zsort='average') 
plt.show() 

enter image description here

maksimum bin yer bulmak istiyoruz.

Sanırım hist.argmax() ile kurtulabilirim. Ama diğerleriyle nasıl baş edeceğimi (1d pozisyonunu 2d'ye çevirmek) yapmıyorum. Ya da daha iyi çözümler varsa?

cevap

2

np.unravel_index `satır, sütun = np.unravel_index (hist.argmax(), hist.shape)` Bu

+1

Yani bu olurdu çözme en deyimsel yoludur? satır geri dönecek, col, ve sonra ben bir değer alabilirim [xedges [row], yedges [col] ' – cqcn1991

1

Düzleştirilmiş dizideki konumu, sütun sayısıyla bölün ve satır alırsınız. Kalan sütun.

row, col = divmod(np.argmax(hist), hist.shape[0])