SOHO-2143

[Michal Kusiak]

Michale mógłbyś pomóc w pewnej sprawie. Chodzi o datę zdjęć na LASCO

np. te zdjęcie

Kod: [Zaznacz]

ftp://sohoftp.nascom.nasa.gov/qkl/lasco/quicklook/level_05/110930/c3/33290016.fts
w nagłówku ma między innymi:


FILENAME   = '33290016.fts'       /                                               
FILEORIG   = '110930_082956.img'  /                                               
DATE       = '2011/09/30 08:38:07.860' /                                           
DATE-OBS   = '2011/09/30'         /                                               
TIME-OBS   = '08:30:06.932'       /                                               


orginalny plik ma nazwę 110930_082956.img co sugerowało by godzinę wykonania zdjęcia na
08:29:56 co jest inne niż przyjęte 08:30:06.932

Czy to oznacza że początek ekspozycji miał miejsce o08:29:56 a koniec o 08:30:06.932 ?

Którą więc wartość uwzględnić w astrometri i czy to wogóle będzie miało znaczący wpływ na wyliczenie orbity ?


PS. Odnośnie zdjęć ze STEREO one mają w nagłówku pozycję sondy

dla    20110930_060901_s4h1A.fts

HCIX_OBS   = 118389458568. /                                               
HCIY_OBS   = 82441408544.4 /                                               
HCIZ_OBS   = -10802185912.3 /                                               
HAEX_OBS   = -51482583113.7 /                                               
HAEY_OBS   = 135199056712. /                                               
HAEZ_OBS   = -309295243.476 /                                               
HEEX_OBS   = -35594055773.9 /                                               
HEEY_OBS   = 140222347211. /                                               
HEEZ_OBS   = -305776906.501 /                                               
HEQX_OBS   = -34561067395.5 /                                               
HEQY_OBS   = 140064921959. /                                               
HEQZ_OBS   = -10802185912.3 /                                               

Ja osobiście korzystam z daty i godziny podawanej w plikach tekstowych na serwerze ze zdjęciami. Ekspozycja zdjęcia LASCO C3 wynosi między 19 a 20 sekund czyli środek ekspozycji wypada w 9-10 sekundzie.
Obierając podczas astrometrii czas można wybrać środek ekspozycji choć osobiście uważam że obranie za moment obserwacji, początku/ końca ekspozycji nie ma dużego znaczenia na wynik ponieważ:

 biorąc pod uwagę średnią szybkość komety Kreutza na dzisiaj  lub obiektów przelatujących w bliskiej odległości od Słońca
