Python code to fetch Carb nebula FITS image and Display

In Astronomy, there will be needs to extract Images from FITS format file for further analysis or research, Below is the Two steps code for "getting image link" and "fetching & displaying image from the link" 

The code uses python libraries such as PYVO,  Astropy, Matplotlib

STEP 1-excute below code to get url of the astronomical object

import pyvo as vo

from astropy.coordinates import SkyCoord

from astropy.units import Quantity

pos = SkyCoord.from_name('Eta Carina') #Carb Nebula

size = Quantity(0.5, unit="deg")

sia_service = vo.dal.SIAService("http://dc.zah.uni-heidelberg.de/hppunion/q/im/siap.xml")

sia_results = sia_service.search(pos=pos, size=size)


sia_results

 -------------- its will list all FITS format image from online repository------------

<Table length=8>
                                     accref                                      ...
                                                                                 ...
                                     object                                      ...
-------------------------------------------------------------------------------- ...
http://dc.zah.uni-heidelberg.de/getproduct/boydende/data/fits/HAR081_006454.fits ...
http://dc.zah.uni-heidelberg.de/getproduct/boydende/data/fits/HAR081_006478.fits ...
http://dc.zah.uni-heidelberg.de/getproduct/boydende/data/fits/HAR081_006507.fits ...
http://dc.zah.uni-heidelberg.de/getproduct/boydende/data/fits/HAR081_006510.fits ...
http://dc.zah.uni-heidelberg.de/getproduct/boydende/data/fits/HAR081_006523.fits ...
http://dc.zah.uni-heidelberg.de/getproduct/boydende/data/fits/HAR081_006481.fits ...
http://dc.zah.uni-heidelberg.de/getproduct/boydende/data/fits/HAR081_006498.fits ...
http://dc.zah.uni-heidelberg.de/getproduct/boydende/data/fits/HAR081_006479.fits ...

-----------------------------------------------------------------

STEP 2-select the image url which you want to display. copy past the url as shown below in bold part.

import matplotlib.pyplot as plt
from astropy.io import fits

# Open the FITS file
fits_file = 'http://dc.zah.uni-heidelberg.de/getproduct/boydende/data/fits/HAR081_006454.fits'  # Replace with your FITS file path
hdulist = fits.open(fits_file)

# Display the structure of the FITS file
hdulist.info()

# Get the image data from the primary HDU (Header Data Unit)
image_data = hdulist[0].data

# Close the FITS file
hdulist.close()

# Plot the image data using matplotlib
plt.imshow(image_data, cmap='gray')
plt.colorbar()
plt.title('FITS Image')
plt.show()

-----output----




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