我如何计算全食的路径？

Question

在观察了2015年3月20日的全食之后，我开始考虑使用pyephem计算全食的路径。

我的一般方法是计算最大日食的时间和地点，然后从那一点开始前后工作。通过从最大的eclipse参数开始，我可以对此进行检查：

所以考虑到时间，月亮和太阳的RA / dec，我将如何计算最大日食的纬度和经度？