OpticsPractice-value/EXP06_221109.py at main · TriFullKim/OpticsPractice-value · GitHub

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
from math import radians, cos, sin
import pandas as pd

rd_exp0601 = pd.read_csv('RAW_DATA/exp06/VALUE221102EXP0601.csv')
rd_exp0602_1 = pd.read_csv('RAW_DATA/exp06/VALUE221102EXP0602.1.csv')
rd_exp0602_2 = pd.read_csv('RAW_DATA/exp06/VALUE221102EXP0602.2.csv')
rd_exp0603 = pd.read_csv('RAW_DATA/exp06/VALUE221102EXP0603.csv')

rd_exp0601 = rd_exp0601.set_index(keys=rd_exp0601.columns[0])
rd_exp0602_1 = rd_exp0602_1.set_index(keys=rd_exp0602_1.columns[0])
rd_exp0602_2 = rd_exp0602_2.set_index(keys=rd_exp0602_2.columns[0])
rd_exp0603 = rd_exp0603.set_index(keys=rd_exp0603.columns[0])


polarAngle = [i*10 for i in range(19)]
luxReferPoint = 3


# EXP0601 ###################
trim_exp0601 = pd.DataFrame()
trim_exp0601 = rd_exp0601 - luxReferPoint

trim_exp0601.loc[:, 'IntensityMean'] = trim_exp0601.mean(axis=1)
for degree in polarAngle:
    trim_exp0601.loc[degree, 'theoricalValue'] = trim_exp0601.loc[0, 'IntensityMean'] * cos(radians(degree)) ** 2
# EXP0601_result ############
result_exp0601 = trim_exp0601.loc[:, 'IntensityMean':'theoricalValue']
# Analysis ##################
# trim_exp0601.loc[:, 'ERROR'] = abs(trim_exp0601.loc[:, 'theoricalValue'] - trim_exp0601.loc[:, 'IntensityMean'])

# END #######################
# rd_exp0601.to_clipboard()
result_exp0601.to_clipboard()


angle0602_1 = ['0deg', '45deg']
# EXP0602_1 ##################
trim_exp0602_1 = pd.DataFrame()
trim_exp0602_1 = rd_exp0602_1 - luxReferPoint
for degree in angle0602_1:
    trim_exp0602_1.loc[:, f'{degree}MEAN'] = trim_exp0602_1.loc[:,f'{degree}1':f'{degree}2'].mean(axis=1)
# EXP0602_1_result ############
result_exp0602_1 = trim_exp0602_1.loc[:,f'{angle0602_1[0]}MEAN':f'{angle0602_1[len(angle0602_1)-1]}MEAN']
# Analysis ##################

# END #########################
# rd_exp0602_1.to_clipboard()
result_exp0602_1.to_clipboard()


angle0602_2 = ['25deg', '55deg', '70deg', '85deg', '115deg']
# EXP0602_2 ###################
trim_exp0602_2 = pd.DataFrame()
trim_exp0602_2 = rd_exp0602_2 - luxReferPoint
for degree in angle0602_2:
    trim_exp0602_2.loc[:, f'{degree}MEAN'] = trim_exp0602_2.loc[:,f'{degree}1':f'{degree}2'].mean(axis=1)
# EXP0602_2_result ############
result_exp0602_2 = trim_exp0602_2.loc[:,f'{angle0602_2[0]}MEAN':f'{angle0602_2[len(angle0602_2)-1]}MEAN']
# Analysis ##################

# END #########################
# rd_exp0602_2.to_clipboard()
result_exp0602_2.to_clipboard()


# EXP0603 ###################
trim_exp0603 = pd.DataFrame()
trim_exp0603= rd_exp0603 - luxReferPoint
# EXP0603_result ############
result_exp0603 = trim_exp0603
# Analysis ##################
trim_exp0603.loc['MIN','1차'] = trim_exp0603.loc[:,'1차'].idxmin()
trim_exp0603.loc['MIN','2차'] = trim_exp0603.loc[:,'2차'].idxmin()

brewsterAngle = trim_exp0603.loc['MIN',:].mean()
n_glass = sin(radians(brewsterAngle))/cos(radians(brewsterAngle)) # tan(brewsterAngle)
# END #######################
rd_exp0603.to_clipboard()