@openglobus/og 0.17.2 → 0.18.1

package/lib/js/Rectangle.d.ts

@@ -0,0 +1,76 @@
+/**
+ * 2D Rectangle class.
+ * @class
+ * @param {number} [left] - Left coordinate. 0 - default.
+ * @param {number} [top] - Top coordinate. 0 - default.
+ * @param {number} [right] - Right coordinate. 0 - default.
+ * @param {number} [bottom] - Bottom coordinate. 0 - default.
+ */
+declare class Rectangle {
+    /**
+     * Left coordinate.
+     * @public
+     * @type {number}
+     */
+    left: number;
+    /**
+     * Right coordinate.
+     * @public
+     * @type {number}
+     */
+    right: number;
+    /**
+     * Top coordinate.
+     * @public
+     * @type {number}
+     */
+    top: number;
+    /**
+     * Top coordinate.
+     * @public
+     * @type {number}
+     */
+    bottom: number;
+    constructor(left?: number, top?: number, right?: number, bottom?: number);
+    set(left?: number, top?: number, right?: number, bottom?: number): void;
+    /**
+     * Clone rectangle object.
+     * @public
+     * @returns {Rectangle}
+     */
+    clone(): Rectangle;
+    /**
+     * Returns rectangle width.
+     * @public
+     * @type {number}
+     */
+    getWidth(): number;
+    /**
+     * Returns rectangle height.
+     * @public
+     * @type {number}
+     */
+    getHeight(): number;
+    /**
+     * Returns rectangle area.
+     * @public
+     * @type {number}
+     */
+    getSquare(): number;
+    /**
+     * Returns rectangle diagonal size.
+     * @public
+     * @type {number}
+     */
+    getDiagonal(): number;
+    /**
+     * Returns true if rectangle fits their size in width and height.
+     * @public
+     * @param {number} width - Width.
+     * @param {number} height - Height.
+     * @type {boolean}
+     */
+    fit(width: number, height: number): boolean;
+    isInside(x: number, y: number): boolean;
+}
+export { Rectangle };

package/lib/js/Rectangle.js

@@ -0,0 +1,77 @@
+/**
+ * 2D Rectangle class.
+ * @class
+ * @param {number} [left] - Left coordinate. 0 - default.
+ * @param {number} [top] - Top coordinate. 0 - default.
+ * @param {number} [right] - Right coordinate. 0 - default.
+ * @param {number} [bottom] - Bottom coordinate. 0 - default.
+ */
+class Rectangle {
+    constructor(left = 0, top = 0, right = 0, bottom = 0) {
+        this.left = left;
+        this.right = right;
+        this.top = top;
+        this.bottom = bottom;
+    }
+    set(left = 0, top = 0, right = 0, bottom = 0) {
+        this.left = left;
+        this.right = right;
+        this.top = top;
+        this.bottom = bottom;
+    }
+    /**
+     * Clone rectangle object.
+     * @public
+     * @returns {Rectangle}
+     */
+    clone() {
+        return new Rectangle(this.left, this.top, this.right, this.bottom);
+    }
+    /**
+     * Returns rectangle width.
+     * @public
+     * @type {number}
+     */
+    getWidth() {
+        return Math.abs(this.right - this.left);
+    }
+    /**
+     * Returns rectangle height.
+     * @public
+     * @type {number}
+     */
+    getHeight() {
+        return Math.abs(this.bottom - this.top);
+    }
+    /**
+     * Returns rectangle area.
+     * @public
+     * @type {number}
+     */
+    getSquare() {
+        return this.getHeight() * this.getWidth();
+    }
+    /**
+     * Returns rectangle diagonal size.
+     * @public
+     * @type {number}
+     */
+    getDiagonal() {
+        let w = this.getWidth(), h = this.getHeight();
+        return Math.sqrt(h * h + w * w);
+    }
+    /**
+     * Returns true if rectangle fits their size in width and height.
+     * @public
+     * @param {number} width - Width.
+     * @param {number} height - Height.
+     * @type {boolean}
+     */
+    fit(width, height) {
+        return this.getWidth() === width && this.getHeight() === height;
+    }
+    isInside(x, y) {
+        return x >= this.left && x <= this.right && y >= this.top && y <= this.bottom;
+    }
+}
+export { Rectangle };

package/lib/js/Stack.d.ts

