Skip to content

Latest commit

 

History

History
498 lines (383 loc) · 16.2 KB

usage.en.md

File metadata and controls

498 lines (383 loc) · 16.2 KB

orb.js(v2)- JavaScript library for astronomical calculations

Please note: There is no compatibility v1 and v2.

Files

orb.v2.js:
all-in-one package

orb-core.v2.js:
Truncated version of orb.js. If you only need calculations for time conversion or coordinate conversion, you can use this .js file

orb-satellite.v2.js:
Truncated version of orb.js. If you only need calculations for artificial satellites or space debris, you can use this .js file

orb-planetary.v2.js:
Truncated version of orb.js. If you only need calculations for planetary objects such as planets or asteroids, you can use this .js file

CDN

<script src="https://cdn.jsdelivr.net/gh/lizard-isana/[email protected]/build/orb.v2.js"></script>

or minified file

<script src="https://cdn.jsdelivr.net/gh/lizard-isana/[email protected]/build/min/orb.v2.min.js"></script>

Example

var date = new Date();

// Position of planets
var mars = new Orb.VSOP("Mars");
var xyz = mars.xyz(date); // ecliptic rectangular coordinates(x, y, z)
var radec = mars.radec(date); // equatorial spherical coordinates(ra, dec, distance)

// Position of the moon
var luna = new Orb.Luna();
var xyz = luna.xyz(date); // Earth centered ecliptic rectangular coordinates (x, y, z)
var radec = luna.radec(date); // equatorial spherical coordinates(ra, dec, distance)

// Apparent position of the Sun
var sun = new Orb.Sun();
var xyz = sun.xyz(date); // Earth centered equatorial rectangular coordinates (x, y, z)
var radec = sun.radec(date); // equatorial spherical coordinates(ra, dec, distance)

// Kepler orbital elements
var asteroid = new Orb.Kepler({
  "gm": 2.9591220828559093*Math.pow(10,-4), //(au^3/d^2) default value is the sun, so you can omit this line.
  "eccentricity":"0.08728849329001058",
  "inclination":"6.812676631845272",
  "longitude_of_ascending_node":"250.5660658100269",
  "argument_of_periapsis":"95.63473165761138",
  "time_of_periapsis":"2456918.756066796",
  "semi_major_axis":"1.001911878091084"
});
var xyz = asteroid.xyz(date); // ecliptic rectangular coordinates(x, y, z, xdot, ydot, zdot)
var radec = asteroid.radec(date); // equatorial spherical coordinates(ra, dec, distance)

// Cartesian state vectors to Kepler orbital elements
var orbital_elements = new Orb.CartesianToKeplerian(xyz)

// Position of artificial satellites from Two Line Elements(TLE)
var tle = {
  first_line:"1 25544U 98067A   15012.59173611  .00014829  00000-0  23845-3 0  7179",
  second_line:"2 25544 051.6466 140.7335 0006107 243.2909 291.5211 15.53213268923827"
}
var satellite = new Orb.SGP4(tle);
var xyz = satellite.xyz(date); // Earth centered equatorial rectangular coordinates (x, y, z, xdot, ydot, zdot)
var latlng = satellite.latlng(date); // geographic spherical coordinates(latitude, longitude, altitude, velocity)

// Azimuth and Elevation
var your_location = {
  "latitude":35.658,
  "longitude":139.741,
  "altitude":0
}

var observe_mars = new Orb.Observation({"observer":your_location,"target":mars});
var horizontal = observe_mars.azel(date); // horizontal coordinates(azimuth, elevation)

var observe_satellite = new Orb.Observation({"observer":your_location,"target":satellite});
var horizontal = observe_satellite.azel(date); // horizontal coordinates(azimuth, elevation)

var sirius = {
  "ra":6.75257,
  "dec":-16.7131,
  "distance":543300
}

var observe_star = new Orb.Observation({
  "observer":your_location ,
  "target":sirius
});
var horizontal = observe_star.azel(date); // horizontal coordinates(azimuth, elevation)

// Time conversion
var time = new Orb.Time(date);
var gmst = time.gmst();
var julian_date = time.jd();
var time_in_day = time.time_in_day();
var day_of_year = time.doy();

// Coordinates conversion
var equatorial_rectangular = Orb.RadecToXYZ(sirius)
var ecliptic_rectangular = Orb.EquatorialToEcliptic({"date":date,"equatorial":equatorial_rectangular})
var equatorial_rectangular = Orb.EclipticToEquatorial({"date":date,"ecliptic":ecliptic_rectangular})
var equatorial_spherical = Orb.XYZtoRadec(equatorial_rectangular)

