【发布时间】:2020-08-19 09:18:55
【问题描述】:
import glob
import pandas as pd
import seaborn as sns
import numpy as np
from scipy import signal
import matplotlib.pyplot as plt
files = glob.glob("Angular_position_*_*.csv")
output = pd.DataFrame()
for f in files:
df = pd.read_csv(f)
time = df.iloc[:,0]
time = time.to_numpy()
ynew = df.iloc[:,1:]
ynew = ynew.to_numpy()
lowPassCutoffFreq = 6.0 # Cut off frequency
Sample_freq = 150; #Target sample frequency
N = 2 # Order of the filter; In this case 2nd order
Wn = lowPassCutoffFreq/(Sample_freq/2) #Normalize frequency
b, a = signal.butter(5, Wn, btype='low',analog=False,output='ba')
#scipy.signal.butter(N, Wn, btype='low', analog=False, output='ba', fs=None)
output = signal.filtfilt(b, a, ynew, axis=0)
np.savetxt("enter directory path/Filtered_files/Filtered_Angular_position_*_*", output, delimiter = ', ', newline = "\n")
我正在尝试读取目录中的所有文件,然后对它们进行低通滤波。之后,结果将一个接一个地保存,但不会保存在一个文件中。结果为每个文件提供了 3 列,理想情况下,我希望它们用标题命名,例如col1,col2,col3。
不使用 glob,我可以单独过滤所有文件,但我有 100 多个这样的文件。
任何帮助将不胜感激。
祝你好运,
【问题讨论】:
标签: python-3.x pandas numpy glob lowpass-filter