我有一个有关天文物体的拟合文件。我可以这样绘制:

from astropy.io import fits
from astropy.wcs import WCS

hdul = fits.open(fitsfilename)[0]

wcs = WCS(hdul.header)
fig = plt.figure(figsize=(12,12))
fig.add_subplot(111, projection=wcs)
plt.imshow(hdul.data)


这有效,并产生一个漂亮的图片:

python - 为什么WCS投影子图上的功能在matplotlib中放在错误的位置?-LMLPHP

我想在此图中添加一些其他功能,但无法正常工作。例如,让我们尝试将圆添加到119°,-67°30'。我通过以下方式扩展代码:

plt.scatter([119],[-67.5],c='r',s=500)


我得到的是:

python - 为什么WCS投影子图上的功能在matplotlib中放在错误的位置?-LMLPHP

确实不是我们想要的,圆大约在118°5',-67°5',而不是应该在的位置(119°,-67°30')。

我出了什么问题,或者什么是解决此问题的好方法?



注意:当我运行wcs = WCS(hdul.header)时,会收到警告:


  警告:Verify警告:验证报告错误:
  [astropy.io.fits.verify]警告:VerifyWarning:卡“ A_2_0”不是
  FITS标准(无效值字符串:'3.29341755408e-05')。固定
  符合FITS标准的“ A_2_0”卡。 [astropy.io.fits.verify]
  警告:Verify警告:注意:astropy.io.fits使用从零开始
  索引。 [astropy.io.fits.verify]警告:Verify警告:卡
  “ A_1_1”不是FITS标准(无效值字符串:
  '1.51709339878e-05')。修复了“ A_1_1”卡符合FITS标准的问题。
  [astropy.io.fits.verify]警告:验证警告:卡“ A_0_2”不是
  FITS标准(无效值字符串:“ 5.17973753556e-06”)。固定
  符合FITS标准的“ A_0_2”卡。 [astropy.io.fits.verify]
  警告:验证警告:卡“ B_2_0”不是FITS标准(无效
  值字符串:“ 2.97627426087e-06”)。固定了“ B_2_0”卡以满足
  FITS标准。 [astropy.io.fits.verify]警告:Verify警告:卡
  “ B_1_1”不是FITS标准(无效值字符串:
  '2.71948126373e-05')。修复了“ B_1_1”卡符合FITS标准的问题。
  [astropy.io.fits.verify]警告:Verify警告:卡“ B_0_2”不是
  FITS标准(无效值字符串:“ 1.66848449653e-05”)。固定
  符合FITS标准的“ B_0_2”卡。 [astropy.io.fits.verify]
  警告:验证警告:卡“ AP_1_0”不是FITS标准(无效
  值字符串:“ 1.79541533196e-06”)。固定了“ AP_1_0”卡以满足
  FITS标准。 [astropy.io.fits.verify]警告:Verify警告:卡
  “ AP_0_1”不是FITS标准(无效值字符串:
  '9.20624843151e-07')。修复了“ AP_0_1”卡符合FITS标准的问题。
  [astropy.io.fits.verify]警告:VerifyWarning:卡“ AP_2_0”不是
  FITS标准(无效值字符串:“-3.29292923201e-05”)。固定
  “ AP_2_0”卡符合FITS标准。 [astropy.io.fits.verify]
  警告:验证警告:卡“ AP_1_1”不是FITS标准(无效
  值字符串:“-1.51738446887e-05”)。修复了“ AP_1_1”卡可以满足
  FITS标准。 [astropy.io.fits.verify]警告:Verify警告:卡
  “ AP_0_2”不是FITS标准(无效值字符串:
  '-5.18321445978e-06')。修复了“ AP_0_2”卡符合FITS标准的问题。
  [astropy.io.fits.verify]警告:VerifyWarning:卡“ BP_1_0”不是
  FITS标准(无效值字符串:'8.99029048217e-07')。固定
  “ BP_1_0”卡符合FITS标准。 [astropy.io.fits.verify]
  警告:验证警告:卡“ BP_0_1”不是FITS标准(无效
  值字符串:“ 1.15967736014e-06”)。修复了“ BP_0_1”卡满足
  FITS标准。 [astropy.io.fits.verify]警告:Verify警告:卡
  “ BP_2_0”不是FITS标准(无效值字符串:
  '-2.97837492348e-06')。修复了“ BP_2_0”卡符合FITS标准的问题。
  [astropy.io.fits.verify]警告:VerifyWarning:卡“ BP_1_1”不是
  FITS标准(无效值字符串:“-2.71998518336e-05”)。固定
  符合FITS标准的“ BP_1_1”卡。 [astropy.io.fits.verify]
  警告:VerifyWarning:卡“ BP_0_2”不是FITS标准(无效
  值字符串:“-1.66872388359e-05”)。修复了“ BP_0_2”卡满足
  FITS标准。 [astropy.io.fits.verify]警告:Verify警告:卡
  “ WCSR_PRJ”不是FITS标准(无效值字符串:“ 3.6679e-07”)。
  修复了“ WCSR_PRJ”卡符合FITS标准的问题。
  [astropy.io.fits.verify]警告:VerifyWarning:卡“ WCSR_PIX”为
  不是FITS标准(无效值字符串:'8.2565e-05')。固定
  “ WCSR_PIX”卡符合FITS标准。 [astropy.io.fits.verify]


因此,这可能是相关的;如何解决它的问题仍然存在。

最佳答案

要在world coordinates中进行绘图,您需要指定transform,例如:

ax = fig.gca()
ax.scatter([34], [3.2], transform=ax.get_transform('world'))


通常,您可以忽略那些FITS标头警告,因为那里没有FITS标头卡与WCS(afaik)相关。

关于python - 为什么WCS投影子图上的功能在matplotlib中放在错误的位置?,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/54109929/

10-09 08:21