@openglobus/og 0.18.0 → 0.18.2

package/{src/astro/astro.ts → lib/js/astro/astro.js}

@@ -1,5 +1,4 @@
-import {JulianDate, ONE_BY_SECONDS_PER_DAY, J2000} from "./jd";
-
+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
@@ -8,41 +7,36 @@ import {JulianDate, ONE_BY_SECONDS_PER_DAY, J2000} from "./jd";
  * @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: JulianDate): JulianDate {
+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/{src/astro/earth.ts → lib/js/astro/earth.js}

@@ -3,108 +3,72 @@ import { Quat } from "../math/Quat";
 import { Vec3 } from "../math/Vec3";
 import * as astro from "./astro";
 import * as jd from "./jd";
-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 function getSunPosition(jDate: JulianDate): Vec3 {
+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
-
+    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
-        )
-    );
-
+    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;