问题标签 [pyfits]

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 投票
2 回答
888 浏览

python - 读取许多文件并创建动画python

我编写了从适合文件创建动画(电影)的脚本。一个文件大小为 2.8 MB,第一个文件大小为 2.8 MB。文件数为 9000。这是代码

我认为我的问题是当我创建数组 img[]...我有 8 GB RAM 并且当 RAM 已满时我的操作系统终止 python 脚本。

我的问题是:如何读取 9000 个文件并创建动画?是否可以创建一些缓冲区或一些并行处理?

有什么建议吗?

0 投票
2 回答
1045 浏览

python - 如何将两个表与pyfits合并?

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

错误是:

0 投票
1 回答
46 浏览

python - 列长操作

首先,抱歉,如果我使用了错误的语言来解释我正在操作的内容,我对 python 还很陌生,而且还远远不够了解它。

我目前正在尝试对一列数据的长度进行操作,但是遇到了一些问题。这些列来自第 7 个 sdss 数据转储中的 .fit 文件。当我运行代码时,根据边界条件打印 x1、x2、x3 的每个值。

从这一点上,我获得了大约 70000、90000 和 30000 的 x1、x2、x3 的值,总共大约 190000。我想包括一条线来计算它们占总值的百分比但是我似乎不能让它工作。

试过了。

这将返回不是对长度进行操作的元素,但是当我尝试将长度定义为一个值时,它返回的答案为 0。

任何帮助将非常感激。

0 投票
0 回答
102 浏览

python - 更新 pyfits bin 表数据

我正在尝试使用 pyfits 更新现有的 fit 表。它适用于表格的某些列,不幸的是不适用于第一列。这是列定义:

以及更新它的简单代码片段:

编写 TEL_ID(以及此处未显示的其他内容)有效,而 EVENT_ID 则无效。我已经尝试过不同的格式(np.int32,int)但总是一样的......

返回 numpy.uint32 (对于未修改的文件)

谢谢你的帮助

编辑:

如果我更改 'EVENT_ID' 的定义,省略 'bscale' 和 'bzero' 值的更新工作。所以似乎无符号整数有问题。

0 投票
2 回答
1603 浏览

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

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

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

这是代码:

这是我得到的错误:

0 投票
2 回答
898 浏览

python - 通过 PyFITS 裁剪 FITS 数据立方体图像

我一直在寻找解决我的问题的方法,但找不到。我有一个 FITS 数据立方体,我需要通过 PyFITS 对其进行裁剪。当我用我的脚本来做这件事时,最后我会得到一个 2-D FITS 图像!第一个维度是能量,第二个和第三个维度分别是经度和纬度。

我的脚本如下:

0 投票
1 回答
504 浏览

python - 保存 2 个带有 3 个标题 astropy 或 pyfits 的文件

我正在尝试制作一个包含两个图像的适合文件,并且我想给它们两个标题。他们两个都应该有一个基本的标题。到目前为止我发现的是这个,但我不知道如何正确添加标题。变量是headermain、header1、header2、data1、data2。我希望它具有这种结构。

我目前有这个,但我找不到如何添加标题。

0 投票
1 回答
1356 浏览

astropy - 将标头复制到新文件中

我有这个脚本可以创建一个具有多个扩展名的文件,但我想将旧文件的标题添加到扩展名中。

现在我试过了:

但两者都给我错误说

ValueError: header 必须是 Header 对象

我怎样才能解决这个问题?

http://pastebin.com/JXki7EPV

0 投票
2 回答
1516 浏览

python - WCS作为使用astropy加载的数据立方体切片的matplotlib投影?

我有一个名为“my_cube.fits”的带有 WCS 的 FITS 文件。该文件在轴 1 和 2(X 和 Y)上具有空间信息,在轴 3 (Z) 上具有光谱信息。当我使用astropy.io.fits加载它时,光谱轴为 0,空间轴为 1 和 2。文件加载如下:

我一直在使用 matplotlib 来显示数据,我想知道如何使用它的 WCS 显示我的数据立方体的单个光谱框架。比方说:

如果我这样做,我有:

0 投票
1 回答
1027 浏览

python - 更改适合图像的投影

这是我在这里的第一个问题,所以我会尽力解释清楚。

我正在处理一些全天适合的天文图像,它们位于银河坐标系中(在 AITOF 投影中)。我想将这些拟合图像转换为 RA、Dec 坐标,以便图像以正确的方式旋转(即银河平面在地图中心将不再水平,而是有点扭曲)。

有谁知道该怎么做?我试图以一种不太优雅的方式做到这一点,我在这里解释:

我打开 fit 文件,以及图像所在的 hdu 扩展名。我阅读了标题的适当关键字我创建了一个具有相同形状的数组,但是每个元素现在都是几个值,这将是每个像素的坐标。然后我将每个元素转换为我想要的新坐标(使用 astropy.coordinates 包),最后我必须移动每个元素以便在新坐标中对它们进行排序(首先它们在银河坐标中排序,现在它们必须按天体分类)。

到目前为止,我的代码非常慢,我相信一定有更好的方法来做我打算做的事情:

在那之后,我不知道如何以正确的方式对数组进行重新排序......

非常感谢你们!