Home > analysis > get_GTsize.m

get_GTsize

PURPOSE ^

SYNOPSIS ^

function [totalSize] = get_GTsize(dataFolder)

DESCRIPTION ^

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function [totalSize] = get_GTsize(dataFolder)
0002 
0003  iPath = 'C:\Research\DataSets\';
0004  iFolder = strcat(iPath,dataFolder);
0005  
0006  
0007  cd (iFolder);
0008  load _tracksGT;
0009 
0010 
0011  if (strcmp(dataFolder,'clpsaline1bl#4')==1)
0012 %     range = [115 135; 140 160 (Gold); 175 210 (Gold) ; 230 245];
0013     ids = [73 66; 78 79; 90 94];
0014     rangeD = [115 135; 140 160; 175 210];
0015 elseif (strcmp(dataFolder,'clpsaline1bl#3')==1)
0016 %     range = [215 225; 235 245; 260 270];
0017     ids = [13 9; 24 25; 31 36];
0018     rangeD = [215 225; 235 245; 260 270];
0019 elseif (strcmp(dataFolder,'clpsaline1bl#1')==1)
0020 %     range = [2 12;35 45; 50 60 fast flowing cell];
0021     ids = [33 32; 5 10; 15 17];
0022     rangeD = [2 12; 35 45; 50 60];
0023     elseif (strcmp(dataFolder,'rhod74-10-')==1)
0024 %     range = [2 22; 50 85; 140 155];
0025     ids = [43 44; 14 15; 63 67];
0026     rangeD = [2 22; 50 85; 140 155];
0027 else
0028     error('not yet supported dataset.');
0029 end
0030 
0031  
0032 totalSize = [];
0033 % for id  = 1: size(ids,2)
0034 %     tr1 = tracksGT(tracksGT(:,1)==ids(id,1),:);
0035 %     tr2 = tracksGT(tracksGT(:,1)==ids(id,2),:);
0036 %
0037 %     tr1 = tr1(tr1(:,2)>=rangeD(id,1) & tr1(:,2)<=rangeD(id,2),:);
0038 %     tr2 = tr2(tr2(:,2)>=rangeD(id,1) & tr2(:,2)<=rangeD(id,2),:);
0039 %
0040 %
0041 %     totalSize = [totalSize; size(tr1,1)+size(tr2,1) ];
0042 %
0043 % end
0044 
0045 
0046 for ra = 1: size(rangeD,2)
0047     tracksR = tracksGT(tracksGT(:,2) >=rangeD(ra,1) & tracksGT(:,2) <= rangeD(ra,2),:);
0048    
0049     idx = unique(tracksR(:,1));
0050    
0051    for id = 1: size(idx,1)
0052       track = tracksR(tracksR(:,1)==idx(id),:);
0053       
0054     
0055        
0056       totalSize = [totalSize size(track,1)];
0057    end
0058     
0059     
0060 end
0061 % nTR = length(totalSize);
0062 % nPO = sum(totalSize);
0063 
0064  
0065  
0066 
0067 
0068 
0069 
0070 
0071 
0072 
0073 
0074 
0075 
0076

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