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°在任意一天加减修正值即可。