@@ -0,0 +1,20 @@
+declare class Node<T> {
+    next: Node<T> | null;
+    prev: Node<T> | null;
+    data: T | null;
+    constructor();
+}
+/**
+ * @class Stack<T>
+ * @param {number} [size=256] - Stack size
+ */
+declare class Stack<T> {
+    protected _current: Node<T>;
+    protected _head: Node<T>;
+    constructor(size?: number);
+    current(): Node<T>;
+    push(data: T): void;
+    pop(): T | null;
+    popPrev(): T | null;
+}
+export { Stack };

package/lib/js/Stack.js

@@ -0,0 +1,42 @@
+class Node {
+    constructor() {
+        this.next = null;
+        this.prev = null;
+        this.data = null;
+    }
+}
+/**
+ * @class Stack<T>
+ * @param {number} [size=256] - Stack size
+ */
+class Stack {
+    constructor(size = 256) {
+        this._current = new Node();
+        this._head = this._current;
+        for (let i = 0; i < size; i++) {
+            let n = new Node();
+            n.prev = this._current;
+            this._current.next = n;
+            this._current = n;
+        }
+        this._current = this._head;
+    }
+    current() {
+        return this._current;
+    }
+    push(data) {
+        this._current = this._current.next;
+        this._current.data = data;
+    }
+    pop() {
+        let res = this._current.data;
+        //this._current.data = null;
+        this._current = this._current.prev;
+        return res;
+    }
+    popPrev() {
+        this._current = this._current.prev;
+        return this._current.data;
+    }
+}
+export { Stack };

package/lib/js/astro/astro.d.ts

@@ -0,0 +1,39 @@
+import { JulianDate } from "./jd";
+/**
+ * Angle between J2000 mean equator and the ecliptic plane.
+ * 23 deg 26' 21".448 (Seidelmann, _Explanatory Supplement to the
+ * Astronomical Almanac_ (1992), eqn 3.222-1.
+ * @const
+ * @type{number}
+ */
+export declare const J2000_OBLIQUITY = 23.4392911;
+/**
+ * IAU 1976 value
+ * @const
+ * @type{number}
+ */
+export declare const AU_TO_METERS = 149597870000;
+/**
+ * Terrestrial and atomic time difference.
+ * @const
+ * @type{number}
+ */
+export declare const TDT_TAI = 32.184;
+/**
+ * Earth gravitational parameter product of gravitational constant G and the mass M of the Earth.
+ * @const
+ * @type{number}
+ */
+export declare const EARTH_GRAVITATIONAL_PARAMETER = 398600435000000;
+/**
+ * Sun gravitational parameter product of gravitational constant G and the mass M of the Sun.
+ * @const
+ * @type{number}
+ */
+export declare const SUN_GRAVITATIONAL_PARAMETER = 132712440018000000000;
+/**
+ * Converts atomic time to barycentric dynamical time.
+ * @param {JulianDate} tai - Atomic time.
+ * @returns {JulianDate} - returns barycentric dynamical time.
+ */
+export declare function TAItoTDB(tai: JulianDate): JulianDate;

package/lib/js/astro/astro.js

@@ -0,0 +1,43 @@
+import { ONE_BY_SECONDS_PER_DAY, J2000 } from "./jd";
+/**
+ * Angle between J2000 mean equator and the ecliptic plane.
+ * 23 deg 26' 21".448 (Seidelmann, _Explanatory Supplement to the
+ * Astronomical Almanac_ (1992), eqn 3.222-1.
+ * @const
+ * @type{number}
+ */
+export const J2000_OBLIQUITY = 23.4392911;
+/**
+ * IAU 1976 value
+ * @const
+ * @type{number}
+ */
+export const AU_TO_METERS = 1.4959787e11;
+/**
+ * Terrestrial and atomic time difference.
+ * @const
+ * @type{number}
+ */
+export const TDT_TAI = 32.184;
+/**
+ * Earth gravitational parameter product of gravitational constant G and the mass M of the Earth.
+ * @const
+ * @type{number}
+ */
+export const EARTH_GRAVITATIONAL_PARAMETER = 3.98600435e14;
+/**
+ * Sun gravitational parameter product of gravitational constant G and the mass M of the Sun.
+ * @const
+ * @type{number}
+ */
+export const SUN_GRAVITATIONAL_PARAMETER = 1.32712440018e20;
+/**
+ * Converts atomic time to barycentric dynamical time.
+ * @param {JulianDate} tai - Atomic time.
+ * @returns {JulianDate} - returns barycentric dynamical time.
+ */
+export function TAItoTDB(tai) {
+    tai += TDT_TAI * ONE_BY_SECONDS_PER_DAY;
+    const g = 6.239996 + 0.0172019696544 * (tai - J2000);
+    return tai + 0.001658 * Math.sin(g + 1.671e-2 * Math.sin(g)) * ONE_BY_SECONDS_PER_DAY;
+}

