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)