I follow privong&#x27;s suggestion for stellarium for celestial body visualization.<p>However, if you want precision data as an input to your code, I would make use of NASA JPL&#x27;s ephemerides system [0]. It is precise and reliable. Another, more pythonic option is pyephem [1], which allows you to &quot;calculate positions of stars and planets&quot; at any given time.<p>Both options can be easily incorporated with astropy and APLpy, suggested by privong, for calculations and plotting.<p>[0] <a href="http:&#x2F;&#x2F;ssd.jpl.nasa.gov&#x2F;?ephemerides" rel="nofollow">http:&#x2F;&#x2F;ssd.jpl.nasa.gov&#x2F;?ephemerides</a><p>[1] <a href="https:&#x2F;&#x2F;pypi.python.org&#x2F;pypi&#x2F;pyephem&#x2F;" rel="nofollow">https:&#x2F;&#x2F;pypi.python.org&#x2F;pypi&#x2F;pyephem&#x2F;</a>