将点列表转换为numpy 2D数组

时间:2011-05-18 15:30:25

标签: python numpy

我正在使用genfromtxt导入一个二维数组,该数组的所有值都列在表单的文本文件中(x和y是整数):

    x1   y1   z1
    x2   y2   z2
    :    :    :

我正在使用下面的for循环,但我很确定必须采用一行方式来完成它。什么是更有效的方式来进行这种转换?

raw = genfromtxt(file,skip_header = 6)

xrange = ( raw[:,0].min() , raw[:,0].max() )
yrange = ( raw[:,1].min() , raw[:,1].max() )

Z = zeros(( xrange[1] - xrange[0] +1 , yrange[1] - yrange[0] +1 ))

for row in raw:
    Z[ row[0]-xrange[0] , row[1]-yrange[0] ] = row[2]

4 个答案:

答案 0 :(得分:3)

您可以使用以下内容替换for循环:

xidx = (raw[:,0]-xrange[0]).astype(int)
yidx = (raw[:,1]-yrange[0]).astype(int)

Z[xidx, yidx] = raw[:,2]

答案 1 :(得分:0)

要从文件导入矩阵,您只需拆分行然后转换为int。

[[int(i) for i in j.split()] for j in open('myfile').readlines()]

当然,我假设你的文件只包含矩阵。

最后,您可以将此二维数组转换为numpy。

答案 2 :(得分:0)

您可以尝试这样的事情:

>>> Z = zeros((3, 3))
>>> test = array([[0, 1, 2], [1, 1, 6], [2, 0, 4]])
>>> Z[test[:, 0:2].T.tolist()]
array([ 0.,  0.,  0.])
>>> Z[test[:, 0:2].T.tolist()] = test[:, 2]
>>> Z
array([[ 0.,  2.,  0.],
       [ 0.,  6.,  0.],
       [ 4.,  0.,  0.]])

在你的情况下:

Z[(raw[:, 0:2] - minimum(raw[:, 0:2], axis=0)).T.tolist()] = raw[:, 2]

答案 3 :(得分:0)

您也可以使用numpy.searchsorted,这也允许非等间距/浮动数据:

raw = genfromtxt(file,skip_header = 6)

xvalues = numpy.sorted(set(raw[:,0]))
xidx = numpy.searchsorted(xvalues, raw[:,0])

yvalues = numpy.sorted(set(raw[:,1]))
yidx = numpy.searchsorted(yvalues, raw[:,1])

Z = numpy.zeros((len(xvalues), len(yvalues)))
Z[xidx, yidx] = raw[:,2]

否则,我会关注Simon的回答。