funhist - 云端在线

程序：

您的姓名

funhist - 创建一列的一维直方图（来自 FITS 二进制表或原始事件文件）
或图像

概要

有趣的 [-n抓-w抓-T] [列] [[lo:hi:]bins]

配置

-n # 按每个 bin 的宽度标准化 bin 值
-w # 在 arg3 中指定 bin 宽度而不是 bin 数量
-T # 以 rdb/starbase 格式输出（制表符分隔符）

商品描述

有趣的 从 FITS 扩展的指定列创建一维直方图
FITS 文件（或来自非 FITS 原始事件文件）的二进制表，或来自 FITS 图像或
数组，并将该直方图写入 ASCII 表。 或者，程序可以执行
图像轴之一的一维投影。

程序的第一个参数是必需的，并指定 Funtools 文件：FITS 表
或图像、原始事件文件或数组。 如果指定了“stdin”，则从
标准输入。 使用 Funtools Bracket Notation 指定 FITS 扩展和过滤器。

对于表，第二个参数也是必需的。 它指定要使用的列
生成直方图。 如果数据文件是图像（或数组）类型，则列是
可选：如果指定了“x”（或“X”）、“y”（或“Y”），则在其上执行投影
分别为 x (dim1) 或 y (dim2) 轴。 （也就是说，这个投影将给出相同的
结果作为对包含等效 x,y 事件行的表执行的直方图。）如果
没有指定列名或为图像指定了“xy”（或“XY”），然后是直方图
对包含在图像像素中的值执行。

后面的参数是可选的，它指定在创建时使用的 bin 数量
直方图，如果需要，还有 bin 值的范围。 对于图像和表格直方图，
范围应指定最小和最大数据值。 对于 x 和 y 上的图像直方图
轴，范围应指定最小和最大图像 bin 值。 如果这个论点是
省略，表的输出箱数根据 TLMIN/TLMAX 计算
标题值（如果这些存在于指定列的表 FITS 标题中）或通过
通过数据计算最小值和最大值。 对于图像，数量
输出箱是根据 DATAMIN/DATAMAX 标头值计算的，或者通过
通过数据计算最小值和最大值。 （请注意，后一种计算可能
如果图像无法放入内存，则失败。）如果数据是浮点数（表或
image) 并且没有指定 bin 的数量，使用任意默认值 128。

对于二进制表处理， -w (bin width) 开关可用于指定宽度
每个 bin 而不是 bin 的数量。 因此：

funhist test.ev 阶段 1:100:5

表示在直方图中使用了 5 个宽度为 20 的 bin，而：

funhist -w test.ev 阶段 1:100:5

表示在直方图中使用了 20 个宽度为 5 的 bin。

数据被分成指定数量的 bin 和生成的一维直方图
（或投影）以 ASCII 表格式输出。 对于表，输出显示
数据的 low_edge（包含）和 hi_edge（不包含）值。 例如，一个 15 行
包含值范围从 -7.5 到 7.5 的“pha”列的表可以这样处理：

[sh] funhist test.ev pha
# 数据文件：/home/eric/data/test.ev
# 列：pha
# 最小、最大、分箱：-7.5 7.5 15

bin 值 lo_edge hi_edge
------ --------- --------------------- -------------- -------
1年 22 -7.50000000 -6.50000000
2年 21 -6.50000000 -5.50000000
3年 20 -5.50000000 -4.50000000
4年 19 -4.50000000 -3.50000000
5年 18 -3.50000000 -2.50000000
6年 17 -2.50000000 -1.50000000
7年 16 -1.50000000 -0.50000000
8 -30
9 16 0.50000000 1.50000000
10 17 1.50000000 2.50000000
11 18 2.50000000 3.50000000
12 19 3.50000000 4.50000000
13 20 4.50000000 5.50000000
14 21 5.50000000 6.50000000
15 22 6.50000000 7.50000000

[sh] funhist test.ev 阶段 1:6
# 数据文件：/home/eric/data/test.ev
# 列：pha
# 最小、最大、分箱：0.5 6.5 6

