def pct_updig(checked, kbouts=2, maxsep=50):
ad = submit_summarize_runs.analysis_dir_from_donebase(
checked.rsplit('.', 1)[0])
#config = eval(open(checked.split('l1800')[0]+'config.dict').read())
newact_tarf = glob(ad + '/newact_ols_postsummary-*.tar')[0]
newact_tarh = tarfile.open(newact_tarf)
newacts = Util.tar2obj_all(newact_tarh)
SHAPE = eval(open(ad + '/SHAPE').read())
cm = vidtools.calc_coordMat(SHAPE)
centroids = [
vidtools.centroid(reduce(lambda x, y: x + y, l), SHAPE, cm)
for l in newacts if len(l) > 0
]
#pdy = config['predug'][1][1]
#return len(filter(lambda x:x<pdy, \
# [y for x,y in centroids]))/float(len(centroids))
k = trial_from_filename(glob(ad + '/*hwc4*hmcdNone*dig_areas.list')[0])
y_offsets = [
p1[1] - p2[1]
for p1, p2 in zip(centroids[:-kbouts], centroids[kbouts:])
if vidtools.hypotenuse(p1, p2) < maxsep
]
return k, len(filter(lambda x: x > 0, y_offsets)) / float(len(y_offsets))
def nearest_point_on_line(lpt1,lpt2,qpt):
h = vidtools.hypotenuse(lpt1,lpt2)
def calc_chord_stats(tun,cm_factor):
h = vidtools.hypotenuse(tun[0],tun[-1])
def calc_chord_stats(tun):
h = vidtools.hypotenuse(tun[0],tun[-1])
ar = numpy.asarray(Image.open(im).convert(mode='L'))
from PIL import Image
ar = numpy.asarray(Image.open(im).convert(mode='L'))
matshow(ar)
fig = plt.figure(1)
def onclick(event):
print 'button=%d, x=%d, y=%d, xdata=%f, ydata=%f'%(
event.button, event.x, event.y, event.xdata, event.ydata)
cid = fig.canvas.mpl_connect('button_press_event', onclick)
fig.canvas.mpl_disconnect(cid)
pylab.ginput(2,0)
p1,p2 = _
from video_analysis import vidtools
from video_analysis import vidtools, viz_vidtools
vidtools.hypotenuse(p1,p2)
p0,p1,p10 = ginput(3,0
)
cm_dists = {}
cm_dists['bot1'] = vidtools.hypotenuse(p0,p1)/1
cm_dists['bot10'] = vidtools.hypotenuse(p0,p10)/1
cm_dists
p0,p1,p10 = ginput(3,0
)
cm_dists['rt1'] = vidtools.hypotenuse(p0,p1)/1
cm_dists['rt10'] = vidtools.hypotenuse(p0,p10)/1
cm_dists
mean([cm_dists['bot1'],cm_dists['bot1'],cm_dists['bot1'],)
cm_dists['bot10'] = vidtools.hypotenuse(p0,p10)/10
cm_dists['rt10'] = vidtools.hypotenuse(p0,p10)/10
cm_dists
def pct_updig(checked,kbouts=2,maxsep=50):
ad = submit_summarize_runs.analysis_dir_from_donebase(checked.rsplit('.',1)[0])
#config = eval(open(checked.split('l1800')[0]+'config.dict').read())
newact_tarf = glob(ad+'/newact_ols_postsummary-*.tar')[0]
newact_tarh = tarfile.open(newact_tarf)
newacts = Util.tar2obj_all(newact_tarh)
SHAPE = eval(open(ad+'/SHAPE').read())
cm = vidtools.calc_coordMat(SHAPE)
centroids = [vidtools.centroid(reduce(lambda x,y:x+y, l),SHAPE,cm) for l in newacts if len(l)>0]
#pdy = config['predug'][1][1]
#return len(filter(lambda x:x<pdy, \
# [y for x,y in centroids]))/float(len(centroids))
k = trial_from_filename(glob(ad+'/*hwc4*hmcdNone*dig_areas.list')[0])
y_offsets = [p1[1]-p2[1] for p1,p2 in zip(centroids[:-kbouts],centroids[kbouts:]) if vidtools.hypotenuse(p1,p2)<maxsep]
return k,len(filter(lambda x:x>0,y_offsets))/float(len(y_offsets))
