Skip to content

Instantly share code, notes, and snippets.

@sxwebdev

sxwebdev/calculate_ds.js

Last active December 6, 2017 17:26
Show Gist options

  • Select an option

    Learn more about clone URLs
  • Save sxwebdev/62d77bcdd188097f4747f0b8e7c613bb to your computer and use it in GitHub Desktop.

Select an option

Learn more about clone URLs
Save sxwebdev/62d77bcdd188097f4747f0b8e7c613bb to your computer and use it in GitHub Desktop.
Download ZIP

Revisions

  1. sxwebdev revised this gist Dec 6, 2017. No changes.
  2. sxwebdev created this gist Dec 6, 2017.
    35 changes: 35 additions & 0 deletions calculate_ds.js
    View file
    Original file line number Diff line number Diff line change
    @@ -0,0 +1,35 @@
    const calculate_ds = (lat1, lon1, lat2, lon2, timestamp1, timestamp2) => {

    lat1 = lat1 * Math.PI / 180.0;
    lon1 = lon1 * Math.PI / 180.0;

    lat2 = lat2 * Math.PI / 180.0;
    lon2 = lon2 * Math.PI / 180.0;

    const r = 6378100;

    const rho1 = r * Math.cos(lat1);
    const z1 = r * Math.sin(lat1);
    const x1 = rho1 * Math.cos(lon1);
    const y1 = rho1 * Math.sin(lon1);

    const rho2 = r * Math.cos(lat2);
    const z2 = r * Math.sin(lat2);
    const x2 = rho2 * Math.cos(lon2);
    const y2 = rho2 * Math.sin(lon2);

    const dot = (x1 * x2 + y1 * y2 + z1 * z2);
    const cos_theta = dot / (r * r);

    const theta = Math.acos(cos_theta);

    const distance = r * theta;
    const time_s = (timestamp2 - timestamp1);
    const speed_mps = distance / time_s;
    const speed_kph = (speed_mps * 3600.0) / 1000.0;

    return { distance, time_s, speed_mps, speed_kph };

    }

    calculate_ds(53.232863,44.899063,53.232858,44.898095,1512512341,1512512346)