Advertisement
Guest User

Untitled

a guest
Oct 25th, 2016
58
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.39 KB | None | 0 0
  1. import csv
  2. import pandas as pd
  3. import numpy as np
  4. from pandas import DataFrame
  5. import pandas.io.data
  6. import scipy.fftpack
  7. import matplotlib.pyplot as plt
  8.  
  9. #reader = csv.DictReader(open('SPNTSTB.txt','rt'), delimiter = 't')
  10. Fname = raw_input("Input File Name (Must be modified '.csv' file): ")
  11. print (Fname)
  12. Newname = Fname + 'Scrb.txt'
  13. print ("Scrubbed File in: "), Newname
  14.  
  15. reader = pd.read_csv(Fname+'.txt', sep = None, engine = 'python')
  16. df = pd.DataFrame(reader)# reader["Gyro_X"].astype(float)
  17. dg = pd.DataFrame(reader)
  18.  
  19.  
  20. Timefld = 0
  21. XCol = 1
  22. YCol = 2
  23. ZCol = 3
  24.  
  25.  
  26. ######################
  27. # Parsing Frame into individual columns
  28. ######################
  29.  
  30. Time_Col = [0]
  31. Samp_X_Col = [0]
  32. Samp2_X_Col = [0]
  33. Samp3_Y_Col = [0]
  34.  
  35. print 'Parsing'
  36. Time_Col = df.iloc[:,[Timefld]]
  37. Samp_X_Col = df.iloc[:,[XCol]]
  38. Samp2_X_Col = df.iloc[:,[YCol]]
  39. Samp3_Y_Col = df.iloc[:,[ZCol]]
  40.  
  41.  
  42. print 'Parsing Complete'
  43.  
  44.  
  45. xSignal = Samp_X_Col
  46.  
  47. # Number of samplepoints
  48. N = 10000
  49. # sample spacing
  50. T = 1.0 / 200.0
  51. yf = scipy.fftpack.rfft(xSignal)
  52. xf = np.linspace(0.0, 1.0/(2.0*T), N/2)
  53.  
  54. fig, ax = plt.subplots()
  55. #ax.plot(xf, 2.0/N * np.abs(yf[:N//2]))
  56. ax.plot(xf, yf[:N//2])
  57. plt.show()
  58.  
  59. -0.659912109
  60. -0.636962891
  61. -0.616210938
  62. -0.598876953
  63. -0.574462891
  64. -0.556396484
  65. -0.532470703
  66. -0.533935547
  67. -0.524902344
  68. -0.511962891
  69. -0.536621094
  70. -0.560302734
  71. -0.558837891
  72. -0.585205078
  73. -0.631591797
  74. -0.627929688
  75. -0.632324219
  76. -0.676269531
  77. -0.677734375
  78. -0.651123047
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement