问题标签 [fits]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
1 回答
1577 浏览

python - 使用标题中的坐标拼接图像数组(FITS 文件)

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

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

  1. 如何调用标题来识别银河纬度?
  2. 以这样的方式拼接数组,使其仅包含纬度范围的值,其他所有值都为 0。
0 投票
1 回答
158 浏览

python - 无法从 SDSS 中读取 WCS 共同添加的 Stripe 82 适合图像

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

我得到错误:

头文件是:

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

0 投票
2 回答
565 浏览

python - 当其中一列是数组时,如何更改 numpy recarray 的 dtype?

之前的帖子dtype中,我已经看到recarray可以使用astype. 但是,我无法使用recarray其中一列中有一个数组的 a 来做到这一点。

recarray来自 FITS 文件记录:

我需要将AC列从转换intfloat所以我尝试过:

而且,虽然看起来dtype确实发生了变化

查看AC中的数据表明它们仍然是整数值。例如,sum计算达到int16变量的限制(所有AC列值都是正数):

有没有办法有效地改变数组数据类型?

0 投票
1 回答
1178 浏览

python - 使用 astropy.io.fits 编写适合文件

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

这是我的代码示例:

我得到的错误是

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

任何帮助,将不胜感激。

0 投票
1 回答
343 浏览

python - Python 和 fit 模块(操作系统接口):“IOError: [Errno 2] No such file or directory”

似乎对于我尝试使用的脚本,我返回了一个有点令人困惑的错误,也就是说,在尝试从旧文件扩展名创建新文件时,它表示它正在尝试创建的文件不存在?谁能解释为什么会发生这种情况以及可能的解决方法?

这是我正在使用的脚本:

运行时:

我收到以下错误:

完整追溯:

0 投票
1 回答
3492 浏览

python - .fits 图像的径向轮廓

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

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

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

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

在此处输入图像描述

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

在此处输入图像描述

0 投票
2 回答
946 浏览

python - 使用另一个图像文件的指定尺寸在astropy中创建一个新的图像文件

我在 astropy 包中遇到适合文件操作的问题,我需要一些帮助。

我基本上想拍摄一张适合文件格式的图像,并创建一个我需要开始输入校正因子的新文件和一个新图像,然后可以将其与校正因子和原始图像一起使用以生成校正图像。这些中的每一个都将具有相同的尺寸。

从此开始:

谁能给我我缺少的明显知识来执行此操作...最后两行只是为了概述缺少的内容。我已经包括了??也许表明我需要其他东西,也许是 fit.writeto() 或类似的东西......

0 投票
1 回答
1250 浏览

python - astropy.io.fits 从具有多个 HDU 的大型拟合文件中读取行

我有一个大约 50GB 的 fit 文件,其中包含多个 HDU,它们都具有相同的格式:一个 (1E5 x 1E6) 数组,其中包含 1E5 个对象和 1E6 个时间戳。HDU 描述了不同的物理属性,例如通量、RA、DEC 等。我只想从每个 HDU 中读取 5 个对象(即 (5 x 1E6) 阵列)。

python 2.7,astropy 1.0.3,linux x86_64

到目前为止,我尝试了很多我发现的建议,但没有任何效果。我最好的方法仍然是:

此代码适用于高达 ~20 GB 的文件,运行良好且快速,但对于较大的文件会耗尽内存(较大的文件仅包含更多对象,而不包含更多时间戳)。我不明白为什么 - astropy.io.fits 本质上使用 mmap 并且据我了解应该只将 (5x1E6) 数组加载到内存中?因此独立于文件大小,我想读出的内容总是相同的大小。

编辑 - 这是错误消息:

编辑 2:谢谢,我现在包含了建议,它使我能够处理高达 50GB 的适合文件。新代码:

没有引起任何改变。

确实不是默认的,需要手动设置。

etc 现在允许我读取 50GB 而不是 20GB 的文件

新问题:任何大于 50GB 的内容仍会导致相同的内存错误 - 但是,现在直接在第一行。

0 投票
0 回答
526 浏览

java - 在 Java 应用程序中打开大文件时出现 Java Map Failed 错误

尝试在 Java 应用程序中打开一个大文件(大约 128KB)会导致 Map Failed 错误:

我尝试在此站点上搜索答案,但我能找到的唯一类似案例是将代码加载到框架中。以这个为例,Java Large File Uploads throws java.io.IOException: Map failed。我正在使用的应用程序是一个 .jar 文件(更多信息在这里 - http://www.star.bris.ac.uk/~mbt/topcat/),我试图加载的文件是 .jar 文件。适合类型。

你能给我一些关于这个问题的建议吗?我正在使用 JRE 构建 1.8.0_73

更新:该错误似乎仅在我增加最大堆大小时出现。保持默认设置时一切正常。

0 投票
2 回答
1603 浏览

python - 我该怎么做才能让我的程序不为 .fits 文件中不存在的 Header 抛出 KeyError?

我正在制作一个程序,在计算机中搜索 .fts 和 .fits 文件,在其中打开文件并检索与标题中特定关键字相对应的信息,并将文件重命名为该关键字。

我遇到了一个问题,我一直收到 KeyError,因为在文件中找不到我正在搜索的标题关键字。有没有解决的办法?即使文件中不存在该关键字,我也希望能够搜索各种关键字并执行某些操作。

这是代码:

这是我得到的错误: