Pare knyfów z IRAFa

Jak używać do kroskorelacji template-ów syntetycznych?

W sieci dostępne są widma syntetyczne w różnej postaci np. fits, txt i inne.

Często jednak fitsy nie są czytane przez IRAFa.

Przykład:

Na stronie http://gaia.esa.int/spectralib mamy możliwosć wygenerować sobie

widmo dla zadanych przez nas parametrach. Otrzymane widmo FITS jest niestety

nieczytelne dla IRAFa. Jednak można sobie z tym poradzić naprzykład wczytując

widmo do GAIA(Starlink) i eksportując do pliku tekstowego (lambda i natężenie).

IRAF może wczytać taki plik komendą rspectext(onedspec).

Przykładowe parametry komendy rspectext:

input = b.dat Input list of text spectra

output = wynik Output list of image spectra

(title = ) Spectrum title

(flux = yes) Flux calibrated?

(dtype = linear) Dispersion type

(crval1 = 7653.0815) Coordinate of first pixel

(cdelt1 = 0.3751) Coordinate interval per pixel

(fd1 = )

(fd2 = )

(mode = ql)

Uzyskujemy plik .imh +.pix. Jednak heder jest pusty. Musimy dodać:

DATE-OBS= '1993-03-17T04:56:38.0'

RA = '12:00:00'

DEC = '12:00:00'

EPOCH = 1993.0

OBSERVAT= 'KPNO'

VHELIO = 0.0

EXPTIME= 0.00000001

Dzięki temu FXCOR przyjmie ten template. Te parametry dają nam heliocentryczną poprawkę

0 z dokładnością 5 znaków po przecinku co oznacza ze nasze Vhelio pozostanie 0.0. Ta barbarzyńska metoda jest opisana w helpie FXCOR.

Filtracja Furierowska dla FXCOR

Filtracja pozwala czasami rozdzielic 2 piki w funkcji CCF jeżeli zlewanie jest spowodowane szumem.

1. Włączamy opcje “filter” (none | both | object | templeit) w “epar fxcor”.

przykładowe ustawienia:

PACKAGE = rv

TASK = filtpars

(f_type = square) Filter window type

(cuton = 1) Cuton component for filter

(cutoff = 200) Cutoff component for filter

(fullon = 0) Component at which (ramp) filter reaches full value

(fulloff= 0) Component at which (ramp) filter reaches zero

(mode = ql)

CDN