Я использую астропию для загрузки и построения графиков из базы данных SDSS.Также я использую функцию Cutout2D, чтобы вырезать квадрат 10x10 arcsec вокруг интересующего меня объекта.Это код:
from astropy.coordinates import SkyCoord
from astroquery.sdss import SDSS
from astropy.wcs import WCS
from astropy import units as u
from astropy.nddata import Cutout2D
import matplotlib.pyplot as plt
import numpy as np
# Coordinates from Skycoord name search
objCoord = SkyCoord(ra=263.7552*u.degree, dec=57.052374*u.degree)
# Get SDSS images
xObj = SDSS.query_region(objCoord, spectro=True)
imgObj = SDSS.get_images(matches=xObj)
image_i = imgObj[0][0]
data_i = image_i.data
# Get image coordinates
wcs = WCS(image_i.header)
# Cut the Field of view 10.0 x 10.0 arcsec
FoV = np.array([10.0, 10.0])
FoV_dimen = u.Quantity((FoV[0], FoV[1]), u.arcsec)
cutout = Cutout2D(data_i, objCoord, FoV_dimen, wcs=wcs)
wcs_cut = cutout.wcs
# Plot the image
fig = plt.figure()
ax = fig.add_subplot(111, projection=wcs_cut)
imgRegion = ax.imshow(cutout.data, origin='lower', vmin=-0.64, vmax=0.45)
ax.scatter(objCoord.ra, objCoord.dec , s=150, linewidths=0.35, edgecolors='black', transform=ax.get_transform('world'), label='SHOC579')
# Format the image axis
ra, dec = ax.coords[0], ax.coords[1]
ra.set_axislabel('Right ascension')
dec.set_axislabel('Declination')
ra.set_ticks(spacing=1 * u.arcsec, color='white', exclude_overlapping=True)
dec.set_ticks(spacing=1 * u.arcsec, color='white', exclude_overlapping=True)
ra.grid(color='red', alpha=0.5, linestyle='solid')
dec.grid(color='black', alpha=0.5, linestyle='solid')
ax.legend()
plt.show()
Он генерирует это изображение:
Несмотря на то, что обе оси изображения были обрезаны с помощьюодинаковая длина и осевая сетка имеют одинаковое разрешение 1 "(1 угловая секунда), прямое восхождение покрывает почти вдвое больший ожидаемый угол.
Однако, когда вы смотрите на изображение небесного сервера SDSS также с 1" (1 угловая секунда)сетка:
Правая сетка восхождения - это та, которая выглядит ближе (около 6 угловых секунд диаметра для объекта)
Isчто-то мне не хватает?