月球/月相算法

时间:2010-03-26 21:16:22

标签: python c algorithm calendar astronomy

有没有人知道在给定日期计算月相或年龄的算法,或者找出给定年份中新月/满月的日期?

谷歌搜索告诉我答案在一些天文学书中,但当我只需要一页时,我真的不想购买整本书。

更新

我应该证明我对google搜索更好一点的说法。我确实找到了只在某个时间段内工作的解决方案(如1900年代);以及基于触发器的解决方案,其计算成本将超出我的预期。

S Lott在他的Python书中有几种算法用于计算给定年份的复活节,大多数算法少于十行代码,有些算法适用于格里高利历中的所有日子。在3月份找到满月是寻找复活节的关键因素,所以我认为应该有一个不需要触发的算法,适用于格里高利历中的所有日期。

8 个答案:

答案 0 :(得分:16)

我暂时将一些代码移植到Python。我本来只是链接到它,但事实证明它在此期间从网上掉了下来,所以我不得不去掉它并再次上传它。请参阅moon.pyJohn Walker's moontool

我无法找到这方面的参考资料,因为它在任何时间跨度都是准确的,但似乎作者非常严谨。这意味着是的,它确实使用了trig,但是我无法想象你将使用它会是什么,这将使它在计算上受到限制。 Python函数调用开销可能超过了trig操作的成本。计算机计算速度非常快。

代码中使用的算法来自以下来源:

Meeus,Jean。天文算法。 Richmond:Willmann-Bell,1991。ISBN 0-943396-35-2。

必备;如果你只买一本书,请确保它是这本书。算法以数学方式呈现,而不是计算机程序,但实现本书中许多算法的源代码可以与QuickBasic,Turbo Pascal或C中的发布者分开订购.Meeus提供了许多有用的调试计算示例你的代码,并经常提出几种算法,在准确性,速度,复杂性和长期(世纪和千年)有效性之间进行不同的权衡。

Duffett-Smith,彼得。实用的天文学与您的计算器。第3版。剑桥:剑桥大学出版社,1981。国际标准书号0-521-28411-2。

尽管标题中有计算器一词;如果您对开发计算行星位置,轨道,日食等的软件感兴趣,这是一个有价值的参考。提供了比Meeus更多的背景信息,这有助于那些不熟悉天文学的人学习常常令人困惑的术语。给出的算法比Meeus提供的算法更简单,更准确,但适用于大多数实际工作。

答案 1 :(得分:11)

答案 2 :(得分:11)

如果你像我一样,你会成为一个细心的程序员。因此,当您看到分散在互联网上的随机代码旨在解决复杂的天文问题时会让您感到紧张,但却无法解释为什么解决方案是正确的。

您认为必须有权威来源,例如 books ,其中包含仔细,完整的解决方案。例如:

  是的,Meeus,Jean。天文算法。 Richmond:Willmann-Bell,1991。ISBN 0-943396-35-2。

     达弗特 - 史密斯,彼得。实用的天文学与您的计算器。第3版。剑桥:   剑桥大学出版社,1981。国际标准书号0-521-28411-2。

您信赖广泛使用且经过良好测试的开源库,这些库可以纠正错误(与静态网页不同)。那么,这是基于PyEphem库的问题的Python解决方案,使用Phases of the Moon接口。

#!/usr/bin/python
import datetime
import ephem

def get_phase_on_day(year,month,day):
  """Returns a floating-point number from 0-1. where 0=new, 0.5=full, 1=new"""
  #Ephem stores its date numbers as floating points, which the following uses
  #to conveniently extract the percent time between one new moon and the next
  #This corresponds (somewhat roughly) to the phase of the moon.

  #Use Year, Month, Day as arguments
  date=ephem.Date(datetime.date(year,month,day))

  nnm = ephem.next_new_moon    (date)
  pnm = ephem.previous_new_moon(date)

  lunation=(date-pnm)/(nnm-pnm)

  #Note that there is a ephem.Moon().phase() command, but this returns the
  #percentage of the moon which is illuminated. This is not really what we want.

  return lunation

