# -*- coding: utf-8 -*-
"""
Created on Fri Dec 04 13:23:33 2015
@author: Kyle Ellefsen
"""
import numpy as np
from plugins.detect_puffs.puff_simulator.puff_simulator import simulate_puffs
from plugins.detect_puffs.threshold_cluster import threshold_cluster
data_window = simulate_puffs(nFrames=1000, nPuffs=20)
data_window = g.m.currentWindow
norm_window = data_window
norm_window.setName('Normalized Window')
binary_window = threshold(1.1, keepSourceWindow=True)
binary_window.setName('Binary Window')
threshold_cluster(binary_window,
data_window,
norm_window,
density_threshold=5.8)
## run through manual steps
## once threshold_cluster has run, see what the percentage of puffs it detected, and the false positive rate
detected_puffs = g.m.puffAnalyzer.puffs.puffs
simulated_puffs = simulate_puffs.puffs
nearest_puffs = []
for puff in detected_puffs:
k = puff.kinetics
detected_puff = np.array([k['t_peak'], k['x'], k['y']])
"""
Created on Fri Feb 06 11:24:36 2015
@author: Kyle Ellefsen
"""
if __name__ == '__main__':
import os, sys
flika_dir = os.path.join(os.path.expanduser('~'), 'Documents', 'GitHub',
'flika')
sys.path.append(flika_dir)
from flika import *
start_flika()
from plugins.detect_puffs.threshold_cluster import threshold_cluster
from plugins.detect_puffs.puff_simulator.puff_simulator import simulate_puffs
simulate_puffs(nFrames=1000, nPuffs=20)
baseline = -5 # This is the mean pixel value in the absence of photons.
subtract(baseline)
data_window = ratio(
0, 30, 'average'
) # ratio(first_frame, nFrames, ratio_type). Now we are in F/F0
data_window.setName('Data Window (F/F0)')
norm_image = data_window.image - 1
norm_window = Window(norm_image)
norm_window.setName('Normalized Window')
blurred_window = gaussian_blur(2, norm_edges=True, keepSourceWindow=True)
blurred_window.setName('Blurred Window')
threshold_cluster(data_window,
blurred_window,
blurred_window,
blur_thresh=.3)
# -*- coding: utf-8 -*-
"""
Created on Fri Dec 04 13:23:33 2015
@author: Kyle Ellefsen
"""
import numpy as np
from plugins.detect_puffs.puff_simulator.puff_simulator import simulate_puffs
from plugins.detect_puffs.threshold_cluster import threshold_cluster
data_window=simulate_puffs(nFrames=1000,nPuffs=20)
data_window=g.m.currentWindow
norm_window=data_window
norm_window.setName('Normalized Window')
binary_window=threshold(1.1, keepSourceWindow=True)
binary_window.setName('Binary Window')
threshold_cluster(binary_window,data_window,norm_window,density_threshold=5.8)
## run through manual steps
## once threshold_cluster has run, see what the percentage of puffs it detected, and the false positive rate
detected_puffs=g.m.puffAnalyzer.puffs.puffs
simulated_puffs=simulate_puffs.puffs
nearest_puffs=[]
for puff in detected_puffs:
k=puff.kinetics
detected_puff=np.array([k['t_peak'],k['x'],k['y']])
difference=np.sqrt(np.sum((simulated_puffs-detected_puff)**2,1))
closest_idx=np.argmin(difference)
# -*- coding: utf-8 -*-
"""
Created on Fri Feb 06 11:24:36 2015
@author: Kyle Ellefsen
"""
if __name__ == '__main__':
import os, sys; flika_dir = os.path.join(os.path.expanduser('~'),'Documents', 'GitHub', 'flika'); sys.path.append(flika_dir); from flika import *; start_flika()
from plugins.detect_puffs.threshold_cluster import threshold_cluster
from plugins.detect_puffs.puff_simulator.puff_simulator import simulate_puffs
simulate_puffs(nFrames=1000,nPuffs=20)
baseline = -5 # This is the mean pixel value in the absence of photons.
subtract(baseline)
data_window=ratio(0, 30, 'average') # ratio(first_frame, nFrames, ratio_type). Now we are in F/F0
data_window.setName('Data Window (F/F0)')
norm_image = data_window.image - 1
norm_window = Window(norm_image)
norm_window.setName('Normalized Window')
blurred_window = gaussian_blur(2, norm_edges=True, keepSourceWindow=True)
blurred_window.setName('Blurred Window')
threshold_cluster(data_window, blurred_window, blurred_window, blur_thresh=.3)
sys.exit(g.app.exec_()) # This is required to run outside of Spyder or PyCharm