orb-core.v2.js

Orb

Returns information of version, author, license

{
  "VERSION":<string>,
  "AUTHOR":<string>,
  "LICENSE":<string>
}

Orb.Const

Various numerical constants for astronomical calculations

}
  "PI","RAD","AU","RE","LD","LY","PC","PARSEC","G": <number>,
  "Planets":<array>, //list of planets,
  "Sun","Mercury","Venus","Earth","Moon","Mars","Jupiter","Saturn","Uranus","Neptune":{
  "radius":<number>, //radius at equator,
  "obliquity":<number>, //obliquity of the planet,
  "mass":<number>, //mass of the planet,
  "gm":<number> //gravitational constant for the planet
}

Orb.Time()

Convert Date to the astronomical time units

// initialize
var date = new Date();
var time = new Orb.Time(date);
// if you omit date object, you get value of now.

var gmst = time.gmst();
var julian_date = time.jd();
var time_in_day = time.time_in_day();
var day_of_year = time.doy();

Orb.Observation()

Calc horizontal coordinates of target from observer.

// initialize
var date = new Date()
var your_location = {
  "latitude":35.658,
  "longitude":139.741,
  "altitude":0
}

// for planets (see Orb.VSOP)
var mars = new Orb.VSOP("Mars");
var obs = new Orb.Observation({"observer":your_location,"target":mars});
var azel = obs.azel(date); // azimuth, elevation -> horizontal coordinates

// for artificial satellites (see Orb.SGP4)
var tle = {
  first_line:"1 25544U 98067A   15012.59173611  .00014829  00000-0  23845-3 0  7179",
  second_line:"2 25544 051.6466 140.7335 0006107 243.2909 291.5211 15.53213268923827"
}
var satellite = new Orb.SGP4(tle);
var obs = new Orb.Observation({"observer":your_location,"target":satellite});
var azel = obs.azel(date); // azimuth, elevation -> horizontal coordinates

// for equatorial sphirical coordinates (RA,Dec)
var obs = new Orb.Observation({
  "observer":your_location ,
  "target":{"ra":0,"dec":0}
});
var azel = obs.azel(date); // azimuth, elevation

// for equatorial rectangular coordinates (x,y,z)
var obs = new Orb.Observation({
  "observer":your_location ,
  "target":{"x":0,"y":0,"z":0}
});

var azel = obs.azel(date); // azimuth, elevation

// if target includes "coorinate_keywords" with "ecliptic", orb.js automatically convert coordinates.
// "coorinate_keywords":"ecliptic rectangular"

    var obs = new Orb.Observation({
      "observer":your_location ,
      "target":{
        "x":0,
        "y":0,
        "z":0,
        "date":new Date(),
        "coorinate_keywords":"ecliptic rectangular"
      }
    });
    //if you omit date object, you get value of the time of observation.

    var azel = obs.azel(date); // azimuth, elevation

Orb.RadecToXYZ()

Convert equatorial spherical coordinates(RA, Dec) to equatorial rectangular coordinates(x,y,z) Please note: "ra"(right ascension) must be hours.

var sirius = {
  "ra":6.75257,
  "dec":-16.7131,
  "distance":543300
}
var xyz = Orb.RadecToXYZ(sirius)

Orb.XYZtoRadec()

Convert equatorial rectangular coordinates(x,y,z) to equatorial spherical coordinates(RA, Dec) Please note: "ra"(right ascension) must be hours.

var sirius = {
  "x": -101858.45670898016,
  "y": 510282.46985869075,
  "z": -156241.94619812947
}

var radec = Orb.XYZtoRadec(sirius)

if "coorinate_keywords" includes "ecliptic" (and optional "date"), this method automatically convert ecliptic coordinate to equatorial coordinate.

var sirius = {
  "x": -101858.45670898016,
  "y": 510282.46985869075,
  "z": -156241.94619812947,
  "date":new Date(),
  "coorinate_keywords":"ecliptic rectangular"
}
//if you omit date object, you get value of now.

var radec = Orb.XYZtoRadec(sirius)

Orb.EquatorialToEcliptic()

Convert equatorial rectangular coordinates(x,y,z) to ecliptic rectangular coordinates(x,y,z)

var sirius = {
  "x": -101858.45670898016,
  "y": 510282.46985869075,
  "z": -156241.94619812947
}
var ecliptic = Orb.EquatorialToEcliptic({
  "equatorial": sirius,
  "date":date
})

Orb.EclipticToEquatorial()

