Fundisp - 云端在线

程序：

您的姓名

Fundisp - 在 Funtools 数据文件中显示数据

概要

基金管理公司 [-f 格式] [-l] [-n] [-T] [列数bitpix=n]

配置

-f # 显示的格式字符串
-l # 将图像显示为包含 X、Y、VAL 列的列表
-n # 不输出标题
-F [c] # 使用指定字符作为列分隔符 (def: space)
-T # 以 rdb/starbase 格式输出（制表符分隔符）

商品描述

基金管理公司 在指定的 FITS 扩展和/或 FITS 的图像部分中显示数据
文件，或在非 FITS 数组或原始事件文件的部分中。

程序的第一个参数指定 FITS 输入图像、数组或原始事件文件
显示。 如果指定了“stdin”，则从标准输入读取数据。 使用 Funtools
用于指定 FITS 扩展、图像部分和过滤器的括号表示法。

如果显示的数据是列（在 FITS 二进制表或原始事件中）
文件），则列出各个行。 可以使用括号表示法添加过滤器。 因此：

[sh] fundisp "test.ev[time-(int)time>.15]"
XY PHA PI 时间 DX DY
------- ------- ------- --------- ---------------- ---- ————————————
10 8 10 8 17.1600 8.50 10.50
9 9 9 9 17.1600 9.50 9.50
10 9 10 9 18.1600 9.50 10.50
10 9 10 9 18.1700 9.50 10.50
8 10 8 10 17.1600 10.50 8.50
9 10 9 10 18.1600 10.50 9.50
9 10 9 10 18.1700 10.50 9.50
10 10 10 10 19.1600 10.50 10.50
10 10 10 10 19.1700 10.50 10.50
10 10 10 10 19.1800 10.50 10.50

[注意：FITS 二进制表测试文件 test.ev，以及 FITS 图像 test.fits，是
包含在 funtools funtest 目录中。]

当使用表显示时 基金管理公司，第二个可选参数可用于
指定要显示的列。 例如：

[sh] fundisp "test.ev[time-(int)time>=.99]" "xy 时间"
XY时间
——————————————————
5 -6 40.99000000
4 -5 59.99000000
-1
-2
-3
-4
-5
-6
-7

特别栏目 $地区 可以指定显示每一行的区域id：

[sh $] fundisp "test.ev[time-(int)time>=.99&&&annulus(0 0 0 10 n=3)]" 'xy time $REGION'
XY 时间区域
—————————————————————————————————
5 -6 40.99000000 3
4 -5 59.99000000 2
-1 0 154.99000000 1
-2 1 168.99000000 1
-3 2 183.99000000 2
-4 3 199.99000000 2
-5 4 216.99000000 2
-6 5 234.99000000 3
-7 6 253.99000000 3

这里只有具有适当小数时间且其位置也在以下之一内的行
显示三个环。

可以使用列前的减号从显示中排除列：

[sh $] fundisp "test.ev[time-(int)time>=.99]" "-time"
XY PHA PI DX DY
-------- -------- -------- ---------- ----------- ----- ------
5 -6 5 -6 5.50 -6.50
4 -5 4 -5 4.50 -5.50
-1 0 -1 0 -1.50 0.50
-2 1 -2 1 -2.50 1.50
-3 2 -3 2 -3.50 2.50
-4 3 -4 3 -4.50 3.50
-5 4 -5 4 -5.50 4.50
-6 5 -6 5 -6.50 5.50
-7 6 -7 6 -7.50 6.50

显示除时间列之外的所有列。

特别栏目 $N 可以指定显示每一行的序号值。 因此，
继续前面的例子：

Fundisp "test.ev[time-(int)time>=.99]" '-time $n'
XY PHA PI DX DY N
------- -------- -------- ---------- ----------- ------ ----- ----------
5 -6 5 -6 5.50 -6.50 337
4 -5 4 -5 4.50 -5.50 356
-1 0 -1 0 -1.50 0.50 451
-2 1 -2 1 -2.50 1.50 465
-3 2 -3 2 -3.50 2.50 480
-4 3 -4 3 -4.50 3.50 496
-5 4 -5 4 -5.50 4.50 513
-6 5 -6 5 -6.50 5.50 531
-7 6 -7 6 -7.50 6.50 550

请注意，列规范用单引号括起来以保护 '$n' 从开始
被壳膨胀。

一般来说，激活和去激活列的规则是：

· 如果只指定了排除列，那么除排除列之外的所有列都将
被激活。

· 如果只指定了包含列，则只激活指定的列。

· 如果包含和排除列的混合被指定，那么除了排除列之外的所有列
列将处于活动状态； 最后一种情况是模棱两可的，规则是任意的。

