d(x1,y1,x2,y2)=r*arccos(sin(x1)*sin(x2)+cos(x1)*cos(x2)*cos(y1-y2)) x1,y1是纬度/经度的弧度单位,r为地球半径 而当y1=y2时,公式就变为: d=r*|x1-x2|