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() '''