Astrometria na danych z Brwinowa

Sprawdzenie astrometrii na danych z Brwinowa :

1. Dane z nocy 20060226 :

Zrobilem astrometrie na danych z nocy 20060226, klatka fitc , zgodnie z informacjami od GW kamera byla wycelowana w rownik dokladnie na
poludnie co oznacza dla BRW ~= WAW (az,h)=(180,38) , co na 2006-02-26 19:48 UT dalo pozycje (RA,DEC)~=(7.7 h, 0 deg )
Zakladajac PIXSCALE = 15um / 85 mm =~ 36.4 , oraz korzystajac z informacji GW ze gwiazdy wedruja po chipie gora-dol zrobilem najpierw
vertical flip :

calcfits canon_best.fit FV canon_best_fv.fit

a nastepnie wykonalem skrypt do photo/astro. Obrot do astrometri wynosil -90 ( -fi -90 ), co oczywiscie wyniklo obortu kamery na
montazu stalym.
Po astrometrii wyszlo (RA,DEC)=(7.75 h , 1.86 deg), POSANGLE=-90.6 , oraz nastepujace parametery kamery :
PIXSCALE=36.00

2/ Zasieg

Skatalogowalem nastepnie ta pojedyncza klatke ast ( jest w nim 17536 gwizd) , do bazy 20060226@heplx43, mozna z niej zobaczyc od razu nastepujace rzeczy :

star_num | max_mag | min_mag |   fov_ra          |   fov_dec
- ---------+---------+---------+-------------+-------------
17164 | 16.5396    | 4.44122   | 20.06450325 | 19.96280479

czyli skatalogowalo sie troche mniej gwiazd niz jest w ast ( bladne pomiary sa odrzucane ) , max_mag=16.53, co jest chyba dosc podejrzane
dlatego sprawdzilem gwiazdy ktore maja odpowidnik w katalogu i to dalo :

star_num | max_mag | min_mag | fov_ra   |   fov_dec
----------+---------+---------+-----------+-------------
    10430 | 14.883 | 4.84314 | 19.948689 | 19.95724773
(1 row)

Co wyglada ze jest bardziej prawdopodobna wartoscia, zrobilem tez rozklady magnitudo , dla wszystkich gwiazd w bazie , oraz tylko
dla gwiazd ktore maja odpowiednik w katalogu, wiecej plotow i eps-y mozna znalezc tutaj
Widac ze zasieg w Brwinowie przy chyba dosc kiepskich warunkach wyniosl na pojedynczej klatce ~11 mag.
( nie bylo flatowania , trzeba zebrac klatki przy stojacej kamerze do zrobienia flata ! )
Dla porownania te same ploty na danych z LCO f=50mm, dla klatki k2a_050510_01000.fitc , wszystkie gwiazdy, tylko znane w katalog cat/act (ASAS)
Baza 20050510 na www ( sorry for old version of interface )
Swoje badania robilem na heplx43:/opt/pi/dev/pisys/daq/ndir/data/20060226/ oraz w podkatalogu f50mm.

FWHM :
FWHM = 3.97811341

3. Dane z nocy 20060307 :

Na poczatku uzylem skryptu fix_radec! zeby dla danych w ktorych nie bylo zapisywanych wspolrzednych do headera ( z test2K2K )
wpisac odpowiednie dane do headera niezbedne dla astrometrii, dodatkowo wylicza on RA,DEC z AZIM,ALT na dana chwile czasu ( czyta
go z headera ).
Chyba cos zostalo przekrecone bo teraz potrzebny byl flip horizontal , fi=-90.6 , zeby poszlo i zrobilem szybka fotometrie
( robi sie w heplx43:/disk01/results/20060307/FH/FH) , natomiast proby robilem w heplx43:/disk01/results/20060307/audela/fits/FH

