astropy WCS信息与fits文件中的头值不同

2024-05-01 21:41:27 发布

您现在位置:Python中文网/ 问答频道 /正文

我在使用fits文件时遇到问题。这个问题与wcs和我的文件头有关,作为参考,fits文件的轴是速度和度数。你知道吗

问题是WCS所说的和我的标题(正确的)所说的不一致。你知道吗

特别是,如果我这样做:

fits.open('file.fits')[0].header['CRVAL2']

,我得到了6012.0

fits.open('pv749290_gu.fits')[0].header['CDELT2']

,我得到4.0分

到目前为止还不错。当我这么做的时候问题就出现了

w = WCS('file.fits')

,因为我得到:

 CRVAL : 0.0  6012000.0 

 CDELT : 2.999833375699044  4000.0  

所以,正如你所看到的,我最初得到的CRVAL2和CDELT2的值突然大了3个数量级,然后这影响了我的图像绘制,因为我使用“w”作为投影来绘制我的轴。 有人能帮我解决这个问题吗?提前谢谢!你知道吗


Tags: 文件标题绘制open速度fileheaderwcs
1条回答
网友
1楼 · 发布于 2024-05-01 21:41:27

所以,如果有人有一天遇到同样的问题:

问题是,为了尽可能清楚地说明问题,我向文件的CUNIT2添加了一个值,即使该关键字最初不在头中。在本例中,我使用hdr['CUNIT2']='KM/S',但是在查看WCS时(文件.fits),CRVAL2的值似乎是m/s而不是km/s,因此我认为WCS的默认值与我给出的单位(?)之间存在某种张力。你知道吗

在任何情况下,通过再次移除标题的CUNIT2标签,并再次读取WCS(文件.fits)WCS和头之间的差异消失了,文件现在有了正确的维度,尽管没有在关键字中指定单位(当然,您可以向CRVAL2添加一个注释,明确指出单位)

相关问题 更多 >