funhist – Online in der Cloud

PROGRAMM:

NAME/FUNKTION

funhist – Erstellen Sie ein 1D-Histogramm einer Spalte (aus einer FITS-Binärtabelle oder einer Rohereignisdatei).
oder ein Bild

ZUSAMMENFASSUNG

lustig [-n⎪-w⎪-T] [Spalte] [[lo:hi:]bins]

OPTIONAL

-n # Bin-Wert um die Breite jedes Bins normalisieren
-w # Geben Sie die Bin-Breite anstelle der Anzahl der Bins in arg3 an
-T # Ausgabe im rdb/starbase-Format (Tabulatortrennzeichen)

BESCHREIBUNG

lustig erstellt ein eindimensionales Histogramm aus den angegebenen Spalten einer FITS-Erweiterung
Binärtabelle einer FITS-Datei (oder aus einer Nicht-FITS-Rohereignisdatei) oder aus einem FITS-Bild oder
Array und schreibt dieses Histogramm als ASCII-Tabelle. Alternativ kann das Programm ausgeführt werden
eine 1D-Projektion einer der Bildachsen.

Das erste Argument des Programms ist erforderlich und gibt die Funtools-Datei an: FITS-Tabelle
oder Bild, Rohereignisdatei oder Array. Wenn „stdin“ angegeben ist, werden Daten aus dem gelesen
Standardeingabe. Verwenden Sie die Funtools-Klammernotation, um FITS-Erweiterungen und Filter anzugeben.

Für eine Tabelle ist auch das zweite Argument erforderlich. Es gibt die Spalte an, in der verwendet werden soll
Erzeugen des Histogramms. Wenn die Datendatei vom Typ Bild (oder Array) ist, ist dies die Spalte
optional: wenn „x“ (oder „X“), „y“ (oder „Y“) angegeben ist, dann wird eine Projektion darüber durchgeführt
die x-Achse (dim1) bzw. die y-Achse (dim2). (Das heißt, diese Projektion wird dasselbe ergeben
Ergebnisse als Histogramm, das für eine Tabelle ausgeführt wird, die die entsprechenden x-, y-Ereigniszeilen enthält.) Wenn
Für das Bild wird kein Spaltenname angegeben oder „xy“ (oder „XY“) angegeben, dann ein Histogramm
wird an den in den Bildpixeln enthaltenen Werten durchgeführt.

Das folgende Argument ist optional und gibt die Anzahl der beim Erstellen zu verwendenden Bins an
das Histogramm und, falls gewünscht, den Bereich der Bin-Werte. Für Bild- und Tabellenhistogramme
Der Bereich sollte die minimalen und maximalen Datenwerte angeben. Für Bildhistogramme auf x und y
Achsen sollte der Bereich die minimalen und maximalen Bild-Bin-Werte angeben. Wenn dieses Argument ist
weggelassen, wird die Anzahl der Ausgabefächer für eine Tabelle entweder aus TLMIN/TLMAX berechnet
Header-Werte (sofern diese im Tabellen-FITS-Header für die angegebene Spalte vorhanden sind) oder nach
Gehen Sie die Daten durch, um den Mindest- und Höchstwert zu berechnen. Für ein Bild ist die Anzahl der
Die Ausgabebins werden entweder aus den DATAMIN/DATAMAX-Headerwerten oder durch Gehen berechnet
Durchsuchen Sie die Daten, um den Mindest- und Höchstwert zu berechnen. (Beachten Sie, dass diese letztere Berechnung möglicherweise
schlägt fehl, wenn das Bild nicht in den Speicher passt.) Wenn es sich bei den Daten um Gleitkommadaten handelt (Tabelle oder
Bild) und die Anzahl der Bins nicht angegeben ist, wird ein willkürlicher Standardwert von 128 verwendet.

Für die Binärtabellenverarbeitung ist die -w Mit dem Schalter (Bin-Breite) kann die Breite angegeben werden
jedes Fach und nicht die Anzahl der Fächer. Daher:

funhist test.ev pha 1:100:5

bedeutet, dass im Histogramm 5 Bins der Breite 20 verwendet werden, während:

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

bedeutet, dass im Histogramm 20 Bins der Breite 5 verwendet werden.

Die Daten werden in die angegebene Anzahl von Bins unterteilt und das resultierende 1D-Histogramm erstellt
(oder Projektion) wird im ASCII-Tabellenformat ausgegeben. Bei einer Tabelle zeigt die Ausgabe die an
low_edge-Werte (inklusive) und hi_edge-Werte (exklusiv) für die Daten. Zum Beispiel eine 15-Reihe
Eine Tabelle mit einer „pha“-Spalte, deren Werte zwischen -7.5 und 7.5 liegen, kann folgendermaßen verarbeitet werden:

[sh] funhist test.ev pha
# Datendatei: /home/eric/data/test.ev
# Spalte: pha
# min, max, Bins: -7.5 7.5 15

Bin-Wert 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 pha 1:6
# Datendatei: /home/eric/data/test.ev
# Spalte: pha
# min, max, Bins: 0.5 6.5 6

Bin-Wert 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
# Datendatei: /home/eric/data/test.ev
# Spalte: pha
# min, max, Bins: 0.5 6.5 3

Bin-Wert lo_edge hi_edge
------ --------- ------ -------------- -------
1 33 0.50000000 2.50000000
2 37 2.50000000 4.50000000
3 41 4.50000000 6.50000000

Für ein Tabellenhistogramm ist das -nDer Schalter (Normalisieren) kann verwendet werden, um den Bin-Wert zu normalisieren
die Breite des Behälters (dh hi_edge-lo_edge):

[sh] funhist -n test.ev pha 1:6:3
# Datendatei: test.ev
# Spalte: pha
# min, max, Bins: 0.5 6.5 3
# Breitennormalisierung (val/(hi_edge-lo_edge)) wird angewendet

Bin-Wert 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

Dies könnte beispielsweise verwendet werden, um eine Lichtkurve mit Werten in der Einheit zu erzeugen
counts/second statt counts.

Bei einem Bildhistogramm zeigt die Ausgabe die niedrigen und hohen Bildwerte (beide inklusive) an.
wird zur Erstellung des Histogramms verwendet. Im folgenden Beispiel hatten beispielsweise 184 Pixel eine
Wert von 1, 31 hatte einen Wert von 2, während nur 2 einen Wert von 3,4,5,6 oder 7 hatte:

[sh] funhist test.fits
# Datendatei: /home/eric/data/test.fits
# min, max, Bins: 1 7 7

Bin-Wert 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

Für die Achsenprojektion eines Bildes zeigt die Ausgabe die unteren und oberen Bildbereiche (beide) an
einschließlich) zur Generierung der Projektion verwendet. Im folgenden Beispiel beispielsweise 21
Zählungen hatten ihren X-Bin-Wert von 2 usw.:

[sh] funhist test.fits x 2:7
# Datendatei: /home/eric/data/test.fits
# Spalte: X
# min, max, Bins: 2 7 6

Bin-Wert 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
# Datendatei: /home/eric/data/test.fits
# Spalte: X
# min, max, Bins: 2 7 2

Bin-Wert lo_bin hi_bin
------ --------------------- --------- -------------------
1 60.00000000 2.00000000 4.00000000
2 51.00000000 5.00000000 7.00000000

Sie können Gnuplot oder andere Plotprogramme verwenden, um die Ergebnisse grafisch darzustellen, indem Sie ein Skript verwenden, z
als:

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

Ähnliche Plotbefehle werden im Skript bereitgestellt funhist.plot:

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

Nutzen Sie Funhist online über die Dienste von onworks.net