Rádiové pozorování meteorů

Popisy projektů, stavební postupy, konstrukce zařízení, návrh antén, přijímačů, recenze, testy.
Odpovědět
Uživatelský avatar
ok1in
Příspěvky: 987
Registrován: 27. 06. 2017, 19:23
Bydliště: Slaný
Věk: 63
Kontaktovat uživatele:

Re: Rádiové pozorování meteorů

#76

Příspěvek od ok1in »

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...
detekce.jpg
detekce.jpg (604.03 KiB) Zobrazeno 7132 x
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...
Uživatelský avatar
ok1in
Příspěvky: 987
Registrován: 27. 06. 2017, 19:23
Bydliště: Slaný
Věk: 63
Kontaktovat uživatele:

Re: Rádiové pozorování meteorů

#77

Příspěvek od ok1in »

Ještě sem dám jeden obrázek - dají se totiž pozorovat i jiné zajímavosti...
cojeste.png
cojeste.png (787.32 KiB) Zobrazeno 7135 x
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.
ladeni.jpg
ladeni.jpg (82.3 KiB) Zobrazeno 7125 x
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)... ;)
Uživatelský avatar
ok1in
Příspěvky: 987
Registrován: 27. 06. 2017, 19:23
Bydliště: Slaný
Věk: 63
Kontaktovat uživatele:

Re: Rádiové pozorování meteorů

#78

Příspěvek od ok1in »

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.
Uživatelský avatar
ok1in
Příspěvky: 987
Registrován: 27. 06. 2017, 19:23
Bydliště: Slaný
Věk: 63
Kontaktovat uživatele:

Re: Rádiové pozorování meteorů

#79

Příspěvek od ok1in »

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... :mrgreen:
Uživatelský avatar
ok1in
Příspěvky: 987
Registrován: 27. 06. 2017, 19:23
Bydliště: Slaný
Věk: 63
Kontaktovat uživatele:

Re: Rádiové pozorování meteorů

#80

Příspěvek od ok1in »

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
.
event20201207_075345_94.jpg
event20201207_075345_94.jpg (194.77 KiB) Zobrazeno 6924 x
Uživatelský avatar
MaG
Příspěvky: 11380
Registrován: 06. 04. 2002, 21:22
Bydliště: Jablonec nad Nisou
Kontaktovat uživatele:

Re: Rádiové pozorování meteorů

#81

Příspěvek od MaG »

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
Uživatelský avatar
ok1in
Příspěvky: 987
Registrován: 27. 06. 2017, 19:23
Bydliště: Slaný
Věk: 63
Kontaktovat uživatele:

Re: Rádiové pozorování meteorů

#82

Příspěvek od ok1in »

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... ;)
Uživatelský avatar
PetrK
Příspěvky: 509
Registrován: 24. 09. 2014, 21:32
Bydliště: Brno
Věk: 48

Re: Rádiové pozorování meteorů

#83

Příspěvek od PetrK »

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)
meteor-radar_GRAVES_2020-12-07_22-22.png
meteor-radar_GRAVES_2020-12-07_22-22.png (195.73 KiB) Zobrazeno 6867 x
[/url]
Dobson 200/1200, SCT 127/1250, AC 80/400
15x70, 10x30, 6x30
Uživatelský avatar
ok1in
Příspěvky: 987
Registrován: 27. 06. 2017, 19:23
Bydliště: Slaný
Věk: 63
Kontaktovat uživatele:

Re: Rádiové pozorování meteorů

#84

Příspěvek od ok1in »

Je to odraz ;) Rušení je ta slabá čára dole...
Uživatelský avatar
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ů

#85

Příspěvek od BiF »

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...
Možná jsem to trochu vytrhl z kontextu, ale - jsou tedy ty "levné" RTL-SDR použitelné, a nebo ne?
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
Uživatelský avatar
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ů

#86

Příspěvek od BiF »

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
Uživatelský avatar
ok1in
Příspěvky: 987
Registrován: 27. 06. 2017, 19:23
Bydliště: Slaný
Věk: 63
Kontaktovat uživatele:

Re: Rádiové pozorování meteorů

#87

Příspěvek od ok1in »

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ě... ;)
Uživatelský avatar
dablik
Příspěvky: 610
Registrován: 17. 09. 2012, 13:17
Bydliště: Zichovec
Věk: 45

Re: Rádiové pozorování meteorů

#88

Příspěvek od dablik »

Tak mne Petr motivoval a zase jsem pripojil yagi-nu :D 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.
na foceni: SV80 APO, WO102 GT, ASI1600MM cool, 8x50 hledacek + QHY5, hodne ruznych krouzku a udelatek :)  
montaz: EQ6R
Uživatelský avatar
ok1in
Příspěvky: 987
Registrován: 27. 06. 2017, 19:23
Bydliště: Slaný
Věk: 63
Kontaktovat uživatele:

Re: Rádiové pozorování meteorů

#89

Příspěvek od ok1in »

squelch musí být vypnutý... :idea:
Uživatelský avatar
PetrK
Příspěvky: 509
Registrován: 24. 09. 2014, 21:32
Bydliště: Brno
Věk: 48

Re: Rádiové pozorování meteorů

#90

Příspěvek od PetrK »

Dvě zrníčka a pěkný kamínek s dopplerem, 8.12.2020 UTC 10:09 ± 2min.

meteor-radar_GRAVES_2020-12-08_11-09.jpg
meteor-radar_GRAVES_2020-12-08_11-09.jpg (261.54 KiB) Zobrazeno 6761 x
[/url]
Dobson 200/1200, SCT 127/1250, AC 80/400
15x70, 10x30, 6x30
Odpovědět