def theta_labeling(raw):
Th_label = []
for part in range(len(All)):
start = raw[part][5]-1
stop = len(raw[part][0])-start
database = analysis.make_twoD(raw[part],start,stop,128)
Th = []
for i in range(len(database)):
RI = analysis.R_I(raw[part],i)
R = RI[0]
I = RI[1]
if analysis.dots(database[i]) >= 85:
error = analysis.for_error(R,I)
if error >= 0.85 and error <= 1:
theta = analysis.ura(R,I)
else:
theta = analysis.circle_fit(R,I)[1]
Th.append(theta)
else:
Th.append(0)
for i in Th:
if i >= raw[part][6]:
Th_label.append(1)
else:
Th_label.append(0)
return Th_label
def theta_labeling(data):
database = analysis.make_twoD(data, start, stop, 128, 80)
Th = []
for i in range(len(database)):
RI = analysis.R_I(data, i)
R = RI[0]
I = RI[1]
if analysis.dots(database[i]) >= 85:
error = analysis.for_error(R, I)
if error >= 0.85 and error <= 1:
print(i, error)
theta = analysis.ura(R, I)
else:
theta = analysis.circle_fit(R, I)[1]
Th.append(theta)
else:
Th.append(0)
angle = analysis.scatter(Th, 10)
plt.plot(angle[0], angle[1])
plt.show()
Th_label = []
for i in Th:
if i >= 220:
Th_label.append(1)
else:
Th_label.append(0)
return Th_label
position = []
print("raw",len(raw))
#package data
data = [freq, Amp, phase]
data1 = [freq1,Amp1,phase1]
data2 = [freq2,Amp2,phase2]
data3 = [freq3,Amp3,phase3]
data4 = [freq4,Amp4,phase4]
All = []
n = 28
start = 24
stop = len(freq)-24
database = analysis.make_twoD(data,start,stop,n)
All.append(database)
start = int(0.025/((max(freq1)-min(freq1))/len(freq1)))
stop = len(freq1)-start
database1 = analysis.make_twoD(data1,start,stop,n,start*2)
All.append(database1)
'''
start = 24
stop = len(freq2)-24
database2 = analysis.make_twoD(data2,start,stop,n)
All.append(database2)
start = 24
stop = len(freq3)-24
database3 = analysis.make_twoD(data3,start,stop,n)
) #data package data = [freq, Amp, phase] data1 = [freq1, Amp1, phase1] data2 = [freq2, Amp2, phase2] data3 = [freq3, Amp3, phase3] data4 = [freq4, Amp4, phase4] subede = [data, data1, data2, data3, data4] #make data in 2D all_database = [] n = 28 start = 24 stop = len(freq) - 24 database = analysis.make_twoD(data, start, stop, n) all_database.append([database, freq, start]) start = int(0.025 / ((max(freq1) - min(freq1)) / len(freq1))) stop = len(freq1) - start database1 = analysis.make_twoD(data1, start, stop, n, start * 2) all_database.append([database1, freq1, start]) start = 24 stop = len(freq2) - 24 database2 = analysis.make_twoD(data2, start, stop, n) all_database.append([database2, freq2, start]) start = 24 stop = len(freq3) - 24 database3 = analysis.make_twoD(data3, start, stop, n)
Amp4 = np.loadtxt(
"./txt_data/S21(B)@fp=4.86-4.92GHz Pp=-5dBm B=1mA IFB=100Hz 601pts/Amp.txt"
)
phase4 = np.loadtxt(
"./txt_data/S21(B)@fp=4.86-4.92GHz Pp=-5dBm B=1mA IFB=100Hz 601pts/phase.txt"
)
data = [freq2, Amp2, phase2]
#print(analysis.fre_point(freq3,750))
RI = analysis.R_I(data, 2500)
plt.plot(RI[0], RI[1])
plt.show()
mat = analysis.make_twoD(data, 0, 4, 28, 50) #data,start,stop,n
for i in mat:
plt.imshow(i, cmap="gray_r")
plt.show()
'''
data = [freq,Amp,phase]
Range = 50#50*0.001/((max(freq)-min(freq))/len(freq))
start = int(Range)
stop = len(freq) - int(Range)
mat = analysis.make_twoD(data,3406,3407,128,Range)#data,start,stop,n
plt.imshow(mat[0],cmap = "gray_r")
plt.show()
'''