У меня есть коллекция изображений RGB Fits, и я хочу нарисовать синюю рамку с WCS на осях.Я нарезал изображение, чтобы взять синий, используя:
from astropy.wcs import WCS
from astropy.io import fits
image1 = fits.open('fits_file.fits')
data1 = image1[0].data
image1_slice_blue = data1[2,:,:]
Когда я строю это, оси соответствуют количеству пикселей, а не координатам WCS.Я попытался взять информацию WCS из файла, используя:
wcs1 = WCS(image1[0].header)
, но у меня есть проблема, потому что NAXIS = 3 (ширина изображения, высота и плоскости), но WCSAXES = 2 (RA и Dec).Ошибка говорит о том, что в основных ключевых словах WCS были обнаружены 3 измерения, и говорит мне выбрать / уменьшить размеры, используя naxis kwarg.Я запустил это снова, но с:
wcs1 = WCS(image1[0].header, naxis=2)
Но это тоже не работает, потому что CTYPE3 не найден.Может кто-нибудь помочь мне решить это?Спасибо
Это изображение, которое я пытался разделить: http://nova.astrometry.net/user_images/2677994#annotated
Я не могу поделиться настоящим файлом, но это файл "new-image.fits" в RHS.