3.1 - Affichage d'un calque de coordonnées via le Header¶
%matplotlib inline
from astropy.wcs import WCS
import numpy as np
from matplotlib.colors import LogNorm
import warnings
warnings.filterwarnings("ignore", category=UserWarning)#delete after also
#warnings.filterwarnings("ignore", category=Warning)#delete after also
#Create WCS object
wcs_for_plot = WCS(image_header)
#prepare plot
fig5 = plt.figure(figsize=(20, 20))
ax = plt.subplot(projection=wcs_for_plot) #add wcs informations on plot
plt.imshow(image_data, origin='lower', cmap='gray', aspect='equal',norm=LogNorm())
plot_title = "Gamma Cassiopae "
plt.title(plot_title ,fontsize=16, loc='left')
#add and configure axis informations
ax = plt.gca()
ax.grid(True, color='white', ls='solid')
ra = ax.coords[0]
dec = ax.coords[1]
ra.set_axislabel('RA (J2000)', minpad=1)
dec.set_axislabel('DEC (J2000)', minpad=1)
ra.set_major_formatter('hh:mm:ss.s')
dec.set_major_formatter('dd:mm:ss.s')
#add other grid as overlay
overlay = ax.get_coords_overlay('icrs') # or galactic, alt-az, ...
overlay.grid(True, color='lime', ls='solid')
overlay[0].set_axislabel('RA (deg)')
overlay[1].set_axislabel('DEC (deg)')
#uncomment the line after for show contour on objects
#ax.contour(image_data, transform=ax.get_transform(wcs_for_plot), levels=np.logspace(-3, 0, 3), colors='orange', alpha=1)
#show
plt.show()
3.2 - Récupération des coordonnées par le service astrometry.net¶
Pour pouvoir récupérer les coordonnées correspondant à l'image, une requête au service d'astrométrie (astrometry.net) est effectuée grâce à la librairie AstroQuery.
Conversion de l'image FITS en JPEG
L'image initiale faisant 80 Mo, il est possible de convertir cette image FITS en jpeg/png, avant de faire appel au service d'astrometrie, afin d'avoir une image plus légère d'environ 100ko et d'accélérer le processus.
#image name (FITS) to convert : image_path
#image name convert in jpeg/png : gamcas_jpeg.jpeg
from matplotlib.colors import LogNorm
%matplotlib inline
fig6 = plt.figure(figsize=(16, 9))
plt.imshow(image_data, cmap='gray', norm=LogNorm())
plt.grid()
plt.axis('off')
plt.savefig("dataset/gamcas_jpeg.jpeg", bbox_inches='tight')
Requêtage de l'API Astrometry.net
Il est nécessaire de préciser la clé d'API lié à votre compte. Cette dernière est disponible sur la page de votre compte, à cette adresse : http://nova.astrometry.net/api_help. L'avertissement présent est normal et demande de préciser notre configuration.
import warnings
warnings.filterwarnings("ignore", category=Warning)
#Import astroquery/astrometry library and World Coordinate System module from astropy
from astroquery.astrometry_net import AstrometryNet
from astropy.wcs import WCS
ast = AstrometryNet()
ast.api_key = '##############' # PUT YOUR API KEY HERE
try_again = True
submission_id = None
while try_again:
try:
if not submission_id:
wcs_header = ast.solve_from_image('dataset/gamcas_jpeg.jpeg',
submission_id=submission_id)
else:
wcs_header = ast.monitor_submission(submission_id,
solve_timeout=120)
except TimeoutError as e:
submission_id = e.args[1]
else:
# got a result, so terminate
try_again = False
#when solve is ok, record data in wcs_gamcas object
if wcs_header:
print('Solve OK - Print wcs result ...')
wcs_gamcas = WCS(wcs_header)
print(wcs_gamcas)
else:
print('nok')
Affichage des coordonnées célestes après le retour de l'API Astrometry.net
Dès lors, les coordonnées peuvent être affichées par-dessus l'image, dans le système de coordonnées souhaité. Plus d'informations sur le système de gestion WCS d'Astropy à ces adresses :
import matplotlib.image as mpimg
#Prepare plot
img = mpimg.imread('dataset/gamcas_jpeg.jpeg') #Image name_file
fig7 = plt.figure(figsize=(20, 20))
ax = plt.subplot(projection=wcs_gamcas) #Add wcs projection with wcs_gamcas file
plt.imshow(img, origin='lower', cmap='gray', aspect='equal',norm=LogNorm())
plt.title('gam Cas - Astrometry.net API',fontsize=16, loc='center')
ax.coords['ra'].set_axislabel('RA')
ax.coords['dec'].set_axislabel('DEC')
overlay = ax.get_coords_overlay('icrs')
overlay.grid(color='red', ls='dotted')
plt.show()
Récupération et affichage des coordonnées sans passer par l'API en ligne de code
Dans l'éventualité où le repérage astrométrique n'est pas effectué par l'API Python via le morceau de code précédent. L'utilisation classique via l'upload de fichier FITS sur le site astrometry.net permet de récupérer le fichier WCS unique, ou l'image FITS incluant l'en-tête WCS. Le code ci-dessous présente l'affichage des coordonnées sur l'image par cette méthode. Le fichier récupéré via le site astrometry.net s'appelle 'wcs.fit'
#Get wcs file path
wcs_file = fits.open('dataset/wcs.fits')
header_wcs_file = wcs_file[0].header
#Create WCS object
wcs_for_plot2 = WCS(header_wcs_file)
#Image plot must be the same than upload image on astrometry.net
import matplotlib.image as mpimg
img = mpimg.imread('dataset/gamcas_jpeg.jpeg')
fig8 = plt.figure(figsize=(20, 20))
ax = plt.subplot(projection=wcs_for_plot2)
plt.imshow(img, origin='lower', cmap='gray', aspect='equal',norm=LogNorm())
plot_title = "Gamma Cassiopae"
plt.title(plot_title ,fontsize=16, loc='left')
ax.coords['ra'].set_axislabel('RA')
ax.coords['dec'].set_axislabel('DEC')
overlay = ax.get_coords_overlay('icrs')
overlay.grid(color='red', ls='dotted')
plt.show()
