return index of all points in a epsilon neighborhood of a vector
[theta,rho,itheta,prho,rhomin] = getdata(S2G)