python - 显示 FITS 文件内容

Question

我有一堆需要可视化的 FITS 文件格式的数据。

一些细节：

• 我使用python和astropy来操作和预览数据；
• 存储在 FITS 文件中的数据基本上是一个 numpy 数组，有 70 行（“订单”），每行 8096 像素（中阶梯状光谱）
• 数据保存为多页 pdf，其中每一页对应于 FITS 观察的一个特定顺序

• 我想按图 1 显示数据：

  • 图 1 对应于每个 FITS 文件中的一个订单；
  • 顶部面板上的“灰色”区域对应于没有数据/观察的区域；
  • 顶部面板上的每条“线”对应于不同的观察结果（x 轴：波长，y 轴：观察日期；z 轴：通量）
  • 红线是可选的
  • 底部面板与上面的数据相同，但所有观测值重叠（x 轴：波长，y 轴：通量）
  • 通量被归一化为面板上的中位数，但拟合值有时会远高于 10^7

现在，我面临以下问题。如果我保存为 pdf（甚至 png 等），我会受到我使用的 dpi 的限制。dpi 越高，我可以更好地预览数据，但由于文件大小，它变得无法使用。但在 dpi 较低的情况下，数据显得模糊。当我使用matplotlib 的显示预览数据时，我可以毫无问题地放大和缩小，但是当我在工作中生成图像时，它变得无法使用，并且无法远程进行。

所以，我的问题是：我可以使用类似于图 1 的文件格式来存储我的数据（理想情况下是多页格式，我可以使用 python 创建，但不限于。），但允许我使用类似于我对matplotlib的节目的“无限”分辨率？“在野外”有几个 FITS 文件查看器，但据我所知，它们不允许按我的意愿预览数据......

标签： pythonmatplotlibpdffits