问题标签 [astropy]

我正在使用 astropy.table.write(filename, path=run_dir) 将 astropy 表写入名为 dat.h5 的文件。但是我收到了文件存在的错误，我在下面显示了它不存在的 pdb 跟踪。到底是怎么回事？

我正在尝试根据从Header. 但是，以我对 Python 的了解和astropy. 我的代码是这样的：

这latitude1将是一个以度为单位的值，在从标题中调用后被识别。所以我需要完成两件事：

1. 如何调用标题来识别银河纬度？
2. 以这样的方式拼接数组，使其仅包含纬度范围的值，其他所有值都为 0。

关于在大表中实现屏蔽的两个相关问题。

我正在处理一个包含 50 多列（从不同来源合并）的大表。列的缺失值为 NaN、0 或空。我想创建一个掩码来查找这些值并为所有值使用相同的缺失值。如果不为每列创建掩码，我无法找到一种方法。

astropy.table.Table 还允许在 Web 浏览器中显示可排序的表格，使用：

将缺失值显示为空白或“--”会很棒。谢谢你的帮助。

我正在为一个项目做一些目录匹配，发现使用astropy与其他软件包进行坐标匹配（例如Topcat）相比，我得到的结果是不同的。特别是，我不明白为什么下面的代码会产生不对称的结果

这打印：10547 10553

这里我只是颠倒目录顺序，得到不同的结果。特别是，对于我的项目，差异要大得多，我不确定问题出在哪里。

谢谢！

我正在使用 Python 2.7.10 和 pyfits 3.3。早些时候，我使用以下代码合并两个表。但是，现在我遇到了一些错误

错误是：

我需要从Distance astropy 类中访问浮点值。

这是一个 MWE：

这会产生一个<class 'astropy.coordinates.distances.Distance'>对象列表：

我需要存储浮点数（而不是对象），但我不知道如何访问它们。由于这个 MWE 是更大代码的一部分，我不能只做d.append(_). 我需要从Distance类生成的对象中访问浮点数。

添加：

我尝试将列表转换为 numpy 数组：

但我得到：

我相信解决这个问题会相对简单，但我似乎无法弄清楚如何对我在 python 中绘制的散点图进行卷积。

我有 2 个数据数组，一个是银河纬度，一个是银河经度，我用锤子投影绘制了它们，以表示银河坐标中恒星的分布。

现在，我想使用 boxcar 平滑来平滑 15 度框的图。我曾尝试将 astropy.convolution 与 convolve 和 Box2DKernel 一起使用，但我似乎无法使其工作。我还查看了来自http://docs.astropy.org/en/stable/convolution/kernels.html的示例， 但我不明白如何将他们的示例转换为我需要做的事情。他们似乎正在绘制一个二维函数并对其进行平滑处理。我可以不卷积一个图并根据它们在图表上的位置将这些点分类吗？我唯一能显示任何东西的东西都会产生一条直线，我不明白为什么。我对python很陌生，所以这给我带来了很多麻烦。

这是我到目前为止的代码：

这将两个数组绘制成一个锤子投影：

这是我对数据进行卷积的尝试：

就像我说的那样，这只是最终给出一些东西而不是错误的众多尝试之一，但我不确定我是否完全正确地使用了该功能。我不确定（或者我不认为）我可以像我那样创建“数据”，但是任何其他数组组合或将每个数组卷积为一维数组也不起作用。

任何想法都会非常有帮助，谢谢。

我有这个 Fortran 程序来计算谱线的等效宽度我希望找到帮助编写 python 代码来执行相同的算法（输入文件包含两列波长和通量）

我有一个包含元素的巨大表数据（或记录数组）：

它们都是整数，并且i是 0 到 100 万（表的大小）之间的随机数。

我还有一个名为的列表Name，其元素是文件的所有名称（总共 900 个名称），例如'/Users/Desktop/Data/spe-3588-55184-0228.jpg'（修改的），都包含三个数字。

现在我想从我tbdata的上面提到的元素中选择那些与 list 名称中的三个数字匹配的数据Name。这是我最初编写的代码：

Python 跑了一夜，仍然没有完成，因为要么数据量很大，要么我的算法太慢......我想知道完成这样一项任务的最快方法是什么？谢谢！

我想我在从 fit 文件中选择数据时搞砸了……基本上我所做的是：

并应用了选择规则：

有没有办法将 A 中的数据重新组合成一个新的 fit 文件？我尝试使用 Pandas 进行操作，但找不到正确的方法...