@math.gl/geospatial 3.6.3 → 4.0.0-alpha.2

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Files changed (48) hide show
  1. package/dist/constants.js.map +1 -0
  2. package/dist/ellipsoid/ellipsoid.d.ts +9 -8
  3. package/dist/ellipsoid/ellipsoid.d.ts.map +1 -1
  4. package/dist/{esm/ellipsoid → ellipsoid}/ellipsoid.js +5 -6
  5. package/dist/ellipsoid/ellipsoid.js.map +1 -0
  6. package/dist/ellipsoid/helpers/ellipsoid-transform.d.ts +4 -4
  7. package/dist/ellipsoid/helpers/ellipsoid-transform.d.ts.map +1 -1
  8. package/dist/{esm/ellipsoid → ellipsoid}/helpers/ellipsoid-transform.js +1 -1
  9. package/dist/ellipsoid/helpers/ellipsoid-transform.js.map +1 -0
  10. package/dist/ellipsoid/helpers/scale-to-geodetic-surface.d.ts +2 -2
  11. package/dist/ellipsoid/helpers/scale-to-geodetic-surface.d.ts.map +1 -1
  12. package/dist/{esm/ellipsoid → ellipsoid}/helpers/scale-to-geodetic-surface.js +1 -1
  13. package/dist/ellipsoid/helpers/scale-to-geodetic-surface.js.map +1 -0
  14. package/dist/index.cjs +407 -0
  15. package/dist/index.d.ts +1 -1
  16. package/dist/index.d.ts.map +1 -1
  17. package/dist/index.js +3 -0
  18. package/dist/index.js.map +1 -0
  19. package/dist/type-utils.d.ts +6 -6
  20. package/dist/type-utils.d.ts.map +1 -1
  21. package/dist/{esm/type-utils.js → type-utils.js} +1 -1
  22. package/dist/type-utils.js.map +1 -0
  23. package/package.json +7 -6
  24. package/src/ellipsoid/ellipsoid.ts +33 -33
  25. package/src/ellipsoid/helpers/ellipsoid-transform.ts +14 -12
  26. package/src/ellipsoid/helpers/scale-to-geodetic-surface.ts +2 -2
  27. package/src/index.ts +1 -1
  28. package/src/type-utils.ts +15 -9
  29. package/dist/es5/constants.js +0 -22
  30. package/dist/es5/constants.js.map +0 -1
  31. package/dist/es5/ellipsoid/ellipsoid.js +0 -210
  32. package/dist/es5/ellipsoid/ellipsoid.js.map +0 -1
  33. package/dist/es5/ellipsoid/helpers/ellipsoid-transform.js +0 -134
  34. package/dist/es5/ellipsoid/helpers/ellipsoid-transform.js.map +0 -1
  35. package/dist/es5/ellipsoid/helpers/scale-to-geodetic-surface.js +0 -74
  36. package/dist/es5/ellipsoid/helpers/scale-to-geodetic-surface.js.map +0 -1
  37. package/dist/es5/index.js +0 -24
  38. package/dist/es5/index.js.map +0 -1
  39. package/dist/es5/type-utils.js +0 -96
  40. package/dist/es5/type-utils.js.map +0 -1
  41. package/dist/esm/constants.js.map +0 -1
  42. package/dist/esm/ellipsoid/ellipsoid.js.map +0 -1
  43. package/dist/esm/ellipsoid/helpers/ellipsoid-transform.js.map +0 -1
  44. package/dist/esm/ellipsoid/helpers/scale-to-geodetic-surface.js.map +0 -1
  45. package/dist/esm/index.js +0 -3
  46. package/dist/esm/index.js.map +0 -1
  47. package/dist/esm/type-utils.js.map +0 -1
  48. /package/dist/{esm/constants.js → constants.js} +0 -0
