@cesium/engine 9.2.0 → 10.1.0

package/Source/Core/Rectangle.js

@@ -285,7 +285,7 @@ Rectangle.fromCartographicArray = function (cartographics, result) {
  * Creates the smallest possible Rectangle that encloses all positions in the provided array.
  *
  * @param {Cartesian3[]} cartesians The list of Cartesian instances.
- * @param {Ellipsoid} [ellipsoid=Ellipsoid.WGS84] The ellipsoid the cartesians are on.
+ * @param {Ellipsoid} [ellipsoid=Ellipsoid.default] The ellipsoid the cartesians are on.
  * @param {Rectangle} [result] The object onto which to store the result, or undefined if a new instance should be created.
  * @returns {Rectangle} The modified result parameter or a new Rectangle instance if none was provided.
  */
@@ -293,7 +293,7 @@ Rectangle.fromCartesianArray = function (cartesians, ellipsoid, result) {
   //>>includeStart('debug', pragmas.debug);
   Check.defined("cartesians", cartesians);
   //>>includeEnd('debug');
-  ellipsoid = defaultValue(ellipsoid, Ellipsoid.WGS84);
+  ellipsoid = defaultValue(ellipsoid, Ellipsoid.default);
 
   let west = Number.MAX_VALUE;
   let east = -Number.MAX_VALUE;
@@ -354,7 +354,7 @@ for (let n = 0; n < fromBoundingSpherePositionsScratch.length; ++n) {
  *
  *
  * @param {BoundingSphere} boundingSphere The bounding sphere.
- * @param {Ellipsoid} [ellipsoid=Ellipsoid.WGS84] The ellipsoid.
+ * @param {Ellipsoid} [ellipsoid=Ellipsoid.default] The ellipsoid.
  * @param {Rectangle} [result] The object onto which to store the result, or undefined if a new instance should be created.
  * @returns {Rectangle} The modified result parameter or a new Rectangle instance if none was provided.
  */
@@ -367,7 +367,7 @@ Rectangle.fromBoundingSphere = function (boundingSphere, ellipsoid, result) {
   const radius = boundingSphere.radius;
 
   if (!defined(ellipsoid)) {
-    ellipsoid = Ellipsoid.WGS84;
+    ellipsoid = Ellipsoid.default;
   }
 
   if (!defined(result)) {
@@ -912,7 +912,7 @@ const subsampleLlaScratch = new Cartographic();
  * for rectangles that cover the poles or cross the equator.
  *
  * @param {Rectangle} rectangle The rectangle to subsample.
- * @param {Ellipsoid} [ellipsoid=Ellipsoid.WGS84] The ellipsoid to use.
+ * @param {Ellipsoid} [ellipsoid=Ellipsoid.default] The ellipsoid to use.
  * @param {number} [surfaceHeight=0.0] The height of the rectangle above the ellipsoid.
  * @param {Cartesian3[]} [result] The array of Cartesians onto which to store the result.
  * @returns {Cartesian3[]} The modified result parameter or a new Array of Cartesians instances if none was provided.
@@ -922,7 +922,7 @@ Rectangle.subsample = function (rectangle, ellipsoid, surfaceHeight, result) {
   Check.typeOf.object("rectangle", rectangle);
   //>>includeEnd('debug');
 
-  ellipsoid = defaultValue(ellipsoid, Ellipsoid.WGS84);
+  ellipsoid = defaultValue(ellipsoid, Ellipsoid.default);
   surfaceHeight = defaultValue(surfaceHeight, 0.0);
 
   if (!defined(result)) {

package/Source/Core/RectangleGeometry.js

@@ -978,7 +978,7 @@ function computeRectangle(rectangle, granularity, rotation, ellipsoid, result) {
  * @param {object} options Object with the following properties:
  * @param {Rectangle} options.rectangle A cartographic rectangle with north, south, east and west properties in radians.
  * @param {VertexFormat} [options.vertexFormat=VertexFormat.DEFAULT] The vertex attributes to be computed.
- * @param {Ellipsoid} [options.ellipsoid=Ellipsoid.WGS84] The ellipsoid on which the rectangle lies.
+ * @param {Ellipsoid} [options.ellipsoid=Ellipsoid.default] The ellipsoid on which the rectangle lies.
  * @param {number} [options.granularity=CesiumMath.RADIANS_PER_DEGREE] The distance, in radians, between each latitude and longitude. Determines the number of positions in the buffer.
  * @param {number} [options.height=0.0] The distance in meters between the rectangle and the ellipsoid surface.
  * @param {number} [options.rotation=0.0] The rotation of the rectangle, in radians. A positive rotation is counter-clockwise.
@@ -998,7 +998,7 @@ function computeRectangle(rectangle, granularity, rotation, ellipsoid, result) {
  * @example
  * // 1. create a rectangle
  * const rectangle = new Cesium.RectangleGeometry({
- *   ellipsoid : Cesium.Ellipsoid.WGS84,
+ *   ellipsoid : Cesium.Ellipsoid.default,
  *   rectangle : Cesium.Rectangle.fromDegrees(-80.0, 39.0, -74.0, 42.0),
  *   height : 10000.0
  * });
@@ -1006,7 +1006,7 @@ function computeRectangle(rectangle, granularity, rotation, ellipsoid, result) {
  *
  * // 2. create an extruded rectangle without a top
  * const rectangle = new Cesium.RectangleGeometry({
- *   ellipsoid : Cesium.Ellipsoid.WGS84,
+ *   ellipsoid : Cesium.Ellipsoid.default,
  *   rectangle : Cesium.Rectangle.fromDegrees(-80.0, 39.0, -74.0, 42.0),
  *   height : 10000.0,
  *   extrudedHeight: 300000
@@ -1037,7 +1037,7 @@ function RectangleGeometry(options) {
     CesiumMath.RADIANS_PER_DEGREE
   );
   this._ellipsoid = Ellipsoid.clone(
-    defaultValue(options.ellipsoid, Ellipsoid.WGS84)
+    defaultValue(options.ellipsoid, Ellipsoid.default)
   );
   this._surfaceHeight = Math.max(height, extrudedHeight);
   this._rotation = defaultValue(options.rotation, 0.0);
@@ -1185,7 +1185,7 @@ RectangleGeometry.unpack = function (array, startingIndex, result) {
  *
  * @param {object} options Object with the following properties:
  * @param {Rectangle} options.rectangle A cartographic rectangle with north, south, east and west properties in radians.
- * @param {Ellipsoid} [options.ellipsoid=Ellipsoid.WGS84] The ellipsoid on which the rectangle lies.
+ * @param {Ellipsoid} [options.ellipsoid=Ellipsoid.default] The ellipsoid on which the rectangle lies.
  * @param {number} [options.granularity=CesiumMath.RADIANS_PER_DEGREE] The distance, in radians, between each latitude and longitude. Determines the number of positions in the buffer.
  * @param {number} [options.rotation=0.0] The rotation of the rectangle, in radians. A positive rotation is counter-clockwise.
  * @param {Rectangle} [result] An object in which to store the result.
@@ -1211,7 +1211,7 @@ RectangleGeometry.computeRectangle = function (options, result) {
     options.granularity,
     CesiumMath.RADIANS_PER_DEGREE
   );
-  const ellipsoid = defaultValue(options.ellipsoid, Ellipsoid.WGS84);
+  const ellipsoid = defaultValue(options.ellipsoid, Ellipsoid.default);
   const rotation = defaultValue(options.rotation, 0.0);
 
   return computeRectangle(rectangle, granularity, rotation, ellipsoid, result);

package/Source/Core/RectangleGeometryLibrary.js

@@ -2,6 +2,7 @@ import Cartesian3 from "./Cartesian3.js";
 import Cartographic from "./Cartographic.js";
 import defined from "./defined.js";
 import DeveloperError from "./DeveloperError.js";
+import Ellipsoid from "./Ellipsoid.js";
 import GeographicProjection from "./GeographicProjection.js";
 import CesiumMath from "./Math.js";
 import Matrix2 from "./Matrix2.js";
@@ -100,6 +101,7 @@ function getRotationOptions(
   const granYSin = granularityY * sinRotation;
   const granXSin = granularityX * sinRotation;
 
+  proj._ellipsoid = Ellipsoid.default;
   nwCartesian = proj.project(nwCorner, nwCartesian);
 
   nwCartesian = Cartesian3.subtract(nwCartesian, centerCartesian, nwCartesian);
@@ -190,6 +192,7 @@ RectangleGeometryLibrary.computeOptions = function (
     if (center.longitude < nwCorner.longitude) {
       center.longitude += CesiumMath.TWO_PI;
     }
+    proj._ellipsoid = Ellipsoid.default;
     centerCartesian = proj.project(center, centerCartesian);
   }
 

package/Source/Core/RectangleOutlineGeometry.js

@@ -247,7 +247,7 @@ function constructExtrudedRectangle(rectangleGeometry, computedOptions) {
  *
  * @param {object} options Object with the following properties:
  * @param {Rectangle} options.rectangle A cartographic rectangle with north, south, east and west properties in radians.
- * @param {Ellipsoid} [options.ellipsoid=Ellipsoid.WGS84] The ellipsoid on which the rectangle lies.
+ * @param {Ellipsoid} [options.ellipsoid=Ellipsoid.default] The ellipsoid on which the rectangle lies.
  * @param {number} [options.granularity=CesiumMath.RADIANS_PER_DEGREE] The distance, in radians, between each latitude and longitude. Determines the number of positions in the buffer.
  * @param {number} [options.height=0.0] The distance in meters between the rectangle and the ellipsoid surface.
  * @param {number} [options.rotation=0.0] The rotation of the rectangle, in radians. A positive rotation is counter-clockwise.
@@ -277,7 +277,7 @@ function RectangleOutlineGeometry(options) {
     options.granularity,
     CesiumMath.RADIANS_PER_DEGREE
   );
-  const ellipsoid = defaultValue(options.ellipsoid, Ellipsoid.WGS84);
+  const ellipsoid = defaultValue(options.ellipsoid, Ellipsoid.default);
   const rotation = defaultValue(options.rotation, 0.0);
 
   //>>includeStart('debug', pragmas.debug);

package/Source/Core/Simon1994PlanetaryPositions.js

@@ -613,7 +613,7 @@ function computeSimonEarth(date, result) {
 }
 
 // Values for the <code>axesTransformation</code> needed for the rotation were found using the STK Components
-// GreographicTransformer on the position of the sun center of mass point and the earth J2000 frame.
+// GeographicTransformer on the position of the sun center of mass point and the earth J2000 frame.
 
 const axesTransformation = new Matrix3(
   1.0000000000000002,
@@ -627,6 +627,7 @@ const axesTransformation = new Matrix3(
   0.9174820620691819
 );
 let translation = new Cartesian3();
+
 /**
  * Computes the position of the Sun in the Earth-centered inertial frame
  *
@@ -662,7 +663,7 @@ Simon1994PlanetaryPositions.computeSunPositionInEarthInertialFrame = function (
 /**
  * Computes the position of the Moon in the Earth-centered inertial frame
  *
- * @param {JulianDate} [julianDate] The time at which to compute the Sun's position, if not provided the current system time is used.
+ * @param {JulianDate} [julianDate] The time at which to compute the Moon's position, if not provided the current system time is used.
  * @param {Cartesian3} [result] The object onto which to store the result.
  * @returns {Cartesian3} Calculated moon position
  */
@@ -679,4 +680,5 @@ Simon1994PlanetaryPositions.computeMoonPositionInEarthInertialFrame = function (
 
   return result;
 };
+
 export default Simon1994PlanetaryPositions;

package/Source/Core/SimplePolylineGeometry.js

@@ -68,7 +68,7 @@ function interpolateColors(p0, p1, color0, color1, minDistance, array, offset) {
  * @param {boolean} [options.colorsPerVertex=false] A boolean that determines whether the colors will be flat across each segment of the line or interpolated across the vertices.
  * @param {ArcType} [options.arcType=ArcType.GEODESIC] The type of line the polyline segments must follow.
  * @param {number} [options.granularity=CesiumMath.RADIANS_PER_DEGREE] The distance, in radians, between each latitude and longitude if options.arcType is not ArcType.NONE. Determines the number of positions in the buffer.
- * @param {Ellipsoid} [options.ellipsoid=Ellipsoid.WGS84] The ellipsoid to be used as a reference.
+ * @param {Ellipsoid} [options.ellipsoid=Ellipsoid.default] The ellipsoid to be used as a reference.
  *
  * @exception {DeveloperError} At least two positions are required.
  * @exception {DeveloperError} colors has an invalid length.
@@ -114,7 +114,7 @@ function SimplePolylineGeometry(options) {
     options.granularity,
     CesiumMath.RADIANS_PER_DEGREE
   );
-  this._ellipsoid = defaultValue(options.ellipsoid, Ellipsoid.WGS84);
+  this._ellipsoid = defaultValue(options.ellipsoid, Ellipsoid.default);
   this._workerName = "createSimplePolylineGeometry";
 
   let numComponents = 1 + positions.length * Cartesian3.packedLength;

package/Source/Core/Stereographic.js

@@ -99,12 +99,12 @@ const scratchCartesian = new Cartesian3();
 /**
  * Computes the latitude based on an ellipsoid.
  *
- * @param {Ellipsoid} [ellipsoid=Ellipsoid.WGS84] The ellipsoid on which to compute the longitude.
+ * @param {Ellipsoid} [ellipsoid=Ellipsoid.default] The ellipsoid on which to compute the longitude.
  * @returns {number} The latitude
  */
 Stereographic.prototype.getLatitude = function (ellipsoid) {
   if (!defined(ellipsoid)) {
-    ellipsoid = Ellipsoid.WGS84;
+    ellipsoid = Ellipsoid.default;
   }
 
   scratchCartographic.latitude = this.conformalLatitude;

package/Source/Core/Transforms.js

@@ -113,7 +113,7 @@ Transforms.localFrameToFixedFrameGenerator = function (firstAxis, secondAxis) {
    * centered at the provided origin to the provided ellipsoid's fixed reference frame.
    * @callback Transforms.LocalFrameToFixedFrame
    * @param {Cartesian3} origin The center point of the local reference frame.
-   * @param {Ellipsoid} [ellipsoid=Ellipsoid.WGS84] The ellipsoid whose fixed frame is used in the transformation.
+   * @param {Ellipsoid} [ellipsoid=Ellipsoid.default] The ellipsoid whose fixed frame is used in the transformation.
    * @param {Matrix4} [result] The object onto which to store the result.
    * @returns {Matrix4} The modified result parameter or a new Matrix4 instance if none was provided.
    */
@@ -199,7 +199,7 @@ Transforms.localFrameToFixedFrameGenerator = function (firstAxis, secondAxis) {
           );
         }
       } else {
-        ellipsoid = defaultValue(ellipsoid, Ellipsoid.WGS84);
+        ellipsoid = defaultValue(ellipsoid, Ellipsoid.default);
         ellipsoid.geodeticSurfaceNormal(origin, scratchCalculateCartesian.up);
 
         const up = scratchCalculateCartesian.up;
@@ -265,7 +265,7 @@ Transforms.localFrameToFixedFrameGenerator = function (firstAxis, secondAxis) {
  *
  * @function
  * @param {Cartesian3} origin The center point of the local reference frame.
- * @param {Ellipsoid} [ellipsoid=Ellipsoid.WGS84] The ellipsoid whose fixed frame is used in the transformation.
+ * @param {Ellipsoid} [ellipsoid=Ellipsoid.default] The ellipsoid whose fixed frame is used in the transformation.
  * @param {Matrix4} [result] The object onto which to store the result.
  * @returns {Matrix4} The modified result parameter or a new Matrix4 instance if none was provided.
  *
@@ -291,7 +291,7 @@ Transforms.eastNorthUpToFixedFrame = Transforms.localFrameToFixedFrameGenerator(
  *
  * @function
  * @param {Cartesian3} origin The center point of the local reference frame.
- * @param {Ellipsoid} [ellipsoid=Ellipsoid.WGS84] The ellipsoid whose fixed frame is used in the transformation.
+ * @param {Ellipsoid} [ellipsoid=Ellipsoid.default] The ellipsoid whose fixed frame is used in the transformation.
  * @param {Matrix4} [result] The object onto which to store the result.
  * @returns {Matrix4} The modified result parameter or a new Matrix4 instance if none was provided.
  *
@@ -317,7 +317,7 @@ Transforms.northEastDownToFixedFrame = Transforms.localFrameToFixedFrameGenerato
  *
  * @function
  * @param {Cartesian3} origin The center point of the local reference frame.
- * @param {Ellipsoid} [ellipsoid=Ellipsoid.WGS84] The ellipsoid whose fixed frame is used in the transformation.
+ * @param {Ellipsoid} [ellipsoid=Ellipsoid.default] The ellipsoid whose fixed frame is used in the transformation.
  * @param {Matrix4} [result] The object onto which to store the result.
  * @returns {Matrix4} The modified result parameter or a new Matrix4 instance if none was provided.
  *
@@ -343,7 +343,7 @@ Transforms.northUpEastToFixedFrame = Transforms.localFrameToFixedFrameGenerator(
  *
  * @function
  * @param {Cartesian3} origin The center point of the local reference frame.
- * @param {Ellipsoid} [ellipsoid=Ellipsoid.WGS84] The ellipsoid whose fixed frame is used in the transformation.
+ * @param {Ellipsoid} [ellipsoid=Ellipsoid.default] The ellipsoid whose fixed frame is used in the transformation.
  * @param {Matrix4} [result] The object onto which to store the result.
  * @returns {Matrix4} The modified result parameter or a new Matrix4 instance if none was provided.
  *
@@ -369,7 +369,7 @@ const scratchHPRMatrix4 = new Matrix4();
  *
  * @param {Cartesian3} origin The center point of the local reference frame.
  * @param {HeadingPitchRoll} headingPitchRoll The heading, pitch, and roll.
- * @param {Ellipsoid} [ellipsoid=Ellipsoid.WGS84] The ellipsoid whose fixed frame is used in the transformation.
+ * @param {Ellipsoid} [ellipsoid=Ellipsoid.default] The ellipsoid whose fixed frame is used in the transformation.
  * @param {Transforms.LocalFrameToFixedFrame} [fixedFrameTransform=Transforms.eastNorthUpToFixedFrame] A 4x4 transformation
  *  matrix from a reference frame to the provided ellipsoid's fixed reference frame
  * @param {Matrix4} [result] The object onto which to store the result.
@@ -424,7 +424,7 @@ const scratchHPRMatrix3 = new Matrix3();
  *
  * @param {Cartesian3} origin The center point of the local reference frame.
  * @param {HeadingPitchRoll} headingPitchRoll The heading, pitch, and roll.
- * @param {Ellipsoid} [ellipsoid=Ellipsoid.WGS84] The ellipsoid whose fixed frame is used in the transformation.
+ * @param {Ellipsoid} [ellipsoid=Ellipsoid.default] The ellipsoid whose fixed frame is used in the transformation.
  * @param {Transforms.LocalFrameToFixedFrame} [fixedFrameTransform=Transforms.eastNorthUpToFixedFrame] A 4x4 transformation
  *  matrix from a reference frame to the provided ellipsoid's fixed reference frame
  * @param {Quaternion} [result] The object onto which to store the result.
@@ -473,7 +473,7 @@ const hprQuaternionScratch = new Quaternion();
  * are above the plane. Negative pitch angles are below the plane. Roll is the first rotation applied about the local east axis.
  *
  * @param {Matrix4} transform The transform
- * @param {Ellipsoid} [ellipsoid=Ellipsoid.WGS84] The ellipsoid whose fixed frame is used in the transformation.
+ * @param {Ellipsoid} [ellipsoid=Ellipsoid.default] The ellipsoid whose fixed frame is used in the transformation.
  * @param {Transforms.LocalFrameToFixedFrame} [fixedFrameTransform=Transforms.eastNorthUpToFixedFrame] A 4x4 transformation
  *  matrix from a reference frame to the provided ellipsoid's fixed reference frame
  * @param {HeadingPitchRoll} [result] The object onto which to store the result.
@@ -489,7 +489,7 @@ Transforms.fixedFrameToHeadingPitchRoll = function (
   Check.defined("transform", transform);
   //>>includeEnd('debug');
 
-  ellipsoid = defaultValue(ellipsoid, Ellipsoid.WGS84);
+  ellipsoid = defaultValue(ellipsoid, Ellipsoid.default);
   fixedFrameTransform = defaultValue(
     fixedFrameTransform,
     Transforms.eastNorthUpToFixedFrame
@@ -538,6 +538,35 @@ const wgs84WRPrecessing = 7.2921158553e-5;
 const twoPiOverSecondsInDay = CesiumMath.TWO_PI / 86400.0;
 let dateInUtc = new JulianDate();
 
+/**
+ * The default function to compute a rotation matrix to transform a point or vector from the International Celestial
+ * Reference Frame (GCRF/ICRF) inertial frame axes to the central body, typically Earth, fixed frame axis at a given
+ * time for use in lighting and transformation from inertial reference frames. This function may return undefined if
+ * the data necessary to do the transformation is not yet loaded.
+ *
+ * @param {JulianDate} date The time at which to compute the rotation matrix.
+ * @param {Matrix3} [result] The object onto which to store the result.  If this parameter is
+ *                  not specified, a new instance is created and returned.
+ * @returns {Matrix3|undefined} The rotation matrix, or undefined if the data necessary to do the
+ *                   transformation is not yet loaded.
+ *
+ * @example
+ * // Set the default ICRF to fixed transformation to that of the Moon.
+ * Cesium.Transforms.computeIcrfToCentralBodyFixedMatrix = Cesium.Transforms.computeIcrfToMoonFixedMatrix;
+ *
+ * @see Transforms.computeIcrfToFixedMatrix
+ * @see Transforms.computeTemeToPseudoFixedMatrix
+ * @see Transforms.computeIcrfToMoonFixedMatrix
+ */
+Transforms.computeIcrfToCentralBodyFixedMatrix = function (date, result) {
+  let transformMatrix = Transforms.computeIcrfToFixedMatrix(date, result);
+  if (!defined(transformMatrix)) {
+    transformMatrix = Transforms.computeTemeToPseudoFixedMatrix(date, result);
+  }
+
+  return transformMatrix;
+};
+
 /**
  * Computes a rotation matrix to transform a point or vector from True Equator Mean Equinox (TEME) axes to the
  * pseudo-fixed axes at a given time.  This method treats the UT1 time standard as equivalent to UTC.
@@ -691,7 +720,7 @@ Transforms.preloadIcrfFixed = function (timeInterval) {
  * @param {JulianDate} date The time at which to compute the rotation matrix.
  * @param {Matrix3} [result] The object onto which to store the result.  If this parameter is
  *                  not specified, a new instance is created and returned.
- * @returns {Matrix3} The rotation matrix, or undefined if the data necessary to do the
+ * @returns {Matrix3|undefined} The rotation matrix, or undefined if the data necessary to do the
  *                   transformation is not yet loaded.
  *
  *
@@ -726,6 +755,112 @@ Transforms.computeIcrfToFixedMatrix = function (date, result) {
   return Matrix3.transpose(fixedToIcrfMtx, result);
 };
 
+const TdtMinusTai = 32.184;
+const J2000d = 2451545;
+const scratchHpr = new HeadingPitchRoll();
+const scratchRotationMatrix = new Matrix3();
+const dateScratch = new JulianDate();
+
+/**
+ * Computes a rotation matrix to transform a point or vector from the Moon-Fixed frame axes
+ * to the International Celestial Reference Frame (GCRF/ICRF) inertial frame axes
+ * at a given time.
+ *
+ * @param {JulianDate} date The time at which to compute the rotation matrix.
+ * @param {Matrix3} [result] The object onto which to store the result.  If this parameter is
+ *                  not specified, a new instance is created and returned.
+ * @returns {Matrix3} The rotation matrix.
+ *
+ * @example
+ * // Transform a point from the Fixed axes to the ICRF axes.
+ * const now = Cesium.JulianDate.now();
+ * const pointInFixed = Cesium.Cartesian3.fromDegrees(0.0, 0.0);
+ * const fixedToIcrf = Cesium.Transforms.computeMoonFixedToIcrfMatrix(now);
+ * let pointInInertial = new Cesium.Cartesian3();
+ * if (Cesium.defined(fixedToIcrf)) {
+ *     pointInInertial = Cesium.Matrix3.multiplyByVector(fixedToIcrf, pointInFixed, pointInInertial);
+ * }
+ */
+Transforms.computeMoonFixedToIcrfMatrix = function (date, result) {
+  //>>includeStart('debug', pragmas.debug);
+  if (!defined(date)) {
+    throw new DeveloperError("date is required.");
+  }
+  //>>includeEnd('debug');
+
+  if (!defined(result)) {
+    result = new Matrix3();
+  }
+
+  // Converts TAI to TT
+  const secondsTT = JulianDate.addSeconds(date, TdtMinusTai, dateScratch);
+
+  // Converts TT to TDB, interval in days since the standard epoch
+  const d = JulianDate.totalDays(secondsTT) - J2000d;
+
+  // Compute the approximate rotation, using https://articles.adsabs.harvard.edu//full/1980CeMec..22..205D/0000209.000.html
+  const e1 = CesiumMath.toRadians(12.112) - CesiumMath.toRadians(0.052992) * d;
+  const e2 = CesiumMath.toRadians(24.224) - CesiumMath.toRadians(0.105984) * d;
+  const e3 = CesiumMath.toRadians(227.645) + CesiumMath.toRadians(13.012) * d;
+  const e4 =
+    CesiumMath.toRadians(261.105) + CesiumMath.toRadians(13.340716) * d;
+  const e5 = CesiumMath.toRadians(358.0) + CesiumMath.toRadians(0.9856) * d;
+
+  scratchHpr.pitch =
+    CesiumMath.toRadians(270.0 - 90) -
+    CesiumMath.toRadians(3.878) * Math.sin(e1) -
+    CesiumMath.toRadians(0.12) * Math.sin(e2) +
+    CesiumMath.toRadians(0.07) * Math.sin(e3) -
+    CesiumMath.toRadians(0.017) * Math.sin(e4);
+  scratchHpr.roll =
+    CesiumMath.toRadians(66.53 - 90) +
+    CesiumMath.toRadians(1.543) * Math.cos(e1) +
+    CesiumMath.toRadians(0.24) * Math.cos(e2) -
+    CesiumMath.toRadians(0.028) * Math.cos(e3) +
+    CesiumMath.toRadians(0.007) * Math.cos(e4);
+  scratchHpr.heading =
+    CesiumMath.toRadians(244.375 - 90) +
+    CesiumMath.toRadians(13.17635831) * d +
+    CesiumMath.toRadians(3.558) * Math.sin(e1) +
+    CesiumMath.toRadians(0.121) * Math.sin(e2) -
+    CesiumMath.toRadians(0.064) * Math.sin(e3) +
+    CesiumMath.toRadians(0.016) * Math.sin(e4) +
+    CesiumMath.toRadians(0.025) * Math.sin(e5);
+  return Matrix3.fromHeadingPitchRoll(scratchHpr, scratchRotationMatrix);
+};
+
+/**
+ * Computes a rotation matrix to transform a point or vector from the International Celestial
+ * Reference Frame (GCRF/ICRF) inertial frame axes to the Moon-Fixed frame axes
+ * at a given time.
+ *
+ * @param {JulianDate} date The time at which to compute the rotation matrix.
+ * @param {Matrix3} [result] The object onto which to store the result.  If this parameter is
+ *                  not specified, a new instance is created and returned.
+ * @returns {Matrix3} The rotation matrix.
+ *
+ * @example
+ * // Set the default ICRF to fixed transformation to that of the Moon.
+ * Cesium.Transforms.computeIcrfToCentralBodyFixedMatrix = Cesium.Transforms.computeIcrfToMoonFixedMatrix;
+ */
+Transforms.computeIcrfToMoonFixedMatrix = function (date, result) {
+  //>>includeStart('debug', pragmas.debug);
+  if (!defined(date)) {
+    throw new DeveloperError("date is required.");
+  }
+  //>>includeEnd('debug');
+  if (!defined(result)) {
+    result = new Matrix3();
+  }
+
+  const fixedToIcrfMtx = Transforms.computeMoonFixedToIcrfMatrix(date, result);
+  if (!defined(fixedToIcrfMtx)) {
+    return undefined;
+  }
+
+  return Matrix3.transpose(fixedToIcrfMtx, result);
+};
+
 const xysScratch = new Iau2006XysSample(0.0, 0.0, 0.0);
 const eopScratch = new EarthOrientationParametersSample(
   0.0,
@@ -747,7 +882,7 @@ const rotation2Scratch = new Matrix3();
  * @param {JulianDate} date The time at which to compute the rotation matrix.
  * @param {Matrix3} [result] The object onto which to store the result.  If this parameter is
  *                  not specified, a new instance is created and returned.
- * @returns {Matrix3} The rotation matrix, or undefined if the data necessary to do the
+ * @returns {Matrix3|undefined} The rotation matrix, or undefined if the data necessary to do the
  *                   transformation is not yet loaded.
  *
  *
@@ -950,7 +1085,7 @@ const upScratch = new Cartesian3();
  *
  * @param {Cartesian3} position The position to transform.
  * @param {Cartesian3} velocity The velocity vector to transform.
- * @param {Ellipsoid} [ellipsoid=Ellipsoid.WGS84] The ellipsoid whose fixed frame is used in the transformation.
+ * @param {Ellipsoid} [ellipsoid=Ellipsoid.default] The ellipsoid whose fixed frame is used in the transformation.
  * @param {Matrix3} [result] The object onto which to store the result.
  * @returns {Matrix3} The modified result parameter or a new Matrix3 instance if none was provided.
  */
@@ -970,10 +1105,10 @@ Transforms.rotationMatrixFromPositionVelocity = function (
   }
   //>>includeEnd('debug');
 
-  const normal = defaultValue(ellipsoid, Ellipsoid.WGS84).geodeticSurfaceNormal(
-    position,
-    normalScratch
-  );
+  const normal = defaultValue(
+    ellipsoid,
+    Ellipsoid.default
+  ).geodeticSurfaceNormal(position, normalScratch);
   let right = Cartesian3.cross(velocity, normal, rightScratch);
 
   if (Cartesian3.equalsEpsilon(right, Cartesian3.ZERO, CesiumMath.EPSILON6)) {
@@ -1073,7 +1208,7 @@ Transforms.basisTo2D = function (projection, matrix, result) {
     );
   }
 
-  // Assuming the instance are positioned in WGS84, invert the WGS84 transform to get the local transform and then convert to 2D
+  // Assuming the instance are positioned on the ellipsoid, invert the ellipsoidal transform to get the local transform and then convert to 2D
   const fromENU = Transforms.eastNorthUpToFixedFrame(
     rtcCenter,
     ellipsoid,
@@ -1091,7 +1226,7 @@ Transforms.basisTo2D = function (projection, matrix, result) {
 /**
  * @private
  */
-Transforms.wgs84To2DModelMatrix = function (projection, center, result) {
+Transforms.ellipsoidTo2DModelMatrix = function (projection, center, result) {
   //>>includeStart('debug', pragmas.debug);
   if (!defined(projection)) {
     throw new DeveloperError("projection is required.");

package/Source/Core/VRTheWorldTerrainProvider.js

@@ -24,7 +24,7 @@ function DataRectangle(rectangle, maxLevel) {
  *
  * Initialization options for the VRTheWorldTerrainProvider constructor
  *
- * @property {Ellipsoid} [ellipsoid] The ellipsoid.  If not specified, the WGS84 ellipsoid is used.
+ * @property {Ellipsoid} [ellipsoid=Ellipsoid.default] The ellipsoid.  If not specified, the default ellipsoid is used.
  * @property {Credit|string} [credit] A credit for the data source, which is displayed on the canvas.
  */
 
@@ -37,7 +37,7 @@ function DataRectangle(rectangle, maxLevel) {
  * @param {VRTheWorldTerrainProvider.ConstructorOptions} options An object describing initialization options
  */
 function TerrainProviderBuilder(options) {
-  this.ellipsoid = defaultValue(options.ellipsoid, Ellipsoid.WGS84);
+  this.ellipsoid = defaultValue(options.ellipsoid, Ellipsoid.default);
   this.tilingScheme = undefined;
   this.heightmapWidth = undefined;
   this.heightmapHeight = undefined;

package/Source/Core/WallGeometry.js

@@ -36,7 +36,7 @@ const scratchNormal = new Cartesian3();
  *        wall at <code>positions</code>. If undefined, the height of each position in used.
  * @param {number[]} [options.minimumHeights] An array parallel to <code>positions</code> that give the minimum height of the
  *        wall at <code>positions</code>. If undefined, the height at each position is 0.0.
- * @param {Ellipsoid} [options.ellipsoid=Ellipsoid.WGS84] The ellipsoid for coordinate manipulation
+ * @param {Ellipsoid} [options.ellipsoid=Ellipsoid.default] The ellipsoid for coordinate manipulation
  * @param {VertexFormat} [options.vertexFormat=VertexFormat.DEFAULT] The vertex attributes to be computed.
  *
  * @exception {DeveloperError} positions length must be greater than or equal to 2.
@@ -95,7 +95,7 @@ function WallGeometry(options) {
     options.granularity,
     CesiumMath.RADIANS_PER_DEGREE
   );
-  const ellipsoid = defaultValue(options.ellipsoid, Ellipsoid.WGS84);
+  const ellipsoid = defaultValue(options.ellipsoid, Ellipsoid.default);
 
   this._positions = wallPositions;
   this._minimumHeights = minimumHeights;
@@ -280,7 +280,7 @@ WallGeometry.unpack = function (array, startingIndex, result) {
  *        wall at <code>positions</code>. If undefined, the height of each position in used.
  * @param {number} [options.minimumHeight] A constant that defines the minimum height of the
  *        wall at <code>positions</code>. If undefined, the height at each position is 0.0.
- * @param {Ellipsoid} [options.ellipsoid=Ellipsoid.WGS84] The ellipsoid for coordinate manipulation
+ * @param {Ellipsoid} [options.ellipsoid=Ellipsoid.default] The ellipsoid for coordinate manipulation
  * @param {VertexFormat} [options.vertexFormat=VertexFormat.DEFAULT] The vertex attributes to be computed.
  * @returns {WallGeometry}
  *

package/Source/Core/WallOutlineGeometry.js

@@ -30,7 +30,7 @@ const scratchCartesian3Position2 = new Cartesian3();
  *        wall at <code>positions</code>. If undefined, the height of each position in used.
  * @param {number[]} [options.minimumHeights] An array parallel to <code>positions</code> that give the minimum height of the
  *        wall at <code>positions</code>. If undefined, the height at each position is 0.0.
- * @param {Ellipsoid} [options.ellipsoid=Ellipsoid.WGS84] The ellipsoid for coordinate manipulation
+ * @param {Ellipsoid} [options.ellipsoid=Ellipsoid.default] The ellipsoid for coordinate manipulation
  *
  * @exception {DeveloperError} positions length must be greater than or equal to 2.
  * @exception {DeveloperError} positions and maximumHeights must have the same length.
@@ -85,7 +85,7 @@ function WallOutlineGeometry(options) {
     options.granularity,
     CesiumMath.RADIANS_PER_DEGREE
   );
-  const ellipsoid = defaultValue(options.ellipsoid, Ellipsoid.WGS84);
+  const ellipsoid = defaultValue(options.ellipsoid, Ellipsoid.default);
 
   this._positions = wallPositions;
   this._minimumHeights = minimumHeights;
@@ -255,7 +255,7 @@ WallOutlineGeometry.unpack = function (array, startingIndex, result) {
  *        wall at <code>positions</code>. If undefined, the height of each position in used.
  * @param {number} [options.minimumHeight] A constant that defines the minimum height of the
  *        wall at <code>positions</code>. If undefined, the height at each position is 0.0.
- * @param {Ellipsoid} [options.ellipsoid=Ellipsoid.WGS84] The ellipsoid for coordinate manipulation
+ * @param {Ellipsoid} [options.ellipsoid=Ellipsoid.default] The ellipsoid for coordinate manipulation
  * @returns {WallOutlineGeometry}
  *
  *

package/Source/Core/WebMercatorTilingScheme.js

@@ -13,8 +13,8 @@ import WebMercatorProjection from "./WebMercatorProjection.js";
  * @constructor
  *
  * @param {object} [options] Object with the following properties:
- * @param {Ellipsoid} [options.ellipsoid=Ellipsoid.WGS84] The ellipsoid whose surface is being tiled. Defaults to
- * the WGS84 ellipsoid.
+ * @param {Ellipsoid} [options.ellipsoid=Ellipsoid.default] The ellipsoid whose surface is being tiled. Defaults to
+ * the default ellipsoid.
  * @param {number} [options.numberOfLevelZeroTilesX=1] The number of tiles in the X direction at level zero of
  *        the tile tree.
  * @param {number} [options.numberOfLevelZeroTilesY=1] The number of tiles in the Y direction at level zero of
@@ -31,7 +31,7 @@ import WebMercatorProjection from "./WebMercatorProjection.js";
 function WebMercatorTilingScheme(options) {
   options = defaultValue(options, defaultValue.EMPTY_OBJECT);
 
-  this._ellipsoid = defaultValue(options.ellipsoid, Ellipsoid.WGS84);
+  this._ellipsoid = defaultValue(options.ellipsoid, Ellipsoid.default);
   this._numberOfLevelZeroTilesX = defaultValue(
     options.numberOfLevelZeroTilesX,
     1

package/Source/Core/createWorldTerrainAsync.js

@@ -1,5 +1,6 @@
 import CesiumTerrainProvider from "./CesiumTerrainProvider.js";
 import defaultValue from "./defaultValue.js";
+import Ellipsoid from "./Ellipsoid.js";
 
 /**
  * Creates a {@link CesiumTerrainProvider} instance for the {@link https://cesium.com/content/#cesium-world-terrain|Cesium World Terrain}.
@@ -43,6 +44,7 @@ function createWorldTerrainAsync(options) {
   return CesiumTerrainProvider.fromIonAssetId(1, {
     requestVertexNormals: defaultValue(options.requestVertexNormals, false),
     requestWaterMask: defaultValue(options.requestWaterMask, false),
+    ellipsoid: Ellipsoid.WGS84,
   });
 }
 export default createWorldTerrainAsync;