funhist - מקוון בענן

תָכְנִית:

שֵׁם

funhist - צור היסטוגרמה 1D של עמודה (מטבלה בינארית FITS או קובץ אירוע גולמי)
או תמונה

תַקצִיר

funhist [-n⎪-w⎪-T] [עמודה] [[lo:hi:]פחים]

אפשרויות

-n # לנרמל את ערך הפח ברוחב של כל פח
-w # ציין את רוחב הפח במקום את מספר הפחים ב-arg3
-T # פלט בפורמט rdb/starbase (מפרידי כרטיסיות)

תיאור

funhist יוצר היסטוגרמה חד מימדית מהעמודות שצוינו של הרחבת FITS
טבלה בינארית של קובץ FITS (או מקובץ אירוע גולמי שאינו FITS), או מתמונת FITS או
מערך, וכותב את ההיסטוגרמה הזו כטבלת ASCII. לחלופין, התוכנית יכולה לבצע
הקרנה 1D של אחד מצירי התמונה.

הארגומנט הראשון לתוכנית נדרש, ומציין את קובץ Funtools: FITS table
או תמונה, קובץ אירוע גולמי או מערך. אם מצוין "stdin", הנתונים נקראים מה-
קלט סטנדרטי. השתמש ב-Funtools Bracket Notation כדי לציין הרחבות FITS ומסננים.

עבור טבלה, נדרש גם הארגומנט השני. זה מציין את העמודה שבה יש להשתמש
יצירת ההיסטוגרמה. אם קובץ הנתונים הוא מסוג image (או מערך), העמודה היא
אופציונלי: אם מצוין "x" (או "X"), "y" (או "Y"), אזי מבוצעת הקרנה מעל
צירי x (dim1) או y (dim2) בהתאמה. (כלומר, השלכה זו תיתן את אותו הדבר
תוצאות כהיסטוגרמה המבוצעת בטבלה המכילה את שורות האירוע המקבילות x,y.) אם
לא צוין שם עמודה או שצוין "xy" (או "XY") עבור התמונה, ולאחר מכן היסטוגרמה
מבוצע על הערכים הכלולים בפיקסלים של התמונה.

הטיעון הבא הוא אופציונלי ומציין את מספר הפחים לשימוש ביצירה
ההיסטוגרמה ואם תרצה, טווח ערכי הפח. להיסטוגרמות של תמונה וטבלה,
הטווח צריך לציין את ערכי הנתונים המינימליים והמקסימליים. עבור היסטוגרמות תמונה על x ו-y
צירים, הטווח צריך לציין את ערכי סל התמונות המינימליים והמקסימליים. אם הטיעון הזה הוא
בהשמטה, מספר פחי הפלט עבור טבלה מחושב מתוך ה-TLMIN/TLMAX
ערכי כותרות (אם אלו קיימים בכותרת FITS בטבלה עבור העמודה שצוינה) או לפי
עוברים על הנתונים כדי לחשב את הערך המינימלי והמקסימלי. עבור תמונה, המספר של
סלי הפלט מחושבים או מתוך ערכי הכותרת DATAMIN/DATAMAX, או על ידי מעבר
דרך הנתונים כדי לחשב ערך מינימלי ומקסימלי. (שים לב שהחישוב האחרון הזה עשוי
נכשל אם התמונה לא יכולה להתאים לזיכרון.) אם הנתונים הם נקודה צפה (טבלה או
image) ומספר הפחים לא צוין, נעשה שימוש ברירת מחדל שרירותית של 128.

עבור עיבוד טבלה בינארית, ה -w ניתן להשתמש במתג (רוחב פח) כדי לציין את הרוחב של
כל סל ולא מספר הפחים. לכן:

funhist test.ev pha 1:100:5

פירוש הדבר ש-5 פחים ברוחב 20 משמשים בהיסטוגרמה, בעוד:

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

פירושו ש-20 פחים ברוחב 5 משמשים בהיסטוגרמה.

הנתונים מחולקים למספר הפחים שצוין ולהיסטוגרמה 1D המתקבלת
(או הקרנה) מופק בפורמט טבלת ASCII. עבור טבלה, הפלט מציג את
ערכי low_edge (כולל) ו-hi_edge (בלעדיים) עבור הנתונים. לדוגמה, 15 שורות
טבלה המכילה עמודת "pha" שהערכים שלה נעים בין -7.5 ל-7.5 ניתנים לעיבוד כך:

[sh] funhist test.ev pha
# קובץ נתונים: /home/eric/data/test.ev
עמודה אחת: pha
# דקות, מקסימום, פחים: -7.5 7.5 15

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

[ש] funhist test.ev פא 1:6
# קובץ נתונים: /home/eric/data/test.ev
עמודה אחת: pha
# דקות, מקסימום, פחים: 0.5 6.5 6

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

[ש] funhist test.ev פא 1:6:3
# קובץ נתונים: /home/eric/data/test.ev
עמודה אחת: pha
# דקות, מקסימום, פחים: 0.5 6.5 3

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

עבור היסטוגרמה טבלה, ה -nניתן להשתמש במתג (נרמול) כדי לנרמל את ערך הפח על ידי
רוחב הפח (כלומר, hi_edge-lo_edge):

[ש] funhist -n test.ev פא 1:6:3
# קובץ נתונים: test.ev
עמודה אחת: pha
# דקות, מקסימום, פחים: 0.5 6.5 3
# נורמליזציה של רוחב (val/(hi_edge-lo_edge)) מוחלת

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

זה יכול לשמש, למשל, כדי לייצר עקומת אור עם ערכים שיש להם יחידות של
ספירות/שנייה במקום ספירות.

עבור היסטוגרמה תמונה, הפלט מציג את ערכי התמונה הנמוכים והגבוהים (כולל שניהם)
משמש ליצירת ההיסטוגרמה. לדוגמה, בדוגמה הבאה, ל-184 פיקסלים היה א
לערך של 1, 31 היה ערך 2, בעוד שרק ל-2 היה ערך של 3,4,5,6 או 7:

[ש] funhist test.fits
# קובץ נתונים: /home/eric/data/test.fits
# דקות, מקסימום, פחים: 1 7 7

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

עבור הקרנת ציר של תמונה, הפלט מציג את סל התמונה הנמוך והגבוה (שניהם
כולל) משמש ליצירת ההקרנה. לדוגמה, בדוגמה הבאה, 21
לספירות היה ערך X bin שלהם 2 וכו':

[sh] funhist test.fits x 2:7
# קובץ נתונים: /home/eric/data/test.fits
עמודה אחת: X
# דקות, מקסימום, פחים: 2 7 6

bin value 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
# קובץ נתונים: /home/eric/data/test.fits
עמודה אחת: X
# דקות, מקסימום, פחים: 2 7 2

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

אתה יכול להשתמש ב-gnuplot או בתוכנות ציור אחרות כדי לתאר את התוצאות, באמצעות סקריפט כזה
כמו:

#!/ Bin / sh
sed -e '1,/---- .*/d
/^$/,$d' ⎪ \
אווק '\
BEGIN{print "set nokey; set title \"funhist\"; set xlabel \"bin\"; set ylabel \"counts\"; plot \"-\" with boxes"} \
{הדפס $3, $2, $4-$3}' ⎪ \
gnuplot -persist - 1>/dev/null 2>&1

פקודות עלילה דומות מסופקות בסקריפט funhist.plot:

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

השתמש ב-funhist באינטרנט באמצעות שירותי onworks.net