{"id":938,"date":"2016-09-14T20:00:14","date_gmt":"2016-09-14T18:00:14","guid":{"rendered":"http:\/\/morony.pl\/?p=938"},"modified":"2016-09-13T23:46:18","modified_gmt":"2016-09-13T21:46:18","slug":"burze-w-polsce-edycja-2014","status":"publish","type":"post","link":"https:\/\/morony.pl\/?p=938","title":{"rendered":"Burze w Polsce, edycja 2014"},"content":{"rendered":"<p>Pod koniec 2013 roku wpad\u0142em na pomys\u0142, aby regularnie zapisywa\u0107 sobie obrazek (map\u0119) przedstawiaj\u0105cy aktualne lokalizacje wy\u0142adowa\u0144 atmosferycznych. Korzysta\u0142em wtedy ze strony <a href=\"http:\/\/www.industry.siemens.com\/services\/global\/de\/blids\/service\/spion\/Seiten\/spion_pl.aspx\">blide.de<\/a>. Oczywi\u015bcie nie mia\u0142em na my\u015bli robienia tego r\u0119cznie, wi\u0119c uruchomi\u0142em na serwerze skrypt, kt\u00f3ry r\u00f3wno co 15 minut zapisywa\u0142 aktualny obrazek. Skrypt dzia\u0142a do dzi\u015b i ca\u0142y czas zapisuje aktualne mapy. Przypomnia\u0142em sobie o nim, gdy na serwerze zacz\u0119\u0142o brakowa\u0107 miejsca na dane. Pobra\u0142em obrazki na sw\u00f3j komputer i zacz\u0105\u0142em zabaw\u0119. Na pierwszy ogie\u0144 poszed\u0142 pierwszy pe\u0142ny rok &#8211; 2014.<\/p>\n<p>Teoretycznie w ci\u0105gu 365 dni powinienem zapisa\u0107 365*24*4 = 35040 obrazk\u00f3w. Faktycznie zapisa\u0142em 34999, czyli ca\u0142kiem nie\u017ale. Pierwsza, oczywista my\u015bl to animacja. Po\u0142\u0105czy\u0142em wszystkie obrazki w  filmik:<\/p>\n<p><center><iframe loading=\"lazy\" width=\"560\" height=\"315\" src=\"https:\/\/www.youtube.com\/embed\/hvQwhKAp37U\" frameborder=\"0\" allowfullscreen><\/iframe><\/center><\/p>\n<p>W kolejnym korku chcia\u0142em podej\u015b\u0107 do tematu bardziej kompleksowo. Jak oszacowa\u0107 ilo\u015b\u0107 wy\u0142adowa\u0144 dla ka\u017cdej z map. Oczywi\u015bcie ilo\u015bci si\u0119 nie da, ale da si\u0119 oszacowa\u0107 obszar obj\u0119ty burz\u0105. Na mapach kolorem czerwonym oznaczone s\u0105 miejsca, gdzie wyst\u0105pi\u0142o co najmniej jedno wy\u0142adowanie w przeci\u0105gu ostatnich 15 minut. Niestety pliki s\u0105 zapisane w mocno skompresowanym jpg, ale wyra\u017anie wida\u0107 co czerwone a co nie. Napisa\u0142em na kolanie szybki skrypt w Matlabie:<\/p>\n<pre style=\"width:780px; overflow:auto;\">\r\nclear;clc\r\nfiles=dir('2014\/blids\/*.png');\r\nS = zeros(length(files), 2);\r\nfor a = 1:length(files)\r\n\tname = ['2014\/blids\/' files(a).name];\r\n\tC = strsplit(name, '\/');\r\n\tC2 = strsplit(C{end},'.');\r\n\tdate = C2{1};\r\n\ttime = C2{2};\r\n\tI = imread(name);\r\n\tidx = (I(:,:,1)>100 & I(:,:,2)<100 &#038; I(:,:,3)<100);\r\n\tS(a,1) = sum(sum(idx));\r\n\tS(a,2) = datenum(str2double(date(1:4)), str2double(date(5:6)), str2double(date(7:8)), str2double(time(1:2)), str2double(time(3:4)), 0);\r\n\t\r\n\tfprintf('%d\\t%d\\n',a,S(a,1));\r\nend\r\n<\/pre>\n<p>Dla ka\u017cdej mapy skrypt obliczy\u0142 i zapisa\u0142 ilo\u015b\u0107 pikseli, w kt\u00f3rych mia\u0142y miejsca wy\u0142adowania. Teraz wystarczy to \u0142adnie nanie\u015b\u0107 na wykres: np. sumarycznie dla ka\u017cdego dnia.<\/p>\n<pre style=\"width:780px; overflow:auto;\">\r\nclf\r\ndays = datenum(2014,01,01):datenum(2014,12,31);\r\nstorm = zeros(size(days));\r\n\r\nfor d = 1:length(days)\r\n\tstorm(d) = sum(S(S(:,2)>=days(d)&S(:,2)<days(d)+1,1));\r\nend\r\n%storm(storm>0)=log10(storm(storm>0));\r\nstairs(days-datenum(2014,01,01),storm, 'LineWidth',2)\r\naxis([0 365 0 2.5e5])\r\nhold on\r\nfor month = 1:12\r\n\tplot(ones(1,2)*datenum(2014,month,1)-datenum(2014,01,01),[0 2.5e5], 'k--')\r\nend\r\n<\/pre>\n<p>A oto efekt:<br \/>\n<center><a href=\"http:\/\/morony.pl\/wp-content\/uploads\/2016\/09\/001.png\"><img loading=\"lazy\" decoding=\"async\" width=\"1024\" height=\"768\" src=\"http:\/\/morony.pl\/wp-content\/uploads\/2016\/09\/001-1024x768.png\" alt=\"001\" style=\"width:600px;height:auto;\" class=\"alignnone size-large wp-image-941\" srcset=\"https:\/\/morony.pl\/wp-content\/uploads\/2016\/09\/001-1024x768.png 1024w, https:\/\/morony.pl\/wp-content\/uploads\/2016\/09\/001-300x225.png 300w, https:\/\/morony.pl\/wp-content\/uploads\/2016\/09\/001-768x576.png 768w\" sizes=\"(max-width: 1024px) 100vw, 1024px\" \/><\/a><\/center><\/p>\n<p>Pionowe linie oddzielaj\u0105 kolejne miesi\u0105ce. Wida\u0107 wzmo\u017con\u0105 aktywno\u015b\u0107 burzow\u0105 od kwietnia do wrze\u015bnia. Maksimum przypada na pocz\u0105tku sierpnia, co wida\u0107 na filmie. <\/p>\n<p>Jak wykorzystaliby\u015bcie takie dane? W wolnej chwili zrobi\u0119 animacj\u0119 i wykres dla 2015 roku. <\/p>\n","protected":false},"excerpt":{"rendered":"<p>Pod koniec 2013 roku wpad\u0142em na pomys\u0142, aby regularnie zapisywa\u0107 sobie obrazek (map\u0119) przedstawiaj\u0105cy aktualne lokalizacje wy\u0142adowa\u0144 atmosferycznych. Korzysta\u0142em wtedy ze strony blide.de. Oczywi\u015bcie nie mia\u0142em na my\u015bli robienia tego r\u0119cznie, wi\u0119c uruchomi\u0142em na serwerze skrypt, kt\u00f3ry r\u00f3wno co 15 minut zapisywa\u0142 aktualny obrazek. Skrypt dzia\u0142a do dzi\u015b i ca\u0142y czas zapisuje aktualne mapy. Przypomnia\u0142em [&hellip;]<\/p>\n","protected":false},"author":3,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[1],"tags":[],"_links":{"self":[{"href":"https:\/\/morony.pl\/index.php?rest_route=\/wp\/v2\/posts\/938"}],"collection":[{"href":"https:\/\/morony.pl\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/morony.pl\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/morony.pl\/index.php?rest_route=\/wp\/v2\/users\/3"}],"replies":[{"embeddable":true,"href":"https:\/\/morony.pl\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=938"}],"version-history":[{"count":4,"href":"https:\/\/morony.pl\/index.php?rest_route=\/wp\/v2\/posts\/938\/revisions"}],"predecessor-version":[{"id":943,"href":"https:\/\/morony.pl\/index.php?rest_route=\/wp\/v2\/posts\/938\/revisions\/943"}],"wp:attachment":[{"href":"https:\/\/morony.pl\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=938"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/morony.pl\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=938"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/morony.pl\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=938"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}