Сохранять координаты WCS при нарезке изображения RGB Fits - PullRequest
0 голосов
/ 21 февраля 2019

У меня есть коллекция изображений 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.

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...