-
Notifications
You must be signed in to change notification settings - Fork 35
Open
Description
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()