@wemap/geo 14.6.0 → 14.7.0

package/dist/cdn/wemap-geo.umd.js

@@ -0,0 +1 @@
+(function(e,t){typeof exports==`object`&&typeof module<`u`?t(exports,require("@wemap/maths")):typeof define==`function`&&define.amd?define([`exports`,`@wemap/maths`],t):(e=typeof globalThis<`u`?globalThis:e||self,t((e.Wemap=e.Wemap||{},e.Wemap.Geo={}),e.Wemap.Maths))})(this,function(e,t){Object.defineProperty(e,Symbol.toStringTag,{value:`Module`});var n=Object.defineProperty,r=(e,t)=>{let r={};for(var i in e)n(r,i,{get:e[i],enumerable:!0});return t||n(r,Symbol.toStringTag,{value:`Module`}),r},i=r({CIRCUMFERENCE:()=>_,EARTH_GRAVITY:()=>s,ECCENTRICITY:()=>d,ECCENTRICITY_2:()=>f,ELLIPSOID_FLATNESS:()=>u,EPS_DEG_MM:()=>c,EPS_MM:()=>l,R_MAJOR:()=>a,R_MAJOR_2:()=>p,R_MAJOR_4:()=>m,R_MINOR:()=>o,R_MINOR_2:()=>h,R_MINOR_4:()=>g,SECOND_ECCENTRICITY_2:()=>v}),a=6378137,o=6356752.314245179,s=9.80665,c=1e-8,l=.001,u=1/298.257223563,d=Math.sqrt(u*(2-u)),f=d*d,p=a*a,m=p*p,h=o*o,g=h*h,_=a*2*Math.PI,v=(p-h)/h,y=class e{static VERIFY_TYPING=!1;static checkType(e){if(e!==null&&!(typeof e==`number`&&!isNaN(e))){if(Array.isArray(e)&&e.length===2){let[t,n]=e;if(typeof t==`number`&&!isNaN(t)&&typeof n==`number`&&!isNaN(n)){if(t>n||t===n)throw Error(`Invalid level range: [${t}, ${n}]`);return}}throw Error(`Unknown level format: ${e}`)}}static isRange(t){return e.VERIFY_TYPING&&e.checkType(t),Array.isArray(t)}static clone(t){return e.VERIFY_TYPING&&e.checkType(t),t===null?null:typeof t==`number`?t:[t[0],t[1]]}static fromString(t){if(t===null)return null;if(typeof t!=`string`||!t.length)throw Error(`argument must be a non empty string, got ${typeof t}`);if(!isNaN(Number(t)))return parseFloat(t);let n=e=>{let n=e.map(e=>Number(e));if(n.some(e=>Number.isNaN(e))||e.some(e=>!e.length))throw Error(`Cannot parse following level: ${t}`);return n},r=t.split(`;`);if(r.length>1){let t=n(r),i=Math.min(...t),a=Math.max(...t);return e.checkType([i,a]),[i,a]}else{let r=t.substring(1).indexOf(`-`)+1;if(r>0){let[i,a]=n([t.substring(0,r),t.substring(r+1)]);return e.checkType([i,a]),[i,a]}}throw Error(`Cannot parse following level: ${t}`)}static parse(t){if(t==null)return null;if(typeof t==`number`){if(Number.isNaN(t))throw Error(`level number must not be NaN`);return t}if(Array.isArray(t)){if(t.length!==2)throw Error(`level range must have exactly 2 elements, got ${t.length}`);let n=[t[0],t[1]];return e.checkType(n),e.clone(n)}if(typeof t==`string`)return e.fromString(t);throw Error(`Cannot parse level from ${typeof t}`)}static contains(t,n){return e.VERIFY_TYPING&&(e.checkType(t),e.checkType(n)),t===n?!0:Array.isArray(t)?Array.isArray(n)?t[0]<=n[0]&&t[1]>=n[1]:n!==null&&t[0]<=n&&t[1]>=n:t===null||n===null?!1:t<=n[0]&&t>=n[1]}static intersection(t,n){if(e.VERIFY_TYPING&&(e.checkType(t),e.checkType(n)),t===null||n===null)return null;if(e.equals(t,n))return e.clone(t);if(typeof t==`number`&&typeof n==`number`)return t===n?t:null;if(Array.isArray(t)&&!Array.isArray(n))return e.contains(t,n)?n:null;if(!Array.isArray(t)&&Array.isArray(n))return e.contains(n,t)?t:null;let r=Math.max(t[0],n[0]),i=Math.min(t[1],n[1]);return i===r?i:i<r?null:[r,i]}static intersect(t,n){return e.VERIFY_TYPING&&(e.checkType(t),e.checkType(n)),t===null&&n===null||e.intersection(t,n)!==null}static union(t,n){if(e.VERIFY_TYPING&&(e.checkType(t),e.checkType(n)),t===n)return e.clone(t);if(n===null||t===null)return null;let r,i;return!Array.isArray(t)&&!Array.isArray(n)?(r=Math.min(t,n),i=Math.max(t,n)):Array.isArray(t)&&!Array.isArray(n)?(r=Math.min(t[0],n),i=Math.max(t[1],n)):!Array.isArray(t)&&Array.isArray(n)?(r=Math.min(n[0],t),i=Math.max(n[1],t)):(r=Math.min(t[0],n[0]),i=Math.max(t[1],n[1])),r===i?r:[r,i]}static multiplyBy(t,n){return e.VERIFY_TYPING&&e.checkType(t),t===null?null:Array.isArray(t)?[t[0]*n,t[1]*n]:t*n}static toString(t){return e.VERIFY_TYPING&&e.checkType(t),t===null?null:Array.isArray(t)?t[0]+`;`+t[1]:String(t)}static equals(t,n){return e.VERIFY_TYPING&&(e.checkType(t),e.checkType(n)),t===n?!0:Array.isArray(t)&&Array.isArray(n)?t[0]===n[0]&&t[1]===n[1]:!1}static diff(t,n){return e.VERIFY_TYPING&&(e.checkType(t),e.checkType(n)),t===null||n===null?null:!Array.isArray(t)&&!Array.isArray(n)?n-t:Array.isArray(t)&&!Array.isArray(n)?t[0]===n?n-t[1]:t[1]===n?n-t[0]:null:Array.isArray(n)&&!Array.isArray(t)?t===n[0]?n[1]-t:t===n[1]?n[0]-t:null:e.equals(t,n)?0:null}};function b(e,n){let r=(0,t.deg2rad)(e),i=(0,t.deg2rad)(n);return{sinLat:Math.sin(r),cosLat:Math.cos(r),sinLng:Math.sin(i),cosLng:Math.cos(i)}}function x(e){return a/Math.sqrt(1-f*e**2)}function S(e){return o/(1-f*e**2)}function C(e,t,n){let{sinLat:r,cosLat:i,sinLng:a,cosLng:o}=b(e,t),s=x(r),c=(s+n)*i;return[c*o,c*a,(s*(1-f)+n)*r]}function w(e){let[n,r,i]=e,s=Math.hypot(n,r),c=Math.atan2(a*i,o*s),l=Math.atan2(i+o*v*Math.sin(c)**3,s-a*f*Math.cos(c)**3),u=Math.sin(l),d=s*Math.cos(l)+i*u-a*Math.sqrt(1-f*u**2);return[(0,t.rad2deg)(l),(0,t.rad2deg)(Math.atan2(r,n)),d]}function T(e,n){let r=(0,t.deg2rad)(e),i=(0,t.deg2rad)(n),a=Math.cos(r);return[a*Math.cos(i),a*Math.sin(i),Math.sin(r)]}function E(e){return[(0,t.rad2deg)(Math.atan2(e[2],Math.hypot(e[0],e[1]))),(0,t.rad2deg)(Math.atan2(e[1],e[0]))]}function D(e,t){return 2*Math.asin(Math.hypot(e[0]-t[0],e[1]-t[1],e[2]-t[2])/2)}var O=class e{_lat;_lng;_alt=null;_level=null;_heightFromFloor=null;_heightFromGround=null;_ecef;_groundEcef=null;autoWrap=!0;constructor(e,t,n=null,r=null){this.lat=e,this.lng=t,this.alt=n,this.level=r,this._ecef=null}get lat(){return this._lat}set lat(e){if(Math.abs(e)<=90)this._lat=e;else throw Error(`lat argument is not in [-90; 90], value is ${e}`);this._ecef=null,this._groundEcef=null}get latitude(){return this._lat}set latitude(e){throw Error(`Please use Coordinates#lat setter instead of Coordinate#latitude`)}get lng(){return this._lng}set lng(e){this._lng=e,this.autoWrap&&this.wrap(),this._ecef=null,this._groundEcef=null}get longitude(){return this._lng}set longitude(e){throw Error(`Please use Coordinates#lng setter instead of Coordinate#longitude`)}get alt(){return this._alt}set alt(e){this._alt=e,this._ecef=null}get level(){return this._level}set level(e){y.checkType(e),this._level=e}get heightFromFloor(){return this._heightFromFloor}set heightFromFloor(e){this._heightFromFloor=e}get heightFromGround(){return this._heightFromGround}set heightFromGround(e){this._heightFromGround=e}clone(){let t=new e(this.lat,this.lng,this.alt);return this.level!==null&&(t.level=y.clone(this.level)),t}wrap(){(this._lng<=-180||this._lng>180)&&(this._lng=(0,t.wrap)(this._lng,-180,180))}static equals(e,t,n=c,r=l){return this.equalsWithoutLevel(e,t,n,r)?e===null||t===null||y.equals(e.level,t.level):!1}equals(t){return e.equals(this,t)}static equalsWithoutLevel(t,n,r=c,i=l){return t===null&&t===n?!0:!(t instanceof e)||!(n instanceof e)?!1:Math.abs(n.lat-t.lat)<r&&Math.abs(n.lng-t.lng)<r&&(t.alt===n.alt||t.alt!==null&&n.alt!==null&&Math.abs(n.alt-t.alt)<i)}equalsWithoutLevel(t,n=c,r=l){return e.equalsWithoutLevel(this,t,n,r)}destinationPoint(e,t,n=null){let r=this.clone();return r.move(e,t,n),r}move(e,t,n=null){let{sinLat:r,cosLat:i,sinLng:a,cosLng:o}=b(this.lat,this.lng),s=S(r),c=e/(2*s),l=2*s*Math.sin(c),u=l*Math.cos(c),d=u*Math.sin(t),p=u*Math.cos(t),m=-l*Math.sin(c),h=x(r),g=h*i,[_,v]=w([g*o-a*d-r*o*p+i*o*m,g*a+o*d-r*a*p+i*a*m,h*(1-f)*r+i*p+r*m]);if(this.lat=_,this.lng=v,n!==null){if(this.alt===null)throw Error(`Point altitude is not defined`);this.alt+=n}return this}distanceTo(e){let[n,r,i]=this.groundEcef,[a,o,s]=e.groundEcef,c=Math.hypot(a-n,o-r,s-i),l=S(Math.sin((0,t.deg2rad)((this.lat+e.lat)/2)));return 2*l*Math.asin(Math.min(1,c/(2*l)))}static distanceBetween(e,t){return e.distanceTo(t)}bearingTo(e){let t=this.groundEcef,n=e.groundEcef,r=n[0]-t[0],i=n[1]-t[1],a=n[2]-t[2],{sinLat:o,cosLat:s,sinLng:c,cosLng:l}=b(this.lat,this.lng),u=-c*r+l*i,d=-o*l*r-o*c*i+s*a;return Math.atan2(u,d)}static bearingTo(e,t){return e.bearingTo(t)}get enuToEcefRotation(){let e=t.Quaternion.fromAxisAngle([0,0,1],Math.PI/2+(0,t.deg2rad)(this.lng)),n=t.Quaternion.fromAxisAngle([1,0,0],Math.PI/2-(0,t.deg2rad)(this.lat));return t.Quaternion.multiply(e,n)}get ecefToEnuRotation(){let e=t.Quaternion.fromAxisAngle([1,0,0],(0,t.deg2rad)(this.lat)-Math.PI/2),n=t.Quaternion.fromAxisAngle([0,0,1],-(0,t.deg2rad)(this.lng)-Math.PI/2);return t.Quaternion.multiply(e,n)}get ecef(){return this._ecef||=C(this.lat,this.lng,this.alt||0),this._ecef}get groundEcef(){return this.alt?(this._groundEcef||=C(this.lat,this.lng,0),this._groundEcef):this.ecef}static fromECEF(n){if(t.Vector3.norm(n)===0)throw Error(`ECEF vector cannot be the earth center`);let[r,i,a]=w(n),o=new e(r,i,a);return o._ecef=n,o}getSegmentProjection(n,r){let i=T(n.lat,n.lng),a=T(r.lat,r.lng),o=T(this.lat,this.lng),s=t.Vector3.cross(i,a);if(t.Vector3.norm(s)===0)return null;let c=t.Vector3.cross(o,s),l=t.Vector3.normalize(t.Vector3.cross(s,c));if(Math.abs(D(i,a)-D(i,l)-D(a,l))>.001/6378137)return null;let[u,d]=E(l),f;return n.alt!==null&&r.alt!==null&&(f=(n.alt+r.alt)/2),new e(u,d,f,y.union(n.level,r.level))}toString(){let e=`[`+this._lat.toFixed(7)+`, `+this._lng.toFixed(7);return this._alt!==null&&(e+=`, `+this._alt.toFixed(2)),this._level!==null&&(e+=`, [`+y.toString(this._level)+`]`),e+=`]`,e}toJson(){return{lat:Number(this.lat.toFixed(8)),lng:Number(this.lng.toFixed(8)),...this.alt!==null&&{alt:Number(this.alt.toFixed(3))},...this.level!==null&&{level:this.level}}}static fromJson(t){return new e(t.lat,t.lng,t.alt,t.level)}toCompressedJson(){return this.level===null?[Number(this.lat.toFixed(8)),Number(this.lng.toFixed(8))]:[Number(this.lat.toFixed(8)),Number(this.lng.toFixed(8)),this.level]}static fromCompressedJson(t){let n=new e(t[0],t[1]);return t.length>2&&(n.level=t[2]),n}},k=class e extends O{_time=null;_accuracy=null;_bearing=null;constructor(e,t,n=null,r=null,i=null,a=null,o=null){super(e,t,n,r),this.time=i,this.accuracy=a,this.bearing=o}get time(){return this._time}set time(e){this._time=e}get accuracy(){return this._accuracy}set accuracy(e){if(e!==null&&e<0)throw Error(`accuracy argument is not a positive number`);this._accuracy=e}get bearing(){return this._bearing}set bearing(e){this._bearing=e===null?null:e%(2*Math.PI)}move(e,t,n=null){return super.move(e,t,n),this}destinationPoint(e,t,n=null){let r=this.clone();return r.move(e,t,n),r}static fromCoordinates(t){return new e(t.lat,t.lng,t.alt,t.level)}clone(){let t=e.fromCoordinates(super.clone());return t.time=this.time,t.accuracy=this.accuracy,t.bearing=this.bearing,t}static equals(t,n,r=c,i=l){return t===null&&t===n?!0:!(t instanceof e)||!(n instanceof e)||!super.equals(t,n,r,i)?!1:t.time===n.time&&t.accuracy===n.accuracy&&t.bearing===n.bearing}equals(t,n=c,r=l){return e.equals(this,t,n,r)}toJson(){return{...super.toJson(),...this.time!==null&&{time:this.time},...this.accuracy!==null&&{accuracy:this.accuracy},...this.bearing!==null&&{bearing:this.bearing}}}static fromJson(t){let n=e.fromCoordinates(O.fromJson(t));return t.time!==void 0&&(n.time=t.time),t.accuracy!==void 0&&(n.accuracy=t.accuracy),t.bearing!==void 0&&(n.bearing=t.bearing),n}},A=r({calcDistance:()=>F,createSegmentsAtLevel:()=>ee,geolocationPositionToUserPosition:()=>P,sampleRoute:()=>j,segmentIntersection:()=>te,simplifyRoute:()=>N,trimRoute:()=>M});function j(e,t=.7,n=0,r=Number.MAX_VALUE){let i=n+r,a=[],o,s=0,c=0,l=!1;for(let r=0;r<e.length-1;r++){let u=e[r],d=e[r+1],f=u.distanceTo(d),p=u.bearingTo(d),m=0;if(!l)if(n<s+f)l=!0,c=n-s;else{s+=f;continue}for(o=Object.assign(u.clone(),{bearing:p});m+c<f&&s+c<=i;){let e=o.destinationPoint(c,p);e.bearing=p,a.push(e),o=e,m+=c,s+=c,c=t}if(s+c>i)break;let h=f-m;s+=h,c-=h}return a}function M(e,t=e[0],n=Number.MAX_VALUE){let r=[],i=null,a,o=0;if(e.length<=1)throw Error(`Route must have at least 2 points`);for(a=1;a<e.length;a++){let n=e[a-1],o=e[a];if(O.equals(t,n)){r.push(n),i=n;break}let s=t.getSegmentProjection(n,o);if(s&&O.equals(t,s)&&!s.equals(o)){r.push(s),i=s;break}}if(!r.length)throw Error(`startPosition is not on the route`);for(;i&&a<e.length;){let t=e[a],s=i.distanceTo(t);if(o+s>=n||Math.abs(o+s-n)<=.001){let e=i.bearingTo(t),a=n-o,s=i.destinationPoint(a,e);r.push(s);break}r.push(t),i=t,o+=s,a++}return r}function N(e,n=(0,t.deg2rad)(5)){let r=e[0].equals(e[e.length-1]),i=e.slice(0,e.length-+!!r),a=i.length;for(let o=+!r;o<a;o++){let r=e[(0,t.positiveMod)(o-1,a)],s=e[o],c=e[(0,t.positiveMod)(o+1,a)],l=r.bearingTo(s),u=s.bearingTo(c);Math.abs(u-l)<n&&(i=i.filter(e=>e!==s))}return r&&i.push(i[0]),i}function P(e){if(e===null)return null;let{latitude:n,longitude:r,accuracy:i,heading:a}=e.coords,o=new k(n,r);return o.time=e.timestamp,o.accuracy=i,o.bearing=a?(0,t.deg2rad)(a):null,o}function F(e){return e.reduce((e,t,n,r)=>e+(n?r[n-1].distanceTo(t):0),0)}function ee(e,t,n=!0){let r=e.length,i=!1,a=[],o=[];for(let s=0;s<r;s++){let r=e[s],c=y.intersect(t,r.level);n&&!c&&i?o.push(r):c&&(i||(o=[],a.push(o),n&&s!==0&&o.push(e[s-1])),o.push(r)),i=c}return a}function te(e,t,n,r){let i=e=>e*Math.PI/180,o=(e.lat+t.lat+n.lat+r.lat)/4,s=a*180/Math.PI,c=a*180/Math.PI*Math.cos(i(o)),l=e=>({x:e.lng*c,y:e.lat*s}),u=l(e),d=l(t),f=l(n),p=l(r),m=d.x-u.x,h=d.y-u.y,g=p.x-f.x,_=p.y-f.y,v=m*_-h*g;if(Math.abs(v)<1e-10)return null;let y=f.x-u.x,b=f.y-u.y,x=(y*_-b*g)/v,S=(y*h-b*m)/v;return x<0||x>1||S<0||S>1?null:new O((u.y+x*h)/s,(u.x+x*m)/c)}var I=class e{northEast;southWest;constructor(e,t){if(this.northEast=e,this.southWest=t,this.northEast&&this.southWest&&this.getNorth()<this.getSouth())throw Error(`Incorrect bounding box`)}get center(){return new O((this.southWest.lat+this.northEast.lat)/2,(this.northEast.lng+this.southWest.lng)/2)}contains(e){return e.lat<=this.northEast.lat&&e.lat>=this.southWest.lat&&e.lng<=this.northEast.lng&&e.lng>=this.southWest.lng}extend(t){let n=this.southWest,r=this.northEast,i,a;if(t instanceof O)i=t,a=t;else if(t instanceof e)i=t.southWest,a=t.northEast;else throw Error(`Unknown parameter`);return this.southWest=new O(Math.min(i.lat,n.lat),Math.min(i.lng,n.lng)),this.northEast=new O(Math.max(a.lat,r.lat),Math.max(a.lng,r.lng)),this}extendsWithMeasure(e){if(typeof e!=`number`)throw Error(`measure is not a number`);return this.northEast=this.northEast.destinationPoint(e,0).move(e,Math.PI/2),this.southWest=this.southWest.clone().destinationPoint(e,-Math.PI/2).destinationPoint(e,Math.PI),this}pad(e){let t=this.southWest,n=this.northEast,r=Math.abs(t.lat-n.lat)*e,i=Math.abs(t.lng-n.lng)*e;return this.southWest=new O(t.lat-r,t.lng-i),this.northEast=new O(n.lat+r,n.lng+i),this}getSouthWest(){return this.southWest}getNorthEast(){return this.northEast}getNorthWest(){return new O(this.getNorth(),this.getWest())}getSouthEast(){return new O(this.getSouth(),this.getEast())}getWest(){return this.southWest.lng}getSouth(){return this.southWest.lat}getEast(){return this.northEast.lng}getNorth(){return this.northEast.lat}static equals(e,t){return O.equals(e.northEast,t.northEast)&&O.equals(e.southWest,t.southWest)}equals(t){return e.equals(this,t)}static fromArray(t){return new e(new O(t[3],t[2]),new O(t[1],t[0]))}static fromCoordinates(t){return t.length===0?null:t.reduce((e,t)=>e.extend(t),new e(t[0],t[0]))}toArray(){return[this.getWest(),this.getSouth(),this.getEast(),this.getNorth()]}},L=class e{x;y;z;time=null;_accuracy=null;_bearing=null;constructor(e,t,n,r=null,i=null,a=null){this.x=e,this.y=t,this.z=n,this.time=r,this.accuracy=i,this.bearing=a}get accuracy(){return this._accuracy}set accuracy(e){if(e!==null&&e<0)throw Error(`accuracy argument is not a positive number`);this._accuracy=e}get bearing(){return this._bearing}set bearing(e){this._bearing=e===null?null:e%(2*Math.PI)}clone(){return new e(this.x,this.y,this.z,this.time,this.accuracy,this.bearing)}static equals(t,n,r=l){return t===null&&t===n?!0:!(t instanceof e)||!(n instanceof e)?!1:Math.abs(n.x-t.x)<r&&Math.abs(n.y-t.y)<r&&Math.abs(n.z-t.z)<r&&t.time===n.time&&t.accuracy===n.accuracy&&t.bearing===n.bearing}equals(t){return e.equals(this,t)}toJson(){return{x:this.x,y:this.y,z:this.z,...this.time!==null&&{time:this.time},...this.accuracy!==null&&{accuracy:this.accuracy},...this.bearing!==null&&{bearing:this.bearing}}}static fromJson(t){return new e(t.x,t.y,t.z,t.time,t.accuracy,t.bearing)}},R=class extends L{},z=11102230246251565e-32,B=134217729,V=3.000000000000001*z;function H(e,t,n,r,i){let a,o,s,c,l=t[0],u=r[0],d=0,f=0;u>l==u>-l?(a=l,l=t[++d]):(a=u,u=r[++f]);let p=0;if(d<e&&f<n)for(u>l==u>-l?(o=l+a,s=a-(o-l),l=t[++d]):(o=u+a,s=a-(o-u),u=r[++f]),a=o,s!==0&&(i[p++]=s);d<e&&f<n;)u>l==u>-l?(o=a+l,c=o-a,s=a-(o-c)+(l-c),l=t[++d]):(o=a+u,c=o-a,s=a-(o-c)+(u-c),u=r[++f]),a=o,s!==0&&(i[p++]=s);for(;d<e;)o=a+l,c=o-a,s=a-(o-c)+(l-c),l=t[++d],a=o,s!==0&&(i[p++]=s);for(;f<n;)o=a+u,c=o-a,s=a-(o-c)+(u-c),u=r[++f],a=o,s!==0&&(i[p++]=s);return(a!==0||p===0)&&(i[p++]=a),p}function U(e,t){let n=t[0];for(let r=1;r<e;r++)n+=t[r];return n}function W(e){return new Float64Array(e)}var G=(3+16*z)*z,K=(2+12*z)*z,q=(9+64*z)*z*z,J=W(4),Y=W(8),X=W(12),Z=W(16),Q=W(4);function ne(e,t,n,r,i,a,o){let s,c,l,u,d,f,p,m,h,g,_,v,y,b,x,S,C,w,T=e-i,E=n-i,D=t-a,O=r-a;b=T*O,f=B*T,p=f-(f-T),m=T-p,f=B*O,h=f-(f-O),g=O-h,x=m*g-(b-p*h-m*h-p*g),S=D*E,f=B*D,p=f-(f-D),m=D-p,f=B*E,h=f-(f-E),g=E-h,C=m*g-(S-p*h-m*h-p*g),_=x-C,d=x-_,J[0]=x-(_+d)+(d-C),v=b+_,d=v-b,y=b-(v-d)+(_-d),_=y-S,d=y-_,J[1]=y-(_+d)+(d-S),w=v+_,d=w-v,J[2]=v-(w-d)+(_-d),J[3]=w;let k=U(4,J),A=K*o;if(k>=A||-k>=A||(d=e-T,s=e-(T+d)+(d-i),d=n-E,l=n-(E+d)+(d-i),d=t-D,c=t-(D+d)+(d-a),d=r-O,u=r-(O+d)+(d-a),s===0&&c===0&&l===0&&u===0)||(A=q*o+V*Math.abs(k),k+=T*u+O*s-(D*l+E*c),k>=A||-k>=A))return k;b=s*O,f=B*s,p=f-(f-s),m=s-p,f=B*O,h=f-(f-O),g=O-h,x=m*g-(b-p*h-m*h-p*g),S=c*E,f=B*c,p=f-(f-c),m=c-p,f=B*E,h=f-(f-E),g=E-h,C=m*g-(S-p*h-m*h-p*g),_=x-C,d=x-_,Q[0]=x-(_+d)+(d-C),v=b+_,d=v-b,y=b-(v-d)+(_-d),_=y-S,d=y-_,Q[1]=y-(_+d)+(d-S),w=v+_,d=w-v,Q[2]=v-(w-d)+(_-d),Q[3]=w;let j=H(4,J,4,Q,Y);b=T*u,f=B*T,p=f-(f-T),m=T-p,f=B*u,h=f-(f-u),g=u-h,x=m*g-(b-p*h-m*h-p*g),S=D*l,f=B*D,p=f-(f-D),m=D-p,f=B*l,h=f-(f-l),g=l-h,C=m*g-(S-p*h-m*h-p*g),_=x-C,d=x-_,Q[0]=x-(_+d)+(d-C),v=b+_,d=v-b,y=b-(v-d)+(_-d),_=y-S,d=y-_,Q[1]=y-(_+d)+(d-S),w=v+_,d=w-v,Q[2]=v-(w-d)+(_-d),Q[3]=w;let M=H(j,Y,4,Q,X);return b=s*u,f=B*s,p=f-(f-s),m=s-p,f=B*u,h=f-(f-u),g=u-h,x=m*g-(b-p*h-m*h-p*g),S=c*l,f=B*c,p=f-(f-c),m=c-p,f=B*l,h=f-(f-l),g=l-h,C=m*g-(S-p*h-m*h-p*g),_=x-C,d=x-_,Q[0]=x-(_+d)+(d-C),v=b+_,d=v-b,y=b-(v-d)+(_-d),_=y-S,d=y-_,Q[1]=y-(_+d)+(d-S),w=v+_,d=w-v,Q[2]=v-(w-d)+(_-d),Q[3]=w,Z[H(M,X,4,Q,Z)-1]}function re(e,t,n,r,i,a){let o=(t-a)*(n-i),s=(e-i)*(r-a),c=o-s,l=Math.abs(o+s);return Math.abs(c)>=G*l?c:-ne(e,t,n,r,i,a,l)}(7+56*z)*z,(3+28*z)*z,(26+288*z)*z*z,W(4),W(4),W(4),W(4),W(4),W(4),W(4),W(4),W(4),W(8),W(8),W(8),W(4),W(8),W(8),W(16),W(12),W(192),W(192),(10+96*z)*z,(4+48*z)*z,(44+576*z)*z*z,W(4),W(4),W(4),W(4),W(4),W(4),W(4),W(4),W(8),W(8),W(8),W(8),W(8),W(8),W(8),W(8),W(8),W(4),W(4),W(4),W(8),W(16),W(16),W(16),W(32),W(32),W(48),W(64),W(1152),W(1152),(16+224*z)*z,(5+72*z)*z,(71+1408*z)*z*z,W(4),W(4),W(4),W(4),W(4),W(4),W(4),W(4),W(4),W(4),W(24),W(24),W(24),W(24),W(24),W(24),W(24),W(24),W(24),W(24),W(1152),W(1152),W(1152),W(1152),W(1152),W(2304),W(2304),W(3456),W(5760),W(8),W(8),W(8),W(16),W(24),W(48),W(48),W(96),W(192),W(384),W(384),W(384),W(768),W(96),W(96),W(96),W(1152);function ie(e,t){var n,r,i=0,a,o,s,c,l,u,d,f=e[0],p=e[1],m=t.length;for(n=0;n<m;n++){r=0;var h=t[n],g=h.length-1;if(u=h[0],u[0]!==h[g][0]&&u[1]!==h[g][1])throw Error(`First and last coordinates in a ring must be the same`);for(o=u[0]-f,s=u[1]-p;r<g;r++){if(d=h[r+1],c=d[0]-f,l=d[1]-p,s===0&&l===0){if(c<=0&&o>=0||o<=0&&c>=0)return 0}else if(l>=0&&s<=0||l<=0&&s>=0){if(a=re(o,c,s,l,0,0),a===0)return 0;(a>0&&l>0&&s<=0||a<0&&l<=0&&s>0)&&i++}u=d,s=l,o=c}}return i%2!=0}function ae(e){if(!e)throw Error(`coord is required`);if(!Array.isArray(e)){if(e.type===`Feature`&&e.geometry!==null&&e.geometry.type===`Point`)return[...e.geometry.coordinates];if(e.type===`Point`)return[...e.coordinates]}if(Array.isArray(e)&&e.length>=2&&!Array.isArray(e[0])&&!Array.isArray(e[1]))return[...e];throw Error(`coord must be GeoJSON Point or an Array of numbers`)}function oe(e){return e.type===`Feature`?e.geometry:e}function se(e,t,n={}){if(!e)throw Error(`point is required`);if(!t)throw Error(`polygon is required`);let r=ae(e),i=oe(t),a=i.type,o=t.bbox,s=i.coordinates;if(o&&ce(r,o)===!1)return!1;a===`Polygon`&&(s=[s]);let c=!1;for(var l=0;l<s.length;++l){let e=ie(r,s[l]);if(e===0)return!n.ignoreBoundary;e&&(c=!0)}return c}function ce(e,t){return t[0]<=e[0]&&t[1]<=e[1]&&t[2]>=e[0]&&t[3]>=e[1]}var le=se,ue=class e{scale=1;heading=0;origin;buildingLevels;constructor(e,t=[]){this.origin=e,this.buildingLevels=t}localToWorld(e){let n=t.Vector3.multiplyScalar(e,this.scale),r=t.Quaternion.fromAxisAngle([0,0,1],-this.heading),i=t.Quaternion.multiply(this.origin.enuToEcefRotation,r),a=t.Quaternion.rotate(i,n),o=t.Vector3.sum(this.origin.ecef,a),s=O.fromECEF(o),c=s.alt;s.heightFromGround=s.alt-this.origin.alt;let l=e=>{s.level=e.id,s.heightFromFloor=c-e.floorAltitude};if(!this.buildingLevels.length)return s;let u=this.buildingLevels.filter(e=>c>=e.floorAltitude&&c<=e.ceilingAltitude);if(!u.length)return s;let d=u[0],f=[s.lng,s.lat],p=u.filter(e=>e.geometries.length);if(!p.length)return l(d),s;let m=p.filter(e=>e.geometries.some(e=>le(f,e)))[0];return m&&l(m),s}worldToLocal(e){let n=t.Quaternion.fromAxisAngle([0,0,1],this.heading),r=t.Quaternion.multiply(n,this.origin.ecefToEnuRotation),i=t.Vector3.subtract(e.ecef,this.origin.ecef),a=t.Vector3.multiplyScalar(i,1/this.scale);return t.Quaternion.rotate(r,a)}toJson(){return{origin:this.origin.toJson(),...this.scale!==1&&{scale:this.scale},...this.heading!==0&&{heading:this.heading}}}static fromJson(t){let n=new e(O.fromJson(t.origin));return n.scale=t.scale===void 0?1:t.scale,n.heading=t.heading===void 0?0:t.heading,n}},$=class e{_quaternion=[1,0,0,0];_heading=null;_eulerAngles=null;_time=null;_accuracy=null;constructor(e,t=null,n=null){this.quaternion=e,this.time=t,this.accuracy=n}static unitary(){return new e([1,0,0,0])}get quaternion(){return this._quaternion}set quaternion(e){if(Math.abs(1-t.Quaternion.norm(e))>1e-4)throw Error(`quaternion is not a unit quaternion`);this._quaternion=e,this._heading=null,this._eulerAngles=null}get time(){return this._time}set time(e){this._time=e}get accuracy(){return this._accuracy}set accuracy(e){if(e!==null&&(e<0||e>Math.PI))throw Error(`accuracy argument (`+e+`) is not in range [0; PI]`);this._accuracy=e}get eulerAngles(){return this._eulerAngles===null&&(this._eulerAngles=t.Rotations.quaternionToEulerZXY(this.quaternion)),this._eulerAngles}get eulerAnglesDegrees(){return this.eulerAngles.map(e=>(0,t.rad2deg)(e))}get heading(){if(this._heading===null){let e=0;typeof window<`u`&&window&&window.orientation&&(e=(0,t.deg2rad)(window.orientation)),this._heading=t.Rotations.getHeadingFromQuaternion(this.quaternion)+e}return this._heading}get headingDegrees(){return(0,t.rad2deg)(this.heading)}static equals(n,r){return n===null&&n===r?!0:!(n instanceof e)||!(r instanceof e)?!1:n===r||t.Quaternion.equals(n.quaternion,r.quaternion)}equals(t){return e.equals(this,t)}toJson(){return this.time===null&&this.accuracy===null?this.quaternion:{q:this.quaternion,...this.time!==null&&{time:this.time},...this.accuracy!==null&&{accuracy:this.accuracy}}}static fromJson(t){return Array.isArray(t)?new e(t,null,null):new e(t.q,t.time,t.accuracy)}clone(){return new e(this.quaternion.slice(0),this.time,this.accuracy)}static diff(n,r){let i=t.Quaternion.multiply(t.Quaternion.inverse(n.quaternion),r.quaternion),a=null;r.time!==null&&n.time!==null&&(a=r.time-n.time);let o=null;return n.accuracy!==null&&r.accuracy!==null&&(o=Math.max(r.accuracy-n.accuracy)),new e(i,a,o)}};e.AbsoluteHeading=class e{heading;time=null;_accuracy=null;constructor(e,t=null,n=null){this.heading=e,this.time=t,this.accuracy=n}get accuracy(){return this._accuracy}set accuracy(e){if(e!==null&&(e<0||e>Math.PI))throw Error(`accuracy argument (`+e+`) is not in range [0; PI]`);this._accuracy=e}toAttitude(){return new $(t.Quaternion.fromAxisAngle([0,0,1],-this.heading),this.time,this.accuracy)}static equals(t,n){return t===null&&t===n?!0:!(t instanceof e)||!(n instanceof e)?!1:Math.abs(t.heading-n.heading)<1e-8}equals(t){return e.equals(this,t)}toJson(){return{heading:this.heading,...this.time!==null&&{time:this.time},...this.accuracy!==null&&{accuracy:this.accuracy}}}static fromJson(t){return new e(t.heading,t.time,t.accuracy)}clone(){return new e(this.heading,this.time,this.accuracy)}},e.Attitude=$,e.BoundingBox=I,Object.defineProperty(e,"Constants",{enumerable:!0,get:function(){return i}}),e.Coordinates=O,e.GeoRef=ue,e.GeoRelativePosition=R,e.Level=y,e.RelativePosition=L,e.UserPosition=k,Object.defineProperty(e,"Utils",{enumerable:!0,get:function(){return A}})});

