astro-viewer 1.3.0

package/lib-esm/shader/ShaderManagerMutliHiPS.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ShaderManagerMutliHiPS.js","sourceRoot":"","sources":["../../src/shader/ShaderManagerMutliHiPS.ts"],"names":[],"mappings":"AAAA,4BAA4B;AAI5B,MAAM,CAAC,OAAO,OAAO,sBAAsB;IACzC,MAAM,CAAC,MAAM;QACX,OAAO;;;;;;;;;;;;;MAaL,CAAC;IACL,CAAC;IAED,MAAM,CAAC,YAAY;QACjB,OAAO;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;MAyCL,CAAC;IACL,CAAC;IAED,MAAM,CAAC,eAAe;QACpB,OAAO;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;MA8FL,CAAC;IACL,CAAC;IAED,MAAM,CAAC,YAAY;QACjB,OAAO;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;MAqIL,CAAC;IACL,CAAC;IAED,MAAM,CAAC,SAAS;QACd,OAAO;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;MA+GL,CAAC;IACL,CAAC;IAED,MAAM,CAAC,aAAa;QAClB,OAAO;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;MA8GL,CAAC;IACL,CAAC;IAED,MAAM,CAAC,UAAU;QACf,OAAO;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;MA8GL,CAAC;IACL,CAAC;IAED,MAAM,CAAC,eAAe;QACpB,OAAO;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;MA8GL,CAAC;IACL,CAAC;CACF"}

package/lib-esm/utils/CoordsType.d.ts

@@ -0,0 +1,10 @@
+/**
+ * Enum for coordinate types.
+ * @author Fabrizio Giordano (Fab77)
+ */
+export declare enum CoordsType {
+    CARTESIAN = "cartesian",
+    SPHERICAL = "spherical",
+    ASTRO = "astro"
+}
+//# sourceMappingURL=CoordsType.d.ts.map

package/lib-esm/utils/CoordsType.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"CoordsType.d.ts","sourceRoot":"","sources":["../../src/utils/CoordsType.ts"],"names":[],"mappings":"AAAA;;;GAGG;AACH,oBAAY,UAAU;IACpB,SAAS,cAAc;IACvB,SAAS,cAAc;IACvB,KAAK,UAAU;CAChB"}

package/lib-esm/utils/CoordsType.js

@@ -0,0 +1,12 @@
+/**
+ * Enum for coordinate types.
+ * @author Fabrizio Giordano (Fab77)
+ */
+export var CoordsType;
+(function (CoordsType) {
+    CoordsType["CARTESIAN"] = "cartesian";
+    CoordsType["SPHERICAL"] = "spherical";
+    CoordsType["ASTRO"] = "astro";
+})(CoordsType || (CoordsType = {}));
+// export default CoordsType;
+//# sourceMappingURL=CoordsType.js.map

package/lib-esm/utils/CoordsType.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"CoordsType.js","sourceRoot":"","sources":["../../src/utils/CoordsType.ts"],"names":[],"mappings":"AAAA;;;GAGG;AACH,MAAM,CAAN,IAAY,UAIX;AAJD,WAAY,UAAU;IACpB,qCAAuB,CAAA;IACvB,qCAAuB,CAAA;IACvB,6BAAe,CAAA;AACjB,CAAC,EAJW,UAAU,KAAV,UAAU,QAIrB;AAED,6BAA6B"}

package/lib-esm/utils/FoVUtils.d.ts

@@ -0,0 +1,46 @@
+/**
+ * @author Fabrizio Giordano (Fab)
+ */
+import { Point } from '../model/Point.js';
+import { ReadonlyMat4 } from 'gl-matrix';
+import Camera from '../Camera.js';
+import { AbstractSkyEntity } from '../model/AbstractSkyEntity.js';
+import { HealpixGrid } from '../model/grid/HealpixGrid.js';
+export declare class FoVUtils {
+    /**
+     * Return the minimum FoV value between `_fovY_deg` and `_fovX_deg`.
+     * (Kept here for parity; this class doesn’t maintain those fields.)
+     */
+    getMinFoV(this: {
+        _fovY_deg: number;
+        _fovX_deg: number;
+    }): number;
+    /**
+     * Compute the FoV polygon as a list of Points (clockwise).
+     * Uses ray picking + frustum planes against a unit sphere.
+     */
+    static getFoVPolygon(camera: Camera, canvas: HTMLCanvasElement, model: AbstractSkyEntity, healpixGrid: HealpixGrid, webgl: WebGL2RenderingContext, pMatrix: ReadonlyMat4): Point[];
+    /**
+     * Ray pick against 8 key screen positions (corners + midpoints).
+     * Returns Points in clockwise order starting from top-left.
+     */
+    static getScreenCornersIntersection(pMatrix: ReadonlyMat4, camera: Camera, canvas: HTMLCanvasElement, healpixGrid: HealpixGrid, webgl: WebGL2RenderingContext): Point[];
+    /** Returns the center point (in J2000) of the current view as a `Point`. */
+    static getCenterJ2000(canvas: HTMLCanvasElement, healpixGrid: HealpixGrid, webgl: WebGL2RenderingContext, camera: Camera, pMatrix: ReadonlyMat4): Point;
+    /** Middle point on the unit sphere along the arc between two 3D points. */
+    static computeMiddlePoint(p1: Point, p2: Point): Point[];
+    /**
+     * Nearest intersection point between a frustum plane and the unit sphere,
+     * using the plane normal.
+     */
+    static getNearestSpherePoint(plane: number[]): Point[];
+    /**
+     * Intersections between a frustum plane and the unit sphere,
+     * computed via two perpendicular planes.
+     * Returns two points (first from `plane4Circle_1`, second from `plane4Circle_2`).
+     */
+    static getFrustumIntersectionWithSphere(_M: ReadonlyMat4, plane4Sphere: number[], plane4Circle_1: number[], plane4Circle_2: number[]): Point[];
+    /** Build ADQL string from an array of Points (ra,dec pairs). */
+    static getAstroFoVPolygon(points: Point[]): string;
+}
+//# sourceMappingURL=FoVUtils.d.ts.map

