นี่คือคำสั่ง funcnts ที่สามารถเรียกใช้ในผู้ให้บริการโฮสต์ฟรีของ OnWorks โดยใช้เวิร์กสเตชันออนไลน์ฟรีของเรา เช่น Ubuntu Online, Fedora Online, โปรแกรมจำลองออนไลน์ของ Windows หรือโปรแกรมจำลองออนไลน์ของ MAC OS
โครงการ:
ชื่อ
funcnts - นับโฟตอนในพื้นที่ที่กำหนดพร้อมการลบ bkgd
เรื่องย่อ
ฟังก์ชั่น [สวิตช์] [source_region] [bkgd_file] [bkgd_region⎪bkgd_value]
OPTIONS
-e "source_exposure[;bkgd_exposure]"
# แหล่งที่มา (bkgd) พอดีกับภาพที่เปิดรับโดยใช้ไฟล์ที่ตรงกัน
-w "source_exposure[;bkgd_exposure]"
# แหล่งที่มา (bkgd) พอดีกับภาพการเปิดรับแสงโดยใช้ WCS transform
-t "source_timecorr[;bkgd_timecorr]"
# ค่าการแก้ไขเวลาต้นทาง (bkgd) หรือชื่อพารามิเตอร์ส่วนหัว
-g # เอาต์พุตโดยใช้รูปแบบ g ที่ดี
-G # เอาต์พุตโดยใช้รูปแบบ %.14g (ความแม่นยำสูงสุด)
-i "[คอลัมน์;]int1;int2..." # ช่วงเวลาตามคอลัมน์
-m # จับคู่แหล่งที่มาแต่ละรายการและภูมิภาค bkgd
-p # เอาต์พุตเป็นพิกเซล แม้ว่าจะมี wcs อยู่
-r # เอาต์พุตรัศมีด้านใน/ด้านนอก (และมุม) สำหรับวงแหวน (และแพนด้า)
-s # ค่าผลลัพธ์รวม
-v "scol[;bcol]" # คอลัมน์ค่า src และ bkgd สำหรับตาราง
-T # เอาต์พุตในรูปแบบ starbase/rdb
-z # พื้นที่เอาต์พุตที่มีพื้นที่เป็นศูนย์
DESCRIPTION
ฟังก์ชั่น นับโฟตอนในพื้นที่ต้นทางที่ระบุและรายงานผลลัพธ์สำหรับแต่ละ
ภาค. ขอบเขตถูกกำหนดโดยใช้กลไกการกรองพื้นที่เชิงพื้นที่ โฟตอนคือ
ยังนับในพื้นที่ bkgd ที่ระบุซึ่งใช้กับไฟล์ข้อมูลเดียวกันหรือต่างกัน
แฟ้มข้อมูล. (หรือ ค่าพื้นหลังคงที่ในจำนวน/พิกเซล**2 สามารถเป็น
ระบุ) ภูมิภาค bkgd ถูกจับคู่แบบหนึ่งต่อหนึ่งกับภูมิภาคต้นทางหรือรวมกลุ่ม
และทำให้เป็นมาตรฐานตามพื้นที่ แล้วลบออกจากจำนวนต้นทางในแต่ละภูมิภาค
ผลลัพธ์ที่แสดงรวมจำนวนลบ bkgd ในแต่ละภูมิภาค เช่นเดียวกับข้อผิดพลาด
ในการนับ พื้นที่ในแต่ละภูมิภาค และความสว่างของพื้นผิว (cnts/area**2)
คำนวณสำหรับแต่ละภูมิภาค
อาร์กิวเมนต์แรกของโปรแกรมระบุ FITS อิมเมจอินพุต อาร์เรย์ หรือไฟล์เหตุการณ์ดิบ
ในการประมวลผล. หากระบุ "stdin" ข้อมูลจะถูกอ่านจากอินพุตมาตรฐาน ใช้ Funtools
เครื่องหมายวงเล็บเพื่อระบุส่วนขยาย ส่วนของรูปภาพ และฟิลเตอร์ของ FITS
อาร์กิวเมนต์ที่สองที่เป็นทางเลือกคือตัวบอกขอบเขตแหล่งที่มา หากไม่มีการระบุภูมิภาค
ใช้ทั้งฟิลด์
อาร์กิวเมนต์พื้นหลังสามารถมีรูปแบบใดรูปแบบหนึ่งจากสองรูปแบบ ขึ้นอยู่กับว่าแยกจากกันหรือไม่
ไฟล์พื้นหลังถูกระบุ หากไฟล์ต้นฉบับจะใช้เป็นพื้นหลังด้วย คำสั่ง
อาร์กิวเมนต์ที่สามสามารถเป็นพื้นที่พื้นหลังหรือค่าคงที่แสดงถึง
พื้นหลัง cnts/พิกเซล อีกทางหนึ่ง อาร์กิวเมนต์ที่สามอาจเป็นไฟล์ข้อมูลแบ็กกราวด์
ในกรณีนี้อาร์กิวเมนต์ที่สี่คือพื้นที่พื้นหลัง ถ้าไม่มีอาร์กิวเมนต์ที่สามคือ
ระบุ ใช้ค่าคงที่ 0 (เช่น ไม่มีพื้นหลัง)
โดยสรุป อาร์กิวเมนต์คำสั่งต่อไปนี้ถูกต้อง:
[sh] funcnts sfile # นับในไฟล์ต้นฉบับ
[sh] funcnts sfile sregion # นับในพื้นที่ต้นทาง
[sh] funcnts sfile sregion bregion # bkgd reg. มาจากไฟล์ต้นฉบับ
[sh] funcnts sfile sregion bvalue # bkgd reg. เป็นค่าคงที่
[sh] funcnts sfile sregion bfile bregion # bkgd reg. มาจากไฟล์แยกต่างหาก
หมายเหตุ: ไม่เหมือนกับโปรแกรม Funtools อื่น ๆ ภูมิภาคต้นทางและพื้นหลังถูกระบุเป็น
แยกอาร์กิวเมนต์บนบรรทัดคำสั่ง แทนที่จะวางไว้ในวงเล็บเป็น part
ของชื่อไฟล์ต้นทางและพื้นหลัง นี่เป็นเพราะว่าภูมิภาคใน funcnts ไม่ได้เรียบง่าย
ใช้เป็นตัวกรองข้อมูล แต่ยังใช้ในการคำนวณพื้นที่ การรับแสง ฯลฯ หากคุณใส่
ภูมิภาคต้นทางภายในวงเล็บ (เช่น ใช้เป็นตัวกรองอย่างง่าย) แทนที่จะระบุ
เป็นอาร์กิวเมนต์ที่สอง โปรแกรมจะนับเฉพาะโฟตอนที่ผ่านตัวกรองภูมิภาคเท่านั้น
อย่างไรก็ตาม การคำนวณพื้นที่จะดำเนินการทั้งสนามตั้งแต่ สนาม() คือ
ภูมิภาคต้นทางเริ่มต้น นี่เป็นพฤติกรรมที่ต้องการน้อยมาก ในทางกลับกันกับ FITS
ตารางไบนารีมักจะเป็นประโยชน์ในการใส่ตัวกรองคอลัมน์ในวงเล็บชื่อไฟล์เพื่อที่
เฉพาะเหตุการณ์ที่ตรงกับตัวกรองคอลัมน์เท่านั้นที่จะถูกนับภายในภูมิภาค
ตัวอย่างเช่น หากต้องการแยกจำนวนภายในรัศมี 22 พิกเซลจากจุดศูนย์กลางของ
เหมาะกับตารางไบนารี snr.ev และลบพื้นหลังที่กำหนดจากภาพเดียวกันภายใน
วงแหวนรัศมี 50-100 พิกเซล:
[sh] funcnts snr.ev "วงกลม(502,512,22)" "วงแหวน(502,512,50,100)"
# แหล่งที่มา
# ไฟล์ข้อมูล: snr.ev
# องศา/พิกเซล: 0.00222222
# พื้นหลัง
# ไฟล์ข้อมูล: snr.ev
# หน่วยคอลัมน์
# พื้นที่: arcsec**2
# surf_bri: cts/arcsec**2
# surf_err: cts/arcsec**2
# ผลลัพธ์การลบพื้นหลัง
reg net_counts พื้นหลังข้อผิดพลาด berror พื้นที่ surf_bri surf_err
---- ------------ --------- ------------ --------- ---- ----- --------- ---------
1 3826.403 66.465 555.597 5.972 96831.98 0.040 0.001
# มีการใช้ส่วนประกอบต้นทางและพื้นหลังต่อไปนี้:
ภูมิภาคต้นทาง
----------------
วงกลม(502,512,22)
reg นับพิกเซล
-
1 4382.000 1513
ภูมิภาคพื้นหลัง
--------------------
วงแหวน(502,512,50,100)
reg นับพิกเซล
-
ทั้งหมด 8656.000 23572
หน่วยพื้นที่สำหรับคอลัมน์ผลลัพธ์ที่ระบุว่า "พื้นที่", "surf_bri" (ความสว่างของพื้นผิว) และ
"surf_err" จะได้รับในอาร์ควินาที (หากข้อมูล WCS ที่เหมาะสมอยู่ใน
ส่วนหัวของไฟล์ข้อมูลหรือหน่วยเป็นพิกเซล หากไฟล์ข้อมูลมีข้อมูล WCS แต่คุณไม่ต้องการ arc-
หน่วยที่สอง ใช้ -p สลับเพื่อบังคับเอาต์พุตเป็นพิกเซล นอกจากนี้ ภูมิภาคที่มีพื้นที่ศูนย์
ปกติจะไม่รวมอยู่ในตารางหลัก (ลบพื้นหลัง) แต่รวมอยู่ด้วย
ในซอร์สรองและตาราง bkgd หากคุณต้องการให้ภูมิภาคเหล่านี้รวมอยู่ใน
ตารางหลัก ใช้ the -z สวิตซ์.
โปรดทราบว่าคำสั่ง sed อย่างง่ายจะแยกผลลัพธ์ที่ลบพื้นหลังออกเพื่อเพิ่มเติม
การวิเคราะห์:
[sh] แมว funcnts.sed
1,/---- .*/ด
/^$/,$d
[sh] sed -f funcnts.sed funcnts.out
1 3826.403 66.465 555.597 5.972 96831.98 0.040 0.001
หากระบุไฟล์ต้นฉบับและไฟล์พื้นหลังแยกกัน ฟังก์ชั่น จะพยายามทำให้เป็นปกติ
พื้นที่พื้นหลังเพื่อให้ขนาดพิกเซลพื้นหลังเท่ากับพิกเซลต้นทาง
ขนาด. การทำให้เป็นมาตรฐานนี้สามารถเกิดขึ้นได้ก็ต่อเมื่อข้อมูล WCS ที่เหมาะสมคือ
มีอยู่ในทั้งสองไฟล์ (เช่น ค่าองศา/พิกเซลใน CDELT) หากไฟล์ใดไฟล์หนึ่งไม่
มีข้อมูลขนาดที่จำเป็น ไม่มีการทำให้เป็นมาตรฐาน ในกรณีนี้,
เป็นความรับผิดชอบของผู้ใช้ที่จะต้องแน่ใจว่าขนาดพิกเซลเท่ากันสำหรับทั้งสอง
ไฟล์
โดยปกติ หากระบุพื้นที่พื้นหลังมากกว่าหนึ่งพื้นที่ ฟังก์ชั่น จะรวมพวกมันทั้งหมด
ให้เป็นพื้นที่เดียวและใช้พื้นที่พื้นหลังนี้เพื่อสร้างพื้นหลังที่หักออก
ผลลัพธ์สำหรับแต่ละภูมิภาคต้นทาง NS -m (จับคู่หลายพื้นหลัง) สวิตช์บอก ฟังก์ชั่น
เพื่อสร้างการติดต่อแบบหนึ่งต่อหนึ่งระหว่างพื้นหลังและภูมิภาคต้นทางแทน
โดยใช้ขอบเขตพื้นหลังเดียวที่รวมกัน ตัวอย่างเช่น กรณีเริ่มต้นคือการรวม2
ขอบเขตพื้นหลังเป็นภูมิภาคเดียว แล้วใช้ภูมิภาคนั้นกับแต่ละแหล่งที่มา
ภูมิภาค:
[sh] funcnts snr.ev "วงแหวน(502,512,0,22,n=2)" "วงแหวน(502,512,50,100,n=2)"
# แหล่งที่มา
# ไฟล์ข้อมูล: snr.ev
# องศา/พิกเซล: 0.00222222
# พื้นหลัง
# ไฟล์ข้อมูล: snr.ev
# หน่วยคอลัมน์
# พื้นที่: arcsec**2
# surf_bri: cts/arcsec**2
# surf_err: cts/arcsec**2
# ผลลัพธ์การลบพื้นหลัง
reg net_counts พื้นหลังข้อผิดพลาด berror พื้นที่ surf_bri surf_err
---- ------------ --------- ------------ --------- ---- ----- --------- ---------
1 3101.029 56.922 136.971 1.472 23872.00 0.130 0.002
2 725.375 34.121 418.625 4.500 72959.99 0.010 0.000
# มีการใช้ส่วนประกอบต้นทางและพื้นหลังต่อไปนี้:
ภูมิภาคต้นทาง
----------------
วงแหวน(502,512,0,22,n=2)
reg นับพิกเซล
-
1 3238.000 373
2 1144.000 1140
ภูมิภาคพื้นหลัง
--------------------
วงแหวน(502,512,50,100,n=2)
reg นับพิกเซล
-
ทั้งหมด 8656.000 23572
โปรดทราบว่ากฎพื้นฐานของตัวกรองภูมิภาค "แต่ละโฟตอนจะถูกนับเพียงครั้งเดียวและไม่มีโฟตอนคือ
นับมากกว่าหนึ่งครั้ง" ยังคงมีผลเมื่อใช้ The -m เพื่อให้ตรงกับพื้นที่พื้นหลัง ที่
คือ หากพื้นที่พื้นหลังสองส่วนทับซ้อนกัน พิกเซลที่ทับซ้อนกันจะถูกนับในหนึ่งเดียว
ของพวกเขา. ในกรณีที่เลวร้ายที่สุด หากพื้นที่พื้นหลังสองส่วนเป็นภูมิภาคเดียวกัน
อันดับแรกจะนับจำนวนและพื้นที่ทั้งหมด และครั้งที่สองจะไม่ได้อะไรเลย
การใช้ -m เปลี่ยนสาเหตุ ฟังก์ชั่น เพื่อใช้แต่ละพื้นที่พื้นหลังทั้งสองอย่างแยกจากกัน
กับภูมิภาคต้นทางสองแห่ง:
[sh] funcnts -m snr.ev "วงแหวน(502,512,0,22,n=2)" "แอน(502,512,50,100,n=2)"
# แหล่งที่มา
# ไฟล์ข้อมูล: snr.ev
# องศา/พิกเซล: 0.00222222
# พื้นหลัง
# ไฟล์ข้อมูล: snr.ev
# หน่วยคอลัมน์
# พื้นที่: arcsec**2
# surf_bri: cts/arcsec**2
# surf_err: cts/arcsec**2
# ผลลัพธ์การลบพื้นหลัง
reg net_counts พื้นหลังข้อผิดพลาด berror พื้นที่ surf_bri surf_err
---- ------------ --------- ------------ --------- ---- ----- --------- ---------
1 3087.015 56.954 150.985 2.395 23872.00 0.129 0.002
2 755.959 34.295 388.041 5.672 72959.99 0.010 0.000
# มีการใช้ส่วนประกอบต้นทางและพื้นหลังต่อไปนี้:
ภูมิภาคต้นทาง
----------------
วงแหวน(502,512,0,22,n=2)
reg นับพิกเซล
-
1 3238.000 373
2 1144.000 1140
ภูมิภาคพื้นหลัง
--------------------
แอน(502,512,50,100,n=2)
reg นับพิกเซล
-
1 3975.000 9820
2 4681.000 13752
โปรดทราบว่าปริมาณทศนิยมส่วนใหญ่จะแสดงโดยใช้รูปแบบ "f" เปลี่ยนได้
นี่เป็นรูปแบบ "g" โดยใช้ -g สวิตซ์. ซึ่งจะมีประโยชน์เมื่อนับในแต่ละพิกเซล
มีขนาดเล็กมากหรือใหญ่มาก หากคุณต้องการความแม่นยำสูงสุดและไม่สนใจเกี่ยวกับ
เรียงแถวอย่างสวยงาม ใช้ -Gซึ่งส่งออกค่าลอยตัวทั้งหมดเป็น %.14g
เมื่อนับโฟตอนโดยใช้วงแหวนและรูปร่างแพนด้า (พายและวงแหวน) มักจะเป็น
มีประโยชน์ในการเข้าถึงรัศมี (และมุมแพนด้า) สำหรับแต่ละภูมิภาคที่แยกจากกัน NS -r
สวิตช์จะเพิ่มคอลัมน์รัศมีและมุมลงในตารางผลลัพธ์:
[sh] funcnts -r snr.ev "วงแหวน(502,512,0,22,n=2)" "แอน(502,512,50,100,n=2)"
# แหล่งที่มา
# ไฟล์ข้อมูล: snr.ev
# องศา/พิกเซล: 0.00222222
# พื้นหลัง
# ไฟล์ข้อมูล: snr.ev
# หน่วยคอลัมน์
# พื้นที่: arcsec**2
# surf_bri: cts/arcsec**2
# surf_err: cts/arcsec**2
#รัศมี: arcsecs
# มุม: องศา
# ผลลัพธ์การลบพื้นหลัง
reg net_counts พื้นหลังข้อผิดพลาด berror พื้นที่ surf_bri surf_err รัศมี1 รัศมี2 มุม1 มุม2
---- ------------ --------- ------------ --------- ---- ----- --------- --------- --------- --------- --------- ---------
1 3101.029 56.922 136.971 1.472 23872.00 0.130 0.002 0.00 88.00 ไม่มี
2 725.375 34.121 418.625 4.500 72959.99 0.010 0.000 88.00 176.00 ไม่มี
# มีการใช้ส่วนประกอบต้นทางและพื้นหลังต่อไปนี้:
ภูมิภาคต้นทาง
----------------
วงแหวน(502,512,0,22,n=2)
reg นับพิกเซล
-
1 3238.000 373
2 1144.000 1140
ภูมิภาคพื้นหลัง
--------------------
แอน(502,512,50,100,n=2)
reg นับพิกเซล
-
ทั้งหมด 8656.000 23572
รัศมีกำหนดเป็นหน่วยพิกเซลหรืออาร์ควินาที (ขึ้นอยู่กับการมีอยู่ของข้อมูล WCS)
ในขณะที่ค่ามุม (เมื่อมี) อยู่ในหน่วยองศา คอลัมน์เหล่านี้สามารถใช้เพื่อพล็อต
โปรไฟล์เรเดียล ตัวอย่างเช่น สคริปต์ functs.พล็อต ในการกระจาย funtools) will
พล็อตโปรไฟล์เรเดียลโดยใช้ gnuplot (เวอร์ชัน 3.7 ขึ้นไป) เวอร์ชั่นย่อของสิ่งนี้
สคริปต์แสดงอยู่ด้านล่าง:
#!bin / sh /
ถ้า [ x"$1" = xgnuplot ]; แล้ว
ถ้า [ x`อันไหน gnuplot 2>/dev/null` = x ]; แล้ว
echo "ข้อผิดพลาด: gnuplot ไม่พร้อมใช้งาน"
ทางออก 1
fi
อ๊าก '
เริ่ม{HEADER=1; ข้อมูล=0; ไฟล์=""; XLABEL="ไม่ทราบ"; YLABEL="ไม่ทราบ"}
ส่วนหัว==1{
ถ้า (1 == "#" && $2 == "ข้อมูล" && $3 == "ไฟล์:" ){
ถ้า ( FILES != "" ) FILES = ไฟล์ ","
ไฟล์ = ไฟล์ $4
}
อื่น if( $1 == "#" && $2 == "รัศมี:" ){
XLABEL = $3
}
อื่น if( $1 == "#" && $2 == "surf_bri:" ){
YLABEL = $3
}
อย่างอื่น if($1 == "----"){
printf "set nokey; ตั้งชื่อ \"funcnts(%s)\"\n", FILES
printf "set xlabel \" radius(%s)\"\n", XLABEL
printf "set ylabel \"surf_bri(%s)\"\n", YLABEL
พิมพ์ "พล็อต \"-\" โดยใช้ 3:4:6:7:8 พร้อมกล่องข้อผิดพลาด"
ส่วนหัว = 0
ข้อมูล = 1
ถัดไป
}
}
ข้อมูล==1{
ถ้า (NF == 12 ){
พิมพ์ $9, $10, ($9+$10)/2, $7, $8, $7-$8, $7+$8, $10-$9
}
อื่น{
ทางออก
}
}
' ⎪ gnuplot -persist - 1>/dev/null 2>&1
เอลฟ์ [ x"$1" = xds9 ]; แล้ว
อ๊าก '
เริ่ม{HEADER=1; ข้อมูล=0; XLABEL="ไม่ทราบ"; YLABEL="ไม่ทราบ"}
ส่วนหัว==1{
ถ้า (1 == "#" && $2 == "ข้อมูล" && $3 == "ไฟล์:" ){
ถ้า ( FILES != "" ) FILES = ไฟล์ ","
ไฟล์ = ไฟล์ $4
}
อื่น if( $1 == "#" && $2 == "รัศมี:" ){
XLABEL = $3
}
อื่น if( $1 == "#" && $2 == "surf_bri:" ){
YLABEL = $3
}
อย่างอื่น if($1 == "----"){
printf "funcnts(%s) รัศมี(%s) surf_bri(%s) 3\n", FILES, XLABEL, YLABEL
ส่วนหัว = 0
ข้อมูล = 1
ถัดไป
}
}
ข้อมูล==1{
ถ้า (NF == 12 ){
พิมพ์ $9, $7, $8
}
อื่น{
ทางออก
}
}
'
อื่น
echo "funcnts -r ... ⎪ funcnts.plot [ds9⎪gnuplot]"
ทางออก 1
fi
ดังนั้นการวิ่ง ฟังก์ชั่น และพล็อตผลลัพธ์โดยใช้ gnuplot (เวอร์ชัน 3.7 ขึ้นไป) ใช้:
funcnts -r snr.ev "วงแหวน(502,512,0,50,n=5)" ... ⎪ funcnts.plot gnuplot
พื้นที่ -s (ผลรวม) สวิตช์สาเหตุ ฟังก์ชั่น เพื่อสร้างตารางเพิ่มเติมของผลรวม (แบบบูรณาการ)
ค่าที่ลบพื้นหลังพร้อมกับตารางเริ่มต้นของแต่ละค่า:
[sh] funcnts -s snr.ev "วงแหวน(502,512,0,50,n=5)" "วงแหวน(502,512,50,100)"
# แหล่งที่มา
# ไฟล์ข้อมูล: snr.ev
# องศา/พิกเซล: 0.00222222
# พื้นหลัง
# ไฟล์ข้อมูล: snr.ev
# หน่วยคอลัมน์
# พื้นที่: arcsec**2
# surf_bri: cts/arcsec**2
# surf_err: cts/arcsec**2
# สรุปผลการลบพื้นหลัง
มากถึง net_counts พื้นหลังข้อผิดพลาด berror พื้นที่ surf_bri surf_err
---- ------------ --------- ------------ --------- ---- ----- --------- ---------
1 2880.999 54.722 112.001 1.204 19520.00 0.148 0.003
2 3776.817 65.254 457.183 4.914 79679.98 0.047 0.001
3 4025.492 71.972 1031.508 11.087 179775.96 0.022 0.000
4 4185.149 80.109 1840.851 19.786 320831.94 0.013 0.000
5 4415.540 90.790 2873.460 30.885 500799.90 0.009 0.000
# ผลลัพธ์การลบพื้นหลัง
reg นับข้อผิดพลาดพื้นหลัง berror พื้นที่ surf_bri surf_err
---- ------------ --------- ------------ --------- ---- ----- --------- ---------
1 2880.999 54.722 112.001 1.204 19520.00 0.148 0.003
2 895.818 35.423 345.182 3.710 60159.99 0.015 0.001
3 248.675 29.345 574.325 6.173 100095.98 0.002 0.000
4 159.657 32.321 809.343 8.699 141055.97 0.001 0.000
5 230.390 37.231 1032.610 11.099 179967.96 0.001 0.000
# มีการใช้ส่วนประกอบต้นทางและพื้นหลังต่อไปนี้:
ภูมิภาคต้นทาง
----------------
วงแหวน(502,512,0,50,n=5)
reg นับพิกเซล sumcnts sumix
-
1 2993.000 305 2993.000 305
2 1241.000 940 4234.000 1245
3 823.000 1564 5057.000 2809
4 969.000 2204 6026.000 5013
5 1263.000 2812 7289.000 7825
ภูมิภาคพื้นหลัง
--------------------
วงแหวน(502,512,50,100)
reg นับพิกเซล
-
ทั้งหมด 8656.000 23572
พื้นที่ -t และ -e สามารถใช้สวิตช์เพื่อปรับจังหวะเวลาและการรับแสงตามลำดับ
ข้อมูล โปรดทราบว่าการแก้ไขเหล่านี้มีขึ้นเพื่อใช้ในเชิงคุณภาพตั้งแต่
การใช้ปัจจัยการแก้ไขที่แม่นยำยิ่งขึ้นนั้นเป็นความพยายามที่ซับซ้อนและขึ้นอยู่กับภารกิจ
อัลกอริทึมสำหรับการใช้การแก้ไขอย่างง่ายเหล่านี้มีดังนี้:
C = จำนวนดิบในภูมิภาคต้นทาง
Ac= พื้นที่ของภูมิภาคต้นทาง
Tc= เวลาเปิดรับแสงสำหรับแหล่งข้อมูล
Ec= การเปิดรับแสงเฉลี่ยในพื้นที่ต้นทาง จากแผนที่การเปิดรับแสง
B= จำนวนดิบในพื้นที่พื้นหลัง
Ab= พื้นที่ของภูมิภาคพื้นหลัง
Tb= เวลา (การรับแสง) สำหรับข้อมูลพื้นหลัง
Eb= การเปิดรับแสงเฉลี่ยในพื้นที่พื้นหลัง จากแผนที่การเปิดรับแสง
จากนั้น จำนวนสุทธิในภูมิภาคต้นทางคือ
สุทธิ= C - B * (Ac*Tc*Ec)/(Ab*Tb*Eb)
ด้วยการเผยแพร่ข้อผิดพลาดมาตรฐานสำหรับ Error on Net อัตราสุทธิจะเป็น
อัตราสุทธิ = สุทธิ/(Ac*Tc*Ec)
ค่าแสงเฉลี่ยในแต่ละภูมิภาคคำนวณโดยการรวมค่าพิกเซลใน
แผนที่แสงสำหรับภูมิภาคที่กำหนดแล้วหารด้วยจำนวนพิกเซลในนั้น
ภาค. แผนที่แสงมักจะถูกสร้างขึ้นที่ปัจจัยการบล็อก > 1 (เช่น บล็อก 4 หมายความว่า
พิกเซลการรับแสงแต่ละพิกเซลมีพิกเซล 4x4 ที่ความละเอียดเต็ม) และ ฟังก์ชั่น จะจัดการกับ
การปิดกั้นโดยอัตโนมัติ ใช้ -e เปลี่ยน คุณสามารถจัดหาทั้งแหล่งที่มาและพื้นหลัง
ไฟล์การเปิดเผย (คั่นด้วย ";") หากคุณมีไฟล์ข้อมูลต้นทางและพื้นหลังแยกจากกัน
หากคุณไม่ได้ให้ไฟล์เปิดเผยพื้นหลังเพื่อไปกับข้อมูลพื้นหลังแยกต่างหาก
ไฟล์, ฟังก์ชั่น ถือว่ามีการใช้การเปิดเผยกับไฟล์ข้อมูลพื้นหลังแล้ว
นอกจากนี้ จะถือว่าข้อผิดพลาดบนพิกเซลในไฟล์ข้อมูลพื้นหลังเป็นศูนย์
NB: -e สวิตช์ถือว่าแผนที่แสงซ้อนทับไฟล์ภาพ เผงยกเว้น
สำหรับปัจจัยการบล็อก แต่ละพิกเซลในภาพจะถูกปรับขนาดตามปัจจัยการบล็อกเพื่อเข้าถึง
พิกเซลที่เกี่ยวข้องในแผนที่การเปิดรับแสง หากแผนที่เปิดรับแสงของคุณไม่ตรงกัน
ด้วยภาพ, do ไม่ ใช้ -e การแก้ไขการรับแสง ในกรณีนี้ก็ยังได้
เพื่อทำการแก้ไขค่าแสง if ทั้งภาพและแผนที่เปิดรับแสงมี WCS . ที่ถูกต้อง
ข้อมูล: ใช้ -w สลับเพื่อให้การแปลงจากพิกเซลภาพเป็นแสง
pixel ใช้ข้อมูล WCS นั่นคือแต่ละพิกเซลในพื้นที่ภาพจะเป็น
เปลี่ยนจากพิกัดภาพเป็นพิกัดท้องฟ้าก่อน จากนั้นเปลี่ยนจากพิกัดท้องฟ้าเป็น
พิกัดการเปิดรับแสง โปรดทราบว่าการใช้ -w สามารถเพิ่มเวลาที่ใช้ในการดำเนินการได้
การแก้ไขการรับแสงอย่างมาก
การแก้ไขเวลาสามารถใช้ได้กับทั้งแหล่งข้อมูลและข้อมูลพื้นหลังโดยใช้ปุ่ม -t สวิตซ์.
ค่าสำหรับการแก้ไขอาจเป็นค่าคงที่ตัวเลขหรือชื่อของส่วนหัวก็ได้
พารามิเตอร์ในไฟล์ต้นทาง (หรือพื้นหลัง):
[sh] funcnts -t 23.4 ... # หมายเลขสำหรับแหล่งที่มา
[sh] funcnts -t "LIVETIME;23.4" ... # param สำหรับแหล่งที่มา ตัวเลขสำหรับ bkgd
เมื่อมีการระบุการแก้ไขเวลาจะนำไปใช้กับการนับสุทธิด้วย (ดู
อัลกอริธึมด้านบน) เพื่อให้หน่วยของความสว่างของพื้นผิวกลายเป็น cnts/area**2/sec
พื้นที่ -i (ช่วง) ใช้สวิตช์เพื่อเรียกใช้ ฟังก์ชั่น ในช่วงเวลาตามคอลัมน์หลาย ๆ ตัวด้วย
ผ่านข้อมูลเพียงครั้งเดียว เทียบเท่าการวิ่ง ฟังก์ชั่น หลายครั้ง
ด้วยตัวกรองคอลัมน์อื่นที่เพิ่มไปยังแหล่งข้อมูลและข้อมูลพื้นหลังในแต่ละครั้ง แต่ละ
ช่วงเต็ม ฟังก์ชั่น สร้างเอาต์พุตโดยใส่อักขระป้อนบรรทัด (^L)
ระหว่างการวิ่งแต่ละครั้ง นอกจากนี้ ผลลัพธ์สำหรับแต่ละช่วงเวลาจะมีช่วงเวลา
ข้อกำหนดในส่วนหัว ช่วงเวลามีประโยชน์มากสำหรับการสร้างความแข็งของเอ็กซ์เรย์
อัตราส่วนได้อย่างมีประสิทธิภาพ แน่นอนว่ารองรับก็ต่อเมื่อมีข้อมูลอินพุตอยู่
ในตาราง
รองรับสองรูปแบบสำหรับข้อกำหนดช่วงเวลา รูปแบบทั่วไปส่วนใหญ่เป็นแบบกึ่ง
รายการตัวกรองที่คั่นด้วยโคลอนเพื่อใช้เป็นช่วงเวลา:
funcnts -i "pha=1:5;pha=6:10;pha=11:15" snr.ev "circle(502,512,22)" ...
ตามแนวคิดแล้วสิ่งนี้จะเทียบเท่ากับการวิ่ง ฟังก์ชั่น สามครั้ง:
funcnts snr.ev'[pha=1:5]' "วงกลม(502,512,22)"
funcnts snr.ev'[pha=6:10]' "วงกลม(502,512,22)"
funcnts snr.ev'[pha=11:15]' "วงกลม(502,512,22)"
อย่างไรก็ตามการใช้ไฟล์ -i สวิตช์จะต้องผ่านข้อมูลเพียงครั้งเดียว
โปรดทราบว่าสามารถใช้ตัวกรองที่ซับซ้อนเพื่อระบุช่วงเวลาได้:
funcnts -i "pha=1:5&&pi=4;pha=6:10&&pi=5;pha=11:15&&pi=6" snr.ev ...
โปรแกรมเพียงแค่รันข้อมูลผ่านตัวกรองแต่ละตัว และสร้างสาม ฟังก์ชั่น
เอาต์พุต คั่นด้วยอักขระป้อนบรรทัด
ในความเป็นจริง ถึงแม้ว่าเจตนาจะสนับสนุนช่วงเวลาสำหรับอัตราส่วนความแข็ง ค่าที่ระบุ
ตัวกรองไม่จำเป็นต้องเป็นช่วงๆ เลย ไม่จำเป็นต้องมีตัวกรอง "ช่วงเวลา" ตัวเดียว
เกี่ยวข้องกับผู้อื่น ตัวอย่างเช่น:
funcnts -i "pha=1:5;pi=6:10;energy=11:15" snr.ev "circle(502,512,22)" ...
เทียบเท่ากับการวิ่ง ฟังก์ชั่น สามครั้งด้วยข้อกำหนดตัวกรองที่ไม่เกี่ยวข้อง
รูปแบบช่วงเวลาที่สองได้รับการสนับสนุนสำหรับกรณีง่าย ๆ ที่ใช้คอลัมน์เดียว
เพื่อระบุช่วงที่เป็นเนื้อเดียวกันหลายช่วงสำหรับคอลัมน์นั้น ในรูปแบบนี้ ชื่อคอลัมน์
ถูกระบุก่อน ตามด้วยช่วงเวลา:
funcnts -i "pha;1:5;6:10;11:15" snr.ev "circle(502,512,22)" ...
นี่เทียบเท่ากับตัวอย่างแรก แต่ต้องการการพิมพ์น้อยกว่า NS ฟังก์ชั่น โครงการ
จะเติม "pha=" นำหน้าแต่ละช่วงที่ระบุ (โปรดทราบว่ารูปแบบนี้
ไม่มีอักขระ "=" ในอาร์กิวเมนต์ของคอลัมน์)
โดยปกติเมื่อ ฟังก์ชั่น ทำงานบนตารางไบนารี FITS (หรือตารางเหตุการณ์ดิบ) one
จำนวนนับรวมจะถูกสะสมสำหรับแต่ละแถว (เหตุการณ์) ที่มีอยู่ภายในภูมิภาคที่กำหนด NS -v
"scol[;bcol]" (คอลัมน์ค่า) สวิตช์จะสะสมการนับโดยใช้ค่าจาก
คอลัมน์ที่ระบุสำหรับเหตุการณ์ที่กำหนด หากระบุคอลัมน์เดียว จะใช้สำหรับ
ทั้งภูมิภาคต้นทางและพื้นหลัง สองคอลัมน์แยกกัน คั่นด้วยเซมิโคลอน
สามารถระบุแหล่งที่มาและพื้นหลังได้ โทเค็นพิเศษ '$none' สามารถใช้เพื่อ
ระบุว่าจะใช้คอลัมน์ค่าสำหรับคอลัมน์หนึ่งแต่ไม่ใช่คอลัมน์อื่น ตัวอย่างเช่น,
'pha;$none' จะใช้คอลัมน์ pha สำหรับแหล่งที่มา แต่ใช้จำนวนอินทิกรัลสำหรับ
พื้นหลัง ในขณะที่ '$none;pha' จะทำการสนทนา หากคอลัมน์ค่าเป็นประเภท
ตรรกะ ดังนั้นค่าที่ใช้จะเป็น 1 สำหรับ T และ 0 สำหรับ F สำหรับคอลัมน์ค่าที่ใช้สำหรับ
ตัวอย่าง เพื่อรวมความน่าจะเป็นแทนการนับจำนวนเต็ม
ถ้า -T ใช้สวิตช์ (ตาราง rdb) เอาต์พุตจะสอดคล้องกับฐานข้อมูล starbase/rdb
รูปแบบ: แท็บจะถูกแทรกระหว่างคอลัมน์แทนที่จะเป็นช่องว่างและการป้อนบรรทัดจะเป็น
แทรกระหว่างตาราง
สุดท้ายโปรดทราบว่า ฟังก์ชั่น เป็นโปรแกรมสร้างภาพ แม้จะรันได้โดยตรงบน FITS
ตารางไบนารี ซึ่งหมายความว่าใช้การกรองภาพกับแถวเพื่อให้แน่ใจว่า
เพื่อให้ได้ผลลัพธ์ที่เหมือนกันโดยไม่คำนึงถึงว่าตารางหรือเทียบเท่า binned
ภาพถูกใช้ อย่างไรก็ตาม ด้วยเหตุนี้ จำนวนการนับที่พบโดยใช้ ฟังก์ชั่น สามารถ
แตกต่างจากจำนวนเหตุการณ์ที่พบโดยใช้โปรแกรมกรองแถว เช่น กองทุน or
สนุกได้ สำหรับข้อมูลเพิ่มเติมเกี่ยวกับความแตกต่างเหล่านี้ ดูการสนทนาของภูมิภาค
ขอบเขต
ใช้ funcnts ออนไลน์โดยใช้บริการ onworks.net