Celestial Nav "sight reductions" w/ python ?

[captcha1]

I bought a used Davis 25 sextant a few days ago, and I'm interested to do "sight reductions" w/ python.

George Kaplan from USNO writes about using "least squares" to calculate sight reductions : https://gkaplan.us/content/nav_algorithms.html & https://aa.usno.navy.mil/publications/reports . I'd like to implement this in python.

Have you seen python code that does sight reductions ?

I see various python software that could be useful, eg : numpy.linalg.lstsq , scipy.linalg.lstsq , scipy.optimize.minimize

I found almanac data in python : pymeeus, pyephem & skyfield.

I also found a recent paper about using Kalman filters : https://doi.org/10.1017/S0373463321000758 ...