package/lib-esm/utils/FoVUtils.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"FoVUtils.d.ts","sourceRoot":"","sources":["../../src/utils/FoVUtils.ts"],"names":[],"mappings":"AACA;;GAEG;AAEH,OAAO,EAAE,KAAK,EAAE,MAAM,mBAAmB,CAAC;AAG1C,OAAO,EAAQ,YAAY,EAAE,MAAM,WAAW,CAAC;AAE/C,OAAO,MAAM,MAAM,cAAc,CAAC;AAClC,OAAO,EAAE,iBAAiB,EAAE,MAAM,+BAA+B,CAAC;AAClE,OAAO,EAAE,WAAW,EAAE,MAAM,8BAA8B,CAAC;AAE3D,qBAAa,QAAQ;IACnB;;;OAGG;IACH,SAAS,CAAC,IAAI,EAAE;QAAE,SAAS,EAAE,MAAM,CAAC;QAAC,SAAS,EAAE,MAAM,CAAA;KAAE,GAAG,MAAM;IAIjE;;;OAGG;IACH,MAAM,CAAC,aAAa,CAElB,MAAM,EAAE,MAAM,EACd,MAAM,EAAE,iBAAiB,EACzB,KAAK,EAAE,iBAAiB,EACxB,WAAW,EAAE,WAAW,EACxB,KAAK,EAAE,sBAAsB,EAC7B,OAAO,EAAE,YAAY,GACpB,KAAK,EAAE;IA+KV;;;OAGG;IACH,MAAM,CAAC,4BAA4B,CACjC,OAAO,EAAE,YAAY,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,iBAAiB,EAAE,WAAW,EAAE,WAAW,EAAE,KAAK,EAAE,sBAAsB,GAAG,KAAK,EAAE;IAmCrI,4EAA4E;IAC5E,MAAM,CAAC,cAAc,CAAC,MAAM,EAAE,iBAAiB,EAC7C,WAAW,EAAE,WAAW,EACxB,KAAK,EAAE,sBAAsB,EAC7B,MAAM,EAAE,MAAM,EAAE,OAAO,EAAE,YAAY,GAAG,KAAK;IAY/C,2EAA2E;IAC3E,MAAM,CAAC,kBAAkB,CAAC,EAAE,EAAE,KAAK,EAAE,EAAE,EAAE,KAAK,GAAG,KAAK,EAAE;IAexD;;;OAGG;IACH,MAAM,CAAC,qBAAqB,CAAC,KAAK,EAAE,MAAM,EAAE,GAAG,KAAK,EAAE;IAmBtD;;;;OAIG;IACH,MAAM,CAAC,gCAAgC,CACrC,EAAE,EAAE,YAAY,EAChB,YAAY,EAAE,MAAM,EAAE,EACtB,cAAc,EAAE,MAAM,EAAE,EACxB,cAAc,EAAE,MAAM,EAAE,GACvB,KAAK,EAAE;IAsDV,gEAAgE;IAChE,MAAM,CAAC,kBAAkB,CAAC,MAAM,EAAE,KAAK,EAAE,GAAG,MAAM;CAGnD"}

package/lib-esm/utils/FoVUtils.js

