将bitstring numpy数组转换为整数base 2的最快方法

时间:2016-05-22 10:15:17

标签: python arrays numpy int bitstring

我有一个由bittrings组成的numpy数组,我打算将bittrings转换为整数base 2,以便执行一些xor按位运算。我可以用python将字符串转换为带有基数2的整数:

int('000011000',2)

我想知道是否有更快更好的方法在numpy中做到这一点。我正在研究的numpy数组的一个例子是这样的:

array([['0001'],
       ['0010']], 
      dtype='|S4')

我希望将其转换为:

array([[1],[2]])

2 个答案:

答案 0 :(得分:2)

可以使用np.fromstring将每个字符串位分离为uint8类型数字,然后使用一些带有矩阵乘法的数学来转换/缩减为十进制格式。因此,使用A作为输入数组,一种方法就是这样 -

# Convert each bit of input string to numerals
str2num = (np.fromstring(A, dtype=np.uint8)-48).reshape(-1,4)

# Setup conversion array for binary number to decimal equivalent
de2bi_convarr = 2**np.arange(3,-1,-1)

# Use matrix multiplication for reducing each row of str2num to a single decimal
out = str2num.dot(de2bi_convarr)

示例运行 -

In [113]: A    # Modified to show more variety
Out[113]: 
array([['0001'],
       ['1001'],
       ['1100'],
       ['0010']], 
      dtype='|S4')

In [114]: str2num = (np.fromstring(A, dtype=np.uint8)-48).reshape(-1,4)

In [115]: str2num
Out[115]: 
array([[0, 0, 0, 1],
       [1, 0, 0, 1],
       [1, 1, 0, 0],
       [0, 0, 1, 0]], dtype=uint8)

In [116]: de2bi_convarr = 2**np.arange(3,-1,-1)

In [117]: de2bi_convarr
Out[117]: array([8, 4, 2, 1])

In [118]: out = str2num.dot(de2bi_convarr)

In [119]: out
Out[119]: array([ 1,  9, 12,  2])

可以建议另一种方法来避免np.fromstring。使用此方法,我们将在开始时转换为int数据类型,然后将每个数字分开,这应该等同于前一个方法中的str2num。其余的代码将保持不变。因此,另一种实现方式是 -

# Convert to int array and thus convert each bit of input string to numerals
str2num = np.remainder(A.astype(np.int)//(10**np.arange(3,-1,-1)),10)

de2bi_convarr = 2**np.arange(3,-1,-1)
out = str2num.dot(de2bi_convarr)

运行时测试

让我们列出迄今为止列出的所有方法来解决问题,包括@Kasramvd's loopy solution

In [198]: # Setup a huge array of such strings
     ...: A = np.array([['0001'],['1001'],['1100'],['0010']],dtype='|S4')
     ...: A = A.repeat(10000,axis=0)


In [199]: def app1(A):             
     ...:     str2num = (np.fromstring(A, dtype=np.uint8)-48).reshape(-1,4)
     ...:     de2bi_convarr = 2**np.arange(3,-1,-1)
     ...:     out = str2num.dot(de2bi_convarr)    
     ...:     return out
     ...: 
     ...: def app2(A):             
     ...:     str2num = np.remainder(A.astype(np.int)//(10**np.arange(3,-1,-1)),10)
     ...:     de2bi_convarr = 2**np.arange(3,-1,-1)
     ...:     out = str2num.dot(de2bi_convarr)    
     ...:     return out
     ...: 

In [200]: %timeit app1(A)
1000 loops, best of 3: 1.46 ms per loop

In [201]: %timeit app2(A)
10 loops, best of 3: 36.6 ms per loop

In [202]: %timeit np.array([[int(i[0], 2)] for i in A]) # @Kasramvd's solution
10 loops, best of 3: 61.6 ms per loop

答案 1 :(得分:0)

由于KISS principle,我想建议使用列表理解的以下方法:

>>> np.array([[int(i[0], 2)] for i in a])
array([[1],
       [2]])