问题标签 [fits]

我有一个 Matlab 脚本，我正在尝试将其转换为 C++（见下文），因为它非常慢。我是 C++ 新手，开始我尝试使用codegen但它不起作用（我收到消息??? 此文本包含非空顶级表达式。它似乎是一个脚本。）

您对如何开始转换代码有什么建议吗？另外，最好的 C++ 函数是什么fitsread？

这是我的代码：

这个是专门为天文界设计的。

我有两个 .fit 格式的表，其中两个表都可能包含一个共同的列/标题（即“SDSS 标识符”）。现在，我想连接这些表，以便它生成一个新表，其中仅包含在两个表中找到的对象。这样，任何仅在一个表中的对象都将被丢弃。

我已经阅读了关于 table join astropy operations tutorials 的文档，但我想知道是否支持 .fit 文件（或者我是否需要 HDU 拟合、csv 或 ascii 等）以及某些不匹配的对象是否包含在最终合并表。

如果有人能就如何实现这一目标向我提出建议，那就太好了。

我有一个 CMB 全天地图的 FITS 文件，其中包含 T、Q、U 值。

使用函数healpy.fitsfunc.getformat(t)，输出是“数据类型t的FITS约定格式字符串”。

对于我的 FITS 文件，我得到“A29”。这是标准化的 FITS 约定吗？

我从一个 FITS-Table 读取数据，并获得一个包含该表的 Struct，其中每个标签代表一列。

有没有办法将数组结构重新格式化为结构数组？那么 Array 中的一个 Struct 代表一个 Row 吗？

@mgalloy 提出的一般解决方案（见下文）：

我有这个代码：

尝试在 Python 2.7 下运行时，一切正常。但是在 Python 3.4 下存在这个问题：

在此警告之后，我无法对我的文件执行操作。任何想法如何处理这个问题？

astropy.io.fits手册指出，我们可以使用超过 8 个字符的标题关键字。在这种情况下，将创建 HIERARCH 卡。该手册还指出，如果我们要存储超过 80 个字符的关键字-值对，将自动创建继续卡。

然而，在实践中，这两个定义似乎只能相互排斥，即我们不能创建包含关键字值对的 FITS 文件，其中关键字长于 8 个字符（即 HIERARCH 关键字）并且值非常长细绳。

一个例子：

在这里header1并且header2将是正确的，但是当调用repr(header3)或尝试保存具有此类标头的 FITS 文件时，ValueError: The keyword TEST TEST with its value is too long会引发错误。

这是 FITS 标准的“意料之外的功能”，即HIERARCH关键字不能在CONTINUE卡片上继续，或者这可能只是一个错误astropy.io.fits？

我目前正在尝试使用 Python 3.4 的 astropy.io.fits 从大型数据 FITS 格式的数据文件中读取数据。虽然我可以使用 memmapper 模式成功打开文件，但我无法访问文件中的数据。这是我的代码

在第二行table = hdu_list[1].data我得到OSError: [WinError 8] Not enough storage is available to process this command

我的想法是，将数据分配给变量table时，会尝试读取整个文件，从而导致存储错误。是否有任何解决方法，或者有一种方法可以一次简单地读取表数据的每一行？

我正在尝试使用下面的代码为 FITS 文件中的像素设置阈值。但是，我收到一条错误消息：

有想法该怎么解决这个吗？

这是代码：

我有一张适合的图像，里面有很多天文物体。我正在尝试从我感兴趣的对象中创建小的 4x4“图章”（对象周围的部分/图像）。我已经计算了原始拟合文件中对象的像素坐标，并创建了一个包含坐标的文档。我知道 imshow() 命令可能是最好的选择，但我对如何使用像素坐标来完成任务感到困惑。

第一个 imshow() 给了我原始的 fit 文档。

如果我想在图像中查看像一颗星这样的东西，同时知道它的像素坐标，我将如何去做。

我的长期目标是能够创建许多“邮票”，因为我有许多不同版本的同一天文图像，因此它们是在不同的滤镜中拍摄的。我希望能够切换出原始拟合文件，并为不同过滤器中的相同对象创建一组这些“标记”。

我正在尝试编写一个脚本，该脚本不仅可以读取 FITS 文件，而且可以打印并将数据保存到表中。到目前为止，我的代码读取文件似乎没有问题；打印它们似乎是问题所在。到目前为止，这是我的代码（当我运行此代码时，我收到错误消息“找不到此格式的有效行”：