Projektteil 2
You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.

EpizentrumRechnung.java 8.1KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179
  1. package de.edotzlaff.detection.detektion.berechnung;
  2. import android.os.Build;
  3. import androidx.annotation.RequiresApi;
  4. import com.google.common.collect.ImmutableSet;
  5. import com.google.common.collect.Sets;
  6. import de.edotzlaff.detection.detektion.Detektionssignal;
  7. import de.edotzlaff.detection.detektion.Epizentrum;
  8. import de.edotzlaff.detection.detektion.berechnung.mathObjekte.Koordinaten;
  9. import de.edotzlaff.detection.detektion.berechnung.mathObjekte.Kreis;
  10. import de.edotzlaff.detection.detektion.berechnung.mathObjekte.Schnittpunkt;
  11. import java.math.BigDecimal;
  12. import java.math.RoundingMode;
  13. import java.time.LocalDateTime;
  14. import java.util.ArrayList;
  15. import java.util.Set;
  16. public class EpizentrumRechnung {
  17. //ToDo Schrittweite und Abstand auf Szenario anpassen
  18. private static final double SCHRITTWEITE = 0.02;
  19. private static final double MAX_ABSTAND = 0.009;
  20. private static final int DOUBLESTELLEN = 8;
  21. @RequiresApi(api = Build.VERSION_CODES.O)
  22. public ArrayList<Endpunkt> initEndpunkte(ArrayList<Detektionssignal> detektionssignal) {
  23. Detektionssignal nahegelegensteMessung = detektionssignal.get(0);
  24. ArrayList<Endpunkt> arrBerechnungsObjekte = new ArrayList<>();
  25. for (int i = 1; i < detektionssignal.size(); i++) {
  26. if (detektionssignal.get(i).getAnkunftsZeit().isBefore(nahegelegensteMessung.getAnkunftsZeit())) {
  27. nahegelegensteMessung = detektionssignal.get(i);
  28. }
  29. }
  30. for (Detektionssignal aufzeichnung : detektionssignal) {
  31. int zeitunterschied = aufzeichnung.getAnkunftsZeit().getHour() * 3600 + aufzeichnung.getAnkunftsZeit().getMinute() * 60 + aufzeichnung.getAnkunftsZeit().getSecond()
  32. - nahegelegensteMessung.getAnkunftsZeit().getHour() * 3600 - nahegelegensteMessung.getAnkunftsZeit().getMinute() * 60 - nahegelegensteMessung.getAnkunftsZeit().getSecond();
  33. double startradius = zeitunterschied * aufzeichnung.getGeschwErdbeben();
  34. Endpunkt tempEndpunkt = new Endpunkt(aufzeichnung.getDeviceId());
  35. Kreis k = new Kreis(aufzeichnung.getKoordinaten().getX(), aufzeichnung.getKoordinaten().getY(), startradius);
  36. tempEndpunkt.setKreis(k);
  37. arrBerechnungsObjekte.add(tempEndpunkt);
  38. }
  39. return arrBerechnungsObjekte;
  40. }
  41. public ArrayList<Schnittpunkt> berechneKreisSchnittpunkte(Kreis kreis1, Kreis kreis2) {
  42. ArrayList<Schnittpunkt> schnittpunkte = new ArrayList<>();
  43. double dx = kreis2.getMittelPunkt().getX() - kreis1.getMittelPunkt().getX();
  44. double dy = kreis2.getMittelPunkt().getY() - kreis1.getMittelPunkt().getY();
  45. double c = Math.sqrt(Math.pow(dx, 2) + Math.pow(dy, 2));
  46. double x = ((Math.pow(kreis1.getRadius(), 2) + Math.pow(c, 2) - Math.pow(kreis2.getRadius(), 2))) / (2 * c);
  47. double y = Math.pow(kreis1.getRadius(), 2) - Math.pow(x, 2);
  48. //Kein Schnittpunkt vorhanden
  49. if (y < 0) {
  50. return schnittpunkte;
  51. }
  52. if (y > 0) {
  53. y = Math.sqrt(y);
  54. }
  55. double ex0 = dx / c;
  56. double ex1 = dy / c;
  57. double ey0 = -ex1;
  58. Schnittpunkt s1 = new Schnittpunkt(kreis1.getMittelPunkt().getX() + x * ex0,
  59. kreis1.getMittelPunkt().getY() + x * ex1);
  60. //Genau ein Schnittpunkt vorhanden
  61. if (y == 0) {
  62. schnittpunkte.add(s1);
  63. return schnittpunkte;
  64. }
  65. Schnittpunkt s2 = new Schnittpunkt(s1.getKoordinaten().getX() - y * ey0, s1.getKoordinaten().getY() - y * ex0);
  66. s1.getKoordinaten().setLocation(s1.getKoordinaten().getX() + y * ey0, s1.getKoordinaten().getY() + y * ex0);
  67. schnittpunkte.add(s1);
  68. schnittpunkte.add(s2);
  69. return schnittpunkte;
  70. }
  71. public boolean genuegendSchnittpunktevorhanden(ArrayList<EndpunktVerbindung> endpunktVerbindungen) {
  72. for (EndpunktVerbindung endpunktVerbindung : endpunktVerbindungen) {
  73. ArrayList<Schnittpunkt> templist;
  74. templist = berechneKreisSchnittpunkte(endpunktVerbindung.getEndpunkt1().getKreis(), endpunktVerbindung.getEndpunkt2().getKreis());
  75. if (templist.isEmpty()) {
  76. return false;
  77. }
  78. endpunktVerbindung.setSchnittpunkte(templist);
  79. }
  80. return true;
  81. }
  82. public Epizentrum ueberpruefeObEpizentrum(ArrayList<EndpunktVerbindung> endpunktVerbindungen) {
  83. int counter;
  84. ArrayList<Schnittpunkt> schnittpunkte = new ArrayList<>();
  85. for (Schnittpunkt s1 : endpunktVerbindungen.get(0).getSchnittpunkte()) {
  86. counter = 0;
  87. for (int i = 1; i < endpunktVerbindungen.size(); i++) {
  88. boolean bereitsSchnittpunktErmittelt = false;
  89. for (Schnittpunkt s2 : endpunktVerbindungen.get(i).getSchnittpunkte()) {
  90. double abstand = Math.sqrt(
  91. Math.pow(s2.getKoordinaten().getX() - s1.getKoordinaten().getX(), 2) +
  92. Math.pow(s2.getKoordinaten().getY() - s1.getKoordinaten().getY(), 2));
  93. if (abstand < MAX_ABSTAND && !bereitsSchnittpunktErmittelt) {
  94. bereitsSchnittpunktErmittelt = true;
  95. counter++;
  96. schnittpunkte.add(s2);
  97. }
  98. if (counter == endpunktVerbindungen.size() - 1) {
  99. schnittpunkte.add(s1);
  100. Epizentrum tempEpizentrum = new Epizentrum();
  101. tempEpizentrum.setKoordinaten(berechneEpizentrumKoordinaten(schnittpunkte));
  102. return tempEpizentrum;
  103. }
  104. }
  105. }
  106. }
  107. return null;
  108. }
  109. public Koordinaten berechneEpizentrumKoordinaten(ArrayList<Schnittpunkt> schnittpunkte) {
  110. double x = 0;
  111. double y = 0;
  112. for (Schnittpunkt s : schnittpunkte) {
  113. x+= s.getKoordinaten().getX();
  114. y+= s.getKoordinaten().getY();
  115. }
  116. return new Koordinaten(new BigDecimal(x/schnittpunkte.size()).setScale(DOUBLESTELLEN, RoundingMode.HALF_UP).doubleValue(),
  117. new BigDecimal(y/schnittpunkte.size()).setScale(DOUBLESTELLEN, RoundingMode.HALF_UP).doubleValue());
  118. }
  119. public ArrayList<EndpunktVerbindung> ermittleEndpunktKonbinationen(ArrayList<Endpunkt> tempBerechungsObjekte) {
  120. ArrayList<EndpunktVerbindung> verbindungen = new ArrayList<>();
  121. Set<Set<Endpunkt>> combinations = Sets.combinations(ImmutableSet.copyOf(tempBerechungsObjekte), 2);
  122. for (Set<Endpunkt> tempEndpunkt : combinations) {
  123. ArrayList<Endpunkt> helpArrayToGetValues = new ArrayList<>(tempEndpunkt);
  124. EndpunktVerbindung endpunktVerbindung = new EndpunktVerbindung(helpArrayToGetValues.get(0), helpArrayToGetValues.get(1));
  125. verbindungen.add(endpunktVerbindung);
  126. }
  127. return verbindungen;
  128. }
  129. @RequiresApi(api = Build.VERSION_CODES.O)
  130. public LocalDateTime berechneEntstehungsZeitpunkt(Detektionssignal detektionssignal, Koordinaten koordEpi) {
  131. double abstand = Math.sqrt(
  132. Math.pow(detektionssignal.getKoordinaten().getX() - koordEpi.getX(), 2) +
  133. Math.pow(detektionssignal.getKoordinaten().getY() - koordEpi.getY(), 2));
  134. int sec = (int) abstand/ detektionssignal.geschwErdbeben;
  135. return detektionssignal.getAnkunftsZeit().minusSeconds(sec);
  136. }
  137. @RequiresApi(api = Build.VERSION_CODES.O)
  138. public LocalDateTime berechneEntstehungsZeitpunktLongLatt(Detektionssignal detektionssignal, Koordinaten koordEpi)
  139. {
  140. float[] results = new float[1];
  141. LocalDateTime.distanceBetween(detektionssignal.getKoordinaten().getX(), detektionssignal.getKoordinaten().getY(), koordEpi.getX(), koordEpi.getY());
  142. float distanceInMeters = results[0];
  143. int sec = (int) distanceInMeters/ detektionssignal.geschwErdbeben;
  144. return detektionssignal.getAnkunftsZeit().minusSeconds(sec);
  145. }
  146. public void erhoheRadius(ArrayList<Endpunkt> tempEndpunktListe) {
  147. for (Endpunkt b : tempEndpunktListe) {
  148. b.erhoeheRadius(SCHRITTWEITE);
  149. }
  150. }
  151. }