NumPy的方式
下面是使用advanced indexing
一個矢量NumPy的方式 -
# Extract array data
In [10]: a = df.values
# Get integer based column IDs
In [11]: col_idx = np.searchsorted(df.columns, columns_to_select)
# Use NumPy's advanced indexing to extract relevant elem per row
In [12]: a[np.arange(len(col_idx)), col_idx]
Out[12]: array([ 10, 2, 3, 400])
如果df
列名不排序,我們需要使用sorter
大吵np.searchsorted
。該代碼提取col_idx
,例如通用df
是:
# https://stackoverflow.com/a/38489403/ @Divakar
def column_index(df, query_cols):
cols = df.columns.values
sidx = np.argsort(cols)
return sidx[np.searchsorted(cols,query_cols,sorter=sidx)]
所以,col_idx
會像這樣獲得 -
col_idx = column_index(df, columns_to_select)
進一步優化
剖析它揭示了瓶頸正在處理字符串np.searchsorted
,通常NumPy的弱點與字符串沒有太大關係。所以,爲了克服這個問題,並使用列名爲單個字母的特殊情況,我們可以快速將它們轉換爲數字,然後將它們提供給searchsorted
以加快處理速度。
因此,獲得基於整數列ID,對於列名的單字母排序的情況的優化版本,將是 -
def column_index_singlechar_sorted(df, query_cols):
c0 = np.fromstring(''.join(df.columns), dtype=np.uint8)
c1 = np.fromstring(''.join(query_cols), dtype=np.uint8)
return np.searchsorted(c0, c1)
這給了我們解決方案的修改版,像這樣 -
計時 -
In [149]: # Setup df with 26 uppercase column letters and many rows
...: import string
...: df = pd.DataFrame(np.random.randint(0,9,(1000000,26)))
...: s = list(string.uppercase[:df.shape[1]])
...: df.columns = s
...: idx = np.random.randint(0,df.shape[1],len(df))
...: columns_to_select = np.take(s, idx).tolist()
# With df.lookup from @MaxU's soln
In [150]: %timeit pd.Series(df.lookup(df.index, columns_to_select))
10 loops, best of 3: 76.7 ms per loop
# With proposed one from this soln
In [151]: %%timeit
...: a = df.values
...: col_idx = column_index_singlechar_sorted(df, columns_to_select)
...: out = pd.Series(a[np.arange(len(col_idx)), col_idx])
10 loops, best of 3: 59 ms per loop
鑑於df.lookup
解決了一般情況,這可能是一個更好的選擇,但其他可能的優化,如這篇文章中所示,也可以很方便!
一個你身後秒。 ;-) – Wen
@溫,是的,我知道這種感覺 - 對不起:) – MaxU
@MaxU這正是我所期待的。謝謝! –