StatisticsCalculator.java 14 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527
  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 langstePeriode = 0;
  185. for(int r=0; r<maxRegenval.size(); r++)
  186. {
  187. if(maxRegenval.get(r) > maxNeerslag)
  188. {
  189. periode = 0;
  190. }
  191. else
  192. {
  193. periode ++;
  194. if(periode > langstePeriode)
  195. {
  196. langstePeriode = periode;
  197. }
  198. }
  199. }
  200. return langstePeriode;
  201. }
  202. public static int[] langsteRegenPeriode(ArrayList<Double> array)
  203. {
  204. ArrayList<Double> rainday = new ArrayList<Double>();
  205. int i = 0;
  206. double longRainday = 0;
  207. // Bereken de Rainrate per dag
  208. for(double rain : array)
  209. {
  210. i++;
  211. if (i% 1440 ==0)
  212. {
  213. rainday.add(longRainday);
  214. longRainday = 0;
  215. }
  216. if (rain > longRainday)
  217. {
  218. longRainday = rain;
  219. }
  220. }
  221. int[] index = new int[2];
  222. int index1 = 0;
  223. int index2 = 0;
  224. int index1_1 = 0;
  225. boolean regen = false;
  226. int p = 0;
  227. int maxDays = 0;
  228. for(int t = 0; t < rainday.size(); t++)
  229. {
  230. if (rainday.get(t) > 0 )
  231. {
  232. p++;
  233. if(!regen)
  234. {
  235. regen = true;
  236. index1_1 = t;
  237. }
  238. else
  239. {
  240. if(p > maxDays)
  241. {
  242. maxDays = p;
  243. index1 = index1_1;
  244. index2 = t-1;
  245. regen=false;
  246. }
  247. p = 0;
  248. }
  249. }
  250. }
  251. index[0] = index1*1440;
  252. if(index2<0)
  253. {
  254. index2 =0;
  255. }
  256. index[1] = index2*1440;
  257. return index;
  258. }
  259. /**
  260. * Berekent de totale regenval
  261. *
  262. * @param array Een ArrayList met alle waarden
  263. *
  264. * @return De gevonden totale regenval
  265. */
  266. public static short maximaleRegenPeriode(ArrayList<Double> array)
  267. {
  268. short regen = 0;
  269. ArrayList<Short> regenPerUur = new ArrayList<Short>();
  270. for(int i=0; i<array.size();i++)
  271. {
  272. if(i%60==0)
  273. {
  274. regen /= 60;
  275. regenPerUur.add(regen);
  276. regen = 0;
  277. }
  278. regen += array.get(i);
  279. }
  280. short totaleRegen = 0;
  281. for(short sh : regenPerUur)
  282. {
  283. totaleRegen += sh;
  284. }
  285. return totaleRegen;
  286. }
  287. /**
  288. * Berekent de langste periode waarin het boven 25 graden celcius is geweest
  289. *
  290. * @param array Een ArrayList met alle waarden
  291. *
  292. * @return De posities van de begin en eindtijd in de array van de periode
  293. */
  294. public static int[] langsteZomersePeriode(ArrayList<Double> array)
  295. {
  296. ArrayList<Double> maxTempDag = new ArrayList<Double>();
  297. int i = 0;
  298. double maxTemp = 0;
  299. //Bereken de maximale temperatuur per dag.
  300. for(double db : array)
  301. {
  302. i++;
  303. if(i%1440==0)
  304. {
  305. maxTempDag.add( (double) maxTemp);
  306. maxTemp = 0;
  307. }
  308. if(db > maxTemp)
  309. {
  310. maxTemp = db;
  311. }
  312. }
  313. //Creer een array met twee indexen
  314. int[] index = new int[2];
  315. int index1 = 0;
  316. int index2 = 0;
  317. //Een tijdelijke index om bij te houden wat de index was in het geval dat de vorige groter was.
  318. int index1_1 = 0;
  319. //Een boolean
  320. boolean zomers = false;
  321. int p = 0;
  322. int maxDays = 0;
  323. //Doorloop je maximale temperaturen en zoek de langste periode
  324. for(int t=0; t<maxTempDag.size(); t++) {
  325. if(maxTempDag.get(t) > 770) {
  326. p++;
  327. if(!zomers)
  328. {
  329. zomers = true;
  330. index1_1 = t;
  331. }
  332. }
  333. else {
  334. if(p > maxDays) {
  335. maxDays = p;
  336. index1 = index1_1;
  337. index2 = t-1;
  338. zomers=false;
  339. }
  340. p=0;
  341. }
  342. }
  343. //Het terugsturen van de gevonden indexen. Je doet ze keer 1440 omdat je bij het berekenen van
  344. //de temperatuur per dag je het in stukken van 1440 samenvoegde.
  345. index[0] = index1*1440;
  346. if(index2<0){index2=0;}
  347. index[1] = index2*1440;
  348. return index;
  349. }
  350. public static Periode langsteHittegolfPeriode(ArrayList<Measurement> array){
  351. int maxDay = 0;
  352. ArrayList<Measurement> maxTempPerDay = new ArrayList<Measurement>();
  353. int vorigeMeting = array.get(0).getDateStamp().getDay();
  354. //calculate the max for every day
  355. for(int i = 0; i < array.size(); i++){
  356. if(array.get(i).getDateStamp().getDay() != vorigeMeting){
  357. maxTempPerDay.add(array.get(maxDay));
  358. maxDay = i;
  359. vorigeMeting = array.get(i).getDateStamp().getDay();
  360. }
  361. if(array.get(i).getRawOutsideTemp() > array.get(maxDay).getRawOutsideTemp()){
  362. maxDay = i;
  363. }
  364. }
  365. int periodLength = 0;
  366. int maxPeriodLength = 5;
  367. int numberOfTropicalDays = 0;
  368. Calendar begin = Calendar.getInstance();
  369. Calendar eind = Calendar.getInstance();
  370. Periode periode = new Periode(begin, eind, "Geen hittegolfperiode gevonden");
  371. for(int i = 0; i < maxTempPerDay.size(); i++){
  372. if( maxTempPerDay.get(i).getRawOutsideTemp() >= 770){ // If the temperature is bigger then 25 degree Celcius (770 degree Fahrenheit) the summer period will start
  373. periodLength++;
  374. if (maxTempPerDay.get(i).getRawOutsideTemp() >= 860){ // If the temperature is bigger then 30 degree Celcius (860 degree Fahrenheit) it will add a tropical day
  375. numberOfTropicalDays++;
  376. }
  377. }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
  378. begin.setTimeInMillis(maxTempPerDay.get(i-periodLength).getDateStamp().getTime());
  379. eind.setTimeInMillis(maxTempPerDay.get(i).getDateStamp().getTime());
  380. maxPeriodLength = periodLength;
  381. periodLength = 0;
  382. numberOfTropicalDays = 0;
  383. periode = new Periode(eind, begin, "Langste hittegolf periode");
  384. }else{
  385. numberOfTropicalDays = 0;
  386. periodLength = 0;
  387. }
  388. }
  389. return periode;
  390. }
  391. public static int[] langsteTempStijgingPeriode(ArrayList<Double> array)
  392. {
  393. ArrayList<Double> maxTempDag = new ArrayList<Double>();
  394. ArrayList<Integer> periodeLengteList = new ArrayList<Integer>();
  395. ArrayList<Integer> periodeStartList = new ArrayList<Integer>();
  396. int[] index = new int[2];
  397. int index1 = 0;
  398. int index2 = 0;
  399. int indexStart = 0;
  400. int indexEind = 1;
  401. int periodeLengte = 0;
  402. int i = 0;
  403. double maxTemp = 0;
  404. //berekent maximale temp. per dag
  405. for(double db : array)
  406. {
  407. i ++;
  408. if(i % 1440 == 0)
  409. {
  410. maxTempDag.add(maxTemp);
  411. maxTemp = 0;
  412. }
  413. if(db > maxTemp)
  414. {
  415. maxTemp = db;
  416. }
  417. }
  418. for(int t = 1; t < maxTempDag.size(); t ++) //zoekt naar periodes van temp.stijgingen
  419. {
  420. if (maxTempDag.get(t) > maxTempDag.get(t - 1))
  421. {
  422. indexEind ++;
  423. }
  424. if (maxTempDag.get(t) <= maxTempDag.get(t - 1))
  425. {
  426. periodeLengte = indexEind - indexStart;
  427. periodeLengteList.add(periodeLengte);
  428. periodeStartList.add(indexStart);
  429. indexStart = indexEind;
  430. }
  431. }
  432. Integer tempMax = Collections.max(periodeLengteList); //onthoudt de grootste periode van een temp.stijging
  433. for (int u = 0; u < periodeLengteList.size(); u ++) //pakt de index van de grootste periode
  434. {
  435. if (tempMax == periodeLengteList.get(u))
  436. {
  437. index1 = u;
  438. index2 = tempMax - u;
  439. }
  440. }
  441. index[0] = index1 * 1440;
  442. index[1] = index2 * 1440;
  443. return index;
  444. }
  445. }