4.1. Xác định cấp sao nhìn thấy của các thiên thé bang phần mềm IRAF
4.1.1. Xác định cấp sao của một ngôi sao
Gói phot của IRAF cho phép ta do cấp sao của các ngôi sao dya trên phương pháp quang trắc Aperture. Đặc biệt, gói này rất thuận tiện cho việc đo cấp sao của chỉ một ngôi sao, rất đễ đàng và tiện lợi.
Bước l: Ta có thé truy cập gói này bằng cách gõ như sau:
digiphot <enter>
apphot <enter>
Bước 2: Trước khi thực hiện đo cap sao, ta can chỉnh sửa một vài thông số của gói này
(Hình 4.1):
epar phot <enter>
PACKAGE = apphot TASK = phot
image
skufile
(coords (output = (plotfilz (datapar=
(centerp=
(fitskup=
(photpar=
(interac=
(radplot=
(icoenan=
(gcomman=
(ucsin (ucsout (cache
mm.
Inage Reduction and Analysis Facility
ẹ W0_hoanchinh. fits The input image(s)
The input sky File(s)
) The input coordinate files(s) (default: image.co default) The output photometry File(s) (default: inase.a
) The output plots xetacode File ) Data dependent parameters
) Centering parameters ) Sky Fitting parameters ) Photometry parameters yes) Interactive mode ?
yes) Plot the radial profiles in interactive node 2 ) Image cursor: [x y wes] key [ond]
) Graphics cursor: lễ y wos) key [ond]
)_.wesin) The input coordinate system (logical, tv.physical )_.wesout) The output coordinate system (logical, tv,physica
)..cache) Cache the input image pixels in maxory ?
).verify) Verify critical parameters in non-interactive mo HERA fon HELP
Hình 4.1: Các thông số của gói phot
37
Ta chi quan tâm đến các thông số: datapars, centerpars, fitskypars và photpars. Gõ
“e" ở mỗi thông số dé chinh sửa các thông số ân bên trong.
Chỉnh sửa các thông số của đatapars (Hình 4.2), centerpars (Hình 4.3), fitskypars
(Hình 4.4), và photpars (Hình 4.5) như trong hình:
mm.
Inage Reduction and Analysis Facility PACKAGE = apphot
TASK = datapars
(scale = § 1.) Inage scale in units per pixel (fuhses£= 4,) FWHM of the PSF in scale units
(emissio= yes) Features are positive ?
(signa = INDEF) Standard deviation of background in counts
(datamin= INDEF) Minimum good data value
(datamax= INDEF) Haxiwun sood data value
(noise = poisson) Noise model
(cedread= ) CCD readout noise image header keyword
(gan = ) CCD gain image header keyword
(readnoi= 6.5) CCD readout noise in electrons
(epadu = 2.23399339522165) Gain in electrons per count
(exposur= ) Exposure time image header keyword
(airnassz ) Rirmass image header Keyword
(filter = ) Filter image header keyword
(obst ime= ) Time of observation image header keyword
(itine = 1.) Exposure tine
(xairnas= INDEF) Airmass
Hore]
SEN for HELP
Hình 4.2: Các thông sé của datapars
Inage Reduction and Analysis Facility PACKAGE = apphot
TREK = centerpas
(calgori= centroid) Centering algorithm
(ebox = §f 11,2) Centering box width in scale units
(cthresh= 0.) Centering threshold in sigma above background
(Rinsrra= 1,) Hinisu signal-to-noise ratio for centering algo (cmaxite= 10) Maximum number of iterations for centering algor
(naxshif= 10.) Maximum center shift in scale units
(clean = no) Supmetry clean before centering 2 (relean = 1,) Cleaning radius in scale units
(relip = 2.) Clipping radius in scale units (kelean = 3.) Rejection Limit in sigma
eee no Mark the computed center on display 2
s g
BERD for HELP
Hình 4.3: Các thông số của centerpars
38
Đi với các thông số cbox, annulus, dannulus, apertures ta xác định như sau:
chox = 2*FWHM
annulus = 4*FWHM
dannulus = 3,5*FWHM apertures = 3*FWHM [3]
Image Reduction and Analysis Facility PACKACE = apphot
TASK = fitskypars
(salgori= median) Sky fitting algorithn
(arewlus= §j 22.4) Inner radius of sky arewlus in scale units (dwnuluz 19.6) Uidth of sky annulus in scale units
(skwualus 0,) User sku value
(snaxite= 10) Haxinus runber of sky Fitting iterations Csloclip= 0.) Lower clipping Factor in percent
(shiclip= 0,) Upper clipping factor in percent
(srreJec= 54) Maximum number of sky fitting rejection iteratio (sloreje= 3.) Lower K-signa rejection Linit in sky signa
(shireje= 3.) Upper K~siqằa rejection limit in sky signa
(khist = 3.) Half uidth of histogram in sky signa
(binsize= 0.10000000149012) Binsize of histogram in sky signa
(smooth = no) Boxcar’ smooth the histogran
(rgow = 0.) Region growing radius in scale units
(Msku = no) Mark sky annuli on the display
(mode = ql)
; ERG for HELP Hình 4.4: Các thông số của fitskypars
[ TRAF |
Image Reduction and Analysis Facility PACKAGE = apphot
TASK = photpars
(weight i= constant) Photometric veighting schewe for uphot (apertur= § 16.8) List of aperture radii in scale units
(zmag z 25.) Zero point of magnitude scale
(nkapert= yes) Draw apertures on the display
(mode = ql)
BERR for HELP
Hình 4.5: Các thông số của photpars 39
Một thông so mà ta can quan tâm nữa là zmag. zmag là điềm gốc tinh cấp sao. Gia trị zmag lả đặc trưng cho từng đêm, nó phụ thuộc vảo điều kiện thời tiết, môi trường ở đêm
đỏ. Vì vậy, mỗi khi đo cấp sao của ảnh chụp ở một đêm bat kỳ, ta phải tính toán lại giá trị
zmag đề điều chỉnh lại trong IRAF.
De tính gia trị zmag, ban dau ta lay một ngôi sao làm chuan. Ta sẽ xác định cap sao chính xác trên lý thuyét mir của sao nảy băng phan mềm Stellarium (Hình 4.6). Sau đó, ta sẽ đo cap sao trên thực nghiệm mrx băng phan mem IRAF khi zmag = 25 (gia tri zmag mặc định của IRAF là 25):
Bước 3: Dé do cấp sao mry, ta truy cập vào gói phot:
phot <enter>
Nhập tên của hình ảnh chứa ngôi sao can do (vi dụ: sa0116863_hoanchinh.fits) rôi nhắn Enter. Lúc nay, con trỏ chuột sẽ chuyển sang nhấp nháy. Đưa chuột qua cửa sô DS9
tới vi trí ngôi sao can đo, nhân phím F. Trên cửa so IRAF sẽ xuat hiện các thông tin của
ngôi sao, trong đó có cấp sao (Hình 4.7).
= w50_hoarchinh.fits The input image(s)
skyfile = The input sky file(s)
(coords } The input coordinate files(s) (default: image.co
(output default) The cutput photometry file(s) (default: image.na
(plotfil= } The cutput plots netacode file
(datapar= } Data dependent pa-axeters (centerp= } Centering parameters
(Fitskyp= ) Sky Fitting parancters
(photpx= } Photometry paraneters
Cinterac= yes) Interactive mode ?
(radplot= yes) Plot the radial profiles in interactive mode ? (icommarr= ) Image cursor: (x y wes) key Cond)
(geonmarr= ) Graphics cursor: tx y ves) key [ond]
(wesin = )_.wesin) The input coordinate system (logical .tv,physical
(wesout = )_.wesout) The output coordinate syste (logical. tv.physica (cache = )_.cache) Cache the input image pixels in xenory ?
me = )_.verify) Yerify critical parameters in non-interactive mo
apphot> phot
The input image(s) (m50_hoanchinh, fits): sa0116863_hoanchinh, fits
Warning: Graphics overlay not available for display dev
ssol16865_hosnchinh. Fite 378.76 243,57 3.437027 ok
Hình 4.7: Cap sao của SAO116863 đo được bằng phan mém IRAF
Bước 4: Sau khi đã có hai giá trị mm va mị+, ta sẽ xác định zmag bằng phương trình sau:
mrw — Mur = 25 — zmag [4]
Bước 5: Sau khi xác định xong giá trị zmag của đêm đó, ta sẽ điều chỉnh lại giá trị này
trong thông số photpars. Sau đó, ta sẽ giừ nguyên giá trị này dé tính toán cấp sao cho
những ngôi sao khác.
Bước 6: Thực hiện lại bước 3 dé do cap sao của những ngôi sao khác.