Welcome to the NavList Message Boards.

NavList:

A Community Devoted to the Preservation and Practice of Celestial Navigation and Other Methods of Traditional Wayfinding

Compose Your Message

Message:αβγ
Message:abc
Add Images & Files
    Name or NavList Code:
    Email:
       
    Reply
    Re: Skyfield
    From: Ken Muldrew
    Date: 2018 Oct 22, 17:01 -0600
    Robin Stuart wrote: " I’ve never understood why such a big deal is made about Python as to me it doesnt really seem to offer anything very much is new and that hasn’t been around in other languages for many years. "

    Two of the main reasons for its popularity among scientific programmers are the abundance of thoroughly tested libraries (such as Skyfield) and the simplicity of an interpreted language (hardware is fast enough now that for one-off programs, runtime is only a minor inconvenience). Even more convenient is the Jupyter notebook interface allowing people to easily share programs + text + graphics that are independent of computer system, setup, etc. (see some notebooks on http://nbviewer.jupyter.org/ for examples). And stack overflow (obviously).

    Ken Muldrew.

    On Mon, Oct 22, 2018 at 4:40 PM Robin Stuart <NoReply_Stuart@fer3.com> wrote:

    I recently found an application that makes it easy to generate positions of celestial bodies to the full accuracy of the JPL Ephemerides without the overhead of figuring out the file format and setting up the interpolation. I assume that there are others on Navlist who might find this useful. It’s called Skyfield https://rhodesmill.org/skyfield/ and it can also load the entire Hipparcos catalog with positions, visual magnitudes, proper motion and parallax which are incorporated when star positions are calculated. It’s written in Python. I’ve never understood why such a big deal is made about Python as to me it doesnt really seem to offer anything very much is new and that hasn’t been around in other languages for many years. Be that as it may, Skyfield makes things easy and generally does a nice job. Calculating the apparent geocentric position of the Moon for the equator of date is as simple as

    MoonRA, MoonDec, MoonDist = earth.at(t).observe(moon).apparent().radec(epoch=t)

    Results agree with those from the https://ssd.jpl.nasa.gov/horizons.cgi and MICA.

    As an exercise I generated a lunar almanac for 24 June, 1915 and compared it to the published Nautical Almanac. The hour here is astronomical time with 0h corresponding to 12h civil time. The time scale is UT1. I’ll leave it to you to decide whether the the NA of the time quotes more significant figures than are justified.

    Robin Stuart



    --
       \----------------------------+---------------------------------+   o_,
     O_/ \    Ken Muldrew, PhD      | Voice: (403) 226-6222           |   <\__/7
     <\__  \  Programming/Electrical| Cell:  (---) too-poor           |     | /
      "\ L  | Arbourealis RDC ltd.  | kbmuldrew@gmail.com             |   / /
       <    +-----------------------+---------------------------------+ / /
                   Morning coffee recapitulate phylogeny               L/

       
    Reply
    Browse Files

    Drop Files

    NavList

    What is NavList?

    Get a NavList ID Code

    Name:
    (please, no nicknames or handles)
    Email:
    Do you want to receive all group messages by email?
    Yes No

    A NavList ID Code guarantees your identity in NavList posts and allows faster posting of messages.

    Retrieve a NavList ID Code

    Enter the email address associated with your NavList messages. Your NavList code will be emailed to you immediately.
    Email:

    Email Settings

    NavList ID Code:

    Custom Index

    Subject:
    Author:
    Start date: (yyyymm dd)
    End date: (yyyymm dd)

    Visit this site
    Visit this site
    Visit this site
    Visit this site
    Visit this site
    Visit this site