除了明确指定列名外，特殊符号 + 和 - 可
用于激活和取消激活 所有 列。 如果您想激活
$REGION 列以及所有其他列。 根据规则，语法“$REGION”
仅激活区域列并停用其余部分。 使用“+ $REGION”激活
所有列以及区域列。

如果显示的数据是图像数据（无论是在 FITS 主图像中，还是在 FITS 图像中
扩展名或数组文件），则生成 mxn 像素显示，其中 m 和 n 是
图像的尺寸。 默认情况下，像素值使用相同的数据类型显示
就像在文件中一样。 但是，对于整数数据，其中 BSCALE 和 BZERO 标头参数是
目前，数据显示为浮点数。 在任何一种情况下，显示数据类型都可以是
使用以下形式的可选第二个参数覆盖：

位像素=n

其中 n 为 8,16,32、32、64、-XNUMX、-XNUMX，分别表示 unsigned char、short、int、float 和 double。

当然，跑步 基金管理公司 除了最小的图像之外的任何东西通常都会导致显示
它的大小使它不可读。 因此，可以使用括号表示法（见下文）来
在生成显示之前对图像应用部分和/或块。 例如：

[sh] fundisp "test.fits[2:6,2:7]" bitpix=-32
2 3 4 5 6
---------- ---------- ---------- ---------- ----------
2：3.00 4.00 5.00 6.00 7.00
3：4.00 5.00 6.00 7.00 8.00
4：5.00 6.00 7.00 8.00 9.00
5：6.00 7.00 8.00 9.00 10.00
6：7.00 8.00 9.00 10.00 11.00
7：8.00 9.00 10.00 11.00 12.00

请注意，可以通过简单地将 FITS 二进制表显示为图像
表通过 有趣的形象 第一：

[sh] ./funimage test.ev stdout 获得fundisp "stdin[2:6,2:7]" bitpix=8
2 3 4 5 6
-------- ------- ------- -------- --------
2：3 4 5 6 7
3：4 5 6 7 8
4：5 6 7 8 9
5：6 7 8 9 10
6：7 8 9 10 11
7：8 9 10 11 12

如果 -l (list) 开关使用，然后图像显示为包含列表
列：X、Y、VAL。 例如：

fundisp -l "test1.fits[2:6,2:7]" bitpix=-32
XY值
---------- ---------- ------------
2 2 6.00
3 2 1.00
4 2 1.00
5 2 1.00
6 2 1.00
2 3 1.00
3 3 5.00
4 3 1.00
5 3 1.00
6 3 1.00
2 4 1.00
3 4 1.00
4 4 4.00
5 4 1.00
6 4 1.00
2 5 1.00
3 5 1.00
4 5 1.00
5 5 3.00
6 5 1.00
2 6 1.00
3 6 1.00
4 6 1.00
5 6 1.00
6 6 2.00
2 7 1.00
3 7 1.00
4 7 1.00
5 7 1.00
6 7 1.00

如果 -n (nohead) 开关被使用，然后没有表头输出。 这个有用，
例如，当fundisp 输出被定向到gnuplot 时。

基金管理公司 程序使用一组默认的显示格式：

数据类型 TFORM 格式
———————————————
双 D "%21.8f"
浮动 E "%11.2f"
int J "%10d"
短我“%8d”
字节 B“%6d”
字符串 A "%12.12s"
位 X“%8x”
逻辑 L“%1x”

因此，1 双和 2 短裤的默认显示给出：

[sh]fundisp snr.ev“时间xy”

时间 XY
--------- -------- --------
79494546.56818075 546 201
79488769.94469175 548 201
...

您可以更改单个列或给定的所有列的显示格式
数据类型通过 -f 开关。 伴随 -f 的格式字符串是一个空格-
关键字=格式值的分隔列表。 关键字值可以是列名
（在这种情况下，关联格式仅与该列有关）或 FITS 表 TFORM
说明符（在这种情况下，格式适用于具有该数据类型的所有列）。 为了
例如，您可以像这样更改所有列的双精度和短格式：

[sh] fundisp -f "D=%22.11f I=%3d" snr.ev "time xy"

时间 XY
--------------- --- ---
79494546.56818075478 546 201
79488769.94469174743 548 201
...

或者，您可以像这样更改时间和 x 列的格式：

[sh] fundisp -f "time=%22.11fx=%3d" snr.ev "time xy"

时间 XY
-------------- --- --------
79494546.56818075478 546 201
79488769.94469174743 548 201
...

