Changeset 804 for trunk/src/series/float_tracking.m

Timestamp:

Jul 30, 2014, 8:38:17 PM (10 years ago)

Author:

sommeria

Message:

float tracking updated

File:

• trunk/src/series/float_tracking.m (modified) (5 diffs)

trunk/src/series/float_tracking.m

@@ -78 +78 @@
     end
     % parameters specific to the function 'particle_tracking'
-%     Par.Nblock=[];%size of image subblocks for background determination, =[]: no sublock
-%     Par.ThreshLum=-2000;% luminosity threshold for particle detection, < 0 for black particles, >0 for white particles
-%   ParamOut.ActionInput=Par;
+        Par.Nblock=10;%size of image subblocks for background determination, =[]: no sublock
+    Par.ThreshLum=210;% luminosity threshold for particle detection, < 0 for black particles, >0 for white particles
+  ParamOut.ActionInput=Par;
     return
 end
@@ -92 +92 @@
     checkrun=0;
 end
+hseries=findobj(allchild(0),'Tag','series');
+RUNHandle=findobj(hseries,'Tag','RUN');%handle of RUN button in GUI series
+WaitbarHandle=findobj(hseries,'Tag','Waitbar');%handle of waitbar in GUI series
 
 %% define the directory for result file
@@ -138 +141 @@
 %filter for particle center of mass(luminosity)
 %Nblock=Param.ActionInput.Nblock;
-%ThreshLum=Param.ActionInput.ThreshLum;% luminosity threshold for particle detection, < 0 for black particles, >0 for white particles
-AbsThreshold=30; %threshold below which a pixel is considered belonging to a float
+ThreshLum=Param.ActionInput.ThreshLum;% luminosity threshold for particle detection, < 0 for black particles, >0 for white particles
+%AbsThreshold=30; %threshold below which a pixel is considered belonging to a float
 SizePart=4;
 %
@@ -164 +167 @@
 for ifile=1:nbfield
     if checkrun
-        if strcmp(get(Param.RUNHandle,'BusyAction'),'queue')
-            update_waitbar(Param.WaitbarHandle,ifile/nbfield)
-        else
-            break% leave the loop if the STOP button is activated on the GUI series
+                update_waitbar(WaitbarHandle,ifile/nbfield)
+        if ~isempty(RUNHandle) &&ishandle(RUNHandle) && ~strcmp(get(RUNHandle,'BusyAction'),'queue')
+            disp('program stopped by user')
+            return
         end
     end
@@ -186 +189 @@
     Mask(:,end-SizePart:end)=0;  
     if ifile ==1
-        [Js,Is]=find(A<AbsThreshold & Mask==1);%indices (I,J) of dark pixels
+        if ThreshLum>0 %brigth particles
+        [Js,Is]=find(A>ThreshLum & Mask==1);%indices (I,J) of dark pixels
+        else %black particle
+               [Js,Is]=find(A<ThreshLum & Mask==1);%indices (I,J) of dark pixels
+        end
     else
         Is=round(Xtime(ifile-1,:));