@@ -0,0 +1,217 @@
+'use strict';
+/**
+ * @author Fabrizio Giordano (Fab)
+ */
+import { Point } from '../model/Point.js';
+import RayPickingUtils from './RayPickingUtils.js';
+import { CoordsType } from './CoordsType.js';
+import { mat4 } from 'gl-matrix';
+export class FoVUtils {
+    /**
+     * Return the minimum FoV value between `_fovY_deg` and `_fovX_deg`.
+     * (Kept here for parity; this class doesn’t maintain those fields.)
+     */
+    getMinFoV() {
+        return this._fovY_deg <= this._fovX_deg ? this._fovY_deg : this._fovX_deg;
+    }
+    /**
+     * Compute the FoV polygon as a list of Points (clockwise).
+     * Uses ray picking + frustum planes against a unit sphere.
+     */
+    static getFoVPolygon(
+    // _pMatrix: ReadonlyMat4 | null,
+    camera, canvas, model, healpixGrid, webgl, pMatrix) {
+        // const pMatrix = (computePerspectiveMatrixSingleton.pMatrix ??
+        //   _pMatrix) as ReadonlyMat4;
+        // const pMatrix = computePerspectiveMatrixSingleton.pMatrix as ReadonlyMat4 
+        const vMatrix = camera.getCameraMatrix();
+        const mMatrix = model.getModelMatrix();
+        const canvasWidth = canvas.clientWidth;
+        const canvasHeight = canvas.clientHeight;
+        let points = [];
+        // First check: does the sphere cover the whole screen?
+        const intersectionWithModel = RayPickingUtils.getIntersectionPointWithSingleModel(0, 0, healpixGrid, webgl, camera, pMatrix);
+        if (intersectionWithModel.length > 0) {
+            // Fully covered → grab corners + midpoints (CASE C)
+            const cornersPoints = FoVUtils.getScreenCornersIntersection(pMatrix, camera, canvas, healpixGrid, webgl);
+            points = cornersPoints;
+        }
+        else {
+            // Partial coverage: build frustum planes
+            let M = mat4.create();
+            M = mat4.multiply(M, vMatrix, mMatrix);
+            M = mat4.multiply(M, pMatrix, M);
+            const topPlane = [M[3] - M[1], M[7] - M[5], M[11] - M[9], M[15] - M[13]]; // m41-m21, ...
+            const bottomPlane = [M[3] + M[1], M[7] + M[5], M[11] + M[9], M[15] + M[13]];
+            const rightPlane = [M[3] - M[0], M[7] - M[4], M[11] - M[8], M[15] - M[12]];
+            const leftPlane = [M[3] + M[0], M[7] + M[4], M[11] + M[8], M[15] + M[12]];
+            const intersectionTopMiddle = RayPickingUtils.getIntersectionPointWithSingleModel(canvasWidth / 2, 0, healpixGrid, webgl, camera, pMatrix);
+            const intersectionRightMiddle = RayPickingUtils.getIntersectionPointWithSingleModel(canvasWidth, canvasHeight / 2, healpixGrid, webgl, camera, pMatrix);
+            // CASE A: zoomed out, hemisphere fully visible
+            if (intersectionTopMiddle.length === 0 &&
+                intersectionRightMiddle.length === 0) {
+                const topPoints = FoVUtils.getNearestSpherePoint(topPlane);
+                const bottomPoints = FoVUtils.getNearestSpherePoint(bottomPlane);
+                const leftPoints = FoVUtils.getNearestSpherePoint(leftPlane);
+                const rightPoints = FoVUtils.getNearestSpherePoint(rightPlane);
+                const middleLeftTop = FoVUtils.computeMiddlePoint(leftPoints[0], topPoints[0])[0];
+                const middleTopRight = FoVUtils.computeMiddlePoint(topPoints[0], rightPoints[0])[0];
+                const middleRightBottom = FoVUtils.computeMiddlePoint(rightPoints[0], bottomPoints[0])[0];
+                const middleBottomLeft = FoVUtils.computeMiddlePoint(bottomPoints[0], leftPoints[0])[0];
+                points.push(topPoints[0], middleTopRight, rightPoints[0], middleRightBottom, bottomPoints[0], middleBottomLeft, leftPoints[0], middleLeftTop);
+            }
+            // CASE E: no intersection on top/bottom planes
+            else if (intersectionTopMiddle.length === 0) {
+                const topPoints = FoVUtils.getNearestSpherePoint(topPlane);
+                const bottomPoints = FoVUtils.getNearestSpherePoint(bottomPlane);
+                const leftPoints = FoVUtils.getFrustumIntersectionWithSphere(M, leftPlane, bottomPlane, topPlane);
+                const rightPoints = FoVUtils.getFrustumIntersectionWithSphere(M, rightPlane, topPlane, bottomPlane);
+                const middleLeftTop = FoVUtils.computeMiddlePoint(leftPoints[1], topPoints[0])[0];
+                const middleTopRight = FoVUtils.computeMiddlePoint(topPoints[0], rightPoints[0])[0];
+                const middleRightBottom = FoVUtils.computeMiddlePoint(rightPoints[1], bottomPoints[0])[0];
+                const middleBottomLeft = FoVUtils.computeMiddlePoint(bottomPoints[0], leftPoints[0])[0];
+                points.push(topPoints[0], middleTopRight, rightPoints[0], rightPoints[1], middleRightBottom, bottomPoints[0], middleBottomLeft, leftPoints[0], leftPoints[1], middleLeftTop);
+            }
+            // CASE D: no intersection on right/left planes
+            else if (intersectionRightMiddle.length === 0) {
+                const topPoints = FoVUtils.getFrustumIntersectionWithSphere(M, topPlane, leftPlane, rightPlane);
+                const bottomPoints = FoVUtils.getFrustumIntersectionWithSphere(M, bottomPlane, rightPlane, leftPlane);
+                const leftPoints = FoVUtils.getNearestSpherePoint(leftPlane);
+                const rightPoints = FoVUtils.getNearestSpherePoint(rightPlane);
+                const middleLeftTop = FoVUtils.computeMiddlePoint(leftPoints[0], topPoints[0])[0];
+                const middleTopRight = FoVUtils.computeMiddlePoint(topPoints[1], rightPoints[0])[0];
+                const middleRightBottom = FoVUtils.computeMiddlePoint(rightPoints[0], bottomPoints[0])[0];
+                const middleBottomLeft = FoVUtils.computeMiddlePoint(bottomPoints[1], leftPoints[0])[0];
+                points.push(topPoints[0], topPoints[1], middleTopRight, rightPoints[0], middleRightBottom, bottomPoints[0], bottomPoints[1], middleBottomLeft, leftPoints[0], middleLeftTop);
+            }
+            // CASE B: all frustum planes intersect
+            else {
+                const topPoints = FoVUtils.getFrustumIntersectionWithSphere(M, topPlane, leftPlane, rightPlane);
+                const bottomPoints = FoVUtils.getFrustumIntersectionWithSphere(M, bottomPlane, rightPlane, leftPlane);
+                const leftPoints = FoVUtils.getFrustumIntersectionWithSphere(M, leftPlane, bottomPlane, topPlane);
+                const rightPoints = FoVUtils.getFrustumIntersectionWithSphere(M, rightPlane, topPlane, bottomPlane);
+                points.push(topPoints[0], topPoints[1], rightPoints[0], rightPoints[1], bottomPoints[0], bottomPoints[1], leftPoints[0], leftPoints[1]);
+            }
+        }
+        return points;
+    }
+    /**
+     * Ray pick against 8 key screen positions (corners + midpoints).
+     * Returns Points in clockwise order starting from top-left.
+     */
+    static getScreenCornersIntersection(pMatrix, camera, canvas, healpixGrid, webgl) {
+        const w = canvas.clientWidth;
+        const h = canvas.clientHeight;
+        const topLeft = RayPickingUtils.getIntersectionPointWithSingleModel(0, 0, healpixGrid, webgl, camera, pMatrix);
+        const middleTop = RayPickingUtils.getIntersectionPointWithSingleModel(w / 2, 0, healpixGrid, webgl, camera, pMatrix);
+        const topRight = RayPickingUtils.getIntersectionPointWithSingleModel(w, 0, healpixGrid, webgl, camera, pMatrix);
+        const middleRight = RayPickingUtils.getIntersectionPointWithSingleModel(w, h / 2, healpixGrid, webgl, camera, pMatrix);
+        const bottomRight = RayPickingUtils.getIntersectionPointWithSingleModel(w, h, healpixGrid, webgl, camera, pMatrix);
+        const middleBottom = RayPickingUtils.getIntersectionPointWithSingleModel(w / 2, h, healpixGrid, webgl, camera, pMatrix);
+        const bottomLeft = RayPickingUtils.getIntersectionPointWithSingleModel(0, h, healpixGrid, webgl, camera, pMatrix);
+        const middleLeft = RayPickingUtils.getIntersectionPointWithSingleModel(0, h / 2, healpixGrid, webgl, camera, pMatrix);
+        const out = [];
+        const pushIf = (ip) => {
+            if (ip.length > 0) {
+                out.push(new Point({ x: ip[0], y: ip[1], z: ip[2] }, CoordsType.CARTESIAN));
+            }
+        };
+        pushIf(topLeft);
+        pushIf(middleTop);
+        pushIf(topRight);
+        pushIf(middleRight);
+        pushIf(bottomRight);
+        pushIf(middleBottom);
+        pushIf(bottomLeft);
+        pushIf(middleLeft);
+        return out;
+    }
+    /** Returns the center point (in J2000) of the current view as a `Point`. */
+    static getCenterJ2000(canvas, healpixGrid, webgl, camera, pMatrix) {
+        const w = canvas.clientWidth;
+        const h = canvas.clientHeight;
+        const center = RayPickingUtils.getIntersectionPointWithSingleModel(w / 2, h / 2, healpixGrid, webgl, camera, pMatrix);
+        if (center.length <= 0)
+            throw Error(`Central point is null`);
+        return new Point({ x: center[0], y: center[1], z: center[2] }, CoordsType.CARTESIAN);
+    }
+    /** Middle point on the unit sphere along the arc between two 3D points. */
+    static computeMiddlePoint(p1, p2) {
+        // midpoint of segment
+        const xm = (p1.x + p2.x) / 2;
+        const ym = (p1.y + p2.y) / 2;
+        const zm = (p1.z + p2.z) / 2;
+        // project the midpoint back to unit sphere
+        const len = Math.hypot(xm, ym, zm) || 1;
+        const x = xm / len;
+        const y = ym / len;
+        const z = zm / len;
+        return [new Point({ x, y, z }, CoordsType.CARTESIAN)];
+    }
+    /**
+     * Nearest intersection point between a frustum plane and the unit sphere,
+     * using the plane normal.
+     */
+    static getNearestSpherePoint(plane) {
+        const [A, B, C, D] = plane;
+        const R = 1;
+        const invLen = 1 / Math.sqrt(A * A + B * B + C * C);
+        const t1 = R * invLen;
+        const t2 = -R * invLen;
+        const P1 = [A * t1, B * t1, C * t1];
+        const P2 = [A * t2, B * t2, C * t2];
+        const den = Math.sqrt(A * A + B * B + C * C) || 1;
+        const dist1 = Math.abs(A * P1[0] + B * P1[1] + C * P1[2] + D) / den;
+        const dist2 = Math.abs(A * P2[0] + B * P2[1] + C * P2[2] + D) / den;
+        const P = dist1 <= dist2 ? P1 : P2;
+        return [new Point({ x: P[0], y: P[1], z: P[2] }, CoordsType.CARTESIAN)];
+    }
+    /**
+     * Intersections between a frustum plane and the unit sphere,
+     * computed via two perpendicular planes.
+     * Returns two points (first from `plane4Circle_1`, second from `plane4Circle_2`).
+     */
+    static getFrustumIntersectionWithSphere(_M, plane4Sphere, plane4Circle_1, plane4Circle_2) {
+        const [A0, B0, C0, D0] = plane4Sphere;
+        // center of the circle (projection of sphere center onto plane)
+        const denom0 = (A0 * A0 + B0 * B0 + C0 * C0) || 1;
+        const x_c = -(A0 * D0) / denom0;
+        const y_c = -(B0 * D0) / denom0;
+        const z_c = -(C0 * D0) / denom0;
+        const d = Math.abs(D0) / Math.sqrt(denom0); // distance from sphere center (0,0,0)
+        const R = 1;
+        const out = [];
+        if (R > d) {
+            const r = Math.sqrt(R * R - d * d);
+            const pick = (plane) => {
+                const [A, B, C, D] = plane;
+                const invLen = 1 / Math.sqrt(A * A + B * B + C * C);
+                const t1 = r * invLen;
+                const t2 = -r * invLen;
+                const P1 = [x_c + A * t1, y_c + B * t1, z_c + C * t1];
+                const P2 = [x_c + A * t2, y_c + B * t2, z_c + C * t2];
+                const den = Math.sqrt(A * A + B * B + C * C) || 1;
+                const dist1 = Math.abs(A * P1[0] + B * P1[1] + C * P1[2] + D) / den;
+                const dist2 = Math.abs(A * P2[0] + B * P2[1] + C * P2[2] + D) / den;
+                return dist1 <= dist2 ? P1 : P2;
+            };
+            const P_intersection_1 = pick(plane4Circle_1);
+            const P_intersection_2 = pick(plane4Circle_2);
+            out.push(new Point({ x: P_intersection_1[0], y: P_intersection_1[1], z: P_intersection_1[2] }, CoordsType.CARTESIAN), new Point({ x: P_intersection_2[0], y: P_intersection_2[1], z: P_intersection_2[2] }, CoordsType.CARTESIAN));
+        }
+        else if (R === d) {
+            // Tangent: both intersections collapse to the circle center on the plane
+            out.push(new Point({ x: x_c, y: y_c, z: z_c }, CoordsType.CARTESIAN), new Point({ x: x_c, y: y_c, z: z_c }, CoordsType.CARTESIAN));
+        }
+        else {
+            // No intersection; return empty to avoid pushing undefined values
+            // console.log('Frustum plane not intersecting the sphere');
+        }
+        return out;
+    }
+    /** Build ADQL string from an array of Points (ra,dec pairs). */
+    static getAstroFoVPolygon(points) {
+        return points.map(p => p.toADQL()).join(',');
+    }
+}
+//# sourceMappingURL=FoVUtils.js.map

