问题标签 [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 回答
782 浏览

python - 使用 pyfits 一次读取 FITS 二进制表一行

我有一个包含二进制表的 60GB FITS 文件。我想一次读取(并处理)该表一行/条目/行/块*。

(*我不确定正确的命名法)

我正在使用 pyfits,我想做的归结为:

该变量entry具有类型<class 'pyfits.fitsrec.FITS_record'>,并且长度等于二进制表中的列数。然而,似乎发生的是整个二进制表在这一行被读入内存:entry = hdulist[1].data[i].

我已经查看了 pyfits 文档,但我找不到任何方法似乎可以从表条目的表条目(或一次少量条目)上从二进制表扩展读取数据。我不想从表中选择某些条目,只需按顺序扫描它们即可。

我想我的问题是:

0) 这一步发生了hdulist[1].data[i]什么?为什么所有内容都被读入内存?(有没有办法解决这个问题?)

1)我错过了什么,pyfits 真的可以做我想做的吗?

2)那里有另一个python库吗?(即在 FITS 扩展中使用二进制表)

3)如果不是,我可以用不同的二进制(或其他压缩/非 ascii)格式(不是 FITS)重写数据并找到其他一些 python 库或模块来做我想做的事吗?

0 投票
1 回答
761 浏览

python - Using python to analyze a FITS image for a star count

I have several FITS images taken with my telescope and I need help with a python code to analyze the images and obtain a star count for each image. I don't have any code yet because I'm extremely new to python. Any help would be greatly appreciated.

0 投票
1 回答
228 浏览

python - 有条件地从 PyFITS 返回的 Numpy 数组中选择值

我在 pyfits 中打开了一个 FITS 文件。HEADER 文件使用 36个列标签读取XTENSION='BINTABLE',如、、、、、等。DIMENSION= 52989R x 36C'ZBEST''ZQUALITY''M_B''UB''PGAL'

现在,我必须从数据中选择“ZQUALITY”大于 2 且“PGAL”等于 3 的对象。然后我必须为符合上述条件的相应对象的“ZBEST”制作直方图。我还必须为这些对象绘制“M_B”与“UB”。

最后,我想将“ZBEST”分成三个切片(zbest < 0.5)、(0.5 < zbest < 1.0)、(zbest > 1.0),并分别绘制直方图和“M_B”与“UB”图。

我坚持选择服从这两个条件的数据。谁能告诉我如何从满足两个条件('ZQUALITY' > 2 & 'PGAL' == 3 )的数据中选择对象?我正在使用 like:data.field[('zquality' > 2) & ('pgal'==3)]但它不起作用。

0 投票
0 回答
826 浏览

matlab - MATLAB 的 fitiolib 中的错误:无法创建命名文件

我需要编写一组图像,这些图像是在 for 循环中计算的。如果我使用 定义路径sprintf,则会收到错误消息

您对如何避免这种情况有什么建议吗?这是计算图像的代码部分:

0 投票
0 回答
174 浏览

python - 如何将 .map 文件转换为 FITS 文件?

我想知道是否有人知道如何将 .map 文件转换为 .fits 文件,以便我可以在 Python 中读取 .fits 文件?

0 投票
1 回答
77 浏览

linux - 在 csh 中为 FITS 文件更改 BUNIT

我正在.csh 中编写代码,并且正在尝试将FITS文件的 bunit 标头从 K (kelvin) 更改为 km/s。我怎样才能做到这一点?

我知道在Python中我会使用 new_fitsfile.header['BUNIT']='km/s',但这在当前的 .csh 代码中不起作用,也不能将其切换为 Python 代码。

0 投票
4 回答
162 浏览

python - 按 FITS 文件中多列中的值聚合表行

我有一个包含许多列的 FITS 文件。该文件的简化示例的一部分如下所示:

如果您注意到这里, A 列和 C列的前四行相同,但B列从 1 到 4 不等。那么A 和 C列的接下来 5 个值相同,但B列从 40 到 44 不等。

我想做的是,编写一个程序来创建这样的文件:

也就是说,B 列应该包含 A 列和 C 列相同的值的数量!

我想知道如何在 Python 中做到这一点。它不一定需要处理一个FITS文件,如果FITS文件中有一些不能使用的例程,我也可以将其转换为ASCII格式。

到目前为止我已经尝试过:

我碰到了一个调用的例程Collections,它有一个调用的子例程Counter,它计算列表中相等的值的数量并返回它们。

我试过了:

但这只会给我A列中相等的值。有人可以告诉我如何使用此例程将这些值也与C列进行比较吗?

0 投票
1 回答
1004 浏览

fits - 使用正确的 wcs 转换使用 aplpy 绘制适合图像

我想用 aplpy 绘制一个 fit 文件,但出现以下错误:

如果我删除标题关键字 CTYPE3、CD3_3、WAT3_001,则不再出现错误,但在绘制 fit 文件时,我发现 WCS 转换没有正确地使用 RA 和 DEC 坐标在整个地方完成。

这是标题:

0 投票
2 回答
1764 浏览

linux - IRAF 中的“像素超出范围”错误

我正在尝试使用 SNOOPY(点扩散函数拟合器)和 IRAF 对拟合图像进行 PSF 拟合。我可以使用 imexam 很好地打开此图像,但是当我选择一个点(星形或其他)时出现错误:

看来我所看到的和 IRAF 看到的(幕后)并不相同。好像有某种坐标偏移或什么的。

如何解决这个问题?

[科学Linux 6、16bit、IRAF v2.16]

0 投票
1 回答
3828 浏览

python - 使用python在fits文件中查找像素的物理坐标

我想从 python 脚本中获取给定像素的物理天空坐标。我想使用 astropy 的 WCS,但我会在 python 中做任何事情。

我已经尝试了这两个代码片段。

问题是我第一次尝试使用 WCS 时出错,结果只是我输入的像素值。