Home > vis > find_neighbor_density.m

find_neighbor_density

PURPOSE ^

SYNOPSIS ^

This is a script file.

DESCRIPTION ^

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

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 % plot out 2 plots with nNeighbor and dNeighbor
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

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