I'm from Spain and I'm currently finishing my degree in Aerospace Engineering. I am starting with all of this. My objective is to obtain the Chebyshev coefficients for the major planets position, relative to a determined timespan.

Im now using JPL Horizon's software and I'm feeling a bit lost. I only know how to get the vectors and Keplerian elements of the bodies I'm trying to study, but not the appropriate Chebyshev coefficients. Im using the telnet interface. What steps should I follow to get that tables?

  • $\begingroup$ Hello @Diego Suárez González and *Welcome to Space! If you have some success feel free to stop by and write a short answer about what worked best for you. Answering your own question is always okay in Stack Exchange, and you can click accept on your own answer as well. That's particularly helpful because it will show future readers how you were able to solve your problem. $\endgroup$ – uhoh Mar 23 '19 at 0:08

The best and most concise answer is here

The "gold standard" is the SPICE toolkit. You can read more about using spice in this excellent answer.

The JPL Development Ephemerides are lists of the Chebyshev coefficients you seek, along with some additional information. They can be downloaded as text files or binary. If you are familliar with Python you can have a look inside the Skyfield package (Github) to see how it's done there.

Here's a quick way to get DE421: https://pypi.org/project/de421/

Here they are in ascii format: ftp://ssd.jpl.nasa.gov/pub/eph/planets/ascii

Further reading:

according to ftp://ssd.jpl.nasa.gov/pub/eph/planets/ascii/ascii_format.txt:


On the first record of an export binary file or in the "GROUP 1050 of the ascii
"header", there are 3 sets of 15 integers each.  (Older versions have only the first 13
integers in each set)

The 15 triplets give information about the location, order and time-coverage of
the chebychev polynomials corresponding to the following 15 items:

       Earth-Moon barycenter
       Moon (geocentric)
       Earth Nutations in longitude and obliquity (IAU 1980 model)
       Lunar mantle libration
       Lunar mantle angular velocity
       TT-TDB (at geocenter)

Word (1,i) is the starting location in each data record of the chebychev 
coefficients belonging to the ith item.  Word (2,i) is the number of chebychev 
coefficients per component of the ith item, and Word (3,i) is the number of 
complete sets of coefficients in each data record for the ith item.

Items not stored on the ascii files have 0 coefficients [Word (3,i)].

Data Records ("GROUP 1070")

These records contain the actual ephemeris data in the form of chebychev 

The first two double precision words in each data record contain

         Julian date of earliest data in record.
         Julian date of latest data in record.

The remaining data are chebychev position coefficients for each component of 
each body on the tape.  The chebychev coefficients for the planets represent 
the solar system barycentric positions of the centers of the planetary systems.

There are three Cartesian components (x, y, z), for each of the items #1-11; 
there are two components for the 12th item, nutations : d(psi) and d(epsilon);
there are three components for the 13th item, librations : phi, theta, psi;
there are three components for the 14th item, mantle omega_x,omega_y,omega_z;
there is one component for the 17th item, TT-TDB.

Planetary positions are stored in units of kilometers (TDB-compatible).
The nutations and librations are stored in units of radians.
The mantle angular velocities are stored in radians/day.
TT-TDB is stored in seconds.
| improve this answer | |
  • $\begingroup$ Thank you for your answer, but i have some doubts remaining. How can I read that .bsp files? As I understood, those files are binary and need to be converted but, is there a "reader" to see the original file structure described in the quote above? or I can only use the spice toolkit to obtain the chebyshev polynomials? $\endgroup$ – Diego Suárez González Apr 3 '19 at 20:22
  • $\begingroup$ @DiegoSuárezGonzález If I remember correctly there are readers in several languages (FORTRAN, C...) available there, and if you look in Github you can find more. As I mentioned in the answer for planet positions (and a lot more) you can also use the Python package Skyfield which downloads and reads them as well. You should do some exploring of these options and see what works best for you. I can't choose for you because I won't know what suits you best. If you have specific problems, ask a new question. There's no limit to the number of good questions you can ask! $\endgroup$ – uhoh Apr 3 '19 at 22:32

You don't get those from HORIZONS. You download them from here. They are best accessed with the SPICE toolkit, or you can roll your own to extract the Чебышёв coefficients from the SPK files.

| improve this answer | |
  • $\begingroup$ It's good you've chimed in; I've added a pointer to this answer from mine. $\endgroup$ – uhoh Mar 23 '19 at 4:56

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.