0

我在处理 fit 文件时遇到问题。问题与我的文件的 wcs 和标题有关，并且对于信息，我的 fit 文件的轴是速度和度数。

问题是 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”作为投影来绘制我的轴。有人可以帮我解决这个问题吗？提前致谢！

4

1 回答 1

1

所以，如果有一天有人遇到同样的问题：

问题是，为了尽可能清楚，我向文件的 CUNIT2 添加了一个值，即使最初该关键字不在标题中。在这种情况下，我使用的是 hdr['CUNIT2']='KM/S'，但是在查看 WCS(file.fits) 时，CRVAL2 的值似乎是 m/s 而不是 km/s，所以我认为 WCS 的默认值和我给出的单位之间存在一些紧张关系（？）。

在任何情况下，通过再次删除标题的 CUNIT2 标签，并再次读取 WCS(file.fits)，WCS 和标题之间的差异消失了，文件现在具有正确的尺寸，尽管单位没有在关键字（当然，您可以在 CRVAL2 中添加注释，明确说明单位。）

于 2018-11-13T09:54:07.307 回答