Dane po odjeciu darka i obrocie znajduja sie tutaj : heplx43:/disk01/results/20060307/FH/FH
dodaje do nich pare brakujacych kluczy w headerze a nastepnie zapuszcze na nich szybka fotometrie i astrometrie , bedzie mozna
je skatalogowac, ale to co ja chce zrobic to raczej wyznaczyc sigma_g/g dla roznych laplacjanow w funkcji jasnosci gwiazd
i wybrac mozliwie najlepszy laplacjan do nowych obiektywow

Poza tym przekrecone klatki beda juz mialy zapisana w fits headerze astrometrie i bedzie mozna ich uzyc do analizy algorytmow z
potwierdzeniem na nastepnej , wzor config jest w /disk01/results/20060307/FH/FH/CCDSINGLE
to jest config ktory czyta astrometrie z fits headera

FWHM :
FWHM = 3.15979362

ZASIEG - DOBRY SPOSOB
Zrobilem katalogowanie tej nocy : heplx49:/opt/pi/dev/pisys/daq/ndir/results/tmp do bazy brw20060307@heplx49 , nastepnie zrobilem
selekcje gwiazd o no_measurements>=300 i sprawdzilem rozklad ich jasnosci , sposob z mag_cat>0 nie jest dobry bo w
katalgou ktorego uzywam do normalizacji nie ma gwiazd o malych jasnosciach !!!

ZASIEG WEDLUG TYCHO,GCVS :
Tycho nie jest dobry do badania zasiegu z powodu malego zasiegu, w momencie gdy bedziemy mieli zasiec > 13 mag to nie bedzie on wtedy dobry
rozklad jasnosci w tycho jest nastepujacy ( log ).
Gcvs jest lepszy ale ma gwiazdy zmienne, wiec niestabilne, ale za to ma szeroki zakres magnitudo , rozklad w gcvs jest taki ( log ).
Wiecej plotow , testy wykonalem w heplx43:/disk01/msok/analysis/zasieg/
Chyba najlepsza metoda pozostaje robienie rozkladu jasnosci gwiazd o duzej liczbie pomiarow - wtedy to na pewno nie sa fluktuacje
i mozna to tak zrobic ( ewentualnie zobaczyc ich rozklad na chipie - czy gdzies sie na brzegach nie grupuja )
Dla porownania mozna zmatchowac z gcvs i zobaczyc jak sie ma rozklad jasnosci tych gwiazd ktore maja odpowiednik w GCVS
Ale tycho nie jest lepszy do tego celu niz cat/act ( asas - astrometry catalog ) bo to chyba jest to samo, a przynajmniej zasieg ma podobny.
To co mozna jeszcze sprobowac zrobic, to sciagnac z SIMBAD-a pelny katalog ( tylko gwiazd ) i na tej podstawie to juz byloby idealnie chyba

ALGORITHM PARAMETERS :
Analysis of algorithms on this data was done by Lech Piotrowski, results are available here .

4. 20060312 - dwie klatki

Na tych dwoch klatkach wykonalem astrometrie :
    fix_radec! - trzeba pamietac zeby miec odpowiedni ccd.cfg z obs_site_brwinow.cfg zamiast LCO ( sprawdzic RA,DEC ~ equator powinien byc )
    piphoto fits_list -list=fits_list -dir=../mag1/ -dark=dark1.fit -tresh=3.5 -flip -force
    ./do_astrometry! mag_list ../ast1/ -90
- czyli fi=-90 , i poszlo , powstaly asty , <FWHM> :
    k2a_060312_00004.ast FWHM = 3.15979362
    k2a_060312_00010.ast FWHM = 3.15979362
- robilem to na swoim komputerze : /opt/pi/dev/pisys/daq/ndir/data/20060312/ast1
- nastepnie katalogowanie skrypt : cat_single_night! :
   nohup ./cat_single_night! "-no_update -skip_field_change -max_ast_err=1. -min_alt=0" > out 2>&1 &
