Home > analysis > regionAnalysis.m

regionAnalysis

PURPOSE ^

SYNOPSIS ^

This is a script file.

DESCRIPTION ^

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 load _param/_wbc;
0002 P = initParam('wbc');
0003 
0004 FP = wbc.experiment.FP;
0005 TP = wbc.experiment.TP;
0006 
0007 fpRegionAverage = zeros(0,1);
0008 tpRegionAverage = zeros(0,1);
0009 fpRegionStd = zeros(0,1);
0010 tpRegionStd = zeros(0,1);
0011 for seq = 1: length(P.IFOLDERS)
0012     fprintf('\nAnalyzing at sequence %s', P.IFOLDERS(seq).Name);
0013     
0014     fpPerSeq = FP(FP(:,1)==seq,:);
0015     frs =  unique(fpPerSeq(:,2));
0016     
0017     for frId  = 1: length(frs)
0018         fr = frs(frId);
0019         im = im_get (P.IFOLDERS(seq).Name, fr,P.EXTENSION);
0020 %         tempDiffName= strcat(P.IFOLDERS(seq).Name,'/TD/TD_',...
0021 %             P.IFOLDERS(seq).Name,num2str(fr),'.png');
0022 %         im = imread(tempDiffName);
0023         
0024         % Read in  the detection mask
0025         maskName = strcat(P.IFOLDERS(seq).Name,'/DT/DT_',...
0026             P.IFOLDERS(seq).Name,num2str(fr),'.png');
0027         mask     = imread(maskName);
0028         [imL, nL] = bwlabel(mask,8);
0029         
0030         % Get FP per frame
0031         fpPerFrame = fpPerSeq(fpPerSeq(:,2) == fr,:);
0032         for fpId = 1 : size(fpPerFrame,1)
0033            fp = fpPerFrame(fpId,:);
0034            fpRegion = imL(round(fp(5)),round(fp(4)));
0035            
0036            regionVals = im(imL==fpRegion);
0037            fpRegionAverage(end+1) = mean(regionVals(:)); 
0038            fpRegionStd(end+1) = std(double(regionVals(:)));
0039         end
0040         
0041         tpPerFrame = TP(TP(:,1)==seq & TP(:,2)==fr,:);
0042         for tpId = 1: size(tpPerFrame,1)
0043            tp = tpPerFrame(tpId,:);
0044            tpRegion = imL(round(tp(5)),round(tp(4)));
0045            
0046            regionVals = im(imL == tpRegion);
0047            tpRegionAverage(end+1) = mean(regionVals(:));                       
0048            tpRegionStd(end+1) = std(double(regionVals(:)));
0049         end
0050         
0051         
0052     end
0053 end
0054 
0055 [tp tpOut] = hist(tpRegionAverage,256);
0056 [fp fpOut] = hist(fpRegionAverage,256);
0057 % S = horzcat(tp',fp');
0058 subplot(2,1,1);
0059 bar(tp);
0060 title('Average Temporal Difference Within Detection Boundary');
0061 legend('TP');
0062 xlabel('Temporal Difference Value');
0063 ylabel('Frequency');
0064 
0065 subplot(2,1,2);
0066 bar(fp,'r');
0067 title('Average Temporal Difference Within Detection Boundary');
0068 legend('FP');
0069 xlabel('Temporal Difference Value');
0070 ylabel('Frequency');
0071 
0072 
0073

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