-
Notifications
You must be signed in to change notification settings - Fork 0
/
Test.py
119 lines (106 loc) · 3.28 KB
/
Test.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
# To add a new cell, type '#%%'
# To add a new markdown cell, type '#%% [markdown]'
#%% [markdown]
# # Normalización de longitudes de onda
#%% [markdown]
# Se utiliza pandas y numpy para crear el Dataframe y los arreglos necesarios, pandas permite seleccionar solo los datos dentro del rango de longitudes de ondas deseadas. Utilizando la siguiente ecuacion es posible normalizar los datos de potencia entre [0-1]
#
# \begin{equation*}
# X' = \frac{X-X_{min}}{X_{max}-X_{min}}
# \end{equation*}
#
#%%
#Importar librerias
import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
import scipy
from scipy import signal
#%%
#Cargar datos del documento
excelfile=pd.read_csv('Files/6602nm.csv')
#Delimitar datos a utilizar por longitud de onda
excelfile=excelfile[(excelfile['Longitud'] > 500) & (excelfile['Longitud'] < 800)]
#Creando arreglos
wave = excelfile['Longitud']
X = excelfile['Amplitud']
#Suevizado de datos
X = signal.savgol_filter(X,53,5)
Xmax = np.amax(X)
Xmin = np.amin(X)
Amplitud = (X-Xmin)/(Xmax-Xmin)
plt.figure(figsize=(10,10))
plt.plot(wave,Amplitud)
plt.minorticks_on()
plt.grid(which='major',linestyle='-')
plt.grid(which='minor',linestyle=':')
plt.xlim(560,760)
plt.ylim(0,1)
plt.xlabel('Longitud de onda', fontsize=18)
plt.ylabel('Amplitud', fontsize=18)
plt.title('Espectro normalizado 660nm', fontsize=18)
plt.savefig('Graphs/6603norm.png')
#%%
excelfile=pd.read_csv('Files/7802nm.csv')
excelfile=excelfile[(excelfile['Longitud'] > 600) & (excelfile['Longitud'] < 900)]
wave = excelfile['Longitud']
X = excelfile['Amplitud']
#Suevizado de datos
X = signal.savgol_filter(X,53,5)
Xmax = np.amax(X)
Xmin = np.amin(X)
Amplitud = (X-Xmin)/(Xmax-Xmin)
plt.figure(figsize=(10,10))
plt.plot(wave,Amplitud)
plt.minorticks_on()
plt.grid(which='major',linestyle='-')
plt.grid(which='minor',linestyle=':')
plt.xlim(680,880)
plt.ylim(0,1.1)
plt.xlabel('Longitud de onda', fontsize=18)
plt.ylabel('Amplitud', fontsize=18)
plt.title('Espectro normalizado 780nm', fontsize=18)
plt.savefig('Graphs/780norm.png')
#%%
excelfile=pd.read_csv('Files/8082nm.csv');
excelfile=excelfile[(excelfile['Longitud'] > 700) & (excelfile['Longitud'] < 950)]
wave = excelfile['Longitud']
X = excelfile['Amplitud']
#Suevizado de datos
X = signal.savgol_filter(X,53,5)
Xmax = np.amax(X)
Xmin = np.amin(X)
Amplitud = (X-Xmin)/(Xmax-Xmin)
plt.figure(figsize=(10,10))
plt.plot(wave,Amplitud)
plt.minorticks_on()
plt.grid(which='major',linestyle='-')
plt.grid(which='minor',linestyle=':')
plt.xlabel('Longitud de onda', fontsize=18)
plt.xlim(708,908)
plt.ylim(0,1.1)
plt.ylabel('Amplitud', fontsize=18)
plt.title('Espectro normalizado 808nm', fontsize=18)
plt.savefig('Graphs/8082norm.png')
#%%
excelfile=pd.read_csv('Files/8502nm.csv')
excelfile=excelfile[(excelfile['Longitud'] > 700) & (excelfile['Longitud'] < 950)]
wave = excelfile['Longitud']
X = excelfile['Amplitud']
#Suevizado de datos
X = signal.savgol_filter(X,53,5)
Xmax = np.amax(X)
Xmin = np.amin(X)
Amplitud = (X-Xmin)/(Xmax-Xmin)
plt.figure(figsize=(10,10))
plt.plot(wave,Amplitud)
plt.minorticks_on()
plt.grid(which='major',linestyle='-')
plt.grid(which='minor',linestyle=':')
plt.xlabel('Longitud de onda', fontsize=18)
plt.xlim(740,940)
plt.ylim(0,1.1)
plt.ylabel('Amplitud', fontsize=18)
plt.title('Espectro normalizado 850nm', fontsize=18)
plt.savefig('Graphs/850norm.png')
#%%