plot

import numpy as np
import matplotlib.pyplot as plt
from astropy.io import fits
from function import listdirfile
path = './J125035+002000_tmp/'
wcslist = listdirfile(path,'.wcs')
ls_ra = []
ls_dec = []
wcslist.sort()
for i in wcslist:
    print(i)
    wcsfile = path +i
    hdul = fits.open(wcsfile)
    hdr = hdul[0].header
    crval1 = hdr['CRVAL1']
    crval2 = hdr['CRVAL2']
    ls_ra.append(crval1)
    ls_dec.append(crval2)
n=np.arange(len(ls_ra))
fig,ax=plt.subplots()
ax.scatter(ls_ra,ls_dec,c='r')
for j,txt in enumerate(n):
    ax.annotate(txt,(ls_ra[j],ls_dec[j]))

plt.plot(ls_ra,ls_dec)
plt.show()