How to calculate planet rotation/orientation
Posted: 24.09.2018, 01:50
I am creating a POV-Ray scene to model the Solar System:
http://lib.povray.org/searchcollection/index2.php ... emOrrery&contributorTag=SharkD
Using this PDF I was able to plot the planet positions:
https://ssd.jpl.nasa.gov/txt/aprx_pos_planets.pdf
However, now I want to rotate/orient each planet properly using rotational elements found here:
https://astropedia.astrogeology.usgs.gov/download/Docs/WGCCRE/WGCCRE2009reprint.pdf
I *think* this is how it works:
1. Position the planet so that the pole axis is aligned with the Z (up) coordinate axis.
2. Position the planet so that the the prime meridian points toward the negative X axis.
3. Rotate around the Z axis by W.
4. Rotate around the X axis by (90-delta_0) degrees.
5. Rotate around the Z axis by (90+alpha_0) degrees.
6. Rotate around the X axis by 23.43928 to get out of the ICRF frame and into the ecliptic frame.
7. Rotate around the Z axis by 90 degrees for reasons unknown to me.
These steps seem to work reliably for Earth and Mars (give or take one hour). But for Saturn they are off by about ~90 degrees around the Z axis (as far as I can tell). Am I following the wrong steps? Where am I making a mistake? Thanks!
http://lib.povray.org/searchcollection/index2.php ... emOrrery&contributorTag=SharkD
Using this PDF I was able to plot the planet positions:
https://ssd.jpl.nasa.gov/txt/aprx_pos_planets.pdf
However, now I want to rotate/orient each planet properly using rotational elements found here:
https://astropedia.astrogeology.usgs.gov/download/Docs/WGCCRE/WGCCRE2009reprint.pdf
I *think* this is how it works:
1. Position the planet so that the pole axis is aligned with the Z (up) coordinate axis.
2. Position the planet so that the the prime meridian points toward the negative X axis.
3. Rotate around the Z axis by W.
4. Rotate around the X axis by (90-delta_0) degrees.
5. Rotate around the Z axis by (90+alpha_0) degrees.
6. Rotate around the X axis by 23.43928 to get out of the ICRF frame and into the ecliptic frame.
7. Rotate around the Z axis by 90 degrees for reasons unknown to me.
These steps seem to work reliably for Earth and Mars (give or take one hour). But for Saturn they are off by about ~90 degrees around the Z axis (as far as I can tell). Am I following the wrong steps? Where am I making a mistake? Thanks!