问题标签 [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.
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 库或模块来做我想做的事吗?
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.
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)]
但它不起作用。
matlab - MATLAB 的 fitiolib 中的错误:无法创建命名文件
我需要编写一组图像,这些图像是在 for 循环中计算的。如果我使用 定义路径sprintf
,则会收到错误消息
您对如何避免这种情况有什么建议吗?这是计算图像的代码部分:
python - 如何将 .map 文件转换为 FITS 文件?
我想知道是否有人知道如何将 .map 文件转换为 .fits 文件,以便我可以在 Python 中读取 .fits 文件?
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列进行比较吗?
fits - 使用正确的 wcs 转换使用 aplpy 绘制适合图像
我想用 aplpy 绘制一个 fit 文件,但出现以下错误:
如果我删除标题关键字 CTYPE3、CD3_3、WAT3_001,则不再出现错误,但在绘制 fit 文件时,我发现 WCS 转换没有正确地使用 RA 和 DEC 坐标在整个地方完成。
这是标题:
linux - IRAF 中的“像素超出范围”错误
我正在尝试使用 SNOOPY(点扩散函数拟合器)和 IRAF 对拟合图像进行 PSF 拟合。我可以使用 imexam 很好地打开此图像,但是当我选择一个点(星形或其他)时出现错误:
看来我所看到的和 IRAF 看到的(幕后)并不相同。好像有某种坐标偏移或什么的。
如何解决这个问题?
[科学Linux 6、16bit、IRAF v2.16]
python - 使用python在fits文件中查找像素的物理坐标
我想从 python 脚本中获取给定像素的物理天空坐标。我想使用 astropy 的 WCS,但我会在 python 中做任何事情。
我已经尝试了这两个代码片段。
问题是我第一次尝试使用 WCS 时出错,结果只是我输入的像素值。