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
0021
0022
0023
0024
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
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
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