Question

基本上我只想使用已经合并到astropy中的pyfits从找到here的r-waveband拟合文件中获取图像数据。这是我的尝试：

from astropy.io import fits
import matplotlib.pyplot as plt

hdulist = fits.open("frame-r-004646-1-0019.fits")
hdulist.info()

imgplot = plt.imshow(hdulist[0].data)
imgplot.set_cmap('binary_r')

plt.show()
hdulist.close()

目前，我只能看到微小的点，而不是jpeg图像中看到的明亮的星星。此外，ImageHDU不包含图像（因为它是1D阵列所以不能包含图像）知道为什么会这样吗？

最后，我试图从（ra，dec）这个图像中识别我的目标星系。拟合文件具有与其相关联的（ra，dec）和每像素0.396弧秒的像素比例。因此，我想知道合并（ra，dec）所涉及的哪个角落，以便我可以在图像中推断目标的位置！

Answer 1

你会想要一些像样的颜色标准化：

imshow(im,vmin=np.percentile(im,5),vmax=np.percentile(im,95),cmap='gray')

其中im = hdulist [0] .data

为了渲染WCS，我建议你使用aplpy，如果你不知道像ds9这样的真实观众。

gc = aplpy.FITSFigure("frame-r-004646-1-0019.fits")
gc.show_grayscale(vmin=np.percentile(im,5),vmax=np.percentile(im,95))

从SDSS DR10的Fits文件中获取图像

1 个答案: