funhist - Online in de cloud

PROGRAMMA:

NAAM

funhist - maak een 1D-histogram van een kolom (vanuit een FITS binaire tabel of een onbewerkt gebeurtenisbestand)
of een afbeelding

KORTE INHOUD

grappig [-n⎪-w⎪-T] [kolom] [[lo:hi:]bakken]

OPTIES

-n # normaliseer de bin-waarde met de breedte van elke bin
-w # specificeer de bakbreedte in plaats van het aantal bakken in arg3
-T # uitvoer in rdb/starbase-formaat (tabbladscheidingstekens)

PRODUCTBESCHRIJVING

grappig creëert een eendimensionaal histogram van de opgegeven kolommen van een FITS-extensie
binaire tabel van een FITS-bestand (of van een niet-FITS onbewerkt gebeurtenisbestand), of van een FITS-afbeelding of
array, en schrijft dat histogram als een ASCII-tabel. Als alternatief kan het programma presteren
een 1D-projectie van een van de beeldassen.

Het eerste argument voor het programma is vereist en specificeert het Funtools-bestand: FITS-tabel
of afbeelding, onbewerkt gebeurtenisbestand of array. Als "stdin" is opgegeven, worden gegevens gelezen uit het
standaard invoer. Gebruik Funtools Bracket Notation om FITS-extensies en filters te specificeren.

Voor een tabel is het tweede argument ook vereist. Het specificeert de kolom waarin moet worden gebruikt
het genereren van het histogram. Als het gegevensbestand van het type image (of array) is, is de kolom dat ook
optioneel: als "x" (of "X"), "y" (of "Y") is opgegeven, wordt er een projectie uitgevoerd
respectievelijk de x (dim1) of y (dim2) assen. (Dat wil zeggen, deze projectie zal hetzelfde opleveren
resultaten als een histogram uitgevoerd op een tabel die de equivalente x,y gebeurtenisrijen bevat.) If
er is geen kolomnaam opgegeven of "xy" (of "XY") is opgegeven voor de afbeelding, en vervolgens een histogram
wordt uitgevoerd op de waarden in de beeldpixels.

Het volgende argument is optioneel en specificeert het aantal bakken dat moet worden gebruikt bij het maken
het histogram en, indien gewenst, het bereik van bin-waarden. Voor beeld- en tabelhistogrammen:
het bereik moet de minimale en maximale gegevenswaarden specificeren. Voor beeldhistogrammen op de x en y
assen, moet het bereik de minimale en maximale afbeeldingsbakwaarden specificeren. Als dit argument dat is
weggelaten, wordt het aantal uitvoerbakken voor een tabel berekend op basis van de TLMIN/TLMAX
headers-waarden (als deze bestaan ​​in de tabel FITS-header voor de opgegeven kolom) of by
het doornemen van de gegevens om de min- en max-waarde te berekenen. Voor een afbeelding is het aantal
uitvoerbakken worden berekend op basis van de DATAMIN/DATAMAX-headerwaarden, of door te gaan
door de gegevens om de min- en max-waarde te berekenen. (Merk op dat deze laatste berekening mogelijk
mislukken als de afbeelding niet in het geheugen past.) Als de gegevens drijvende komma zijn (tabel of
afbeelding) en het aantal bakken is niet gespecificeerd, wordt een willekeurige standaard van 128 gebruikt.

Voor binaire tabelverwerking wordt de -w (bakbreedte) schakelaar kan worden gebruikt om de breedte van te specificeren
elke bak in plaats van het aantal bakken. Dus:

funhist test.ev pha 1:100:5

betekent dat er 5 bakken met een breedte van 20 worden gebruikt in het histogram, terwijl:

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

betekent dat er 20 bakken met breedte 5 worden gebruikt in het histogram.

De gegevens worden opgedeeld in het opgegeven aantal bakken en het resulterende 1D-histogram
(of projectie) wordt uitgevoerd in ASCII-tabelformaat. Voor een tabel geeft de uitvoer de
low_edge (inclusief) en hi_edge (exclusief) waarden voor de gegevens. Bijvoorbeeld een 15-rij
tabel met een "pha"-kolom waarvan de waarden variëren van -7.5 tot 7.5 en kan als volgt worden verwerkt:

[sh] funhistische test.ev pha
# gegevensbestand: /home/eric/data/test.ev
# kolom: fa
# min,max,bakken: -7.5 7.5 15

bin-waarde 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
# gegevensbestand: /home/eric/data/test.ev
# kolom: fa
# min,max,bakken: 0.5 6.5 6

bin-waarde 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
# gegevensbestand: /home/eric/data/test.ev
# kolom: fa
# min,max,bakken: 0.5 6.5 3

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

Voor een tabelhistogram wordt de -n(normaliseren) schakelaar kan worden gebruikt om de bin-waarde te normaliseren
de breedte van de bak (dat wil zeggen, hi_edge-lo_edge):

[sh] funhist -n test.ev pha 1:6:3
# gegevensbestand: test.ev
# kolom: fa
# min,max,bakken: 0.5 6.5 3
# breedtenormalisatie (val/(hi_edge-lo_edge)) wordt toegepast

bin-waarde 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

Dit zou bijvoorbeeld kunnen worden gebruikt om een ​​lichtcurve te produceren met waarden met eenheden van
tellingen/seconde in plaats van tellingen.

Voor een beeldhistogram geeft de uitvoer de lage en hoge beeldwaarden weer (beide inclusief)
gebruikt om het histogram te genereren. In het volgende voorbeeld hadden 184 pixels bijvoorbeeld een
waarde van 1, 31 had een waarde van 2, terwijl slechts 2 een waarde van 3,4,5,6 of 7 had:

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

bin-waarde 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

Voor de asprojectie van een afbeelding geeft de uitvoer de lage en hoge afbeeldingsvakken weer (beide
inclusief) gebruikt om de projectie te genereren. In het volgende voorbeeld bijvoorbeeld 21
tellingen hadden een X-bakwaarde van 2, enz.:

[sh] funhist test.past x 2:7
# gegevensbestand: /home/eric/data/test.fits
# kolom: X
# min,max,bakken: 2 7 6

bin-waarde 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
# gegevensbestand: /home/eric/data/test.fits
# kolom: X
# min,max,bakken: 2 7 2

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

U kunt gnuplot of andere plotprogramma's gebruiken om de resultaten in een grafiek weer te geven, met behulp van een script zoals
als:

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

Soortgelijke plotcommando's worden in het script geleverd funhistisch.plot:

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

Gebruik funhist online met behulp van onworks.net-services