def get_moons_in_year(year):
  """Returns a list of the full and new moons in a year. The list contains tuples
of either the form (DATE,'full') or the form (DATE,'new')"""
  moons=[]

  date=ephem.Date(datetime.date(year,01,01))
  while date.datetime().year==year:
    date=ephem.next_full_moon(date)
    moons.append( (date,'full') )

  date=ephem.Date(datetime.date(year,01,01))
  while date.datetime().year==year:
    date=ephem.next_new_moon(date)
    moons.append( (date,'new') )

  #Note that previous_first_quarter_moon() and previous_last_quarter_moon()
  #are also methods

  moons.sort(key=lambda x: x[0])

  return moons

print get_phase_on_day(2013,1,1)

print get_moons_in_year(2013)

返回

0.632652265318

[(2013/1/11 19:43:37, 'new'), (2013/1/27 04:38:22, 'full'), (2013/2/10 07:20:06, 'new'), (2013/2/25 20:26:03, 'full'), (2013/3/11 19:51:00, 'new'), (2013/3/27 09:27:18, 'full'), (2013/4/10 09:35:17, 'new'), (2013/4/25 19:57:06, 'full'), (2013/5/10 00:28:22, 'new'), (2013/5/25 04:24:55, 'full'), (2013/6/8 15:56:19, 'new'), (2013/6/23 11:32:15, 'full'), (2013/7/8 07:14:16, 'new'), (2013/7/22 18:15:31, 'full'), (2013/8/6 21:50:40, 'new'), (2013/8/21 01:44:35, 'full'), (2013/9/5 11:36:07, 'new'), (2013/9/19 11:12:49, 'full'), (2013/10/5 00:34:31, 'new'), (2013/10/18 23:37:39, 'full'), (2013/11/3 12:49:57, 'new'), (2013/11/17 15:15:44, 'full'), (2013/12/3 00:22:22, 'new'), (2013/12/17 09:28:05, 'full'), (2014/1/1 11:14:10, 'new'), (2014/1/16 04:52:10, 'full')]

答案 3 :(得分:7)

另外,pyephem — scientific-grade astronomy routines [PyPI],这是一个Python包,但有computational guts in C,而

  

精度< 0.05“从-1369到+2950   使用表查找技术来限制对三角函数的调用。

答案 4 :(得分:2)

Pyephem默认使用协调通用(UTC)时间。我想要一个程序来生成一个在太平洋时区准确的满月列表。下面的代码将计算给定年份的满月,然后使用ephem.localtime()方法对其进行调整,以校准到所需的时区。它似乎也适当地考虑了夏令时。感谢理查德,这段代码与他所写的相似。

#!/usr/bin/python
import datetime
import ephem
import os
import time

# Set time zone to pacific
os.environ['TZ'] = 'US/Pacific'
time.tzset()

print("Time zone calibrated to", os.environ['TZ'])

def get_full_moons_in_year(year):
    """
    Generate a list of full moons for a given year calibrated to the local time zone
    :param year: year to determine the list of full moons
    :return: list of dates as strings in the format YYYY-mm-dd
    """
    moons = []

    date = ephem.Date(datetime.date(year - 1, 12, 31))
    end_date = ephem.Date(datetime.date(year + 1, 1, 1))

    while date <= end_date:
        date = ephem.next_full_moon(date)

        # Convert the moon dates to the local time zone, add to list if moon date still falls in desired year
        local_date = ephem.localtime(date)
        if local_date.year == year:
            # Append the date as a string to the list for easier comparison later
            moons.append(local_date.strftime("%Y-%m-%d"))

    return moons

moons = get_full_moons_in_year(2015)
print(moons)

上面的代码将返回:

Time zone calibrated to US/Pacific
['2015-01-04', '2015-02-03', '2015-03-05', '2015-04-04', '2015-05-03', '2015-06-02', '2015-07-01', '2015-07-31', '2015-08-29', '2015-09-27', '2015-10-27', '2015-11-25', '2015-12-25']

答案 5 :(得分:0)

我知道你正在寻找Python,但如果你能理解C#,那里有一个名为Chronos XP的开源项目,它可以很好地完成这项工作。

答案 6 :(得分:-1)

快速谷歌透露this

答案 7 :(得分:-1)

如果您不需要高精度,您可以随时(ab)使用月球(或阴阳极)日历类(例如,Microsoft .NET中的HijriCalendarChineseLunisolarCalendar)来计算(任何日期的月相,作为日历的“日期”属性,是月球(或月球)日历日,始终对应于月相(例如,第1天是新月,第15天是满月等。)