package/lib-esm/utils/FoVUtils.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"FoVUtils.js","sourceRoot":"","sources":["../../src/utils/FoVUtils.ts"],"names":[],"mappings":"AAAA,YAAY,CAAC;AACb;;GAEG;AAEH,OAAO,EAAE,KAAK,EAAE,MAAM,mBAAmB,CAAC;AAC1C,OAAO,eAAe,MAAM,sBAAsB,CAAC;AACnD,OAAO,EAAE,UAAU,EAAE,MAAM,iBAAiB,CAAC;AAC7C,OAAO,EAAE,IAAI,EAAgB,MAAM,WAAW,CAAC;AAM/C,MAAM,OAAO,QAAQ;IACnB;;;OAGG;IACH,SAAS;QACP,OAAO,IAAI,CAAC,SAAS,IAAI,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,SAAS,CAAC,CAAC,CAAC,IAAI,CAAC,SAAS,CAAC;IAC5E,CAAC;IAED;;;OAGG;IACH,MAAM,CAAC,aAAa;IAClB,iCAAiC;IACjC,MAAc,EACd,MAAyB,EACzB,KAAwB,EACxB,WAAwB,EACxB,KAA6B,EAC7B,OAAqB;QAErB,gEAAgE;QAChE,+BAA+B;QAC/B,6EAA6E;QAC7E,MAAM,OAAO,GAAG,MAAM,CAAC,eAAe,EAAE,CAAC;QACzC,MAAM,OAAO,GAAG,KAAK,CAAC,cAAc,EAAE,CAAC;QACvC,MAAM,WAAW,GAAG,MAAM,CAAC,WAAW,CAAC;QACvC,MAAM,YAAY,GAAG,MAAM,CAAC,YAAY,CAAC;QAEzC,IAAI,MAAM,GAAY,EAAE,CAAC;QAEzB,uDAAuD;QACvD,MAAM,qBAAqB,GAAG,eAAe,CAAC,mCAAmC,CAAC,CAAC,EAAE,CAAC,EAAE,WAAW,EAAE,KAAK,EAAE,MAAM,EAAE,OAAO,CAAC,CAAC;QAE7H,IAAI,qBAAqB,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;YACrC,oDAAoD;YACpD,MAAM,aAAa,GAAG,QAAQ,CAAC,4BAA4B,CACzD,OAAO,EACP,MAAM,EACN,MAAM,EAAE,WAAW,EAAE,KAAK,CAC3B,CAAC;YACF,MAAM,GAAG,aAAa,CAAC;QACzB,CAAC;aAAM,CAAC;YACN,yCAAyC;YACzC,IAAI,CAAC,GAAG,IAAI,CAAC,MAAM,EAAE,CAAC;YACtB,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,EAAE,OAAO,EAAE,OAAO,CAAC,CAAC;YACvC,CAAC,GAAG,IAAI,CAAC,QAAQ,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,CAAC,CAAC;YAEjC,MAAM,QAAQ,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAI,eAAe;YAC5F,MAAM,WAAW,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YAC5E,MAAM,UAAU,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YAC3E,MAAM,SAAS,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YAE1E,MAAM,qBAAqB,GAAG,eAAe,CAAC,mCAAmC,CAAC,WAAW,GAAG,CAAC,EAAE,CAAC,EAAE,WAAW,EAAE,KAAK,EAAE,MAAM,EAAE,OAAO,CAAC,CAAC;YAC3I,MAAM,uBAAuB,GAAG,eAAe,CAAC,mCAAmC,CAAC,WAAW,EAAE,YAAY,GAAG,CAAC,EAAE,WAAW,EAAE,KAAK,EAAE,MAAM,EAAE,OAAO,CAAC,CAAC;YAExJ,+CAA+C;YAC/C,IACE,qBAAqB,CAAC,MAAM,KAAK,CAAC;gBAClC,uBAAuB,CAAC,MAAM,KAAK,CAAC,EACpC,CAAC;gBACD,MAAM,SAAS,GAAG,QAAQ,CAAC,qBAAqB,CAAC,QAAQ,CAAC,CAAC;gBAC3D,MAAM,YAAY,GAAG,QAAQ,CAAC,qBAAqB,CAAC,WAAW,CAAC,CAAC;gBACjE,MAAM,UAAU,GAAG,QAAQ,CAAC,qBAAqB,CAAC,SAAS,CAAC,CAAC;gBAC7D,MAAM,WAAW,GAAG,QAAQ,CAAC,qBAAqB,CAAC,UAAU,CAAC,CAAC;gBAE/D,MAAM,aAAa,GAAG,QAAQ,CAAC,kBAAkB,CAAC,UAAU,CAAC,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAClF,MAAM,cAAc,GAAG,QAAQ,CAAC,kBAAkB,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACpF,MAAM,iBAAiB,GAAG,QAAQ,CAAC,kBAAkB,CAAC,WAAW,CAAC,CAAC,CAAC,EAAE,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC1F,MAAM,gBAAgB,GAAG,QAAQ,CAAC,kBAAkB,CAAC,YAAY,CAAC,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAExF,MAAM,CAAC,IAAI,CACT,SAAS,CAAC,CAAC,CAAC,EACZ,cAAc,EACd,WAAW,CAAC,CAAC,CAAC,EACd,iBAAiB,EACjB,YAAY,CAAC,CAAC,CAAC,EACf,gBAAgB,EAChB,UAAU,CAAC,CAAC,CAAC,EACb,aAAa,CACd,CAAC;YACJ,CAAC;YACD,+CAA+C;iBAC1C,IAAI,qBAAqB,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;gBAC5C,MAAM,SAAS,GAAG,QAAQ,CAAC,qBAAqB,CAAC,QAAQ,CAAC,CAAC;gBAC3D,MAAM,YAAY,GAAG,QAAQ,CAAC,qBAAqB,CAAC,WAAW,CAAC,CAAC;gBACjE,MAAM,UAAU,GAAG,QAAQ,CAAC,gCAAgC,CAC1D,CAAC,EACD,SAAS,EACT,WAAW,EACX,QAAQ,CACT,CAAC;gBACF,MAAM,WAAW,GAAG,QAAQ,CAAC,gCAAgC,CAC3D,CAAC,EACD,UAAU,EACV,QAAQ,EACR,WAAW,CACZ,CAAC;gBAEF,MAAM,aAAa,GAAG,QAAQ,CAAC,kBAAkB,CAAC,UAAU,CAAC,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAClF,MAAM,cAAc,GAAG,QAAQ,CAAC,kBAAkB,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACpF,MAAM,iBAAiB,GAAG,QAAQ,CAAC,kBAAkB,CAAC,WAAW,CAAC,CAAC,CAAC,EAAE,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC1F,MAAM,gBAAgB,GAAG,QAAQ,CAAC,kBAAkB,CAAC,YAAY,CAAC,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAExF,MAAM,CAAC,IAAI,CACT,SAAS,CAAC,CAAC,CAAC,EACZ,cAAc,EACd,WAAW,CAAC,CAAC,CAAC,EACd,WAAW,CAAC,CAAC,CAAC,EACd,iBAAiB,EACjB,YAAY,CAAC,CAAC,CAAC,EACf,gBAAgB,EAChB,UAAU,CAAC,CAAC,CAAC,EACb,UAAU,CAAC,CAAC,CAAC,EACb,aAAa,CACd,CAAC;YACJ,CAAC;YACD,+CAA+C;iBAC1C,IAAI,uBAAuB,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;gBAC9C,MAAM,SAAS,GAAG,QAAQ,CAAC,gCAAgC,CACzD,CAAC,EACD,QAAQ,EACR,SAAS,EACT,UAAU,CACX,CAAC;gBACF,MAAM,YAAY,GAAG,QAAQ,CAAC,gCAAgC,CAC5D,CAAC,EACD,WAAW,EACX,UAAU,EACV,SAAS,CACV,CAAC;gBACF,MAAM,UAAU,GAAG,QAAQ,CAAC,qBAAqB,CAAC,SAAS,CAAC,CAAC;gBAC7D,MAAM,WAAW,GAAG,QAAQ,CAAC,qBAAqB,CAAC,UAAU,CAAC,CAAC;gBAE/D,MAAM,aAAa,GAAG,QAAQ,CAAC,kBAAkB,CAAC,UAAU,CAAC,CAAC,CAAC,EAAE,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAClF,MAAM,cAAc,GAAG,QAAQ,CAAC,kBAAkB,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACpF,MAAM,iBAAiB,GAAG,QAAQ,CAAC,kBAAkB,CAAC,WAAW,CAAC,CAAC,CAAC,EAAE,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC1F,MAAM,gBAAgB,GAAG,QAAQ,CAAC,kBAAkB,CAAC,YAAY,CAAC,CAAC,CAAC,EAAE,UAAU,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAExF,MAAM,CAAC,IAAI,CACT,SAAS,CAAC,CAAC,CAAC,EACZ,SAAS,CAAC,CAAC,CAAC,EACZ,cAAc,EACd,WAAW,CAAC,CAAC,CAAC,EACd,iBAAiB,EACjB,YAAY,CAAC,CAAC,CAAC,EACf,YAAY,CAAC,CAAC,CAAC,EACf,gBAAgB,EAChB,UAAU,CAAC,CAAC,CAAC,EACb,aAAa,CACd,CAAC;YACJ,CAAC;YACD,uCAAuC;iBAClC,CAAC;gBACJ,MAAM,SAAS,GAAG,QAAQ,CAAC,gCAAgC,CACzD,CAAC,EACD,QAAQ,EACR,SAAS,EACT,UAAU,CACX,CAAC;gBACF,MAAM,YAAY,GAAG,QAAQ,CAAC,gCAAgC,CAC5D,CAAC,EACD,WAAW,EACX,UAAU,EACV,SAAS,CACV,CAAC;gBACF,MAAM,UAAU,GAAG,QAAQ,CAAC,gCAAgC,CAC1D,CAAC,EACD,SAAS,EACT,WAAW,EACX,QAAQ,CACT,CAAC;gBACF,MAAM,WAAW,GAAG,QAAQ,CAAC,gCAAgC,CAC3D,CAAC,EACD,UAAU,EACV,QAAQ,EACR,WAAW,CACZ,CAAC;gBAEF,MAAM,CAAC,IAAI,CACT,SAAS,CAAC,CAAC,CAAC,EACZ,SAAS,CAAC,CAAC,CAAC,EACZ,WAAW,CAAC,CAAC,CAAC,EACd,WAAW,CAAC,CAAC,CAAC,EACd,YAAY,CAAC,CAAC,CAAC,EACf,YAAY,CAAC,CAAC,CAAC,EACf,UAAU,CAAC,CAAC,CAAC,EACb,UAAU,CAAC,CAAC,CAAC,CACd,CAAC;YACJ,CAAC;QACH,CAAC;QAED,OAAO,MAAM,CAAC;IAChB,CAAC;IAED;;;OAGG;IACH,MAAM,CAAC,4BAA4B,CACjC,OAAqB,EAAE,MAAc,EAAE,MAAyB,EAAE,WAAwB,EAAE,KAA6B;QACzH,MAAM,CAAC,GAAG,MAAM,CAAC,WAAW,CAAC;QAC7B,MAAM,CAAC,GAAG,MAAM,CAAC,YAAY,CAAC;QAE9B,MAAM,OAAO,GAAG,eAAe,CAAC,mCAAmC,CAAC,CAAC,EAAE,CAAC,EAAE,WAAW,EAAE,KAAK,EAAE,MAAM,EAAE,OAAO,CAAC,CAAC;QAC/G,MAAM,SAAS,GAAG,eAAe,CAAC,mCAAmC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,WAAW,EAAE,KAAK,EAAE,MAAM,EAAE,OAAO,CAAC,CAAC;QACrH,MAAM,QAAQ,GAAG,eAAe,CAAC,mCAAmC,CAAC,CAAC,EAAE,CAAC,EAAE,WAAW,EAAE,KAAK,EAAE,MAAM,EAAE,OAAO,CAAC,CAAC;QAEhH,MAAM,WAAW,GAAG,eAAe,CAAC,mCAAmC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,WAAW,EAAE,KAAK,EAAE,MAAM,EAAE,OAAO,CAAC,CAAC;QAEvH,MAAM,WAAW,GAAG,eAAe,CAAC,mCAAmC,CAAC,CAAC,EAAE,CAAC,EAAE,WAAW,EAAE,KAAK,EAAE,MAAM,EAAE,OAAO,CAAC,CAAC;QACnH,MAAM,YAAY,GAAG,eAAe,CAAC,mCAAmC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,WAAW,EAAE,KAAK,EAAE,MAAM,EAAE,OAAO,CAAC,CAAC;QACxH,MAAM,UAAU,GAAG,eAAe,CAAC,mCAAmC,CAAC,CAAC,EAAE,CAAC,EAAE,WAAW,EAAE,KAAK,EAAE,MAAM,EAAE,OAAO,CAAC,CAAC;QAElH,MAAM,UAAU,GAAG,eAAe,CAAC,mCAAmC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,WAAW,EAAE,KAAK,EAAE,MAAM,EAAE,OAAO,CAAC,CAAC;QAEtH,MAAM,GAAG,GAAY,EAAE,CAAC;QACxB,MAAM,MAAM,GAAG,CAAC,EAAY,EAAE,EAAE;YAC9B,IAAI,EAAE,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;gBAClB,GAAG,CAAC,IAAI,CAAC,IAAI,KAAK,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,UAAU,CAAC,SAAS,CAAC,CAAC,CAAC;YAC9E,CAAC;QACH,CAAC,CAAC;QAEF,MAAM,CAAC,OAAO,CAAC,CAAC;QAChB,MAAM,CAAC,SAAS,CAAC,CAAC;QAClB,MAAM,CAAC,QAAQ,CAAC,CAAC;QACjB,MAAM,CAAC,WAAW,CAAC,CAAC;QACpB,MAAM,CAAC,WAAW,CAAC,CAAC;QACpB,MAAM,CAAC,YAAY,CAAC,CAAC;QACrB,MAAM,CAAC,UAAU,CAAC,CAAC;QACnB,MAAM,CAAC,UAAU,CAAC,CAAC;QAEnB,OAAO,GAAG,CAAC;IACb,CAAC;IAED,4EAA4E;IAC5E,MAAM,CAAC,cAAc,CAAC,MAAyB,EAC7C,WAAwB,EACxB,KAA6B,EAC7B,MAAc,EAAE,OAAqB;QACrC,MAAM,CAAC,GAAG,MAAM,CAAC,WAAW,CAAC;QAC7B,MAAM,CAAC,GAAG,MAAM,CAAC,YAAY,CAAC;QAE9B,MAAM,MAAM,GAAG,eAAe,CAAC,mCAAmC,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,WAAW,EAAE,KAAK,EAAE,MAAM,EAAE,OAAO,CAAC,CAAC;QACtH,IAAI,MAAM,CAAC,MAAM,IAAI,CAAC;YAAE,MAAM,KAAK,CAAC,uBAAuB,CAAC,CAAA;QAC5D,OAAO,IAAI,KAAK,CACd,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,EAC5C,UAAU,CAAC,SAAS,CACrB,CAAC;IACJ,CAAC;IAED,2EAA2E;IAC3E,MAAM,CAAC,kBAAkB,CAAC,EAAS,EAAE,EAAS;QAC5C,sBAAsB;QACtB,MAAM,EAAE,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QAC7B,MAAM,EAAE,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QAC7B,MAAM,EAAE,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QAE7B,2CAA2C;QAC3C,MAAM,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,IAAI,CAAC,CAAC;QACxC,MAAM,CAAC,GAAG,EAAE,GAAG,GAAG,CAAC;QACnB,MAAM,CAAC,GAAG,EAAE,GAAG,GAAG,CAAC;QACnB,MAAM,CAAC,GAAG,EAAE,GAAG,GAAG,CAAC;QAEnB,OAAO,CAAC,IAAI,KAAK,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,UAAU,CAAC,SAAS,CAAC,CAAC,CAAC;IACxD,CAAC;IAED;;;OAGG;IACH,MAAM,CAAC,qBAAqB,CAAC,KAAe;QAC1C,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC;QAE3B,MAAM,CAAC,GAAG,CAAC,CAAC;QACZ,MAAM,MAAM,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QACpD,MAAM,EAAE,GAAG,CAAC,GAAG,MAAM,CAAC;QACtB,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC;QAEvB,MAAM,EAAE,GAA6B,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC;QAC9D,MAAM,EAAE,GAA6B,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,EAAE,CAAC,CAAC;QAE9D,MAAM,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC;QAClD,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;QACpE,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;QAEpE,MAAM,CAAC,GAAG,KAAK,IAAI,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;QACnC,OAAO,CAAC,IAAI,KAAK,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,EAAE,UAAU,CAAC,SAAS,CAAC,CAAC,CAAC;IAC1E,CAAC;IAED;;;;OAIG;IACH,MAAM,CAAC,gCAAgC,CACrC,EAAgB,EAChB,YAAsB,EACtB,cAAwB,EACxB,cAAwB;QAExB,MAAM,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,GAAG,YAAY,CAAC;QAEtC,gEAAgE;QAChE,MAAM,MAAM,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,IAAI,CAAC,CAAC;QAClD,MAAM,GAAG,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,MAAM,CAAC;QAChC,MAAM,GAAG,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,MAAM,CAAC;QAChC,MAAM,GAAG,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,MAAM,CAAC;QAEhC,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,CAAC,CAAC,CAAC,sCAAsC;QAClF,MAAM,CAAC,GAAG,CAAC,CAAC;QAEZ,MAAM,GAAG,GAAY,EAAE,CAAC;QAExB,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;YACV,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;YAEnC,MAAM,IAAI,GAAG,CAAC,KAAe,EAA4B,EAAE;gBACzD,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,GAAG,KAAK,CAAC;gBAC3B,MAAM,MAAM,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;gBACpD,MAAM,EAAE,GAAG,CAAC,GAAG,MAAM,CAAC;gBACtB,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC;gBAEvB,MAAM,EAAE,GAA6B,CAAC,GAAG,GAAG,CAAC,GAAG,EAAE,EAAE,GAAG,GAAG,CAAC,GAAG,EAAE,EAAE,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;gBAChF,MAAM,EAAE,GAA6B,CAAC,GAAG,GAAG,CAAC,GAAG,EAAE,EAAE,GAAG,GAAG,CAAC,GAAG,EAAE,EAAE,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC;gBAEhF,MAAM,GAAG,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC;gBAClD,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;gBACpE,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC;gBAEpE,OAAO,KAAK,IAAI,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC;YAClC,CAAC,CAAC;YAEF,MAAM,gBAAgB,GAAG,IAAI,CAAC,cAAc,CAAC,CAAC;YAC9C,MAAM,gBAAgB,GAAG,IAAI,CAAC,cAAc,CAAC,CAAC;YAE9C,GAAG,CAAC,IAAI,CACN,IAAI,KAAK,CAAC,EAAE,CAAC,EAAE,gBAAgB,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,gBAAgB,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,gBAAgB,CAAC,CAAC,CAAC,EAAE,EAAE,UAAU,CAAC,SAAS,CAAC,EAC3G,IAAI,KAAK,CAAC,EAAE,CAAC,EAAE,gBAAgB,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,gBAAgB,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,gBAAgB,CAAC,CAAC,CAAC,EAAE,EAAE,UAAU,CAAC,SAAS,CAAC,CAC5G,CAAC;QACJ,CAAC;aAAM,IAAI,CAAC,KAAK,CAAC,EAAE,CAAC;YACnB,yEAAyE;YACzE,GAAG,CAAC,IAAI,CACN,IAAI,KAAK,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,EAAE,UAAU,CAAC,SAAS,CAAC,EAC3D,IAAI,KAAK,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE,EAAE,UAAU,CAAC,SAAS,CAAC,CAC5D,CAAC;QACJ,CAAC;aAAM,CAAC;YACN,kEAAkE;YAClE,4DAA4D;QAC9D,CAAC;QAED,OAAO,GAAG,CAAC;IACb,CAAC;IAED,gEAAgE;IAChE,MAAM,CAAC,kBAAkB,CAAC,MAAe;QACvC,OAAO,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,MAAM,EAAE,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;IAC/C,CAAC;CACF"}

