问题标签 [photutils]

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 回答
147 浏览

python - Photutils/IRAF Starfinder:给出总数

我一直在使用来自 photutils 包的 IRAF Starfinder 来识别来自 fit 图像的来源,或者更确切地说来自从 fit 图像中提取的 2D 数组。来源的识别似乎工作得很好,我得到了一个包含所有常用参数的天文表。但是,我希望简单地获取源的总计数(总计数表示源半径内像素值的总和)。我意识到这里也没有定义半径,但是 FWHM 的某些因素在这里就足够了。通量参数似乎是我想要的最接近的东西,但在文档中它将其描述为对象通量,计算为卷积图像中的峰值密度除以检测阈值,这不是我所追求的。是否有捷径可寻?我对这些软件包和这种分析形式很陌生,所以如果这是一个明显/愚蠢的问题,我深表歉意。感谢您的时间!

0 投票
1 回答
62 浏览

python - 从具有过饱和像素的图像中精确检测光源中心

我有一个 16 位图像,其主要对象是 6 个 LED,两个反射方块,图像的其余部分非常暗。这两个正方形与两个 LED 部分重叠。大约 50% 的可见 LED 过饱和(它们的中心),像素强度为 65532。在中心附近,靠近边缘,可以看到溢出的光,我认为这可以用来精确确定每个 LED 的质心引领。我正在尝试通过使用 astropy photutils 包(https://photutils.readthedocs.io/en/stable/index.html),特别是它的方法:source detection()、centroids()、aperture photometry... 两个主要任务是如何区分与提到的反射正方形重叠的 LED 以及如何精确确定每个 LED 的中心。我试图增加阈值,这部分有效,正方形不再可见,但我也失去了 LED 的边缘,这导致它只有过饱和的中心。
有人熟悉这个问题吗?

0 投票
0 回答
34 浏览

python - Find_Peaks:图像数据的无效形状 (4951,)。可能的解决方案?

大家好,当我尝试通过 imshow 显示时遇到了这个错误:Invalid shape (4951,) for image data

该数据是我通过运行图像得到的像素“peak_vals”输出photutils.find_peaks()。原始形状是 (5820,)。我很确定这个错误是由于形状不规则而发生的,但我不确定是否可以将其重塑为正确的尺寸。

所以我的问题是:

  1. 如果有一种可能的重塑方法,它是什么?

  2. 如果不是,我如何在我正在使用的图像中找到连接的像素(或超像素)?我最初的方法是将 peak_vals 数据(像素强度数据)拟合到椭圆,并将数据/像素过滤到适合椭圆内的那些。2 次迭代导致数据重塑(从 (5820,) 到 (4951,)。

0 投票
1 回答
39 浏览

python - 如何为长轴和短轴的值列表绘制椭圆孔径?

我想在 UV FITS 文件上绘制不同区域的椭圆孔,用于长轴和短轴和 theta 的不同值,长轴与正 x 方向的角度。我有一组值来定义所需的光圈中心、相应的主轴和次轴以及 theta 值。我附上我写的代码。但是,在定义光圈时,它不接受数组作为长轴、短轴或 theta 的合格值。

0 投票
0 回答
63 浏览

python - 来自 photoutils 孔径光度计的更多详细信息

我在 SAOimageDS9 中手动绘制圆形区域,并使用 photutils python 包在给定坐标和半径的情况下对 DS9 的区域进行圆形光圈测光。

我能够获得总通量测量值,但想知道是否有办法使用孔径光度法获得更多数据值,例如孔径内的像素数和 DS9 的平均/中值/标准偏差。

如果有人知道一种直接从 DS9 区域分析中自动检索这些值的方法,那也将非常有帮助!

0 投票
2 回答
114 浏览

python - 如何将 QTable 保存到文件(.txt、.csv 等)

我正在尝试将源自椭圆等光度拟合的等光度列表的 QTable 对象保存到任何文件类型(.txt、.csv 等),让我只需将其加载回脚本而无需清理它或任何一种。

  • QTable 创建如下:
  • 我的尝试:

在这里使用 json 不起作用。出现以下错误:

这里的任何人都有过 photutils 数据处理或将 isoophote 保存到文件的经验吗?这真的应该是一件非常简单的事情,只是一种方式,这样我就不必每次想要使用等光度拟合的结果时都重新运行整个脚本。我的整个数据集包含 26 张图像,这意味着大约需要 2.5 小时的计算时间,而不会在两者之间进行保存。

提前致谢!

0 投票
0 回答
25 浏览

python - 使用 jupyter 笔记本时遇到阈值问题

您好,我正在尝试对我的 Fit 文件进行阈值检测,但是当我运行它时,我的内核死了。我调试它并停止我的内核以查看卡在哪里。它通常是我代码的第二行。我正在尝试做类似的事情:https ://photutils.readthedocs.io/en/latest/segmentation.html

我的 fit 文件是 24000x24000 像素可能是它太多并且阈值未检测到还是我的 sigma 值错误?image_data1 是它读取我的 fit 文件的方式。有什么建议吗?

0 投票
0 回答
23 浏览

astropy - 如何告诉`photoutils`仅绘制满足条件的孔径?

我正在按照文档中photutils示例来检测图像中的来源:

哪个生产

在此处输入图像描述

我添加了最后四行,目的是只在最亮的星星周围绘制光圈。但是,for循环不会改变图像。我明白为什么(它会多次绘制所有光圈,对于 mag < -2 的 4 颗恒星中的每颗都绘制一次),但是我如何更改它以仅为这些恒星绘制它们?

0 投票
1 回答
50 浏览

python - 如何为已安装模块中发生的迭代创建进度条

我的目标是为已安装模块内发生的迭代创建一个进度条。

为了在用户定义的函数中为迭代创建一个进度条,我将一个tqdm.notebook.tqdm_notebook对象作为可迭代对象传递:

哪个有效:

在此处输入图像描述

但是,当我尝试对已安装模块内的 for 循环执行相同操作时,此操作会失败。在 Astropy 的 Photutils 模块中,有for label in labels一行(此处),我可以传递标签对象。

可重现的示例(主要基于- 在安装 photutils: 后工作pip install photutils):

这有效:

试图传递tqdm.notebook.tqdm_notebook对象来创建进度条:

我得到一个AttributeError: 'int' object has no attribute '_comparable'. 完整追溯:

一种解决方法是修改 Photutils 并tqdm在其中使用(我在这个 fork上做过,它有效),但这似乎有点过头了,我希望有一种更简单的方法来做到这一点。