def runfromfft_graph(f, psd): peaklist = aa.peak_map(psd) peaksdex = aa.peak_assign(peaklist) output = guesswork(f, psd, peaklist, peaksdex) ax = gf.init_image(ylabel='POWER') gf.std_graph(ax, output[:, 0], output[:, 1]) ax2 = gf.get_twinx(ax) gf.std_graph(ax2, output[:, 0], output[:, 2], c='r') plt.show() print(aa.print_identify(aa.identify(output)))
def splinters(fname, numparts): diskdata = np.load(fname) length = len(diskdata) - numparts length = int(length / numparts) i = 0 datas = [] while (i < numparts): newdata = diskdata[i * length:(i + 1) * length] datas.append(newdata) i += 1 result = [] for each in datas: harmonic = harmonics.run(each) identify = aa.identify(harmonic) result.append(identify) return result
def analyze(recording, all_output): print('\nAnalyzing.. ') output = harmonics.run(recording) aa.print_identify(aa.identify(output)) all_output.append(output)
# data = np.load('%s%d.npy' % (fname, index)) # print('%dTH ' % index, end=' ') # # fftfreq.run('%s%d.npy' % (fname, index)) # if index == 50: # break # intervaldata = harmonics.run(data) # intervaldata = aa.identify(intervaldata) # print(intervaldata) # result.append(intervaldata) # index += 1 gen = aa.split_seconds(fname) result = [] for each in gen: intervaldata = harmonics.run(each) intervaldata = aa.identify(intervaldata) result.append(intervaldata) result = np.asarray(result) print(result) def plotdots(ax, i, detect): colored = 'r' if detect: colored = 'b' gf.std_graph(ax, [i, i + 1], [0, 0], c=colored, lw=3) print(np.average(result[:, 2]) - 120) print(np.average(result[:, 3]))
def runfromfft(f, psd): peaklist = aa.peak_map(psd) peaksdex = aa.peak_assign(peaklist) output = guesswork(f, psd, peaklist, peaksdex) print(aa.print_identify(aa.identify(output))) return aa.identify(output)