def process_images():
xVmin = -0.0805
yVmin = -0.0730
xVmax = -0.0399
yVmax = -0.0317
x_factor = (xVmax-xVmin)/120
y_factor = (yVmax-yVmin)/120
print(x_factor)
print(y_factor)
path = 'Z:\\Lab\\Cantilever\\Measurements\\20150807_PushingDiamonds\\Run4\\*.csv'
files = glob.glob(path)
files.sort()
im_list = list()
for f in files:
im_list.append(pd.read_csv(f, header=-1))
print('finished reading')
cur_x = 0
cur_y = 0
n_list = []
d_list = []
offset_x = 0
offset_y = 0
image1 = np.array(im_list[0])
for n in xrange(0, len(im_list) - 1):
if n % 2 == 0:
offset_x = cur_x
offset_y = cur_y
image1 = np.array(im_list[n])
print(n)
print(offset_x)
n_list.append(n)
image2 = np.array(im_list[n+1])
x, y, corr, old_im = track.corr_NVs_no_subset(image1, image2)
#print(x)
#print(y)
#plt.imshow(corr)
#plt.show()
cur_x = x * x_factor + offset_x
cur_y = y * y_factor + offset_y
d_list.append(np.sqrt(cur_x**2 + cur_y**2))
dirpath = 'Z:\\Lab\\Cantilever\\Measurements\\20150807_PushingDiamonds\\Run4'
tag = 'processed'
rw.save_data([n_list, d_list], dirpath, tag)
plt.plot(n_list, d_list)
plt.show()
def process_images():
xVmin = -0.0805
yVmin = -0.0730
xVmax = -0.0399
yVmax = -0.0317
x_factor = (xVmax-xVmin)/120
y_factor = (yVmax-yVmin)/120
print(x_factor)
print(y_factor)
path = 'Z:\\Lab\\Cantilever\\Measurements\\20150807_PushingDiamonds\\Run4\\*.csv'
files = glob.glob(path)
files.sort()
im_list = list()
for f in files:
im_list.append(pd.read_csv(f, header=-1))
print('finished reading')
cur_x = 0
cur_y = 0
n_list = []
d_list = []
offset_x = 0
offset_y = 0
image1 = np.array(im_list[0])
for n in xrange(0, len(im_list) - 1):
if n % 2 == 0:
offset_x = cur_x
offset_y = cur_y
image1 = np.array(im_list[n])
print(n)
print(offset_x)
n_list.append(n)
image2 = np.array(im_list[n+1])
x, y, corr, old_im = track.corr_NVs_no_subset(image1, image2)
#print(x)
#print(y)
#plt.imshow(corr)
#plt.show()
cur_x = x * x_factor + offset_x
cur_y = y * y_factor + offset_y
d_list.append(np.sqrt(cur_x**2 + cur_y**2))
dirpath = 'Z:\\Lab\\Cantilever\\Measurements\\20150807_PushingDiamonds\\Run4'
tag = 'processed'
rw.save_data([n_list, d_list], dirpath, tag)
plt.plot(n_list, d_list)
plt.show()
recording = True
i = 0
while recording:
fo = find_peak(zi, **find_peak_settings)
max_frequencies.append([time.strftime("%Y-%m-%d_%H-%M-%S"), fo])
find_peak_settings.update({
'f_min': fo - f_range / 2,
'f_max': fo + f_range / 2
})
i += 1
if i > 10:
recording = False
print(i, fo)
# save data
rw.save_data(max_frequencies, dirpath, tag + '_frequencies')
rw.save_json(zi_settings,
dirpath + rw.date_prefix() + tag + '_zi_settings.json')
rw.save_json(find_peak_settings,
dirpath + rw.date_prefix() + tag + '_peak_settings.json')
plt.show()
## ==================== TEST SCRIPT T FIND PEAK ========================
# dirpath = 'Z:\\Lab\\Cantilever\\Measurements\\20160302_Resonator_2.2\\20160303_Cold_Resonator_2.2\\data\\'
# tag = 'resonator_2.2'
#
#
# zi = ZI.ZIHF2(**zi_settings)
max_frequencies = []
f_range = find_peak_settings['f_max'] - find_peak_settings['f_min']
recording = True
i = 0
while recording:
fo = find_peak(zi, **find_peak_settings)
max_frequencies.append([time.strftime("%Y-%m-%d_%H-%M-%S"), fo])
find_peak_settings.update({'f_min': fo-f_range/2, 'f_max' : fo+f_range/2})
i+=1
if i>10:
recording = False
print(i, fo)
# save data
rw.save_data(max_frequencies, dirpath, tag + '_frequencies')
rw.save_json(zi_settings, dirpath + rw.date_prefix() + tag + '_zi_settings.json')
rw.save_json(find_peak_settings, dirpath + rw.date_prefix() + tag + '_peak_settings.json')
plt.show()
## ==================== TEST SCRIPT T FIND PEAK ========================
# dirpath = 'Z:\\Lab\\Cantilever\\Measurements\\20160302_Resonator_2.2\\20160303_Cold_Resonator_2.2\\data\\'
# tag = 'resonator_2.2'
#
#
# zi = ZI.ZIHF2(**zi_settings)
#