使用Python ephem库计算月相-CSDN博客

本文链接：https://blog.csdn.net/weixin_28835551/article/details/114465187

该博客介绍了如何使用Python的ephem库来计算月相。通过导入ephem模块并选择Moon对象，可以创建一个仅关注月亮的实例，或者与其他天体一起考虑。示例代码展示了如何初始化和处理数据，以便在Python环境中进行天文学计算。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

# 需要导入模块: import ephem [as 别名]

# 或者: from ephem import Moon [as 别名]

def __init__(self, sun_only = False, moon_only = False):

if sun_only:

self.ephemerides = [ephem.Sun()]

self.data = None

elif moon_only:

self.ephemerides = [ephem.Moon()]

self.data = None

else:

self.ephemerides = [ephem.Venus(), ephem.Mars(), ephem.Jupiter(), ephem.Saturn(), ephem.Moon(), ephem.Sun()]

self.data = ascii.read(Configuration.star_catalog_file, guess=False, format='fixed_width_no_header', names=('HR', 'Name', 'DM', 'HD', 'SAO', 'FK5', 'IRflag', 'r_IRflag', 'Multiple', 'ADS', 'ADScomp', 'VarID', 'RAh1900', 'RAm1900', 'RAs1900', 'DE-1900', 'DEd1900', 'DEm1900', 'DEs1900', 'RAh', 'RAm', 'RAs', 'DE-', 'DEd', 'DEm', 'DEs', 'GLON', 'GLAT', 'Vmag', 'n_Vmag', 'u_Vmag', 'B-V', 'u_B-V', 'U-B', 'u_U-B', 'R-I', 'n_R-I', 'SpType', 'n_SpType', 'pmRA', 'pmDE', 'n_Parallax', 'Parallax', 'RadVel', 'n_RadVel', 'l_RotVel', 'RotVel', 'u_RotVel', 'Dmag', 'Sep', 'MultID', 'MultCnt', 'NoteFlag'), col_starts=(0, 4, 14, 25, 31, 37, 41, 42, 43, 44, 49, 51, 60, 62, 64, 68, 69, 71, 73, 75, 77, 79, 83, 84, 86, 88, 90, 96, 102, 107, 108, 109, 114, 115, 120, 121, 126, 127, 147, 148, 154, 160, 161, 166, 170, 174, 176, 179, 180, 184, 190, 194, 196), col_ends=(3, 13, 24, 30, 36, 40, 41, 42, 43, 48, 50, 59, 61, 63, 67, 68, 70, 72, 74, 76, 78, 82, 83, 85, 87, 89, 95, 101, 106, 107, 108, 113, 114, 119, 120, 125, 126, 146, 147, 153, 159, 160, 165, 169, 173, 175, 178, 179, 183, 189, 193, 195, 196))

# removed masked rows

self.data = self.data[:][~np.ma.getmaskarray(self.data['DE-'])]