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
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
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
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