@@ -1 +0,0 @@
1
- {"version":3,"sources":["../../../src/ellipsoid/ellipsoid.ts"],"names":["scratchVector","Vector3","scratchNormal","scratchK","scratchPosition","scratchHeight","scratchCartesian","wgs84","Ellipsoid","x","y","z","_MathUtils","EPSILON1","radii","radiiSquared","radiiToTheFourth","oneOverRadii","oneOverRadiiSquared","minimumRadius","Math","min","maximumRadius","max","squaredXOverSquaredZ","Object","freeze","right","Boolean","equals","toString","cartographic","result","normal","k","height","geodeticSurfaceNormalCartographic","copy","scale","gamma","sqrt","dot","add","to","cartesian","from","point","scaleToGeodeticSurface","undefined","geodeticSurfaceNormal","h","subtract","longitude","atan2","latitude","asin","sign","vec3","length","origin","Matrix4","firstAxis","secondAxis","thirdAxis","normalize","cartographicVectorRadians","cosLatitude","cos","set","sin","positionX","positionY","positionZ","beta","multiplyScalar","position","buffer","EPSILON15","abs","WGS84_RADIUS_X","WGS84_RADIUS_Y","WGS84_RADIUS_Z"],"mappings":";;;;;;;;;;;;;;;;;;;AAIA;;AAUA;;AAEA;;AACA;;AAEA;;AACA;;;;;;AAEA,IAAMA,aAAa,GAAG,IAAIC,aAAJ,EAAtB;AACA,IAAMC,aAAa,GAAG,IAAID,aAAJ,EAAtB;AACA,IAAME,QAAQ,GAAG,IAAIF,aAAJ,EAAjB;AACA,IAAMG,eAAe,GAAG,IAAIH,aAAJ,EAAxB;AACA,IAAMI,aAAa,GAAG,IAAIJ,aAAJ,EAAtB;AACA,IAAMK,gBAAgB,GAAG,IAAIL,aAAJ,EAAzB;AAEA,IAAIM,KAAJ;;IAOqBC,S;AAkBnB,uBAAuC;AAAA,QAA3BC,CAA2B,uEAAvB,GAAuB;AAAA,QAAlBC,CAAkB,uEAAd,GAAc;AAAA,QAATC,CAAS,uEAAL,GAAK;AAAA;AAAA;AAAA;AAAA;AAAA;AAAA;AAAA;AAAA;AAAA,kEAPGC,iBAAWC,QAOd;AAAA;AACrC,sBAAOJ,CAAC,IAAI,GAAZ;AACA,sBAAOC,CAAC,IAAI,GAAZ;AACA,sBAAOC,CAAC,IAAI,GAAZ;AAEA,SAAKG,KAAL,GAAa,IAAIb,aAAJ,CAAYQ,CAAZ,EAAeC,CAAf,EAAkBC,CAAlB,CAAb;AAEA,SAAKI,YAAL,GAAoB,IAAId,aAAJ,CAAYQ,CAAC,GAAGA,CAAhB,EAAmBC,CAAC,GAAGA,CAAvB,EAA0BC,CAAC,GAAGA,CAA9B,CAApB;AAEA,SAAKK,gBAAL,GAAwB,IAAIf,aAAJ,CAAYQ,CAAC,GAAGA,CAAJ,GAAQA,CAAR,GAAYA,CAAxB,EAA2BC,CAAC,GAAGA,CAAJ,GAAQA,CAAR,GAAYA,CAAvC,EAA0CC,CAAC,GAAGA,CAAJ,GAAQA,CAAR,GAAYA,CAAtD,CAAxB;AAEA,SAAKM,YAAL,GAAoB,IAAIhB,aAAJ,CAClBQ,CAAC,KAAK,GAAN,GAAY,GAAZ,GAAkB,MAAMA,CADN,EAElBC,CAAC,KAAK,GAAN,GAAY,GAAZ,GAAkB,MAAMA,CAFN,EAGlBC,CAAC,KAAK,GAAN,GAAY,GAAZ,GAAkB,MAAMA,CAHN,CAApB;AAMA,SAAKO,mBAAL,GAA2B,IAAIjB,aAAJ,CACzBQ,CAAC,KAAK,GAAN,GAAY,GAAZ,GAAkB,OAAOA,CAAC,GAAGA,CAAX,CADO,EAEzBC,CAAC,KAAK,GAAN,GAAY,GAAZ,GAAkB,OAAOA,CAAC,GAAGA,CAAX,CAFO,EAGzBC,CAAC,KAAK,GAAN,GAAY,GAAZ,GAAkB,OAAOA,CAAC,GAAGA,CAAX,CAHO,CAA3B;AAMA,SAAKQ,aAAL,GAAqBC,IAAI,CAACC,GAAL,CAASZ,CAAT,EAAYC,CAAZ,EAAeC,CAAf,CAArB;AAEA,SAAKW,aAAL,GAAqBF,IAAI,CAACG,GAAL,CAASd,CAAT,EAAYC,CAAZ,EAAeC,CAAf,CAArB;;AAEA,QAAI,KAAKI,YAAL,CAAkBJ,CAAlB,KAAwB,CAA5B,EAA+B;AAC7B,WAAKa,oBAAL,GAA4B,KAAKT,YAAL,CAAkBN,CAAlB,GAAsB,KAAKM,YAAL,CAAkBJ,CAApE;AACD;;AAEDc,IAAAA,MAAM,CAACC,MAAP,CAAc,IAAd;AACD;;;;WAGD,gBAAOC,KAAP,EAAkC;AAChC,aAAO,SAASA,KAAT,IAAkBC,OAAO,CAACD,KAAK,IAAI,KAAKb,KAAL,CAAWe,MAAX,CAAkBF,KAAK,CAACb,KAAxB,CAAV,CAAhC;AACD;;;WAGD,oBAAmB;AACjB,aAAO,KAAKA,KAAL,CAAWgB,QAAX,EAAP;AACD;;;WAMD,iCAAwBC,YAAxB,EAA0D;AAAA,UAApBC,MAAoB,uEAAX,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAAW;AACxD,UAAMC,MAAM,GAAG/B,aAAf;AACA,UAAMgC,CAAC,GAAG/B,QAAV;;AAEA,uDAAqB4B,YAArB;AAAA,UAAWI,MAAX;;AACA,WAAKC,iCAAL,CAAuCL,YAAvC,EAAqDE,MAArD;AACAC,MAAAA,CAAC,CAACG,IAAF,CAAO,KAAKtB,YAAZ,EAA0BuB,KAA1B,CAAgCL,MAAhC;AAEA,UAAMM,KAAK,GAAGnB,IAAI,CAACoB,IAAL,CAAUP,MAAM,CAACQ,GAAP,CAAWP,CAAX,CAAV,CAAd;AACAA,MAAAA,CAAC,CAACI,KAAF,CAAQ,IAAIC,KAAZ;AAEAN,MAAAA,MAAM,CAACK,KAAP,CAAaH,MAAb;AAEAD,MAAAA,CAAC,CAACQ,GAAF,CAAMT,MAAN;AAEA,aAAOC,CAAC,CAACS,EAAF,CAAKX,MAAL,CAAP;AACD;;;WAOD,iCAAwBY,SAAxB,EAAuD;AAAA,UAApBZ,MAAoB,uEAAX,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAAW;AACrD1B,MAAAA,gBAAgB,CAACuC,IAAjB,CAAsBD,SAAtB;AACA,UAAME,KAAK,GAAG,KAAKC,sBAAL,CAA4BzC,gBAA5B,EAA8CF,eAA9C,CAAd;;AAEA,UAAI,CAAC0C,KAAL,EAAY;AACV,eAAOE,SAAP;AACD;;AAED,UAAMf,MAAM,GAAG,KAAKgB,qBAAL,CAA2BH,KAA3B,EAAkC5C,aAAlC,CAAf;AAEA,UAAMgD,CAAC,GAAG7C,aAAV;AACA6C,MAAAA,CAAC,CAACb,IAAF,CAAO/B,gBAAP,EAAyB6C,QAAzB,CAAkCL,KAAlC;AAEA,UAAMM,SAAS,GAAGhC,IAAI,CAACiC,KAAL,CAAWpB,MAAM,CAACvB,CAAlB,EAAqBuB,MAAM,CAACxB,CAA5B,CAAlB;AACA,UAAM6C,QAAQ,GAAGlC,IAAI,CAACmC,IAAL,CAAUtB,MAAM,CAACtB,CAAjB,CAAjB;AACA,UAAMwB,MAAM,GAAGf,IAAI,CAACoC,IAAL,CAAUC,IAAI,CAAChB,GAAL,CAASS,CAAT,EAAY5C,gBAAZ,CAAV,IAA2CmD,IAAI,CAACC,MAAL,CAAYR,CAAZ,CAA1D;AAEA,aAAO,0CAA0B,CAACE,SAAD,EAAYE,QAAZ,EAAsBnB,MAAtB,CAA1B,EAAyDH,MAAzD,CAAP;AACD;;;WAOD,iCAAwB2B,MAAxB,EAAwD;AAAA,UAAxB3B,MAAwB,uEAAf,IAAI4B,aAAJ,EAAe;AACtD,aAAO,iCAAuB,IAAvB,EAA6B,MAA7B,EAAqC,OAArC,EAA8C,IAA9C,EAAoDD,MAApD,EAA4D3B,MAA5D,CAAP;AACD;;;WAsBD,gCAAuB6B,SAAvB,EAAkCC,UAAlC,EAA8CC,SAA9C,EAAyDJ,MAAzD,EAAyF;AAAA,UAAxB3B,MAAwB,uEAAf,IAAI4B,aAAJ,EAAe;AACvF,aAAO,iCAAuB,IAAvB,EAA6BC,SAA7B,EAAwCC,UAAxC,EAAoDC,SAApD,EAA+DJ,MAA/D,EAAuE3B,MAAvE,CAAP;AACD;;;WAMD,iCAAwBY,SAAxB,EAAuD;AAAA,UAApBZ,MAAoB,uEAAX,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAAW;AACrD,aAAOhC,aAAa,CAAC6C,IAAd,CAAmBD,SAAnB,EAA8BoB,SAA9B,GAA0CrB,EAA1C,CAA6CX,MAA7C,CAAP;AACD;;;WAKD,2CAAkCD,YAAlC,EAAoE;AAAA,UAApBC,MAAoB,uEAAX,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAAW;AAClE,UAAMiC,yBAAyB,GAAG,0CAA0BlC,YAA1B,CAAlC;AAEA,UAAMqB,SAAS,GAAGa,yBAAyB,CAAC,CAAD,CAA3C;AACA,UAAMX,QAAQ,GAAGW,yBAAyB,CAAC,CAAD,CAA1C;AAEA,UAAMC,WAAW,GAAG9C,IAAI,CAAC+C,GAAL,CAASb,QAAT,CAApB;AAEAtD,MAAAA,aAAa,CACVoE,GADH,CACOF,WAAW,GAAG9C,IAAI,CAAC+C,GAAL,CAASf,SAAT,CADrB,EAC0Cc,WAAW,GAAG9C,IAAI,CAACiD,GAAL,CAASjB,SAAT,CADxD,EAC6EhC,IAAI,CAACiD,GAAL,CAASf,QAAT,CAD7E,EAEGU,SAFH;AAIA,aAAOhE,aAAa,CAAC2C,EAAd,CAAiBX,MAAjB,CAAP;AACD;;;WAKD,+BAAsBY,SAAtB,EAAqD;AAAA,UAApBZ,MAAoB,uEAAX,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAAW;AACnD,aAAOhC,aAAa,CAAC6C,IAAd,CAAmBD,SAAnB,EAA8BN,KAA9B,CAAoC,KAAKpB,mBAAzC,EAA8D8C,SAA9D,GAA0ErB,EAA1E,CAA6EX,MAA7E,CAAP;AACD;;;WAKD,gCAAuBY,SAAvB,EAA4CZ,MAA5C,EAAyE;AACvE,aAAO,sCAAuBY,SAAvB,EAAkC,IAAlC,EAAwCZ,MAAxC,CAAP;AACD;;;WAID,kCAAyBY,SAAzB,EAAsF;AAAA,UAAxCZ,MAAwC,uEAArB,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAAqB;AACpF5B,MAAAA,eAAe,CAACyC,IAAhB,CAAqBD,SAArB;AAEA,UAAM0B,SAAS,GAAGlE,eAAe,CAACK,CAAlC;AACA,UAAM8D,SAAS,GAAGnE,eAAe,CAACM,CAAlC;AACA,UAAM8D,SAAS,GAAGpE,eAAe,CAACO,CAAlC;AACA,UAAMO,mBAAmB,GAAG,KAAKA,mBAAjC;AAEA,UAAMuD,IAAI,GACR,MACArD,IAAI,CAACoB,IAAL,CACE8B,SAAS,GAAGA,SAAZ,GAAwBpD,mBAAmB,CAACT,CAA5C,GACE8D,SAAS,GAAGA,SAAZ,GAAwBrD,mBAAmB,CAACR,CAD9C,GAEE8D,SAAS,GAAGA,SAAZ,GAAwBtD,mBAAmB,CAACP,CAHhD,CAFF;AAQA,aAAOP,eAAe,CAACsE,cAAhB,CAA+BD,IAA/B,EAAqC9B,EAArC,CAAwCX,MAAxC,CAAP;AACD;;;WAID,wCAA+B2C,QAA/B,EAA2F;AAAA,UAAxC3C,MAAwC,uEAArB,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAAqB;AACzF,aAAO5B,eAAe,CAACyC,IAAhB,CAAqB8B,QAArB,EAA+BrC,KAA/B,CAAqC,KAAKrB,YAA1C,EAAwD0B,EAAxD,CAA2DX,MAA3D,CAAP;AACD;;;WAID,0CAAiC2C,QAAjC,EAA6F;AAAA,UAAxC3C,MAAwC,uEAArB,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAAqB;AAC3F,aAAO5B,eAAe,CAACyC,IAAhB,CAAqB8B,QAArB,EAA+BrC,KAA/B,CAAqC,KAAKxB,KAA1C,EAAiD6B,EAAjD,CAAoDX,MAApD,CAAP;AACD;;;WAGD,+CACE2C,QADF,EAIY;AAAA,UAFVC,MAEU,uEAFO,CAEP;AAAA,UADV5C,MACU,uEADS,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CACT;AAEV,wBAAO,kBAAO,KAAKlB,KAAL,CAAWL,CAAlB,EAAqB,KAAKK,KAAL,CAAWJ,CAAhC,EAAmCE,iBAAWiE,SAA9C,CAAP;AACA,wBAAO,KAAK/D,KAAL,CAAWH,CAAX,GAAe,CAAtB;AAEAP,MAAAA,eAAe,CAACyC,IAAhB,CAAqB8B,QAArB;AACA,UAAMhE,CAAC,GAAGP,eAAe,CAACO,CAAhB,IAAqB,IAAI,KAAKa,oBAA9B,CAAV;;AAEA,UAAIJ,IAAI,CAAC0D,GAAL,CAASnE,CAAT,KAAe,KAAKG,KAAL,CAAWH,CAAX,GAAeiE,MAAlC,EAA0C;AACxC,eAAO5B,SAAP;AACD;;AAED,aAAO5C,eAAe,CAACgE,GAAhB,CAAoB,GAApB,EAAyB,GAAzB,EAA8BzD,CAA9B,EAAiCgC,EAAjC,CAAoCX,MAApC,CAAP;AACD;;;;;;8BAzOkBxB,S,WAEgB,IAAIA,SAAJ,CAAcuE,yBAAd,EAA8BC,yBAA9B,EAA8CC,yBAA9C,C","sourcesContent":["// This file is derived from the Cesium math library under Apache 2 license\n// See LICENSE.md and https://github.com/AnalyticalGraphicsInc/cesium/blob/master/LICENSE.md\n\n/* eslint-disable */\nimport {\n Vector3,\n Matrix4,\n toRadians,\n toDegrees,\n assert,\n equals,\n _MathUtils,\n NumericArray\n} from '@math.gl/core';\nimport * as vec3 from 'gl-matrix/vec3';\n\nimport {WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z} from '../constants';\nimport {fromCartographicToRadians, toCartographicFromRadians} from '../type-utils';\n\nimport scaleToGeodeticSurface from './helpers/scale-to-geodetic-surface';\nimport localFrameToFixedFrame from './helpers/ellipsoid-transform';\n\nconst scratchVector = new Vector3();\nconst scratchNormal = new Vector3();\nconst scratchK = new Vector3();\nconst scratchPosition = new Vector3();\nconst scratchHeight = new Vector3();\nconst scratchCartesian = new Vector3();\n\nlet wgs84;\n\n/**\n * A quadratic surface defined in Cartesian coordinates by the equation\n * `(x / a)^2 + (y / b)^2 + (z / c)^2 = 1`. Primarily used\n * to represent the shape of planetary bodies.\n */\nexport default class Ellipsoid {\n /** An Ellipsoid instance initialized to the WGS84 standard. */\n static readonly WGS84: Ellipsoid = new Ellipsoid(WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z);\n\n readonly radii: Vector3;\n readonly radiiSquared: Vector3;\n readonly radiiToTheFourth: Vector3;\n readonly oneOverRadii: Vector3;\n readonly oneOverRadiiSquared: Vector3;\n readonly minimumRadius: number;\n readonly maximumRadius: number;\n readonly centerToleranceSquared: number = _MathUtils.EPSILON1;\n readonly squaredXOverSquaredZ: number;\n\n /** Creates an Ellipsoid from a Cartesian specifying the radii in x, y, and z directions. */\n constructor(x: number, y: number, z: number);\n constructor();\n\n constructor(x = 0.0, y = 0.0, z = 0.0) {\n assert(x >= 0.0);\n assert(y >= 0.0);\n assert(z >= 0.0);\n\n this.radii = new Vector3(x, y, z);\n\n this.radiiSquared = new Vector3(x * x, y * y, z * z);\n\n this.radiiToTheFourth = new Vector3(x * x * x * x, y * y * y * y, z * z * z * z);\n\n this.oneOverRadii = new Vector3(\n x === 0.0 ? 0.0 : 1.0 / x,\n y === 0.0 ? 0.0 : 1.0 / y,\n z === 0.0 ? 0.0 : 1.0 / z\n );\n\n this.oneOverRadiiSquared = new Vector3(\n x === 0.0 ? 0.0 : 1.0 / (x * x),\n y === 0.0 ? 0.0 : 1.0 / (y * y),\n z === 0.0 ? 0.0 : 1.0 / (z * z)\n );\n\n this.minimumRadius = Math.min(x, y, z);\n\n this.maximumRadius = Math.max(x, y, z);\n\n if (this.radiiSquared.z !== 0) {\n this.squaredXOverSquaredZ = this.radiiSquared.x / this.radiiSquared.z;\n }\n\n Object.freeze(this);\n }\n\n /** Compares this Ellipsoid against the provided Ellipsoid componentwise */\n equals(right: Ellipsoid): boolean {\n return this === right || Boolean(right && this.radii.equals(right.radii));\n }\n\n /** Creates a string representing this Ellipsoid in the format '(radii.x, radii.y, radii.z)'. */\n toString(): string {\n return this.radii.toString();\n }\n\n /** Converts the provided cartographic to Cartesian representation. */\n cartographicToCartesian(cartographic: number[], result: Vector3): Vector3;\n cartographicToCartesian(cartographic: number[], result?: number[]): number[];\n\n cartographicToCartesian(cartographic, result = [0, 0, 0]) {\n const normal = scratchNormal;\n const k = scratchK;\n\n const [, , height] = cartographic;\n this.geodeticSurfaceNormalCartographic(cartographic, normal);\n k.copy(this.radiiSquared).scale(normal);\n\n const gamma = Math.sqrt(normal.dot(k));\n k.scale(1 / gamma);\n\n normal.scale(height);\n\n k.add(normal);\n\n return k.to(result);\n }\n\n /** Converts the provided cartesian to cartographic (lng/lat/z) representation.\n * The cartesian is undefined at the center of the ellipsoid. */\n cartesianToCartographic(cartesian: number[], result: Vector3): Vector3;\n cartesianToCartographic(cartesian: number[], result?: number[]): number[];\n\n cartesianToCartographic(cartesian, result = [0, 0, 0]) {\n scratchCartesian.from(cartesian);\n const point = this.scaleToGeodeticSurface(scratchCartesian, scratchPosition);\n\n if (!point) {\n return undefined;\n }\n\n const normal = this.geodeticSurfaceNormal(point, scratchNormal);\n\n const h = scratchHeight;\n h.copy(scratchCartesian).subtract(point);\n\n const longitude = Math.atan2(normal.y, normal.x);\n const latitude = Math.asin(normal.z);\n const height = Math.sign(vec3.dot(h, scratchCartesian)) * vec3.length(h);\n\n return toCartographicFromRadians([longitude, latitude, height], result);\n }\n\n /** Computes a 4x4 transformation matrix from a reference frame with an east-north-up axes\n * centered at the provided origin to the provided ellipsoid's fixed reference frame. */\n eastNorthUpToFixedFrame(origin: number[], result?: Matrix4): Matrix4;\n eastNorthUpToFixedFrame(origin: number[], result: number[]): number[];\n\n eastNorthUpToFixedFrame(origin, result = new Matrix4()) {\n return localFrameToFixedFrame(this, 'east', 'north', 'up', origin, result);\n }\n\n /** Computes a 4x4 transformation matrix from a reference frame centered at\n * the provided origin to the ellipsoid's fixed reference frame.\n */\n localFrameToFixedFrame(\n firstAxis: string,\n secondAxis: string,\n thirdAxis: string,\n origin: Readonly<NumericArray>,\n result?: Matrix4\n ): Matrix4;\n localFrameToFixedFrame<Matrix4T>(\n firstAxis: string,\n secondAxis: string,\n thirdAxis: string,\n origin: Readonly<NumericArray>,\n result: number[]\n ): number[];\n\n // Computes a 4x4 transformation matrix from a reference frame centered at\n // the provided origin to the ellipsoid's fixed reference frame.\n localFrameToFixedFrame(firstAxis, secondAxis, thirdAxis, origin, result = new Matrix4()) {\n return localFrameToFixedFrame(this, firstAxis, secondAxis, thirdAxis, origin, result);\n }\n\n /** Computes the unit vector directed from the center of this ellipsoid toward\n * the provided Cartesian position. */\n geocentricSurfaceNormal(cartesian: number[], result?: number[]): number[];\n geocentricSurfaceNormal<NumArray>(cartesian: number[], result: NumArray): NumArray;\n geocentricSurfaceNormal(cartesian, result = [0, 0, 0]) {\n return scratchVector.from(cartesian).normalize().to(result);\n }\n\n /** Computes the normal of the plane tangent to the surface of the ellipsoid at provided position. */\n geodeticSurfaceNormalCartographic<NumArray>(cartographic: number[], result: NumArray): NumArray;\n geodeticSurfaceNormalCartographic(cartographic: number[]): number[];\n geodeticSurfaceNormalCartographic(cartographic, result = [0, 0, 0]) {\n const cartographicVectorRadians = fromCartographicToRadians(cartographic);\n\n const longitude = cartographicVectorRadians[0];\n const latitude = cartographicVectorRadians[1];\n\n const cosLatitude = Math.cos(latitude);\n\n scratchVector\n .set(cosLatitude * Math.cos(longitude), cosLatitude * Math.sin(longitude), Math.sin(latitude))\n .normalize();\n\n return scratchVector.to(result);\n }\n\n /** Computes the normal of the plane tangent to the surface of the ellipsoid at the provided position. */\n geodeticSurfaceNormal<NumArrayT>(cartesian: number[], result: NumArrayT): NumArrayT;\n geodeticSurfaceNormal(cartesian: number[]): number[];\n geodeticSurfaceNormal(cartesian, result = [0, 0, 0]) {\n return scratchVector.from(cartesian).scale(this.oneOverRadiiSquared).normalize().to(result);\n }\n\n /** Scales the provided Cartesian position along the geodetic surface normal\n * so that it is on the surface of this ellipsoid. If the position is\n * at the center of the ellipsoid, this function returns undefined. */\n scaleToGeodeticSurface(cartesian: number[], result?: number[]): number[] {\n return scaleToGeodeticSurface(cartesian, this, result);\n }\n\n /** Scales the provided Cartesian position along the geocentric surface normal\n * so that it is on the surface of this ellipsoid. */\n scaleToGeocentricSurface(cartesian: number[], result: number[] = [0, 0, 0]): number[] {\n scratchPosition.from(cartesian);\n\n const positionX = scratchPosition.x;\n const positionY = scratchPosition.y;\n const positionZ = scratchPosition.z;\n const oneOverRadiiSquared = this.oneOverRadiiSquared;\n\n const beta =\n 1.0 /\n Math.sqrt(\n positionX * positionX * oneOverRadiiSquared.x +\n positionY * positionY * oneOverRadiiSquared.y +\n positionZ * positionZ * oneOverRadiiSquared.z\n );\n\n return scratchPosition.multiplyScalar(beta).to(result);\n }\n\n /** Transforms a Cartesian X, Y, Z position to the ellipsoid-scaled space by multiplying\n * its components by the result of `Ellipsoid#oneOverRadii` */\n transformPositionToScaledSpace(position: number[], result: number[] = [0, 0, 0]): number[] {\n return scratchPosition.from(position).scale(this.oneOverRadii).to(result);\n }\n\n /** Transforms a Cartesian X, Y, Z position from the ellipsoid-scaled space by multiplying\n * its components by the result of `Ellipsoid#radii`. */\n transformPositionFromScaledSpace(position: number[], result: number[] = [0, 0, 0]): number[] {\n return scratchPosition.from(position).scale(this.radii).to(result);\n }\n\n /** Computes a point which is the intersection of the surface normal with the z-axis. */\n getSurfaceNormalIntersectionWithZAxis(\n position: number[],\n buffer: number = 0,\n result: number[] = [0, 0, 0]\n ): number[] {\n // Ellipsoid must be an ellipsoid of revolution (radii.x == radii.y)\n assert(equals(this.radii.x, this.radii.y, _MathUtils.EPSILON15));\n assert(this.radii.z > 0);\n\n scratchPosition.from(position);\n const z = scratchPosition.z * (1 - this.squaredXOverSquaredZ);\n\n if (Math.abs(z) >= this.radii.z - buffer) {\n return undefined;\n }\n\n return scratchPosition.set(0.0, 0.0, z).to(result);\n }\n}\n"],"file":"ellipsoid.js"}
@@ -1,134 +0,0 @@
1
- "use strict";
2
-
3
- Object.defineProperty(exports, "__esModule", {
4
- value: true
5
- });
6
- exports.default = localFrameToFixedFrame;
7
-
8
- var _core = require("@math.gl/core");
9
-
10
- var EPSILON14 = 1e-14;
11
- var scratchOrigin = new _core.Vector3();
12
- var VECTOR_PRODUCT_LOCAL_FRAME = {
13
- up: {
14
- south: 'east',
15
- north: 'west',
16
- west: 'south',
17
- east: 'north'
18
- },
19
- down: {
20
- south: 'west',
21
- north: 'east',
22
- west: 'north',
23
- east: 'south'
24
- },
25
- south: {
26
- up: 'west',
27
- down: 'east',
28
- west: 'down',
29
- east: 'up'
30
- },
31
- north: {
32
- up: 'east',
33
- down: 'west',
34
- west: 'up',
35
- east: 'down'
36
- },
37
- west: {
38
- up: 'north',
39
- down: 'south',
40
- north: 'down',
41
- south: 'up'
42
- },
43
- east: {
44
- up: 'south',
45
- down: 'north',
46
- north: 'up',
47
- south: 'down'
48
- }
49
- };
50
- var degeneratePositionLocalFrame = {
51
- north: [-1, 0, 0],
52
- east: [0, 1, 0],
53
- up: [0, 0, 1],
54
- south: [1, 0, 0],
55
- west: [0, -1, 0],
56
- down: [0, 0, -1]
57
- };
58
- var scratchAxisVectors = {
59
- east: new _core.Vector3(),
60
- north: new _core.Vector3(),
61
- up: new _core.Vector3(),
62
- west: new _core.Vector3(),
63
- south: new _core.Vector3(),
64
- down: new _core.Vector3()
65
- };
66
- var scratchVector1 = new _core.Vector3();
67
- var scratchVector2 = new _core.Vector3();
68
- var scratchVector3 = new _core.Vector3();
69
-
70
- function localFrameToFixedFrame(ellipsoid, firstAxis, secondAxis, thirdAxis, cartesianOrigin, result) {
71
- var thirdAxisInferred = VECTOR_PRODUCT_LOCAL_FRAME[firstAxis] && VECTOR_PRODUCT_LOCAL_FRAME[firstAxis][secondAxis];
72
- (0, _core.assert)(thirdAxisInferred && (!thirdAxis || thirdAxis === thirdAxisInferred));
73
- var firstAxisVector;
74
- var secondAxisVector;
75
- var thirdAxisVector;
76
- var origin = scratchOrigin.copy(cartesianOrigin);
77
- var atPole = (0, _core.equals)(origin.x, 0.0, EPSILON14) && (0, _core.equals)(origin.y, 0.0, EPSILON14);
78
-
79
- if (atPole) {
80
- var sign = Math.sign(origin.z);
81
- firstAxisVector = scratchVector1.fromArray(degeneratePositionLocalFrame[firstAxis]);
82
-
83
- if (firstAxis !== 'east' && firstAxis !== 'west') {
84
- firstAxisVector.scale(sign);
85
- }
86
-
87
- secondAxisVector = scratchVector2.fromArray(degeneratePositionLocalFrame[secondAxis]);
88
-
89
- if (secondAxis !== 'east' && secondAxis !== 'west') {
90
- secondAxisVector.scale(sign);
91
- }
92
-
93
- thirdAxisVector = scratchVector3.fromArray(degeneratePositionLocalFrame[thirdAxis]);
94
-
95
- if (thirdAxis !== 'east' && thirdAxis !== 'west') {
96
- thirdAxisVector.scale(sign);
97
- }
98
- } else {
99
- var up = scratchAxisVectors.up,
100
- east = scratchAxisVectors.east,
101
- north = scratchAxisVectors.north;
102
- east.set(-origin.y, origin.x, 0.0).normalize();
103
- ellipsoid.geodeticSurfaceNormal(origin, up);
104
- north.copy(up).cross(east);
105
- var down = scratchAxisVectors.down,
106
- west = scratchAxisVectors.west,
107
- south = scratchAxisVectors.south;
108
- down.copy(up).scale(-1);
109
- west.copy(east).scale(-1);
110
- south.copy(north).scale(-1);
111
- firstAxisVector = scratchAxisVectors[firstAxis];
112
- secondAxisVector = scratchAxisVectors[secondAxis];
113
- thirdAxisVector = scratchAxisVectors[thirdAxis];
114
- }
115
-
116
- result[0] = firstAxisVector.x;
117
- result[1] = firstAxisVector.y;
118
- result[2] = firstAxisVector.z;
119
- result[3] = 0.0;
120
- result[4] = secondAxisVector.x;
121
- result[5] = secondAxisVector.y;
122
- result[6] = secondAxisVector.z;
123
- result[7] = 0.0;
124
- result[8] = thirdAxisVector.x;
125
- result[9] = thirdAxisVector.y;
126
- result[10] = thirdAxisVector.z;
127
- result[11] = 0.0;
128
- result[12] = origin.x;
129
- result[13] = origin.y;
130
- result[14] = origin.z;
131
- result[15] = 1.0;
132
- return result;
133
- }
134
- //# sourceMappingURL=ellipsoid-transform.js.map
@@ -1 +0,0 @@
1
- {"version":3,"sources":["../../../../src/ellipsoid/helpers/ellipsoid-transform.ts"],"names":["EPSILON14","scratchOrigin","Vector3","VECTOR_PRODUCT_LOCAL_FRAME","up","south","north","west","east","down","degeneratePositionLocalFrame","scratchAxisVectors","scratchVector1","scratchVector2","scratchVector3","localFrameToFixedFrame","ellipsoid","firstAxis","secondAxis","thirdAxis","cartesianOrigin","result","thirdAxisInferred","firstAxisVector","secondAxisVector","thirdAxisVector","origin","copy","atPole","x","y","sign","Math","z","fromArray","scale","set","normalize","geodeticSurfaceNormal","cross"],"mappings":";;;;;;;AAAA;;AAIA,IAAMA,SAAS,GAAG,KAAlB;AAEA,IAAMC,aAAa,GAAG,IAAIC,aAAJ,EAAtB;AAGA,IAAMC,0BAA0B,GAAG;AACjCC,EAAAA,EAAE,EAAE;AACFC,IAAAA,KAAK,EAAE,MADL;AAEFC,IAAAA,KAAK,EAAE,MAFL;AAGFC,IAAAA,IAAI,EAAE,OAHJ;AAIFC,IAAAA,IAAI,EAAE;AAJJ,GAD6B;AAOjCC,EAAAA,IAAI,EAAE;AACJJ,IAAAA,KAAK,EAAE,MADH;AAEJC,IAAAA,KAAK,EAAE,MAFH;AAGJC,IAAAA,IAAI,EAAE,OAHF;AAIJC,IAAAA,IAAI,EAAE;AAJF,GAP2B;AAajCH,EAAAA,KAAK,EAAE;AACLD,IAAAA,EAAE,EAAE,MADC;AAELK,IAAAA,IAAI,EAAE,MAFD;AAGLF,IAAAA,IAAI,EAAE,MAHD;AAILC,IAAAA,IAAI,EAAE;AAJD,GAb0B;AAmBjCF,EAAAA,KAAK,EAAE;AACLF,IAAAA,EAAE,EAAE,MADC;AAELK,IAAAA,IAAI,EAAE,MAFD;AAGLF,IAAAA,IAAI,EAAE,IAHD;AAILC,IAAAA,IAAI,EAAE;AAJD,GAnB0B;AAyBjCD,EAAAA,IAAI,EAAE;AACJH,IAAAA,EAAE,EAAE,OADA;AAEJK,IAAAA,IAAI,EAAE,OAFF;AAGJH,IAAAA,KAAK,EAAE,MAHH;AAIJD,IAAAA,KAAK,EAAE;AAJH,GAzB2B;AA+BjCG,EAAAA,IAAI,EAAE;AACJJ,IAAAA,EAAE,EAAE,OADA;AAEJK,IAAAA,IAAI,EAAE,OAFF;AAGJH,IAAAA,KAAK,EAAE,IAHH;AAIJD,IAAAA,KAAK,EAAE;AAJH;AA/B2B,CAAnC;AAuCA,IAAMK,4BAA4B,GAAG;AACnCJ,EAAAA,KAAK,EAAE,CAAC,CAAC,CAAF,EAAK,CAAL,EAAQ,CAAR,CAD4B;AAEnCE,EAAAA,IAAI,EAAE,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAF6B;AAGnCJ,EAAAA,EAAE,EAAE,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAH+B;AAInCC,EAAAA,KAAK,EAAE,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAJ4B;AAKnCE,EAAAA,IAAI,EAAE,CAAC,CAAD,EAAI,CAAC,CAAL,EAAQ,CAAR,CAL6B;AAMnCE,EAAAA,IAAI,EAAE,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAC,CAAR;AAN6B,CAArC;AASA,IAAME,kBAAkB,GAAG;AACzBH,EAAAA,IAAI,EAAE,IAAIN,aAAJ,EADmB;AAEzBI,EAAAA,KAAK,EAAE,IAAIJ,aAAJ,EAFkB;AAGzBE,EAAAA,EAAE,EAAE,IAAIF,aAAJ,EAHqB;AAIzBK,EAAAA,IAAI,EAAE,IAAIL,aAAJ,EAJmB;AAKzBG,EAAAA,KAAK,EAAE,IAAIH,aAAJ,EALkB;AAMzBO,EAAAA,IAAI,EAAE,IAAIP,aAAJ;AANmB,CAA3B;AASA,IAAMU,cAAc,GAAG,IAAIV,aAAJ,EAAvB;AACA,IAAMW,cAAc,GAAG,IAAIX,aAAJ,EAAvB;AACA,IAAMY,cAAc,GAAG,IAAIZ,aAAJ,EAAvB;;AAOe,SAASa,sBAAT,CACbC,SADa,EAEbC,SAFa,EAGbC,UAHa,EAIbC,SAJa,EAKbC,eALa,EAMbC,MANa,EAOH;AACV,MAAMC,iBAAiB,GACrBnB,0BAA0B,CAACc,SAAD,CAA1B,IACCd,0BAA0B,CAACc,SAAD,CAA1B,CAAsCC,UAAtC,CAFH;AAIA,oBAAOI,iBAAiB,KAAK,CAACH,SAAD,IAAcA,SAAS,KAAKG,iBAAjC,CAAxB;AAEA,MAAIC,eAAJ;AACA,MAAIC,gBAAJ;AACA,MAAIC,eAAJ;AAEA,MAAMC,MAAM,GAAGzB,aAAa,CAAC0B,IAAd,CAAmBP,eAAnB,CAAf;AAGA,MAAMQ,MAAM,GAAG,kBAAcF,MAAM,CAACG,CAArB,EAAwB,GAAxB,EAA6B7B,SAA7B,KAA2C,kBAAc0B,MAAM,CAACI,CAArB,EAAwB,GAAxB,EAA6B9B,SAA7B,CAA1D;;AAEA,MAAI4B,MAAJ,EAAY;AAEV,QAAMG,IAAI,GAAGC,IAAI,CAACD,IAAL,CAAUL,MAAM,CAACO,CAAjB,CAAb;AAEAV,IAAAA,eAAe,GAAGX,cAAc,CAACsB,SAAf,CAAyBxB,4BAA4B,CAACO,SAAD,CAArD,CAAlB;;AACA,QAAIA,SAAS,KAAK,MAAd,IAAwBA,SAAS,KAAK,MAA1C,EAAkD;AAChDM,MAAAA,eAAe,CAACY,KAAhB,CAAsBJ,IAAtB;AACD;;AAEDP,IAAAA,gBAAgB,GAAGX,cAAc,CAACqB,SAAf,CAAyBxB,4BAA4B,CAACQ,UAAD,CAArD,CAAnB;;AACA,QAAIA,UAAU,KAAK,MAAf,IAAyBA,UAAU,KAAK,MAA5C,EAAoD;AAClDM,MAAAA,gBAAgB,CAACW,KAAjB,CAAuBJ,IAAvB;AACD;;AAEDN,IAAAA,eAAe,GAAGX,cAAc,CAACoB,SAAf,CAAyBxB,4BAA4B,CAACS,SAAD,CAArD,CAAlB;;AACA,QAAIA,SAAS,KAAK,MAAd,IAAwBA,SAAS,KAAK,MAA1C,EAAkD;AAChDM,MAAAA,eAAe,CAACU,KAAhB,CAAsBJ,IAAtB;AACD;AACF,GAlBD,MAkBO;AAEL,QAAO3B,EAAP,GAA0BO,kBAA1B,CAAOP,EAAP;AAAA,QAAWI,IAAX,GAA0BG,kBAA1B,CAAWH,IAAX;AAAA,QAAiBF,KAAjB,GAA0BK,kBAA1B,CAAiBL,KAAjB;AAEAE,IAAAA,IAAI,CAAC4B,GAAL,CAAS,CAACV,MAAM,CAACI,CAAjB,EAAoBJ,MAAM,CAACG,CAA3B,EAA8B,GAA9B,EAAmCQ,SAAnC;AACArB,IAAAA,SAAS,CAACsB,qBAAV,CAAgCZ,MAAhC,EAAwCtB,EAAxC;AACAE,IAAAA,KAAK,CAACqB,IAAN,CAAWvB,EAAX,EAAemC,KAAf,CAAqB/B,IAArB;AAEA,QAAOC,IAAP,GAA4BE,kBAA5B,CAAOF,IAAP;AAAA,QAAaF,IAAb,GAA4BI,kBAA5B,CAAaJ,IAAb;AAAA,QAAmBF,KAAnB,GAA4BM,kBAA5B,CAAmBN,KAAnB;AAEAI,IAAAA,IAAI,CAACkB,IAAL,CAAUvB,EAAV,EAAc+B,KAAd,CAAoB,CAAC,CAArB;AACA5B,IAAAA,IAAI,CAACoB,IAAL,CAAUnB,IAAV,EAAgB2B,KAAhB,CAAsB,CAAC,CAAvB;AACA9B,IAAAA,KAAK,CAACsB,IAAN,CAAWrB,KAAX,EAAkB6B,KAAlB,CAAwB,CAAC,CAAzB;AAGAZ,IAAAA,eAAe,GAAGZ,kBAAkB,CAACM,SAAD,CAApC;AACAO,IAAAA,gBAAgB,GAAGb,kBAAkB,CAACO,UAAD,CAArC;AACAO,IAAAA,eAAe,GAAGd,kBAAkB,CAACQ,SAAD,CAApC;AACD;;AAGDE,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYE,eAAe,CAACM,CAA5B;AACAR,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYE,eAAe,CAACO,CAA5B;AACAT,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYE,eAAe,CAACU,CAA5B;AACAZ,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAY,GAAZ;AACAA,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYG,gBAAgB,CAACK,CAA7B;AACAR,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYG,gBAAgB,CAACM,CAA7B;AACAT,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYG,gBAAgB,CAACS,CAA7B;AACAZ,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAY,GAAZ;AACAA,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYI,eAAe,CAACI,CAA5B;AACAR,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYI,eAAe,CAACK,CAA5B;AACAT,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAaI,eAAe,CAACQ,CAA7B;AACAZ,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAa,GAAb;AACAA,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAaK,MAAM,CAACG,CAApB;AACAR,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAaK,MAAM,CAACI,CAApB;AACAT,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAaK,MAAM,CAACO,CAApB;AACAZ,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAa,GAAb;AACA,SAAOA,MAAP;AACD","sourcesContent":["import {Vector3, assert, equals as equalsEpsilon} from '@math.gl/core';\n\nimport type Ellipsoid from '../ellipsoid';\n\nconst EPSILON14 = 1e-14;\n\nconst scratchOrigin = new Vector3();\n\n// Caclulate third axis from given two axii\nconst VECTOR_PRODUCT_LOCAL_FRAME = {\n up: {\n south: 'east',\n north: 'west',\n west: 'south',\n east: 'north'\n },\n down: {\n south: 'west',\n north: 'east',\n west: 'north',\n east: 'south'\n },\n south: {\n up: 'west',\n down: 'east',\n west: 'down',\n east: 'up'\n },\n north: {\n up: 'east',\n down: 'west',\n west: 'up',\n east: 'down'\n },\n west: {\n up: 'north',\n down: 'south',\n north: 'down',\n south: 'up'\n },\n east: {\n up: 'south',\n down: 'north',\n north: 'up',\n south: 'down'\n }\n} as const;\n\nconst degeneratePositionLocalFrame = {\n north: [-1, 0, 0],\n east: [0, 1, 0],\n up: [0, 0, 1],\n south: [1, 0, 0],\n west: [0, -1, 0],\n down: [0, 0, -1]\n} as const;\n\nconst scratchAxisVectors = {\n east: new Vector3(),\n north: new Vector3(),\n up: new Vector3(),\n west: new Vector3(),\n south: new Vector3(),\n down: new Vector3()\n};\n\nconst scratchVector1 = new Vector3();\nconst scratchVector2 = new Vector3();\nconst scratchVector3 = new Vector3();\n\ntype Axis = 'up' | 'down' | 'north' | 'east' | 'south' | 'west';\n\n// Computes a 4x4 transformation matrix from a reference frame\n// centered at the provided origin to the provided ellipsoid's fixed reference frame.\n// eslint-disable-next-line max-statements, max-params, complexity\nexport default function localFrameToFixedFrame(\n ellipsoid: Ellipsoid,\n firstAxis: Axis,\n secondAxis: Axis,\n thirdAxis: Axis,\n cartesianOrigin: number[],\n result: number[]\n): number[] {\n const thirdAxisInferred =\n VECTOR_PRODUCT_LOCAL_FRAME[firstAxis] &&\n (VECTOR_PRODUCT_LOCAL_FRAME[firstAxis][secondAxis] as Axis);\n // firstAxis and secondAxis must be east, north, up, west, south or down.');\n assert(thirdAxisInferred && (!thirdAxis || thirdAxis === thirdAxisInferred));\n\n let firstAxisVector: Vector3;\n let secondAxisVector: Vector3;\n let thirdAxisVector: Vector3;\n\n const origin = scratchOrigin.copy(cartesianOrigin);\n\n // If x and y are zero, assume origin is at a pole, which is a special case.\n const atPole = equalsEpsilon(origin.x, 0.0, EPSILON14) && equalsEpsilon(origin.y, 0.0, EPSILON14);\n\n if (atPole) {\n // Look up axis value and adjust\n const sign = Math.sign(origin.z);\n\n firstAxisVector = scratchVector1.fromArray(degeneratePositionLocalFrame[firstAxis]);\n if (firstAxis !== 'east' && firstAxis !== 'west') {\n firstAxisVector.scale(sign);\n }\n\n secondAxisVector = scratchVector2.fromArray(degeneratePositionLocalFrame[secondAxis]);\n if (secondAxis !== 'east' && secondAxis !== 'west') {\n secondAxisVector.scale(sign);\n }\n\n thirdAxisVector = scratchVector3.fromArray(degeneratePositionLocalFrame[thirdAxis]);\n if (thirdAxis !== 'east' && thirdAxis !== 'west') {\n thirdAxisVector.scale(sign);\n }\n } else {\n // Calculate all axis\n const {up, east, north} = scratchAxisVectors;\n\n east.set(-origin.y, origin.x, 0.0).normalize();\n ellipsoid.geodeticSurfaceNormal(origin, up);\n north.copy(up).cross(east);\n\n const {down, west, south} = scratchAxisVectors;\n\n down.copy(up).scale(-1);\n west.copy(east).scale(-1);\n south.copy(north).scale(-1);\n\n // Pick three axis based on desired orientation\n firstAxisVector = scratchAxisVectors[firstAxis];\n secondAxisVector = scratchAxisVectors[secondAxis];\n thirdAxisVector = scratchAxisVectors[thirdAxis];\n }\n\n // TODO - assuming the result is column-major\n result[0] = firstAxisVector.x;\n result[1] = firstAxisVector.y;\n result[2] = firstAxisVector.z;\n result[3] = 0.0;\n result[4] = secondAxisVector.x;\n result[5] = secondAxisVector.y;\n result[6] = secondAxisVector.z;\n result[7] = 0.0;\n result[8] = thirdAxisVector.x;\n result[9] = thirdAxisVector.y;\n result[10] = thirdAxisVector.z;\n result[11] = 0.0;\n result[12] = origin.x;\n result[13] = origin.y;\n result[14] = origin.z;\n result[15] = 1.0;\n return result;\n}\n"],"file":"ellipsoid-transform.js"}
@@ -1,74 +0,0 @@
1
- "use strict";
2
-
3
- Object.defineProperty(exports, "__esModule", {
4
- value: true
5
- });
6
- exports.default = scaleToGeodeticSurface;
7
-
8
- var _core = require("@math.gl/core");
9
-
10
- var scratchVector = new _core.Vector3();
11
- var scaleToGeodeticSurfaceIntersection = new _core.Vector3();
12
- var scaleToGeodeticSurfaceGradient = new _core.Vector3();
13
-
14
- function scaleToGeodeticSurface(cartesian, ellipsoid) {
15
- var result = arguments.length > 2 && arguments[2] !== undefined ? arguments[2] : [];
16
- var oneOverRadii = ellipsoid.oneOverRadii,
17
- oneOverRadiiSquared = ellipsoid.oneOverRadiiSquared,
18
- centerToleranceSquared = ellipsoid.centerToleranceSquared;
19
- scratchVector.from(cartesian);
20
- var positionX = scratchVector.x;
21
- var positionY = scratchVector.y;
22
- var positionZ = scratchVector.z;
23
- var oneOverRadiiX = oneOverRadii.x;
24
- var oneOverRadiiY = oneOverRadii.y;
25
- var oneOverRadiiZ = oneOverRadii.z;
26
- var x2 = positionX * positionX * oneOverRadiiX * oneOverRadiiX;
27
- var y2 = positionY * positionY * oneOverRadiiY * oneOverRadiiY;
28
- var z2 = positionZ * positionZ * oneOverRadiiZ * oneOverRadiiZ;
29
- var squaredNorm = x2 + y2 + z2;
30
- var ratio = Math.sqrt(1.0 / squaredNorm);
31
-
32
- if (!Number.isFinite(ratio)) {
33
- return undefined;
34
- }
35
-
36
- var intersection = scaleToGeodeticSurfaceIntersection;
37
- intersection.copy(cartesian).scale(ratio);
38
-
39
- if (squaredNorm < centerToleranceSquared) {
40
- return intersection.to(result);
41
- }
42
-
43
- var oneOverRadiiSquaredX = oneOverRadiiSquared.x;
44
- var oneOverRadiiSquaredY = oneOverRadiiSquared.y;
45
- var oneOverRadiiSquaredZ = oneOverRadiiSquared.z;
46
- var gradient = scaleToGeodeticSurfaceGradient;
47
- gradient.set(intersection.x * oneOverRadiiSquaredX * 2.0, intersection.y * oneOverRadiiSquaredY * 2.0, intersection.z * oneOverRadiiSquaredZ * 2.0);
48
- var lambda = (1.0 - ratio) * scratchVector.len() / (0.5 * gradient.len());
49
- var correction = 0.0;
50
- var xMultiplier;
51
- var yMultiplier;
52
- var zMultiplier;
53
- var func;
54
-
55
- do {
56
- lambda -= correction;
57
- xMultiplier = 1.0 / (1.0 + lambda * oneOverRadiiSquaredX);
58
- yMultiplier = 1.0 / (1.0 + lambda * oneOverRadiiSquaredY);
59
- zMultiplier = 1.0 / (1.0 + lambda * oneOverRadiiSquaredZ);
60
- var xMultiplier2 = xMultiplier * xMultiplier;
61
- var yMultiplier2 = yMultiplier * yMultiplier;
62
- var zMultiplier2 = zMultiplier * zMultiplier;
63
- var xMultiplier3 = xMultiplier2 * xMultiplier;
64
- var yMultiplier3 = yMultiplier2 * yMultiplier;
65
- var zMultiplier3 = zMultiplier2 * zMultiplier;
66
- func = x2 * xMultiplier2 + y2 * yMultiplier2 + z2 * zMultiplier2 - 1.0;
67
- var denominator = x2 * xMultiplier3 * oneOverRadiiSquaredX + y2 * yMultiplier3 * oneOverRadiiSquaredY + z2 * zMultiplier3 * oneOverRadiiSquaredZ;
68
- var derivative = -2.0 * denominator;
69
- correction = func / derivative;
70
- } while (Math.abs(func) > _core._MathUtils.EPSILON12);
71
-
72
- return scratchVector.scale([xMultiplier, yMultiplier, zMultiplier]).to(result);
73
- }
74
- //# sourceMappingURL=scale-to-geodetic-surface.js.map
@@ -1 +0,0 @@
1
- {"version":3,"sources":["../../../../src/ellipsoid/helpers/scale-to-geodetic-surface.ts"],"names":["scratchVector","Vector3","scaleToGeodeticSurfaceIntersection","scaleToGeodeticSurfaceGradient","scaleToGeodeticSurface","cartesian","ellipsoid","result","oneOverRadii","oneOverRadiiSquared","centerToleranceSquared","from","positionX","x","positionY","y","positionZ","z","oneOverRadiiX","oneOverRadiiY","oneOverRadiiZ","x2","y2","z2","squaredNorm","ratio","Math","sqrt","Number","isFinite","undefined","intersection","copy","scale","to","oneOverRadiiSquaredX","oneOverRadiiSquaredY","oneOverRadiiSquaredZ","gradient","set","lambda","len","correction","xMultiplier","yMultiplier","zMultiplier","func","xMultiplier2","yMultiplier2","zMultiplier2","xMultiplier3","yMultiplier3","zMultiplier3","denominator","derivative","abs","_MathUtils","EPSILON12"],"mappings":";;;;;;;AACA;;AAGA,IAAMA,aAAa,GAAG,IAAIC,aAAJ,EAAtB;AACA,IAAMC,kCAAkC,GAAG,IAAID,aAAJ,EAA3C;AACA,IAAME,8BAA8B,GAAG,IAAIF,aAAJ,EAAvC;;AAKe,SAASG,sBAAT,CACbC,SADa,EAEbC,SAFa,EAIH;AAAA,MADVC,MACU,uEADS,EACT;AACV,MAAOC,YAAP,GAAoEF,SAApE,CAAOE,YAAP;AAAA,MAAqBC,mBAArB,GAAoEH,SAApE,CAAqBG,mBAArB;AAAA,MAA0CC,sBAA1C,GAAoEJ,SAApE,CAA0CI,sBAA1C;AAEAV,EAAAA,aAAa,CAACW,IAAd,CAAmBN,SAAnB;AAEA,MAAMO,SAAS,GAAGZ,aAAa,CAACa,CAAhC;AACA,MAAMC,SAAS,GAAGd,aAAa,CAACe,CAAhC;AACA,MAAMC,SAAS,GAAGhB,aAAa,CAACiB,CAAhC;AAEA,MAAMC,aAAa,GAAGV,YAAY,CAACK,CAAnC;AACA,MAAMM,aAAa,GAAGX,YAAY,CAACO,CAAnC;AACA,MAAMK,aAAa,GAAGZ,YAAY,CAACS,CAAnC;AAEA,MAAMI,EAAE,GAAGT,SAAS,GAAGA,SAAZ,GAAwBM,aAAxB,GAAwCA,aAAnD;AACA,MAAMI,EAAE,GAAGR,SAAS,GAAGA,SAAZ,GAAwBK,aAAxB,GAAwCA,aAAnD;AACA,MAAMI,EAAE,GAAGP,SAAS,GAAGA,SAAZ,GAAwBI,aAAxB,GAAwCA,aAAnD;AAGA,MAAMI,WAAW,GAAGH,EAAE,GAAGC,EAAL,GAAUC,EAA9B;AACA,MAAME,KAAK,GAAGC,IAAI,CAACC,IAAL,CAAU,MAAMH,WAAhB,CAAd;;AAGA,MAAI,CAACI,MAAM,CAACC,QAAP,CAAgBJ,KAAhB,CAAL,EAA6B;AAC3B,WAAOK,SAAP;AACD;;AAGD,MAAMC,YAAY,GAAG7B,kCAArB;AACA6B,EAAAA,YAAY,CAACC,IAAb,CAAkB3B,SAAlB,EAA6B4B,KAA7B,CAAmCR,KAAnC;;AAGA,MAAID,WAAW,GAAGd,sBAAlB,EAA0C;AACxC,WAAOqB,YAAY,CAACG,EAAb,CAAgB3B,MAAhB,CAAP;AACD;;AAED,MAAM4B,oBAAoB,GAAG1B,mBAAmB,CAACI,CAAjD;AACA,MAAMuB,oBAAoB,GAAG3B,mBAAmB,CAACM,CAAjD;AACA,MAAMsB,oBAAoB,GAAG5B,mBAAmB,CAACQ,CAAjD;AAIA,MAAMqB,QAAQ,GAAGnC,8BAAjB;AACAmC,EAAAA,QAAQ,CAACC,GAAT,CACER,YAAY,CAAClB,CAAb,GAAiBsB,oBAAjB,GAAwC,GAD1C,EAEEJ,YAAY,CAAChB,CAAb,GAAiBqB,oBAAjB,GAAwC,GAF1C,EAGEL,YAAY,CAACd,CAAb,GAAiBoB,oBAAjB,GAAwC,GAH1C;AAOA,MAAIG,MAAM,GAAI,CAAC,MAAMf,KAAP,IAAgBzB,aAAa,CAACyC,GAAd,EAAjB,IAAyC,MAAMH,QAAQ,CAACG,GAAT,EAA/C,CAAb;AACA,MAAIC,UAAU,GAAG,GAAjB;AAEA,MAAIC,WAAJ;AACA,MAAIC,WAAJ;AACA,MAAIC,WAAJ;AACA,MAAIC,IAAJ;;AAEA,KAAG;AACDN,IAAAA,MAAM,IAAIE,UAAV;AAEAC,IAAAA,WAAW,GAAG,OAAO,MAAMH,MAAM,GAAGL,oBAAtB,CAAd;AACAS,IAAAA,WAAW,GAAG,OAAO,MAAMJ,MAAM,GAAGJ,oBAAtB,CAAd;AACAS,IAAAA,WAAW,GAAG,OAAO,MAAML,MAAM,GAAGH,oBAAtB,CAAd;AAEA,QAAMU,YAAY,GAAGJ,WAAW,GAAGA,WAAnC;AACA,QAAMK,YAAY,GAAGJ,WAAW,GAAGA,WAAnC;AACA,QAAMK,YAAY,GAAGJ,WAAW,GAAGA,WAAnC;AAEA,QAAMK,YAAY,GAAGH,YAAY,GAAGJ,WAApC;AACA,QAAMQ,YAAY,GAAGH,YAAY,GAAGJ,WAApC;AACA,QAAMQ,YAAY,GAAGH,YAAY,GAAGJ,WAApC;AAEAC,IAAAA,IAAI,GAAGzB,EAAE,GAAG0B,YAAL,GAAoBzB,EAAE,GAAG0B,YAAzB,GAAwCzB,EAAE,GAAG0B,YAA7C,GAA4D,GAAnE;AAIA,QAAMI,WAAW,GACfhC,EAAE,GAAG6B,YAAL,GAAoBf,oBAApB,GACAb,EAAE,GAAG6B,YAAL,GAAoBf,oBADpB,GAEAb,EAAE,GAAG6B,YAAL,GAAoBf,oBAHtB;AAKA,QAAMiB,UAAU,GAAG,CAAC,GAAD,GAAOD,WAA1B;AAEAX,IAAAA,UAAU,GAAGI,IAAI,GAAGQ,UAApB;AACD,GA3BD,QA2BS5B,IAAI,CAAC6B,GAAL,CAAST,IAAT,IAAiBU,iBAAWC,SA3BrC;;AA6BA,SAAOzD,aAAa,CAACiC,KAAd,CAAoB,CAACU,WAAD,EAAcC,WAAd,EAA2BC,WAA3B,CAApB,EAA6DX,EAA7D,CAAgE3B,MAAhE,CAAP;AACD","sourcesContent":["/* eslint-disable */\nimport {Vector3, _MathUtils} from '@math.gl/core';\nimport type Ellipsoid from '../ellipsoid';\n\nconst scratchVector = new Vector3();\nconst scaleToGeodeticSurfaceIntersection = new Vector3();\nconst scaleToGeodeticSurfaceGradient = new Vector3();\n\n// Scales the provided Cartesian position along the geodetic surface normal\n// so that it is on the surface of this ellipsoid. If the position is\n// at the center of the ellipsoid, this function returns undefined.\nexport default function scaleToGeodeticSurface(\n cartesian: number[],\n ellipsoid: Ellipsoid,\n result: number[] = []\n): number[] {\n const {oneOverRadii, oneOverRadiiSquared, centerToleranceSquared} = ellipsoid;\n\n scratchVector.from(cartesian);\n\n const positionX = scratchVector.x;\n const positionY = scratchVector.y;\n const positionZ = scratchVector.z;\n\n const oneOverRadiiX = oneOverRadii.x;\n const oneOverRadiiY = oneOverRadii.y;\n const oneOverRadiiZ = oneOverRadii.z;\n\n const x2 = positionX * positionX * oneOverRadiiX * oneOverRadiiX;\n const y2 = positionY * positionY * oneOverRadiiY * oneOverRadiiY;\n const z2 = positionZ * positionZ * oneOverRadiiZ * oneOverRadiiZ;\n\n // Compute the squared ellipsoid norm.\n const squaredNorm = x2 + y2 + z2;\n const ratio = Math.sqrt(1.0 / squaredNorm);\n\n // When very close to center or at center\n if (!Number.isFinite(ratio)) {\n return undefined;\n }\n\n // As an initial approximation, assume that the radial intersection is the projection point.\n const intersection = scaleToGeodeticSurfaceIntersection;\n intersection.copy(cartesian).scale(ratio);\n\n // If the position is near the center, the iteration will not converge.\n if (squaredNorm < centerToleranceSquared) {\n return intersection.to(result);\n }\n\n const oneOverRadiiSquaredX = oneOverRadiiSquared.x;\n const oneOverRadiiSquaredY = oneOverRadiiSquared.y;\n const oneOverRadiiSquaredZ = oneOverRadiiSquared.z;\n\n // Use the gradient at the intersection point in place of the true unit normal.\n // The difference in magnitude will be absorbed in the multiplier.\n const gradient = scaleToGeodeticSurfaceGradient;\n gradient.set(\n intersection.x * oneOverRadiiSquaredX * 2.0,\n intersection.y * oneOverRadiiSquaredY * 2.0,\n intersection.z * oneOverRadiiSquaredZ * 2.0\n );\n\n // Compute the initial guess at the normal vector multiplier, lambda.\n let lambda = ((1.0 - ratio) * scratchVector.len()) / (0.5 * gradient.len());\n let correction = 0.0;\n\n let xMultiplier;\n let yMultiplier;\n let zMultiplier;\n let func;\n\n do {\n lambda -= correction;\n\n xMultiplier = 1.0 / (1.0 + lambda * oneOverRadiiSquaredX);\n yMultiplier = 1.0 / (1.0 + lambda * oneOverRadiiSquaredY);\n zMultiplier = 1.0 / (1.0 + lambda * oneOverRadiiSquaredZ);\n\n const xMultiplier2 = xMultiplier * xMultiplier;\n const yMultiplier2 = yMultiplier * yMultiplier;\n const zMultiplier2 = zMultiplier * zMultiplier;\n\n const xMultiplier3 = xMultiplier2 * xMultiplier;\n const yMultiplier3 = yMultiplier2 * yMultiplier;\n const zMultiplier3 = zMultiplier2 * zMultiplier;\n\n func = x2 * xMultiplier2 + y2 * yMultiplier2 + z2 * zMultiplier2 - 1.0;\n\n // \"denominator\" here refers to the use of this expression in the velocity and acceleration\n // computations in the sections to follow.\n const denominator =\n x2 * xMultiplier3 * oneOverRadiiSquaredX +\n y2 * yMultiplier3 * oneOverRadiiSquaredY +\n z2 * zMultiplier3 * oneOverRadiiSquaredZ;\n\n const derivative = -2.0 * denominator;\n\n correction = func / derivative;\n } while (Math.abs(func) > _MathUtils.EPSILON12);\n\n return scratchVector.scale([xMultiplier, yMultiplier, zMultiplier]).to(result);\n}\n"],"file":"scale-to-geodetic-surface.js"}
package/dist/es5/index.js DELETED
@@ -1,24 +0,0 @@
1
- "use strict";
2
-
3
- var _interopRequireDefault = require("@babel/runtime/helpers/interopRequireDefault");
4
-
5
- Object.defineProperty(exports, "__esModule", {
6
- value: true
7
- });
8
- Object.defineProperty(exports, "Ellipsoid", {
9
- enumerable: true,
10
- get: function get() {
11
- return _ellipsoid.default;
12
- }
13
- });
14
- Object.defineProperty(exports, "isWGS84", {
15
- enumerable: true,
16
- get: function get() {
17
- return _typeUtils.isWGS84;
18
- }
19
- });
20
-
21
- var _ellipsoid = _interopRequireDefault(require("./ellipsoid/ellipsoid"));
22
-
23
- var _typeUtils = require("./type-utils");
24
- //# sourceMappingURL=index.js.map
@@ -1 +0,0 @@
1
- {"version":3,"sources":["../../src/index.ts"],"names":[],"mappings":";;;;;;;;;;;;;;;;;;;;AAAA;;AACA","sourcesContent":["export {default as Ellipsoid} from './ellipsoid/ellipsoid';\nexport {isWGS84} from './type-utils';\n"],"file":"index.js"}
@@ -1,96 +0,0 @@
1
- "use strict";
2
-
3
- Object.defineProperty(exports, "__esModule", {
4
- value: true
5
- });
6
- exports.fromCartographic = fromCartographic;
7
- exports.fromCartographicToRadians = fromCartographicToRadians;
8
- exports.fromCartographicToDegrees = fromCartographicToDegrees;
9
- exports.toCartographic = toCartographic;
10
- exports.toCartographicFromRadians = toCartographicFromRadians;
11
- exports.toCartographicFromDegrees = toCartographicFromDegrees;
12
- exports.isWGS84 = isWGS84;
13
-
14
- var _core = require("@math.gl/core");
15
-
16
- var _constants = require("./constants");
17
-
18
- function identity(x) {
19
- return x;
20
- }
21
-
22
- var scratchVector = new _core.Vector3();
23
-
24
- function fromCartographic(cartographic) {
25
- var result = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : [];
26
- var map = arguments.length > 2 && arguments[2] !== undefined ? arguments[2] : identity;
27
-
28
- if ('longitude' in cartographic) {
29
- result[0] = map(cartographic.longitude);
30
- result[1] = map(cartographic.latitude);
31
- result[2] = cartographic.height;
32
- } else if ('x' in cartographic) {
33
- result[0] = map(cartographic.x);
34
- result[1] = map(cartographic.y);
35
- result[2] = cartographic.z;
36
- } else {
37
- result[0] = map(cartographic[0]);
38
- result[1] = map(cartographic[1]);
39
- result[2] = cartographic[2];
40
- }
41
-
42
- return result;
43
- }
44
-
45
- function fromCartographicToRadians(cartographic) {
46
- var vector = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : [];
47
- return fromCartographic(cartographic, vector, _core.config._cartographicRadians ? identity : _core.toRadians);
48
- }
49
-
50
- function fromCartographicToDegrees(cartographic) {
51
- var vector = arguments.length > 1 && arguments[1] !== undefined ? arguments[1] : [];
52
- return fromCartographic(cartographic, vector, _core.config._cartographicRadians ? _core.toDegrees : identity);
53
- }
54
-
55
- function toCartographic(vector, cartographic) {
56
- var map = arguments.length > 2 && arguments[2] !== undefined ? arguments[2] : identity;
57
-
58
- if ('longitude' in cartographic) {
59
- cartographic.longitude = map(vector[0]);
60
- cartographic.latitude = map(vector[1]);
61
- cartographic.height = vector[2];
62
- } else if ('x' in cartographic) {
63
- cartographic.x = map(vector[0]);
64
- cartographic.y = map(vector[1]);
65
- cartographic.z = vector[2];
66
- } else {
67
- cartographic[0] = map(vector[0]);
68
- cartographic[1] = map(vector[1]);
69
- cartographic[2] = vector[2];
70
- }
71
-
72
- return cartographic;
73
- }
74
-
75
- function toCartographicFromRadians(vector, cartographic) {
76
- return toCartographic(vector, cartographic, _core.config._cartographicRadians ? identity : _core.toDegrees);
77
- }
78
-
79
- function toCartographicFromDegrees(vector, cartographic) {
80
- return toCartographic(vector, cartographic, _core.config._cartographicRadians ? _core.toRadians : identity);
81
- }
82
-
83
- function isWGS84(vector) {
84
- if (!vector) {
85
- return false;
86
- }
87
-
88
- scratchVector.from(vector);
89
- var oneOverRadiiSquared = _constants.WGS84_CONSTANTS.oneOverRadiiSquared,
90
- centerToleranceSquared = _constants.WGS84_CONSTANTS.centerToleranceSquared;
91
- var x2 = vector[0] * vector[0] * oneOverRadiiSquared[0];
92
- var y2 = vector[1] * vector[1] * oneOverRadiiSquared[1];
93
- var z2 = vector[2] * vector[2] * oneOverRadiiSquared[2];
94
- return Math.abs(x2 + y2 + z2 - 1) < centerToleranceSquared;
95
- }
96
- //# sourceMappingURL=type-utils.js.map
@@ -1 +0,0 @@
1
- {"version":3,"sources":["../../src/type-utils.ts"],"names":["identity","x","scratchVector","Vector3","fromCartographic","cartographic","result","map","longitude","latitude","height","y","z","fromCartographicToRadians","vector","config","_cartographicRadians","toRadians","fromCartographicToDegrees","toDegrees","toCartographic","toCartographicFromRadians","toCartographicFromDegrees","isWGS84","from","oneOverRadiiSquared","WGS84_CONSTANTS","centerToleranceSquared","x2","y2","z2","Math","abs"],"mappings":";;;;;;;;;;;;;AAIA;;AACA;;AAgBA,SAASA,QAAT,CAAkBC,CAAlB,EAAqC;AACnC,SAAOA,CAAP;AACD;;AAED,IAAMC,aAAa,GAAG,IAAIC,aAAJ,EAAtB;;AAQO,SAASC,gBAAT,CACLC,YADK,EAIK;AAAA,MAFVC,MAEU,uEAFD,EAEC;AAAA,MADVC,GACU,uEADJP,QACI;;AACV,MAAI,eAAeK,YAAnB,EAAiC;AAC/BC,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYC,GAAG,CAACF,YAAY,CAACG,SAAd,CAAf;AACAF,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYC,GAAG,CAACF,YAAY,CAACI,QAAd,CAAf;AACAH,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYD,YAAY,CAACK,MAAzB;AACD,GAJD,MAIO,IAAI,OAAOL,YAAX,EAAyB;AAC9BC,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYC,GAAG,CAACF,YAAY,CAACJ,CAAd,CAAf;AACAK,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYC,GAAG,CAACF,YAAY,CAACM,CAAd,CAAf;AACAL,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYD,YAAY,CAACO,CAAzB;AACD,GAJM,MAIA;AACLN,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYC,GAAG,CAACF,YAAY,CAAC,CAAD,CAAb,CAAf;AACAC,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYC,GAAG,CAACF,YAAY,CAAC,CAAD,CAAb,CAAf;AACAC,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYD,YAAY,CAAC,CAAD,CAAxB;AACD;;AACD,SAAOC,MAAP;AACD;;AAOM,SAASO,yBAAT,CACLR,YADK,EAGK;AAAA,MADVS,MACU,uEADD,EACC;AACV,SAAOV,gBAAgB,CAACC,YAAD,EAAeS,MAAf,EAAuBC,aAAOC,oBAAP,GAA8BhB,QAA9B,GAAyCiB,eAAhE,CAAvB;AACD;;AAOM,SAASC,yBAAT,CACLb,YADK,EAGK;AAAA,MADVS,MACU,uEADD,EACC;AACV,SAAOV,gBAAgB,CAACC,YAAD,EAAeS,MAAf,EAAuBC,aAAOC,oBAAP,GAA8BG,eAA9B,GAA0CnB,QAAjE,CAAvB;AACD;;AAEM,SAASoB,cAAT,CACLN,MADK,EAELT,YAFK,EAIF;AAAA,MADHE,GACG,uEAD0BP,QAC1B;;AACH,MAAI,eAAeK,YAAnB,EAAiC;AAC/BA,IAAAA,YAAY,CAACG,SAAb,GAAyBD,GAAG,CAACO,MAAM,CAAC,CAAD,CAAP,CAA5B;AACAT,IAAAA,YAAY,CAACI,QAAb,GAAwBF,GAAG,CAACO,MAAM,CAAC,CAAD,CAAP,CAA3B;AACAT,IAAAA,YAAY,CAACK,MAAb,GAAsBI,MAAM,CAAC,CAAD,CAA5B;AACD,GAJD,MAIO,IAAI,OAAOT,YAAX,EAAyB;AAC9BA,IAAAA,YAAY,CAACJ,CAAb,GAAiBM,GAAG,CAACO,MAAM,CAAC,CAAD,CAAP,CAApB;AACAT,IAAAA,YAAY,CAACM,CAAb,GAAiBJ,GAAG,CAACO,MAAM,CAAC,CAAD,CAAP,CAApB;AACAT,IAAAA,YAAY,CAACO,CAAb,GAAiBE,MAAM,CAAC,CAAD,CAAvB;AACD,GAJM,MAIA;AACLT,IAAAA,YAAY,CAAC,CAAD,CAAZ,GAAkBE,GAAG,CAACO,MAAM,CAAC,CAAD,CAAP,CAArB;AACAT,IAAAA,YAAY,CAAC,CAAD,CAAZ,GAAkBE,GAAG,CAACO,MAAM,CAAC,CAAD,CAAP,CAArB;AACAT,IAAAA,YAAY,CAAC,CAAD,CAAZ,GAAkBS,MAAM,CAAC,CAAD,CAAxB;AACD;;AACD,SAAOT,YAAP;AACD;;AAEM,SAASgB,yBAAT,CACLP,MADK,EAELT,YAFK,EAGF;AACH,SAAOe,cAAc,CAACN,MAAD,EAAST,YAAT,EAAuBU,aAAOC,oBAAP,GAA8BhB,QAA9B,GAAyCmB,eAAhE,CAArB;AACD;;AAEM,SAASG,yBAAT,CACLR,MADK,EAELT,YAFK,EAGF;AACH,SAAOe,cAAc,CAACN,MAAD,EAAST,YAAT,EAAuBU,aAAOC,oBAAP,GAA8BC,eAA9B,GAA0CjB,QAAjE,CAArB;AACD;;AAGM,SAASuB,OAAT,CAAiBT,MAAjB,EAA0D;AAC/D,MAAI,CAACA,MAAL,EAAa;AACX,WAAO,KAAP;AACD;;AACDZ,EAAAA,aAAa,CAACsB,IAAd,CAAmBV,MAAnB;AACA,MAAOW,mBAAP,GAAsDC,0BAAtD,CAAOD,mBAAP;AAAA,MAA4BE,sBAA5B,GAAsDD,0BAAtD,CAA4BC,sBAA5B;AACA,MAAMC,EAAE,GAAGd,MAAM,CAAC,CAAD,CAAN,GAAYA,MAAM,CAAC,CAAD,CAAlB,GAAwBW,mBAAmB,CAAC,CAAD,CAAtD;AACA,MAAMI,EAAE,GAAGf,MAAM,CAAC,CAAD,CAAN,GAAYA,MAAM,CAAC,CAAD,CAAlB,GAAwBW,mBAAmB,CAAC,CAAD,CAAtD;AACA,MAAMK,EAAE,GAAGhB,MAAM,CAAC,CAAD,CAAN,GAAYA,MAAM,CAAC,CAAD,CAAlB,GAAwBW,mBAAmB,CAAC,CAAD,CAAtD;AACA,SAAOM,IAAI,CAACC,GAAL,CAASJ,EAAE,GAAGC,EAAL,GAAUC,EAAV,GAAe,CAAxB,IAA6BH,sBAApC;AACD","sourcesContent":["// This file is derived from the Cesium math library under Apache 2 license\n// See LICENSE.md and https://github.com/AnalyticalGraphicsInc/cesium/blob/master/LICENSE.md\n\nimport type {NumericArray} from '@math.gl/core';\nimport {Vector3, toRadians, toDegrees, config} from '@math.gl/core';\nimport {WGS84_CONSTANTS} from './constants';\n\ntype LngLatHeightObject = {\n longitude: number;\n latitude: number;\n height: number;\n};\n\ntype XYZObject = {\n x: number;\n y: number;\n z: number;\n};\n\ntype Cartographic = LngLatHeightObject | XYZObject | NumericArray;\n\nfunction identity(x: number): number {\n return x;\n}\n\nconst scratchVector = new Vector3();\n\nexport function fromCartographic(cartographic: Cartographic): number[];\nexport function fromCartographic<NumArrayT>(\n cartographic: Cartographic,\n result: NumArrayT,\n map?: (x: number) => number\n): NumArrayT;\nexport function fromCartographic(\n cartographic: Cartographic,\n result = [] as number[],\n map = identity\n): number[] {\n if ('longitude' in cartographic) {\n result[0] = map(cartographic.longitude);\n result[1] = map(cartographic.latitude);\n result[2] = cartographic.height;\n } else if ('x' in cartographic) {\n result[0] = map(cartographic.x);\n result[1] = map(cartographic.y);\n result[2] = cartographic.z;\n } else {\n result[0] = map(cartographic[0]);\n result[1] = map(cartographic[1]);\n result[2] = cartographic[2];\n }\n return result;\n}\n\nexport function fromCartographicToRadians(cartographic: Cartographic, result?: number[]): number[];\nexport function fromCartographicToRadians<TArray>(\n cartographic: Cartographic,\n result: TArray\n): TArray;\nexport function fromCartographicToRadians(\n cartographic: Cartographic,\n vector = [] as number[]\n): number[] {\n return fromCartographic(cartographic, vector, config._cartographicRadians ? identity : toRadians);\n}\n\nexport function fromCartographicToDegrees(cartographic: Cartographic, result?: number[]): number[];\nexport function fromCartographicToDegrees<TArray>(\n cartographic: Cartographic,\n result: TArray\n): TArray;\nexport function fromCartographicToDegrees(\n cartographic: Cartographic,\n vector = [] as number[]\n): number[] {\n return fromCartographic(cartographic, vector, config._cartographicRadians ? toDegrees : identity);\n}\n\nexport function toCartographic<T extends Cartographic>(\n vector: Readonly<NumericArray>,\n cartographic: T,\n map: (x: number) => number = identity\n): T {\n if ('longitude' in cartographic) {\n cartographic.longitude = map(vector[0]);\n cartographic.latitude = map(vector[1]);\n cartographic.height = vector[2];\n } else if ('x' in cartographic) {\n cartographic.x = map(vector[0]);\n cartographic.y = map(vector[1]);\n cartographic.z = vector[2];\n } else {\n cartographic[0] = map(vector[0]);\n cartographic[1] = map(vector[1]);\n cartographic[2] = vector[2];\n }\n return cartographic;\n}\n\nexport function toCartographicFromRadians<T extends Cartographic>(\n vector: Readonly<NumericArray>,\n cartographic: T\n): T {\n return toCartographic(vector, cartographic, config._cartographicRadians ? identity : toDegrees);\n}\n\nexport function toCartographicFromDegrees<T extends Cartographic>(\n vector: Readonly<NumericArray>,\n cartographic: T\n): T {\n return toCartographic(vector, cartographic, config._cartographicRadians ? toRadians : identity);\n}\n\n// Estimates if a vector is close to the surface of the WGS84 Ellipsoid\nexport function isWGS84(vector: Readonly<NumericArray>): boolean {\n if (!vector) {\n return false;\n }\n scratchVector.from(vector);\n const {oneOverRadiiSquared, centerToleranceSquared} = WGS84_CONSTANTS;\n const x2 = vector[0] * vector[0] * oneOverRadiiSquared[0];\n const y2 = vector[1] * vector[1] * oneOverRadiiSquared[1];\n const z2 = vector[2] * vector[2] * oneOverRadiiSquared[2];\n return Math.abs(x2 + y2 + z2 - 1) < centerToleranceSquared;\n}\n\n/*\n\nexport function fromCartographic(cartographic: Cartographic, result?: number[]): number[];\nexport function fromCartographic(cartographic: Cartographic, result: TypedArray): TypedArray;\nexport function fromCartographicToRadians(cartographic: Cartographic, result?: number[]): number[];\nexport function fromCartographicToRadians(\n cartographic: Cartographic,\n result: TypedArray\n): TypedArray;\nexport function fromCartographicToDegrees(cartographic: Cartographic, result?: number[]): number[];\nexport function fromCartographicToDegrees(\n cartographic: Cartographic,\n result: TypedArray\n): TypedArray;\n\nexport function toCartographic(vector: number[] | TypedArray, result: Cartographic): number[];\nexport function toCartographicFromRadians(\n vector: number[] | TypedArray,\n result: Cartographic\n): number[];\nexport function toCartographicFromDegrees(\n vector: number[] | TypedArray,\n result: Cartographic\n): number[];\n\n// Estimates if a vector is close to the surface of the WGS84 Ellipsoid\nexport function isWGS84(vector: number[] | TypedArray): boolean;\n*/\n"],"file":"type-utils.js"}
@@ -1 +0,0 @@
1
- {"version":3,"sources":["../../src/constants.ts"],"names":["WGS84_RADIUS_X","WGS84_RADIUS_Y","WGS84_RADIUS_Z","WGS84_CONSTANTS","radii","radiiSquared","oneOverRadii","oneOverRadiiSquared","maximumRadius","Math","max","centerToleranceSquared"],"mappings":"AAGA,OAAO,MAAMA,cAAc,GAAG,SAAvB;AACP,OAAO,MAAMC,cAAc,GAAG,SAAvB;AACP,OAAO,MAAMC,cAAc,GAAG,kBAAvB;AAIP,OAAO,MAAMC,eAAe,GAAG;AAC7BC,EAAAA,KAAK,EAAE,CAACJ,cAAD,EAAiBC,cAAjB,EAAiCC,cAAjC,CADsB;AAE7BG,EAAAA,YAAY,EAAE,CACZL,cAAc,GAAGA,cADL,EAEZC,cAAc,GAAGA,cAFL,EAGZC,cAAc,GAAGA,cAHL,CAFe;AAO7BI,EAAAA,YAAY,EAAE,CAAC,MAAMN,cAAP,EAAuB,MAAMC,cAA7B,EAA6C,MAAMC,cAAnD,CAPe;AAQ7BK,EAAAA,mBAAmB,EAAE,CACnB,OAAOP,cAAc,GAAGA,cAAxB,CADmB,EAEnB,OAAOC,cAAc,GAAGA,cAAxB,CAFmB,EAGnB,OAAOC,cAAc,GAAGA,cAAxB,CAHmB,CARQ;AAa7BM,EAAAA,aAAa,EAAEC,IAAI,CAACC,GAAL,CAASV,cAAT,EAAyBC,cAAzB,EAAyCC,cAAzC,CAbc;AAc7BS,EAAAA,sBAAsB,EAAE;AAdK,CAAxB","sourcesContent":["// This file is derived from the Cesium math library under Apache 2 license\n// See LICENSE.md and https://github.com/AnalyticalGraphicsInc/cesium/blob/master/LICENSE.md\n\nexport const WGS84_RADIUS_X = 6378137.0;\nexport const WGS84_RADIUS_Y = 6378137.0;\nexport const WGS84_RADIUS_Z = 6356752.3142451793;\n\n// Pre-calculated ellipsoid defaults to avoid utils depending on `ellipsoid.js`\n\nexport const WGS84_CONSTANTS = {\n radii: [WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z],\n radiiSquared: [\n WGS84_RADIUS_X * WGS84_RADIUS_X,\n WGS84_RADIUS_Y * WGS84_RADIUS_Y,\n WGS84_RADIUS_Z * WGS84_RADIUS_Z\n ],\n oneOverRadii: [1.0 / WGS84_RADIUS_X, 1.0 / WGS84_RADIUS_Y, 1.0 / WGS84_RADIUS_Z],\n oneOverRadiiSquared: [\n 1.0 / (WGS84_RADIUS_X * WGS84_RADIUS_X),\n 1.0 / (WGS84_RADIUS_Y * WGS84_RADIUS_Y),\n 1.0 / (WGS84_RADIUS_Z * WGS84_RADIUS_Z)\n ],\n maximumRadius: Math.max(WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z),\n centerToleranceSquared: 1e-1 // EPSILON1;\n};\n"],"file":"constants.js"}
@@ -1 +0,0 @@
1
- {"version":3,"sources":["../../../src/ellipsoid/ellipsoid.ts"],"names":["Vector3","Matrix4","assert","equals","_MathUtils","vec3","WGS84_RADIUS_X","WGS84_RADIUS_Y","WGS84_RADIUS_Z","fromCartographicToRadians","toCartographicFromRadians","scaleToGeodeticSurface","localFrameToFixedFrame","scratchVector","scratchNormal","scratchK","scratchPosition","scratchHeight","scratchCartesian","wgs84","Ellipsoid","constructor","x","y","z","EPSILON1","radii","radiiSquared","radiiToTheFourth","oneOverRadii","oneOverRadiiSquared","minimumRadius","Math","min","maximumRadius","max","squaredXOverSquaredZ","Object","freeze","right","Boolean","toString","cartographicToCartesian","cartographic","result","normal","k","height","geodeticSurfaceNormalCartographic","copy","scale","gamma","sqrt","dot","add","to","cartesianToCartographic","cartesian","from","point","undefined","geodeticSurfaceNormal","h","subtract","longitude","atan2","latitude","asin","sign","length","eastNorthUpToFixedFrame","origin","firstAxis","secondAxis","thirdAxis","geocentricSurfaceNormal","normalize","cartographicVectorRadians","cosLatitude","cos","set","sin","scaleToGeocentricSurface","positionX","positionY","positionZ","beta","multiplyScalar","transformPositionToScaledSpace","position","transformPositionFromScaledSpace","getSurfaceNormalIntersectionWithZAxis","buffer","EPSILON15","abs"],"mappings":";AAIA,SACEA,OADF,EAEEC,OAFF,EAKEC,MALF,EAMEC,MANF,EAOEC,UAPF,QASO,eATP;AAUA,OAAO,KAAKC,IAAZ,MAAsB,gBAAtB;AAEA,SAAQC,cAAR,EAAwBC,cAAxB,EAAwCC,cAAxC,QAA6D,cAA7D;AACA,SAAQC,yBAAR,EAAmCC,yBAAnC,QAAmE,eAAnE;AAEA,OAAOC,sBAAP,MAAmC,qCAAnC;AACA,OAAOC,sBAAP,MAAmC,+BAAnC;AAEA,MAAMC,aAAa,GAAG,IAAIb,OAAJ,EAAtB;AACA,MAAMc,aAAa,GAAG,IAAId,OAAJ,EAAtB;AACA,MAAMe,QAAQ,GAAG,IAAIf,OAAJ,EAAjB;AACA,MAAMgB,eAAe,GAAG,IAAIhB,OAAJ,EAAxB;AACA,MAAMiB,aAAa,GAAG,IAAIjB,OAAJ,EAAtB;AACA,MAAMkB,gBAAgB,GAAG,IAAIlB,OAAJ,EAAzB;AAEA,IAAImB,KAAJ;AAOA,eAAe,MAAMC,SAAN,CAAgB;AAkB7BC,EAAAA,WAAW,CAACC,CAAC,GAAG,GAAL,EAAUC,CAAC,GAAG,GAAd,EAAmBC,CAAC,GAAG,GAAvB,EAA4B;AAAA;;AAAA;;AAAA;;AAAA;;AAAA;;AAAA;;AAAA;;AAAA,oDAPGpB,UAAU,CAACqB,QAOd;;AAAA;;AACrCvB,IAAAA,MAAM,CAACoB,CAAC,IAAI,GAAN,CAAN;AACApB,IAAAA,MAAM,CAACqB,CAAC,IAAI,GAAN,CAAN;AACArB,IAAAA,MAAM,CAACsB,CAAC,IAAI,GAAN,CAAN;AAEA,SAAKE,KAAL,GAAa,IAAI1B,OAAJ,CAAYsB,CAAZ,EAAeC,CAAf,EAAkBC,CAAlB,CAAb;AAEA,SAAKG,YAAL,GAAoB,IAAI3B,OAAJ,CAAYsB,CAAC,GAAGA,CAAhB,EAAmBC,CAAC,GAAGA,CAAvB,EAA0BC,CAAC,GAAGA,CAA9B,CAApB;AAEA,SAAKI,gBAAL,GAAwB,IAAI5B,OAAJ,CAAYsB,CAAC,GAAGA,CAAJ,GAAQA,CAAR,GAAYA,CAAxB,EAA2BC,CAAC,GAAGA,CAAJ,GAAQA,CAAR,GAAYA,CAAvC,EAA0CC,CAAC,GAAGA,CAAJ,GAAQA,CAAR,GAAYA,CAAtD,CAAxB;AAEA,SAAKK,YAAL,GAAoB,IAAI7B,OAAJ,CAClBsB,CAAC,KAAK,GAAN,GAAY,GAAZ,GAAkB,MAAMA,CADN,EAElBC,CAAC,KAAK,GAAN,GAAY,GAAZ,GAAkB,MAAMA,CAFN,EAGlBC,CAAC,KAAK,GAAN,GAAY,GAAZ,GAAkB,MAAMA,CAHN,CAApB;AAMA,SAAKM,mBAAL,GAA2B,IAAI9B,OAAJ,CACzBsB,CAAC,KAAK,GAAN,GAAY,GAAZ,GAAkB,OAAOA,CAAC,GAAGA,CAAX,CADO,EAEzBC,CAAC,KAAK,GAAN,GAAY,GAAZ,GAAkB,OAAOA,CAAC,GAAGA,CAAX,CAFO,EAGzBC,CAAC,KAAK,GAAN,GAAY,GAAZ,GAAkB,OAAOA,CAAC,GAAGA,CAAX,CAHO,CAA3B;AAMA,SAAKO,aAAL,GAAqBC,IAAI,CAACC,GAAL,CAASX,CAAT,EAAYC,CAAZ,EAAeC,CAAf,CAArB;AAEA,SAAKU,aAAL,GAAqBF,IAAI,CAACG,GAAL,CAASb,CAAT,EAAYC,CAAZ,EAAeC,CAAf,CAArB;;AAEA,QAAI,KAAKG,YAAL,CAAkBH,CAAlB,KAAwB,CAA5B,EAA+B;AAC7B,WAAKY,oBAAL,GAA4B,KAAKT,YAAL,CAAkBL,CAAlB,GAAsB,KAAKK,YAAL,CAAkBH,CAApE;AACD;;AAEDa,IAAAA,MAAM,CAACC,MAAP,CAAc,IAAd;AACD;;AAGDnC,EAAAA,MAAM,CAACoC,KAAD,EAA4B;AAChC,WAAO,SAASA,KAAT,IAAkBC,OAAO,CAACD,KAAK,IAAI,KAAKb,KAAL,CAAWvB,MAAX,CAAkBoC,KAAK,CAACb,KAAxB,CAAV,CAAhC;AACD;;AAGDe,EAAAA,QAAQ,GAAW;AACjB,WAAO,KAAKf,KAAL,CAAWe,QAAX,EAAP;AACD;;AAMDC,EAAAA,uBAAuB,CAACC,YAAD,EAAeC,MAAM,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAAxB,EAAmC;AACxD,UAAMC,MAAM,GAAG/B,aAAf;AACA,UAAMgC,CAAC,GAAG/B,QAAV;AAEA,UAAM,IAAKgC,MAAL,IAAeJ,YAArB;AACA,SAAKK,iCAAL,CAAuCL,YAAvC,EAAqDE,MAArD;AACAC,IAAAA,CAAC,CAACG,IAAF,CAAO,KAAKtB,YAAZ,EAA0BuB,KAA1B,CAAgCL,MAAhC;AAEA,UAAMM,KAAK,GAAGnB,IAAI,CAACoB,IAAL,CAAUP,MAAM,CAACQ,GAAP,CAAWP,CAAX,CAAV,CAAd;AACAA,IAAAA,CAAC,CAACI,KAAF,CAAQ,IAAIC,KAAZ;AAEAN,IAAAA,MAAM,CAACK,KAAP,CAAaH,MAAb;AAEAD,IAAAA,CAAC,CAACQ,GAAF,CAAMT,MAAN;AAEA,WAAOC,CAAC,CAACS,EAAF,CAAKX,MAAL,CAAP;AACD;;AAODY,EAAAA,uBAAuB,CAACC,SAAD,EAAYb,MAAM,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAArB,EAAgC;AACrD1B,IAAAA,gBAAgB,CAACwC,IAAjB,CAAsBD,SAAtB;AACA,UAAME,KAAK,GAAG,KAAKhD,sBAAL,CAA4BO,gBAA5B,EAA8CF,eAA9C,CAAd;;AAEA,QAAI,CAAC2C,KAAL,EAAY;AACV,aAAOC,SAAP;AACD;;AAED,UAAMf,MAAM,GAAG,KAAKgB,qBAAL,CAA2BF,KAA3B,EAAkC7C,aAAlC,CAAf;AAEA,UAAMgD,CAAC,GAAG7C,aAAV;AACA6C,IAAAA,CAAC,CAACb,IAAF,CAAO/B,gBAAP,EAAyB6C,QAAzB,CAAkCJ,KAAlC;AAEA,UAAMK,SAAS,GAAGhC,IAAI,CAACiC,KAAL,CAAWpB,MAAM,CAACtB,CAAlB,EAAqBsB,MAAM,CAACvB,CAA5B,CAAlB;AACA,UAAM4C,QAAQ,GAAGlC,IAAI,CAACmC,IAAL,CAAUtB,MAAM,CAACrB,CAAjB,CAAjB;AACA,UAAMuB,MAAM,GAAGf,IAAI,CAACoC,IAAL,CAAU/D,IAAI,CAACgD,GAAL,CAASS,CAAT,EAAY5C,gBAAZ,CAAV,IAA2Cb,IAAI,CAACgE,MAAL,CAAYP,CAAZ,CAA1D;AAEA,WAAOpD,yBAAyB,CAAC,CAACsD,SAAD,EAAYE,QAAZ,EAAsBnB,MAAtB,CAAD,EAAgCH,MAAhC,CAAhC;AACD;;AAOD0B,EAAAA,uBAAuB,CAACC,MAAD,EAAS3B,MAAM,GAAG,IAAI3C,OAAJ,EAAlB,EAAiC;AACtD,WAAOW,sBAAsB,CAAC,IAAD,EAAO,MAAP,EAAe,OAAf,EAAwB,IAAxB,EAA8B2D,MAA9B,EAAsC3B,MAAtC,CAA7B;AACD;;AAsBDhC,EAAAA,sBAAsB,CAAC4D,SAAD,EAAYC,UAAZ,EAAwBC,SAAxB,EAAmCH,MAAnC,EAA2C3B,MAAM,GAAG,IAAI3C,OAAJ,EAApD,EAAmE;AACvF,WAAOW,sBAAsB,CAAC,IAAD,EAAO4D,SAAP,EAAkBC,UAAlB,EAA8BC,SAA9B,EAAyCH,MAAzC,EAAiD3B,MAAjD,CAA7B;AACD;;AAMD+B,EAAAA,uBAAuB,CAAClB,SAAD,EAAYb,MAAM,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAArB,EAAgC;AACrD,WAAO/B,aAAa,CAAC6C,IAAd,CAAmBD,SAAnB,EAA8BmB,SAA9B,GAA0CrB,EAA1C,CAA6CX,MAA7C,CAAP;AACD;;AAKDI,EAAAA,iCAAiC,CAACL,YAAD,EAAeC,MAAM,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAAxB,EAAmC;AAClE,UAAMiC,yBAAyB,GAAGpE,yBAAyB,CAACkC,YAAD,CAA3D;AAEA,UAAMqB,SAAS,GAAGa,yBAAyB,CAAC,CAAD,CAA3C;AACA,UAAMX,QAAQ,GAAGW,yBAAyB,CAAC,CAAD,CAA1C;AAEA,UAAMC,WAAW,GAAG9C,IAAI,CAAC+C,GAAL,CAASb,QAAT,CAApB;AAEArD,IAAAA,aAAa,CACVmE,GADH,CACOF,WAAW,GAAG9C,IAAI,CAAC+C,GAAL,CAASf,SAAT,CADrB,EAC0Cc,WAAW,GAAG9C,IAAI,CAACiD,GAAL,CAASjB,SAAT,CADxD,EAC6EhC,IAAI,CAACiD,GAAL,CAASf,QAAT,CAD7E,EAEGU,SAFH;AAIA,WAAO/D,aAAa,CAAC0C,EAAd,CAAiBX,MAAjB,CAAP;AACD;;AAKDiB,EAAAA,qBAAqB,CAACJ,SAAD,EAAYb,MAAM,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAArB,EAAgC;AACnD,WAAO/B,aAAa,CAAC6C,IAAd,CAAmBD,SAAnB,EAA8BP,KAA9B,CAAoC,KAAKpB,mBAAzC,EAA8D8C,SAA9D,GAA0ErB,EAA1E,CAA6EX,MAA7E,CAAP;AACD;;AAKDjC,EAAAA,sBAAsB,CAAC8C,SAAD,EAAsBb,MAAtB,EAAmD;AACvE,WAAOjC,sBAAsB,CAAC8C,SAAD,EAAY,IAAZ,EAAkBb,MAAlB,CAA7B;AACD;;AAIDsC,EAAAA,wBAAwB,CAACzB,SAAD,EAAsBb,MAAgB,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAAzC,EAA8D;AACpF5B,IAAAA,eAAe,CAAC0C,IAAhB,CAAqBD,SAArB;AAEA,UAAM0B,SAAS,GAAGnE,eAAe,CAACM,CAAlC;AACA,UAAM8D,SAAS,GAAGpE,eAAe,CAACO,CAAlC;AACA,UAAM8D,SAAS,GAAGrE,eAAe,CAACQ,CAAlC;AACA,UAAMM,mBAAmB,GAAG,KAAKA,mBAAjC;AAEA,UAAMwD,IAAI,GACR,MACAtD,IAAI,CAACoB,IAAL,CACE+B,SAAS,GAAGA,SAAZ,GAAwBrD,mBAAmB,CAACR,CAA5C,GACE8D,SAAS,GAAGA,SAAZ,GAAwBtD,mBAAmB,CAACP,CAD9C,GAEE8D,SAAS,GAAGA,SAAZ,GAAwBvD,mBAAmB,CAACN,CAHhD,CAFF;AAQA,WAAOR,eAAe,CAACuE,cAAhB,CAA+BD,IAA/B,EAAqC/B,EAArC,CAAwCX,MAAxC,CAAP;AACD;;AAID4C,EAAAA,8BAA8B,CAACC,QAAD,EAAqB7C,MAAgB,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAAxC,EAA6D;AACzF,WAAO5B,eAAe,CAAC0C,IAAhB,CAAqB+B,QAArB,EAA+BvC,KAA/B,CAAqC,KAAKrB,YAA1C,EAAwD0B,EAAxD,CAA2DX,MAA3D,CAAP;AACD;;AAID8C,EAAAA,gCAAgC,CAACD,QAAD,EAAqB7C,MAAgB,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAAxC,EAA6D;AAC3F,WAAO5B,eAAe,CAAC0C,IAAhB,CAAqB+B,QAArB,EAA+BvC,KAA/B,CAAqC,KAAKxB,KAA1C,EAAiD6B,EAAjD,CAAoDX,MAApD,CAAP;AACD;;AAGD+C,EAAAA,qCAAqC,CACnCF,QADmC,EAEnCG,MAAc,GAAG,CAFkB,EAGnChD,MAAgB,GAAG,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAHgB,EAIzB;AAEV1C,IAAAA,MAAM,CAACC,MAAM,CAAC,KAAKuB,KAAL,CAAWJ,CAAZ,EAAe,KAAKI,KAAL,CAAWH,CAA1B,EAA6BnB,UAAU,CAACyF,SAAxC,CAAP,CAAN;AACA3F,IAAAA,MAAM,CAAC,KAAKwB,KAAL,CAAWF,CAAX,GAAe,CAAhB,CAAN;AAEAR,IAAAA,eAAe,CAAC0C,IAAhB,CAAqB+B,QAArB;AACA,UAAMjE,CAAC,GAAGR,eAAe,CAACQ,CAAhB,IAAqB,IAAI,KAAKY,oBAA9B,CAAV;;AAEA,QAAIJ,IAAI,CAAC8D,GAAL,CAAStE,CAAT,KAAe,KAAKE,KAAL,CAAWF,CAAX,GAAeoE,MAAlC,EAA0C;AACxC,aAAOhC,SAAP;AACD;;AAED,WAAO5C,eAAe,CAACgE,GAAhB,CAAoB,GAApB,EAAyB,GAAzB,EAA8BxD,CAA9B,EAAiC+B,EAAjC,CAAoCX,MAApC,CAAP;AACD;;AAzO4B;;gBAAVxB,S,WAEgB,IAAIA,SAAJ,CAAcd,cAAd,EAA8BC,cAA9B,EAA8CC,cAA9C,C","sourcesContent":["// This file is derived from the Cesium math library under Apache 2 license\n// See LICENSE.md and https://github.com/AnalyticalGraphicsInc/cesium/blob/master/LICENSE.md\n\n/* eslint-disable */\nimport {\n Vector3,\n Matrix4,\n toRadians,\n toDegrees,\n assert,\n equals,\n _MathUtils,\n NumericArray\n} from '@math.gl/core';\nimport * as vec3 from 'gl-matrix/vec3';\n\nimport {WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z} from '../constants';\nimport {fromCartographicToRadians, toCartographicFromRadians} from '../type-utils';\n\nimport scaleToGeodeticSurface from './helpers/scale-to-geodetic-surface';\nimport localFrameToFixedFrame from './helpers/ellipsoid-transform';\n\nconst scratchVector = new Vector3();\nconst scratchNormal = new Vector3();\nconst scratchK = new Vector3();\nconst scratchPosition = new Vector3();\nconst scratchHeight = new Vector3();\nconst scratchCartesian = new Vector3();\n\nlet wgs84;\n\n/**\n * A quadratic surface defined in Cartesian coordinates by the equation\n * `(x / a)^2 + (y / b)^2 + (z / c)^2 = 1`. Primarily used\n * to represent the shape of planetary bodies.\n */\nexport default class Ellipsoid {\n /** An Ellipsoid instance initialized to the WGS84 standard. */\n static readonly WGS84: Ellipsoid = new Ellipsoid(WGS84_RADIUS_X, WGS84_RADIUS_Y, WGS84_RADIUS_Z);\n\n readonly radii: Vector3;\n readonly radiiSquared: Vector3;\n readonly radiiToTheFourth: Vector3;\n readonly oneOverRadii: Vector3;\n readonly oneOverRadiiSquared: Vector3;\n readonly minimumRadius: number;\n readonly maximumRadius: number;\n readonly centerToleranceSquared: number = _MathUtils.EPSILON1;\n readonly squaredXOverSquaredZ: number;\n\n /** Creates an Ellipsoid from a Cartesian specifying the radii in x, y, and z directions. */\n constructor(x: number, y: number, z: number);\n constructor();\n\n constructor(x = 0.0, y = 0.0, z = 0.0) {\n assert(x >= 0.0);\n assert(y >= 0.0);\n assert(z >= 0.0);\n\n this.radii = new Vector3(x, y, z);\n\n this.radiiSquared = new Vector3(x * x, y * y, z * z);\n\n this.radiiToTheFourth = new Vector3(x * x * x * x, y * y * y * y, z * z * z * z);\n\n this.oneOverRadii = new Vector3(\n x === 0.0 ? 0.0 : 1.0 / x,\n y === 0.0 ? 0.0 : 1.0 / y,\n z === 0.0 ? 0.0 : 1.0 / z\n );\n\n this.oneOverRadiiSquared = new Vector3(\n x === 0.0 ? 0.0 : 1.0 / (x * x),\n y === 0.0 ? 0.0 : 1.0 / (y * y),\n z === 0.0 ? 0.0 : 1.0 / (z * z)\n );\n\n this.minimumRadius = Math.min(x, y, z);\n\n this.maximumRadius = Math.max(x, y, z);\n\n if (this.radiiSquared.z !== 0) {\n this.squaredXOverSquaredZ = this.radiiSquared.x / this.radiiSquared.z;\n }\n\n Object.freeze(this);\n }\n\n /** Compares this Ellipsoid against the provided Ellipsoid componentwise */\n equals(right: Ellipsoid): boolean {\n return this === right || Boolean(right && this.radii.equals(right.radii));\n }\n\n /** Creates a string representing this Ellipsoid in the format '(radii.x, radii.y, radii.z)'. */\n toString(): string {\n return this.radii.toString();\n }\n\n /** Converts the provided cartographic to Cartesian representation. */\n cartographicToCartesian(cartographic: number[], result: Vector3): Vector3;\n cartographicToCartesian(cartographic: number[], result?: number[]): number[];\n\n cartographicToCartesian(cartographic, result = [0, 0, 0]) {\n const normal = scratchNormal;\n const k = scratchK;\n\n const [, , height] = cartographic;\n this.geodeticSurfaceNormalCartographic(cartographic, normal);\n k.copy(this.radiiSquared).scale(normal);\n\n const gamma = Math.sqrt(normal.dot(k));\n k.scale(1 / gamma);\n\n normal.scale(height);\n\n k.add(normal);\n\n return k.to(result);\n }\n\n /** Converts the provided cartesian to cartographic (lng/lat/z) representation.\n * The cartesian is undefined at the center of the ellipsoid. */\n cartesianToCartographic(cartesian: number[], result: Vector3): Vector3;\n cartesianToCartographic(cartesian: number[], result?: number[]): number[];\n\n cartesianToCartographic(cartesian, result = [0, 0, 0]) {\n scratchCartesian.from(cartesian);\n const point = this.scaleToGeodeticSurface(scratchCartesian, scratchPosition);\n\n if (!point) {\n return undefined;\n }\n\n const normal = this.geodeticSurfaceNormal(point, scratchNormal);\n\n const h = scratchHeight;\n h.copy(scratchCartesian).subtract(point);\n\n const longitude = Math.atan2(normal.y, normal.x);\n const latitude = Math.asin(normal.z);\n const height = Math.sign(vec3.dot(h, scratchCartesian)) * vec3.length(h);\n\n return toCartographicFromRadians([longitude, latitude, height], result);\n }\n\n /** Computes a 4x4 transformation matrix from a reference frame with an east-north-up axes\n * centered at the provided origin to the provided ellipsoid's fixed reference frame. */\n eastNorthUpToFixedFrame(origin: number[], result?: Matrix4): Matrix4;\n eastNorthUpToFixedFrame(origin: number[], result: number[]): number[];\n\n eastNorthUpToFixedFrame(origin, result = new Matrix4()) {\n return localFrameToFixedFrame(this, 'east', 'north', 'up', origin, result);\n }\n\n /** Computes a 4x4 transformation matrix from a reference frame centered at\n * the provided origin to the ellipsoid's fixed reference frame.\n */\n localFrameToFixedFrame(\n firstAxis: string,\n secondAxis: string,\n thirdAxis: string,\n origin: Readonly<NumericArray>,\n result?: Matrix4\n ): Matrix4;\n localFrameToFixedFrame<Matrix4T>(\n firstAxis: string,\n secondAxis: string,\n thirdAxis: string,\n origin: Readonly<NumericArray>,\n result: number[]\n ): number[];\n\n // Computes a 4x4 transformation matrix from a reference frame centered at\n // the provided origin to the ellipsoid's fixed reference frame.\n localFrameToFixedFrame(firstAxis, secondAxis, thirdAxis, origin, result = new Matrix4()) {\n return localFrameToFixedFrame(this, firstAxis, secondAxis, thirdAxis, origin, result);\n }\n\n /** Computes the unit vector directed from the center of this ellipsoid toward\n * the provided Cartesian position. */\n geocentricSurfaceNormal(cartesian: number[], result?: number[]): number[];\n geocentricSurfaceNormal<NumArray>(cartesian: number[], result: NumArray): NumArray;\n geocentricSurfaceNormal(cartesian, result = [0, 0, 0]) {\n return scratchVector.from(cartesian).normalize().to(result);\n }\n\n /** Computes the normal of the plane tangent to the surface of the ellipsoid at provided position. */\n geodeticSurfaceNormalCartographic<NumArray>(cartographic: number[], result: NumArray): NumArray;\n geodeticSurfaceNormalCartographic(cartographic: number[]): number[];\n geodeticSurfaceNormalCartographic(cartographic, result = [0, 0, 0]) {\n const cartographicVectorRadians = fromCartographicToRadians(cartographic);\n\n const longitude = cartographicVectorRadians[0];\n const latitude = cartographicVectorRadians[1];\n\n const cosLatitude = Math.cos(latitude);\n\n scratchVector\n .set(cosLatitude * Math.cos(longitude), cosLatitude * Math.sin(longitude), Math.sin(latitude))\n .normalize();\n\n return scratchVector.to(result);\n }\n\n /** Computes the normal of the plane tangent to the surface of the ellipsoid at the provided position. */\n geodeticSurfaceNormal<NumArrayT>(cartesian: number[], result: NumArrayT): NumArrayT;\n geodeticSurfaceNormal(cartesian: number[]): number[];\n geodeticSurfaceNormal(cartesian, result = [0, 0, 0]) {\n return scratchVector.from(cartesian).scale(this.oneOverRadiiSquared).normalize().to(result);\n }\n\n /** Scales the provided Cartesian position along the geodetic surface normal\n * so that it is on the surface of this ellipsoid. If the position is\n * at the center of the ellipsoid, this function returns undefined. */\n scaleToGeodeticSurface(cartesian: number[], result?: number[]): number[] {\n return scaleToGeodeticSurface(cartesian, this, result);\n }\n\n /** Scales the provided Cartesian position along the geocentric surface normal\n * so that it is on the surface of this ellipsoid. */\n scaleToGeocentricSurface(cartesian: number[], result: number[] = [0, 0, 0]): number[] {\n scratchPosition.from(cartesian);\n\n const positionX = scratchPosition.x;\n const positionY = scratchPosition.y;\n const positionZ = scratchPosition.z;\n const oneOverRadiiSquared = this.oneOverRadiiSquared;\n\n const beta =\n 1.0 /\n Math.sqrt(\n positionX * positionX * oneOverRadiiSquared.x +\n positionY * positionY * oneOverRadiiSquared.y +\n positionZ * positionZ * oneOverRadiiSquared.z\n );\n\n return scratchPosition.multiplyScalar(beta).to(result);\n }\n\n /** Transforms a Cartesian X, Y, Z position to the ellipsoid-scaled space by multiplying\n * its components by the result of `Ellipsoid#oneOverRadii` */\n transformPositionToScaledSpace(position: number[], result: number[] = [0, 0, 0]): number[] {\n return scratchPosition.from(position).scale(this.oneOverRadii).to(result);\n }\n\n /** Transforms a Cartesian X, Y, Z position from the ellipsoid-scaled space by multiplying\n * its components by the result of `Ellipsoid#radii`. */\n transformPositionFromScaledSpace(position: number[], result: number[] = [0, 0, 0]): number[] {\n return scratchPosition.from(position).scale(this.radii).to(result);\n }\n\n /** Computes a point which is the intersection of the surface normal with the z-axis. */\n getSurfaceNormalIntersectionWithZAxis(\n position: number[],\n buffer: number = 0,\n result: number[] = [0, 0, 0]\n ): number[] {\n // Ellipsoid must be an ellipsoid of revolution (radii.x == radii.y)\n assert(equals(this.radii.x, this.radii.y, _MathUtils.EPSILON15));\n assert(this.radii.z > 0);\n\n scratchPosition.from(position);\n const z = scratchPosition.z * (1 - this.squaredXOverSquaredZ);\n\n if (Math.abs(z) >= this.radii.z - buffer) {\n return undefined;\n }\n\n return scratchPosition.set(0.0, 0.0, z).to(result);\n }\n}\n"],"file":"ellipsoid.js"}
@@ -1 +0,0 @@
1
- {"version":3,"sources":["../../../../src/ellipsoid/helpers/ellipsoid-transform.ts"],"names":["Vector3","assert","equals","equalsEpsilon","EPSILON14","scratchOrigin","VECTOR_PRODUCT_LOCAL_FRAME","up","south","north","west","east","down","degeneratePositionLocalFrame","scratchAxisVectors","scratchVector1","scratchVector2","scratchVector3","localFrameToFixedFrame","ellipsoid","firstAxis","secondAxis","thirdAxis","cartesianOrigin","result","thirdAxisInferred","firstAxisVector","secondAxisVector","thirdAxisVector","origin","copy","atPole","x","y","sign","Math","z","fromArray","scale","set","normalize","geodeticSurfaceNormal","cross"],"mappings":"AAAA,SAAQA,OAAR,EAAiBC,MAAjB,EAAyBC,MAAM,IAAIC,aAAnC,QAAuD,eAAvD;AAIA,MAAMC,SAAS,GAAG,KAAlB;AAEA,MAAMC,aAAa,GAAG,IAAIL,OAAJ,EAAtB;AAGA,MAAMM,0BAA0B,GAAG;AACjCC,EAAAA,EAAE,EAAE;AACFC,IAAAA,KAAK,EAAE,MADL;AAEFC,IAAAA,KAAK,EAAE,MAFL;AAGFC,IAAAA,IAAI,EAAE,OAHJ;AAIFC,IAAAA,IAAI,EAAE;AAJJ,GAD6B;AAOjCC,EAAAA,IAAI,EAAE;AACJJ,IAAAA,KAAK,EAAE,MADH;AAEJC,IAAAA,KAAK,EAAE,MAFH;AAGJC,IAAAA,IAAI,EAAE,OAHF;AAIJC,IAAAA,IAAI,EAAE;AAJF,GAP2B;AAajCH,EAAAA,KAAK,EAAE;AACLD,IAAAA,EAAE,EAAE,MADC;AAELK,IAAAA,IAAI,EAAE,MAFD;AAGLF,IAAAA,IAAI,EAAE,MAHD;AAILC,IAAAA,IAAI,EAAE;AAJD,GAb0B;AAmBjCF,EAAAA,KAAK,EAAE;AACLF,IAAAA,EAAE,EAAE,MADC;AAELK,IAAAA,IAAI,EAAE,MAFD;AAGLF,IAAAA,IAAI,EAAE,IAHD;AAILC,IAAAA,IAAI,EAAE;AAJD,GAnB0B;AAyBjCD,EAAAA,IAAI,EAAE;AACJH,IAAAA,EAAE,EAAE,OADA;AAEJK,IAAAA,IAAI,EAAE,OAFF;AAGJH,IAAAA,KAAK,EAAE,MAHH;AAIJD,IAAAA,KAAK,EAAE;AAJH,GAzB2B;AA+BjCG,EAAAA,IAAI,EAAE;AACJJ,IAAAA,EAAE,EAAE,OADA;AAEJK,IAAAA,IAAI,EAAE,OAFF;AAGJH,IAAAA,KAAK,EAAE,IAHH;AAIJD,IAAAA,KAAK,EAAE;AAJH;AA/B2B,CAAnC;AAuCA,MAAMK,4BAA4B,GAAG;AACnCJ,EAAAA,KAAK,EAAE,CAAC,CAAC,CAAF,EAAK,CAAL,EAAQ,CAAR,CAD4B;AAEnCE,EAAAA,IAAI,EAAE,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAF6B;AAGnCJ,EAAAA,EAAE,EAAE,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAH+B;AAInCC,EAAAA,KAAK,EAAE,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAP,CAJ4B;AAKnCE,EAAAA,IAAI,EAAE,CAAC,CAAD,EAAI,CAAC,CAAL,EAAQ,CAAR,CAL6B;AAMnCE,EAAAA,IAAI,EAAE,CAAC,CAAD,EAAI,CAAJ,EAAO,CAAC,CAAR;AAN6B,CAArC;AASA,MAAME,kBAAkB,GAAG;AACzBH,EAAAA,IAAI,EAAE,IAAIX,OAAJ,EADmB;AAEzBS,EAAAA,KAAK,EAAE,IAAIT,OAAJ,EAFkB;AAGzBO,EAAAA,EAAE,EAAE,IAAIP,OAAJ,EAHqB;AAIzBU,EAAAA,IAAI,EAAE,IAAIV,OAAJ,EAJmB;AAKzBQ,EAAAA,KAAK,EAAE,IAAIR,OAAJ,EALkB;AAMzBY,EAAAA,IAAI,EAAE,IAAIZ,OAAJ;AANmB,CAA3B;AASA,MAAMe,cAAc,GAAG,IAAIf,OAAJ,EAAvB;AACA,MAAMgB,cAAc,GAAG,IAAIhB,OAAJ,EAAvB;AACA,MAAMiB,cAAc,GAAG,IAAIjB,OAAJ,EAAvB;AAOA,eAAe,SAASkB,sBAAT,CACbC,SADa,EAEbC,SAFa,EAGbC,UAHa,EAIbC,SAJa,EAKbC,eALa,EAMbC,MANa,EAOH;AACV,QAAMC,iBAAiB,GACrBnB,0BAA0B,CAACc,SAAD,CAA1B,IACCd,0BAA0B,CAACc,SAAD,CAA1B,CAAsCC,UAAtC,CAFH;AAIApB,EAAAA,MAAM,CAACwB,iBAAiB,KAAK,CAACH,SAAD,IAAcA,SAAS,KAAKG,iBAAjC,CAAlB,CAAN;AAEA,MAAIC,eAAJ;AACA,MAAIC,gBAAJ;AACA,MAAIC,eAAJ;AAEA,QAAMC,MAAM,GAAGxB,aAAa,CAACyB,IAAd,CAAmBP,eAAnB,CAAf;AAGA,QAAMQ,MAAM,GAAG5B,aAAa,CAAC0B,MAAM,CAACG,CAAR,EAAW,GAAX,EAAgB5B,SAAhB,CAAb,IAA2CD,aAAa,CAAC0B,MAAM,CAACI,CAAR,EAAW,GAAX,EAAgB7B,SAAhB,CAAvE;;AAEA,MAAI2B,MAAJ,EAAY;AAEV,UAAMG,IAAI,GAAGC,IAAI,CAACD,IAAL,CAAUL,MAAM,CAACO,CAAjB,CAAb;AAEAV,IAAAA,eAAe,GAAGX,cAAc,CAACsB,SAAf,CAAyBxB,4BAA4B,CAACO,SAAD,CAArD,CAAlB;;AACA,QAAIA,SAAS,KAAK,MAAd,IAAwBA,SAAS,KAAK,MAA1C,EAAkD;AAChDM,MAAAA,eAAe,CAACY,KAAhB,CAAsBJ,IAAtB;AACD;;AAEDP,IAAAA,gBAAgB,GAAGX,cAAc,CAACqB,SAAf,CAAyBxB,4BAA4B,CAACQ,UAAD,CAArD,CAAnB;;AACA,QAAIA,UAAU,KAAK,MAAf,IAAyBA,UAAU,KAAK,MAA5C,EAAoD;AAClDM,MAAAA,gBAAgB,CAACW,KAAjB,CAAuBJ,IAAvB;AACD;;AAEDN,IAAAA,eAAe,GAAGX,cAAc,CAACoB,SAAf,CAAyBxB,4BAA4B,CAACS,SAAD,CAArD,CAAlB;;AACA,QAAIA,SAAS,KAAK,MAAd,IAAwBA,SAAS,KAAK,MAA1C,EAAkD;AAChDM,MAAAA,eAAe,CAACU,KAAhB,CAAsBJ,IAAtB;AACD;AACF,GAlBD,MAkBO;AAEL,UAAM;AAAC3B,MAAAA,EAAD;AAAKI,MAAAA,IAAL;AAAWF,MAAAA;AAAX,QAAoBK,kBAA1B;AAEAH,IAAAA,IAAI,CAAC4B,GAAL,CAAS,CAACV,MAAM,CAACI,CAAjB,EAAoBJ,MAAM,CAACG,CAA3B,EAA8B,GAA9B,EAAmCQ,SAAnC;AACArB,IAAAA,SAAS,CAACsB,qBAAV,CAAgCZ,MAAhC,EAAwCtB,EAAxC;AACAE,IAAAA,KAAK,CAACqB,IAAN,CAAWvB,EAAX,EAAemC,KAAf,CAAqB/B,IAArB;AAEA,UAAM;AAACC,MAAAA,IAAD;AAAOF,MAAAA,IAAP;AAAaF,MAAAA;AAAb,QAAsBM,kBAA5B;AAEAF,IAAAA,IAAI,CAACkB,IAAL,CAAUvB,EAAV,EAAc+B,KAAd,CAAoB,CAAC,CAArB;AACA5B,IAAAA,IAAI,CAACoB,IAAL,CAAUnB,IAAV,EAAgB2B,KAAhB,CAAsB,CAAC,CAAvB;AACA9B,IAAAA,KAAK,CAACsB,IAAN,CAAWrB,KAAX,EAAkB6B,KAAlB,CAAwB,CAAC,CAAzB;AAGAZ,IAAAA,eAAe,GAAGZ,kBAAkB,CAACM,SAAD,CAApC;AACAO,IAAAA,gBAAgB,GAAGb,kBAAkB,CAACO,UAAD,CAArC;AACAO,IAAAA,eAAe,GAAGd,kBAAkB,CAACQ,SAAD,CAApC;AACD;;AAGDE,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYE,eAAe,CAACM,CAA5B;AACAR,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYE,eAAe,CAACO,CAA5B;AACAT,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYE,eAAe,CAACU,CAA5B;AACAZ,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAY,GAAZ;AACAA,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYG,gBAAgB,CAACK,CAA7B;AACAR,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYG,gBAAgB,CAACM,CAA7B;AACAT,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYG,gBAAgB,CAACS,CAA7B;AACAZ,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAY,GAAZ;AACAA,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYI,eAAe,CAACI,CAA5B;AACAR,EAAAA,MAAM,CAAC,CAAD,CAAN,GAAYI,eAAe,CAACK,CAA5B;AACAT,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAaI,eAAe,CAACQ,CAA7B;AACAZ,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAa,GAAb;AACAA,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAaK,MAAM,CAACG,CAApB;AACAR,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAaK,MAAM,CAACI,CAApB;AACAT,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAaK,MAAM,CAACO,CAApB;AACAZ,EAAAA,MAAM,CAAC,EAAD,CAAN,GAAa,GAAb;AACA,SAAOA,MAAP;AACD","sourcesContent":["import {Vector3, assert, equals as equalsEpsilon} from '@math.gl/core';\n\nimport type Ellipsoid from '../ellipsoid';\n\nconst EPSILON14 = 1e-14;\n\nconst scratchOrigin = new Vector3();\n\n// Caclulate third axis from given two axii\nconst VECTOR_PRODUCT_LOCAL_FRAME = {\n up: {\n south: 'east',\n north: 'west',\n west: 'south',\n east: 'north'\n },\n down: {\n south: 'west',\n north: 'east',\n west: 'north',\n east: 'south'\n },\n south: {\n up: 'west',\n down: 'east',\n west: 'down',\n east: 'up'\n },\n north: {\n up: 'east',\n down: 'west',\n west: 'up',\n east: 'down'\n },\n west: {\n up: 'north',\n down: 'south',\n north: 'down',\n south: 'up'\n },\n east: {\n up: 'south',\n down: 'north',\n north: 'up',\n south: 'down'\n }\n} as const;\n\nconst degeneratePositionLocalFrame = {\n north: [-1, 0, 0],\n east: [0, 1, 0],\n up: [0, 0, 1],\n south: [1, 0, 0],\n west: [0, -1, 0],\n down: [0, 0, -1]\n} as const;\n\nconst scratchAxisVectors = {\n east: new Vector3(),\n north: new Vector3(),\n up: new Vector3(),\n west: new Vector3(),\n south: new Vector3(),\n down: new Vector3()\n};\n\nconst scratchVector1 = new Vector3();\nconst scratchVector2 = new Vector3();\nconst scratchVector3 = new Vector3();\n\ntype Axis = 'up' | 'down' | 'north' | 'east' | 'south' | 'west';\n\n// Computes a 4x4 transformation matrix from a reference frame\n// centered at the provided origin to the provided ellipsoid's fixed reference frame.\n// eslint-disable-next-line max-statements, max-params, complexity\nexport default function localFrameToFixedFrame(\n ellipsoid: Ellipsoid,\n firstAxis: Axis,\n secondAxis: Axis,\n thirdAxis: Axis,\n cartesianOrigin: number[],\n result: number[]\n): number[] {\n const thirdAxisInferred =\n VECTOR_PRODUCT_LOCAL_FRAME[firstAxis] &&\n (VECTOR_PRODUCT_LOCAL_FRAME[firstAxis][secondAxis] as Axis);\n // firstAxis and secondAxis must be east, north, up, west, south or down.');\n assert(thirdAxisInferred && (!thirdAxis || thirdAxis === thirdAxisInferred));\n\n let firstAxisVector: Vector3;\n let secondAxisVector: Vector3;\n let thirdAxisVector: Vector3;\n\n const origin = scratchOrigin.copy(cartesianOrigin);\n\n // If x and y are zero, assume origin is at a pole, which is a special case.\n const atPole = equalsEpsilon(origin.x, 0.0, EPSILON14) && equalsEpsilon(origin.y, 0.0, EPSILON14);\n\n if (atPole) {\n // Look up axis value and adjust\n const sign = Math.sign(origin.z);\n\n firstAxisVector = scratchVector1.fromArray(degeneratePositionLocalFrame[firstAxis]);\n if (firstAxis !== 'east' && firstAxis !== 'west') {\n firstAxisVector.scale(sign);\n }\n\n secondAxisVector = scratchVector2.fromArray(degeneratePositionLocalFrame[secondAxis]);\n if (secondAxis !== 'east' && secondAxis !== 'west') {\n secondAxisVector.scale(sign);\n }\n\n thirdAxisVector = scratchVector3.fromArray(degeneratePositionLocalFrame[thirdAxis]);\n if (thirdAxis !== 'east' && thirdAxis !== 'west') {\n thirdAxisVector.scale(sign);\n }\n } else {\n // Calculate all axis\n const {up, east, north} = scratchAxisVectors;\n\n east.set(-origin.y, origin.x, 0.0).normalize();\n ellipsoid.geodeticSurfaceNormal(origin, up);\n north.copy(up).cross(east);\n\n const {down, west, south} = scratchAxisVectors;\n\n down.copy(up).scale(-1);\n west.copy(east).scale(-1);\n south.copy(north).scale(-1);\n\n // Pick three axis based on desired orientation\n firstAxisVector = scratchAxisVectors[firstAxis];\n secondAxisVector = scratchAxisVectors[secondAxis];\n thirdAxisVector = scratchAxisVectors[thirdAxis];\n }\n\n // TODO - assuming the result is column-major\n result[0] = firstAxisVector.x;\n result[1] = firstAxisVector.y;\n result[2] = firstAxisVector.z;\n result[3] = 0.0;\n result[4] = secondAxisVector.x;\n result[5] = secondAxisVector.y;\n result[6] = secondAxisVector.z;\n result[7] = 0.0;\n result[8] = thirdAxisVector.x;\n result[9] = thirdAxisVector.y;\n result[10] = thirdAxisVector.z;\n result[11] = 0.0;\n result[12] = origin.x;\n result[13] = origin.y;\n result[14] = origin.z;\n result[15] = 1.0;\n return result;\n}\n"],"file":"ellipsoid-transform.js"}
@@ -1 +0,0 @@
1
- {"version":3,"sources":["../../../../src/ellipsoid/helpers/scale-to-geodetic-surface.ts"],"names":["Vector3","_MathUtils","scratchVector","scaleToGeodeticSurfaceIntersection","scaleToGeodeticSurfaceGradient","scaleToGeodeticSurface","cartesian","ellipsoid","result","oneOverRadii","oneOverRadiiSquared","centerToleranceSquared","from","positionX","x","positionY","y","positionZ","z","oneOverRadiiX","oneOverRadiiY","oneOverRadiiZ","x2","y2","z2","squaredNorm","ratio","Math","sqrt","Number","isFinite","undefined","intersection","copy","scale","to","oneOverRadiiSquaredX","oneOverRadiiSquaredY","oneOverRadiiSquaredZ","gradient","set","lambda","len","correction","xMultiplier","yMultiplier","zMultiplier","func","xMultiplier2","yMultiplier2","zMultiplier2","xMultiplier3","yMultiplier3","zMultiplier3","denominator","derivative","abs","EPSILON12"],"mappings":"AACA,SAAQA,OAAR,EAAiBC,UAAjB,QAAkC,eAAlC;AAGA,MAAMC,aAAa,GAAG,IAAIF,OAAJ,EAAtB;AACA,MAAMG,kCAAkC,GAAG,IAAIH,OAAJ,EAA3C;AACA,MAAMI,8BAA8B,GAAG,IAAIJ,OAAJ,EAAvC;AAKA,eAAe,SAASK,sBAAT,CACbC,SADa,EAEbC,SAFa,EAGbC,MAAgB,GAAG,EAHN,EAIH;AACV,QAAM;AAACC,IAAAA,YAAD;AAAeC,IAAAA,mBAAf;AAAoCC,IAAAA;AAApC,MAA8DJ,SAApE;AAEAL,EAAAA,aAAa,CAACU,IAAd,CAAmBN,SAAnB;AAEA,QAAMO,SAAS,GAAGX,aAAa,CAACY,CAAhC;AACA,QAAMC,SAAS,GAAGb,aAAa,CAACc,CAAhC;AACA,QAAMC,SAAS,GAAGf,aAAa,CAACgB,CAAhC;AAEA,QAAMC,aAAa,GAAGV,YAAY,CAACK,CAAnC;AACA,QAAMM,aAAa,GAAGX,YAAY,CAACO,CAAnC;AACA,QAAMK,aAAa,GAAGZ,YAAY,CAACS,CAAnC;AAEA,QAAMI,EAAE,GAAGT,SAAS,GAAGA,SAAZ,GAAwBM,aAAxB,GAAwCA,aAAnD;AACA,QAAMI,EAAE,GAAGR,SAAS,GAAGA,SAAZ,GAAwBK,aAAxB,GAAwCA,aAAnD;AACA,QAAMI,EAAE,GAAGP,SAAS,GAAGA,SAAZ,GAAwBI,aAAxB,GAAwCA,aAAnD;AAGA,QAAMI,WAAW,GAAGH,EAAE,GAAGC,EAAL,GAAUC,EAA9B;AACA,QAAME,KAAK,GAAGC,IAAI,CAACC,IAAL,CAAU,MAAMH,WAAhB,CAAd;;AAGA,MAAI,CAACI,MAAM,CAACC,QAAP,CAAgBJ,KAAhB,CAAL,EAA6B;AAC3B,WAAOK,SAAP;AACD;;AAGD,QAAMC,YAAY,GAAG7B,kCAArB;AACA6B,EAAAA,YAAY,CAACC,IAAb,CAAkB3B,SAAlB,EAA6B4B,KAA7B,CAAmCR,KAAnC;;AAGA,MAAID,WAAW,GAAGd,sBAAlB,EAA0C;AACxC,WAAOqB,YAAY,CAACG,EAAb,CAAgB3B,MAAhB,CAAP;AACD;;AAED,QAAM4B,oBAAoB,GAAG1B,mBAAmB,CAACI,CAAjD;AACA,QAAMuB,oBAAoB,GAAG3B,mBAAmB,CAACM,CAAjD;AACA,QAAMsB,oBAAoB,GAAG5B,mBAAmB,CAACQ,CAAjD;AAIA,QAAMqB,QAAQ,GAAGnC,8BAAjB;AACAmC,EAAAA,QAAQ,CAACC,GAAT,CACER,YAAY,CAAClB,CAAb,GAAiBsB,oBAAjB,GAAwC,GAD1C,EAEEJ,YAAY,CAAChB,CAAb,GAAiBqB,oBAAjB,GAAwC,GAF1C,EAGEL,YAAY,CAACd,CAAb,GAAiBoB,oBAAjB,GAAwC,GAH1C;AAOA,MAAIG,MAAM,GAAI,CAAC,MAAMf,KAAP,IAAgBxB,aAAa,CAACwC,GAAd,EAAjB,IAAyC,MAAMH,QAAQ,CAACG,GAAT,EAA/C,CAAb;AACA,MAAIC,UAAU,GAAG,GAAjB;AAEA,MAAIC,WAAJ;AACA,MAAIC,WAAJ;AACA,MAAIC,WAAJ;AACA,MAAIC,IAAJ;;AAEA,KAAG;AACDN,IAAAA,MAAM,IAAIE,UAAV;AAEAC,IAAAA,WAAW,GAAG,OAAO,MAAMH,MAAM,GAAGL,oBAAtB,CAAd;AACAS,IAAAA,WAAW,GAAG,OAAO,MAAMJ,MAAM,GAAGJ,oBAAtB,CAAd;AACAS,IAAAA,WAAW,GAAG,OAAO,MAAML,MAAM,GAAGH,oBAAtB,CAAd;AAEA,UAAMU,YAAY,GAAGJ,WAAW,GAAGA,WAAnC;AACA,UAAMK,YAAY,GAAGJ,WAAW,GAAGA,WAAnC;AACA,UAAMK,YAAY,GAAGJ,WAAW,GAAGA,WAAnC;AAEA,UAAMK,YAAY,GAAGH,YAAY,GAAGJ,WAApC;AACA,UAAMQ,YAAY,GAAGH,YAAY,GAAGJ,WAApC;AACA,UAAMQ,YAAY,GAAGH,YAAY,GAAGJ,WAApC;AAEAC,IAAAA,IAAI,GAAGzB,EAAE,GAAG0B,YAAL,GAAoBzB,EAAE,GAAG0B,YAAzB,GAAwCzB,EAAE,GAAG0B,YAA7C,GAA4D,GAAnE;AAIA,UAAMI,WAAW,GACfhC,EAAE,GAAG6B,YAAL,GAAoBf,oBAApB,GACAb,EAAE,GAAG6B,YAAL,GAAoBf,oBADpB,GAEAb,EAAE,GAAG6B,YAAL,GAAoBf,oBAHtB;AAKA,UAAMiB,UAAU,GAAG,CAAC,GAAD,GAAOD,WAA1B;AAEAX,IAAAA,UAAU,GAAGI,IAAI,GAAGQ,UAApB;AACD,GA3BD,QA2BS5B,IAAI,CAAC6B,GAAL,CAAST,IAAT,IAAiB9C,UAAU,CAACwD,SA3BrC;;AA6BA,SAAOvD,aAAa,CAACgC,KAAd,CAAoB,CAACU,WAAD,EAAcC,WAAd,EAA2BC,WAA3B,CAApB,EAA6DX,EAA7D,CAAgE3B,MAAhE,CAAP;AACD","sourcesContent":["/* eslint-disable */\nimport {Vector3, _MathUtils} from '@math.gl/core';\nimport type Ellipsoid from '../ellipsoid';\n\nconst scratchVector = new Vector3();\nconst scaleToGeodeticSurfaceIntersection = new Vector3();\nconst scaleToGeodeticSurfaceGradient = new Vector3();\n\n// Scales the provided Cartesian position along the geodetic surface normal\n// so that it is on the surface of this ellipsoid. If the position is\n// at the center of the ellipsoid, this function returns undefined.\nexport default function scaleToGeodeticSurface(\n cartesian: number[],\n ellipsoid: Ellipsoid,\n result: number[] = []\n): number[] {\n const {oneOverRadii, oneOverRadiiSquared, centerToleranceSquared} = ellipsoid;\n\n scratchVector.from(cartesian);\n\n const positionX = scratchVector.x;\n const positionY = scratchVector.y;\n const positionZ = scratchVector.z;\n\n const oneOverRadiiX = oneOverRadii.x;\n const oneOverRadiiY = oneOverRadii.y;\n const oneOverRadiiZ = oneOverRadii.z;\n\n const x2 = positionX * positionX * oneOverRadiiX * oneOverRadiiX;\n const y2 = positionY * positionY * oneOverRadiiY * oneOverRadiiY;\n const z2 = positionZ * positionZ * oneOverRadiiZ * oneOverRadiiZ;\n\n // Compute the squared ellipsoid norm.\n const squaredNorm = x2 + y2 + z2;\n const ratio = Math.sqrt(1.0 / squaredNorm);\n\n // When very close to center or at center\n if (!Number.isFinite(ratio)) {\n return undefined;\n }\n\n // As an initial approximation, assume that the radial intersection is the projection point.\n const intersection = scaleToGeodeticSurfaceIntersection;\n intersection.copy(cartesian).scale(ratio);\n\n // If the position is near the center, the iteration will not converge.\n if (squaredNorm < centerToleranceSquared) {\n return intersection.to(result);\n }\n\n const oneOverRadiiSquaredX = oneOverRadiiSquared.x;\n const oneOverRadiiSquaredY = oneOverRadiiSquared.y;\n const oneOverRadiiSquaredZ = oneOverRadiiSquared.z;\n\n // Use the gradient at the intersection point in place of the true unit normal.\n // The difference in magnitude will be absorbed in the multiplier.\n const gradient = scaleToGeodeticSurfaceGradient;\n gradient.set(\n intersection.x * oneOverRadiiSquaredX * 2.0,\n intersection.y * oneOverRadiiSquaredY * 2.0,\n intersection.z * oneOverRadiiSquaredZ * 2.0\n );\n\n // Compute the initial guess at the normal vector multiplier, lambda.\n let lambda = ((1.0 - ratio) * scratchVector.len()) / (0.5 * gradient.len());\n let correction = 0.0;\n\n let xMultiplier;\n let yMultiplier;\n let zMultiplier;\n let func;\n\n do {\n lambda -= correction;\n\n xMultiplier = 1.0 / (1.0 + lambda * oneOverRadiiSquaredX);\n yMultiplier = 1.0 / (1.0 + lambda * oneOverRadiiSquaredY);\n zMultiplier = 1.0 / (1.0 + lambda * oneOverRadiiSquaredZ);\n\n const xMultiplier2 = xMultiplier * xMultiplier;\n const yMultiplier2 = yMultiplier * yMultiplier;\n const zMultiplier2 = zMultiplier * zMultiplier;\n\n const xMultiplier3 = xMultiplier2 * xMultiplier;\n const yMultiplier3 = yMultiplier2 * yMultiplier;\n const zMultiplier3 = zMultiplier2 * zMultiplier;\n\n func = x2 * xMultiplier2 + y2 * yMultiplier2 + z2 * zMultiplier2 - 1.0;\n\n // \"denominator\" here refers to the use of this expression in the velocity and acceleration\n // computations in the sections to follow.\n const denominator =\n x2 * xMultiplier3 * oneOverRadiiSquaredX +\n y2 * yMultiplier3 * oneOverRadiiSquaredY +\n z2 * zMultiplier3 * oneOverRadiiSquaredZ;\n\n const derivative = -2.0 * denominator;\n\n correction = func / derivative;\n } while (Math.abs(func) > _MathUtils.EPSILON12);\n\n return scratchVector.scale([xMultiplier, yMultiplier, zMultiplier]).to(result);\n}\n"],"file":"scale-to-geodetic-surface.js"}
package/dist/esm/index.js DELETED
@@ -1,3 +0,0 @@
1
- export { default as Ellipsoid } from './ellipsoid/ellipsoid';
2
- export { isWGS84 } from './type-utils';
3
- //# sourceMappingURL=index.js.map
@@ -1 +0,0 @@
1
- {"version":3,"sources":["../../src/index.ts"],"names":["default","Ellipsoid","isWGS84"],"mappings":"AAAA,SAAQA,OAAO,IAAIC,SAAnB,QAAmC,uBAAnC;AACA,SAAQC,OAAR,QAAsB,cAAtB","sourcesContent":["export {default as Ellipsoid} from './ellipsoid/ellipsoid';\nexport {isWGS84} from './type-utils';\n"],"file":"index.js"}
@@ -1 +0,0 @@
1
- {"version":3,"sources":["../../src/type-utils.ts"],"names":["Vector3","toRadians","toDegrees","config","WGS84_CONSTANTS","identity","x","scratchVector","fromCartographic","cartographic","result","map","longitude","latitude","height","y","z","fromCartographicToRadians","vector","_cartographicRadians","fromCartographicToDegrees","toCartographic","toCartographicFromRadians","toCartographicFromDegrees","isWGS84","from","oneOverRadiiSquared","centerToleranceSquared","x2","y2","z2","Math","abs"],"mappings":"AAIA,SAAQA,OAAR,EAAiBC,SAAjB,EAA4BC,SAA5B,EAAuCC,MAAvC,QAAoD,eAApD;AACA,SAAQC,eAAR,QAA8B,aAA9B;;AAgBA,SAASC,QAAT,CAAkBC,CAAlB,EAAqC;AACnC,SAAOA,CAAP;AACD;;AAED,MAAMC,aAAa,GAAG,IAAIP,OAAJ,EAAtB;AAQA,OAAO,SAASQ,gBAAT,CACLC,YADK,EAELC,MAAM,GAAG,EAFJ,EAGLC,GAAG,GAAGN,QAHD,EAIK;AACV,MAAI,eAAeI,YAAnB,EAAiC;AAC/BC,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYC,GAAG,CAACF,YAAY,CAACG,SAAd,CAAf;AACAF,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYC,GAAG,CAACF,YAAY,CAACI,QAAd,CAAf;AACAH,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYD,YAAY,CAACK,MAAzB;AACD,GAJD,MAIO,IAAI,OAAOL,YAAX,EAAyB;AAC9BC,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYC,GAAG,CAACF,YAAY,CAACH,CAAd,CAAf;AACAI,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYC,GAAG,CAACF,YAAY,CAACM,CAAd,CAAf;AACAL,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYD,YAAY,CAACO,CAAzB;AACD,GAJM,MAIA;AACLN,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYC,GAAG,CAACF,YAAY,CAAC,CAAD,CAAb,CAAf;AACAC,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYC,GAAG,CAACF,YAAY,CAAC,CAAD,CAAb,CAAf;AACAC,IAAAA,MAAM,CAAC,CAAD,CAAN,GAAYD,YAAY,CAAC,CAAD,CAAxB;AACD;;AACD,SAAOC,MAAP;AACD;AAOD,OAAO,SAASO,yBAAT,CACLR,YADK,EAELS,MAAM,GAAG,EAFJ,EAGK;AACV,SAAOV,gBAAgB,CAACC,YAAD,EAAeS,MAAf,EAAuBf,MAAM,CAACgB,oBAAP,GAA8Bd,QAA9B,GAAyCJ,SAAhE,CAAvB;AACD;AAOD,OAAO,SAASmB,yBAAT,CACLX,YADK,EAELS,MAAM,GAAG,EAFJ,EAGK;AACV,SAAOV,gBAAgB,CAACC,YAAD,EAAeS,MAAf,EAAuBf,MAAM,CAACgB,oBAAP,GAA8BjB,SAA9B,GAA0CG,QAAjE,CAAvB;AACD;AAED,OAAO,SAASgB,cAAT,CACLH,MADK,EAELT,YAFK,EAGLE,GAA0B,GAAGN,QAHxB,EAIF;AACH,MAAI,eAAeI,YAAnB,EAAiC;AAC/BA,IAAAA,YAAY,CAACG,SAAb,GAAyBD,GAAG,CAACO,MAAM,CAAC,CAAD,CAAP,CAA5B;AACAT,IAAAA,YAAY,CAACI,QAAb,GAAwBF,GAAG,CAACO,MAAM,CAAC,CAAD,CAAP,CAA3B;AACAT,IAAAA,YAAY,CAACK,MAAb,GAAsBI,MAAM,CAAC,CAAD,CAA5B;AACD,GAJD,MAIO,IAAI,OAAOT,YAAX,EAAyB;AAC9BA,IAAAA,YAAY,CAACH,CAAb,GAAiBK,GAAG,CAACO,MAAM,CAAC,CAAD,CAAP,CAApB;AACAT,IAAAA,YAAY,CAACM,CAAb,GAAiBJ,GAAG,CAACO,MAAM,CAAC,CAAD,CAAP,CAApB;AACAT,IAAAA,YAAY,CAACO,CAAb,GAAiBE,MAAM,CAAC,CAAD,CAAvB;AACD,GAJM,MAIA;AACLT,IAAAA,YAAY,CAAC,CAAD,CAAZ,GAAkBE,GAAG,CAACO,MAAM,CAAC,CAAD,CAAP,CAArB;AACAT,IAAAA,YAAY,CAAC,CAAD,CAAZ,GAAkBE,GAAG,CAACO,MAAM,CAAC,CAAD,CAAP,CAArB;AACAT,IAAAA,YAAY,CAAC,CAAD,CAAZ,GAAkBS,MAAM,CAAC,CAAD,CAAxB;AACD;;AACD,SAAOT,YAAP;AACD;AAED,OAAO,SAASa,yBAAT,CACLJ,MADK,EAELT,YAFK,EAGF;AACH,SAAOY,cAAc,CAACH,MAAD,EAAST,YAAT,EAAuBN,MAAM,CAACgB,oBAAP,GAA8Bd,QAA9B,GAAyCH,SAAhE,CAArB;AACD;AAED,OAAO,SAASqB,yBAAT,CACLL,MADK,EAELT,YAFK,EAGF;AACH,SAAOY,cAAc,CAACH,MAAD,EAAST,YAAT,EAAuBN,MAAM,CAACgB,oBAAP,GAA8BlB,SAA9B,GAA0CI,QAAjE,CAArB;AACD;AAGD,OAAO,SAASmB,OAAT,CAAiBN,MAAjB,EAA0D;AAC/D,MAAI,CAACA,MAAL,EAAa;AACX,WAAO,KAAP;AACD;;AACDX,EAAAA,aAAa,CAACkB,IAAd,CAAmBP,MAAnB;AACA,QAAM;AAACQ,IAAAA,mBAAD;AAAsBC,IAAAA;AAAtB,MAAgDvB,eAAtD;AACA,QAAMwB,EAAE,GAAGV,MAAM,CAAC,CAAD,CAAN,GAAYA,MAAM,CAAC,CAAD,CAAlB,GAAwBQ,mBAAmB,CAAC,CAAD,CAAtD;AACA,QAAMG,EAAE,GAAGX,MAAM,CAAC,CAAD,CAAN,GAAYA,MAAM,CAAC,CAAD,CAAlB,GAAwBQ,mBAAmB,CAAC,CAAD,CAAtD;AACA,QAAMI,EAAE,GAAGZ,MAAM,CAAC,CAAD,CAAN,GAAYA,MAAM,CAAC,CAAD,CAAlB,GAAwBQ,mBAAmB,CAAC,CAAD,CAAtD;AACA,SAAOK,IAAI,CAACC,GAAL,CAASJ,EAAE,GAAGC,EAAL,GAAUC,EAAV,GAAe,CAAxB,IAA6BH,sBAApC;AACD","sourcesContent":["// This file is derived from the Cesium math library under Apache 2 license\n// See LICENSE.md and https://github.com/AnalyticalGraphicsInc/cesium/blob/master/LICENSE.md\n\nimport type {NumericArray} from '@math.gl/core';\nimport {Vector3, toRadians, toDegrees, config} from '@math.gl/core';\nimport {WGS84_CONSTANTS} from './constants';\n\ntype LngLatHeightObject = {\n longitude: number;\n latitude: number;\n height: number;\n};\n\ntype XYZObject = {\n x: number;\n y: number;\n z: number;\n};\n\ntype Cartographic = LngLatHeightObject | XYZObject | NumericArray;\n\nfunction identity(x: number): number {\n return x;\n}\n\nconst scratchVector = new Vector3();\n\nexport function fromCartographic(cartographic: Cartographic): number[];\nexport function fromCartographic<NumArrayT>(\n cartographic: Cartographic,\n result: NumArrayT,\n map?: (x: number) => number\n): NumArrayT;\nexport function fromCartographic(\n cartographic: Cartographic,\n result = [] as number[],\n map = identity\n): number[] {\n if ('longitude' in cartographic) {\n result[0] = map(cartographic.longitude);\n result[1] = map(cartographic.latitude);\n result[2] = cartographic.height;\n } else if ('x' in cartographic) {\n result[0] = map(cartographic.x);\n result[1] = map(cartographic.y);\n result[2] = cartographic.z;\n } else {\n result[0] = map(cartographic[0]);\n result[1] = map(cartographic[1]);\n result[2] = cartographic[2];\n }\n return result;\n}\n\nexport function fromCartographicToRadians(cartographic: Cartographic, result?: number[]): number[];\nexport function fromCartographicToRadians<TArray>(\n cartographic: Cartographic,\n result: TArray\n): TArray;\nexport function fromCartographicToRadians(\n cartographic: Cartographic,\n vector = [] as number[]\n): number[] {\n return fromCartographic(cartographic, vector, config._cartographicRadians ? identity : toRadians);\n}\n\nexport function fromCartographicToDegrees(cartographic: Cartographic, result?: number[]): number[];\nexport function fromCartographicToDegrees<TArray>(\n cartographic: Cartographic,\n result: TArray\n): TArray;\nexport function fromCartographicToDegrees(\n cartographic: Cartographic,\n vector = [] as number[]\n): number[] {\n return fromCartographic(cartographic, vector, config._cartographicRadians ? toDegrees : identity);\n}\n\nexport function toCartographic<T extends Cartographic>(\n vector: Readonly<NumericArray>,\n cartographic: T,\n map: (x: number) => number = identity\n): T {\n if ('longitude' in cartographic) {\n cartographic.longitude = map(vector[0]);\n cartographic.latitude = map(vector[1]);\n cartographic.height = vector[2];\n } else if ('x' in cartographic) {\n cartographic.x = map(vector[0]);\n cartographic.y = map(vector[1]);\n cartographic.z = vector[2];\n } else {\n cartographic[0] = map(vector[0]);\n cartographic[1] = map(vector[1]);\n cartographic[2] = vector[2];\n }\n return cartographic;\n}\n\nexport function toCartographicFromRadians<T extends Cartographic>(\n vector: Readonly<NumericArray>,\n cartographic: T\n): T {\n return toCartographic(vector, cartographic, config._cartographicRadians ? identity : toDegrees);\n}\n\nexport function toCartographicFromDegrees<T extends Cartographic>(\n vector: Readonly<NumericArray>,\n cartographic: T\n): T {\n return toCartographic(vector, cartographic, config._cartographicRadians ? toRadians : identity);\n}\n\n// Estimates if a vector is close to the surface of the WGS84 Ellipsoid\nexport function isWGS84(vector: Readonly<NumericArray>): boolean {\n if (!vector) {\n return false;\n }\n scratchVector.from(vector);\n const {oneOverRadiiSquared, centerToleranceSquared} = WGS84_CONSTANTS;\n const x2 = vector[0] * vector[0] * oneOverRadiiSquared[0];\n const y2 = vector[1] * vector[1] * oneOverRadiiSquared[1];\n const z2 = vector[2] * vector[2] * oneOverRadiiSquared[2];\n return Math.abs(x2 + y2 + z2 - 1) < centerToleranceSquared;\n}\n\n/*\n\nexport function fromCartographic(cartographic: Cartographic, result?: number[]): number[];\nexport function fromCartographic(cartographic: Cartographic, result: TypedArray): TypedArray;\nexport function fromCartographicToRadians(cartographic: Cartographic, result?: number[]): number[];\nexport function fromCartographicToRadians(\n cartographic: Cartographic,\n result: TypedArray\n): TypedArray;\nexport function fromCartographicToDegrees(cartographic: Cartographic, result?: number[]): number[];\nexport function fromCartographicToDegrees(\n cartographic: Cartographic,\n result: TypedArray\n): TypedArray;\n\nexport function toCartographic(vector: number[] | TypedArray, result: Cartographic): number[];\nexport function toCartographicFromRadians(\n vector: number[] | TypedArray,\n result: Cartographic\n): number[];\nexport function toCartographicFromDegrees(\n vector: number[] | TypedArray,\n result: Cartographic\n): number[];\n\n// Estimates if a vector is close to the surface of the WGS84 Ellipsoid\nexport function isWGS84(vector: number[] | TypedArray): boolean;\n*/\n"],"file":"type-utils.js"}
File without changes