Question

我的第一篇帖子。

所以我将数据加载到名为 f1_data 的变量中，然后将其传递给 pm.removeDC（）函数进行一些信号处理，并保持结果变成了同一个变量。但是，我想只用原来的 f1_data 替换第8列，我称之为 raw_data ，我无法弄明白为什么它没有＆＃39;工作。这是功能。帮助任何人？

内部文件pm.py

def removeDC(data):    
    # define the filter
    butter_order = 2
    hp_cutoff_Hz = 1.0
    b, a = signal.butter(butter_order, hp_cutoff_Hz/(fs_Hz / 2.0), 'highpass')

    for i in range(1,9):
        data[:,i] = signal.lfilter(b, a, data[:,i], 0)

    return (data)

def get_epoch1(data, t_sec, epoch, f_tup, col):
#f_tup = (f_wdir, f_name, f_columns, out_save, out_dir, out_number, fig_width)
    f_name = f_tup[1]
    fig_width = f_tup[6]

    epoch_boolvector = (t_sec >= epoch[0][0]) & (t_sec <= epoch[0][1])
    epoch_timescale = t_sec[epoch_boolvector]
    epoch_data = data[epoch_boolvector]

    plt.figure(figsize=(fig_width,8), dpi=96)
    plt.plot(epoch_timescale, epoch_data[:,col]);
    plt.xlim(epoch_timescale[0], epoch_timescale[-1])
    plt.show()
    return (epoch_boolvector, epoch_timescale, epoch_data)

在主文件中

#load the whole data
(f1_data, f1_data_indices, f1_timescale) = pm.load_data(f1_wdir, f1_name)

raw_data = f1_data[:] #create copy of f1_data

(f1ep1_boolvector, f1ep1_timescale, f1ep1_data) = pm.get_epoch1(f1_data, f1_timescale, f1_epochs[1], f1_tup, 8)

#--- filter data to remove DC (1Hz)
f1_data = pm.removeDC(f1_data)


# replace only channel 8 with original data
f1_data[:,8] = raw_data[:,8]
(f1ep2_boolvector, f1ep2_timescale, f1ep2_data) = pm.get_epoch1(f1_data, f1_timescale, f1_epochs[1], f1_tup, 8)

Answer 1

解决方案是import copy并使用copy.deepcopy函数。有关详细信息，请查看此链接： docs.python.org/2/library/copy.html

当我raw_data = f1_data[:]后，我pm.removeDC()：

 raw_data is f1_data: False
(raw_data == f1_data).all(): True

但是当我raw_data = copy.deepcopy(f1_data)后，pm.removeDC()：

 raw_data is f1_data: False  
(raw_data == f1_data).all(): False

Python中奇怪的变量行为

1 个答案: