Polskie Forum Astronautyczne
Astronautyka => Słońce => Komety SOHO i STEREO => Wątek zaczęty przez: Michal Kusiak w Października 01, 2011, 20:22
-
Cześć,
Wczorajszy poranek przejdzie do mojej pamięci jako jeden z najbardziej wyjątkowych.
30 września podczas porannej transmisji danych z sondy SOHO pobrałem pierwsze zdjęcia uzyskane przez koronograf LASCO C3.
Już na pierwszym zdjęciu z rozpoczętej transmisji pojawił się dreszczyk, ponieważ tuż przy brzegu pola widzenia (elongacja 8 stopni) znajdowała się jasna kometa mająca wówczas 6 magnitudo i krótki warkocz o długości kilku minut łuku.
W napięciu oczekiwałem na kolejne dwie fotografie i przyszedł ten moment. :) Bach! O godzinie 10:38:43 (CEST) zgłosiłem kometę mając wymagane 3 pozycje do zgloszenia. Okazało się jednak, że oprócz mnie obiekt zgłosili trzej kolejni łowcy Bo Zhou-Chiny, Szymon Liwo i Zhijan Xu-Chiny i co ciekawe wszystkie zgłoszenia spłynęły w czasie zaledwie 9 sekund, co pozwoliło na przyznanie współodkrycia każdemu ze zgłaszających.
Wstępna analiza nowoodkrytej komety wskazywała na maksimum blasku około +1 mag, rankiem lub w okolicy południa 1 października, a policzona przeze mnie orbita ewidentnie pokazała że mamy do czynienia z kometą Kreutza której peryhelium (~0.0057 AU) wypada 1 października około 23 CEST.
Jak się okazało byłem w błędzie co do jasności, co mnie bardzo ucieszyło. Dlaczego? Otóż kometa podczas zbliżania się do Słońca prawdopodobnie doznała niewielkiego wybuchu, co przeniosło się na uwolnienie znacznie większej ilości pyłu i ewentualnych molekuł (ewentualnych, ponieważ w tak bliskiej okolicy Słońca dochodzi do silnej fotodysocjacji która powoduje b. szybki ich rozpad). Kometa już na zdjęciach porannych uzyskała 1.0 mag, wybuch spowodował jednak wzrost blasku do conajmniej -1 magnitudo, co uczyniło SOHO-2143 (bo tak została oznaczona) najjaśniejszą kometą Kreutza zaobserwowaną przez cały okres misji sondy SOHO.
(http://img593.imageshack.us/img593/3285/201110011442c3.gif)
Nietypowe przepalenie komety pochodzi w większości od pyłu emitowanego przez drobne fragmenty pozostałych po rozpadzie jądra
(http://img829.imageshack.us/img829/3575/1448.gif)
SOHO-2143 na C2 (14:48UT). Zgrubienie w warkoczu to najprowdopodobniej chmura pyłu i drobnych fragmentów komety
Liczę na to że to nie ostatni tak jasny obiekt.
Roczna liczba odkryć komet z grupy Kreutza systematycznie wzrasta, co mam nadzieję, sugeruje przeloty kolejnych tak jasnych, lub jaśniejszych fragmentów.
Michał
-
Gratulacje! Świetny obiekt. :)
-
Ano, bardzo ładny ślad. Gratulacje.
-
Ależ jasna ! Gratuluję.
Ciekawy jestem zdjęć ze STEREO, czy też ją mają ?
-
Na STEREO też jest,
Alan Watson zaobserwował kometę na zdjęciach HI-1a z 28 września.
Poza tym jutro / pojutrze powinny pojawić się zdjęcia z COR2 w pełnym formacie i tam też na pewno będzie ją widać.
MK
-
Na STEREO też jest,
Alan Watson zaobserwował kometę na zdjęciach HI-1a z 28 września.
Poza tym jutro / pojutrze powinny pojawić się zdjęcia z COR2 w pełnym formacie i tam też na pewno będzie ją widać.
MK
Dzięki :)
Widzę ją w prawym dolnym rogu na HI-1A ale z 29 września
http://stereo-ssc.nascom.nasa.gov/browse/2011/09/29/ahead_20110929_hi1_512.mpg
http://stereo-ssc.nascom.nasa.gov/browse/2011/09/29/index.shtml
PS. Skąd pobierasz aktualniejsze zdjęcia z C2 i C3 niż te na
http://sohodata.nascom.nasa.gov/cgi-bin/data_query
-
Widać ją nieźle! :)
Zdjęcia C3,C2 pobieram z poniższych źródeł:
http://lasco-www.nrl.navy.mil/javagif/gifs/
albo w postaci FTS:
ftp://sohoftp.nascom.nasa.gov/qkl/lasco/quicklook/level_05/
-
Gratulacje :)
-
Jestem pod wrażeniem, gratulacje!
-
Dziś pojawiła się reszta zdjęć z LASCO z widoczną kometą (wczoraj była przerwa w dostawie zdjęć)
Zrobiłem więc animację
http://www.youtube.com/watch?v=a0D0ttMTN4s
-
Świetna robota Adam!
W ruchu widać pięknie jak dynamicznie przebiegał rozpad jądra.
-
...
Poza tym jutro / pojutrze powinny pojawić się zdjęcia z COR2 w pełnym formacie i tam też na pewno będzie ją widać.
No i są już. Animacja na stronie STEREO jest mocno skompresowana więc utworzyłem po swojemu ::)
http://www.youtube.com/watch?v=LLyV9Kd22sc
Chociaż i tu Youtube skompresował po swojemu :-\
-
Niemniej wszyło bardzo fajnie ^^
-
Zabawne, że zaraz po zniknięciu komety za przesłoną, następuje CME z drugiej strony tarczy :).
-
Ciekaw jestem czy kometa załapała się na zdjęciach z SDO AIA 171 podobnie jak ta z 5/6 lipca 2011 ???
http://sdo.gsfc.nasa.gov/gallery/potw.php?v=item&id=58
Ponżej screen z SOHO dla porównania
-
Da się jakkolwiek oszacować masę i rozmiary tego obiektu?
-
Czy zna ktoś rzetelne źródło ze zdjęciami AIA171? Z tego co pamiętam zdjęcia wykonywane są co ok. 5 minut, ale niestety na GSFC dane są niepełne.
Policzyłem orbitę dla komety SOHO-2143 używając LASCO C3 oraz STEREO: COR2B i HI1A. Warto byłoby zatem spróbować odszukać ten obiekt, choć zakładam, że szansa jest nieduża, ponieważ jeszcze przed wejściem w koronograf C2 (elongacja ok 2 stopnie) kometa zaczęła dość szybko się rozpadać i prawdopodobnie jądro było już w postaci chmury drobnych fragmentów.
Próby odnalezienia tej komety rozpocznie niebawem także Alan Hale, bazując na poniższym przeliczeniu.
Orbita odrzuca kilka pozycji które w znaczący sposób zmieniały wynik, co więcej zastosowałem przeliczenie wagowe przyjmując wagę pozycji z COR 2B=3.7, LASCO C3=1.0, HI1a=0.8 (waga=rozdzielczość piksela LASCO/rozdzielczości piksela instrumentu)
Orbital elements:
SOHO2143
Perihelion 2011 Oct 1.852270 TT = 20:27:16 (JD 2455836.352270)
Epoch 2011 Sep 30.0 TT = JDT 2455834.5 Earth MOID: 0.5607
M 359.985 (2000.0) P Q
n 0.007755 Peri. 84.912 0.172652 -0.983140
a 25.27709 Node 5.969 -0.961986 -0.181440
e 0.99979 Incl. 144.614 0.211597 -0.022689
P 127.08 q 784011.7km Q 50.54894
From 50 observations 2011 Sept. 29-Oct. 1; RMS error 16.948 arcseconds
11 09 29.56182 C49 19 56 04.9 -27 51 11 34.0+ 13.7+ .84558 .18615
11 09 29.61737 C49 19 55 52.7 -27 42 47 1.4+ 24.7+ .84776 .18306
11 09 29.64515 X C49 19 55 48.0 -27 38 03 4.7+ 61.7+ .84885 .18150
11 09 29.70071 X C49 19 55 37.2 -27 29 44 4.9- 66.1+ .85106 .17837
11 09 29.75626 C49 19 55 24.8 -27 22 07 32.6- 27.4+ .85329 .17521
11 09 29.78404 X C49 19 55 22.2 -27 17 21 2.7+ 65.1+ .85442 .17362
11 09 29.81182 C49 19 55 17.0 -27 13 29 4.4+ 48.5+ .85555 .17202
11 09 29.89515 C49 19 55 00.6 -27 00 58 3.6+ 51.9+ .85897 .16718
11 09 29.97848 X C49 19 54 42.5 -26 48 08 11.6- 71.1+ .86244 .16228
11 09 30.00626 X C49 19 54 36.1 -26 43 54 19.6- 74.1+ .86361 .16063
11 09 30.26734 C50 05 50 53.66 +19 49 23.0 4.2+ 6.7- 1.1979 .14480
11 09 30.27091 X 249 11 52 20.7 -04 26 10 7.1- 50.8- .95971 .14449
11 09 30.32090 249 11 53 17.8 -04 24 16 6.7+ 2.3- .95951 .14134
11 09 30.33757 249 11 53 35.1 -04 24 10 16.7- 18.5- .95944 .14028
11 09 30.35067 C50 05 51 08.75 +19 57 34.4 1.2+ 4.2- 1.1934 .13953
11 09 30.35425 249 11 53 55.3 -04 23 44 2.6+ 15.0- .95938 .13921
11 09 30.36280 249 11 54 03.5 -04 23 24 20.2- 6.6- .95935 .13867
11 09 30.37924 249 11 54 22.1 -04 23 07 22.3- 12.1- .95930 .13761
11 09 30.39590 249 11 54 41.2 -04 22 59 21.5- 27.1- .95924 .13654
11 09 30.43400 C50 05 51 24.49 +20 05 57.7 2.2+ 2.8- 1.1887 .13416
11 09 30.47567 C50 05 51 32.42 +20 10 15.3 1.4+ 1.4- 1.1864 .13144
11 09 30.47924 249 11 56 18.1 -04 20 46 12.0- 12.7- .95901 .13112
11 09 30.49591 249 11 56 38.8 -04 20 14 6.6+ 5.1- .95898 .13002
11 09 30.51734 C50 05 51 40.47 +20 14 35.8 .78+ .43- 1.1840 .12869
11 09 30.52090 249 11 57 08.1 -04 19 31 6.8+ .89+ .95893 .12836
11 09 30.53817 C50 05 51 44.53 +20 16 45.9 .37+ 1.8- 1.1828 .12730
11 09 30.55900 C50 05 51 48.63 +20 18 59.7 .03+ .51+ 1.1816 .12591
11 09 30.56315 249 11 57 57.4 -04 18 37 1.8- 9.0- .95885 .12554
11 09 30.57984 C50 05 51 52.76 +20 21 13.3 .25- .04+ 1.1804 .12451
11 09 30.60067 C50 05 51 56.94 +20 23 29.2 .34- 1.7+ 1.1791 .12310
11 09 30.62924 249 11 59 16.6 -04 16 52 1.8+ 7.4- .95878 .12107
11 09 30.64234 C50 05 52 05.35 +20 28 01.6 1.0- 2.5+ 1.1767 .12026
11 09 30.64591 249 11 59 36.7 -04 16 21 1.8+ 3.1- .95877 .11993
11 09 30.66257 249 11 59 57.1 -04 15 55 5.2+ 4.1- .95876 .11878
11 09 30.66317 C50 05 52 09.61 +20 30 20.2 1.4- 3.7+ 1.1754 .11883
11 09 30.68400 C50 05 52 14.01 +20 32 38.1 .21- 3.2+ 1.1742 .11738
11 09 30.68757 249 12 00 27.0 -04 15 17 2.5- 7.2- .95876 .11705
11 09 30.70424 249 12 00 47.6 -04 14 47 .76+ 4.9- .95876 .11589
11 09 30.72567 C50 05 52 22.65 +20 37 20.4 1.7- 5.1+ 1.1717 .11448
11 09 30.72925 249 12 01 17.8 -04 14 03 7.3- 3.1- .95877 .11414
11 09 30.74590 249 12 01 39.2 -04 13 36 5.2+ 4.6- .95878 .11297
11 09 30.76734 C50 05 52 31.55 +20 42 05.5 1.8- 5.1+ 1.1691 .11154
11 09 30.77091 249 12 02 10.0 -04 12 47 1.3+ 1.1+ .95880 .11120
11 09 30.78758 249 12 02 30.9 -04 12 20 2.6+ 1.2- .95882 .11001
11 09 30.81257 249 12 03 02.0 -04 11 38 1.4- 3.8- .95886 .10821
11 09 30.86257 249 12 04 05.9 -04 09 59 6.0+ 3.9+ .95896 .10458
11 09 30.88306 249 12 04 32.1 -04 09 21 5.5+ 3.5+ .95901 .10308
11 09 30.88758 249 12 04 37.6 -04 09 17 .86+ 1.1- .95902 .10274
11 09 30.89590 249 12 04 48.6 -04 09 03 5.2+ 2.8- .95904 .10213
11 09 30.90424 249 12 04 59.4 -04 08 38 5.7+ 6.2+ .95907 .10151
11 09 30.91257 249 12 05 10.3 -04 08 26 7.6+ 2.2+ .95909 .10089
11 09 30.92093 249 12 05 21.0 -04 08 09 5.5+ 3.1+ .95912 .10026
11 09 30.92924 249 12 05 31.5 -04 07 56 .94+ .06- .95915 .09964
11 09 30.93757 249 12 05 43.0 -04 07 34 10.6+ 5.7+ .95918 .09902
11 10 01.18757 249 12 11 21.5 -03 58 34 14.0- 5.1+ .96071 .07927
11 10 01.22925 X 249 12 12 19.7 -03 56 50 39.0- 7.5+ .96111 .07575
11 10 01.27091 249 12 13 21.9 -03 55 14 20.2- 2.1- .96156 .07216
11 10 01.31257 X 249 12 14 20.9 -03 53 30 67.3- 8.0- .96208 .06849
Station data:
(249) SOHO.
(C49) STEREO-A.
(C50) STEREO-B.
-
Da się jakkolwiek oszacować masę i rozmiary tego obiektu?
Według pracy doktorskiej M.Knight "Dissertation" można oszacować rozmiar. Według jego przeliczeń największe fragmenty Kreutza obserwowane na SOHO miały średnicę do ok. 80 metrów.
-
Według pracy doktorskiej M.Knight "Dissertation" można oszacować rozmiar. Według jego przeliczeń największe fragmenty Kreutza obserwowane na SOHO miały średnicę do ok. 80 metrów.
A to oznacza marne szanse na misje bezzałogową. ;)
-
Czy zna ktoś rzetelne źródło ze zdjęciami AIA171? Z tego co pamiętam zdjęcia wykonywane są co ok. 5 minut, ale niestety na GSFC dane są niepełne.
Ja korzystam z :
http://aia.lmsal.com
a konkretnie ze zdjęć JPEG2000
http://jsoc.stanford.edu/data/aia/images
Ciekawą stroną jest też ta:
Guide to SDO Data Analysis
https://www.lmsal.com/sdodocs/doc/dcur/SDOD0060.zip/zip/entry/index.html
-
Michale mógłbyś opisać jak wyznaczasz pozycje RA i Dec dla SOHO (LASCO C2 i C3), STEREO (COR2, HI1)
by potem móc wyznaczyć orbitę komety ?
-
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.
-
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.
XYZ chyba można wygenerować w JPL HORIZON lub Solexem ale jak wprowadzić położenie sondy do astrometrii ??
Mógłbyś to rozwinąć.
-
Zerknij na poniższą ścieżkę:
WATSON02 S2010 03 11.57100 23 55 46.3 -04 54 44 EC038249
WATSON02 s2010 03 11.57100 2 +0.00873605 -0.00528850 -0.00257565 EC038249
WATSON02 S2010 03 11.59598 23 55 28.3 -04 53 01 EC038249
WATSON02 s2010 03 11.59598 2 +0.00873752 -0.00528521 -0.00257343 EC038249
WATSON02 S2010 03 11.61265 23 55 17.6 -04 51 26 EC038249
WATSON02 s2010 03 11.61265 2 +0.00873825 -0.00528356 -0.00257231 EC038249
Miałem na myśli:
Jak widzisz mamy po dwa wiersze dla danej daty/momentu obserwacji, pierwsza określa wyznaczone pozycje. Druga z kolei określa wektory X,Y,Z.
Żeby program obliczający orbitę odpowiednio zinterpretował je w obliczeniach trzeba dodać
- przed wierszem określającym pozycję komety/planetoidy duże "S"
- przed wierszem z wektorem małe "s".
Program interpretuje wówczas jako obserwacje wykonane przez satelitę, określając przy tym wektor zaczepiony w odpowiednim punkcie odniesienia.
Ja korzystam ze źródeł z wektorami znajdującymi się na serwerach:
ftp://sohoftp.nascom.nasa.gov/pub/data/ancillary/orbit/predictive/
http://www.srl.caltech.edu/STEREO2/Position/
-
Ja korzystam ze źródeł z wektorami znajdującymi się na serwerach:
ftp://sohoftp.nascom.nasa.gov/pub/data/ancillary/orbit/predictive/
http://www.srl.caltech.edu/STEREO2/Position/
Dane z tej strony http://www.srl.caltech.edu/STEREO2/Position
rozumiem.
Ale z ftp://sohoftp.nascom.nasa.gov/pub/data/ancillary/orbit/predictive/nie mogę rozgryźć ...
-
W plikach dat masz nastepujące dane które bierzesz pod uwagę:
Date Time DOY X Y Z
21-Sep-2011 00:00:00.000 2011 264 00000000 -1279359.16436 535138.16950 282857.06930
Dane X,Y,Z masz podane w kilometrach.
Źródło:
ftp://sohoftp.nascom.nasa.gov/pub/data/ancillary/orbit/predictive/SO_OR_PRE_20110921_V01.DAT
-
W plikach dat masz nastepujące dane które bierzesz pod uwagę:
Date Time DOY X Y Z
21-Sep-2011 00:00:00.000 2011 264 00000000 -1279359.16436 535138.16950 282857.06930
Dane X,Y,Z masz podane w kilometrach.
Źródło:
ftp://sohoftp.nascom.nasa.gov/pub/data/ancillary/orbit/predictive/SO_OR_PRE_20110921_V01.DAT
Ale np. weźmy godzinę 12:00:00
21-Sep-2011 12:00:00.000 2011 264 43200000 -1282984.51389 517815.42649 277247.15270
rozumiem że 264 jest dniem w roku ale zamiast 43200000 nie powinno być 50000000 ?
Co właściwie oznacza w tym przypadku 43200000 ?
Aha i zobaczyłem że są też ***V02.DAT
tam są inne liczby x y z i pozostałe
Czym więc różni się SO_OR_PRE_20110921_V01.DAT od SO_OR_PRE_20110921_V02.DAT
-
różnice między v01 a v02 są zwykle delikatne (nie przekraczają 10 kilometrów co przeklada sie na blad pomiaru pozycji mieszczący się dokladnosci przyrządu), ale jeśli jest v02 dostępny wówczas wybierz go jako za właściwy gdyż zawiera korektę i jesli tylko można, warto pozbywać się systematycznych błędów.
43200000 oznacza: 43200,000 sekundę dnia ;)
-
Michale, dzięki za wyjaśnienia.
Ja mam pewien pomysł na wyznaczanie astrometrii komet z SOHO i STEREO...
Mianowicie do Guide8 można załadować pozycje sond właśnie w oparciu o XYZ. Sam mam tak
kilka sond i księżyców Saturna zaimplementowanych. Jest też chyba program do konwersji
plików SPICE (*.bsp) do formatu XYZ
W Guide mogę za pomocą " Shift + < " znaleźć się na wybranej sondzi, planetoidzie etc. i
spojrzeć w przestrzeń z tej perspektywy.
Często korzystam z tego analizując zdjęcia z HI 1
Więc nachodzi mi na myśl poprostu zaznaczanie kursorem miejsca gdzie znajduje się kometa w danym czasie i odczytanie
RA i Dec z pozycji sondy. Może nawet da się wgrać do Guide ów FITSy jako tło.
Nie wiem jak z dokładnością tej metody, jak będę miał czas to zmierzę się z tym dla HI 1-A i LASCO C3 dla komety w tym poście.
Wyliczę orbitę FIND_ORBitem
-
Mały przykład. Poniżej zdjęcie z SOHO 30.09.2011 godz. 08:30
i ten sam widok z SOHO wygenerowany Guidem
-
Możesz spróbować, choć osobiście obawiam się, że astrometria wykonana taka metodą będzie miała dość słabą jakość. Przede wszystkim, Guide nie jest przeznaczony do wykonywania precyzyjnych pomiarów pozycji.
Pozycja powinna być mierzona do wartości subpikselowych. Co więcej tak jak wspomniałem zdjęcia pochodzące z koronografów LASCO są obarczone niesymetryczną dystorsją, stąd możesz mieć problem z nałożeniem na siebie mapy i zdjęcia.
Ciekaw jestem, jaki uzyskasz rezultat, zatem oczekuję na wieści. :)
-
Problemem może być zbyt duża jasność komety (zwłaszcza w miarę zbliżania się do Słońca)
I tu nawet FITS nie pomoże.
Dlatego najpierw spróbuję LASCO C3 kilka pozycji zaraz po wejściu w pole widzenia...
-
Oczywiście. Przepalenie komety bądź jasny warkocz mogą wpływać na wyznaczone pozycje.
W przypadku tak dużego przepalenia jak u SOHO2143 niestety wyznaczenie pozycji jest już raczej niemożliwe.
Godzinę temu,, znalazłem kolejnego ciekawego Kreutza na zdjęciach C3 tuż przy brzegu pola widzenia (z jasnością 9-9.5 mag). Niezależne zgłoszenie nadesłał Sergei Schmalz (21 sekund później). Kometa jest duża! Spodziewam się jednak, że jutro około 17UT gdy kometa będąc 2-3 stopnie od Słońca, uzyska około 4 może 3 mag. Konfiguracja orbity jest w tej chwili dość niekorzystna (duży kąt pozycyjny) co przekłada się na zjawisko "rozpraszania w tył" (backward scattering), które w rezultacie przyciemnia kometę o jakieś 1.5-2 mag.
Zobaczymy zresztą co z kometą się stanie... warto jutro monitorować zdjęcia.
-
Z 7 pozycji i tylko LASCO C3
*******************************************************************************
2455834.854236111 = A.D. 2011-Sep-30 08:30:06.0000 (CT)
X =-8.760517076678498E-03 Y = 1.530257377144018E-03 Z = 5.110797637666262E-04
2455835.063148148 = A.D. 2011-Sep-30 13:30:56.0000 (CT)
X =-8.758079709826233E-03 Y = 1.470425508106596E-03 Z = 5.149809080274550E-04
2455835.154224537 = A.D. 2011-Sep-30 15:42:05.0000 (CT)
X =-8.756915411014149E-03 Y = 1.444315743935875E-03 Z = 5.166714730108776E-04
2455835.229236111 = A.D. 2011-Sep-30 17:30:06.0000 (CT)
X =-8.755910247609503E-03 Y = 1.422800136411107E-03 Z = 5.180591472693902E-04
2455835.529236111 = A.D. 2011-Oct-01 00:42:06.0000 (CT)
X =-8.751474299211769E-03 Y = 1.336653914930227E-03 Z = 5.235659050753672E-04
2455835.695891204 = A.D. 2011-Oct-01 04:42:05.0000 (CT)
X =-8.748723956085032E-03 Y = 1.288734755455538E-03 Z = 5.265946657921168E-04
2455835.937569445 = A.D. 2011-Oct-01 10:30:06.0000 (CT)
X =-8.744374921359654E-03 Y = 1.219169132553775E-03 Z = 5.309474615634443E-04
*******************************************************************************
ASTROMETRY:
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
I za nic nie mogę obliczyć orbity...
Albo źle wyznaczone pozycje, albo źle wygenerowana orbita SOHO (XYZ) w HORIZONS
albo jedno i drugie.
Bardzo też możliwe że za mało wyznaczonych pozycji..
Może warto zamienić XYZ z HORIZON na te z
ftp://sohoftp.nascom.nasa.gov/pub/data/ancillary/orbit/predictive/
-
Ja uzyskałem coś takiego:
Orbital elements:
SOHO2143
Perihelion 2011 Oct 1.802442 TT; Constraint: e=1,q=0.005
Epoch 2011 Sep 30.0 TT = JDT 2455834.5 Earth MOID: 0.5413
M 359.999 (2000.0) P Q
n 0.000043 Peri. 74.412 0.369477 -0.926000
a 802.9261 Node 7.596 -0.901635 -0.377435
e 0.99999 Incl. 144.092 0.224811 0.008125
P 22751 q 747871.5km Q 1605.847
From 4 observations 2011 Sept. 30-Oct. 1 (15.2 hr); RMS error 73.266 arcseconds
11 09 30.35424 X 249 11 54 05.4 -04 22 51 122- 33.1+ .98847 .13632
11 09 30.56315 249 11 58 02.6 -04 19 38 63.8- 79.8- .98504 .12247
11 09 30.65422 249 11 59 54.4 -04 15 03 27.2+ 48.4+ .98363 .11619
11 09 30.72924 249 12 01 26.0 -04 12 50 63.7+ 53.4+ .98252 .11089
11 10 01.02924 X 249 12 07 36.9 -03 59 47 59.3- 248+ .97859 .08835
11 10 01.19589 249 12 11 30.1 -03 57 40 23.7- 22.0- .97690 .07458
11 10 01.43757 X 249 12 17 31.2 -03 44 21 229- 63.1+ .97556 .05214
Station data:
(249) SOHO.
Jak widzisz, błędy są dość spore i pomimo odrzucenia najgorszych obserwacji przekraczają katowy rozmiar 1 piksela.
-
A jednak coś uzyskałem pozostawiając 5 i 4 obserwacje...
Jako stałą podałem tylko e=1
SOHO2143
Perihelion 2011 Oct 2.239246 TT; Constraint: e=1
Epoch 2011 Oct 1.0 TT = JDT 2455835.5 Earth MOID: 0.4050
q 0.031589451 (2000.0) P Q
Peri. 59.895694 0.593017895 -0.801093425
Node 7.912140 -0.774616360 -0.595094297
e 1.0 Incl. 143.896824 0.219770947 0.064125666
From 5 observations 2011 Sept. 30-Oct. 1 (20.2 hr); RMS error 53.237 arcseconds
11 09 30.35424 249 11 54 05.4 -04 22 51 11.6+ 30.8+ .94143 .14321
11 09 30.56315 249 11 58 02.6 -04 19 38 40.6- 93.5- .93869 .13117
11 09 30.65422 249 11 59 54.4 -04 15 03 10.6+ 31.3+ .93766 .12581
11 09 30.72924 249 12 01 26.0 -04 12 50 19.9+ 34.6+ .93690 .12133
11 10 01.02924 X 249 12 07 36.9 -03 59 47 128- 239+ .93483 .10285
11 10 01.19589 249 12 11 30.1 -03 57 40 1.8- 3.2- .93453 .09218
Station data:
(249) SOHO.
i z 4 obserwacji
SOHO2143
Perihelion 2011 Oct 2.279271 TT; Constraint: e=1
Epoch 2011 Oct 1.0 TT = JDT 2455835.5 Earth MOID: 0.3909
q 0.034010081 (2000.0) P Q
Peri. 58.906815 0.604658280 -0.792580327
Node 7.704059 -0.766621290 -0.605961008
e 1.0 Incl. 144.013678 0.216055921 0.068027053
From 4 observations 2011 Sept. 30-Oct. 1 (20.2 hr); RMS error 7.629 arcseconds
11 09 30.35424 249 11 54 05.4 -04 22 51 1.7+ 7.5- .93810 .14415
11 09 30.56315 X 249 11 58 02.6 -04 19 38 55.8- 125- .93542 .13228
11 09 30.65422 249 11 59 54.4 -04 15 03 5.8- 3.2+ .93443 .12700
11 09 30.72924 249 12 01 26.0 -04 12 50 3.2+ 9.5+ .93370 .12259
11 10 01.02924 X 249 12 07 36.9 -03 59 47 137- 227+ .93175 .10447
11 10 01.19589 249 12 11 30.1 -03 57 40 .92+ 5.4- .93152 .09403
Station data:
(249) SOHO.
Spróbuję dodać astrometrię z STEREO-A HI1
-
Przyjrzałem się Twoim wynikom i faktycznie orbita zbliżona jest do grupy Kreutza. Fizycznie/dynamicznie jednak nie odpowiada wymaganiom. Mądrzejsi ode mnie opisali niegdyś genezę powstania grupy Kreutza. Mam nadzieję, że również Tobie się przyda trochę wiedzy z tego tematu.
Skrótowo mówiąc, różnice w wartościach peryhelium i podział na dwie subgrupy jest efektem rozpadu macierzystej komety/superfragmentów w wyniku sił pływowych i dalszych rozpadów zachodzących już niepływowo (np. rotacyjnie). W momencie rozpadu fragmenty dostają dodatkowy zastrzyk energii i pędu. Jak może słyszałeś elementy orbity dzielą się na dynamiczne, geometryczne i czasowe.
Do dynamicznych należą właśnie e oraz q stąd na nie rozpad działa najbardziej. Zmiana czasowego bywa następstwem dynamicznego poprzez ewentualną ingerencję kształtu orbity i wydłużenie bądź skrócenie okresu obiegu.
Różnice w elementach geometrycznych są jakby to nazwać "artefaktem oddziaływań Jowisza na powstałe fragmenty". Słońce to dyrygent, a Jowisz odgrywa "pierwsze skrzypce" :), a skutkiem tego wszystkiego jest dość ciekawa ewolucja grupy i bardzo ciekawa "stosunkowo liniowa" zależność pomiędzy parametrami Peri i Node.
(http://img218.imageshack.us/img218/16/kreutzrb7.gif)
Głównym podziałem na dwie subgrupy jest granica w wartości peryhelium (q~0.007)
Zauważ, że na poniższym obrazku próbka danych wyraźnie wskazuje że koncentracja fragmentów z grupy Kreutza przeważa w subgrupie 1, której peryhelium jest mniejsze od 0.007 AU.
(http://img149.imageshack.us/img149/4327/nodell1.gif)
-
Popełniłem błąd w liczeniu...
trzeba umieć obsługiwać FIND_ORB, chodzi o to że najpierw trzeba podać stałe i wyliczać orbitę wykluczając
stopniowo obserwacje.
Czyli najpierw sobie ustalam
e=1,q=0.005,i=144
a potem gdy orbita jest w miarę podobna do rodziny Kreutza mogę pozbyć się jakiejś stałej zostawiając jedną.
I tak uzyskałem taką orbitę, dość dobrze pasującą do grafiki przedstawionej przez Ciebie Michale
SOHO2143
Perihelion 2011 Oct 1.849601 TT; Constraint: e=1,
Epoch 2011 Sep 30.0 TT = JDT 2455834.5 Earth MOID: 0.5748
M 359.99998 (2000.0) P Q
n 0.00001053 Peri. 97.83164 -0.03228073 -0.99665889
a2061.33306 Node 7.35972 -0.97746009 0.01581077
e 0.9999987 Incl. 144.14785 0.20863776 -0.08013160
P 93588 q 386824.291km Q 4122.66354
From 4 observations 2011 Sept. 30-Oct. 1 (20.2 hr); RMS error 11.738 arcseconds
gdybym zebrał więcej astrometrii mogłoby być dużo dokładniej..
-
Michale mógłbyś pomóc w pewnej sprawie. Chodzi o datę zdjęć na LASCO
np. te zdjęcie
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 /
-
Witam Wszystkich serdecznie,
Ostatnio próbuję obliczać orbity komet SOHO. Na razie z mizernym skutkiem. Mam kilka pytań:
Co oznacza liczba "2" przed wartościami wektorów XYZ sondy ?
Czy są inne programy do liczenia orbit oprócz Find Orb ?
Jak obliczyć orbitę ręcznie mając 3 pozycje? Jakich użyć wzorów ?
-
Kolejne pytanie: co oznacza P, Q i trzy wartości wartości pod nimi. Czasami w ostatnim wersie podaje się jeszcze takie parametry:
P 3.80 H 19.7 G 0.15 U 9. Byłbym wdzięczny jakby ktoś mi to wyjaśnił.
-
Witam Wszystkich serdecznie,
Ostatnio próbuję obliczać orbity komet SOHO. Na razie z mizernym skutkiem. Mam kilka pytań:
Co oznacza liczba "2" przed wartościami wektorów XYZ sondy ?
Czy są inne programy do liczenia orbit oprócz Find Orb ?
Jak obliczyć orbitę ręcznie mając 3 pozycje? Jakich użyć wzorów ?
Zerknij na http://www.minorplanetcenter.net/iau/info/SatelliteObs.html
2 przed XYZ oznacza paralaksę typu 2 czyli odległości są w AU, paralaksa typu 1 ma odległości w km
---
Wiem że SOLEX 11.0 Light ma program ExOrb do liczenia orbit
http://chemistry.unina.it/~alvitagl/solex/
-
Kolejne pytanie: co oznacza P, Q i trzy wartości wartości pod nimi. Czasami w ostatnim wersie podaje się jeszcze takie parametry:
P 3.80 H 19.7 G 0.15 U 9. Byłbym wdzięczny jakby ktoś mi to wyjaśnił.
to są wektory, xyz dla P i xyz dla Q ale nie wiem co P i Q oznacza
Przykład:
http://www.minorplanetcenter.net/db_search/show_object?object_id=136199
-
to są wektory, xyz dla P i xyz dla Q ale nie wiem co P i Q oznacza
Coś znalazłem:
CHAPTER 10 COMPUTATION OF AN EPHEMERIS
www.astro.uvic.ca/~tatum/celmechs/celm10.pdf
więcej na:
http://www.astro.uvic.ca/~tatum/index.html
-
Dzięki Adam.
-
Michale mógłbyś pomóc w pewnej sprawie. Chodzi o datę zdjęć na LASCO
np. te zdjęcie
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").
-
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..
-
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
-
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. :)
-
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
-
Dzięki Michał. Spróbuję coś policzyć.
-
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.
-
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..
-
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!
-
Dzięki Michał. Bacznie obserwuję jak moja kometa leci ku Słońcu. Jest już jaśniejsza i ma warkocz.
-
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 !
-
Dzięki Kanarkusmaximus. To już dziewiąta kometa na moim koncie ;D
-
Dzięki Kanarkusmaximus. To już dziewiąta kometa na moim koncie ;D
Łał, gratulacje!! :)
-
Na chwilę obecną:
(http://img337.imageshack.us/img337/2599/kometa.gif)
-
Kometa Szymona na C2
(http://img21.imageshack.us/img21/4274/201110302036c2.gif)