function pts = im2pts(filename)
% function pts = im2pts(filename)
% This function probably needs to be tuned for different images.
% The defaults work pretty well for LINK[1-3].JPG

[i,j]=find(edge(medfilt2(rgb2gray(imread(filename)), [5 5])));
pts = [i j];