bin 值 lo_edge hi_edge
------ --------- --------------------- -------------- -------
1 16 0.50000000 1.50000000
2 17 1.50000000 2.50000000
3 18 2.50000000 3.50000000
4 19 3.50000000 4.50000000
5 20 4.50000000 5.50000000
6 21 5.50000000 6.50000000

[sh] funhist test.ev 阶段 1:6:3
# 数据文件：/home/eric/data/test.ev
# 列：pha
# 最小、最大、分箱：0.5 6.5 3

bin 值 lo_edge hi_edge
------ --------- --------------------- -------------- -------
1 33 0.50000000 2.50000000
2 37 2.50000000 4.50000000
3 41 4.50000000 6.50000000

对于表直方图， -n(normalize) 开关可用于通过以下方式对 bin 值进行归一化
bin 的宽度（即 hi_edge-lo_edge）：

[sh] funhist -n test.ev pha 1:6:3
# 数据文件：test.ev
# 列：pha
# 最小、最大、分箱：0.5 6.5 3
# 应用宽度归一化 (val/(hi_edge-lo_edge))

bin 值 lo_edge hi_edge
—————————————————————————————— -------------------
1 16.50000000 0.50000000 2.50000000
2 6.16666667 2.50000000 4.50000000
3 4.10000000 4.50000000 6.50000000

例如，这可以用于生成具有单位为的值的光曲线
计数/秒而不是计数。

对于图像直方图，输出显示图像的低值和高值（包括两者）
用于生成直方图。 例如，在下面的例子中，184 个像素有一个
值为 1、31 的值为 2，而只有 2 的值为 3,4,5,6、7、XNUMX、XNUMX 或 XNUMX：

[sh] funhist test.fits
# 数据文件：/home/eric/data/test.fits
# 最小、最大、分箱：1 7 7

bin 值 lo_val hi_val
—————————————————————————————— -------------------
1 184.00000000 1.00000000 1.00000000
2 31.00000000 2.00000000 2.00000000
3 2.00000000 3.00000000 3.00000000
4 2.00000000 4.00000000 4.00000000
5 2.00000000 5.00000000 5.00000000
6 2.00000000 6.00000000 6.00000000
7 2.00000000 7.00000000 7.00000000

对于图像的轴投影，输出显示低图像箱和高图像箱（均
包括）用于生成投影。 例如，在下面的例子中，21
counts 的 X bin 值为 2，依此类推：

[sh] funhist test.fits x 2:7
# 数据文件：/home/eric/data/test.fits
# 列：X
# 最小、最大、分箱：2 7 6

bin 值 lo_bin hi_bin
—————————————————————————————— -------------------
1 21.00000000 2.00000000 2.00000000
2 20.00000000 3.00000000 3.00000000
3 19.00000000 4.00000000 4.00000000
4 18.00000000 5.00000000 5.00000000
5 17.00000000 6.00000000 6.00000000
6 16.00000000 7.00000000 7.00000000

[sh] funhist test.fits x 2:7:2
# 数据文件：/home/eric/data/test.fits
# 列：X
# 最小、最大、分箱：2 7 2

bin 值 lo_bin hi_bin
—————————————————————————————— -------------------
1 60.00000000 2.00000000 4.00000000
2 51.00000000 5.00000000 7.00000000

您可以使用 gnuplot 或其他绘图程序来绘制结果，使用脚本，例如
如：

#!/ bin / sh的
sed -e '1,/---- .*/d
/^$/,$d' 获得 \
awk '\
BEGIN{print "set nokey; set title \"funhist\"; set xlabel \"bin\"; set ylabel \"counts\"; plot \"-\" with boxes"} \
{print $3, $2, $4-$3}'获取\
gnuplot -persist - 1>/dev/null 2>&1

脚本中提供了类似的绘图命令 有趣的情节:

funhist test.ev pha ... 玩 funhist.plot gnuplot

使用 onworks.net 服务在线使用 funhist