SciPy Medfilt仅返回0值

时间:2016-10-04 18:53:30

标签: python numpy scipy

您好我正在编写一个代码,用于计算数据集中给定窗口大小的中值。我正在使用SciPy的medfilt。我不明白为什么返回的中值数组全为零。我已经改变了内核大小,但这并没有影响任何东西,我想知道我的数组的形状是否会影响medfilt。这是我的代码:

    from __future__ import division
    import numpy as np
    import matplotlib.pyplot as plt
    import scipy as sp
    from scipy import signal

    filename = "hbond.txt"
    hbond_val = []
    with open(filename,'r') as f:
        for line in f:
            f_line = np.array(line.split())
            hbond_val.append(f_line)
    bond_array = np.asarray(hbond_val)
    bond_array_float = bond_array.astype(float)
    bond_1d = np.reshape(bond_array_float,50001)
    #print bond_1d.shape
    median = sp.signal.medfilt(bond_array_float,101)
    #plt.plot(range(len(bond_array)),bond_array,'b')
    #plt.plot(range(len(median)),median,'r')
    #plt.show()
    print median #median returns array full of zeros

0 个答案:

没有答案