代码之家 › 专栏 › 技术社区 › hamster on wheels

快速将阵列分配给等长的n个存储单元

discretization binning numpy python

hamster on wheels · 技术社区 · 7 年前

0.0 到 10.0

我想分配 arr 到 5

等长是指bin间隔为 [0.0, 2.0), [2.0, 4.0), [4.0, 6.0), [6.0, 8.0), [8.0, 10.0] .

问题是最后一个间隔与其他间隔不同。

测试:

import numpy as np
# Things we know and can pre-calculate
n_bins = 5
minimal = 0.0  
maximal = 10.0
reciprocal_bin_length = n_bins / (maximal - minimal)

# Let's say the stream gives 1001 numbers every time.
data = np.arange(1001)/100

norm_data = (data - minimal) * reciprocal_bin_length
norm_data = norm_data.astype(int)
print(norm_data.max())
print(norm_data.min())

5
0

箱子索引应该是0、1、2、3或4,但不是5。

2 回复 | 直到 7 年前

willeM_ Van Onsem 7 年前

穷人的解决方案 “可能是计算 minimum 在您的阵列之间 norm_data 和 nbins-1 :

norm_data = np.minimum(norm_data,nbins-1)

所以所有 5 4 120.0 也将在箱子4)中结束。

hamster on wheels 7 年前

import numpy as np
# Things we know and can pre-calculate
n_bins = 5
minimal = 0.0  
maximal = 10.0
approx = 1.001  # <-- this is new
reciprocal_bin_length = n_bins / (maximal*approx - minimal)

# Let's say the stream gives 1001 numbers every time.
data = np.arange(1001)/100

# can use numexpr for speed.
norm_data = (data - minimal) * reciprocal_bin_length
norm_data = norm_data.astype(int)
print(norm_data.max())
print(norm_data.min())

推荐文章

Redshoe · 当执行numpy.genfromtxt()时,python是否会读取文件的所有行?

4 月前

unfolx · numpy数组不等式的执行时间

5 月前

mchaudh4 · 用numpy表示三对角矩阵

5 月前

Geremia · 2D NumPy数组+1D数组?

5 月前

LMC · Numpy数组布尔索引以获取包含元素

5 月前

HJA24 · 根据条件用值正向填充Numpy矩阵/掩码

6 月前

alpelito7 · 如何只获取numpy数组中每个递增值的第一次出现?

6 月前

Amarth Gûl · 找到一组向量的最近收敛点

6 月前

Mr. W · numpy.divide是函数、类还是其他什么?

6 月前

Mr. W · 为什么numpy.array在编辑内部数据时如此缓慢?

6 月前