package/lib/js/astro/earth.d.ts

@@ -0,0 +1,8 @@
+import { Vec3 } from "../math/Vec3";
+import { JulianDate } from "./jd";
+/**
+ * Returns Sun position in the geocentric coordinate system by the time.
+ * @param {JulianDate} jDate - Julian date time.
+ * @returns {Vec3} - Sun geocentric coordinates.
+ */
+export declare function getSunPosition(jDate: JulianDate): Vec3;

package/lib/js/astro/earth.js

@@ -0,0 +1,75 @@
+import * as math from "../math";
+import { Quat } from "../math/Quat";
+import { Vec3 } from "../math/Vec3";
+import * as astro from "./astro";
+import * as jd from "./jd";
+/**
+ * Returns Sun position in the geocentric coordinate system by the time.
+ * @param {JulianDate} jDate - Julian date time.
+ * @returns {Vec3} - Sun geocentric coordinates.
+ */
+export function getSunPosition(jDate) {
+    // http://stjarnhimlen.se/comp/tutorial.html
+    // a  Mean distance, or semi-major axis
+    // e  Eccentricity
+    // T  Time at perihelion
+    // q  Perihelion distance  = a * (1 - e)
+    // Q  Aphelion distance    = a * (1 + e)
+    // i  Inclination, i.e. the "tilt" of the orbit relative to the
+    //    ecliptic.  The inclination varies from 0 to 180 degrees. If
+    //    the inclination is larger than 90 degrees, the planet is in
+    //    a retrogade orbit, i.e. it moves "backwards".  The most
+    //    well-known celestial body with retrogade motion is Comet Halley.
+    // N  (usually written as "Capital Omega") Longitude of Ascending
+    //    Node. This is the angle, along the ecliptic, from the Vernal
+    //    Point to the Ascending Node, which is the intersection between
+    //    the orbit and the ecliptic, where the planet moves from south
+    //    of to north of the ecliptic, i.e. from negative to positive
+    //    latitudes.
+    // w  (usually written as "small Omega") The angle from the Ascending
+    //    node to the Perihelion, along the orbit.
+    // P  Orbital period       = 365.256898326 * a**1.5/sqrt(1+m) days,
+    //    where m = the mass of the planet in solar masses (0 for
+    //    comets and asteroids). sqrt() is the square root function.
+    // n  Daily motion         = 360_deg / P    degrees/day
+    // t  Some epoch as a day count, e.g. Julian Day Number. The Time
+    //    at Perihelion, T, should then be expressed as the same day count.
+    // t - T   Time since Perihelion, usually in days
+    // M  Mean Anomaly         = n * (t - T)  =  (t - T) * 360_deg / P
+    //    Mean Anomaly is 0 at perihelion and 180 degrees at aphelion
+    // L  Mean Longitude       = M + w + N
+    // E  Eccentric anomaly, defined by Kepler's equation:   M = E - e * sin(E)
+    //    An auxiliary angle to compute the position in an elliptic orbit
+    // v  True anomaly: the angle from perihelion to the planet, as seen
+    //    from the Sun
+    // r  Heliocentric distance: the planet's distance from the Sun.
+    // x,y,z  Rectangular coordinates. Used e.g. when a heliocentric
+    //        position (seen from the Sun) should be converted to a
+    //        corresponding geocentric position (seen from the Earth).
+    var d = jDate - jd.J2000;
+    var w = 282.9404 + 4.70935e-5 * d; // longitude of perihelion
+    // var a = 1.000000; // mean distance, a.u.
+    var e = 0.016709 - 1.151e-9 * d; // eccentricity
+    var M = math.rev(356.047 + 0.9856002585 * d); // mean anomaly
+    var oblecl = astro.J2000_OBLIQUITY - 3.563e-7 * d; // obliquity of the ecliptic
+    // var L = math.rev(w + M); // Sun's mean longitude
+    var E = M + math.DEGREES * e * Math.sin(M * math.RADIANS) * (1 + e * Math.cos(M * math.RADIANS)); // eccentric anomaly
+    // Sun rectangular coordinates, where the X axis points towards the perihelion
+    var x = Math.cos(E * math.RADIANS) - e;
+    var y = Math.sin(E * math.RADIANS) * Math.sqrt(1 - e * e);
+    var r = Math.sqrt(x * x + y * y); // distance
+    var v = Math.atan2(y, x) * math.DEGREES; // true anomaly
+    var lon = math.rev(v + w); // longitude of the Sun
+    // the Sun's ecliptic rectangular coordinates
+    x = r * Math.cos(lon * math.RADIANS);
+    y = r * Math.sin(lon * math.RADIANS);
+    // We use oblecl, and rotate these coordinates
+    var xequat = x;
+    var yequat = y * Math.cos(oblecl * math.RADIANS);
+    var zequat = y * Math.sin(oblecl * math.RADIANS);
+    var theta = math.TWO_PI * ((d * 24.0) / 23.9344694 - 259.853 / 360.0); // Siderial spin time
+    return Quat.zRotation(-theta).mulVec3(new Vec3(-xequat * astro.AU_TO_METERS, -yequat * astro.AU_TO_METERS, zequat * astro.AU_TO_METERS));
+    // Convert to RA and Decl
+    // var RA = Math.atan2(yequat, xequat) * math.DEGREES;
+    // var Decl = Math.atan2(zequat, Math.sqrt(xequat * xequat + yequat * yequat)) * math.DEGREES;
+}

