Невозможно создать и наложить контур из подгонки изображения с помощью aplpy - PullRequest
Новая
       29

Невозможно создать и наложить контур из подгонки изображения с помощью aplpy

0 голосов
/ 11 апреля 2019

Я пытаюсь наложить контур из файла подгонки, содержащего радиоизображение на рентгеновском изображении, используя функцию "show_contour", предоставляемую библиотекой aplpy.Мой код очень прост.

"program1" 

import aplpy 
import astropy.io.fits as fits
import astropy.wcs as wcs
hdul = fits.open('x-ray_image.fits')
fig = aplpy.FITSFigure(hdul[1])
fig.show_contour('radio_image.fits')

Я запускаю программу1, и возникает ошибка."сообщение об ошибке" 

WCS has more than 2 dimensions, so ``slice`` should be set

Я смотрю wcs radio_image.fits 

file_radio = fits.open('radio_image.fits')
w = wcs.WCS(file_radio[0])
print(w)

"сообщение о выходе" 

Number of WCS axes: 4
CTYPE : 'RA---SIN'  'DEC--SIN'  'FREQ'  'STOKES'
CRVAL : xxx.xxxxxxx  yy.yyyyyy  f  s  
CRPIX : 514.0  306.0  1.0  1.0  
NAXIS : 1090  634  1  1

x, y, f иs на самом деле числа.Я понимаю, что radio_image.fits имеет кубические данные.Поэтому я изменил program1, как показано ниже

"program2" 

import aplpy 
import astropy.io.fits as fits
import astropy.wcs as wcs
hdul = fits.open('x-ray_image.fits')
fig = aplpy.FITSFigure(hdul[1])
fig.show_contour('radio_image.fits',dimension=[0,1],slices=[f,s])

Но такая же ошибка возникает 

WCS has more than 2 dimensions, so ``slice`` should be set

Пожалуйста, скажите мне, как исправить мой код.

...