Convert ecliptic rectangular coordinates(x,y,z) to equatorial rectangular coordinates(x,y,z)

 var sirius = {
   "x":-101858.45670898016,
   "y":510282.46985869075,
   "z":-156241.94619812947
 }
var equatorial = Orb.EclipticToEquatorial({
  "ecliptic": sirius,
  "date":date
})

Orb.NutationAndObliquity()

Calc nutation and obliquity of Earth

var date = new Date()
var nao = Orb.NutationAndObliquity(date)
var obliquity = nao.obliquity()
var nutation = nao.nutation()

Orb.Obliquity()

Calc obliquity of Earth. This is wrapper function of Orb.NutationAndObliquity()

var date = new Date()
var obliquity = Orb.Obliquity(date)

Orb.RoundAngle()

var degree = 400;
var round = Orb.RoundAngle(degree) // returns 40

var degree = -10;
var round = Orb.RoundAngle(degree) // returns 350

Orb.ZeroFill()

var num = 1;
var fill = Orb.ZeroFill(num,2) // returns 01
var fill = Orb.ZeroFill(num,4) // returns 0001

orb-satellite.v2.js

Orb.SGP4()

Calc position of an Earth-orbiting object from Two Line Elements(TLE).

//initialize
var tle = {
  first_line:"1 25544U 98067A   15012.59173611  .00014829  00000-0  23845-3 0  7179",
  second_line:"2 25544 051.6466 140.7335 0006107 243.2909 291.5211 15.53213268923827"
}
var satellite = new Orb.SGP4(tle);
var xyz = satellite.xyz(date); // equatorial rectangular coordinates (Earth centered): x, y, z, xdot, ydot, zdot
var latlng = satellite.latlng(date); // geographic spherical coordinates :  latitude, longitude, altitude, velocity

ORBIT MEAN-ELEMENTS MESSAGE(OMM) format is also acceptable
cf. ORBIT DATA MESSAGES(CCSDS 502.0-B-2) https://public.ccsds.org/pubs/502x0b2c1.pdf

const omm = {
  "CCSDS_OMM_VERS": "2.0",
  "COMMENT": "GENERATED VIA SPACE-TRACK.ORG API",
  "CREATION_DATE": "2020-06-24 23:30:28",
  "ORIGINATOR": "18 SPCS",
  "OBJECT_NAME": "ISS (ZARYA)",
  "OBJECT_ID": "1998-067A",
  "CENTER_NAME": "EARTH",
  "REF_FRAME": "TEME",
  "TIME_SYSTEM": "UTC",
  "MEAN_ELEMENT_THEORY": "SGP4",
  "EPOCH": "2020-06-24T23:30:28",
  "MEAN_MOTION": "15.49457702",
  "ECCENTRICITY": "0.0002538",
  "INCLINATION": "51.6446",
  "RA_OF_ASC_NODE": "309.8972",
  "ARG_OF_PERICENTER": "80.4322",
  "MEAN_ANOMALY": "66.556",
  "EPHEMERIS_TYPE": "0",
  "CLASSIFICATION_TYPE": "U",
  "NORAD_CAT_ID": "25544",
  "ELEMENT_SET_NO": "999",
  "REV_AT_EPOCH": "23322",
  "BSTAR": "0.000017286",
  "MEAN_MOTION_DOT": "0.00000516",
  "MEAN_MOTION_DDOT": "0",
  "USER_DEFINED_TLE_LINE0": "0 ISS (ZARYA)",
  "USER_DEFINED_TLE_LINE1": "1 25544U 98067A   20176.97949255  .00000516  00000-0  17286-4 0  9996",
  "USER_DEFINED_TLE_LINE2": "2 25544  51.6446 309.8972 0002538  80.4322  66.5560 15.49457702233227"
}
const satellite = new Orb.SGP4(omm);

orb-planetary.v2.js (require orb-core.v2.js)

Orb.Kepler()

Orb.KeplerianToCartesian() is same as above

//initialize
var asteroid = new Orb.Kepler({
  "gm": 2.9591220828559093E-4; //in au^3/d^2
  "eccentricity":"0.08728849329001058",
  "inclination":"6.812676631845272", //in degree
  "longitude_of_ascending_node":"250.5660658100269", //in degree
  "argument_of_periapsis":"95.63473165761138", //in degree
  "time_of_periapsis":"2456918.756066796", //in jd
  "semi_major_axis":"1.001911878091084" //in au
});

var xyz = asteroid.xyz(date); // x, y, z, xdot, ydot, zdot -> ecliptic rectangular coordinates in au and au/d
var radec = asteroid.radec(date); // ra, dec, distance -> equatorial spherical coordinates, distance in au

