"""Example conversion from pixel coordinates to ra/dec."""
from astropy.io import fits
from astropy.wcs import WCS

x, y = 250, 250

f = fits.open('f242_e0_c98.fits')
mywcs = WCS(f[0].header)
ra, dec = mywcs.all_pix2world([[x, y]], 0)[0]