Python Code to Plot position of planets on RA-DEC Grid

 Below is Python  Output & Code to "Plot position of planets on RA-DEC Grid".


Change the Setting parameters - Longitude, Latitude and Time.

The comments "#" gives short description of the code. You may need to install the required libraries.

***************************************************************************

from astropy.coordinates import EarthLocation, AltAz, get_sun, get_moon, solar_system_ephemeris, get_body

from astropy.time import Time

import astropy.units as u

import matplotlib.pyplot as plt


def plot_planet_positions(date_time, latitude, longitude):

    # Define observer's location

    observer_location = EarthLocation(lat=latitude*u.deg, lon=longitude*u.deg, height=0*u.m)

    

    # Define the time of observation

    observing_time = Time(date_time)

    

    # List of planets

    planets = ['mercury', 'venus', 'earth', 'mars', 'jupiter', 'saturn', 'uranus', 'neptune']

    

    # Lists to store RA and DEC values for each planet

    ra_values = []

    dec_values = []

    

    # Calculate RA and DEC for each planet

    with solar_system_ephemeris.set('builtin'):

        for planet_name in planets:

            planet = get_body(planet_name, observing_time, observer_location)

            ra_values.append(planet.ra.deg)

            dec_values.append(planet.dec.deg)

    

    # Plotting

    plt.figure(figsize=(10, 5))

    plt.scatter(ra_values, dec_values, color='blue')

    for i, planet_name in enumerate(planets):

        plt.text(ra_values[i], dec_values[i], planet_name.capitalize(), fontsize=9, ha='right')

    

    plt.xlabel('Right Ascension (degrees)')

    plt.ylabel('Declination (degrees)')

    plt.title('Positions of Planets on RA-DEC Grid')

    plt.grid(True)

    plt.gca().invert_xaxis()  # Invert RA axis to match astronomical convention

    plt.show()


# Example usage:

date_time = '2024-07-17 20:00:00'  # UTC time

latitude = 37.7749  # San Francisco, CA latitude

longitude = -122.4194  # San Francisco, CA longitude


plot_planet_positions(date_time, latitude, longitude)


***************************************************************************

Credits- ChatGPT 

Comments

Post a Comment

Popular posts from this blog

Distances used in Astronomy.

 There are various distances used to measure in Space/Astronomy. They are Mile, Kilometer, AU, Light year, Parsec. Mile and Kilometer - it is used for near by objects or Moon. AU(Astronomical Unit) -it is average distance between Sun & Earth. Probably used within solar system limits. Light year -Distance travelled by light in 1 year. It is used for larger distance beyond solar system. The distance between "Proxima Centauri" (the nearest star from earth) is 4.246 light years. Parsec - it is a unit of distance equal to 3.26 light-years. Cosmic distance ladder (extragalactic distance scale) - It is a technique of several methods used by Astronomers to determine the distance of cosmological bodies beyond our own galaxy, which are not easily determined using traditional methods. Learn more from below links- https://en.wikipedia.org/wiki/Cosmic_distance_ladder#Galactic_distance_indicators https://solarsystem.nasa.gov/news/1230/cosmic-distances/ https://en.wikipedia.org/wiki/As
Read more

Astronomy Python - Webb's First Deep Field - Marking the Edges of the Galaxies

Image
  In Astronomy, the image taken by telescopes of deep space field, galaxies, planets may appear to be complex - it may be difficult to find boundaries of objects when objects are too near. Or  simply one would like to know the shape of the object. Edge detection technique is used in image processing for finding the boundaries of objects within the image. Edge detection would be very useful while analyzing such difficult images. Below is the sample image taken by James Webb telescope. Image Credits: NASA, ESA, CSA, STScI Below is grey image(left) and edge detection image(right) An edge detection operato r used in the above example is called as Canny edge detector Below is the python code used. ********************************************************************** import cv2 import numpy as np import matplotlib.pyplot as plt # read the image image = cv2.imread("E:\\astronomy_related\\edge_detection\\webb_first_deep_field.jpeg") # convert it to grayscale gray = cv2.cvtColor(ima
Read more

Astronomy Python - Jupiter Planet - Cropping Red Gaint Spot

Image
 In Astronomy, there will be often need to crop(cut) certain part of image.  By using cropping technique, it becomes easy to highlight import part of the image, as shown in the below example,  Red Giant Spot is cropped and then shown in another window (top right). Image credit -Nasa Jupiter Below is the python code used. When code is executed original image appears and then selected the part of the image which you want to cut. ****************************************************************** import cv2 import numpy as np    # Read image image = cv2.imread("E:\\astronomy_related\\image_crop\\jupiter.jpg")    # Select ROI r = cv2.selectROI("select the area", image)    # Crop image cropped_image = image[int(r[1]):int(r[1]+r[3]),                        int(r[0]):int(r[0]+r[2])]    # Display cropped image cv2.imshow("Cropped image", cropped_image) cv2.waitKey(0) ********************************************************* Credits code credit- https://www.geeksfo
Read more