问题标签 [astropy]

我尝试使用 satrapy 和以下代码从 FITS 文件中读取世界坐标系（WCS）：

我得到错误：

头文件是：

我已经将它与头文件中的其他 WCS 进行了比较，它看起来确实不同，但我确定如何解决问题。我也不确定如何使用RADESYSa而不是RADECSYS. 我似乎找不到任何文档。任何帮助，将不胜感激。

在以前的天文学版本中，可以通过以下方式处理不确定性的传播：

对 NDData 的更改似乎已删除此功能。我得到“AttributeError：'NDData'对象没有属性'add'”，我在文档中找不到任何有用的建议。现在如何处理错误传播？

我正在使用直接卷积算法来计算此图像之间的卷积：

这个内核：

我正在使用astropy中的实现进行直接卷积。

这将导致以下卷积，将所有设置（包括边界处理）保留为默认值，即 astropy.convolution.convolve(image,kernel)：

这个卷积有一些令人费解的伪影。特别是，在距边缘约 50 个像素的偏移处有一个“正方形”图案。在我看来，这是由于内核的范围所致；尽管内核的正式大小为 249x249，但大多数信息显然包含在大约 100 像素的半径内——这意味着当内核应用于边缘时，我们可能会遇到麻烦。

这让我想到了我的问题：

1. 这个假设是否正确——它确实是一个边缘问题？
2. 我将如何解决这个问题？我不知道如何证明使用不同的边缘处理（零填充、插值、换行......）我确信不同的情况需要不同的解决方案，但我不确定如何决定这个......
3. 只是...试图了解使用直接算法和 FFT 卷积之间的区别。如果内核和图像同样大，则 FT 卷积不需要零填充，不会出现边缘效应。对于直接方法，您会不经意间进行一些边缘处理……那么结果是否相等？因为原则上只有他们的表现应该不同，对吧？

我正在尝试使用 astropy.io 将数据附加到适合文件。

这是我的代码示例：

我得到的错误是

1. 这可能是我正在使用的 astropy.io 版本的问题吗？
2. 有没有更简单的方法可以使用 astropy.io 向 fit 文件添加扩展名或列？

任何帮助，将不胜感激。

我有两个具有赤道坐标(ra, dec)且(ra_0, dec_0)位于距离r和的源r_0，我需要计算它们之间的 3D 距离。

据我所知，我使用了两种应该给出相同结果的方法，但事实并非如此。

第一种方法是应用astropy的separator_3d函数。第二种方法是使用用球坐标给出两个源之间距离的表达式：

如此处所示。

在下面的 MCVE 中，返回的值为：

这两个值不应该相等吗？

MCVE：

I have tried debugging my code and I've realised that ultimately it breaks down when I try to save my AltAz coordinates into a .csv file because its not a numpy array, its a SkyCoord object. Could someone suggest a simple of way of converting a large table of Equatorial coordinates to AltAz or how I can get my code to save to file.

我在这里看到了几个关于 PyEphem 的答案，以及它如何产生日落/日出时间，但是如果我能找到一个仅使用 Astropy 包的解决方案，那对我来说会更有用。目前我找到的最接近的是astropy.coordinates包中的(get_sun)函数。这个日落/日出功能是否已内置到 Astropy 中？

我一直在尝试使用我在网上找到的修改后的脚本来绘制拟合图像的径向轮廓。我总是得到与预期完全不同的 y 轴单位。我什至不确定y轴单位是什么。我附上了拟合文件和我不断获得的轮廓以及我使用另一个程序绘制的正确径向轮廓。

我对python很陌生，所以我不知道为什么会这样。任何解决此问题的帮助将不胜感激。

这是我一直在使用的代码：

带有错误 y 轴单位的配置文件

使用 Celtech Aperture Photometry Tool 创建的具有预期正确单位的配置文件。

我正在尝试将 ascii 表读入 Python 中的 Numpy/Pandas/Astropy 数组/数据帧/表。表中的每一行如下所示：

问题是列之间没有明确的分隔符/分隔符，因此对于某些行，两列之间没有空格，如下所示：

从网页上它说这些被称为“卡片图像”。表格格式的信息描述如下：

目录数据文件由 80 个字符的卡片图像组成，每条谱线一张卡片图像。每张卡片图像的格式为：FREQ, ERR, LGINT, DR, ELO, GUP, TAG, QNFMT, QN', QN" (F13.4,F8.4, F8.4, I2,F10.4, I3, I7、I4、6I2、6I2)

我真的很想要一种只使用上面给出的格式说明符的方法。我唯一发现的是 Numpy 的 genfromtxt 函数。但是，以下不起作用。

任何人都知道我如何使用给定的每一列的格式规范将该表读入 Python 中？

有没有办法将 numpy 屏蔽数组中的屏蔽值替换为 null 或 None 值？这是我尝试过但不起作用的方法。

具体来说

生产

当我将数据帧输出为 csv 文件时，我想用 '' 替换掩码值 '--'。一种解决方法是将输出的 csv 文件读回 python 并将 '--' 替换为 ''，但这是一个可怕的解决方案。必须有更好的解决方案。我不希望在 csv 文件中将掩码值打印为“--”。