我在使用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”作为投影来绘制我的轴。 有人能帮我解决这个问题吗?提前谢谢!你知道吗
所以,如果有人有一天遇到同样的问题:
问题是,为了尽可能清楚地说明问题,我向文件的CUNIT2添加了一个值,即使该关键字最初不在头中。在本例中,我使用hdr['CUNIT2']='KM/S',但是在查看WCS时(文件.fits),CRVAL2的值似乎是m/s而不是km/s,因此我认为WCS的默认值与我给出的单位(?)之间存在某种张力。你知道吗
在任何情况下,通过再次移除标题的CUNIT2标签,并再次读取WCS(文件.fits)WCS和头之间的差异消失了,文件现在有了正确的维度,尽管没有在关键字中指定单位(当然,您可以向CRVAL2添加一个注释,明确指出单位)
相关问题 更多 >
编程相关推荐