StatisticsCalculator.java 14 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531
  1. import java.util.ArrayList;
  2. import java.util.Calendar;
  3. import java.util.Collections;
  4. import java.util.Map;
  5. import java.util.HashMap;
  6. /**
  7. * Berekent allerlei statistische waarden voor de gegeven data
  8. *
  9. * @author Projectgroep B5
  10. */
  11. public class StatisticsCalculator {
  12. /**
  13. * Berekent de maximale waarde
  14. *
  15. * @param array Een ArrayList met alle waarden
  16. *
  17. * @return De gevonden maximale waarde
  18. */
  19. public static double max(ArrayList<Double> array)
  20. {
  21. double max = 0;
  22. for(double waarde : array){
  23. if(waarde > max){
  24. max = waarde;
  25. }
  26. }
  27. return max;
  28. }
  29. /**
  30. * Berekent de minimale waarde
  31. *
  32. * @param array Een ArrayList met alle waarden
  33. *
  34. * @return De gevonden minimale waarde
  35. */
  36. public static double min(ArrayList<Double> array)
  37. {
  38. double min = array.get(0);
  39. for(double waarde : array){
  40. if(waarde < min){
  41. min = waarde;
  42. }
  43. }
  44. return min;
  45. }
  46. /**
  47. * Berekent de gemiddelde waarde
  48. *
  49. * @param array Een ArrayList met alle waarden
  50. *
  51. * @return De gevonden gemiddelde waarde
  52. */
  53. public static double avg(ArrayList<Double> array)
  54. {
  55. double avg = 0;
  56. for(double waarde : array){
  57. avg += waarde;
  58. }
  59. avg /= array.size();
  60. return avg;
  61. }
  62. /**
  63. * Berekent de mediaan
  64. *
  65. * @param array Een ArrayList met alle waarden
  66. *
  67. * @return De gevonden mediaan
  68. */
  69. public static double median(ArrayList<Double> array2){
  70. ArrayList<Double> array = new ArrayList<Double>();
  71. for(double db : array2)
  72. {
  73. array.add(db);
  74. }
  75. Collections.sort(array); //sort the array
  76. double median = 0;
  77. int middle = array.size()/2; //calculate the middle of the array
  78. if (array.size()%2 == 1) { //check if the array is even or uneven
  79. median = array.get(middle);
  80. } else {
  81. median = (array.get(middle-1) + array.get(middle+1)) / 2;
  82. }
  83. return median;
  84. }
  85. /**
  86. * Berekent de modus
  87. *
  88. * @param array Een ArrayList met alle waarden
  89. *
  90. * @return De gevonden modus
  91. */
  92. public static double modus(ArrayList<Double> array){
  93. Map<Double,Integer> map = new HashMap<Double,Integer>();
  94. for(double number : array)
  95. {
  96. if(map.containsKey(number)){
  97. map.put(number,(int)map.get(number)+1);
  98. }else{
  99. map.put(number, 1);
  100. }
  101. }
  102. double maxVal = 0;
  103. double maxKey = 0;
  104. for (Map.Entry<Double, Integer> entry : map.entrySet()) {
  105. if (entry.getValue() > maxVal) {
  106. maxVal = entry.getValue();
  107. maxKey = entry.getKey();
  108. }
  109. }
  110. return maxKey;
  111. }
  112. /**
  113. * Berekent de afwijking
  114. *
  115. * @param array Een ArrayList met alle waarden
  116. *
  117. * @return De gevonden afwijking
  118. */
  119. public static double afwijking(ArrayList<Double> array){
  120. double mediaan = StatisticsCalculator.median(array);
  121. double afwijking = 0;
  122. for(double m :array){
  123. afwijking += (mediaan-m)*(mediaan-m);
  124. }
  125. afwijking /= array.size();
  126. afwijking = Math.sqrt(afwijking);
  127. return afwijking;
  128. }
  129. public static double graadDagen(ArrayList<Double> array)
  130. {
  131. ArrayList<Double> avgTempDag = new ArrayList<Double>();
  132. int i = 0;
  133. double avgTemp = 0;
  134. //Breken de gemiddelde temperatuur per dag.
  135. for(double db : array)
  136. {
  137. i++;
  138. if(i%1440==0)
  139. {
  140. avgTempDag.add(avgTemp);
  141. avgTemp = 0;
  142. }
  143. avgTemp += db;
  144. }
  145. avgTemp /= array.size();
  146. int graaddagen = 0;
  147. int periodeBegin = 0;
  148. int periodeEind = 1439;
  149. for(int j = 0; j < avgTempDag.size(); j++)
  150. {
  151. if(avgTemp < 18)
  152. {
  153. graaddagen += 18 - avgTemp;
  154. }
  155. }
  156. periodeBegin += 1400;
  157. periodeEind += 1400;
  158. //graaddagen afronden
  159. return graaddagen;
  160. }
  161. public static int[] langsteDroogstePeriode(ArrayList<Double> array)
  162. {
  163. return langsteDroogstePeriodeMetMax(array, 0);
  164. }
  165. public static int[] langsteDroogstePeriodeMetMax(ArrayList<Double> array, int maxNeerslag)
  166. {
  167. int i = 0;
  168. double maximaleRegenval = 0;
  169. ArrayList<Double> maxRegenval = new ArrayList<Double>();
  170. for(double rainAmount : array)
  171. {
  172. i++;
  173. if(i%1440==0)
  174. {
  175. maxRegenval.add(maximaleRegenval);
  176. maximaleRegenval = 0;
  177. }
  178. if(rainAmount > maximaleRegenval)
  179. {
  180. maximaleRegenval = rainAmount;
  181. }
  182. }
  183. int periode = 0;
  184. int aantalDagen = 0;
  185. int langstePeriode = 0;
  186. ArrayList<Integer> periodes = new ArrayList<Integer>();
  187. for(int r=0; r<maxRegenval.size(); r++)
  188. {
  189. if(maxRegenval.get(r) > maxNeerslag)
  190. {
  191. r++;
  192. periodes.add(periode);
  193. periode = 0;
  194. }
  195. else
  196. {
  197. periode ++;
  198. if(periode > langstePeriode)
  199. {
  200. langstePeriode = periode;
  201. }
  202. }
  203. }
  204. return langstePeriode;
  205. }
  206. public static int[] langsteRegenPeriode(ArrayList<Double> array)
  207. {
  208. ArrayList<Double> rainday = new ArrayList<Double>();
  209. int i = 0;
  210. double longRainday = 0;
  211. // Bereken de Rainrate per dag
  212. for(double rain : array)
  213. {
  214. i++;
  215. if (i% 1440 ==0)
  216. {
  217. rainday.add(longRainday);
  218. longRainday = 0;
  219. }
  220. if (rain > longRainday)
  221. {
  222. longRainday = rain;
  223. }
  224. }
  225. int[] index = new int[2];
  226. int index1 = 0;
  227. int index2 = 0;
  228. int index1_1 = 0;
  229. boolean regen = false;
  230. int p = 0;
  231. int maxDays = 0;
  232. for(int t = 0; t < rainday.size(); t++)
  233. {
  234. if (rainday.get(t) > 0 )
  235. {
  236. p++;
  237. if(!regen)
  238. {
  239. regen = true;
  240. index1_1 = t;
  241. }
  242. else
  243. {
  244. if(p > maxDays)
  245. {
  246. maxDays = p;
  247. index1 = index1_1;
  248. index2 = t-1;
  249. regen=false;
  250. }
  251. p = 0;
  252. }
  253. }
  254. }
  255. index[0] = index1*1440;
  256. if(index2<0)
  257. {
  258. index2 =0;
  259. }
  260. index[1] = index2*1440;
  261. return index;
  262. }
  263. /**
  264. * Berekent de totale regenval
  265. *
  266. * @param array Een ArrayList met alle waarden
  267. *
  268. * @return De gevonden totale regenval
  269. */
  270. public static short maximaleRegenPeriode(ArrayList<Double> array)
  271. {
  272. short regen = 0;
  273. ArrayList<Short> regenPerUur = new ArrayList<Short>();
  274. for(int i=0; i<array.size();i++)
  275. {
  276. if(i%60==0)
  277. {
  278. regen /= 60;
  279. regenPerUur.add(regen);
  280. regen = 0;
  281. }
  282. regen += array.get(i);
  283. }
  284. short totaleRegen = 0;
  285. for(short sh : regenPerUur)
  286. {
  287. totaleRegen += sh;
  288. }
  289. return totaleRegen;
  290. }
  291. /**
  292. * Berekent de langste periode waarin het boven 25 graden celcius is geweest
  293. *
  294. * @param array Een ArrayList met alle waarden
  295. *
  296. * @return De posities van de begin en eindtijd in de array van de periode
  297. */
  298. public static int[] langsteZomersePeriode(ArrayList<Double> array)
  299. {
  300. ArrayList<Double> maxTempDag = new ArrayList<Double>();
  301. int i = 0;
  302. double maxTemp = 0;
  303. //Bereken de maximale temperatuur per dag.
  304. for(double db : array)
  305. {
  306. i++;
  307. if(i%1440==0)
  308. {
  309. maxTempDag.add( (double) maxTemp);
  310. maxTemp = 0;
  311. }
  312. if(db > maxTemp)
  313. {
  314. maxTemp = db;
  315. }
  316. }
  317. //Creer een array met twee indexen
  318. int[] index = new int[2];
  319. int index1 = 0;
  320. int index2 = 0;
  321. //Een tijdelijke index om bij te houden wat de index was in het geval dat de vorige groter was.
  322. int index1_1 = 0;
  323. //Een boolean
  324. boolean zomers = false;
  325. int p = 0;
  326. int maxDays = 0;
  327. //Doorloop je maximale temperaturen en zoek de langste periode
  328. for(int t=0; t<maxTempDag.size(); t++) {
  329. if(maxTempDag.get(t) > 770) {
  330. p++;
  331. if(!zomers)
  332. {
  333. zomers = true;
  334. index1_1 = t;
  335. }
  336. }
  337. else {
  338. if(p > maxDays) {
  339. maxDays = p;
  340. index1 = index1_1;
  341. index2 = t-1;
  342. zomers=false;
  343. }
  344. p=0;
  345. }
  346. }
  347. //Het terugsturen van de gevonden indexen. Je doet ze keer 1440 omdat je bij het berekenen van
  348. //de temperatuur per dag je het in stukken van 1440 samenvoegde.
  349. index[0] = index1*1440;
  350. if(index2<0){index2=0;}
  351. index[1] = index2*1440;
  352. return index;
  353. }
  354. public static Periode langsteHittegolfPeriode(ArrayList<Measurement> array){
  355. double maxDay = 0;
  356. ArrayList<Double> maxTempPerDay = new ArrayList<Double>();
  357. //calculate the max for every day
  358. for(int i = 0; i < array.size(); i++){
  359. if((i % 1440) == 0){
  360. maxTempPerDay.add(maxDay);
  361. maxDay = 0;
  362. }
  363. if(array.get(i).getRawOutsideTemp() > maxDay){
  364. maxDay = array.get(i).getRawOutsideTemp();
  365. }
  366. }
  367. int periodLength = 0;
  368. int maxPeriodLength = 5;
  369. int numberOfTropicalDays = 0;
  370. Calendar begin = Calendar.getInstance();
  371. Calendar eind = Calendar.getInstance();
  372. Periode periode = new Periode(begin, eind, "Geen hittegolfperiode gevonden");
  373. for(int i = 0; i < maxTempPerDay.size(); i++){
  374. if( maxTempPerDay.get(i) >= 770){ // If the temperature is bigger then 25 degree Celcius (770 degree Fahrenheit) the summer period will start
  375. periodLength++;
  376. if (maxTempPerDay.get(i) >= 860){ // If the temperature is bigger then 30 degree Celcius (860 degree Fahrenheit) it will add a tropical day
  377. numberOfTropicalDays++;
  378. }
  379. }else if(periodLength > maxPeriodLength && numberOfTropicalDays >= 3){ //period has ended. If it is bigger then 5 days AND there are 3 or more tropical days, the period is a heat wave
  380. maxPeriodLength = periodLength;
  381. periodLength = 0;
  382. numberOfTropicalDays = 0;
  383. eind.setTimeInMillis(array.get((i-1)*1440).getDateStamp().getTime());
  384. periode = new Periode(eind, begin, "Langste hittegolf periode");
  385. }else{
  386. numberOfTropicalDays = 0;
  387. periodLength = 0;
  388. }
  389. if(periodLength == 5){
  390. begin.setTimeInMillis(array.get((i-5)*1440).getDateStamp().getTime());
  391. }
  392. }
  393. return periode;
  394. }
  395. public static int[] langsteTempStijgingPeriode(ArrayList<Double> array)
  396. {
  397. ArrayList<Double> maxTempDag = new ArrayList<Double>();
  398. ArrayList<Integer> periodeLengteList = new ArrayList<Integer>();
  399. ArrayList<Integer> periodeStartList = new ArrayList<Integer>();
  400. int[] index = new int[2];
  401. int index1 = 0;
  402. int index2 = 0;
  403. int indexStart = 0;
  404. int indexEind = 1;
  405. int periodeLengte = 0;
  406. int i = 0;
  407. double maxTemp = 0;
  408. //berekent maximale temp. per dag
  409. for(double db : array)
  410. {
  411. i ++;
  412. if(i % 1440 == 0)
  413. {
  414. maxTempDag.add(maxTemp);
  415. maxTemp = 0;
  416. }
  417. if(db > maxTemp)
  418. {
  419. maxTemp = db;
  420. }
  421. }
  422. for(int t = 1; t < maxTempDag.size(); t ++) //zoekt naar periodes van temp.stijgingen
  423. {
  424. if (maxTempDag.get(t) > maxTempDag.get(t - 1))
  425. {
  426. indexEind ++;
  427. }
  428. if (maxTempDag.get(t) <= maxTempDag.get(t - 1))
  429. {
  430. periodeLengte = indexEind - indexStart;
  431. periodeLengteList.add(periodeLengte);
  432. periodeStartList.add(indexStart);
  433. indexStart = indexEind;
  434. }
  435. }
  436. Integer tempMax = Collections.max(periodeLengteList); //onthoudt de grootste periode van een temp.stijging
  437. for (int u = 0; u < periodeLengteList.size(); u ++) //pakt de index van de grootste periode
  438. {
  439. if (tempMax == periodeLengteList.get(u))
  440. {
  441. index1 = u;
  442. index2 = tempMax - u;
  443. }
  444. }
  445. index[0] = index1 * 1440;
  446. index[1] = index2 * 1440;
  447. return index;
  448. }
  449. }