package/lib/js/astro/jd.d.ts

@@ -0,0 +1,234 @@
+export type JulianDate = number;
+/**
+ * Seconds in millisecond.
+ * @const
+ */
+export declare const SECONDS_PER_MILLISECOND = 0.001;
+/**
+ * Milliseconds in second.
+ * @const
+ */
+export declare const MILLISECONDS_PER_SECOND = 1000;
+/**
+ * Seconds in minute.
+ * @const
+ */
+export declare const SECONDS_PER_MINUTE = 60;
+/**
+ * One by seconds in minute.
+ * @const
+ */
+export declare const ONE_BY_SECONDS_PER_MINUTE: number;
+/**
+ * Minutes in hour.
+ * @const
+ */
+export declare const MINUTES_PER_HOUR = 60;
+/**
+ * Hours in day.
+ * @const
+ */
+export declare const HOURS_PER_DAY = 24;
+/**
+ * One by hours in day.
+ * @const
+ */
+export declare const ONE_BY_HOURS_PER_DAY: number;
+/**
+ * Seconds in hour.
+ * @const
+ */
+export declare const SECONDS_PER_HOUR = 3600;
+/**
+ * One by seconds in hour.
+ * @const
+ */
+export declare const ONE_BY_SECONDS_PER_HOUR: number;
+/**
+ * Seconds in 12 hours.
+ * @const
+ */
+export declare const SECONDS_PER_12_HOURS: number;
+/**
+ * Minutes in day.
+ * @const
+ */
+export declare const MINUTES_PER_DAY = 1440;
+/**
+ * One by minutes in day.
+ * @const
+ */
+export declare const ONE_BY_MINUTES_PER_DAY: number;
+/**
+ * Seconds in day.
+ * @const
+ */
+export declare const SECONDS_PER_DAY = 86400;
+/**
+ * Milliseconds in day.
+ * @const
+ */
+export declare const MILLISECONDS_PER_DAY = 86400000;
+/**
+ * One by milliseconds in day.
+ * @const
+ */
+export declare const ONE_BY_MILLISECONDS_PER_DAY: number;
+/**
+ * One by seconds in day.
+ * @const
+ */
+export declare const ONE_BY_SECONDS_PER_DAY: number;
+/**
+ * Days in julian century.
+ * @const
+ */
+export declare const DAYS_PER_JULIAN_CENTURY = 36525;
+/**
+ * Days in julian year.
+ * @const
+ */
+export declare const DAYS_PER_JULIAN_YEAR = 365.25;
+/**
+ * Seconds in picoseconds.
+ * @const
+ */
+export declare const PICOSECOND = 1e-9;
+/**
+ * Modified julian date difference.
+ * @const
+ */
+export declare const MODIFIED_JULIAN_DATE_DIFFERENCE = 2400000.5;
+/**
+ * Julian date of 2000 year. Epoch.
+ * @const
+ */
+export declare const J2000 = 2451545;
+/**
+ * Returns julian days from Epoch.
+ * @param {JulianDate} jd - Julian date.
+ * @returns {number} Days from epoch
+ */
+export declare function T(jd: JulianDate): number;
+/**
+ * Gets the date's julian day.
+ * @param {number} year - Year.
+ * @param {number} month - Month.
+ * @param {number} day - Day.
+ * @returns {number} Day number
+ */
+export declare function getDayNumber(year: number, month: number, day: number): number;
+/**
+ * Converts javascript date to the universal(UTC) julian date.
+ * @param {Date} date - Date.
+ * @returns {JulianDate} UTC julian date
+ */
+export declare function DateToUTC(date: Date): JulianDate;
+/**
+ * Converts javascript date to the atomic(TAI) julian date.
+ * @param {Date} date - Date.
+ * @returns {JulianDate} TAI julian date
+ */
+export declare function DateToTAI(date: Date): JulianDate;
+/**
+ * Converts coordinated universal(UTC) julian date to atomic(TAI) julian date.
+ * @param {JulianDate} jd - UTC julian date.
+ * @returns {JulianDate} TAI julian date
+ */
+export declare function UTCtoTAI(jd: JulianDate): JulianDate;
+/**
+ * Converts atomic julian date(TAI) to the coordinated universal(UTC) julian date.
+ * @param {JulianDate} tai - TAI julian date.
+ * @returns {JulianDate | undefined} UTC julian date
+ */
+export declare function TAItoUTC(tai: JulianDate): JulianDate | undefined;
+/**
+ * Converts UTC julian date to the javascript date object.
+ * @param {JulianDate} utc - UTC julian date.
+ * @returns {Date} JavaScript Date object
+ */
+export declare function UTCtoDate(utc: JulianDate): Date;
+/**
+ * Converts TAI julian date to the javascript date object.
+ * @param {JulianDate} tai - TAI julian date.
+ * @returns {Date} JavaScript Date object
+ */
+export declare function TAItoDate(tai: JulianDate): Date;
+/**
+ * Adds milliseconds to the julian date.
+ * @param {JulianDate} jd - Julian date.
+ * @param {number} milliseconds - Milliseconds to add.
+ * @returns {JulianDate} Julian date
+ */
+export declare function addMilliseconds(jd: JulianDate, milliseconds: number): JulianDate;
+/**
+ * Adds seconds to the julian date.
+ * @param {JulianDate} jd - Julian date.
+ * @param {number} seconds - Seconds to add.
+ * @returns {JulianDate} Julian date
+ */
+export declare function addSeconds(jd: JulianDate, seconds: number): JulianDate;
+/**
+ * Adds hours to the julian date.
+ * @param {JulianDate} jd - Julian date.
+ * @param {number} hours - Hours to add.
+ * @returns {JulianDate} Julian date
+ */
+export declare function addHours(jd: JulianDate, hours: number): JulianDate;
+/**
+ * Adds minutes to the julian date.
+ * @param {JulianDate} jd - Julian date.
+ * @param {number} minutes - Minutes to add.
+ * @returns {JulianDate} Julian date
+ */
+export declare function addMinutes(jd: JulianDate, minutes: number): JulianDate;
+/**
+ * Adds days to the julian date.
+ * @param {JulianDate} jd - Julian date.
+ * @param {number} days - Days to add.
+ * @returns {JulianDate} Julian date
+ */
+export declare function addDays(jd: JulianDate, days: number): JulianDate;
+/**
+ * Gets milliseconds of a julian date.
+ * @param {JulianDate} jd - julian date.
+ * @returns {number} Milliseconds
+ */
+export declare function getMilliseconds(jd: JulianDate): number;
+/**
+ * Gets seconds of a julian date.
+ * @param {JulianDate} jd - julian date.
+ * @returns {number} Seconds
+ */
+export declare function getSeconds(jd: JulianDate): number;
+/**
+ * Gets hours of a julian date.
+ * @param {JulianDate} jd - julian date.
+ * @returns {number} Hours
+ */
+export declare function getHours(jd: JulianDate): number;
+/**
+ * Gets minutes of a julian date.
+ * @param {JulianDate} jd - julian date.
+ * @returns {number} Minutes
+ */
+export declare function getMinutes(jd: JulianDate): number;
+/**
+ * Gets days of a julian date.
+ * @param {JulianDate} jd - julian date.
+ * @returns {number} Days
+ */
+export declare function getDays(jd: JulianDate): number;
+/**
+ * Returns days in seconds.
+ * @param {number} s - Seconds.
+ * @returns {number} Days
+ */
+export declare function secondsToDays(s: number): number;
+/**
+ * Returns seconds in days.
+ * @param {number} d - Days.
+ * @returns {number} Seconds
+ */
+export declare function daysToSeconds(d: number): number;
+export declare const J2000TAI: number;