def astropymethod2(img, hduid=1):
    hdu = fits.open(img)
w = WCS(hdu[hduid].header)
lon, lat = w.wcs_pix2world(100.0, 100.0, 1)
print(lon, lat)
Найдите физические координаты пикселя в подходящем файле с помощью python