package/dist/index.js

@@ -21,8 +21,9 @@ var c = Object.defineProperty, l = (e, t) => {
 	R_MAJOR_4: () => b,
 	R_MINOR: () => f,
 	R_MINOR_2: () => x,
-	R_MINOR_4: () => S
-}), d = 6378137, f = 6356752.3142, p = 9.80665, m = 1e-8, h = .001, g = (d - f) / d, _ = Math.sqrt(g * (2 - g)), v = _ * _, y = d * d, b = y * y, x = f * f, S = x * x, C = d * 2 * Math.PI, w = class e {
+	R_MINOR_4: () => S,
+	SECOND_ECCENTRICITY_2: () => w
+}), d = 6378137, f = 6356752.314245179, p = 9.80665, m = 1e-8, h = .001, g = 1 / 298.257223563, _ = Math.sqrt(g * (2 - g)), v = _ * _, y = d * d, b = y * y, x = f * f, S = x * x, C = d * 2 * Math.PI, w = (y - x) / x, T = class e {
 	static VERIFY_TYPING = !1;
 	static checkType(e) {
 		if (e !== null && !(typeof e == "number" && !isNaN(e))) {
@@ -110,7 +111,55 @@ var c = Object.defineProperty, l = (e, t) => {
 	static diff(t, n) {
 		return e.VERIFY_TYPING && (e.checkType(t), e.checkType(n)), t === null || n === null ? null : !Array.isArray(t) && !Array.isArray(n) ? n - t : Array.isArray(t) && !Array.isArray(n) ? t[0] === n ? n - t[1] : t[1] === n ? n - t[0] : null : Array.isArray(n) && !Array.isArray(t) ? t === n[0] ? n[1] - t : t === n[1] ? n[0] - t : null : e.equals(t, n) ? 0 : null;
 	}
-}, T = class t {
+};
+//#endregion
+//#region src/coordinates/Coordinates.ts
+function E(e, t) {
+	let n = r(e), i = r(t);
+	return {
+		sinLat: Math.sin(n),
+		cosLat: Math.cos(n),
+		sinLng: Math.sin(i),
+		cosLng: Math.cos(i)
+	};
+}
+function D(e) {
+	return d / Math.sqrt(1 - v * e ** 2);
+}
+function O(e) {
+	return f / (1 - v * e ** 2);
+}
+function k(e, t, n) {
+	let { sinLat: r, cosLat: i, sinLng: a, cosLng: o } = E(e, t), s = D(r), c = (s + n) * i;
+	return [
+		c * o,
+		c * a,
+		(s * (1 - v) + n) * r
+	];
+}
+function A(e) {
+	let [t, n, r] = e, i = Math.hypot(t, n), o = Math.atan2(d * r, f * i), s = Math.atan2(r + f * w * Math.sin(o) ** 3, i - d * v * Math.cos(o) ** 3), c = Math.sin(s), l = i * Math.cos(s) + r * c - d * Math.sqrt(1 - v * c ** 2);
+	return [
+		a(s),
+		a(Math.atan2(n, t)),
+		l
+	];
+}
+function j(e, t) {
+	let n = r(e), i = r(t), a = Math.cos(n);
+	return [
+		a * Math.cos(i),
+		a * Math.sin(i),
+		Math.sin(n)
+	];
+}
+function M(e) {
+	return [a(Math.atan2(e[2], Math.hypot(e[0], e[1]))), a(Math.atan2(e[1], e[0]))];
+}
+function N(e, t) {
+	return 2 * Math.asin(Math.hypot(e[0] - t[0], e[1] - t[1], e[2] - t[2]) / 2);
+}
+var P = class t {
 	_lat;
 	_lng;
 	_alt = null;
@@ -118,6 +167,7 @@ var c = Object.defineProperty, l = (e, t) => {
 	_heightFromFloor = null;
 	_heightFromGround = null;
 	_ecef;
+	_groundEcef = null;
 	autoWrap = !0;
 	constructor(e, t, n = null, r = null) {
 		this.lat = e, this.lng = t, this.alt = n, this.level = r, this._ecef = null;
@@ -128,7 +178,7 @@ var c = Object.defineProperty, l = (e, t) => {
 	set lat(e) {
 		if (Math.abs(e) <= 90) this._lat = e;
 		else throw Error(`lat argument is not in [-90; 90], value is ${e}`);
-		this._ecef = null;
+		this._ecef = null, this._groundEcef = null;
 	}
 	get latitude() {
 		return this._lat;
@@ -140,7 +190,7 @@ var c = Object.defineProperty, l = (e, t) => {
 		return this._lng;
 	}
 	set lng(e) {
-		this._lng = e, this.autoWrap && this.wrap(), this._ecef = null;
+		this._lng = e, this.autoWrap && this.wrap(), this._ecef = null, this._groundEcef = null;
 	}
 	get longitude() {
 		return this._lng;
@@ -158,7 +208,7 @@ var c = Object.defineProperty, l = (e, t) => {
 		return this._level;
 	}
 	set level(e) {
-		w.checkType(e), this._level = e;
+		T.checkType(e), this._level = e;
 	}
 	get heightFromFloor() {
 		return this._heightFromFloor;
@@ -174,13 +224,13 @@ var c = Object.defineProperty, l = (e, t) => {
 	}
 	clone() {
 		let e = new t(this.lat, this.lng, this.alt);
-		return this.level !== null && (e.level = w.clone(this.level)), e;
+		return this.level !== null && (e.level = T.clone(this.level)), e;
 	}
 	wrap() {
 		(this._lng <= -180 || this._lng > 180) && (this._lng = o(this._lng, -180, 180));
 	}
 	static equals(e, t, n = m, r = h) {
-		return this.equalsWithoutLevel(e, t, n, r) ? e === null || t === null || w.equals(e.level, t.level) : !1;
+		return this.equalsWithoutLevel(e, t, n, r) ? e === null || t === null || T.equals(e.level, t.level) : !1;
 	}
 	equals(e) {
 		return t.equals(this, e);
@@ -196,23 +246,27 @@ var c = Object.defineProperty, l = (e, t) => {
 		return r.move(e, t, n), r;
 	}
 	move(e, t, n = null) {
-		let i = e / d, o = Math.cos(i), s = Math.sin(i), c = r(this.lat), l = r(this.lng), u = Math.asin(Math.sin(c) * o + Math.cos(c) * s * Math.cos(t)), f = l + Math.atan2(Math.sin(t) * s * Math.cos(c), o - Math.sin(c) * Math.sin(u));
-		if (this.lat = a(u), this.lng = a(f), n !== null) {
+		let { sinLat: r, cosLat: i, sinLng: a, cosLng: o } = E(this.lat, this.lng), s = O(r), c = e / (2 * s), l = 2 * s * Math.sin(c), u = l * Math.cos(c), d = u * Math.sin(t), f = u * Math.cos(t), p = -l * Math.sin(c), m = D(r), h = m * i, [g, _] = A([
+			h * o - a * d - r * o * f + i * o * p,
+			h * a + o * d - r * a * f + i * a * p,
+			m * (1 - v) * r + i * f + r * p
+		]);
+		if (this.lat = g, this.lng = _, n !== null) {
 			if (this.alt === null) throw Error("Point altitude is not defined");
 			this.alt += n;
 		}
 		return this;
 	}
 	distanceTo(e) {
-		let t = this.lat, n = this.lng, i = e.lat, a = e.lng, o = r(i - t), s = r(a - n), c = Math.sin(s / 2), l = Math.sin(o / 2), u = r(t), f = Math.cos(u), p = r(i), m = Math.cos(p), h = l * l + f * m * c * c, g = Math.sqrt(h), _ = Math.sqrt(1 - h);
-		return 2 * Math.atan2(g, _) * d;
+		let [t, n, i] = this.groundEcef, [a, o, s] = e.groundEcef, c = Math.hypot(a - t, o - n, s - i), l = O(Math.sin(r((this.lat + e.lat) / 2)));
+		return 2 * l * Math.asin(Math.min(1, c / (2 * l)));
 	}
 	static distanceBetween(e, t) {
 		return e.distanceTo(t);
 	}
 	bearingTo(e) {
-		let t = r(this.lat), n = r(e.lat), i = r(e.lng - this.lng);
-		return Math.atan2(Math.sin(i) * Math.cos(n), Math.cos(t) * Math.sin(n) - Math.sin(t) * Math.cos(n) * Math.cos(i));
+		let t = this.groundEcef, n = e.groundEcef, r = n[0] - t[0], i = n[1] - t[1], a = n[2] - t[2], { sinLat: o, cosLat: s, sinLng: c, cosLng: l } = E(this.lat, this.lng), u = -c * r + l * i, d = -o * l * r - o * c * i + s * a;
+		return Math.atan2(u, d);
 	}
 	static bearingTo(e, t) {
 		return e.bearingTo(t);
@@ -242,33 +296,27 @@ var c = Object.defineProperty, l = (e, t) => {
 		return e.multiply(t, n);
 	}
 	get ecef() {
-		if (!this._ecef) {
-			let e = r(this.lat), t = r(this.lng), n = this.alt || 0, i = (d + n) * Math.cos(e) * Math.cos(t), a = (d + n) * Math.cos(e) * Math.sin(t), o = (d + n) * Math.sin(e);
-			this._ecef = [
-				i,
-				a,
-				o
-			];
-		}
-		return this._ecef;
+		return this._ecef ||= k(this.lat, this.lng, this.alt || 0), this._ecef;
+	}
+	get groundEcef() {
+		return this.alt ? (this._groundEcef ||= k(this.lat, this.lng, 0), this._groundEcef) : this.ecef;
 	}
 	static fromECEF(e) {
-		let n = e[0], r = e[1], i = e[2], o = Math.sqrt(n ** 2 + r ** 2), s = Math.atan2(r, n), c = Math.atan2(i, o), l = o / Math.cos(c) - d;
-		s %= 2 * Math.PI;
-		let u = new t(a(c), a(s), l);
-		return u._ecef = e, u;
+		if (n.norm(e) === 0) throw Error("ECEF vector cannot be the earth center");
+		let [r, i, a] = A(e), o = new t(r, i, a);
+		return o._ecef = e, o;
 	}
 	getSegmentProjection(e, r) {
-		let i = n.normalize(e.ecef), a = n.normalize(r.ecef), o = n.normalize(this.ecef), s = n.cross(i, a);
+		let i = j(e.lat, e.lng), a = j(r.lat, r.lng), o = j(this.lat, this.lng), s = n.cross(i, a);
 		if (n.norm(s) === 0) return null;
-		let c = n.cross(o, s), l = n.normalize(n.cross(s, c)), u = n.multiplyScalar(l, d), f = t.fromECEF(u), p;
-		e.alt !== null && r.alt !== null && (p = (e.alt + r.alt) / 2);
-		let m = new t(f.lat, f.lng, p, w.union(e.level, r.level));
-		return Math.abs(e.distanceTo(r) - e.distanceTo(m) - r.distanceTo(m)) > .001 ? null : m;
+		let c = n.cross(o, s), l = n.normalize(n.cross(s, c));
+		if (Math.abs(N(i, a) - N(i, l) - N(a, l)) > .001 / 6378137) return null;
+		let [u, d] = M(l), f;
+		return e.alt !== null && r.alt !== null && (f = (e.alt + r.alt) / 2), new t(u, d, f, T.union(e.level, r.level));
 	}
 	toString() {
 		let e = "[" + this._lat.toFixed(7) + ", " + this._lng.toFixed(7);
-		return this._alt !== null && (e += ", " + this._alt.toFixed(2)), this._level !== null && (e += ", [" + w.toString(this._level) + "]"), e += "]", e;
+		return this._alt !== null && (e += ", " + this._alt.toFixed(2)), this._level !== null && (e += ", [" + T.toString(this._level) + "]"), e += "]", e;
 	}
 	toJson() {
 		return {
@@ -292,7 +340,7 @@ var c = Object.defineProperty, l = (e, t) => {
 		let n = new t(e[0], e[1]);
 		return e.length > 2 && (n.level = e[2]), n;
 	}
-}, E = class e extends T {
+}, F = class e extends P {
 	_time = null;
 	_accuracy = null;
 	_bearing = null;
@@ -347,19 +395,19 @@ var c = Object.defineProperty, l = (e, t) => {
 		};
 	}
 	static fromJson(t) {
-		let n = e.fromCoordinates(T.fromJson(t));
+		let n = e.fromCoordinates(P.fromJson(t));
 		return t.time !== void 0 && (n.time = t.time), t.accuracy !== void 0 && (n.accuracy = t.accuracy), t.bearing !== void 0 && (n.bearing = t.bearing), n;
 	}
-}, D = /* @__PURE__ */ l({
-	calcDistance: () => M,
-	createSegmentsAtLevel: () => N,
-	geolocationPositionToUserPosition: () => j,
-	sampleRoute: () => O,
-	segmentIntersection: () => P,
-	simplifyRoute: () => A,
-	trimRoute: () => k
+}, I = /* @__PURE__ */ l({
+	calcDistance: () => V,
+	createSegmentsAtLevel: () => H,
+	geolocationPositionToUserPosition: () => B,
+	sampleRoute: () => L,
+	segmentIntersection: () => U,
+	simplifyRoute: () => z,
+	trimRoute: () => R
 });
-function O(e, t = .7, n = 0, r = Number.MAX_VALUE) {
+function L(e, t = .7, n = 0, r = Number.MAX_VALUE) {
 	let i = n + r, a = [], o, s = 0, c = 0, l = !1;
 	for (let r = 0; r < e.length - 1; r++) {
 		let u = e[r], d = e[r + 1], f = u.distanceTo(d), p = u.bearingTo(d), m = 0;
@@ -378,17 +426,17 @@ function O(e, t = .7, n = 0, r = Number.MAX_VALUE) {
 	}
 	return a;
 }
-function k(e, t = e[0], n = Number.MAX_VALUE) {
+function R(e, t = e[0], n = Number.MAX_VALUE) {
 	let r = [], i = null, a, o = 0;
 	if (e.length <= 1) throw Error("Route must have at least 2 points");
 	for (a = 1; a < e.length; a++) {
 		let n = e[a - 1], o = e[a];
-		if (T.equals(t, n)) {
+		if (P.equals(t, n)) {
 			r.push(n), i = n;
 			break;
 		}
 		let s = t.getSegmentProjection(n, o);
-		if (s && T.equals(t, s) && !s.equals(o)) {
+		if (s && P.equals(t, s) && !s.equals(o)) {
 			r.push(s), i = s;
 			break;
 		}
@@ -405,7 +453,7 @@ function k(e, t = e[0], n = Number.MAX_VALUE) {
 	}
 	return r;
 }
-function A(e, t = r(5)) {
+function z(e, t = r(5)) {
 	let n = e[0].equals(e[e.length - 1]), a = e.slice(0, e.length - +!!n), o = a.length;
 	for (let r = +!n; r < o; r++) {
 		let n = e[i(r - 1, o)], s = e[r], c = e[i(r + 1, o)], l = n.bearingTo(s), u = s.bearingTo(c);
@@ -413,51 +461,51 @@ function A(e, t = r(5)) {
 	}
 	return n && a.push(a[0]), a;
 }
-function j(e) {
+function B(e) {
 	if (e === null) return null;
-	let { latitude: t, longitude: n, accuracy: i, heading: a } = e.coords, o = new E(t, n);
+	let { latitude: t, longitude: n, accuracy: i, heading: a } = e.coords, o = new F(t, n);
 	return o.time = e.timestamp, o.accuracy = i, o.bearing = a ? r(a) : null, o;
 }
-function M(e) {
+function V(e) {
 	return e.reduce((e, t, n, r) => e + (n ? r[n - 1].distanceTo(t) : 0), 0);
 }
-function N(e, t, n = !0) {
+function H(e, t, n = !0) {
 	let r = e.length, i = !1, a = [], o = [];
 	for (let s = 0; s < r; s++) {
-		let r = e[s], c = w.intersect(t, r.level);
+		let r = e[s], c = T.intersect(t, r.level);
 		n && !c && i ? o.push(r) : c && (i || (o = [], a.push(o), n && s !== 0 && o.push(e[s - 1])), o.push(r)), i = c;
 	}
 	return a;
 }
-function P(e, t, n, r) {
+function U(e, t, n, r) {
 	let i = (e) => e * Math.PI / 180, a = (e.lat + t.lat + n.lat + r.lat) / 4, o = d * 180 / Math.PI, s = d * 180 / Math.PI * Math.cos(i(a)), c = (e) => ({
 		x: e.lng * s,
 		y: e.lat * o
 	}), l = c(e), u = c(t), f = c(n), p = c(r), m = u.x - l.x, h = u.y - l.y, g = p.x - f.x, _ = p.y - f.y, v = m * _ - h * g;
 	if (Math.abs(v) < 1e-10) return null;
 	let y = f.x - l.x, b = f.y - l.y, x = (y * _ - b * g) / v, S = (y * h - b * m) / v;
-	return x < 0 || x > 1 || S < 0 || S > 1 ? null : new T((l.y + x * h) / o, (l.x + x * m) / s);
+	return x < 0 || x > 1 || S < 0 || S > 1 ? null : new P((l.y + x * h) / o, (l.x + x * m) / s);
 }
 //#endregion
 //#region src/coordinates/BoundingBox.ts
-var F = class e {
+var W = class e {
 	northEast;
 	southWest;
 	constructor(e, t) {
 		if (this.northEast = e, this.southWest = t, this.northEast && this.southWest && this.getNorth() < this.getSouth()) throw Error("Incorrect bounding box");
 	}
 	get center() {
-		return new T((this.southWest.lat + this.northEast.lat) / 2, (this.northEast.lng + this.southWest.lng) / 2);
+		return new P((this.southWest.lat + this.northEast.lat) / 2, (this.northEast.lng + this.southWest.lng) / 2);
 	}
 	contains(e) {
 		return e.lat <= this.northEast.lat && e.lat >= this.southWest.lat && e.lng <= this.northEast.lng && e.lng >= this.southWest.lng;
 	}
 	extend(t) {
 		let n = this.southWest, r = this.northEast, i, a;
-		if (t instanceof T) i = t, a = t;
+		if (t instanceof P) i = t, a = t;
 		else if (t instanceof e) i = t.southWest, a = t.northEast;
 		else throw Error("Unknown parameter");
-		return this.southWest = new T(Math.min(i.lat, n.lat), Math.min(i.lng, n.lng)), this.northEast = new T(Math.max(a.lat, r.lat), Math.max(a.lng, r.lng)), this;
+		return this.southWest = new P(Math.min(i.lat, n.lat), Math.min(i.lng, n.lng)), this.northEast = new P(Math.max(a.lat, r.lat), Math.max(a.lng, r.lng)), this;
 	}
 	extendsWithMeasure(e) {
 		if (typeof e != "number") throw Error("measure is not a number");
@@ -465,7 +513,7 @@ var F = class e {
 	}
 	pad(e) {
 		let t = this.southWest, n = this.northEast, r = Math.abs(t.lat - n.lat) * e, i = Math.abs(t.lng - n.lng) * e;
-		return this.southWest = new T(t.lat - r, t.lng - i), this.northEast = new T(n.lat + r, n.lng + i), this;
+		return this.southWest = new P(t.lat - r, t.lng - i), this.northEast = new P(n.lat + r, n.lng + i), this;
 	}
 	getSouthWest() {
 		return this.southWest;
@@ -474,10 +522,10 @@ var F = class e {
 		return this.northEast;
 	}
 	getNorthWest() {
-		return new T(this.getNorth(), this.getWest());
+		return new P(this.getNorth(), this.getWest());
 	}
 	getSouthEast() {
-		return new T(this.getSouth(), this.getEast());
+		return new P(this.getSouth(), this.getEast());
 	}
 	getWest() {
 		return this.southWest.lng;
@@ -492,13 +540,13 @@ var F = class e {
 		return this.northEast.lat;
 	}
 	static equals(e, t) {
-		return T.equals(e.northEast, t.northEast) && T.equals(e.southWest, t.southWest);
+		return P.equals(e.northEast, t.northEast) && P.equals(e.southWest, t.southWest);
 	}
 	equals(t) {
 		return e.equals(this, t);
 	}
 	static fromArray(t) {
-		return new e(new T(t[3], t[2]), new T(t[1], t[0]));
+		return new e(new P(t[3], t[2]), new P(t[1], t[0]));
 	}
 	static fromCoordinates(t) {
 		return t.length === 0 ? null : t.reduce((e, t) => e.extend(t), new e(t[0], t[0]));
@@ -511,7 +559,7 @@ var F = class e {
 			this.getNorth()
 		];
 	}
-}, I = class e {
+}, G = class e {
 	x;
 	y;
 	z;
@@ -556,7 +604,7 @@ var F = class e {
 	static fromJson(t) {
 		return new e(t.x, t.y, t.z, t.time, t.accuracy, t.bearing);
 	}
-}, L = class extends I {}, R = class t {
+}, K = class extends G {}, q = class t {
 	scale = 1;
 	heading = 0;
 	origin;
@@ -569,7 +617,7 @@ var F = class e {
 			0,
 			0,
 			1
-		], this.heading), a = e.multiply(i, this.origin.enuToEcefRotation), o = e.rotate(a, r), c = n.sum(this.origin.ecef, o), l = T.fromECEF(c), u = l.alt;
+		], -this.heading), a = e.multiply(this.origin.enuToEcefRotation, i), o = e.rotate(a, r), c = n.sum(this.origin.ecef, o), l = P.fromECEF(c), u = l.alt;
 		l.heightFromGround = l.alt - this.origin.alt;
 		let d = (e) => {
 			l.level = e.id, l.heightFromFloor = u - e.floorAltitude;
@@ -587,7 +635,7 @@ var F = class e {
 			0,
 			0,
 			1
-		], -this.heading), i = e.multiply(this.origin.ecefToEnuRotation, r), a = n.subtract(t.ecef, this.origin.ecef), o = n.multiplyScalar(a, 1 / this.scale);
+		], this.heading), i = e.multiply(r, this.origin.ecefToEnuRotation), a = n.subtract(t.ecef, this.origin.ecef), o = n.multiplyScalar(a, 1 / this.scale);
 		return e.rotate(i, o);
 	}
 	toJson() {
@@ -598,10 +646,10 @@ var F = class e {
 		};
 	}
 	static fromJson(e) {
-		let n = new t(T.fromJson(e.origin));
+		let n = new t(P.fromJson(e.origin));
 		return n.scale = e.scale === void 0 ? 1 : e.scale, n.heading = e.heading === void 0 ? 0 : e.heading, n;
 	}
-}, z = class n {
+}, J = class n {
 	_quaternion = [
 		1,
 		0,
@@ -684,7 +732,7 @@ var F = class e {
 		let o = null;
 		return t.accuracy !== null && r.accuracy !== null && (o = Math.max(r.accuracy - t.accuracy)), new n(i, a, o);
 	}
-}, B = class t {
+}, Y = class t {
 	heading;
 	time = null;
 	_accuracy = null;
@@ -699,7 +747,7 @@ var F = class e {
 		this._accuracy = e;
 	}
 	toAttitude() {
-		return new z(e.fromAxisAngle([
+		return new J(e.fromAxisAngle([
 			0,
 			0,
 			1
@@ -726,4 +774,4 @@ var F = class e {
 	}
 };
 //#endregion
-export { B as AbsoluteHeading, z as Attitude, F as BoundingBox, u as Constants, T as Coordinates, R as GeoRef, L as GeoRelativePosition, w as Level, I as RelativePosition, E as UserPosition, D as Utils };
+export { Y as AbsoluteHeading, J as Attitude, W as BoundingBox, u as Constants, P as Coordinates, q as GeoRef, K as GeoRelativePosition, T as Level, G as RelativePosition, F as UserPosition, I as Utils };

package/dist/src/Constants.d.ts

@@ -1,5 +1,9 @@
 export declare const R_MAJOR = 6378137;
-export declare const R_MINOR = 6356752.3142;
+/**
+ * R_MAJOR * (1 - ELLIPSOID_FLATNESS), to the last bit. Kept as a literal so it stays a
+ * literal type for consumers, but the flattening below is the primitive one.
+ */
+export declare const R_MINOR = 6356752.314245179;
 export declare const EARTH_GRAVITY = 9.80665;
 /**
  * latitude and longitude epsilon in degrees
@@ -10,6 +14,7 @@ export declare const EPS_DEG_MM = 1e-8;
  * epsilon in meters which corresponds to 1 millimeter
  */
 export declare const EPS_MM = 0.001;
+/** WGS84 defines the flattening, and derives the semi-minor axis from it. */
 export declare const ELLIPSOID_FLATNESS: number;
 export declare const ECCENTRICITY: number;
 export declare const ECCENTRICITY_2: number;
@@ -18,4 +23,9 @@ export declare const R_MAJOR_4: number;
 export declare const R_MINOR_2: number;
 export declare const R_MINOR_4: number;
 export declare const CIRCUMFERENCE: number;
+/**
+ * Second (polar) eccentricity squared, e'^2 = (a^2 - b^2) / b^2.
+ * Used by the geodetic <-> ECEF conversion.
+ */
+export declare const SECOND_ECCENTRICITY_2: number;
 //# sourceMappingURL=Constants.d.ts.map

package/dist/src/Constants.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"Constants.d.ts","sourceRoot":"","sources":["../../src/Constants.ts"],"names":[],"mappings":"AAAA,eAAO,MAAM,OAAO,UAAY,CAAC;AACjC,eAAO,MAAM,OAAO,eAAe,CAAC;AACpC,eAAO,MAAM,aAAa,UAAU,CAAC;AAErC;;;GAGG;AACH,eAAO,MAAM,UAAU,OAAO,CAAC;AAE/B;;GAEG;AACH,eAAO,MAAM,MAAM,QAAO,CAAC;AAG3B,eAAO,MAAM,kBAAkB,QAAgC,CAAC;AAEhE,eAAO,MAAM,YAAY,QAA2D,CAAC;AACrF,eAAO,MAAM,cAAc,QAA8B,CAAC;AAE1D,eAAO,MAAM,SAAS,QAAoB,CAAC;AAC3C,eAAO,MAAM,SAAS,QAAwB,CAAC;AAC/C,eAAO,MAAM,SAAS,QAAoB,CAAC;AAC3C,eAAO,MAAM,SAAS,QAAwB,CAAC;AAC/C,eAAO,MAAM,aAAa,QAAwB,CAAC"}
+{"version":3,"file":"Constants.d.ts","sourceRoot":"","sources":["../../src/Constants.ts"],"names":[],"mappings":"AAAA,eAAO,MAAM,OAAO,UAAY,CAAC;AACjC;;;GAGG;AACH,eAAO,MAAM,OAAO,oBAAoB,CAAC;AACzC,eAAO,MAAM,aAAa,UAAU,CAAC;AAErC;;;GAGG;AACH,eAAO,MAAM,UAAU,OAAO,CAAC;AAE/B;;GAEG;AACH,eAAO,MAAM,MAAM,QAAO,CAAC;AAG3B,6EAA6E;AAC7E,eAAO,MAAM,kBAAkB,QAAoB,CAAC;AAEpD,eAAO,MAAM,YAAY,QAA2D,CAAC;AACrF,eAAO,MAAM,cAAc,QAA8B,CAAC;AAE1D,eAAO,MAAM,SAAS,QAAoB,CAAC;AAC3C,eAAO,MAAM,SAAS,QAAwB,CAAC;AAC/C,eAAO,MAAM,SAAS,QAAoB,CAAC;AAC3C,eAAO,MAAM,SAAS,QAAwB,CAAC;AAC/C,eAAO,MAAM,aAAa,QAAwB,CAAC;AAEnD;;;GAGG;AACH,eAAO,MAAM,qBAAqB,QAAsC,CAAC"}

package/dist/src/Utils.d.ts

@@ -24,9 +24,19 @@ export declare function calcDistance(coords: Coordinates[]): number;
 export declare function createSegmentsAtLevel(itineraryCoords: Coordinates[], segmentsLevel: Level_t, useMultiLevelSegments?: boolean): Coordinates[][];
 /**
  * Computes the intersection point of two WGS84 segments.
- * Uses a local planar projection in meters (valid approximation for short-to-medium distances).
  * Returns null if segments are parallel or do not intersect.
  * Note: it does not use level information.
+ *
+ * The projection below is linear, so it maps lines to lines: the intersection point comes out
+ * the same for any kLat / kLng, spherical or ellipsoidal, because the factors cancel on the
+ * way out. They are kept because they make the code read as metres, and because the parallel
+ * test is not scale invariant: denom is in projected units, so its 1e-10 threshold means
+ * square metres. Make that epsilon relative before changing the factors.
+ *
+ * What the result does approximate is treating the segments as straight in lat/lng rather than
+ * as geodesics, which is off by 0.6mm on 100m segments and grows with the square of their
+ * length (5.6cm at 1km, 5.6m at 10km). Only Edge.intersection uses it, on graph edges, so
+ * that is well inside the noise.
  */
 export declare function segmentIntersection(a1: Coordinates, a2: Coordinates, b1: Coordinates, b2: Coordinates): Coordinates | null;
 //# sourceMappingURL=Utils.d.ts.map

package/dist/src/Utils.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"Utils.d.ts","sourceRoot":"","sources":["../../src/Utils.ts"],"names":[],"mappings":"AAKA,OAAO,WAAW,MAAM,2BAA2B,CAAC;AACpD,OAAO,YAAY,MAAM,4BAA4B,CAAC;AAEtD,OAAO,KAAK,EAAE,OAAO,EAAE,MAAM,SAAS,CAAC;AAEvC,MAAM,MAAM,WAAW,GAAG,WAAW,GAAG;IAAE,OAAO,EAAE,MAAM,CAAA;CAAE,CAAC;AAE5D;;;;;GAKG;AACH,wBAAgB,WAAW,CACvB,KAAK,EAAE,WAAW,EAAE,EACpB,QAAQ,SAAM,EACd,aAAa,SAAI,EACjB,MAAM,SAAmB,iBAwD5B;AAED;;;;;GAKG;AACH,wBAAgB,SAAS,CAAC,KAAK,EAAE,WAAW,EAAE,EAAE,aAAa,cAAW,EAAE,MAAM,SAAmB,iBAqDlG;AAGD,wBAAgB,aAAa,CAAC,MAAM,EAAE,WAAW,EAAE,EAAE,cAAc,SAAa,iBA0B/E;AAED,wBAAgB,iCAAiC,CAAC,mBAAmB,EAAE,mBAAmB,uBAYzF;AAED,wBAAgB,YAAY,CAAC,MAAM,EAAE,WAAW,EAAE,UAGjD;AAGD,wBAAgB,qBAAqB,CAAC,eAAe,EAAE,WAAW,EAAE,EAAE,aAAa,EAAE,OAAO,EAAE,qBAAqB,UAAO,mBA2BzH;AAED;;;;;GAKG;AACH,wBAAgB,mBAAmB,CAC/B,EAAE,EAAE,WAAW,EAAE,EAAE,EAAE,WAAW,EAChC,EAAE,EAAE,WAAW,EAAE,EAAE,EAAE,WAAW,GAC7B,WAAW,GAAG,IAAI,CAmCxB"}
+{"version":3,"file":"Utils.d.ts","sourceRoot":"","sources":["../../src/Utils.ts"],"names":[],"mappings":"AAKA,OAAO,WAAW,MAAM,2BAA2B,CAAC;AACpD,OAAO,YAAY,MAAM,4BAA4B,CAAC;AAEtD,OAAO,KAAK,EAAE,OAAO,EAAE,MAAM,SAAS,CAAC;AAEvC,MAAM,MAAM,WAAW,GAAG,WAAW,GAAG;IAAE,OAAO,EAAE,MAAM,CAAA;CAAE,CAAC;AAE5D;;;;;GAKG;AACH,wBAAgB,WAAW,CACvB,KAAK,EAAE,WAAW,EAAE,EACpB,QAAQ,SAAM,EACd,aAAa,SAAI,EACjB,MAAM,SAAmB,iBAwD5B;AAED;;;;;GAKG;AACH,wBAAgB,SAAS,CAAC,KAAK,EAAE,WAAW,EAAE,EAAE,aAAa,cAAW,EAAE,MAAM,SAAmB,iBAqDlG;AAGD,wBAAgB,aAAa,CAAC,MAAM,EAAE,WAAW,EAAE,EAAE,cAAc,SAAa,iBA0B/E;AAED,wBAAgB,iCAAiC,CAAC,mBAAmB,EAAE,mBAAmB,uBAYzF;AAED,wBAAgB,YAAY,CAAC,MAAM,EAAE,WAAW,EAAE,UAGjD;AAGD,wBAAgB,qBAAqB,CAAC,eAAe,EAAE,WAAW,EAAE,EAAE,aAAa,EAAE,OAAO,EAAE,qBAAqB,UAAO,mBA2BzH;AAED;;;;;;;;;;;;;;;GAeG;AACH,wBAAgB,mBAAmB,CAC/B,EAAE,EAAE,WAAW,EAAE,EAAE,EAAE,WAAW,EAChC,EAAE,EAAE,WAAW,EAAE,EAAE,EAAE,WAAW,GAC7B,WAAW,GAAG,IAAI,CAmCxB"}

package/dist/src/coordinates/Coordinates.d.ts

@@ -7,9 +7,20 @@ import { CoordinatesCompressedJson, CoordinatesJson, Level_t } from '../types';
  * Basic geo methods are directly accessibles from here:
  * distanceTo, bearingTo, toEcef...
  *
- * /!\ This class has been adapted to use earth as a sphere and not as an ellipsoid
- * /!\ So, this class does not stricly represent WGS84 coordinates anymore
- * /!\ This modifications have been made for computational improvements.
+ * Coordinates are WGS84 and every transformation here uses the ellipsoid, so the local up
+ * axis is the ellipsoid normal and the ECEF / ENU geometry reads `alt` as a height above the
+ * ellipsoid.
+ *
+ * What `alt` actually carries is not settled across the SDK, and that predates this class:
+ * the position providers put a height above the ground in it (Constants.DEFAULT_ALTITUDE is
+ * 1.6m, a phone in a hand), while GeoRef.localToWorld compares it against the absolute
+ * floorAltitude / ceilingAltitude of a building level. Both cannot be right. Only differences
+ * matter to the geometry here, which is why the ambiguity has never surfaced as a bug, but
+ * level detection does depend on the answer.
+ *
+ * (Until 2026-08 this class modelled the earth as a sphere of radius R_MAJOR for speed. That
+ * cost ~0.35% of anisotropic scale error at mid-latitude and tilted the local vertical by
+ * 11.5', i.e. ~33cm of altitude error per 100m of horizontal travel.)
  */
 declare class Coordinates {
     private _lat;
@@ -19,6 +30,7 @@ declare class Coordinates {
     private _heightFromFloor;
     private _heightFromGround;
     private _ecef;
+    private _groundEcef;
     autoWrap: boolean;
     constructor(lat: number, lng: number, alt?: number | null, level?: Level_t | null);
     get lat(): number;
@@ -32,6 +44,7 @@ declare class Coordinates {
     /**
      * alt does not denote the altitude of a point but its height from
      *  the "level" field (if defined) or from the ground
+     * /!\ This is the providers' convention, not GeoRef's: see the class header.
      */
     get alt(): number | null;
     set alt(alt: number | null);
@@ -55,33 +68,62 @@ declare class Coordinates {
      */
     destinationPoint(distance: number, bearing: number, elevation?: number | null): Coordinates;
     /**
-     * Source: http://www.movable-type.co.uk/scripts/latlong.html#destPoint
+     * Moves the point by a ground distance (in meters) along a bearing, on the WGS84 ellipsoid.
+     * The displacement is built in the local ENU frame as the chord of the arc travelled, then
+     * converted back through ECEF.
+     * ponytail: micrometer-accurate at working range (7um at 10km), then degrades as
+     * distance^3 / radius^2 (7mm at 100km, ~20m at 1200km). Use Karney's geodesics if
+     * continental distances ever matter.
      * @throws {Error} if elevation is defined and point altitude is not defined
      */
     move(distance: number, bearing: number, elevation?: number | null): this;
     /**
-     * Returns a distance between two points in meters
+     * Returns the ground distance between two points in meters, on the WGS84 ellipsoid.
+     * Altitudes are ignored, as many callers compare positions across levels.
+     *
+     * The chord between the two ECEF points is exact; converting it back to an arc length on
+     * the locally fitting sphere is what costs accuracy, and it costs very little: 3nm at
+     * 340m, 3um at 14km, 1.2m at 557km.
+     *
+     * The one case where this is worse than the old haversine is a path along the equator,
+     * where the arc radius is R_MAJOR but the Gaussian radius is R_MINOR: 264m over 3340km,
+     * 9.3km over 10000km, against an exact haversine there. Everywhere else, and at every
+     * distance this SDK works at, it is orders of magnitude better.
      */
     distanceTo(location2: Coordinates): number;
     static distanceBetween(point1: Coordinates, point2: Coordinates): number;
+    /**
+     * Azimuth to another point, in radians, clockwise from north.
+     * Taken in the local ENU frame, so it matches the true geodesic azimuth within 0.02" even
+     * at 100km.
+     */
     bearingTo(location2: Coordinates): number;
     static bearingTo(point1: Coordinates, point2: Coordinates): number;
     /**
-     * ECEF Transformations
-     * Here we used a light version of ECEF considering earth
-     * as a sphere instead of an ellipse
+     * ECEF Transformations (WGS84 ellipsoid)
      */
     get enuToEcefRotation(): Quaternion_t;
     get ecefToEnuRotation(): Quaternion_t;
     /**
      * https://gist.github.com/klucar/1536194
-     * Adapted for spherical formula
      */
     get ecef(): Vector3_t;
+    /**
+     * ECEF at altitude 0, for the ground-level geometry of distanceTo, move and bearingTo.
+     * The cached ecef already is that vector whenever alt is null or 0, which is both the
+     * common case and the hot one: the routing graph measures the same vertices thousands of
+     * times, and OsrmRemoteRouter walks hundreds of millions of pairs on a long route.
+     */
+    private get groundEcef();
     static fromECEF(ecef: Vector3_t): Coordinates;
     /**
      * https://stackoverflow.com/questions/1299567/how-to-calculate-distance-from-a-point-to-a-line-segment-on-a-sphere
-     * Adapted to ECEF
+     *
+     * Great-circle geometry, so it stays on the sphere: normalizing an ellipsoidal ECEF
+     * vector would yield a geocentric direction and shift the projection by tens of meters.
+     * The closure test that rejects a point off the segment measures angles between the same
+     * unit vectors, for the same reason: mixing in the geodesic distanceTo would compare a
+     * great circle with a geodesic and reject valid projections on long segments.
      */
     getSegmentProjection(p1: Coordinates, p2: Coordinates): Coordinates | null;
     /**

package/dist/src/coordinates/Coordinates.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"Coordinates.d.ts","sourceRoot":"","sources":["../../../src/coordinates/Coordinates.ts"],"names":[],"mappings":"AAAA,OAAO,EAC0C,KAAK,SAAS,EAAE,KAAK,YAAY,EACjF,MAAM,cAAc,CAAC;AAGtB,OAAO,KAAK,EAAE,yBAAyB,EAAE,eAAe,EAAE,OAAO,EAAE,MAAM,UAAU,CAAC;AAGpF;;;;;;;;;;GAUG;AACH,cAAM,WAAW;IAEb,OAAO,CAAC,IAAI,CAAU;IACtB,OAAO,CAAC,IAAI,CAAU;IACtB,OAAO,CAAC,IAAI,CAAuB;IACnC,OAAO,CAAC,MAAM,CAAwB;IACtC,OAAO,CAAC,gBAAgB,CAAuB;IAC/C,OAAO,CAAC,iBAAiB,CAAuB;IAEhD,OAAO,CAAC,KAAK,CAAmB;IAEhC,QAAQ,UAAQ;gBAEJ,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,GAAE,MAAM,GAAG,IAAW,EAAE,KAAK,GAAE,OAAO,GAAG,IAAW;IAQ7F,IAAI,GAAG,IAIM,MAAM,CAFlB;IAED,IAAI,GAAG,CAAC,GAAG,EAAE,MAAM,EAOlB;IAED,IAAI,QAAQ,WAEX;IAED,IAAI,QAAQ,CAAC,CAAC,QAAA,EAEb;IAED,IAAI,GAAG,WAEN;IAED,IAAI,GAAG,CAAC,GAAG,QAAA,EAMV;IAED,IAAI,SAAS,WAEZ;IAED,IAAI,SAAS,CAAC,CAAC,QAAA,EAEd;IAED;;;OAGG;IACH,IAAI,GAAG,kBAEN;IAED,IAAI,GAAG,CAAC,GAAG,eAAA,EAGV;IAED,IAAI,KAAK,YAER;IAED,IAAI,KAAK,CAAC,KAAK,SAAA,EAGd;IAED,IAAI,eAAe,kBAElB;IAED,IAAI,eAAe,CAAC,eAAe,eAAA,EAElC;IAED,IAAI,gBAAgB,kBAEnB;IAED,IAAI,gBAAgB,CAAC,gBAAgB,eAAA,EAEpC;IAED;;OAEG;IACH,KAAK;IAQL,IAAI;IAOJ,MAAM,CAAC,MAAM,CAAC,IAAI,EAAE,WAAW,GAAG,IAAI,EAAE,IAAI,EAAE,WAAW,GAAG,IAAI,EAAE,GAAG,SAAa,EAAE,MAAM,SAAS;IAUnG,MAAM,CAAC,KAAK,EAAE,WAAW;IAIzB,MAAM,CAAC,kBAAkB,CAAC,IAAI,EAAE,WAAW,GAAG,IAAI,EAAE,IAAI,EAAE,WAAW,GAAG,IAAI,EAAE,GAAG,SAAa,EAAE,MAAM,SAAS;IAiB/G,kBAAkB,CAAC,KAAK,EAAE,WAAW,EAAE,GAAG,SAAa,EAAE,MAAM,SAAS;IAIxE;;OAEG;IACH,gBAAgB,CAAC,QAAQ,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,SAAS,GAAE,MAAM,GAAG,IAAW;IAMnF;;;OAGG;IACH,IAAI,CAAC,QAAQ,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,SAAS,GAAE,MAAM,GAAG,IAAW;IA+BvE;;OAEG;IACH,UAAU,CAAC,SAAS,EAAE,WAAW;IA0BjC,MAAM,CAAC,eAAe,CAAC,MAAM,EAAE,WAAW,EAAE,MAAM,EAAE,WAAW;IAI/D,SAAS,CAAC,SAAS,EAAE,WAAW;IAShC,MAAM,CAAC,SAAS,CAAC,MAAM,EAAE,WAAW,EAAE,MAAM,EAAE,WAAW;IAKzD;;;;OAIG;IAEH,IAAI,iBAAiB,IAAI,YAAY,CAIpC;IAED,IAAI,iBAAiB,IAAI,YAAY,CAIpC;IAED;;;OAGG;IACH,IAAI,IAAI,IAAI,SAAS,CAepB;IAGD,MAAM,CAAC,QAAQ,CAAC,IAAI,EAAE,SAAS,GAAG,WAAW;IAoB7C;;;OAGG;IACH,oBAAoB,CAAC,EAAE,EAAE,WAAW,EAAE,EAAE,EAAE,WAAW,GAAG,WAAW,GAAG,IAAI;IAmC1E;;OAEG;IAEH,QAAQ;IAYR,MAAM,IAAI,eAAe;IASzB,MAAM,CAAC,QAAQ,CAAC,IAAI,EAAE,eAAe;IAKrC,gBAAgB,IAAI,yBAAyB;IAW7C,MAAM,CAAC,kBAAkB,CAAC,IAAI,EAAE,yBAAyB;CAO5D;AAED,eAAe,WAAW,CAAC"}
+{"version":3,"file":"Coordinates.d.ts","sourceRoot":"","sources":["../../../src/coordinates/Coordinates.ts"],"names":[],"mappings":"AAAA,OAAO,EAC0C,KAAK,SAAS,EAAE,KAAK,YAAY,EACjF,MAAM,cAAc,CAAC;AAKtB,OAAO,KAAK,EAAE,yBAAyB,EAAE,eAAe,EAAE,OAAO,EAAE,MAAM,UAAU,CAAC;AAgHpF;;;;;;;;;;;;;;;;;;;;;GAqBG;AACH,cAAM,WAAW;IAEb,OAAO,CAAC,IAAI,CAAU;IACtB,OAAO,CAAC,IAAI,CAAU;IACtB,OAAO,CAAC,IAAI,CAAuB;IACnC,OAAO,CAAC,MAAM,CAAwB;IACtC,OAAO,CAAC,gBAAgB,CAAuB;IAC/C,OAAO,CAAC,iBAAiB,CAAuB;IAEhD,OAAO,CAAC,KAAK,CAAmB;IAChC,OAAO,CAAC,WAAW,CAA0B;IAE7C,QAAQ,UAAQ;gBAEJ,GAAG,EAAE,MAAM,EAAE,GAAG,EAAE,MAAM,EAAE,GAAG,GAAE,MAAM,GAAG,IAAW,EAAE,KAAK,GAAE,OAAO,GAAG,IAAW;IAQ7F,IAAI,GAAG,IAIM,MAAM,CAFlB;IAED,IAAI,GAAG,CAAC,GAAG,EAAE,MAAM,EAQlB;IAED,IAAI,QAAQ,WAEX;IAED,IAAI,QAAQ,CAAC,CAAC,QAAA,EAEb;IAED,IAAI,GAAG,WAEN;IAED,IAAI,GAAG,CAAC,GAAG,QAAA,EAOV;IAED,IAAI,SAAS,WAEZ;IAED,IAAI,SAAS,CAAC,CAAC,QAAA,EAEd;IAED;;;;OAIG;IACH,IAAI,GAAG,kBAEN;IAED,IAAI,GAAG,CAAC,GAAG,eAAA,EAGV;IAED,IAAI,KAAK,YAER;IAED,IAAI,KAAK,CAAC,KAAK,SAAA,EAGd;IAED,IAAI,eAAe,kBAElB;IAED,IAAI,eAAe,CAAC,eAAe,eAAA,EAElC;IAED,IAAI,gBAAgB,kBAEnB;IAED,IAAI,gBAAgB,CAAC,gBAAgB,eAAA,EAEpC;IAED;;OAEG;IACH,KAAK;IAQL,IAAI;IAOJ,MAAM,CAAC,MAAM,CAAC,IAAI,EAAE,WAAW,GAAG,IAAI,EAAE,IAAI,EAAE,WAAW,GAAG,IAAI,EAAE,GAAG,SAAa,EAAE,MAAM,SAAS;IAUnG,MAAM,CAAC,KAAK,EAAE,WAAW;IAIzB,MAAM,CAAC,kBAAkB,CAAC,IAAI,EAAE,WAAW,GAAG,IAAI,EAAE,IAAI,EAAE,WAAW,GAAG,IAAI,EAAE,GAAG,SAAa,EAAE,MAAM,SAAS;IAiB/G,kBAAkB,CAAC,KAAK,EAAE,WAAW,EAAE,GAAG,SAAa,EAAE,MAAM,SAAS;IAIxE;;OAEG;IACH,gBAAgB,CAAC,QAAQ,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,SAAS,GAAE,MAAM,GAAG,IAAW;IAMnF;;;;;;;;OAQG;IACH,IAAI,CAAC,QAAQ,EAAE,MAAM,EAAE,OAAO,EAAE,MAAM,EAAE,SAAS,GAAE,MAAM,GAAG,IAAW;IAuCvE;;;;;;;;;;;;OAYG;IACH,UAAU,CAAC,SAAS,EAAE,WAAW;IAYjC,MAAM,CAAC,eAAe,CAAC,MAAM,EAAE,WAAW,EAAE,MAAM,EAAE,WAAW;IAI/D;;;;OAIG;IACH,SAAS,CAAC,SAAS,EAAE,WAAW;IAahC,MAAM,CAAC,SAAS,CAAC,MAAM,EAAE,WAAW,EAAE,MAAM,EAAE,WAAW;IAKzD;;OAEG;IAEH,IAAI,iBAAiB,IAAI,YAAY,CAIpC;IAED,IAAI,iBAAiB,IAAI,YAAY,CAIpC;IAED;;OAEG;IACH,IAAI,IAAI,IAAI,SAAS,CAOpB;IAGD;;;;;OAKG;IACH,OAAO,KAAK,UAAU,GAerB;IAGD,MAAM,CAAC,QAAQ,CAAC,IAAI,EAAE,SAAS,GAAG,WAAW;IAc7C;;;;;;;;OAQG;IACH,oBAAoB,CAAC,EAAE,EAAE,WAAW,EAAE,EAAE,EAAE,WAAW,GAAG,WAAW,GAAG,IAAI;IAmC1E;;OAEG;IAEH,QAAQ;IAYR,MAAM,IAAI,eAAe;IASzB,MAAM,CAAC,QAAQ,CAAC,IAAI,EAAE,eAAe;IAKrC,gBAAgB,IAAI,yBAAyB;IAW7C,MAAM,CAAC,kBAAkB,CAAC,IAAI,EAAE,yBAAyB;CAO5D;AAED,eAAe,WAAW,CAAC"}

package/dist/src/coordinates/GeoRef.d.ts

@@ -3,6 +3,12 @@ import { BuildingLevel, GeoRefJson } from '../types';
 import { default as Coordinates } from './Coordinates';
 declare class GeoRef {
     scale: number;
+    /**
+     * Orientation of the local frame, in radians, clockwise from north, around the local
+     * vertical. Clockwise because that is the convention everything user-facing uses here
+     * (Coordinates#bearingTo, AbsoluteHeading, Rotations.getHeadingFromQuaternion), which is
+     * why the rotations below negate it to work in the counter-clockwise ENU frame.
+     */
     heading: number;
     origin: Coordinates;
     buildingLevels: BuildingLevel[];

package/dist/src/coordinates/GeoRef.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"GeoRef.d.ts","sourceRoot":"","sources":["../../../src/coordinates/GeoRef.ts"],"names":[],"mappings":"AAEA,OAAO,EAAuB,KAAK,SAAS,EAAE,MAAM,cAAc,CAAC;AACnE,OAAO,KAAK,EAAE,aAAa,EAAE,UAAU,EAAE,MAAM,UAAU,CAAC;AAE1D,OAAO,WAAW,MAAM,eAAe,CAAC;AAExC,cAAM,MAAM;IAER,KAAK,SAAK;IACV,OAAO,SAAK;IACZ,MAAM,EAAE,WAAW,CAAC;IACpB,cAAc,EAAE,aAAa,EAAE,CAAC;gBAEpB,MAAM,EAAE,WAAW,EAAE,cAAc,GAAE,aAAa,EAAO;IAKrE;;OAEG;IACH,YAAY,CAAC,aAAa,EAAE,SAAS;IAyFrC;;OAEG;IACH,YAAY,CAAC,MAAM,EAAE,WAAW;IAQhC,MAAM,IAAI,UAAU;IAQpB,MAAM,CAAC,QAAQ,CAAC,IAAI,EAAE,UAAU;CAMnC;AAED,eAAe,MAAM,CAAC"}
+{"version":3,"file":"GeoRef.d.ts","sourceRoot":"","sources":["../../../src/coordinates/GeoRef.ts"],"names":[],"mappings":"AAEA,OAAO,EAAuB,KAAK,SAAS,EAAE,MAAM,cAAc,CAAC;AACnE,OAAO,KAAK,EAAE,aAAa,EAAE,UAAU,EAAE,MAAM,UAAU,CAAC;AAE1D,OAAO,WAAW,MAAM,eAAe,CAAC;AAExC,cAAM,MAAM;IAER,KAAK,SAAK;IAEV;;;;;OAKG;IACH,OAAO,SAAK;IAEZ,MAAM,EAAE,WAAW,CAAC;IACpB,cAAc,EAAE,aAAa,EAAE,CAAC;gBAEpB,MAAM,EAAE,WAAW,EAAE,cAAc,GAAE,aAAa,EAAO;IAKrE;;OAEG;IACH,YAAY,CAAC,aAAa,EAAE,SAAS;IA8FrC;;OAEG;IACH,YAAY,CAAC,MAAM,EAAE,WAAW;IAQhC,MAAM,IAAI,UAAU;IAQpB,MAAM,CAAC,QAAQ,CAAC,IAAI,EAAE,UAAU;CAMnC;AAED,eAAe,MAAM,CAAC"}

package/dist/vite.config.cdn.d.ts

@@ -0,0 +1,3 @@
+declare const _default: import('vite').UserConfig;
+export default _default;
+//# sourceMappingURL=vite.config.cdn.d.ts.map

package/dist/vite.config.cdn.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"vite.config.cdn.d.ts","sourceRoot":"","sources":["../vite.config.cdn.ts"],"names":[],"mappings":";AAMA,wBAAuE"}

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@wemap/geo",
-  "version": "14.6.0",
+  "version": "14.7.0",
   "description": "Geospatial primitives (coordinates, levels, bounding boxes) for the Wemap SDK.",
   "license": "UNLICENSED",
   "author": "Wemap SAS",
@@ -38,8 +38,8 @@
   ],
   "dependencies": {
     "@turf/boolean-point-in-polygon": "^7.3.4",
-    "@wemap/maths": "14.6.0",
-    "@wemap/utils": "14.6.0"
+    "@wemap/maths": "14.7.0",
+    "@wemap/utils": "14.7.0"
   },
   "devDependencies": {
     "@types/geojson": "^7946.0.14",