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

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

Offline Michal Kusiak

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


http://cometman.tk

Offline Limax7

  • Weteran
  • *****
  • Wiadomości: 1130
  • .: I love comets :.
    • Link do Facebook:
Odp: SOHO-2143
« Odpowiedź #31 dnia: Październik 18, 2011, 12:33 »
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...
Adam Hurcewicz
Uniwersał 150/900
Nikon D3200 + Nikkor AF-S 18-105mm + Sigma APO 70-300mm

Białystok, Polska
http://adamhurcewicz.wordpress.com

Offline Michal Kusiak

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



« Ostatnia zmiana: Październik 18, 2011, 13:43 wysłana przez Michal Kusiak »
Pozdrawiam,
Michał Kusiak


http://cometman.tk

Offline Limax7

  • Weteran
  • *****
  • Wiadomości: 1130
  • .: I love comets :.
    • Link do Facebook:
Odp: SOHO-2143
« Odpowiedź #33 dnia: Październik 18, 2011, 14:55 »
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/
« Ostatnia zmiana: Październik 20, 2011, 12:02 wysłana przez Limax7 »
Adam Hurcewicz
Uniwersał 150/900
Nikon D3200 + Nikkor AF-S 18-105mm + Sigma APO 70-300mm

Białystok, Polska
http://adamhurcewicz.wordpress.com

Polskie Forum Astronautyczne

Odp: SOHO-2143
« Odpowiedź #33 dnia: Październik 18, 2011, 14:55 »

Offline Michal Kusiak

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


http://cometman.tk

Offline Limax7

  • Weteran
  • *****
  • Wiadomości: 1130
  • .: I love comets :.
    • Link do Facebook:
Odp: SOHO-2143
« Odpowiedź #35 dnia: Październik 18, 2011, 19:25 »
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
« Ostatnia zmiana: Październik 18, 2011, 19:37 wysłana przez Limax7 »
Adam Hurcewicz
Uniwersał 150/900
Nikon D3200 + Nikkor AF-S 18-105mm + Sigma APO 70-300mm

Białystok, Polska
http://adamhurcewicz.wordpress.com

Offline Michal Kusiak

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



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.



« Ostatnia zmiana: Październik 18, 2011, 19:58 wysłana przez Michal Kusiak »
Pozdrawiam,
Michał Kusiak


http://cometman.tk

Offline Limax7

  • Weteran
  • *****
  • Wiadomości: 1130
  • .: I love comets :.
    • Link do Facebook:
Odp: SOHO-2143
« Odpowiedź #37 dnia: Październik 18, 2011, 23:02 »
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..
« Ostatnia zmiana: Październik 18, 2011, 23:05 wysłana przez Limax7 »
Adam Hurcewicz
Uniwersał 150/900
Nikon D3200 + Nikkor AF-S 18-105mm + Sigma APO 70-300mm

Białystok, Polska
http://adamhurcewicz.wordpress.com

Offline Limax7

  • Weteran
  • *****
  • Wiadomości: 1130
  • .: I love comets :.
    • Link do Facebook:
Odp: SOHO-2143
« Odpowiedź #38 dnia: Październik 19, 2011, 09:30 »
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 /                                               
Adam Hurcewicz
Uniwersał 150/900
Nikon D3200 + Nikkor AF-S 18-105mm + Sigma APO 70-300mm

Białystok, Polska
http://adamhurcewicz.wordpress.com

Offline Szymon

  • Nowy
  • *
  • Wiadomości: 38
Odp: SOHO-2143
« Odpowiedź #39 dnia: Październik 19, 2011, 17:05 »
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 ?
Pozdrawiam,
Szymon Liwo

Offline Szymon

  • Nowy
  • *
  • Wiadomości: 38
Odp: SOHO-2143
« Odpowiedź #40 dnia: Październik 19, 2011, 18:54 »
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ł.
Pozdrawiam,
Szymon Liwo

Offline Limax7

  • Weteran
  • *****
  • Wiadomości: 1130
  • .: I love comets :.
    • Link do Facebook:
Odp: SOHO-2143
« Odpowiedź #41 dnia: Październik 20, 2011, 08:57 »
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/
Adam Hurcewicz
Uniwersał 150/900
Nikon D3200 + Nikkor AF-S 18-105mm + Sigma APO 70-300mm

Białystok, Polska
http://adamhurcewicz.wordpress.com

Offline Limax7

  • Weteran
  • *****
  • Wiadomości: 1130
  • .: I love comets :.
    • Link do Facebook:
Odp: SOHO-2143
« Odpowiedź #42 dnia: Październik 20, 2011, 09:15 »
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

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

Białystok, Polska
http://adamhurcewicz.wordpress.com

Offline Limax7

  • Weteran
  • *****
  • Wiadomości: 1130
  • .: I love comets :.
    • Link do Facebook:
Odp: SOHO-2143
« Odpowiedź #43 dnia: Październik 20, 2011, 09:30 »
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
« Ostatnia zmiana: Październik 20, 2011, 09:32 wysłana przez Limax7 »
Adam Hurcewicz
Uniwersał 150/900
Nikon D3200 + Nikkor AF-S 18-105mm + Sigma APO 70-300mm

Białystok, Polska
http://adamhurcewicz.wordpress.com

Offline Szymon

  • Nowy
  • *
  • Wiadomości: 38
Odp: SOHO-2143
« Odpowiedź #44 dnia: Październik 20, 2011, 10:22 »
Dzięki Adam.
Pozdrawiam,
Szymon Liwo

Polskie Forum Astronautyczne

Odp: SOHO-2143
« Odpowiedź #44 dnia: Październik 20, 2011, 10:22 »