2014-10-05 126 views
3

我使用numpy和scipy從3D座標信息生成密度圖。我可以通過生成KDE用下面的代碼Python - 獲取最密集點的座標

xyz = np.vstack([x,y,z]) 
kde = stats.gaussian_kde(xyz) 
density = kde(xyz) 

成功生成數據的密度圖,但我怎麼能利用這些信息來發現具有最大密度的三維點相關聯的座標?

我已經試過

max(density) 

它返回一個值,然後我就可以找到與

density.argmax(axis=0) 

但後來我打了一個空白的指數,因爲我似乎無法使用索引從xyz獲取關聯的座標,我不確定這是否是正確的方法。

回答

2

here,我可以用

xyz.T[np.argmax(density)] 

返回我的數據

最密集的點的三維座標