(5.5' kątowej / godzinę) = 330" na godzinę    / :3600 sek
v=0,09167" na sekundę

czyli różnica w położeniu pomiędzy środkiem i początkiem(lub koncem eksp) =~ 0,92"
roznica miedzy poczatkiem a koncem ekspozycji= ~1,84"

Zatem mieści się ona w błędzie dopasowania siatki (dopasowanie dobrze przeprowadzone daje błąd rzędu 1/6 czasem ok 1/10 piksela =5-10").

[Michal Kusiak]

Szymon,
Krótkie wyjaśnienie odnośnie P i Q masz również tutaj:
http://www.minorplanetcenter.net/iau/info/OrbElsExplanation.html

Linki które podał Adam, prezentują podstawy obliczeń elementów orbity w sposób analityczny.
Ażeby przeprowadzić dokładne obliczenia, obecnie trzeba skorzystać z programu komputerowego, gdyż na ruch ciała wpływa więcej ciał niż Słońce i wszelkie poprawki można jedynie przybliżyć w sposób numeryczny.

PS.Nie mogłem odpisać wcześniej, bo byłem w trasie..

[Limax7]

Ja osobiście korzystam z daty i godziny podawanej w plikach tekstowych na serwerze ze zdjęciami. Ekspozycja zdjęcia LASCO C3 wynosi między 19 a 20 sekund czyli środek ekspozycji wypada w 9-10 sekundzie.
Obierając podczas astrometrii czas można wybrać środek ekspozycji choć osobiście uważam że obranie za moment obserwacji, początku/ końca ekspozycji nie ma dużego znaczenia na wynik ponieważ:

 biorąc pod uwagę średnią szybkość komety Kreutza na dzisiaj  lub obiektów przelatujących w bliskiej odległości od Słońca
(5.5' kątowej / godzinę) = 330" na godzinę    / :3600 sek
v=0,09167" na sekundę

czyli różnica w położeniu pomiędzy środkiem i początkiem(lub koncem eksp) =~ 0,92"
roznica miedzy poczatkiem a koncem ekspozycji= ~1,84"

Zatem mieści się ona w błędzie dopasowania siatki (dopasowanie dobrze przeprowadzone daje błąd rzędu 1/6 czasem ok 1/10 piksela =5-10").

Dzięki.

Właśnie skończyłem astrometrię z STEREO-A HI1
XYZ podałem w km właśnie ze strony którą podałeś.
Z JPL HORIZON miałem problem.


ASTROMETRY:
    SOHO2143 4S2011 09 29.34990 19 56 34.8  -28 23 35                        C49
    SOHO2143 4s2011 09 29.34990 1 -198355500  +111291500  +47915470          C49
    SOHO2143 4S2011 09 29.43323 19 56 22.0  -28 11 44                        C49
    SOHO2143 4s2011 09 29.43323 1 -198535100  +111028700  +47801430          C49
    SOHO2143 4S2011 09 29.51656 19 56 10.0  -27 58 26                        C49
    SOHO2143 4s2011 09 29.51656 1 -198714300  +110765700  +47687270          C49
    SOHO2143 4S2011 09 29.59990 19 55 52.4  -27 46 04                        C49
    SOHO2143 4s2011 09 29.59990 1 -198893100  +110502400  +47573010          C49
    SOHO2143 4S2011 09 29.68323 19 55 37.4  -27 34 05                        C49
    SOHO2143 4s2011 09 29.68323 1 -199071400  +110238900  +47458640          C49
    SOHO2143 4S2011 09 29.76656 19 55 24.5  -27 21 59                        C49
    SOHO2143 4s2011 09 29.76656 1 -199249300  +109975100  +47344170          C49
    SOHO2143 4S2011 09 29.84990 19 55 08.4  -27 08 36                        C49
    SOHO2143 4s2011 09 29.84990 1 -199426900  +109711000  +47229580          C49
    SOHO2143 4S2011 09 29.93323 19 54 48.6  -26 55 50                        C49
    SOHO2143 4s2011 09 29.93323 1 -199603900  +109446700  +47114890          C49
    SOHO2143 4S2011 09 30.01656 19 54 34.0  -26 44 51                        C49
    SOHO2143 4s2011 09 30.01656 1 -199780600  +109182200  +47000080          C49
    SOHO2143 4S2011 09 30.09990 19 54 13.7  -26 30 35                        C49
    SOHO2143 4s2011 09 30.09990 1 -199956900  +108917400  +46885170          C49
    SOHO2143 4S2011 09 30.18323 19 53 50.7  -26 12 41                        C49
    SOHO2143 4s2011 09 30.18323 1 -200132700  +108652400  +46770150          C49
    SOHO2143 4S2011 09 30.26656 19 53 40.6  -26 03 33                        C49
    SOHO2143 4s2011 09 30.26656 1 -200308100  +108387100  +46655030          C49
    SOHO2143  S2011 09 30.35424 11 54 05.4  -04 22 51                        249
    SOHO2143  s2011 09 30.35424 2 -0.00876052 +0.00153026 +0.00051108        249
    SOHO2143  S2011 09 30.56315 11 58 02.6  -04 19 38                        249
    SOHO2143  s2011 09 30.56315 2 -0.00875808 +0.00147043 +0.00051498        249
    SOHO2143  S2011 09 30.65422 11 59 54.4  -04 15 03                        249
    SOHO2143  s2011 09 30.65422 2 -0.00875692 +0.00144432 +0.00051667        249
    SOHO2143  S2011 09 30.72924 12 01 26.0  -04 12 50                        249
    SOHO2143  s2011 09 30.72924 2 -0.00875591 +0.00142280 +0.00051806        249
    SOHO2143  S2011 10 01.02924 12 07 36.9  -03 59 47                        249
    SOHO2143  s2011 10 01.02924 2 -0.00875147 +0.00133665 +0.00052357        249
    SOHO2143  S2011 10 01.19589 12 11 30.1  -03 57 40                        249
    SOHO2143  s2011 10 01.19589 2 -0.00874872 +0.00128873 +0.00052659        249
    SOHO2143  S2011 10 01.43757 12 17 31.2  -03 44 21                        249
    SOHO2143  s2011 10 01.43757 2 -0.00874437 +0.00121917 +0.00053094        249


Po stopniowym wykluczaniu najsłabszych pomiarów
wychodzi że kometa miała okres obiegu ok 2.71-2.82 lat
Nie dodawałem żadnych stałych

Myślę, że wykonałem dobre pomiary choć trwało to dość długo
(ps. moje szkraby mi przeszkadzały chcąc bawić się z tatą )

Taki wynik ze wszystkimi 19 pomiarami SOHO i STEREO-A


SOHO2143
   Perihelion 2011 Oct 1.857517 TT = 20:34:49 (JD 2455836.357517)
Epoch 2011 Sep 29.0 TT = JDT 2455833.5   Earth MOID: 0.5898
M 359.12781              (2000.0)            P               Q
n   0.30522308     Peri.   90.68938      0.14928546     -0.98143837     
a   2.18469658     Node    11.60115     -0.96018222     -0.17296349     
e   0.9977780      Incl.  143.22746      0.23614392     -0.08283929     
P   3.23            q 726182.908km  Q 4.36453894
From 19 observations 2011 Sept. 29-Oct. 1;   RMS error 122.785 arcseconds


a taki po wykluczeniu 4 pomiarów pozostawiając 15


SOHO2143
   Perihelion 2011 Oct 1.853577 TT = 20:29:09 (JD 2455836.353577)
Epoch 2011 Sep 29.0 TT = JDT 2455833.5   Earth MOID: 0.5826
M 359.00431              (2000.0)            P               Q
n   0.34892579     Peri.   90.45496      0.14192257     -0.98360771     
a   1.99823854     Node    10.75019     -0.96209970     -0.16350109     
e   0.9976719      Incl.  143.39023      0.23285647     -0.07604777     
P   2.82            q 695929.281km  Q 3.99182508
From 15 observations 2011 Sept. 29-Oct. 1;   RMS error 46.139 arcseconds

[Michal Kusiak]

Część astrometrii z HI1a wygląda fajnie, gdyż błędy systematyczne wynikające z niedokładności optycznej przyrządu są znacząco mniejsze niż w przypadku LASCO.
Zatem Twój sposób jest na dobrej drodze do udoskonalenia.

[Limax7]

Już wiem gdzie miałem błąd generując XYZ w HORIZON

Trzeba było wybrać opcję w
reference plane :  Earth mean ...


Spróbuję więc na sam koniec zmienić to w astrometrii dla SOHO i STEREO
bo narazie jest mieszanka cześć z HORIZON część ze strony Stereo

[Szymon]

Dzięki Michał. Spróbuję coś policzyć.

[Limax7]

Jeszcze trochę pomęczę tę orbitę...

Zmieniłem odległości XYZ na tylko i WYŁĄCZNIE te z JPL HORIZON i jako jednostkę wziąłem km

oto astrometria:


ASTROMETRY:
    SOHO2143 4S2011 09 29.34990 19 56 34.8  -28 23 35                        C49
    SOHO2143 4s2011 09 29.34990 1 -198317013  +111310258  +47923906.6        C49
    SOHO2143 4S2011 09 29.43323 19 56 22.0  -28 11 44                        C49
    SOHO2143 4s2011 09 29.43323 1 -198496651  +111047538  +47809889.1        C49
    SOHO2143 4S2011 09 29.51656 19 56 10.0  -27 58 26                        C49
    SOHO2143 4s2011 09 29.51656 1 -198675869  +110784563  +47695762.1        C49
    SOHO2143 4S2011 09 29.59990 19 55 52.4  -27 46 04                        C49
    SOHO2143 4s2011 09 29.59990 1 -198854669  +110521334  +47581525.9        C49
    SOHO2143 4S2011 09 29.68323 19 55 37.4  -27 34 05                        C49
    SOHO2143 4s2011 09 29.68323 1 -199033049  +110257852  +47467180.7        C49
    SOHO2143 4S2011 09 29.76656 19 55 24.5  -27 21 59                        C49
    SOHO2143 4s2011 09 29.76656 1 -199211010  +109994118  +47352726.9        C49
    SOHO2143 4S2011 09 29.84990 19 55 08.4  -27 08 36                        C49
    SOHO2143 4s2011 09 29.84990 1 -199388550  +109730131  +47238164.5        C49
    SOHO2143 4S2011 09 29.93323 19 54 48.6  -26 55 50                        C49
    SOHO2143 4s2011 09 29.93323 1 -199565670  +109465892  +47123494.0        C49
    SOHO2143 4S2011 09 30.01656 19 54 34.0  -26 44 51                        C49
    SOHO2143 4s2011 09 30.01656 1 -199742369  +109201403  +47008715.5        C49
    SOHO2143 4S2011 09 30.09990 19 54 13.7  -26 30 35                        C49
    SOHO2143 4s2011 09 30.09990 1 -199918646  +108936662  +46893829.4        C49
    SOHO2143 4S2011 09 30.18323 19 53 50.7  -26 12 41                        C49
    SOHO2143 4s2011 09 30.18323 1 -200094502  +108671673  +46778835.7        C49
    SOHO2143 4S2011 09 30.26656 19 53 40.6  -26 03 33                        C49
    SOHO2143 4s2011 09 30.26656 1 -200269935  +108406434  +46663734.9        C49
    SOHO2143  S2011 09 30.35424 11 54 05.4  -04 22 51                        249
    SOHO2143  s2011 09 30.35424 1 -1310572.19 +179738.475 +161259.218        249
    SOHO2143  S2011 09 30.56315 11 58 02.6  -04 19 38                        249
    SOHO2143  s2011 09 30.56315 1 -1310209.07 +171317.393 +158242.567        249
    SOHO2143  S2011 09 30.65422 11 59 54.4  -04 15 03                        249
    SOHO2143  s2011 09 30.65422 1 -1310033.99 +167609.171 +156912.297        249
    SOHO2143  S2011 09 30.72924 12 01 26.0  -04 12 50                        249
    SOHO2143  s2011 09 30.72924 1 -1309883.80 +164573.896 +155822.572        249
    SOHO2143  S2011 10 01.02924 12 07 36.9  -03 59 47                        249
    SOHO2143  s2011 10 01.02924 1 -1309220.81 +152422.034 +15145.1984        249
    SOHO2143  S2011 10 01.19589 12 11 30.1  -03 57 40                        249
    SOHO2143  s2011 10 01.19589 1 -1308809.68 +145664.092 +149015.938        249
    SOHO2143  S2011 10 01.43757 12 17 31.2  -03 44 21                        249
    SOHO2143  s2011 10 01.43757 1 -1308159.60 +135857.171 +145473.816        249



a orbita z 15 pozycji zmieniła się deczko ...


SOHO2143
   Perihelion 2011 Oct 1.852545 TT = 20:27:39 (JD 2455836.352545)
Epoch 2011 Sep 30.0 TT = JDT 2455834.5   Earth MOID: 0.5901
M 359.37774              (2000.0)            P               Q
n   0.33588996     Peri.   92.85355      0.14070082     -0.97946020     
a   2.04961063     Node    13.77430     -0.95941355     -0.17089962     
e   0.9976956      Incl.  142.65499      0.24439497     -0.10700942     
P   2.93            q 706551.749km  Q 4.09449827
From 15 observations 2011 Sept. 29-Oct. 1;   RMS error 47.883 arcseconds


11 09 29.34990    C49 19 56 34.8    -28 23 35       6.7-  16.9-  .84002 .19638
11 09 29.43323    C49 19 56 22.0    -28 11 44       1.7+  42.2-  .84304 .19195
11 09 29.51656    C49 19 56 10.0    -27 58 26      27.0+  17.5+  .84610 .18746
11 09 29.59990    C49 19 55 52.4    -27 46 04      15.3-  18.9+  .84921 .18292
11 09 29.68323    C49 19 55 37.4    -27 34 05      16.2-   5.1-  .85237 .17833
11 09 29.76656    C49 19 55 24.5    -27 21 59      18.2+  24.9-  .85557 .17367
11 09 29.84990    C49 19 55 08.4    -27 08 36      17.8+  29.2+  .85882 .16894
11 09 29.93323    C49 19 54 48.6    -26 55 50      23.8-  43.0+  .86212 .16415
11 09 30.01656    C49 19 54 34.0    -26 44 51      12.9+  53.8-  .86547 .15928
11 09 30.09990    C49 19 54 13.7    -26 30 35      17.2-  42.3+  .86889 .15434
11 09 30.18323  X C49 19 53 50.7    -26 12 41      73.9-   352+  .87237 .14932
11 09 30.26656  X C49 19 53 40.6    -26 03 33      54.0+   131+  .87591 .14422
11 09 30.35424  X 249 11 54 05.4    -04 22 51      20.8+   118+  .95893 .13877
11 09 30.56315    249 11 58 02.6    -04 19 38       6.4-  38.5-  .95857 .12527
11 09 30.65422    249 11 59 54.4    -04 15 03      56.2+  71.7+  .95855 .11915
11 09 30.72924  X 249 12 01 26.0    -04 12 50      74.1+  62.6+  .95861 .11399
11 10 01.02924    249 12 07 36.9    -03 59 47      66.9-  14.4+  .95958 .09207
11 10 01.19589    249 12 11 30.1    -03 57 40      13.9+  75.0-  .96095 .07873
11 10 01.43757    249 12 17 31.2    -03 44 21       4.2+  18.4+  .96435 .05715

Station data:
(249) SOHO.
(C49) STEREO-A.

[Limax7]

Astrometrię robię w "Astrometrica Software" Herberta Raaba i (jeśli są kłopoty z wyznaczeniem odpowiedniej siatki gwiazd) w takim prywatnym programiku napisanym przez kolegę.
Sęk leży w tym ażeby tak jak w przypadku normalnych zdjęć, odpowiednio wyskalować i nałożyć siatkę na zdjęcia LASCO, HI1a lub COR2 (koniecznie ftsy). Odpowiednio -> program musi dopasować minimum 25 gwiazd ze zdjęcia i błąd nałożenia nie może być większy niż 10 sekund łuku (1/6 piksela C3). W przypadku C3 ja staram się uzyskiwać siatkę z minimum 40-45 gwiazdami referencyjnymi, a jeśli się to uda, znacząco poprawia dokładność wyznaczonej pozycji. Warto zauważyć że LASCO ma wyraźną niesymetryczną dystorsję ( zakrzywione pole widzenia) stąd najlepiej nakładanie siatki wykonywać na gwiazdach znajdujących się w "ćwiartce" zdjęcia gdzie znajduje się nasz obiekt

W przypadku COR2 i HI (tak samo dla LASCO) trzeba za pomocą jakiegoś oprogramowania (chocby jpl.horizons) odnalezc obszar, który akurat znajduje się w polu widzenia instrumentów. Potem odszukać jakiś punkt odniesienia o znanej współrzędnej (jasna gwiazda lub seria gwiazd) i nałożyć siatkę.

Gdy mamy współrzędne trzeba następnie wyznaczyć wektory X,Y,Z dla sond, z których posiadamy wyznaczone pozycje naszej komety. Wprowadzić je do danych z astrometrii i w Find Orbie, Solex/ExOrb policzyc orbite.

Zabawa bardzo czasochłonna, często nieustanna walka z szumem i właściwą kalibracją zdjęcia a potem z odpowiednim dopasowaniem.

Wczoraj wieczorem zmierzyłem się z Astrometricą. Ohh ciężka sprawa! Najtrudniej skonfigurować program pod względem ogniskowej i PA oraz znaleźć gwiazdy referencyjne i je dopasować. Nie wiem też jaki katalog gwiazd wybrać do pobierania i jaki zakres gwiazd wybrać do pobrania..

[Michal Kusiak]

Szymon Liwo zgłosił dzisiaj kometę która jutro może okazać się interesującym obiektem.
Obecnie w elongacji 6.5 stopnia z jasnością około 8.5 magnitudo. Jutro spodziewam się że w polu C3  (jeszcze przed wejściem w zasięg C2) uzyska między 4 a 5 magnitudo. Z kolei w polu LASCO C2 może okazać się nieco jaśniejsza.
Brawo Szymon!

[Szymon]

Dzięki Michał. Bacznie obserwuję jak moja kometa leci ku Słońcu. Jest już jaśniejsza i ma warkocz.

[kanarkusmaximus]

Dzięki Michał. Bacznie obserwuję jak moja kometa leci ku Słońcu. Jest już jaśniejsza i ma warkocz.

O! Brawa i gratulacje dla Szymona !

[Szymon]

Dzięki Kanarkusmaximus. To już dziewiąta kometa na moim koncie 

[Matias]

Dzięki Kanarkusmaximus. To już dziewiąta kometa na moim koncie 

Łał, gratulacje!!

[Michal Kusiak]

Na chwilę obecną:

[Rafal Reszelewski]

Kometa Szymona na C2