wolffd@0: % Consider matching sources to detections
wolffd@0: 
wolffd@0: %  s1 d2  
wolffd@0: %         s2 d3
wolffd@0: %  d1
wolffd@0: 
wolffd@0: a  = optimalMatching([52;0.01])
wolffd@0: 
wolffd@0: % sources(:,i) = [x y] coords
wolffd@0: sources = [0.1 0.7; 0.6 0.4]';
wolffd@0: detections = [0.2 0.2; 0.2 0.8; 0.7 0.1]';
wolffd@0: dst = sqdist(sources, detections)
wolffd@0: 
wolffd@0: % a = [2 3] which means s1-d2, s2-d3
wolffd@0: a = optimalMatching(dst)
wolffd@0: 
wolffd@0: % a = [0 1 2] which means d1-0, d2-s1, d3-s2
wolffd@0: a = optimalMatching(dst')