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 9.2KB

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