python - 使用 numpy.genfromtxt 填充缺失值

Question

尽管前面的问题给出了建议：

-9999 作为 numpy.genfromtxt() 的缺失值

使用genfromtxt导入numpy中缺失值的csv数据

我仍然无法处理以缺失值结尾的文本文件，

一个.txt：

1 2 3
4 5 6
7 8

我尝试了多种选项的安排，missing_values但filling_values无法使其正常工作：

import numpy as np

sol = np.genfromtxt("a.txt", 
                    dtype=float,
                    invalid_raise=False, 
                    missing_values=None,
                    usemask=True,
                    filling_values=0.0)
print sol

我想得到的是：

[[1.0 2.0 3.0]
 [4.0 5.0 6.0]
 [7.0 8.0 0.0]]

但相反，我得到：

/usr/local/lib/python2.7/dist-packages/numpy/lib/npyio.py:1641: ConversionWarning: Some errors were detected !
    Line #3 (got 2 columns instead of 3)
  warnings.warn(errmsg, ConversionWarning)
[[1.0 2.0 3.0]
 [4.0 5.0 6.0]]

Answer 1

使用熊猫：

import pandas as pd

df = pd.read_table('data', sep='\s+', header=None)
df.fillna(0, inplace=True)
print(df)
#    0  1  2
# 0  1  2  3
# 1  4  5  6
# 2  7  8  0

pandas.read_table用 s 替换缺失的数据NaN。您可以使用 将这些NaNs 替换为其他值df.fillna。

df是一个pandas.DataFrame。您可以使用以下命令访问底层 NumPy 数组df.values：

print(df.values)
# [[ 1.  2.  3.]
#  [ 4.  5.  6.]
#  [ 7.  8.  0.]]

Answer 2

问题是 numpy 不喜欢参差不齐的数组。由于文件最后一行的第三个位置没有字符，所以 genfromtxt 甚至都不知道它是要解析的东西，更不用说如何处理它了。如果缺失值有一个填充符（任何填充符），例如：

1 2 3
4 5 6
7 8 ''

然后你就可以：

sol = np.genfromtxt("a.txt",
                dtype=float,
                invalid_raise=False,
                missing_values='',
                usemask=False,
                filling_values=0.0)

和：溶胶

array([[  1.,   2.,   3.],
       [  4.,   5.,   6.],
       [  7.,   8.,  nan]])

不幸的是，如果不能选择统一文件的列，您可能会被逐行解析所困。

另一种可能性是如果所有“短”行都在末尾......在这种情况下，您可以利用“usecols”标志来解析所有统一的列，然后使用skip_footer标志来做同样的事情对于剩余的列，同时跳过不可用的列：

sol = np.genfromtxt("a.txt",
                dtype=float,
                invalid_raise=False,
                usemask=False,
                filling_values=0.0,
                usecols=(0,1))
sol
array([[ 1.,  2.],
   [ 4.,  5.],
   [ 7.,  8.]])

sol2 = np.genfromtxt("a.txt",
                dtype=float,
                invalid_raise=False,
                usemask=False,
                filling_values=0.0,
                usecols=(2,),
                skip_footer=1)
sol2
array([ 3.,  6.])

然后从那里组合数组添加填充值：

sol2=np.append(sol2, 0.0)
sol2=sol2.reshape(3,1)
sol=np.hstack([sol,sol2])
sol
array([[ 1.,  2.,  3.],
   [ 4.,  5.,  6.],
   [ 7.,  8.,  0.]])

Answer 3

根据我的经验，最好的方法是手动解析，这个功能对我有用，它可能很慢，但通常足够快。

def manual_parsing(filename,delim,dtype):
    out = list()
    lengths = list()
    with open(filename,'r') as ins:
        for line in ins:
            l = line.split(delim)
            out.append(l)
            lengths.append(len(l))
    lim = np.max(lengths)
    for l in out:
        while len(l)<lim:
            l.append("nan")
    return np.array(out,dtype=dtype)