-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathgrabTL.m
52 lines (47 loc) · 1.3 KB
/
grabTL.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
function [tres, TresL]=grabT(Drp1,StartFrame, EndFrame, Timescale,FissionOrReversal, MinDiameter);
% allows the user to select where Tf is for fission events and Tres is for
% reversals
%
% INPUTS:
% Drp1 = Drp1 signal over time
% StartFrame = first frame
% EndFrame = last frame
% Timescale = time between frames in seconds
% FissionOrReversal = 1 if fission, 0 if reversal
%
% OUTPUTS:
% Tres = Tf for fission or [start finish] Tres for reversals
%
figure
yyaxis left
% Drp1
plot(Drp1,'Linewidth',2); hold on; plot(Drp1, 'd');
xlabel('Time before fission [s]')
ylabel('Normalized Drp1 intensity')
title('Residency time')
yyaxis right
% MinDiameter
plot(MinDiameter,'Linewidth',2)
ylabel('MinDiameter')
% if FissionOrReversal ==1
% [tres,~]=ginput(1);
% % Times=(StartFrame-EndFrame:0)*Timescale; Times=Times';
% Times=StartFrame:1:EndFrame;Times=Times';
% ID=knnsearch(Times,tres);
% tres=Times(ID);
% TresL=EndFrame-tres;
% else
% [tres,~]=ginput(2);
% % Times=(StartFrame-EndFrame:0)*Timescale; Times=Times';
% Times=StartFrame:1:EndFrame;Times=Times';
%
% ID=knnsearch(Times,tres);
% tres=Times(ID);
% TresL=tres(2)-tres(1);
%
% end
pause;
TresL=1;%TresL*Timescale;
tres=1'%tres*Timescale;
close;
end