1、计算地球上经纬度之间的距离d，已知地球上两点的经度、纬度：（X1,Y1), (X2,Y2)，其中X1,X2为经度，Y1,Y2为纬度。

2、计算程序需要转化为弧度（*3.1415926/180）地球半径为R=6371.0 km，则两点距离d=R*arcos[cos(Y1)*cos(Y2)*cos(X1-X2)+sin(Y1)*sin(Y2)]。

3、假设在中国某地，杆影最短时是中午13点20分，且杆长与影长之比为1，则可知该地是北纬45°（tgα=1)，东经100°（从120°里1小时减15°，4分钟减1°）杆长与影长之比需查表求α。再算两至日经度的算法不变 纬度在北半球冬至α+23.5°，夏至α-23.5°在任意一天加减修正值即可。