package/lib-esm/utils/GeomUtils.d.ts

@@ -0,0 +1,36 @@
+import { Point } from "../model/Point.js";
+import Point2D from "../model/Point2D.js";
+type ProjectionFlag = 0 | 1 | 2;
+export interface SelectionObj {
+    poly4selection: Point2D[][];
+    flag: ProjectionFlag;
+    maxx: number;
+    maxy: number;
+    minx: number;
+    miny: number;
+}
+declare class GeomUtils {
+    static orthodromicDistance(p1: Point, p2: Point): number;
+    /**
+     * Decide the 2D projection strategy and pre-project polygons for point-in-polygon tests.
+     * Returns the projected polygons + bbox + a flag describing the projection used:
+     * 0 → all points in same hemisphere with |Dec| > 10 → stereographic-like projection using x,y from 3D
+     * 1 → all points in equatorial belt (|Dec| < 10) → use RA/Dec directly
+     * 2 → equatorial belt and polygon crosses RA=0 → shift RA>180 by -360
+     */
+    static computeSelectionObject(polygons: Point[][]): SelectionObj;
+    /** Stereographic projection from 3D point on unit sphere onto plane */
+    static stereographic(point: Point): {
+        x: number;
+        y: number;
+    };
+    static projectIn2D(point: Point): Point2D;
+    /**
+     * Robust point-in-polygon (ray casting) using the precomputed selection object.
+     * Works with any of the three flags (0,1,2).
+     */
+    static checkPointInsidePolygon5(selectionObj: SelectionObj, point: Point): boolean;
+    static checkPointInsidePolygon4(polygons: Point[][], point: Point): boolean;
+}
+export default GeomUtils;
+//# sourceMappingURL=GeomUtils.d.ts.map

