Parameters to be added at the log.

[Bmorgado19]

The log was intended to be a text to help the user to save all the parameters that may be needed for future analysis, without having to redo any step. Here are some suggestions to be added in future versions:

• Apparent polar radius: polar_radius = equatorial_radius * (1 - oblateness)
• Apparent equivalent radius: equivalent_radius = np.sqrt(equatorial_radius * polar_radius)
• Solar separation angle: Occultation.predict['S-G-T']
• Lunar separation angle: Occultation.predict['M-G-T']
• Height of the observation for each station (Observer)
• Lunar separation angle for each station: M-O-T
• Solar separation angle for each station: S-O-T
• Geometric Albedo
• Chord size

Please add your suggestion at this issue

[altairgomes]

I think it's also important to add the height above the horizon each observation took place. This could be the mid-instant of occultation for positive chords and meantime of observation for negative ones

[Bmorgado19]

• Height of the observation for each station (Observer)

Please add your suggestion at this issue

[Bmorgado19]

The geometric albedo can be calculated using:
geometric_albedo = (10**(0.4*(H_sun - H_body))) * ((u.au.to('km')/equivalent_radius)**2), where H_sun = -26.74 and H_body = Ephem.H. See Sicardy et al. (2011, suplementary material) and references therein.

[gugabrossi]

The geometric albedo (p) can be calculated using the equation:
p = (10**(0.4*(H_sun - H_body))) * ((u.au.to('km')/equivalent_radius)**2)
where H_sun is the apparent of the sun at 1 au

H_sun-V=-26.74 and H_sun-R=-27.10
(See Wilmer, 2018 - AJ for236:47 (14pp) for a whole set of Sun's absolute magnitude - m_sun - in different filters, where one can derive its apparent magnitude using the equation
H_sun = m_sun - 2.5 log[ (d/10)^2 ]
where d is 1 au in parsec)

The geometric albedo (p) can also be given by a simpler relation:
D_km = (1329/ √ p )* 10^(−0.2*H)
where D_km is the object diameter; H is the object absolute magnitude