Sporting hospitals by distance

dist/lib/general/AeroflyMissionDescription.js

@@ -121,7 +121,7 @@ export default class AeroflyMissionDescription {
 
     for (const cp of this.#mission.checkpoints) {
       if (lastCp !== null) {
-        distance += this.getDistanceBetweenCheckpoints(lastCp, cp);
+        distance += AeroflyMissionDescription.getDistanceBetweenCheckpoints(lastCp, cp);
       }
 
       lastCp = cp;
@@ -181,7 +181,7 @@ export default class AeroflyMissionDescription {
    * @param {AeroflyMissionCheckpoint} cp
    * @returns {number} distance in meters
    */
-  getDistanceBetweenCheckpoints(lastCp, cp) {
+  static getDistanceBetweenCheckpoints(lastCp, cp) {
     const lat1 = (lastCp.latitude / 180) * Math.PI;
     const lon1 = (lastCp.longitude / 180) * Math.PI;
     const lat2 = (cp.latitude / 180) * Math.PI;

dist/lib/hems/Configuration.js

@@ -81,7 +81,8 @@ export class Configuration {
      * @type {string}
      */
     this.geoJsonFile =
-      positionals[0] ?? path.join(path.dirname(fileURLToPath(import.meta.url)), "../../data/hems/san_francisco.geojson");
+      positionals[0] ??
+      path.join(path.dirname(fileURLToPath(import.meta.url)), "../../data/hems/san_francisco.geojson");
     if (!path.isAbsolute(this.geoJsonFile)) {
       this.geoJsonFile = path.join(process.cwd(), this.geoJsonFile);
     }

dist/lib/hems/GeoJsonLocations.js

@@ -38,7 +38,6 @@ export class GeoJsonLocations {
       throw Error("Missing FeatureCollection with features in GeoJSON file");
     }
 
-
     const pointFeatures = featureCollection.features.filter((f) => {
       return f.type === "Feature" && f.geometry.type === "Point";
     });

dist/lib/hems/Scenario.js

@@ -43,7 +43,9 @@ export class Scenario {
 
     const isTransfer = this.configuration.canTransfer && locations.hospitals.length > 1 && Math.random() <= 0.1;
     const waypoint1 = this.#getRandLocation(isTransfer ? locations.hospitals : locations.other);
-    const waypoint2 = this.#getRandLocation(locations.hospitals, waypoint1);
+    const waypoint2 = isTransfer
+      ? this.#getRandLocation(locations.hospitals, waypoint1)
+      : this.#getNearestLocation(locations.hospitals, waypoint1);
 
     const conditions = new AeroflyMissionConditions({
       time,
@@ -135,6 +137,51 @@ export class Scenario {
     return location;
   }
 
+  /**
+   * @param {import('./GeoJsonLocations.js').GeoJsonFeature[]} locations
+   * @param {import('./GeoJsonLocations.js').GeoJsonFeature} location
+   * @returns {import('./GeoJsonLocations.js').GeoJsonFeature}
+   */
+  #getNearestLocation(locations, location) {
+    let distance = 9999;
+    let nearestLocation = locations[0];
+    for (const testLocation of locations) {
+      const testDistance = this.#getDistanceBetweenLocations(testLocation, location);
+      if (testDistance < distance) {
+        nearestLocation = testLocation;
+        distance = testDistance;
+      }
+    }
+
+    return nearestLocation;
+  }
+
+  /**
+   *
+   * @param {import('./GeoJsonLocations.js').GeoJsonFeature} lastCp
+   * @param {import('./GeoJsonLocations.js').GeoJsonFeature} cp
+   * @returns {number} distance in meters
+   */
+  #getDistanceBetweenLocations(lastCp, cp) {
+    const lat1 = (lastCp.geometry.coordinates[1] / 180) * Math.PI;
+    const lon1 = (lastCp.geometry.coordinates[0] / 180) * Math.PI;
+    const lat2 = (cp.geometry.coordinates[1] / 180) * Math.PI;
+    const lon2 = (cp.geometry.coordinates[0] / 180) * Math.PI;
+
+    const dLon = lon2 - lon1;
+    const dLat = lat2 - lat1;
+
+    //const y = Math.sin(dLon) * Math.cos(lat2);
+    //const x = Math.cos(lat1) * Math.sin(lat2) - Math.sin(lat1) * Math.cos(lat2) * Math.cos(dLon);
+    //const bearing = ((Math.atan2(y, x) * 180) / Math.PI + 360) % 360;
+
+    const a =
+      Math.sin(dLat / 2) * Math.sin(dLat / 2) +
+      Math.sin(dLon / 2) * Math.sin(dLon / 2) * Math.cos(lat1) * Math.cos(lat2);
+    const c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1 - a));
+    return 6_371_000 * c;
+  }
+
   /**
    *
    * @param {import('./GeoJsonLocations.js').GeoJsonFeature} location