如果允许中心线误差30至50公里，以上算法完全可以用来扩展计算地方食。

不过还需要一些坐标变换。

这个算法中，实际上已经给出了太阳和月亮的坐标，以及太阳月亮坐标在短时间内（如几个小时）的变化速度参数、日月直径参数等。这些参数称为“日食根数”。有了日食根数就可以计算日食了。此外，还需计算恒星时，它反应地球自转的位置。

 

此外，计算日食应在“方位坐标中计算”比较方便，如贝塞尔坐标就是一种方位坐标，再如站心坐标系（地平方位坐标）等。有些参数需在赤道坐标中计算，然后再转到方位坐标系中。

方位坐标的主轴一般为月影轴，也可以使用地日连线或观测点到太阳的连线。

由于涉及坐标变换，所以计算“是否可视”会比较麻烦一些。或许可以找到避免烦琐坐标变换的日食计算方法，但我没目前还没有想出来。好象chw以前曾说过他测算过日食（可能是半年前的帖子），据他说误差为几分钟，不然向他请教一些，说不定可以事半功倍。

感謝樓主的寶貴建議 . ^_^...

首先感谢许先生的代码和算法。
但是个人认为：
2010-07-11 05:37:32 T ac=1
这个2010-07-11不是初一啊，怎么会有日食呢？
是不是应该是
2010-07-12啊？

本人知识有限，但是认为这个时间应该是以北京时间表示的，那么就应该日食发生在初一【不管东经120°能不能看得到】呀。

许先生很久没来论坛了。

 

你所说的问题，其实是同一时刻的不同表达形式的问题。

 

许先生的软件，时间是用“力学时”表示的。

另外，从寿星天文历（5。04版）查得本次日食的信息：

本次日食概述(力学时)
偏食始： 2010-07-11 17:10:44 -161°16',- 11°38'
中心始： 2010-07-11 18:17:57 -171°03',- 26°53'
视午食： 2010-07-11 19:52:01 -116°20',- 22°28'
中心终： 2010-07-11 20:51:07 - 70°47',- 50°55'
偏食终： 2010-07-11 21:58:20 - 75°24',- 36°47'

 

换算成北京时间就是7月12日了。

哦，原来如此。是不是这样换算：
2010-07-11 19:52:01 +（8:00:00 - 00:01:07）
=2010-07-11 27:50:54
=2010-07-12 03:50:54
谢谢回答了。