Autor Wątek: SOHO-2143  (Przeczytany 8909 razy)

0 użytkowników i 1 Gość przegląda ten wątek.

Offline Szaniu

  • Weteran
  • *****
  • Wiadomości: 1543
  • Szepty Ziemi.
    • ExoplanetSpot
Odp: SOHO-2143
« Odpowiedź #15 dnia: Październik 06, 2011, 15:28 »
Da się jakkolwiek oszacować masę i rozmiary tego obiektu?

Offline Michal Kusiak

  • Junior
  • **
  • Wiadomości: 69
  • podczas gry w kometkę lubię żuć gumę orbit
Odp: SOHO-2143
« Odpowiedź #16 dnia: Październik 16, 2011, 20:20 »
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.

« Ostatnia zmiana: Październik 16, 2011, 20:31 wysłana przez Michal Kusiak »
Pozdrawiam,
Michał Kusiak


http://cometman.tk

Offline Michal Kusiak

  • Junior
  • **
  • Wiadomości: 69
  • podczas gry w kometkę lubię żuć gumę orbit
Odp: SOHO-2143
« Odpowiedź #17 dnia: Październik 16, 2011, 20:39 »
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.
Pozdrawiam,
Michał Kusiak


http://cometman.tk

Offline kanarkusmaximus

  • Administrator
  • *****
  • Wiadomości: 15584
  • Ja z tym nie mam nic wspólnego!
    • Kosmonauta.net
Odp: SOHO-2143
« Odpowiedź #18 dnia: Październik 16, 2011, 21:03 »
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ą. ;)

Offline Limax7

  • Weteran
  • *****
  • Wiadomości: 539
  • .: I love comets :.
    • adamhurcewicz.wordpress.com
Odp: SOHO-2143
« Odpowiedź #19 dnia: Październik 17, 2011, 12:21 »
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
Adam Hurcewicz
Uniwersał 150/900
aparat cyfrowy Nikon D7100 + Nikkor AF-S 18-105mm + Sigma APO 70-300mm

Białystok +53°07'09.68" N, +23°08'15.25" E, Alt=164m
http://adamhurcewicz.wordpress.com

Offline Limax7

  • Weteran
  • *****
  • Wiadomości: 539
  • .: I love comets :.
    • adamhurcewicz.wordpress.com
Odp: SOHO-2143
« Odpowiedź #20 dnia: Październik 17, 2011, 12:26 »
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 ?
« Ostatnia zmiana: Październik 21, 2011, 10:42 wysłana przez Limax7 »
Adam Hurcewicz
Uniwersał 150/900
aparat cyfrowy Nikon D7100 + Nikkor AF-S 18-105mm + Sigma APO 70-300mm

Białystok +53°07'09.68" N, +23°08'15.25" E, Alt=164m
http://adamhurcewicz.wordpress.com

Offline Michal Kusiak

  • Junior
  • **
  • Wiadomości: 69
  • podczas gry w kometkę lubię żuć gumę orbit
Odp: SOHO-2143
« Odpowiedź #21 dnia: Październik 17, 2011, 12:55 »
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.
« Ostatnia zmiana: Październik 17, 2011, 13:01 wysłana przez Michal Kusiak »
Pozdrawiam,
Michał Kusiak


http://cometman.tk

Offline Limax7

  • Weteran
  • *****
  • Wiadomości: 539
  • .: I love comets :.
    • adamhurcewicz.wordpress.com
Odp: SOHO-2143
« Odpowiedź #22 dnia: Październik 17, 2011, 15:25 »
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ąć.
Adam Hurcewicz
Uniwersał 150/900
aparat cyfrowy Nikon D7100 + Nikkor AF-S 18-105mm + Sigma APO 70-300mm

Białystok +53°07'09.68" N, +23°08'15.25" E, Alt=164m
http://adamhurcewicz.wordpress.com

Offline Michal Kusiak

  • Junior
  • **
  • Wiadomości: 69
  • podczas gry w kometkę lubię żuć gumę orbit
Odp: SOHO-2143
« Odpowiedź #23 dnia: Październik 17, 2011, 17:11 »
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/   
« Ostatnia zmiana: Październik 17, 2011, 17:17 wysłana przez Michal Kusiak »
Pozdrawiam,
Michał Kusiak


http://cometman.tk

Offline Limax7

  • Weteran
  • *****
  • Wiadomości: 539
  • .: I love comets :.
    • adamhurcewicz.wordpress.com
Odp: SOHO-2143
« Odpowiedź #24 dnia: Październik 17, 2011, 19:29 »
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źć ...
Adam Hurcewicz
Uniwersał 150/900
aparat cyfrowy Nikon D7100 + Nikkor AF-S 18-105mm + Sigma APO 70-300mm

Białystok +53°07'09.68" N, +23°08'15.25" E, Alt=164m
http://adamhurcewicz.wordpress.com

Offline Michal Kusiak

  • Junior
  • **
  • Wiadomości: 69
  • podczas gry w kometkę lubię żuć gumę orbit
Odp: SOHO-2143
« Odpowiedź #25 dnia: Październik 17, 2011, 19:36 »
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

« Ostatnia zmiana: Październik 17, 2011, 20:33 wysłana przez Michal Kusiak »
Pozdrawiam,
Michał Kusiak


http://cometman.tk

Offline Limax7

  • Weteran
  • *****
  • Wiadomości: 539
  • .: I love comets :.
    • adamhurcewicz.wordpress.com
Odp: SOHO-2143
« Odpowiedź #26 dnia: Październik 17, 2011, 21:25 »
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
Adam Hurcewicz
Uniwersał 150/900
aparat cyfrowy Nikon D7100 + Nikkor AF-S 18-105mm + Sigma APO 70-300mm

Białystok +53°07'09.68" N, +23°08'15.25" E, Alt=164m
http://adamhurcewicz.wordpress.com

Offline Michal Kusiak

  • Junior
  • **
  • Wiadomości: 69
  • podczas gry w kometkę lubię żuć gumę orbit
Odp: SOHO-2143
« Odpowiedź #27 dnia: Październik 17, 2011, 22:32 »
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 ;)

« Ostatnia zmiana: Październik 17, 2011, 22:44 wysłana przez Michal Kusiak »
Pozdrawiam,
Michał Kusiak


http://cometman.tk

Offline Limax7

  • Weteran
  • *****
  • Wiadomości: 539
  • .: I love comets :.
    • adamhurcewicz.wordpress.com
Odp: SOHO-2143
« Odpowiedź #28 dnia: Październik 17, 2011, 23:38 »
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
Adam Hurcewicz
Uniwersał 150/900
aparat cyfrowy Nikon D7100 + Nikkor AF-S 18-105mm + Sigma APO 70-300mm

Białystok +53°07'09.68" N, +23°08'15.25" E, Alt=164m
http://adamhurcewicz.wordpress.com

Offline Limax7

  • Weteran
  • *****
  • Wiadomości: 539
  • .: I love comets :.
    • adamhurcewicz.wordpress.com
Odp: SOHO-2143
« Odpowiedź #29 dnia: Październik 18, 2011, 00:54 »
Mały przykład. Poniżej zdjęcie z SOHO 30.09.2011 godz. 08:30
i ten sam widok z SOHO wygenerowany Guidem

Adam Hurcewicz
Uniwersał 150/900
aparat cyfrowy Nikon D7100 + Nikkor AF-S 18-105mm + Sigma APO 70-300mm

Białystok +53°07'09.68" N, +23°08'15.25" E, Alt=164m
http://adamhurcewicz.wordpress.com