Skip to content

Improve HG, HG12, etc. documentation #386

@mkelley

Description

@mkelley

The documentation, particularly https://sbpy.readthedocs.io/en/stable/sbpy/photometry.html, but also the class docstrings, seems to be missing two key concepts regarding the disk integrated phase function models:

  • That the models are based on the astropy modeling framework.
    • The users should be directed to the astropy modeling pages for more help on using these models.
    • One especially important point that should be documented up front: calling the object with a phase angle produces the phase corrected absolute magnitude.
  • That the models do not include rh and delta corrections, i.e., they do not produce apparent magnitudes.

These two points could be illustrated with an example that plots apparent magnitude versus time. Perhaps a simplified version of the following:

import numpy as np
import matplotlib.pyplot as plt
from astropy.time import Time
import astropy.units as u
from astropy.visualization import time_support
from sbpy.photometry import HG1G2
from sbpy.data import Ephem

time_support()  # allow plotting Time objects with matplotlib

epochs={
    "start": Time("2022-09-26 23:14"),
    "stop": Time("2023-03-01"),
    "step": 1 * u.day,
}
eph = Ephem.from_horizons("didymos", epochs=epochs, location="500")
m_HG1G2 = HG1G2(H=18.11, G1=0.84, G2=0.05)

m = m_HG1G2(eph["phase"]) + 5 * np.log10(eph["rh"] / u.au * eph["delta"] / u.au)

plt.clf()
ax = plt.gca()
plt.plot(eph["date"], m)
ax.invert_yaxis()

Metadata

Metadata

Assignees

Type

No type

Projects

No projects

Milestone

Relationships

None yet

Development

No branches or pull requests

Issue actions