ASTROMETRY
1. Nieliniowosc :
Klatka k2a_061123_001.txt
Wzor inicjalizcja :
param->px[0] =
param->py[0] = 0.;
param->px[1] =
pixscale * cos(rd);
param->px[2] =
-pixscale * sin(rd);
param->py[1] =
pixscale * sin(rd);
param->py[2] =
pixscale * cos(rd);
for(i=3; i<21; i++)
param->px[i]=param->py[i]=0.;
I checked what is the linear and nonlinear part of the correction. It
was calculated as :
linx.fit = ABS(linX) / (
ABS(linX) + ABS(nonlinX) )
nonlinx.fit = ABS(nonlinX) / (
ABS(linX) + ABS(nonlinX) )
liny.fit = ABS(linY) / (
ABS(linY) + ABS(nonlinY) )
nonliny.fit = ABS(nonlinY) / (
ABS(linY) + ABS(nonlinY) )
such plots can be see
here
( use Audela or other FITS soft to open them and multiply by constant
value 100 ) . It can be seen that
nonlinear part is of order of 1-2% in the corners and
practicaly 0 in the center.
2. Skakanie wycinkow
- skakanie klatek przy wycinaniu - czemu i jak to sie dzieje ???
/pi20/msok/data/20060521
k2a_060521_00600.fitc RA = 19.34132532 DEC = 9.90468906
find_frames_list! 19.34132532 9.90468906 20060522_000000 5
- - /pi20/msok/data/ > out 2>&1 &
mam przyklad :
file:///pi20/msok/data/20060521/Event_19.34132532_9.90468906_20060522_000000/normal_aver.html
skrypt doit! - robi astro na wszystkich klatkach i
zobaczymy czy
przestanie skakac , jak nie to przyczyna jest inna !
mimo astrometrii na kazdej klatce nadal troche skacze ...
ale chyba mniej skok dopiero na Frame 12, 2006-05-22,
09:06:16.500 / = k2a_060521_00612.fitc
a wczesniej juz na Frame 10, 2006-05-22, 09:05:48.550 / ,
piad2xy k2a_060521_00612.fitc
19.34132532 9.90468906
wyglada na to ze to jest kwestia dokladnosci astrometrii
,tzn trzeba
sprobowac poptrzyc na konkretna gwiazde - jak dobrze
mierzymy jej pozycje
mozna to zrobic na 20061127fast@heplx49 - tutaj to
jest nawet 114 arcsec !!!
a na takiej nocnej astro tak samo - czyli moze trzeba
lepszego bledu astro wymagac ?
Wziasc jakas gwiazde na k2a_060521_00612.fitc sprawdzic jej
pozycje w
SIMBAD i zobaczyc jak sie wylicza jej pozycja i to jest ta
kwestia pewnie
ze taka jest dokladnosc astro i pewnie to trzeba
ewentualnie poprawic ...
gwiazda : 19.676278,10.47075122
sprawdzic czy z nowa astro skacze - jak tak sprawdzic ad2xy!
skryptem czy
to widac w x,y wyliczonym z ra,dec ???
Wyglada na to ze to jest sprawa zaokraglen + dryf, wylicza sie
pozycja
gwiazdy i jest ciag :
bash-2.05b$ ad2xy k2a_060521_00602.ast
19.676278 10.47075122
537.166797 1083.463282
bash-2.05b$ ad2xy k2a_060521_00603.ast
19.676278 10.47075122
537.143834 1083.489035
bash-2.05b$ ad2xy k2a_060521_00604.ast
19.676278 10.47075122
537.145634 1083.525949
bash-2.05b$ ad2xy k2a_060521_00605.ast
19.676278 10.47075122
537.086322 1083.608147
bash-2.05b$ ls
W pewnym
momencie pozycja przeskakuje o pixel i wycinamy juz w innej pozycji
oczywsicie tak naprawde nie mamy dokladnosci astrometrii z dokaldnoscia
do 0.1
pixela wiec czasem przeskok nie jest.
OSTATECZNIE :
1/ astrometria normalnie nie jest
robiona na kazdej klatce to jest pewnie glowne zrodlo przeskokow
2/ nawet jak jest na kazdej klatce to
przeskoki sa spowodowane tym ze montaz ma dryf, a astrometria nie ma
takiej dokladnosci zeby trafic w moment kiedy i
wlasciwie nie jest
mozliwe zeby to wyeliminowac
3. Liczenie AstroFailed ( pole m_nFailedAstrometryCount )
To pole bylo do tej pory flaga w CCDPipeline, natomiast oczywiscie w
sytuacji gdy :
- astro zbieglo flaga m_pAsasTransform->m_bTransformOK=1
- astro nie zbieglo , flaga w strukturze
m_pAsasTransform->m_bTransformOK=1, ale
m_pAsasTransformAsynchro->m_bTransformOK=0 , wtedy jeszcze
to 0 nie bylo kopiowane do m_pAsasTransform bo nie ustawiala sie flaga
m_bAstrometryIsReady=TRUE
- stad counter bledow m_nFailedAstrometryCount sie zwiekszal,
natomiast flaga byla ciagle stara !!!!
To co zrobilem to po kazdej astrometri ustawia sie flaga
m_bAstrometryIsReady=TRUE ( takze gdy astro failed i trzeba je
powtarzac w petli do gCCDParams.m_nASASAstrometryReTry razy
Zatem teraz juz flaga m_pAsasTransform->m_bTransformOK powinna
zostac zupdatowana, dodatkowo counter umiescilem w tej strukturze
rozniez czyli beda one
synchronizowane razem , a wiec konsystentnie.
Zapiski z badan i zmian :
blad - czemu astro=1 a jest ze
ostatnio 3 razy bylo failed ????
ASTROMETRY = 1
FAILED_ASTRO= 3
/var/www/html/pi0/status/arch/20070102/daq.status
dodalem nowe tracy ...
czemu jest ASTRO-OK -
kiedy ret=0 ???
/pi20/msok/data/20070102/run_daq_20070102.out
WYSZEDL PRZY OKAZJI
POWAZNIEJSZY BLAD !
ZE JAK ASTRO RAZ ZBIEGA A
POTEM PRZESTAJE TO W daq.status juz tego nie
widac - a to dlatego ze
update struktury w glownym watku nie jest robiony !
robia taka poprawke ze po
retry ustawi sie to pole ze
mozna zupdatowac strukture
- i przynajmniej bedzie to robione !
oraz pzeniose countery do
struktry zeby one byly razem z flaga AstroOK
synchronizowane !
TEST : puscilem dobre pole
na symulatorze , a potem ustawilem w gdb :
set
m_PipelineCfg.m_fASASAstrometryFi=180
i astro przestala zbiegac, a co
za tym idzie bylo tak jak trzeba
najpierw zbiegala , a
potem przestala
jeszcze raz to
przetestowac czy sie counter bledow nie zeruje ??!!!
W miare juz jest OK, mozna
wrzucic na pi2
5. Astrometria w czasie scanu i jak za malo matchy
- Pokazal sie problem ze jesli astrometria tuz przed scanem sie nie
udala, to potem probowal ja powtarzac w czasie scanu az do skutku, w
tej chwili to zmienilem tak
ze jak astro jest wylaczona to i retry sie nie bedzie sie w kolko robil.
- Poza tym nie zawsze byla sprawdzana liczba matchy czy jest
dostatecznie duza zeby zaakceptowac ta astrometria ( ten warunek tez
przenioslem tak ze teraz sie to sprawdza zawsze )
- Po trzecie , jesli jestesmy w trakcie scanu i z jakis wzgledow
jednak astro sie zrobila to jesli liczba gwiazd nie przekracza limitu
to nie ma zadnej manipulacji modem migawki
na czas SCAN-u !
Luzne notatki z poprawek :
- TESTOWAC , opisac na WWW :
podobno scan 20070122 zrobil sie przy
otwartej migawce - czesciowo
k2a_070122_01456.fitc SHUTMODE = NORMAL
k2a_070122_01457.fitc SHUTMODE = NORMAL
k2a_070122_01458.fitc SHUTMODE = OPENED
k2a_070122_01459.fitc SHUTMODE = OPENED
w trakcie scanu robila sie astro - czemu ???
Robila sie astro ( ktora nie zbiegla przed
scanem ) i zrobila :
Changing shutter mode
to opened ...
Changed shutter mode to : 1
hmm czyli juz jasne, ale jak to ulepszyc ...
linia 213601 w /pi20/msok/data/20070122/log
20070123 tez to samo sie zrobilo
Rozwiazanie problemu - jak sie
astrometria robila w czasie scanu
bo poprzednia nie zbiegla - a w tle sie
powtarza do skutku to
wtedy otwiera migawke jak widzi ze jest duzo
gwiazd
sprawdzic - zeby jak astro juz wylaczona to
zeby nie powtarzal w
czasie scanu do skutku !!!
- TESTOWAC : sprawdzic to pmc - tabele , cos tam LWP
znalazl dziwnego
teraz jak za malo matchy to nie uznaje
takiej astrometrii
a wczesniej uznawal i chyba przez to
bylo zle
Astro kolo bieguna !!!