方法如下：
1、先算两分日
比如在中国某地，杆影最短时是中午13点20分，且杆长与影长之比为1，则可知该地是北纬45°（tgα=1)，东经100°（从120°里1小时减15°，4分钟减1°）杆长与影长之比需查表求α，这里用了特殊角。
2、再算两至日经度的算法不变 纬度在北半球冬至α+23.5°，夏至α-23.5°在任意一天加减修正值即可。
3、修正值算法：就是距两分或两至日的天数差乘以94/365. 比如2013年2月17日，2013年3月22日春分差33天，即太阳直射点在南纬
33×94/365=8.5°
所以今天正午时得到的纬度是（arctgα+8.5）°
tgα= 杆长/影长