funhist - Online en la nube

PROGRAMA:

NOMBRE

funhist: crea un histograma 1D de una columna (a partir de una tabla binaria FITS o un archivo de eventos sin formato)
o una imagen

SINOPSIS

divertido [-n⎪-w⎪-T] [columna] [[lo: hi:] bins]

CAMPUS

-n # normaliza el valor del contenedor por el ancho de cada contenedor
-w # especifica el ancho del contenedor en lugar del número de contenedores en arg3
-T # salida en formato rdb / starbase (separadores de pestañas)

DESCRIPCIÓN

divertido crea un histograma unidimensional a partir de las columnas especificadas de una extensión FITS
tabla binaria de un archivo FITS (o de un archivo de eventos sin procesar que no es FITS), o de una imagen FITS o
matriz y escribe ese histograma como una tabla ASCII. Alternativamente, el programa puede realizar
una proyección 1D de uno de los ejes de la imagen.

El primer argumento del programa es obligatorio y especifica el archivo Funtools: tabla FITS
o imagen, archivo de eventos sin formato o matriz. Si se especifica "stdin", los datos se leen desde el
entrada estándar. Utilice la notación de corchetes de Funtools para especificar extensiones y filtros FITS.

Para una tabla, también se requiere el segundo argumento. Especifica la columna a utilizar en
generando el histograma. Si el archivo de datos es de tipo imagen (o matriz), la columna es
opcional: si se especifica "x" (o "X"), "y" (o "Y"), entonces se realiza una proyección sobre
los ejes x (dim1) o y (dim2), respectivamente. (Es decir, esta proyección dará el mismo
resultados como un histograma realizado en una tabla que contiene las filas de eventos x, y equivalentes).
no se especifica ningún nombre de columna o se especifica "xy" (o "XY") para la imagen, luego un histograma
se realiza sobre los valores contenidos en los píxeles de la imagen.

El argumento que sigue es opcional y especifica el número de bins que se utilizarán para crear
el histograma y, si se desea, el rango de valores de intervalo. Para histogramas de imágenes y tablas,
el rango debe especificar los valores de datos mínimo y máximo. Para histogramas de imágenes en xey
ejes, el rango debe especificar los valores mínimos y máximos del contenedor de imágenes. Si este argumento es
omitido, el número de contenedores de salida para una tabla se calcula a partir de TLMIN / TLMAX
valores de encabezados (si existen en el encabezado de la tabla FITS para la columna especificada) o por
revisando los datos para calcular el valor mínimo y máximo. Para una imagen, el número de
Los contenedores de salida se calculan a partir de los valores del encabezado DATAMIN / DATAMAX, o yendo
a través de los datos para calcular el valor mínimo y máximo. (Tenga en cuenta que este ltimo clculo
fallará si la imagen no cabe en la memoria.) Si los datos son de punto flotante (tabla o
image) y no se especifica el número de bins, se utiliza un valor predeterminado arbitrario de 128.

Para el procesamiento de tablas binarias, el -w (ancho del contenedor) se puede utilizar para especificar el ancho de
cada contenedor en lugar del número de contenedores. Por lo tanto:

funhist test.ev pha 1: 100: 5

significa que se utilizan 5 contenedores de ancho 20 en el histograma, mientras que:

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

significa que se utilizan 20 contenedores de ancho 5 en el histograma.

Los datos se dividen en el número especificado de contenedores y el histograma 1D resultante
(o proyección) se emite en formato de tabla ASCII. Para una tabla, la salida muestra el
valores low_edge (inclusive) y hi_edge (exclusivo) para los datos. Por ejemplo, una fila de 15
tabla que contiene una columna "pha" cuyos valores oscilan entre -7.5 y 7.5 se pueden procesar así:

[sh] prueba funhist.ev pha
# archivo de datos: /home/eric/data/test.ev
# columna: pha
# min, max, bins: -7.5 7.5 15

bin valor 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
# archivo de datos: /home/eric/data/test.ev
# columna: pha
# mín., máx., contenedores: 0.5 6.5 6

bin valor 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
# archivo de datos: /home/eric/data/test.ev
# columna: pha
# mín., máx., contenedores: 0.5 6.5 3

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

Para un histograma de tabla, el -nEl interruptor (normalizar) se puede utilizar para normalizar el valor del intervalo mediante
el ancho del contenedor (es decir, hi_edge-lo_edge):

[sh] funhist -n test.ev pha 1: 6: 3
# archivo de datos: test.ev
# columna: pha
# mín., máx., contenedores: 0.5 6.5 3
# se aplica la normalización de ancho (val / (hi_edge-lo_edge))

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

Esto podría usarse, por ejemplo, para producir una curva de luz con valores que tengan unidades de
recuentos / segundo en lugar de recuentos.

Para un histograma de imagen, la salida muestra los valores de imagen bajo y alto (ambos incluidos)
utilizado para generar el histograma. Por ejemplo, en el siguiente ejemplo, 184 píxeles tenían un
valor de 1, 31 tenía un valor de 2, mientras que solo 2 tenía un valor de 3,4,5,6, 7, XNUMX, XNUMX o XNUMX:

[sh] prueba funhist.ajustes
# archivo de datos: /home/eric/data/test.fits
# mín., máx., contenedores: 1 7 7

valor de 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

Para la proyección del eje de una imagen, la salida muestra las bandejas de imagen baja y alta (ambas
inclusive) utilizado para generar la proyección. Por ejemplo, en el siguiente ejemplo, 21
los recuentos tenían su valor X bin de 2, etc .:

[sh] prueba funhist. Se adapta x 2: 7
# archivo de datos: /home/eric/data/test.fits
# columna: X
# mín., máx., contenedores: 2 7 6

valor de 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] prueba funhist. Se adapta x 2: 7: 2
# archivo de datos: /home/eric/data/test.fits
# columna: X
# mín., máx., contenedores: 2 7 2

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

Puede utilizar gnuplot u otros programas de trazado para graficar los resultados, utilizando un script como
como:

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

Se proporcionan comandos de trazado similares en el script funhist.plot:

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

Utilice funhist en línea utilizando los servicios de onworks.net