Rádiové pozorování meteorů
- ok1in
- Příspěvky: 1000
- Registrován: 27. 06. 2017, 19:23
- Bydliště: Slaný
- Věk: 63
- Kontaktovat uživatele:
Re: Rádiové pozorování meteorů
Osobně používám Spectrum Lab, protože jde o velmi povedený program s množstvím možností (nejen pro meteory) a velmi mocným skriptovacím jazykem. Ten totiž v mém případě plně automatizuje detekci a ukládá mi požadovaná data. Vše se skrývá v možnostech pro periodické, časové nebo podmíněné akce programu. Toto je jeho velmi silná stránka. Obrázky z detekcí, které jsem tady nahrával, jsou generované právě skriptem, který ke každé detekci ukládá i audio záznam, který se může kdykoliv později znovu použít a zpracovávat. Takhle vypadají některé detekce...
Co je na tom obrázku?
V horní části je spektrum signálu s dvěmi křivkami - zelená je realtime stopa, červená maximum průběhu události.
Spodní část je tzv. "vodopád", kde je vidět průběh spektra dané události v čase. U každé detekce je na jejím konci časová značka v UTC, velmi snadno se pak párují eventy se zvukovým záznamem, který ji má v názvu souboru (stejně jako obrázek). Stupnice vpravo dole symbolizuje intenzitu odrazu signálu - čímje stopa jasnější, tím silnější je ionizace a také následný odraz signálu.
Vlevo je detekce série krátkých ech, jde o malé meteority, které šly hezky krátce po sobě. Uprostřed dva větší kusy, první s výrazným počátečním dopplerovým posuvem odraženého signálu, druhý relativně dlouhý pomalý rozpad meteoru. Vpravo je ukázka pořádného šutru, dobře znatelný i když menší doppler při začátku a konce detekce a velmi zřetelným rozpadem (výrazná změna intenzity a kmitočtu doppleru v signálu odrazu).
Skript, který vše ovládá sem později nalinkuju, je delší a do zprávy jej nelze uložit...
Co je na tom obrázku?
V horní části je spektrum signálu s dvěmi křivkami - zelená je realtime stopa, červená maximum průběhu události.
Spodní část je tzv. "vodopád", kde je vidět průběh spektra dané události v čase. U každé detekce je na jejím konci časová značka v UTC, velmi snadno se pak párují eventy se zvukovým záznamem, který ji má v názvu souboru (stejně jako obrázek). Stupnice vpravo dole symbolizuje intenzitu odrazu signálu - čímje stopa jasnější, tím silnější je ionizace a také následný odraz signálu.
Vlevo je detekce série krátkých ech, jde o malé meteority, které šly hezky krátce po sobě. Uprostřed dva větší kusy, první s výrazným počátečním dopplerovým posuvem odraženého signálu, druhý relativně dlouhý pomalý rozpad meteoru. Vpravo je ukázka pořádného šutru, dobře znatelný i když menší doppler při začátku a konce detekce a velmi zřetelným rozpadem (výrazná změna intenzity a kmitočtu doppleru v signálu odrazu).
Skript, který vše ovládá sem později nalinkuju, je delší a do zprávy jej nelze uložit...
Petr, www.ok1in.cz
- ok1in
- Příspěvky: 1000
- Registrován: 27. 06. 2017, 19:23
- Bydliště: Slaný
- Věk: 63
- Kontaktovat uživatele:
Re: Rádiové pozorování meteorů
Ještě sem dám jeden obrázek - dají se totiž pozorovat i jiné zajímavosti...
Na obrázku vlevo je silný odraz začínající záporným dopplerem. Ale je tam ještě jedna zajímavost - ty slabé šikmé čáry pod ním jsou totiž odrazy od satelitů. Podle formace to vypadá na vláček jednoho vesmírného bordeláře (neověřoval jsem co v tu dobu bylo v místě odrazu). Pravý obrázek pak ukazuje dvě slabé čáry s odrazem od letadel.
A jedna veledůležitá poznámka na závěr - přijímač musí být naladěný mimo kmitočet vysílače, podle toho, které postranní pásmo u SSB signálu přijímáte! Při USB musí být níže, při LSB výše. Já jsem mimo o 900 Hz od nominálu, zázněj bez doppleru je tedy v audio spektru přesně 900Hz.
Vše ostatní je posunuto. Proto jsou ty stopy od letadel přesně na 900Hz, ale všechny meteory posunuty o doppler. Pokud by jste byli naladěni přesně na kmitočet vysílače (tady 143 050 MHz) bude zázněj nulový, nebo velmi malý a program nebude nic detekovat.
Tož tolik, protože jsem typus technikus, tak dotazy na to co jsem blbě vysvětlil, se budu snažit nějak lidsky zodpovědět (pokud tedy budu vědět)...
Na obrázku vlevo je silný odraz začínající záporným dopplerem. Ale je tam ještě jedna zajímavost - ty slabé šikmé čáry pod ním jsou totiž odrazy od satelitů. Podle formace to vypadá na vláček jednoho vesmírného bordeláře (neověřoval jsem co v tu dobu bylo v místě odrazu). Pravý obrázek pak ukazuje dvě slabé čáry s odrazem od letadel.
A jedna veledůležitá poznámka na závěr - přijímač musí být naladěný mimo kmitočet vysílače, podle toho, které postranní pásmo u SSB signálu přijímáte! Při USB musí být níže, při LSB výše. Já jsem mimo o 900 Hz od nominálu, zázněj bez doppleru je tedy v audio spektru přesně 900Hz.
Vše ostatní je posunuto. Proto jsou ty stopy od letadel přesně na 900Hz, ale všechny meteory posunuty o doppler. Pokud by jste byli naladěni přesně na kmitočet vysílače (tady 143 050 MHz) bude zázněj nulový, nebo velmi malý a program nebude nic detekovat.
Tož tolik, protože jsem typus technikus, tak dotazy na to co jsem blbě vysvětlil, se budu snažit nějak lidsky zodpovědět (pokud tedy budu vědět)...
Petr, www.ok1in.cz
- ok1in
- Příspěvky: 1000
- Registrován: 27. 06. 2017, 19:23
- Bydliště: Slaný
- Věk: 63
- Kontaktovat uživatele:
Re: Rádiové pozorování meteorů
Skript pro detekci
Kód: Vybrat vše
; Exported "Conditional Actions" for Spectrum Lab
if( never ) then # ====== User definable options ======
if( never ) then # SNR_Threshold - Minimum SNR in dB to trigger an event, increase if noise is being interpreted as meteors, reduce if meteors are not being logged
if( never ) then # Meteor_Gap - the period that the signal has to drop below the threshold to be considered the end of a meteor event (in seconds)
if( never ) then # Log_Threshold - the minimum duration of a meteor event that will be logged (in seconds)
if( never ) then # Waterfall_Length - the width or height of the waterfall display in seconds, subtract 20% to be sure of screenshot containing start and end of events
if( never ) then # Log_Path - the location to log files for meteor events (remember to double escape backslashes, including trailing backslash)
if( never ) then # Capture_Path - the location to save screen captures of meteor events (format as above)
if( never ) then # Sound_Path - the location to save .WAV sound capture of meteor events (format as above)
if( never ) then # Sound_Capture - set to 1 to capture meteor sounds, 0 to disable sound capture
if( initialising ) then SNR_Threshold = 16 : Meteor_Gap = 2 : Log_Threshold = 0.05 : Waterfall_Length = 20 : Log_Path = "c:\\Astro\\Radio\\Spectrum\\_data\\Captured_LOG\\" : Capture_Path = "c:\\Astro\\Radio\\Spectrum\\_data\\Captured_PIC\\" : Sound_Path = "c:\\Astro\\Radio\\Spectrum\\_data\\Captured_SND\\" : Sound_Capture = 1
if( never ) then # ====== End Of User definable options ======
if( never ) then ##### Initialise Variables #######
if( initialising ) then Current_Noise = noise(cfg.SpecFreqMin, cfg.SpecFreqMax) * 6 : Current_Signal = peak_a(cfg.SpecFreqMin, cfg.SpecFreqMax): Current_Signal = 0 : Current_SNR = 0 : Current_Peak_Frequency = 0 : Current_Time = 0
if( initialising ) then Meteor_Start_Time = 0 : Meteor_End_Time = 0 : Daily_Meteor_Count = 0 : Midnight_Flag = 0 : Hourly_Meteor_Count = 0 : Hourly_Meteor_Duration = 0 : Hourly_Meteor_Longest = 0 : Hourly_Flag = 0 : Logging = 1 : Capture_Time = 0
if( initialising ) then Meteor_Max_Signal = 0 : Meteor_Max_Noise = 0 : Meteor_Max_SNR = 0: Meteor_Max_Peak = 0 : Meteor_Count = 0
if( initialising ) then Intermediate_Capture = Waterfall_Length : State = "Waiting" : Queued = 0
if( never ) then ##### Measure current signal, exponential weighted average of noise and peak frequency after each FFT calculation #######
if( new_spectrum ) then Current_Noise = (Current_Noise / 6) * 5 + noise(cfg.SpecFreqMin, cfg.SpecFreqMax) : Current_Signal = peak_a(cfg.SpecFreqMin, cfg.SpecFreqMax) : Current_Peak_Frequency = peak_f(cfg.SpecFreqMin, cfg.SpecFreqMax) : Current_SNR = Current_Signal - (Current_Noise / 6) : Current_Time = time
if( never ) then ##### Signal above threshold, no meteor in progress so start of new meteor event ######
if( State="Waiting" & Current_SNR>=SNR_Threshold ) then State = "Meteor" : Meteor_Start_Time = Current_Time : rec.filename = Sound_Path + "event" + str("YYYYMMDD", now) + "_" + str("hhmmss", now) + "_" + str(Daily_Meteor_Count + 1) + ".wav"
if( never ) then ##### SNR is higher than previous max SNR, record new max (subtraction gives SNR in dB) #####
if( State="Meteor" & Current_SNR>=Meteor_Max_SNR ) then Meteor_Max_Signal = Current_Signal : Meteor_Max_Noise = (Current_Noise / 6) : Meteor_Max_Peak = Current_Peak_Frequency : Meteor_Max_SNR = Current_SNR : rec.trigger = Sound_Capture
if( never ) then ##### Meteor in progress, signal below threshold so start timing gap in case of short break in meteor signal ######
if( State="Meteor" & Current_SNR<SNR_Threshold ) then State = "Gap" : Meteor_End_Time = Current_Time
if( never ) then ##### Timing gap, signal back above threshold so resume current meteor event #######
if( State="Gap" & Current_SNR>=SNR_Threshold ) then State = "Meteor"
if( never ) then ##### Meteor hasn't finished, but start or middle of event needs to be captured before it scrolls off ######
if( (State="Gap" | State="Meteor") & (Current_Time - Meteor_Start_Time)>=Intermediate_Capture & Logging ) then Intermediate_Capture = Intermediate_Capture + Waterfall_Length : capture(Capture_Path + "event" + str("YYYYMMDD", now) + "_" + str("hhmmss", now) + "_" + str(Daily_Meteor_Count + 1) + ".jpg", 100)
if( never ) then ##### Timing gap, signal below threshold and gap duration exceeded so log meteor event ######
if( State="Gap" & Current_SNR<SNR_Threshold & (time - Meteor_End_Time)>=Meteor_Gap ) then State = "Log" : rec.trigger = 0
if( never ) then ##### Log last meteor event ######
if( State="Log" & (Meteor_End_Time - Meteor_Start_Time)>=Log_Threshold & Logging ) then Daily_Meteor_Count = Daily_Meteor_Count + 1 : Hourly_Meteor_Count = Hourly_Meteor_Count + 1 : Hourly_Meteor_Duration = Hourly_Meteor_Duration + (Meteor_End_Time - Meteor_Start_Time)
if( continuation ) then sp.print("--- Det:", str(Daily_Meteor_Count),"=", str(Meteor_End_Time - Meteor_Start_Time) + "s,", str("hh : mm : ss", Meteor_Start_Time))
if( continuation ) then Queued = Queued + 1 : queue_event(now + max(1, (Waterfall_Length - (Meteor_End_Time - Meteor_Start_Time))), 0)
if( continuation ) then fopen(Log_Path + "event_log_" + str("YYYYMM", Meteor_Start_Time) + ".csv",a,r)
if( continuation ) then fp(str("YYYY/MM/DD",Meteor_Start_Time) + "," + str("hh:mm:ss",Meteor_Start_Time) + "," + str(Daily_Meteor_Count) + "," + str(Meteor_Max_Signal) + "," + str(Meteor_Max_Noise) + "," + str(Meteor_Max_SNR) + "," + str(Meteor_Max_Peak) + "," + str(Meteor_End_Time - Meteor_Start_Time))
if( continuation ) then fclose
if( never ) then ##### Reset and wait for next meteor event ########
if( State="Log" & Hourly_Meteor_Longest<(Meteor_End_Time - Meteor_Start_Time) ) then Hourly_Meteor_Longest = (Meteor_End_Time - Meteor_Start_Time)
if( State="Log" ) then Intermediate_Capture = Waterfall_Length : Meteor_Max_Signal = 0 : Meteor_Max_Noise = 0 : Meteor_Max_SNR = 0: Meteor_Max_Peak = 0
if( continuation ) then State = "Waiting"
if( never ) then ##### Capture a screenshot of the display at end of event (see also Screen Capture tab) ########
if( queued_event ) then capture(Capture_Path + "event" + str("YYYYMMDD", now) + "_" + str("hhmmss", now) + "_" + str(Daily_Meteor_Count - Queued + 1) + ".jpg", 100)
if( continuation ) then Queued = Queued - 1
if( never ) then ##### Increment the hourly flag (this will happen multiple times during the last five seconds of the hour) ######
if( str("mmss", Current_Time)="5955" ) then Hourly_Flag = Hourly_Flag + 1
if( never ) then # If the hourly flag has just been triggered, log and then reset the meteor count to zero for the new hour, if the hourly flag has been triggered already, do nothing
if( Hourly_Flag=1 ) then fopen2(Log_Path + "hourly_log_" + str("YYYYMM", Current_Time) + ".csv",a,r)
if( continuation ) then fp2(str("YYYY/MM/DD", Current_Time) + "," + str("hh", Current_Time) + "," + str(Hourly_Meteor_Count))
if( continuation ) then fclose2
if( continuation ) then fopen4(Log_Path + "RMOB-" + str("YYYYMM", Current_Time) + ".dat",a,r)
if( continuation ) then fp4(str("YYYYMMDDhh", Current_Time) + "," + str("hh", Current_Time) + "," + str(Hourly_Meteor_Count))
if( continuation ) then fclose4
if( continuation ) then fopen5(Log_Path + "RMOB_Dur-" + str("YYYYMM", Current_Time) + ".dat",a,r)
if( continuation ) then fp5(str("YYYYMMDDhh", Current_Time) + "," + str("hh", Current_Time) + "," + str(Hourly_Meteor_Count) + "," + str(Hourly_Meteor_Duration) + "," + str(Hourly_Meteor_Longest) + "," + str(Current_Noise))
if( continuation ) then fclose5
if( continuation ) then Hourly_Meteor_Count = 0 : Hourly_Meteor_Duration = 0 : Hourly_Meteor_Longest = 0
if( never ) then # If the current time is on the hour, reset the hourly flag ready to trigger again at the end of the next hour
if( str("mmss", Current_Time)="0000" ) then Hourly_Flag = 0
if( never ) then ##### Increment the midnight flag (this will happen multiple times during the last fice seconds of the day) ######
if( str("hhmmss", Current_Time)="235955" ) then Midnight_Flag = Midnight_Flag + 1
if( never ) then # If the midnight flag has just been triggered, log and then reset the meteor count to zero for the new day, if the midnight flag has been triggered already, do nothing
if( Midnight_Flag=1 ) then fopen3(Log_Path + "daily_log_" + str("YYYYMM", Current_Time) + ".csv",a,r)
if( continuation ) then fp3(str("YYYY/MM/DD", Current_Time) + "," + str(Daily_Meteor_Count))
if( continuation ) then fclose3
if( continuation ) then Daily_Meteor_Count = 0
if( never ) then # If the current time is midnight, reset the midnight flag ready to trigger again at the end of the next day
if( str("hhmmss", Current_Time)="000000" ) then Midnight_Flag = 0
if( never ) then ##### ######
if( never ) then # TODO: upload to web site? auto image display online, auto generate colourgramme via web.
Petr, www.ok1in.cz
- ok1in
- Příspěvky: 1000
- Registrován: 27. 06. 2017, 19:23
- Bydliště: Slaný
- Věk: 63
- Kontaktovat uživatele:
Re: Rádiové pozorování meteorů
Nastavení programu SpecLab
Důležité je nastavit v konfiguraci programu a skriptu správné cesty pro ukládání souborů a zejména pak hodnotu "SNR_Threshold". Její hodnota musí být o ca 1-2 db vyšší, než kolik zobrazuje program na liště vlevo - SNR při live detekci. Čím vyšší bude rozdíl, tím horší bude detekce slabých odrazů, těsně nad šumem.
Tak a to je snad vše, doufám, že jsem vám to neznechutil...
Důležité je nastavit v konfiguraci programu a skriptu správné cesty pro ukládání souborů a zejména pak hodnotu "SNR_Threshold". Její hodnota musí být o ca 1-2 db vyšší, než kolik zobrazuje program na liště vlevo - SNR při live detekci. Čím vyšší bude rozdíl, tím horší bude detekce slabých odrazů, těsně nad šumem.
Tak a to je snad vše, doufám, že jsem vám to neznechutil...
Petr, www.ok1in.cz
- ok1in
- Příspěvky: 1000
- Registrován: 27. 06. 2017, 19:23
- Bydliště: Slaný
- Věk: 63
- Kontaktovat uživatele:
Re: Rádiové pozorování meteorů
Jeden dnešní kousek, stojící za pozornost, celková délka jevu ca 27 sec. - byly tam ještě dozněje slabých odrazů po rozpadu.
Vložené údaje ve spektru na konci odrazu:
---pořadí detekce ve dni
= délka úkazu + čas v UTC
.
Vložené údaje ve spektru na konci odrazu:
---pořadí detekce ve dni
= délka úkazu + čas v UTC
.
Petr, www.ok1in.cz
- MaG
- Příspěvky: 11419
- Registrován: 06. 04. 2002, 21:22
- Bydliště: Jablonec nad Nisou
- Kontaktovat uživatele:
Re: Rádiové pozorování meteorů
Díky moc za skvělý popis. Přál bych si, aby to navnadilo nové fandy do radioastronomie i rádia obecně. Rád bych to nějak zpopularizoval i na astro.cz, ale na tohle nemám, leda převzít zdejší text, pokud bys souhlasil.
Martin Gembec, Astronomické události
EQ6 OnStep, WO FLT98+WO AFR-IV 0.8x, Orion Optics CT8 200/900 s Paracorrem, ZWO ASI294MC Pro, Canon 6Dmod, 30Dmod
EQ6 OnStep, WO FLT98+WO AFR-IV 0.8x, Orion Optics CT8 200/900 s Paracorrem, ZWO ASI294MC Pro, Canon 6Dmod, 30Dmod
- ok1in
- Příspěvky: 1000
- Registrován: 27. 06. 2017, 19:23
- Bydliště: Slaný
- Věk: 63
- Kontaktovat uživatele:
Re: Rádiové pozorování meteorů
Martine, nemám s tím žádný problém - jen by to chtělo, kdyby to raději ještě nekdo zkouknul a zkritizoval. Je to jen hrubý náhled do oblasti, kde je mnohem více možností, popsal jsem jen to nejjednodušší, co může podle mě každý zvládnout sám. Opravdu nejsem žádný renomovaný odborník pro tuto oblast, jen se tím už mnoho let amatérsky zabývám. To jen aby jsi u nějakých opravdových odborníků nenarazil, že publikuješ bludy...
Petr, www.ok1in.cz
Re: Rádiové pozorování meteorů
Snad mám taky jedno malé zrníčko. Pořízeno v 22 hodin a několik minut. Může mi to prosím někdo okomentovat, jestli je to meteor, nebo rušení a jak vylepšit nastavení?
(na kliknutí větší náhled)
[/url]
(na kliknutí větší náhled)
[/url]
Dobson 200/1200, SCT 127/1250, AC 80/400
15x70, 10x30, 6x30
15x70, 10x30, 6x30
- BiF
- Příspěvky: 1143
- Registrován: 10. 11. 2002, 02:00
- Bydliště: Trhové Sviny
- Věk: 52
- Kontaktovat uživatele:
Re: Rádiové pozorování meteorů
Možná jsem to trochu vytrhl z kontextu, ale - jsou tedy ty "levné" RTL-SDR použitelné, a nebo ne?ok1in píše: Druhá možnost je použít vlastní přijímač typu RTL-SDR (velmi levné s technickými omezeními)...
... a prakticky se nedá k ničemu použít...
A jaký případně vybrat? Na e-bay se dá usb RTL-SDR pořídit od 400,- Kč, u nás v rozsahu 500-2000 Kč. Jaký je v nich rozdíl (kromě rozsahu frekvencí)?
Bude platit to, že ten dražší bude lepší?
Newtony 300/1200, 200/860, 100/600; SkyMaster 15x70; paralaktická montáž vlastní výroby; QHY163M; ASI120M; Canon EOS 300D; Philips ToUcam Pro; http://bilek.astronomie.cz/ , email: frantabilek(at)gmail(dot)com
- BiF
- Příspěvky: 1143
- Registrován: 10. 11. 2002, 02:00
- Bydliště: Trhové Sviny
- Věk: 52
- Kontaktovat uživatele:
Re: Rádiové pozorování meteorů
např. toto: 100KHz-1.7GHz Full Band UV HF RTL-SDR USB Tuner Receiver/ R820T+8232 Ham Radio
nebo toto (12-Bit; není to ale RTL-SCR ale nějaký SDR-RE-Play)?: Wideband Full Featured 12bit SDR Receiver SDRPLAY RSP1 RSP2 RTL-SDR HackRF Upgrade AM FM HF SSB CW receiver Full band HAM Radio
nebo ještě něco jiného???
nebo toto (12-Bit; není to ale RTL-SCR ale nějaký SDR-RE-Play)?: Wideband Full Featured 12bit SDR Receiver SDRPLAY RSP1 RSP2 RTL-SDR HackRF Upgrade AM FM HF SSB CW receiver Full band HAM Radio
nebo ještě něco jiného???
Newtony 300/1200, 200/860, 100/600; SkyMaster 15x70; paralaktická montáž vlastní výroby; QHY163M; ASI120M; Canon EOS 300D; Philips ToUcam Pro; http://bilek.astronomie.cz/ , email: frantabilek(at)gmail(dot)com
- ok1in
- Příspěvky: 1000
- Registrován: 27. 06. 2017, 19:23
- Bydliště: Slaný
- Věk: 63
- Kontaktovat uživatele:
Re: Rádiové pozorování meteorů
Večer k tomu něco připíšu, jsem na servisu a tak přes den není moc času reagovat. Jen rychlé info, jednodenní test klasického RTL-SDR za 300 dopadl dost špatně...
Petr, www.ok1in.cz
Re: Rádiové pozorování meteorů
Tak mne Petr motivoval a zase jsem pripojil yagi-nu To je proste hruza. Neco sem nahodi a ja to musim hned zkusit.
Nejspis ji budu muset otocit spravnym smerem a pohrat si s detekci. Zatim slysim jenom hezky sum Nebo taky neco udelat se squelch.
Nejspis ji budu muset otocit spravnym smerem a pohrat si s detekci. Zatim slysim jenom hezky sum Nebo taky neco udelat se squelch.
na foceni: SV80 APO, WO102 GT, ASI1600MM cool, 8x50 hledacek + QHY5, hodne ruznych krouzku a udelatek
montaz: EQ6R
montaz: EQ6R
Re: Rádiové pozorování meteorů
Dobson 200/1200, SCT 127/1250, AC 80/400
15x70, 10x30, 6x30
15x70, 10x30, 6x30