package/lib-esm/utils/GeomUtils.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"GeomUtils.d.ts","sourceRoot":"","sources":["../../src/utils/GeomUtils.ts"],"names":[],"mappings":"AAAA,OAAO,EAAC,KAAK,EAAC,MAAM,mBAAmB,CAAC;AACxC,OAAO,OAAO,MAAM,qBAAqB,CAAC;AAG1C,KAAK,cAAc,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;AAEhC,MAAM,WAAW,YAAY;IAC3B,cAAc,EAAE,OAAO,EAAE,EAAE,CAAC;IAC5B,IAAI,EAAE,cAAc,CAAC;IACrB,IAAI,EAAE,MAAM,CAAC;IACb,IAAI,EAAE,MAAM,CAAC;IACb,IAAI,EAAE,MAAM,CAAC;IACb,IAAI,EAAE,MAAM,CAAC;CACd;AAED,cAAM,SAAS;IAEb,MAAM,CAAC,mBAAmB,CAAC,EAAE,EAAE,KAAK,EAAE,EAAE,EAAE,KAAK,GAAG,MAAM;IAQxD;;;;;;OAMG;IACH,MAAM,CAAC,sBAAsB,CAAC,QAAQ,EAAE,KAAK,EAAE,EAAE,GAAG,YAAY;IAwIhE,uEAAuE;IACvE,MAAM,CAAC,aAAa,CAAC,KAAK,EAAE,KAAK,GAAG;QAAE,CAAC,EAAE,MAAM,CAAC;QAAC,CAAC,EAAE,MAAM,CAAA;KAAE;IAU5D,MAAM,CAAC,WAAW,CAAC,KAAK,EAAE,KAAK,GAAG,OAAO;IAKzC;;;OAGG;IACH,MAAM,CAAC,wBAAwB,CAAC,YAAY,EAAE,YAAY,EAAE,KAAK,EAAE,KAAK,GAAG,OAAO;IAyElF,MAAM,CAAC,wBAAwB,CAAC,QAAQ,EAAE,KAAK,EAAE,EAAE,EAAE,KAAK,EAAE,KAAK,GAAG,OAAO;CAmD5E;AAED,eAAe,SAAS,CAAC"}

