function out=flyBurHandData(data,numFlies,roi) %flyBurHandData(data,roi) for i=1:numFlies inx=data(:,2*i+1); iny=data(:,2*i+2); width=roi; out={}; out.x=inx; out.y=iny; out.r=sqrt((inx-width/2).^2+(iny-width/2).^2); out.theta=atan2(iny-width/2,inx-width/2); out.direction=zeros(size(inx,1),1); out.speed=zeros(size(inx,1),1); out.turning=zeros(size(inx,1),1); for i=1:size(inx,1) if mod(i,1000) ==1 i end if i > 1 out.direction(i)=atan2(iny(i)-iny(i-1),inx(i)-inx(i-1)); out.speed(i)=sqrt((iny(i)-iny(i-1))^2+(inx(i)-inx(i-1))^2); if i>2 out.turning(i)=out.direction(i)-out.direction(i-1); if out.turning(i) >pi() out.turning(i)=out.turning(i)-2*pi(); end if out.turning(i) < -pi() out.turning(i)=out.turning(i)+2*pi(); end end end end end out.speed(out.speed>12)=NaN;