0001 P = initParam();
0002
0003 pos = SS(1).Rich.tracksGT;
0004 distanceTable = Inf*ones(size(pos,1),size(pos,1));
0005
0006 for r = 1: size(pos,1)
0007 for c = 1 : size(pos,1)
0008 p1 = pos(r,3:4);
0009 p2 = pos(c,3:4);
0010 v = p2 - p1;
0011 distanceTable(r,c) = norm(v);
0012 end
0013 end
0014
0015 distanceTable(distanceTable > P.THR_DIS) = Inf;
0016 nNeighbor = zeros(size(distanceTable,1),1);
0017 dNeighbor = zeros(size(distanceTable,1),1);
0018
0019 for c = 1 :size(distanceTable,2)
0020 nNeighbor(c) = sum(distanceTable(:,c)~=Inf);
0021 for r = 1: size(distanceTable,1)
0022 if distanceTable(r,c)~=Inf
0023 dNeighbor(c) = dNeighbor(c) + distanceTable(r,c);
0024 end
0025 end
0026 end
0027
0028
0029 figure(1);imshow(zeros(1000,1000));hold on;
0030 for i = 1: size(pos,1)
0031 colorMapping = jet(max(nNeighbor(:)+1));
0032 colorV = colorMapping(nNeighbor(i),:);
0033 plot(pos(i,3),pos(i,4),'.','Color', colorV);
0034
0035 end
0036