Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- import csv
- import pandas as pd
- import numpy as np
- from pandas import DataFrame
- import pandas.io.data
- import scipy.fftpack
- import matplotlib.pyplot as plt
- #reader = csv.DictReader(open('SPNTSTB.txt','rt'), delimiter = 't')
- Fname = raw_input("Input File Name (Must be modified '.csv' file): ")
- print (Fname)
- Newname = Fname + 'Scrb.txt'
- print ("Scrubbed File in: "), Newname
- reader = pd.read_csv(Fname+'.txt', sep = None, engine = 'python')
- df = pd.DataFrame(reader)# reader["Gyro_X"].astype(float)
- dg = pd.DataFrame(reader)
- Timefld = 0
- XCol = 1
- YCol = 2
- ZCol = 3
- ######################
- # Parsing Frame into individual columns
- ######################
- Time_Col = [0]
- Samp_X_Col = [0]
- Samp2_X_Col = [0]
- Samp3_Y_Col = [0]
- print 'Parsing'
- Time_Col = df.iloc[:,[Timefld]]
- Samp_X_Col = df.iloc[:,[XCol]]
- Samp2_X_Col = df.iloc[:,[YCol]]
- Samp3_Y_Col = df.iloc[:,[ZCol]]
- print 'Parsing Complete'
- xSignal = Samp_X_Col
- # Number of samplepoints
- N = 10000
- # sample spacing
- T = 1.0 / 200.0
- yf = scipy.fftpack.rfft(xSignal)
- xf = np.linspace(0.0, 1.0/(2.0*T), N/2)
- fig, ax = plt.subplots()
- #ax.plot(xf, 2.0/N * np.abs(yf[:N//2]))
- ax.plot(xf, yf[:N//2])
- plt.show()
- -0.659912109
- -0.636962891
- -0.616210938
- -0.598876953
- -0.574462891
- -0.556396484
- -0.532470703
- -0.533935547
- -0.524902344
- -0.511962891
- -0.536621094
- -0.560302734
- -0.558837891
- -0.585205078
- -0.631591797
- -0.627929688
- -0.632324219
- -0.676269531
- -0.677734375
- -0.651123047
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement