funhist - En ligne dans le Cloud

PROGRAMME:

Nom

funhist - crée un histogramme 1D d'une colonne (à partir d'une table binaire FITS ou d'un fichier d'événement brut)
ou une image

SYNOPSIS

funhiste [-n⎪-w⎪-T] [colonne] [[lo:hi:]bins]

OPTIONS

-n # normalise la valeur du bac par la largeur de chaque bac
-w # spécifie la largeur du bac au lieu du nombre de bacs dans arg3
-T # sortie au format rdb/starbase (séparateurs de tabulation)

DESCRIPTION

funhiste crée un histogramme unidimensionnel à partir des colonnes spécifiées d'une extension FITS
table binaire d'un fichier FITS (ou d'un fichier d'événement brut non FITS), ou d'une image FITS ou
tableau et écrit cet histogramme sous forme de table ASCII. Alternativement, le programme peut effectuer
une projection 1D d'un des axes de l'image.

Le premier argument du programme est obligatoire, et spécifie le fichier Funtools : table FITS
ou image, fichier d'événement brut ou tableau. Si "stdin" est spécifié, les données sont lues à partir du
entrée standard. Utilisez Funtools Bracket Notation pour spécifier les extensions et les filtres FITS.

Pour une table, le deuxième argument est également requis. Il spécifie la colonne à utiliser dans
génération de l'histogramme. Si le fichier de données est de type image (ou tableau), la colonne est
facultatif : si "x" (ou "X"), "y" (ou "Y") est spécifié, alors une projection est effectuée sur
les axes x (dim1) ou y (dim2), respectivement. (C'est-à-dire que cette projection donnera le même
résultats sous la forme d'un histogramme exécuté sur une table contenant les lignes d'événements x,y équivalentes.) Si
aucun nom de colonne n'est spécifié ou "xy" (ou "XY") est spécifié pour l'image, puis un histogramme
est effectuée sur les valeurs contenues dans les pixels de l'image.

L'argument qui suit est facultatif et spécifie le nombre de casiers à utiliser pour créer
l'histogramme et, si vous le souhaitez, la plage de valeurs de bac. Pour les histogrammes d'images et de tableaux,
la plage doit spécifier les valeurs de données min et max. Pour les histogrammes d'images sur les x et y
axes, la plage doit spécifier les valeurs minimale et maximale de l'image bin. Si cet argument est
omis, le nombre de bacs de sortie pour une table est calculé soit à partir du TLMIN/TLMAX
valeurs des en-têtes (si elles existent dans l'en-tête FITS du tableau pour la colonne spécifiée) ou par
en parcourant les données pour calculer la valeur min et max. Pour une image, le nombre de
bacs de sortie est calculé soit à partir des valeurs d'en-tête DATAMIN/DATAMAX, soit en allant
à travers les données pour calculer la valeur min et max. (Notez que ce dernier calcul peut
échoue si l'image ne tient pas en mémoire.) Si les données sont en virgule flottante (tableau ou
image) et que le nombre de bacs n'est pas spécifié, une valeur par défaut arbitraire de 128 est utilisée.

Pour le traitement des tables binaires, le -w (largeur du bac) peut être utilisé pour spécifier la largeur de
chaque bac plutôt que le nombre de bacs. Ainsi:

funhist test.ev pha 1:100:5

signifie que 5 cases de largeur 20 sont utilisées dans l'histogramme, tandis que :

funhist -w test.ev pha 1:100:5

signifie que 20 cases de largeur 5 sont utilisées dans l'histogramme.

Les données sont divisées en le nombre spécifié de bacs et l'histogramme 1D résultant
(ou projection) est sorti au format de table ASCII. Pour un tableau, la sortie affiche le
valeurs low_edge (inclusive) et hi_edge (exclusive) pour les données. Par exemple, une rangée de 15
table contenant une colonne « pha » dont les valeurs vont de -7.5 à 7.5 peut être traitée ainsi :

[sh] funhist test.ev pha
# fichier de données : /home/eric/data/test.ev
# colonne : pha
# min, max, bacs : -7.5 7.5 15

valeur 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-0.50000000 0.50000000
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 pha 1:6
# fichier de données : /home/eric/data/test.ev
# colonne : pha
# min, max, bacs : 0.5 6.5 6

valeur 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 pha 1:6:3
# fichier de données : /home/eric/data/test.ev
# colonne : pha
# min, max, bacs : 0.5 6.5 3

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

Pour un histogramme de table, le -n(normaliser) peut être utilisé pour normaliser la valeur du bac en
la largeur du bac (c'est-à-dire hi_edge-lo_edge) :

[sh] funhist -n test.ev pha 1:6:3
# fichier de données : test.ev
# colonne : pha
# min, max, bacs : 0.5 6.5 3
# la normalisation de la largeur (val/(hi_edge-lo_edge)) est appliquée

valeur 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

Cela pourrait être utilisé, par exemple, pour produire une courbe de lumière avec des valeurs ayant des unités de
comptes/seconde au lieu de comptes.

Pour un histogramme d'image, la sortie affiche les valeurs d'image basse et haute (toutes deux incluses)
utilisé pour générer l'histogramme. Par exemple, dans l'exemple suivant, 184 pixels avaient un
valeur de 1, 31 avait une valeur de 2, tandis que seulement 2 avait une valeur de 3,4,5,6 ou 7 :

[sh] funhist test.fits
# fichier de données : /home/eric/data/test.fits
# min, max, bacs : 1 7 7

valeur 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

Pour la projection d'axe d'une image, la sortie affiche les bacs d'image bas et haut (les deux
inclus) utilisé pour générer la projection. Par exemple, dans l'exemple suivant, 21
les comptes avaient leur valeur X bin de 2, etc. :

[sh] test funhist.fits x 2:7
# fichier de données : /home/eric/data/test.fits
# colonne : X
# min, max, bacs : 2 7 6

valeur 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] test funhist.fits x 2:7:2
# fichier de données : /home/eric/data/test.fits
# colonne : X
# min, max, bacs : 2 7 2

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

Vous pouvez utiliser gnuplot ou d'autres programmes de traçage pour représenter graphiquement les résultats, en utilisant un script tel
comme:

#!/ 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

Des commandes de tracé similaires sont fournies dans le script funhist.intrigue:

funhist test.ev pha ... ⎪ funhist.plot gnuplot

Utilisez funhist en ligne en utilisant les services onworks.net