Home > analysis > WBCanalysis2.m

WBCanalysis2

PURPOSE ^

SYNOPSIS ^

function WBCanalysis2(remote)

DESCRIPTION ^

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function WBCanalysis2(remote)
0002 
0003 if (remote==0)
0004     iPath = 'C:\Research\DataSets\LILLYwbcrollingapc\';
0005 else
0006     iPath = '/v1/rich/DataSets/LILLYwbcrollingapc/';
0007 end
0008 
0009 cd(iPath);
0010 dataFolders = dir(iPath);
0011 
0012 fid = fopen('AlternativeAnalysis.txt','w');
0013 
0014 
0015 for fo = 3: length(dataFolders)
0016     if (dataFolders(fo).isdir ==1)
0017         dataFolder = dataFolders(fo).name;
0018         
0019         try
0020 
0021             load(strcat(dataFolder,'/','_frames.mat'));
0022             load(strcat(dataFolder,'/','_tracks.mat'));
0023             % No of static cells by detection BG
0024             load(strcat(dataFolder,'/','_nStatic.mat'));
0025 
0026             noFrames = length(frames);
0027             nCells = [];
0028             m_nCells = 0;
0029             v = [];
0030             
0031             for fr = 1: noFrames
0032                 nCells = [nCells; frames(fr).nCells];
0033             end
0034 
0035             m_nCells = mean(nCells);
0036             
0037             
0038             % Velocity
0039             if (~isempty(tracksR))
0040                 nCells = tracksR(end,1);
0041                 for tr = 1:nCells
0042                     track = tracksR((tracksR(:,1)==tr),:);
0043                     if size(track,1) > 3
0044 
0045                         nFrames = size(track,1);
0046                         % Calculate velocity
0047                         v = [];
0048                         for fr = 1:nFrames-1
0049                             vTemp = round(sqrt((track(fr+1,3)-track(fr,3))^2 + (track(fr+1,4)-track(fr,4))^2));
0050                             v = [v; vTemp];
0051                         end
0052                     end
0053                 end
0054                 if (size(v,1)~=0)
0055                     m_v = mean(v);
0056                     sd_v = std(v);
0057                 else
0058                     m_v = 0;
0059                     sd_v =0;
0060                 end
0061 
0062                 csvFile = strcat(dataFolder,'/',dataFolder,'.alt.csv');
0063                 oData = [m_nCells nStatic m_v sd_v];
0064                 csvwrite(csvFile,oData);
0065                 disp(strcat('Collect report in folder:',dataFolder,'...'));
0066                 
0067             else
0068                 csvFile = strcat(dataFolder,'/',dataFolder,'.alt.csv');
0069                 oData = [m_nCells nStatic 0 0];
0070                 csvwrite(csvFile,oData);
0071                 disp(strcat('Collect report in folder:',dataFolder,'...'));
0072             end
0073             fprintf(fid,'%s %7.2f %4.0f %7.2f %7.2f\n',... 
0074                 dataFolder,oData(1),oData(2),oData(3),oData(4));
0075         catch
0076             disp(['Report file in ', dataFolder ,' is not available']);
0077         end
0078         
0079     end
0080 end
0081 
0082 fclose(fid);
0083 
0084 
0085 
0086 
0087 
0088 end

Generated on Thu 17-Mar-2011 14:45:51 by m2html © 2005