请注意，如果列与 TFORM 中的一个具有相同的名称，则存在潜在冲突
说明符。 在上面的例子中，表中的“X”列与
X（位）数据类型。 为了解决这个冲突，格式字符串被处理成
使用区分大小写的比较首先检查 TFORM 数据类型说明符。 如果
指定的格式值不是大写的 TFORM 值，则不区分大小写的检查是
对列名进行。 这意味着，在上面的例子中，“X=%3d”将指
X（位）数据类型，而 "x=%3d" 将引用 X 列：

[sh] fundisp -f "X=%3d" snr.ev "xy"

XY
—————————
546 201
548 201
...

[sh] fundisp -f "x=%3d" snr.ev "xy"

XY
--- --------
546 201
548 201
...

因此，作为一项规则，最好始终以小写和 TFORM 指定列名
大写的数据类型。

-f [格式] 将更改单次执行 fundisp 的格式。 你也可以使用
此 有趣的格式 环境变量更改所有fundisp 调用的格式。
此环境变量值的格式与用于 -f
转变。 在个别情况下，可以通过使用 -f [格式]
开关。

警告：另请注意，匹配格式是用户的责任
列数据类型的说明符正确。 还要注意的是，为了保持视觉
名称和列之间的对齐，列名称将被截断（在左侧），如果
格式宽度小于名称的长度。 但是，不执行截断
如果输出是 RDB 格式（使用 -T 开关）。

[支持但不推荐使用旧式格式字符串。 它由以空格分隔的
所有数据类型的 C 格式语句，按以下顺序指定：

double float int 短字节字符串位。

这个列表的顺序是基于人们通常想要的假设
更改浮动格式。

如果输入“-”而不是给定数据类型的格式语句，则默认格式
用来。 此外，格式字符串可以在不指定所有格式的情况下终止，并且
默认值将用于列表的其余部分。 请注意，您必须提供最小字段
宽度，即“%6d”和“%-6d”是合法的，“%d”是不合法的。

通过使用 -f [format]，您可以像这样更改双格式和短格式：

[sh]fundisp -f "22.11f - - 3d" snr.ev "time xy"

时间 XY
--------------- --- ---
79494546.56818075478 546 201
79488769.94469174743 548 201
...

注意：此格式已弃用，将在未来版本中删除。]

-F[c] 开关可用于指定（单字符）列分隔符（其中
默认为空格）。 请注意，列格式几乎肯定也会向
将各个列填充到所需的宽度。 这些可以用一个程序删除，例如
sed，以生成未对齐的列为代价。 例如：

Fundisp -F',' snr.ev'[cir 512 512 .1]'
X、Y、PHA、PI、时间、DX、DY
512，512，6，7，79493997.45854475，578，574
512，512，8，9，79494575.58943175，579，573
512，512，5，6，79493631.03866175，578，575
512，512，5，5，79493290.86521725，578，575
512，512，8，9，79493432.00990875，579，573

fundisp -F',' snr.ev'[cir 512 512 .1]' 开启 sed 's/ *, */,/g'
X、Y、PHA、PI、时间、DX、DY
512,512,6,7,79493997.45854475,578,574
512,512,8,9,79494575.58943175,579,573
512,512,5,6,79493631.03866175,578,575
512,512,5,5,79493290.86521725,578,575
512,512,8,9,79493432.00990875,579,573

Fundisp -f "x=%3d y=%3d pi=%1d pha=%1d 时间=%20.11f dx=%3d dy=%3d" -F',' snr.ev'[cir 512
512 .1]' 玩 sed 's/ *, */,/g'
X,Y,A,I,TIME,DX,DY ---,---,-,-,--------------------,---, ---
512,512,6,7,79493997.45854474604,578,574 512,512,8,9,79494575.58943174779,579,573
512,512,5,6,79493631.03866174817,578,575 512,512,5,5,79493290.86521725357,578,575
512,512,8,9,79493432.00990875065,579,573

如果 -T (rdb table) 开关被使用，输出将符合starbase/rdb 数据库
格式：制表符将插入列之间而不是空格。 这种格式不是
显示图像像素时可用（与 -l 转变）。

最后，请注意 基金管理公司 可用于从辅助表创建列过滤器
在 FITS 文件中。 例如，下面的 shell 代码将生成一个良好的时间间隔
(GTI) 过滤器，用于包含标准 GTI 扩展名的 X 射线数据文件：

#!/ bin / sh的
sed '1,/---- .*/d
/^$/,$d' 玩 awk 'tot>0{printf "泰国"};{printftime="$1":"$2; 总++}'

如果此脚本放置在名为“mkgti”的文件中，则可以在诸如以下命令中使用它：

fundisp foo.fits"[GTI]" 获得 mkgti > gti.filter

然后可以在各种 funtools 程序中使用生成的过滤器文件：

funcnts foo.fits"[@gti.filter]" ...

仅处理良好时间间隔内的事件。

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