Please note, Parabolic(eccentricity<0) and Hyperbolic(eccentricity=0) orbits need the "periapsis_distance" insted of "semi_major_axis".

"gm" has default value "2.9591220828559093E-4" so you can omit "gm" for the Sun.

Orb.CartesianToKeplerian()

Orb.Cartesian() is same as above var date = new Date(); date.setTime(Date.UTC(2017,0,1,0,0,0))

//initialize
var orbital_elements = new Orb.CartesianToKeplerian({
  "gm": 2.9591220828559093*Math.pow(10,-4), //au^3/d^2
  "date":date, // or epoch: 2457754.5,
  "x": 0.0830237594569403,
  "y": -3.1268511124864538,
  "z": 4.499475953917434,
  "xdot": -0.002473803722068218,
  "ydot": 0.009696903602742064,
  "zdot": -0.015396150337498575
})

// returns keplerian orbital elements

orbital_elements = {
  epoch: 2457754.5,
  semi_major_axis: -1.2911100416899044,
  eccentricity: 1.197095803399395,
  inclination: 122.60100166153163,
  longitude_of_ascending_node: 24.598195301118153,
  true_anomaly: 221.3987070271094,
  mean_anomaly: 191.04774671275925,
  mean_motion: 0.6718296791172355,
  time_of_periapsis: 2458005.9807823156,
  "argument_of_periapsis": 241.53003549784427,
  "periapsis_distance": 0.2544723709438983
}

"gm" has default value "2.9591220828559093E-4" so you can omit "gm" for the Sun.

Orb.VSOP()

["Mercury","Venus","Earth","Moon","Mars","Jupiter","Saturn","Uranus","Neptune"]

//initialize
var mars = new Orb.VSOP("Mars");
// var mars = new Orb.Mars();
// this is same as above

var xyz = mars.xyz(date); // x, y, z -> ecliptic rectangular coordinates
var radec = mars.radec(date); // ra, dec, distance -> equatorial spherical coordinates

Orb.Luna()

Position and Phase of the moon seen from the Earth. Please note: xyz() returns Earth centered "ecliptic" rectangular coordinates not "equatorial"

//initialize
var luna = new Orb.Luna();
// var luna = new Orb.Moon();
// this is same as above

var xyz = luna.xyz(date); // x, y, z -> ecliptic rectangular coordinates (Earth centered)
var radec = luna.radec(date); // ra, dec, distance -> equatorial spherical coordinates

Orb.Sun()

position of the sun seen from the Earth. Please note: xyz() returns Earth centered "equatorial" rectangular coordinates not "ecliptic"

//initialize
var sun = new Orb.Sun();

var xyz = sun.xyz(date); // x, y, z -> equatorial rectangular coordinates (Earth centered)
var radec = sun.radec(date); // ra, dec, distance -> equatorial spherical coordinates

sapplemental

orb-data-handler.js

Orb.DataLoader()

var data = Orb.DataLoader({
  format:"text", // or "json" or "xml"
  path: <path_to_data>,
  id:"id",
  ajax:true, // or "false"
  callback:function(data,id){
  // do something
  }
})

orb-date-handler.js

Orb.DigitsToDate()

var date = Orb.DigitsToDate("20170101120000"); // return date object 2017.01.01 12:00:00 UTC

Orb.DateToDigits()

var date = new Date();
date.setTime(Date.UTC(2017,0,1,12,0,0))    
var digits = Orb.DateToDigits(date); //return "20170101120000"

Orb.StringToDate()

var str = "2017-01-01T12:00:00"    
//date separator should be "." or "-", time and date separator should be " " or "T" and time separator should be ":"

var date = Orb.StringToDate(str)

Orb.FormatUTCDate()

var date = new Date();
date.setTime(Date.UTC(2017,0,1,12,0,0))  
var str = Orb.FormatUTCDate(date) // return "2017.01.01 12:00:00"

Orb.FormatLocalDate()

var date = new Date();
date.setTime(Date.UTC(2017,0,1,12,0,0))  
var str = Orb.FormatUTCDate(date) // return "2017.01.01 21:00:00" in timezone GMT+9:00

Reference

  • Bretagnon, P.; Francou, G. "Planetary theories in rectangular and spherical variables - VSOP 87 solutions". Astronomy & Astrophysics,1988
  • Jean Meeus. Astronomical Algorithms second edition. Willmann-Bell, 1999
  • 長澤 工. 天体の位置計算 増補版. 地人書館, 1985

License

Copyright (c) 2012-2017 Isana Kashiwai
Licensed under the MIT license.

Administrator

Isana Kashiwai email: isana.k at gmail.com