- selekcja gwiazd znanaych w katalogu :
    psql 20060312 -U pidb_user -c "select magnitude from stars where mag_cat>0" > mag_cat.txt
    .x histofile.C("mag_cat.txt",0,0,4,13,100,0,"Magnitude","Star#     ")
otrzymany rozklad jasnosci. Widac ze tym razem zasieg wynosil okolo 10.5 mag. Czyli o okolo 0.5 mag gorzej niz 20060226 ?.

5. Analiza danych w poszukiwaniu blyskow

Zostala wykonana na danych z 20060307 na heplx43:/disk01/results/20060307/FH/FH/CCDSINGLE , wyniki sa dostepne tutaj ( new ) oraz
z wlaczonym odrzucaniem duzych gwiazd ( <4 mag ) blizej niz 600 arcsec.
NOTE : in order to view fits parts with events! script flip must be set to 0 ( CCD_ASAS_REVERSE_FOR_TRANSFORM=0 ) , otherwise
(RA,DEC) coordinates shown in the program will not be correct

6. Analiza klatek

Wzrokowo wygladaja dobrze, ale na wycinkach widac ze duze gwiazdy wygladaja dosc dziwnie i maja dziwne promieniste obwodki ( przyklad lub tutaj )
Inny przyklad dziwnych obwodek .
Jeszcze sie im poprzygladam

7. Dane 20060323 - k2b

przygotowanie darka :
cd /disk01/results/20060323/DARK/
prep_dark!
poprawienie (Lat,Long) w fits header i policzenie (RA,DEC) z (AZ,H,TIME) :
fix_radec! k2b
ustawienie w ccd_pipeline0.cfg :
     CCD_ASAS_REVERSE_FOR_TRANSFORM=2
     CCD_ASAS_ASTROMETRY_FI=90
     CCD_USE_ASAS_ASTROMETRY_FROM_FITS=0

NOTE : when testing astrometry on this data I decided to change code in DeviceBase.cpp :
    AddHeader( ROTATE , 1 );
to :
AddHeader( ROTATE , 0 );
it must be verified that this is good change ... but as comment in code says astrometry for ROTATE=1 adds 180 to fi and
for fi=-90 converts fi =90 -> fi=-90 ... this is quite strange and problematic so I changed it ( on 20060327 )
zapuszczenie puszkiwania blyskow z potwierdzeniem na nastepnej klatce w : heplx43:/opt/pi/dev/pisys/daq/ndir/results/20060323 , wyniki
analizy blyskow z potwierdzeniem na nastepnej klatce dla nocy 20060323

8. Dla porownania wykonalem katalogowanie starych danych 20040425 z BRW

Zrobilem ich katalogowanie ( heplx43:/disk01/msok/tmp/BRW/fits ) i sprawdzilem zasieg ktory wygladal nastepujaco.
Stad wniosek ze zasieg w Brwinowie na starych obiektywach wynosil ~10.4 mag
Ponizsza tabelka podumowuje zasiegi w BRW i LCO na pojedynczych klatkach :

OBSERWATORIUM	OBIEKTYW	ZASIEG	DANE
LCO	carl zeiss 50mm	11.5	fast@heplx48 noce : 20050506 , 20050509, 20050510, 20050512, 20050513
BRW	carl zeiss 50mm	10.5	brw2004@heplx43 - 20040425
LCO	cannon 85mm	?	?
BRW	cannon 85mm	11.0	20060226@heplx43 , brw20060307@heplx49

Ponizsza tabelka podumowuje zasiegi w BRW i LCO na klatkach posumowanych po 20 :

OBSERWATORIUM	OBIEKTYW	ZASIEG	DANE
LCO	carl zeiss 50mm	12.5	2004_2005@heplx49 - asas-pipeline data 2004-2005
BRW	carl zeiss 50mm	?	?
LCO	cannon 85mm	?	?
BRW	cannon 85mm	?	?