Rev Author Line No. Line
3053 kaklik 1 ;Josef Szylar ver 1.4 - 8.2.2013
2 ;8.2.2013 oprava t0_1 na t0_2 => moc rychlé děje
3 ;oprava dur_count1 na dur_count_1
4 ;nově záznam o generování screenshotu bez meteoru
5 ;oprava timer3.expired na timer3.expired(1)
6  
7 ;PODMÍNKA JEN 80 ZNAKŮ
8 ;PŘÍKAZY JEN 255 ZNAKŮ
9  
10 ;Časovače
11 ;timer0 - časovač maximální mezery mezi 2 meteory K_max_time_btw_met
12 ;timer3 - časovač pro odpočet screenshotu bez meteoru K_max_time_to_scr
13 ;timer4 - časovač pro odečet screenshotu s meteorem po odvalení několika málo sekund aby nebyl meteor moc nahoře
14  
15 ;Proměnné:
16 ;n=0 noise background šum pozadí
17 ;n_1=0 šum v okamžiku detekce meteoru
18 ;f=0 frequency of maximum amplitude frekvence na které se nachází maximum
19 ;aver_amp=0 average amplitude +- 100 Hz průměrná aplituda v rozsahu +- 100Hz od maxima
20 ;aver_amp_1=-999 je average amplitude aktualizovaná v průběhu měření meteoru
21 ;aver_amp_max=-999 je maximální hodnota z hodnot v průběhu měření meteoru
22 ;mag=0 magnitude
23 ;hour_count=0 count of meteors in hour
24 ;strenght_count=0 strenght (souborné hodinové počítadlo délky aktivity meteorů počítadlo 50ms) do RMOB_Dur
25 ;dur_count=0 maximální délka meteoru (count po dobu existence meteoru)
26 ;dur_count_1=0 pomocný k uchování hodnoty dur_count
27 ;dur_count_strenght=0 cout trvání meteoru po celou dobu
28 ;t0=0 průběžný čas a datum (aktualizace každých 50ms) v textové formě "YYYYMMDDhh"
29 ;t0_1 uložený čas pro jméno datového souboru k odeslaní ftp po uplynulé hodině
30 ;t1=0 čas start aktuální meteor v UNIX timestamp formě
31 ;t1_1=0 uložený čas pro zápis do RMOB_dur souboru value:event
32 ;t1_2=0 uložený čas pro odesílání "no meteor" screenshotu
33 ;t2=0 průběžný čas hour (hodnota průběžné hodiny) v textové formě "hh"
34 ;t3=0 průběžný čas hodiny minuty a vteřiny v textové formě "mmss"
35 ;t4=0 čas meteoru s nejdelší dobou aktivity (násobek K_int_cond * dur_count)
36 ;id_met=0 identifikátor meteoru do screenshotu
37  
38 ;Konstanty:
39 ;K_station_name="JOSY_A" station name (devices code - 7 characters maximum, only "A-Z" and "_" - get your unique station code on http://home.robozor.cz)
40 ;K_path_scr="/media/sd/capture/" vložit lokální adresář pro screenshoty
41 ;K_path_audio="/media/sd/capture_audio/" vložit lokální adresář pro audio soubory
42 ;K_path_RMOB="/media/sd/RMOB/" vložit lokální adresář textové soubory RMOB
43 ;K_path_SDR="/media/sd/SDR_data/" vložit lokální adresář textové soubory RMOB
44 ;K_max_time_btw_met=1 maximální doba mezi 2 meteory jsou 1 vteřina (nižší se počítá jako 1 meteor)
45 ;K_max_time_to_scr=100 maximální doba bez meteoru, která pokud vyprší, udělá screenshot (max doba za kterou odroluje waterfall)
46 ;K_init_cond=50 interval pro nastavení vyhodnocení scriptu 50ms Evaluate SpectrumLab condition interval (50ms)
47 ;K_time_to_start_up_waterf=4 vteřiny od okraje screenshot Kolik sekund se má počkat do screenshotu po detekci metreoru (jak daleko od okraje)
48 ;K=7 - citlivost detekce meteoru (odstup od šumu v db)
49 ;K_min_dur_count_bolid=50 detekce bolidu pro audio (2,5 vteřiny)
50  
51 ;Události
52 ;A_init (inicializace)
53 ;A_measurement (kdykoliv po 50 ms)
54 ;A_still_detect (událost trvání meteoru nebo bolidu) - průběžně
55 ;A_start_detect (událost detekce začátku meteoru nebo bolidu)
56 ;A_start_detect_bolid (událost 2,5 vteřiny od detekce meteoru)
57 ;A_end_detect (událost detekce konce meteoru nebo bolidu)
58 ;A_end_detect_end (konec cyklu události detekce konce meteoru nebo bolidu)
59 ;A_write_waterf (událost záznamu popisu detekce do waterfallu)
60 ;A_update_SDR_data_file (událost generování souboru SDR screenshot při záznamu meteoru)
61 ;A_update_RMOB_Dur_file (událost generování souborů RBOB)
62 ;A_met_screenshot (událost generování screenshotu waterfallu pokud je zaznamenán meteor)
63 ;A_no_met_screenshot (událost generování screenshotu waterfallu pokud není dlouho zaznamenán meteor)
64 ;A_once_hour_20s_before (událost 1x za hodinu před ukončením hodiny 20 vteřin)
65 ;A_once_hour
66 ;A_once_hour_10s_before (událost 1x za hodinu před ukončením hodiny 10 vteřin)
67 ;A_once_hour_30s_after (událost 1x za hodinu po ukončení hodiny 30 vteřin)
68  
69  
70 ;-----------------A_init - Inicializace proměnných-------------------------------
71 if( initialising ) then A_init=1:A_still_detect=0:n=0:n_1=0:f=0:f_1=0:aver_amp=0:a=0:mag=0:mag_max=0:hour_count=0:strenght_count=0:strenght_count_1=0:dur_count=0:dur_count_1=0:aver_amp_1=-999:aver_amp_max=-999:t0=0:t0_1=0:t1=0:t1_1=0:t1_2=0:t2=0:t3=0:t4=0
72 if( A_init = 1 ) then A_init=2:id_met="no":id_met2="no":K_station_name="JOSY_A":K_path_scr="capture\\":K_path_audio="capture_audio\\":K_path_RMOB="C:\\Spectrum\\":K_path_SDR="C:\\Spectrum\\":K_path="Z:\\sd\\"
73 if( A_init = 2 ) then A_init=3:K_max_time_btw_met=2:K_max_time_to_scr=100:K_init_cond=50:K_time_to_start_up_waterf=2:K=7:K_min_dur_count_bolid=50
74 if( A_init = 3 ) then A_init=4:A_measurement=0:A_start_detect=0:A_end_detect=0:A_end_detect_end=0:A_write_waterf=0:A_update_SDR_data_file=0:A_update_RMOB_Dur_file=0:A_met_screenshot=0
75 if( A_init = 4 ) then A_init=5:A_no_met_screenshot=0:A_once_hour_20s_before=0:A_once_hour=0:A_once_hour_10s_before=0:A_once_hour_30s_after=0:A_once_hour_40s_after=0:timer3.restart(K_max_time_to_scr)
76 if( A_init = 5 ) then A_init=0:
77  
78 ;----------------A_measurement - (po 50 ms)-------------------------------------
79 ;n=noise(9800,10200) aktualizace a výpočet hodnoty šumu (záleží na nastavení FFT)
80 ;f=peak_f(10300,10900) max intenzita signálu v daném intervalu frekvencí
81 ;aver_amp=avrg(f-100,f+100) vrací prům. intenzitu signálu 100 Hz kolem maxima
82 ;t0=str("YYYYMMDDhh",now) uloží aktuální datum
83 ;t1=now uloží UNIX time stamp na microsekundy (tečka oddělovač)
84 ;t2=str("hh",now) uloží aktuální hodinu
85 ;t3=str("mmss",now) uloží aktuální minutu a vteřinu
86 ;-------------------------------------------------------------------------------
87 if( always ) then A_measurement=1:n=noise(9800,10200):f=peak_f(10300,10900):aver_amp=avrg(f-100,f+100):a=(n+K)
88 if( A_measurement=1 ) then A_measurement=0:t0=str("YYYYMMDDhh",now):t1=round(now*1000):t2=str("hh",now):t3=str("mmss",now)
89  
90 ;-----------------A_still_detect (událost trvání meteoru nebo bolidu)-----------
91 if( aver_amp>a ) then A_still_detect=A_still_detect+1:timer0.restart(K_max_time_btw_met):timer3.restart(K_max_time_to_scr):dur_count=dur_count+1:strenght_count=strenght_count+1:aver_amp_1=aver_amp
92  
93 ;-----------------A_aver_amp (událost aktualizace maximální hodnoty aver_amp)---
94 if( aver_amp_1>aver_amp_max ) then aver_amp_max=aver_amp_1
95  
96 ;-------A_aver_amp (událost aktualizace maximální hodnoty strenght_count_1)-----
97 if( strenght_count>strenght_count_1 ) then strenght_count_1=strenght_count
98  
99 ;-----------------A_start_detect (událost začátku meteoru)----------------------
100 if( A_still_detect=1 ) then A_start_detect=1:t1_1=t1:id_met=(K_station_name+str(t1_1)):sp.print("_ __t+"+str("ss",time))
101 if( A_start_detect=1 ) then A_start_detect=2:n_1=n:f_1=f
102 if( A_start_detect=3 ) then A_start_detect=0:
103  
104 ;------------A_end_detect (událost detekce konce meteoru nebo bolidu)-----------
105 if( timer0.expired(1) ) then A_end_detect=1:id_met2=id_met:hour_count=hour_count+1:timer4.restart(K_time_to_start_up_waterf):dur_count_1=dur_count:
106 if( A_end_detect=1 ) then A_end_detect=2:sp.print("_ __"+id_met2+" HCount"+str(hour_count)+" nb"+str(round(n_1))+" f"+str(f_1)+" mag"+str(round(n_1-aver_amp_max+10)))
107 if( A_end_detect=2 ) then A_end_detect=3:A_update_SDR_data_file=1
108 if( A_end_detect=3 ) then A_end_detect=0
109  
110 ;-----------------A_update_SDR_data_file----------------------------------------
111 ;SDR datový soubor záznamu radiometeorů (projekt astrozor.cz)
112 ;Název souboru:JOSY_A20130123.dat
113 if( A_update_SDR_data_file=1 ) then A_update_SDR_data_file=2:fopen1(K_path_SDR+K_station_name+t0+".dat",a)
114 if( A_update_SDR_data_file=2 ) then A_update_SDR_data_file=3:fp1(id_met2+" ; "+str(round(n_1))+" ; "+str(f_1)+" ; "+str(aver_amp_max)+" ; "+str(round(n_1-aver_amp_max+10))+" ; "+str(dur_count_1)+" ; "+str(K_init_cond)+" ; "+str(K_min_dur_count_bolid))
115 if( A_update_SDR_data_file=3 ) then A_update_SDR_data_file=4:fclose1
116 if( A_update_SDR_data_file=4 ) then A_update_SDR_data_file=5:A_update_RMOB_Dur_file=1
117 if( A_update_SDR_data_file=5 ) then A_update_SDR_data_file=0
118  
119 ;-----A_update_RMOB_Dur_file (událost generování souborů RMOB-Dur)--------------
120 ;RMOB_Dur - soubor pro záznam jednotlivých radiometeorů (projekt www.rmob.org)
121 ;Název souboru:RMOB_Dur-201301.dat,kde 2013 je rok a 01 je měsíc
122 if( A_update_RMOB_Dur_file=1 ) then A_update_RMOB_Dur_file=2:fopen2(K_path_RMOB+"RMOB_Dur-"+str("YYYYMM",now)+".dat",a)
123 if( A_update_RMOB_Dur_file=2 ) then A_update_RMOB_Dur_file=3:fp2(t0,", count=",hour_count,", strenght=",strenght_count,", maxdur=",dur_count_1,", noise=",n_1,", event=",str("mm",t1)):fclose2
124 if( A_update_RMOB_Dur_file=3 ) then A_update_RMOB_Dur_file=4:A_end_detect_end=1
125 if( A_update_RMOB_Dur_file=4 ) then A_update_RMOB_Dur_file=0:
126  
127 ;--------------A_end_detect_end-------------------------------------------------
128 ;zametení
129 if( A_end_detect_end=1 ) then A_end_detect_end=2:dur_count=0:dur_count_1=0:aver_amp_1=-999:aver_amp_max=-999:rec.trigger=0:A_still_detect=0
130 if( A_end_detect_end=2 ) then A_end_detect_end=3:REM exec(K_path+K_path_audio+"ftp_up.bat "+id_met2+".wav")
131 if( A_end_detect_end=3 ) then A_end_detect_end=4:REM exec(K_path+K_path_audio+"ftp_up.bat "+id_met2+".aux")
132 if( A_end_detect_end=4 ) then A_end_detect_end=0
133  
134 ;-----------------A_start_detect_bolid (událost start bolidu)-------------------
135 if( dur_count=K_min_dur_count_bolid ) then rec.filename=K_path+K_path_audio+id_met2+".wav":rec.trigger=1
136  
137 ;---A_met_screenshot (událost generování screenshotu waterfallu pokud je zaznamenán meteor)---
138 if( timer4.expired(1) ) then A_met_screenshot=1:capture (K_path+K_path_scr+id_met2+".jpg")
139 if( A_met_screenshot=1 ) then A_met_screenshot=2:REM exec(K_path+K_path_scr+"ftp_up.bat "+id_met2+".jpg")
140 if( A_met_screenshot=2 ) then A_met_screenshot=0:
141  
142 ;---A_no_met_screenshot (událost generování screenshotu waterfallu pokud není dlouho zaznamenán meteor)---
143 if( timer3.expired(1) ) then A_no_met_screenshot=1:t1_2=t1:capture (K_path+K_path_scr+K_station_name+str(t1_2)+"no_met.jpg")
144 if( A_no_met_screenshot=1 ) then A_no_met_screenshot=2:REM exec(K_path+K_path_scr+"ftp_up.bat "+K_station_name+str(t1_2)+"no_met.jpg")
145 if( A_no_met_screenshot=2 ) then A_no_met_screenshot=3:fopen1(K_path_SDR+K_station_name+t0+".dat",a)
146 if( A_no_met_screenshot=3 ) then A_no_met_screenshot=4:fp1(K_station_name+str(t1_2)+"no_met ; ; ; ; ; ; ; ")
147 if( A_no_met_screenshot=4 ) then A_no_met_screenshot=5:fclose1
148 if( A_no_met_screenshot=5 ) then A_no_met_screenshot=6:t1_2=0
149 if( A_no_met_screenshot=6 ) then A_no_met_screenshot=0:timer3.restart(K_max_time_to_scr)
150  
151 ;-A_once_hour_20s_before (událost 1x za hodinu před ukončením hodiny 20 vteřin)-
152 ;RMOB - soubor pro měření hodinových četností radiometeorů (projekt www.rmob.org)
153 ;Název souboru:RMOB-201301.dat, kde RMOB je zkratka pro Radio Meteor Observing Bulletin, 2013 je rok a 01 je měsíc
154 if( val(t3,"####")=5940 ) then A_once_hour_20s_before=A_once_hour_20s_before+1
155 if( A_once_hour_20s_before=1 ) then A_once_hour=1:fopen4(K_path_RMOB+"RMOB-"+str("YYYYMM",now)+".dat",a):fp4(t0,",",t2,",",hour_count):fclose4:t0_1=t0
156 if( A_once_hour=1 ) then A_once_hour=2:sp.print("Last hour=",hour_count):hour_count=0:strenght_count=0:strenght_count_1=0
157 if( A_once_hour=2 ) then A_once_hour=3:
158 if( A_once_hour=3 ) then A_once_hour=0:
159  
160 ;-A_once_hour_10s_before (událost 1x za hodinu před ukončením hodiny 10 vteřin)-
161 if( val(t3,"####")=5950 ) then A_once_hour_10s_before=1:A_once_hour_20s_before=0:
162 if( A_once_hour_10s_before=1 ) then A_once_hour_10s_before=0
163  
164 ;-A_once_hour_30s_after (událost 1x za hodinu po ukončení hodiny 30 vteřin)-
165 if( val(t3,"####")=0030 ) then A_once_hour_30s_after=A_once_hour_30s_after+1
166 if( A_once_hour_30s_after=1 ) then A_once_hour=1:exec(K_path_SDR+"ftp_up.bat "+K_station_name+t0_1+".dat")
167 if( A_once_hour=1 ) then A_once_hour=0
168  
169 ;-A_once_hour_40s_after (událost 1x za hodinu po ukončení hodiny 40 vteřin)-
170 if( val(t3,"####")=0040 ) then A_once_hour_40s_after=1:A_once_hour_30s_after=0
171 if( A_once_hour_40s_after=1 ) then A_once_hour_40s_after=2:t0_1=0
172 if( A_once_hour_40s_after=2 ) then A_once_hour_40s_after=0
173  
174  
175 ;Nahrada casoveho intervalu timery
176 ;timer5
177 ;timer6
178 ;timer7
3052 kaklik 179