package/lib-esm/utils/GeomUtils.js

@@ -0,0 +1,267 @@
+import { Point } from "../model/Point.js";
+import Point2D from "../model/Point2D.js";
+import { CoordsType } from "./CoordsType.js";
+class GeomUtils {
+    // Orthodromic (great-circle) distance in radians
+    static orthodromicDistance(p1, p2) {
+        return Math.acos(Math.sin(p1.decDeg * Math.PI / 180) * Math.sin(p2.decDeg * Math.PI / 180) +
+            Math.cos(p1.decDeg * Math.PI / 180) * Math.cos(p2.decDeg * Math.PI / 180) *
+                Math.cos((p2.raDeg - p1.raDeg) * Math.PI / 180));
+    }
+    /**
+     * Decide the 2D projection strategy and pre-project polygons for point-in-polygon tests.
+     * Returns the projected polygons + bbox + a flag describing the projection used:
+     * 0 → all points in same hemisphere with |Dec| > 10 → stereographic-like projection using x,y from 3D
+     * 1 → all points in equatorial belt (|Dec| < 10) → use RA/Dec directly
+     * 2 → equatorial belt and polygon crosses RA=0 → shift RA>180 by -360
+     */
+    static computeSelectionObject(polygons) {
+        let poly4selection = [];
+        let flag = 0;
+        let maxx;
+        let maxy;
+        let minx;
+        let miny;
+        const DEC_THRESHOLD = 10;
+        //  1 → northern hemisphere (Dec > +10), -1 → southern (Dec < -10), 0 → equatorial belt
+        let hemisphere = 0;
+        if (polygons[0][0].decDeg >= DEC_THRESHOLD) {
+            hemisphere = 1;
+        }
+        else if (polygons[0][0].decDeg <= -DEC_THRESHOLD) {
+            hemisphere = -1;
+        }
+        else {
+            flag = 1;
+        }
+        // Case flag = 0 → stereographic-like projection using x,y,z from 3D point
+        if (flag === 0) {
+            const first = GeomUtils.projectIn2D(polygons[0][0]);
+            maxx = minx = first.x;
+            maxy = miny = first.y;
+            for (const currpoly of polygons) {
+                const selpoly = [];
+                for (const point of currpoly) {
+                    // If a point violates the hemisphere constraint, fall back to belt logic
+                    if ((point.decDeg > hemisphere * DEC_THRESHOLD && hemisphere === -1) ||
+                        (point.decDeg < hemisphere * DEC_THRESHOLD && hemisphere === 1)) {
+                        flag = 1;
+                        poly4selection = [];
+                        break;
+                    }
+                    const p = GeomUtils.projectIn2D(point);
+                    selpoly.push(p);
+                    if (p.x > maxx)
+                        maxx = p.x;
+                    if (p.y > maxy)
+                        maxy = p.y;
+                    if (p.x < minx)
+                        minx = p.x;
+                    if (p.y < miny)
+                        miny = p.y;
+                }
+                poly4selection.push(selpoly);
+            }
+        }
+        if (flag === 0) {
+            return {
+                poly4selection,
+                flag,
+                maxx: maxx,
+                maxy: maxy,
+                minx: minx,
+                miny: miny,
+            };
+        }
+        // Case flag = 1 or 2 → work directly in (RA,Dec)
+        const RA_THRESHOLD = 180;
+        let belowThreshold = polygons[0][0].raDeg < RA_THRESHOLD;
+        maxx = minx = polygons[0][0].raDeg;
+        maxy = miny = polygons[0][0].decDeg;
+        for (const currpoly of polygons) {
+            const selpoly = [];
+            for (const point of currpoly) {
+                const p = new Point2D(point.raDeg, point.decDeg);
+                selpoly.push(p);
+                if (point.raDeg > maxx)
+                    maxx = point.raDeg;
+                if (point.decDeg > maxy)
+                    maxy = point.decDeg;
+                if (point.raDeg < minx)
+                    minx = point.raDeg;
+                if (point.decDeg < miny)
+                    miny = point.decDeg;
+                // Detect crossing of RA=0 meridian
+                if ((point.raDeg >= RA_THRESHOLD && belowThreshold) ||
+                    (point.raDeg <= RA_THRESHOLD && !belowThreshold)) {
+                    flag = 2;
+                    poly4selection = [];
+                    break;
+                }
+            }
+            poly4selection.push(selpoly);
+        }
+        if (flag === 1) {
+            return {
+                poly4selection,
+                flag,
+                maxx,
+                maxy,
+                minx,
+                miny,
+            };
+        }
+        // Case flag = 2 → shift RA>180 by -360 to unwrap around RA=0
+        let startRA = polygons[0][0].raDeg;
+        maxx = startRA >= RA_THRESHOLD ? startRA - 360 : startRA;
+        maxy = polygons[0][0].decDeg;
+        minx = maxx;
+        miny = maxy;
+        for (const currpoly of polygons) {
+            const selpoly = [];
+            for (const point of currpoly) {
+                const curra = point.raDeg >= RA_THRESHOLD ? point.raDeg - 360 : point.raDeg;
+                if (curra > maxx)
+                    maxx = curra;
+                if (point.decDeg > maxy)
+                    maxy = point.decDeg;
+                if (curra < minx)
+                    minx = curra;
+                if (point.decDeg < miny)
+                    miny = point.decDeg;
+                selpoly.push(new Point2D(curra, point.decDeg));
+            }
+            poly4selection.push(selpoly);
+        }
+        return {
+            poly4selection,
+            flag,
+            maxx,
+            maxy,
+            minx,
+            miny,
+        };
+    }
+    /** Stereographic projection from 3D point on unit sphere onto plane */
+    static stereographic(point) {
+        const x = Number(point.xyz[0]);
+        const y = Number(point.xyz[1]);
+        const z = Number(point.xyz[2]);
+        return {
+            x: (2 * x) / (1 - z),
+            y: (2 * y) / (1 - z),
+        };
+    }
+    static projectIn2D(point) {
+        const p = GeomUtils.stereographic(point);
+        return new Point2D(p.x, p.y);
+    }
+    /**
+     * Robust point-in-polygon (ray casting) using the precomputed selection object.
+     * Works with any of the three flags (0,1,2).
+     */
+    static checkPointInsidePolygon5(selectionObj, point) {
+        let p0;
+        if (selectionObj.flag === 0) {
+            p0 = GeomUtils.projectIn2D(point);
+        }
+        else if (selectionObj.flag === 1) {
+            p0 = new Point2D(point.raDeg, point.decDeg);
+        }
+        else {
+            const RA_THRESHOLD = 180;
+            const raShifted = point.raDeg >= RA_THRESHOLD ? point.raDeg - 360 : point.raDeg;
+            p0 = new Point2D(raShifted, point.decDeg);
+        }
+        const p1 = new Point2D(p0.x, p0.y + 2 * Math.abs(selectionObj.maxy - selectionObj.miny));
+        // quick reject by bbox
+        if (p0.x > selectionObj.maxx ||
+            p0.x < selectionObj.minx ||
+            p0.y > selectionObj.maxy ||
+            p0.y < selectionObj.miny) {
+            return false;
+        }
+        // Ray casting against each sub-polygon
+        for (const currpoly of selectionObj.poly4selection) {
+            let intersections = 0;
+            for (let i = 0; i < currpoly.length - 1; i++) {
+                const p2 = currpoly[i];
+                const p3 = currpoly[i + 1];
+                const denominator = (p3.y - p2.y) * (p1.x - p0.x) - (p3.x - p2.x) * (p1.y - p0.y);
+                const numerator01 = (p3.x - p2.x) * (p0.y - p2.y) - (p3.y - p2.y) * (p0.x - p2.x);
+                const numerator23 = (p1.x - p0.x) * (p0.y - p2.y) - (p1.y - p0.y) * (p0.x - p2.x);
+                if (denominator !== 0) {
+                    const lamda01 = numerator01 / denominator;
+                    const lambda23 = numerator23 / denominator;
+                    if (lamda01 >= 0 && lamda01 <= 1 && lambda23 >= 0 && lambda23 <= 1) {
+                        intersections++;
+                    }
+                }
+            }
+            // close the polygon: last with first
+            {
+                const p2 = currpoly[currpoly.length - 1];
+                const p3 = currpoly[0];
+                const denominator = (p3.y - p2.y) * (p1.x - p0.x) - (p3.x - p2.x) * (p1.y - p0.y);
+                const numerator01 = (p3.x - p2.x) * (p0.y - p2.y) - (p3.y - p2.y) * (p0.x - p2.x);
+                const numerator23 = (p1.x - p0.x) * (p0.y - p2.y) - (p1.y - p0.y) * (p0.x - p2.x);
+                if (denominator !== 0) {
+                    const lamda01 = numerator01 / denominator;
+                    const lambda23 = numerator23 / denominator;
+                    if (lamda01 >= 0 && lamda01 <= 1 && lambda23 >= 0 && lambda23 <= 1) {
+                        intersections++;
+                    }
+                }
+            }
+            if (intersections % 2 === 1) {
+                return true; // inside this subpolygon
+            }
+        }
+        return false;
+    }
+    // Legacy version kept for reference; now typed and using getters
+    static checkPointInsidePolygon4(polygons, point) {
+        const p0 = GeomUtils.projectIn2D(point);
+        let maxdist = point.raDeg + 15;
+        if (maxdist > 360)
+            maxdist = point.raDeg - 15;
+        const p1point = new Point({ raDeg: maxdist, decDeg: point.decDeg }, CoordsType.ASTRO);
+        const p1 = GeomUtils.projectIn2D(p1point);
+        for (const currpoly of polygons) {
+            let intersections = 0;
+            for (let i = 0; i < currpoly.length - 1; i++) {
+                const p2 = GeomUtils.projectIn2D(currpoly[i]);
+                const p3 = GeomUtils.projectIn2D(currpoly[i + 1]);
+                const denominator = (p3.y - p2.y) * (p1.x - p0.x) - (p3.x - p2.x) * (p1.y - p0.y);
+                const numerator01 = (p3.x - p2.x) * (p0.y - p2.y) - (p3.y - p2.y) * (p0.x - p2.x);
+                const numerator23 = (p1.x - p0.x) * (p0.y - p2.y) - (p1.y - p0.y) * (p0.x - p2.x);
+                if (denominator !== 0) {
+                    const lamda01 = numerator01 / denominator;
+                    const lambda23 = numerator23 / denominator;
+                    if (lamda01 >= 0 && lamda01 <= 1 && lambda23 >= 0 && lambda23 <= 1) {
+                        intersections++;
+                    }
+                }
+            }
+            {
+                const p2 = GeomUtils.projectIn2D(currpoly[currpoly.length - 1]);
+                const p3 = GeomUtils.projectIn2D(currpoly[0]);
+                const denominator = (p3.y - p2.y) * (p1.x - p0.x) - (p3.x - p2.x) * (p1.y - p0.y);
+                const numerator01 = (p3.x - p2.x) * (p0.y - p2.y) - (p3.y - p2.y) * (p0.x - p2.x);
+                const numerator23 = (p1.x - p0.x) * (p0.y - p2.y) - (p1.y - p0.y) * (p0.x - p2.x);
+                if (denominator !== 0) {
+                    const lamda01 = numerator01 / denominator;
+                    const lambda23 = numerator23 / denominator;
+                    if (lamda01 >= 0 && lamda01 <= 1 && lambda23 >= 0 && lambda23 <= 1) {
+                        intersections++;
+                    }
+                }
+            }
+            if (intersections % 2 === 1)
+                return true;
+        }
+        return false;
+    }
+}
+export default GeomUtils;
+//# sourceMappingURL=GeomUtils.js.map