English: Distance from Earth center to Moon center, data from JPL Horizons On-Line Ephemeris System
https://ssd.jpl.nasa.gov/horizons.cgi . The overlaying yellow disks represents the time of full moons.
Drawn by the following python3 code:
#Download the moon position in geocentric coordinates from https://ssd.jpl.nasa.gov/horizons.cgi with
# Ephemeris Type : VECTORS
# Target Body : Moon [Luna] [301]
# Coordinate Origin : Geocentric [500]
# Time Span : Start=2020-10-01, Stop=2022-10-01, Step=1 d
# Table Settings : defaults
# Display/Output : download/save (plain text file)
# and save to moon.txt
AU = 149597870.700
import re, math
f = open("moon.txt","r")
o = []
for line in f:
if line[0:2] == " X":
t = (re.sub('=+','=',line.strip().replace(" ","="))).split("=")
r = AU*math.sqrt(float(t[1])**2 + float(t[3])**2 + float(t[5])**2)
o.append(r)
f.close()
import matplotlib.pyplot as plt
plt.rcParams["figure.figsize"] = (14,4)
fig, ax = plt.subplots()
ax.plot(o[:570],'C1') # using only first 570 days in data
# full moon days counting from 2020-10-01, data from
# https://lunaf.com/lunar-calendar/2020/
# https://lunaf.com/lunar-calendar/2021/
# https://lunaf.com/lunar-calendar/2022/
fullmoonday = [1,31,61,91,120,149,179,209,238,268,297,326,356,385,415,445,474,504,534,563]
fullmoondistance = []
for i in range(len(fullmoonday)):
fullmoondistance.append(o[fullmoonday[i]])
ax.scatter(fullmoonday,fullmoondistance,color='y')
ax.set_box_aspect(0.2)
ax.set_ylabel('khoảng cách (km)')
plt.xticks([93,456],['1 tháng 1 năm 2021', '1 tháng 1 năm 2022'])
plt.ylim([350000, 415000])
plt.show()
#Then export the image to SVG
#Then use Inkscape to edit (adding annotations)