Source: Core/EllipseGeometry.js

  1. /*global define*/
  2. define([
  3.         './BoundingSphere',
  4.         './Cartesian2',
  5.         './Cartesian3',
  6.         './Cartographic',
  7.         './ComponentDatatype',
  8.         './defaultValue',
  9.         './defined',
  10.         './defineProperties',
  11.         './DeveloperError',
  12.         './EllipseGeometryLibrary',
  13.         './Ellipsoid',
  14.         './GeographicProjection',
  15.         './Geometry',
  16.         './GeometryAttribute',
  17.         './GeometryAttributes',
  18.         './GeometryInstance',
  19.         './GeometryPipeline',
  20.         './IndexDatatype',
  21.         './Math',
  22.         './Matrix3',
  23.         './Matrix4',
  24.         './PrimitiveType',
  25.         './Quaternion',
  26.         './Rectangle',
  27.         './Transforms',
  28.         './VertexFormat'
  29.     ], function(
  30.         BoundingSphere,
  31.         Cartesian2,
  32.         Cartesian3,
  33.         Cartographic,
  34.         ComponentDatatype,
  35.         defaultValue,
  36.         defined,
  37.         defineProperties,
  38.         DeveloperError,
  39.         EllipseGeometryLibrary,
  40.         Ellipsoid,
  41.         GeographicProjection,
  42.         Geometry,
  43.         GeometryAttribute,
  44.         GeometryAttributes,
  45.         GeometryInstance,
  46.         GeometryPipeline,
  47.         IndexDatatype,
  48.         CesiumMath,
  49.         Matrix3,
  50.         Matrix4,
  51.         PrimitiveType,
  52.         Quaternion,
  53.         Rectangle,
  54.         Transforms,
  55.         VertexFormat) {
  56.     'use strict';
  58.     var scratchCartesian1 = new Cartesian3();
  59.     var scratchCartesian2 = new Cartesian3();
  60.     var scratchCartesian3 = new Cartesian3();
  61.     var scratchCartesian4 = new Cartesian3();
  62.     var texCoordScratch = new Cartesian2();
  63.     var textureMatrixScratch = new Matrix3();
  64.     var quaternionScratch = new Quaternion();
  66.     var scratchNormal = new Cartesian3();
  67.     var scratchTangent = new Cartesian3();
  68.     var scratchBinormal = new Cartesian3();
  70.     var scratchCartographic = new Cartographic();
  71.     var projectedCenterScratch = new Cartesian3();
  73.     var scratchMinTexCoord = new Cartesian2();
  74.     var scratchMaxTexCoord = new Cartesian2();
  76.     function computeTopBottomAttributes(positions, options, extrude) {
  77.         var vertexFormat = options.vertexFormat;
  78.         var center = options.center;
  79.         var semiMajorAxis = options.semiMajorAxis;
  80.         var semiMinorAxis = options.semiMinorAxis;
  81.         var ellipsoid = options.ellipsoid;
  82.         var stRotation = options.stRotation;
  83.         var size = (extrude) ? positions.length / 3 * 2 : positions.length / 3;
  85.         var textureCoordinates = (vertexFormat.st) ? new Float32Array(size * 2) : undefined;
  86.         var normals = (vertexFormat.normal) ? new Float32Array(size * 3) : undefined;
  87.         var tangents = (vertexFormat.tangent) ? new Float32Array(size * 3) : undefined;
  88.         var binormals = (vertexFormat.binormal) ? new Float32Array(size * 3) : undefined;
  90.         var textureCoordIndex = 0;
  92.         // Raise positions to a height above the ellipsoid and compute the
  93.         // texture coordinates, normals, tangents, and binormals.
  94.         var normal = scratchNormal;
  95.         var tangent = scratchTangent;
  96.         var binormal = scratchBinormal;
  98.         var projection = new GeographicProjection(ellipsoid);
  99.         var projectedCenter = projection.project(ellipsoid.cartesianToCartographic(center, scratchCartographic), projectedCenterScratch);
  101.         var geodeticNormal = ellipsoid.scaleToGeodeticSurface(center, scratchCartesian1);
  102.         ellipsoid.geodeticSurfaceNormal(geodeticNormal, geodeticNormal);
  103.         var rotation = Quaternion.fromAxisAngle(geodeticNormal, stRotation, quaternionScratch);
  104.         var textureMatrix = Matrix3.fromQuaternion(rotation, textureMatrixScratch);
  106.         var minTexCoord = Cartesian2.fromElements(Number.POSITIVE_INFINITY, Number.POSITIVE_INFINITY, scratchMinTexCoord);
  107.         var maxTexCoord = Cartesian2.fromElements(Number.NEGATIVE_INFINITY, Number.NEGATIVE_INFINITY, scratchMaxTexCoord);
  109.         var length = positions.length;
  110.         var bottomOffset = (extrude) ? length : 0;
  111.         var stOffset = bottomOffset / 3 * 2;
  112.         for (var i = 0; i < length; i += 3) {
  113.             var i1 = i + 1;
  114.             var i2 = i + 2;
  115.             var position = Cartesian3.fromArray(positions, i, scratchCartesian1);
  117.             if (vertexFormat.st) {
  118.                 var rotatedPoint = Matrix3.multiplyByVector(textureMatrix, position, scratchCartesian2);
  119.                 var projectedPoint = projection.project(ellipsoid.cartesianToCartographic(rotatedPoint, scratchCartographic), scratchCartesian3);
  120.                 Cartesian3.subtract(projectedPoint, projectedCenter, projectedPoint);
  122.                 texCoordScratch.x = (projectedPoint.x + semiMajorAxis) / (2.0 * semiMajorAxis);
  123.                 texCoordScratch.y = (projectedPoint.y + semiMinorAxis) / (2.0 * semiMinorAxis);
  125.                 minTexCoord.x = Math.min(texCoordScratch.x, minTexCoord.x);
  126.                 minTexCoord.y = Math.min(texCoordScratch.y, minTexCoord.y);
  127.                 maxTexCoord.x = Math.max(texCoordScratch.x, maxTexCoord.x);
  128.                 maxTexCoord.y = Math.max(texCoordScratch.y, maxTexCoord.y);
  130.                 if (extrude) {
  131.                     textureCoordinates[textureCoordIndex + stOffset] = texCoordScratch.x;
  132.                     textureCoordinates[textureCoordIndex + 1 + stOffset] = texCoordScratch.y;
  133.                 }
  135.                 textureCoordinates[textureCoordIndex++] = texCoordScratch.x;
  136.                 textureCoordinates[textureCoordIndex++] = texCoordScratch.y;
  137.             }
  139.             normal = ellipsoid.geodeticSurfaceNormal(position, normal);
  141.             if (vertexFormat.normal || vertexFormat.tangent || vertexFormat.binormal) {
  142.                 if (vertexFormat.tangent || vertexFormat.binormal) {
  143.                     tangent = Cartesian3.normalize(Cartesian3.cross(Cartesian3.UNIT_Z, normal, tangent), tangent);
  144.                     Matrix3.multiplyByVector(textureMatrix, tangent, tangent);
  145.                 }
  146.                 if (vertexFormat.normal) {
  147.                     normals[i] = normal.x;
  148.                     normals[i1] = normal.y;
  149.                     normals[i2] = normal.z;
  150.                     if (extrude) {
  151.                         normals[i + bottomOffset] = -normal.x;
  152.                         normals[i1 + bottomOffset] = -normal.y;
  153.                         normals[i2 + bottomOffset] = -normal.z;
  154.                     }
  155.                 }
  157.                 if (vertexFormat.tangent) {
  158.                     tangents[i] = tangent.x;
  159.                     tangents[i1] = tangent.y;
  160.                     tangents[i2] = tangent.z;
  161.                     if (extrude) {
  162.                         tangents[i + bottomOffset] = -tangent.x;
  163.                         tangents[i1 + bottomOffset] = -tangent.y;
  164.                         tangents[i2 + bottomOffset] = -tangent.z;
  165.                     }
  166.                 }
  168.                 if (vertexFormat.binormal) {
  169.                     binormal = Cartesian3.normalize(Cartesian3.cross(normal, tangent, binormal), binormal);
  170.                     binormals[i] = binormal.x;
  171.                     binormals[i1] = binormal.y;
  172.                     binormals[i2] = binormal.z;
  173.                     if (extrude) {
  174.                         binormals[i + bottomOffset] = binormal.x;
  175.                         binormals[i1 + bottomOffset] = binormal.y;
  176.                         binormals[i2 + bottomOffset] = binormal.z;
  177.                     }
  178.                 }
  179.             }
  180.         }
  182.         if (vertexFormat.st) {
  183.             length = textureCoordinates.length;
  184.             for (var k = 0; k < length; k += 2) {
  185.                 textureCoordinates[k] = (textureCoordinates[k] - minTexCoord.x) / (maxTexCoord.x - minTexCoord.x);
  186.                 textureCoordinates[k + 1] = (textureCoordinates[k + 1] - minTexCoord.y) / (maxTexCoord.y - minTexCoord.y);
  187.             }
  188.         }
  190.         var attributes = new GeometryAttributes();
  192.         if (vertexFormat.position) {
  193.             var finalPositions = EllipseGeometryLibrary.raisePositionsToHeight(positions, options, extrude);
  194.             attributes.position = new GeometryAttribute({
  195.                 componentDatatype : ComponentDatatype.DOUBLE,
  196.                 componentsPerAttribute : 3,
  197.                 values : finalPositions
  198.             });
  199.         }
  201.         if (vertexFormat.st) {
  202.             attributes.st = new GeometryAttribute({
  203.                 componentDatatype : ComponentDatatype.FLOAT,
  204.                 componentsPerAttribute : 2,
  205.                 values : textureCoordinates
  206.             });
  207.         }
  209.         if (vertexFormat.normal) {
  210.             attributes.normal = new GeometryAttribute({
  211.                 componentDatatype : ComponentDatatype.FLOAT,
  212.                 componentsPerAttribute : 3,
  213.                 values : normals
  214.             });
  215.         }
  217.         if (vertexFormat.tangent) {
  218.             attributes.tangent = new GeometryAttribute({
  219.                 componentDatatype : ComponentDatatype.FLOAT,
  220.                 componentsPerAttribute : 3,
  221.                 values : tangents
  222.             });
  223.         }
  225.         if (vertexFormat.binormal) {
  226.             attributes.binormal = new GeometryAttribute({
  227.                 componentDatatype : ComponentDatatype.FLOAT,
  228.                 componentsPerAttribute : 3,
  229.                 values : binormals
  230.             });
  231.         }
  232.         return attributes;
  233.     }
  235.     function topIndices(numPts) {
  236.         // numTriangles in half = 3 + 8 + 12 + ... = -1 + 4 + (4 + 4) + (4 + 4 + 4) + ... = -1 + 4 * (1 + 2 + 3 + ...)
  237.         //              = -1 + 4 * ((n * ( n + 1)) / 2)
  238.         // total triangles = 2 * numTrangles in half
  239.         // indices = total triangles * 3;
  240.         // Substitute numPts for n above
  242.         var indices = new Array(12 * (numPts * ( numPts + 1)) - 6);
  243.         var indicesIndex = 0;
  244.         var prevIndex;
  245.         var numInterior;
  246.         var positionIndex;
  247.         var i;
  248.         var j;
  249.         // Indices triangles to the 'right' of the north vector
  251.         prevIndex = 0;
  252.         positionIndex = 1;
  253.         for (i = 0; i < 3; i++) {
  254.             indices[indicesIndex++] = positionIndex++;
  255.             indices[indicesIndex++] = prevIndex;
  256.             indices[indicesIndex++] = positionIndex;
  257.         }
  259.         for (i = 2; i < numPts + 1; ++i) {
  260.             positionIndex = i * (i + 1) - 1;
  261.             prevIndex = (i - 1) * i - 1;
  263.             indices[indicesIndex++] = positionIndex++;
  264.             indices[indicesIndex++] = prevIndex;
  265.             indices[indicesIndex++] = positionIndex;
  267.             numInterior = 2 * i;
  268.             for (j = 0; j < numInterior - 1; ++j) {
  270.                 indices[indicesIndex++] = positionIndex;
  271.                 indices[indicesIndex++] = prevIndex++;
  272.                 indices[indicesIndex++] = prevIndex;
  274.                 indices[indicesIndex++] = positionIndex++;
  275.                 indices[indicesIndex++] = prevIndex;
  276.                 indices[indicesIndex++] = positionIndex;
  277.             }
  279.             indices[indicesIndex++] = positionIndex++;
  280.             indices[indicesIndex++] = prevIndex;
  281.             indices[indicesIndex++] = positionIndex;
  282.         }
  284.         // Indices for center column of triangles
  285.         numInterior = numPts * 2;
  286.         ++positionIndex;
  287.         ++prevIndex;
  288.         for (i = 0; i < numInterior - 1; ++i) {
  289.             indices[indicesIndex++] = positionIndex;
  290.             indices[indicesIndex++] = prevIndex++;
  291.             indices[indicesIndex++] = prevIndex;
  293.             indices[indicesIndex++] = positionIndex++;
  294.             indices[indicesIndex++] = prevIndex;
  295.             indices[indicesIndex++] = positionIndex;
  296.         }
  298.         indices[indicesIndex++] = positionIndex;
  299.         indices[indicesIndex++] = prevIndex++;
  300.         indices[indicesIndex++] = prevIndex;
  302.         indices[indicesIndex++] = positionIndex++;
  303.         indices[indicesIndex++] = prevIndex++;
  304.         indices[indicesIndex++] = prevIndex;
  306.         // Reverse the process creating indices to the 'left' of the north vector
  307.         ++prevIndex;
  308.         for (i = numPts - 1; i > 1; --i) {
  309.             indices[indicesIndex++] = prevIndex++;
  310.             indices[indicesIndex++] = prevIndex;
  311.             indices[indicesIndex++] = positionIndex;
  313.             numInterior = 2 * i;
  314.             for (j = 0; j < numInterior - 1; ++j) {
  315.                 indices[indicesIndex++] = positionIndex;
  316.                 indices[indicesIndex++] = prevIndex++;
  317.                 indices[indicesIndex++] = prevIndex;
  319.                 indices[indicesIndex++] = positionIndex++;
  320.                 indices[indicesIndex++] = prevIndex;
  321.                 indices[indicesIndex++] = positionIndex;
  322.             }
  324.             indices[indicesIndex++] = prevIndex++;
  325.             indices[indicesIndex++] = prevIndex++;
  326.             indices[indicesIndex++] = positionIndex++;
  327.         }
  329.         for (i = 0; i < 3; i++) {
  330.             indices[indicesIndex++] = prevIndex++;
  331.             indices[indicesIndex++] = prevIndex;
  332.             indices[indicesIndex++] = positionIndex;
  333.         }
  334.         return indices;
  335.     }
  337.     var boundingSphereCenter = new Cartesian3();
  339.     function computeEllipse(options) {
  340.         var center = options.center;
  341.         boundingSphereCenter = Cartesian3.multiplyByScalar(options.ellipsoid.geodeticSurfaceNormal(center, boundingSphereCenter), options.height, boundingSphereCenter);
  342.         boundingSphereCenter = Cartesian3.add(center, boundingSphereCenter, boundingSphereCenter);
  343.         var boundingSphere = new BoundingSphere(boundingSphereCenter, options.semiMajorAxis);
  344.         var cep = EllipseGeometryLibrary.computeEllipsePositions(options, true, false);
  345.         var positions = cep.positions;
  346.         var numPts = cep.numPts;
  347.         var attributes = computeTopBottomAttributes(positions, options, false);
  348.         var indices = topIndices(numPts);
  349.         indices = IndexDatatype.createTypedArray(positions.length / 3, indices);
  350.         return {
  351.             boundingSphere : boundingSphere,
  352.             attributes : attributes,
  353.             indices : indices
  354.         };
  355.     }
  357.     function computeWallAttributes(positions, options) {
  358.         var vertexFormat = options.vertexFormat;
  359.         var center = options.center;
  360.         var semiMajorAxis = options.semiMajorAxis;
  361.         var semiMinorAxis = options.semiMinorAxis;
  362.         var ellipsoid = options.ellipsoid;
  363.         var height = options.height;
  364.         var extrudedHeight = options.extrudedHeight;
  365.         var stRotation = options.stRotation;
  366.         var size = positions.length / 3 * 2;
  368.         var finalPositions = new Float64Array(size * 3);
  369.         var textureCoordinates = (vertexFormat.st) ? new Float32Array(size * 2) : undefined;
  370.         var normals = (vertexFormat.normal) ? new Float32Array(size * 3) : undefined;
  371.         var tangents = (vertexFormat.tangent) ? new Float32Array(size * 3) : undefined;
  372.         var binormals = (vertexFormat.binormal) ? new Float32Array(size * 3) : undefined;
  374.         var textureCoordIndex = 0;
  376.         // Raise positions to a height above the ellipsoid and compute the
  377.         // texture coordinates, normals, tangents, and binormals.
  378.         var normal = scratchNormal;
  379.         var tangent = scratchTangent;
  380.         var binormal = scratchBinormal;
  382.         var projection = new GeographicProjection(ellipsoid);
  383.         var projectedCenter = projection.project(ellipsoid.cartesianToCartographic(center, scratchCartographic), projectedCenterScratch);
  385.         var geodeticNormal = ellipsoid.scaleToGeodeticSurface(center, scratchCartesian1);
  386.         ellipsoid.geodeticSurfaceNormal(geodeticNormal, geodeticNormal);
  387.         var rotation = Quaternion.fromAxisAngle(geodeticNormal, stRotation, quaternionScratch);
  388.         var textureMatrix = Matrix3.fromQuaternion(rotation, textureMatrixScratch);
  390.         var minTexCoord = Cartesian2.fromElements(Number.POSITIVE_INFINITY, Number.POSITIVE_INFINITY, scratchMinTexCoord);
  391.         var maxTexCoord = Cartesian2.fromElements(Number.NEGATIVE_INFINITY, Number.NEGATIVE_INFINITY, scratchMaxTexCoord);
  393.         var length = positions.length;
  394.         var stOffset = length / 3 * 2;
  395.         for (var i = 0; i < length; i += 3) {
  396.             var i1 = i + 1;
  397.             var i2 = i + 2;
  398.             var position = Cartesian3.fromArray(positions, i, scratchCartesian1);
  399.             var extrudedPosition;
  401.             if (vertexFormat.st) {
  402.                 var rotatedPoint = Matrix3.multiplyByVector(textureMatrix, position, scratchCartesian2);
  403.                 var projectedPoint = projection.project(ellipsoid.cartesianToCartographic(rotatedPoint, scratchCartographic), scratchCartesian3);
  404.                 Cartesian3.subtract(projectedPoint, projectedCenter, projectedPoint);
  406.                 texCoordScratch.x = (projectedPoint.x + semiMajorAxis) / (2.0 * semiMajorAxis);
  407.                 texCoordScratch.y = (projectedPoint.y + semiMinorAxis) / (2.0 * semiMinorAxis);
  409.                 minTexCoord.x = Math.min(texCoordScratch.x, minTexCoord.x);
  410.                 minTexCoord.y = Math.min(texCoordScratch.y, minTexCoord.y);
  411.                 maxTexCoord.x = Math.max(texCoordScratch.x, maxTexCoord.x);
  412.                 maxTexCoord.y = Math.max(texCoordScratch.y, maxTexCoord.y);
  414.                 textureCoordinates[textureCoordIndex + stOffset] = texCoordScratch.x;
  415.                 textureCoordinates[textureCoordIndex + 1 + stOffset] = texCoordScratch.y;
  417.                 textureCoordinates[textureCoordIndex++] = texCoordScratch.x;
  418.                 textureCoordinates[textureCoordIndex++] = texCoordScratch.y;
  419.             }
  421.             position = ellipsoid.scaleToGeodeticSurface(position, position);
  422.             extrudedPosition = Cartesian3.clone(position, scratchCartesian2);
  423.             normal = ellipsoid.geodeticSurfaceNormal(position, normal);
  424.             var scaledNormal = Cartesian3.multiplyByScalar(normal, height, scratchCartesian4);
  425.             position = Cartesian3.add(position, scaledNormal, position);
  426.             scaledNormal = Cartesian3.multiplyByScalar(normal, extrudedHeight, scaledNormal);
  427.             extrudedPosition = Cartesian3.add(extrudedPosition, scaledNormal, extrudedPosition);
  429.             if (vertexFormat.position) {
  430.                 finalPositions[i + length] = extrudedPosition.x;
  431.                 finalPositions[i1 + length] = extrudedPosition.y;
  432.                 finalPositions[i2 + length] = extrudedPosition.z;
  434.                 finalPositions[i] = position.x;
  435.                 finalPositions[i1] = position.y;
  436.                 finalPositions[i2] = position.z;
  437.             }
  439.             if (vertexFormat.normal || vertexFormat.tangent || vertexFormat.binormal) {
  441.                 binormal = Cartesian3.clone(normal, binormal);
  442.                 var next = Cartesian3.fromArray(positions, (i + 3) % length, scratchCartesian4);
  443.                 Cartesian3.subtract(next, position, next);
  444.                 var bottom = Cartesian3.subtract(extrudedPosition, position, scratchCartesian3);
  446.                 normal = Cartesian3.normalize(Cartesian3.cross(bottom, next, normal), normal);
  448.                 if (vertexFormat.normal) {
  449.                     normals[i] = normal.x;
  450.                     normals[i1] = normal.y;
  451.                     normals[i2] = normal.z;
  453.                     normals[i + length] = normal.x;
  454.                     normals[i1 + length] = normal.y;
  455.                     normals[i2 + length] = normal.z;
  456.                 }
  458.                 if (vertexFormat.tangent) {
  459.                     tangent = Cartesian3.normalize(Cartesian3.cross(binormal, normal, tangent), tangent);
  460.                     tangents[i] = tangent.x;
  461.                     tangents[i1] = tangent.y;
  462.                     tangents[i2] = tangent.z;
  464.                     tangents[i + length] = tangent.x;
  465.                     tangents[i + 1 + length] = tangent.y;
  466.                     tangents[i + 2 + length] = tangent.z;
  467.                 }
  469.                 if (vertexFormat.binormal) {
  470.                     binormals[i] = binormal.x;
  471.                     binormals[i1] = binormal.y;
  472.                     binormals[i2] = binormal.z;
  474.                     binormals[i + length] = binormal.x;
  475.                     binormals[i1 + length] = binormal.y;
  476.                     binormals[i2 + length] = binormal.z;
  477.                 }
  478.             }
  479.         }
  481.         if (vertexFormat.st) {
  482.             length = textureCoordinates.length;
  483.             for (var k = 0; k < length; k += 2) {
  484.                 textureCoordinates[k] = (textureCoordinates[k] - minTexCoord.x) / (maxTexCoord.x - minTexCoord.x);
  485.                 textureCoordinates[k + 1] = (textureCoordinates[k + 1] - minTexCoord.y) / (maxTexCoord.y - minTexCoord.y);
  486.             }
  487.         }
  489.         var attributes = new GeometryAttributes();
  491.         if (vertexFormat.position) {
  492.             attributes.position = new GeometryAttribute({
  493.                 componentDatatype : ComponentDatatype.DOUBLE,
  494.                 componentsPerAttribute : 3,
  495.                 values : finalPositions
  496.             });
  497.         }
  499.         if (vertexFormat.st) {
  500.             attributes.st = new GeometryAttribute({
  501.                 componentDatatype : ComponentDatatype.FLOAT,
  502.                 componentsPerAttribute : 2,
  503.                 values : textureCoordinates
  504.             });
  505.         }
  507.         if (vertexFormat.normal) {
  508.             attributes.normal = new GeometryAttribute({
  509.                 componentDatatype : ComponentDatatype.FLOAT,
  510.                 componentsPerAttribute : 3,
  511.                 values : normals
  512.             });
  513.         }
  515.         if (vertexFormat.tangent) {
  516.             attributes.tangent = new GeometryAttribute({
  517.                 componentDatatype : ComponentDatatype.FLOAT,
  518.                 componentsPerAttribute : 3,
  519.                 values : tangents
  520.             });
  521.         }
  523.         if (vertexFormat.binormal) {
  524.             attributes.binormal = new GeometryAttribute({
  525.                 componentDatatype : ComponentDatatype.FLOAT,
  526.                 componentsPerAttribute : 3,
  527.                 values : binormals
  528.             });
  529.         }
  530.         return attributes;
  531.     }
  533.     function computeWallIndices(positions) {
  534.         var length = positions.length / 3;
  535.         var indices = IndexDatatype.createTypedArray(length, length * 6);
  536.         var index = 0;
  537.         for (var i = 0; i < length; i++) {
  538.             var UL = i;
  539.             var LL = i + length;
  540.             var UR = (UL + 1) % length;
  541.             var LR = UR + length;
  542.             indices[index++] = UL;
  543.             indices[index++] = LL;
  544.             indices[index++] = UR;
  545.             indices[index++] = UR;
  546.             indices[index++] = LL;
  547.             indices[index++] = LR;
  548.         }
  550.         return indices;
  551.     }
  553.     var topBoundingSphere = new BoundingSphere();
  554.     var bottomBoundingSphere = new BoundingSphere();
  556.     function computeExtrudedEllipse(options) {
  557.         var center = options.center;
  558.         var ellipsoid = options.ellipsoid;
  559.         var semiMajorAxis = options.semiMajorAxis;
  560.         var scaledNormal = Cartesian3.multiplyByScalar(ellipsoid.geodeticSurfaceNormal(center, scratchCartesian1), options.height, scratchCartesian1);
  561.         topBoundingSphere.center = Cartesian3.add(center, scaledNormal, topBoundingSphere.center);
  562.         topBoundingSphere.radius = semiMajorAxis;
  564.         scaledNormal = Cartesian3.multiplyByScalar(ellipsoid.geodeticSurfaceNormal(center, scaledNormal), options.extrudedHeight, scaledNormal);
  565.         bottomBoundingSphere.center = Cartesian3.add(center, scaledNormal, bottomBoundingSphere.center);
  566.         bottomBoundingSphere.radius = semiMajorAxis;
  568.         var cep = EllipseGeometryLibrary.computeEllipsePositions(options, true, true);
  569.         var positions = cep.positions;
  570.         var numPts = cep.numPts;
  571.         var outerPositions = cep.outerPositions;
  572.         var boundingSphere = BoundingSphere.union(topBoundingSphere, bottomBoundingSphere);
  573.         var topBottomAttributes = computeTopBottomAttributes(positions, options, true);
  574.         var indices = topIndices(numPts);
  575.         var length = indices.length;
  576.         indices.length = length * 2;
  577.         var posLength = positions.length / 3;
  578.         for (var i = 0; i < length; i += 3) {
  579.             indices[i + length] = indices[i + 2] + posLength;
  580.             indices[i + 1 + length] = indices[i + 1] + posLength;
  581.             indices[i + 2 + length] = indices[i] + posLength;
  582.         }
  584.         var topBottomIndices = IndexDatatype.createTypedArray(posLength * 2 / 3, indices);
  586.         var topBottomGeo = new Geometry({
  587.             attributes : topBottomAttributes,
  588.             indices : topBottomIndices,
  589.             primitiveType : PrimitiveType.TRIANGLES
  590.         });
  592.         var wallAttributes = computeWallAttributes(outerPositions, options);
  593.         indices = computeWallIndices(outerPositions);
  594.         var wallIndices = IndexDatatype.createTypedArray(outerPositions.length * 2 / 3, indices);
  596.         var wallGeo = new Geometry({
  597.             attributes : wallAttributes,
  598.             indices : wallIndices,
  599.             primitiveType : PrimitiveType.TRIANGLES
  600.         });
  602.         var geo = GeometryPipeline.combineInstances([
  603.             new GeometryInstance({
  604.                 geometry : topBottomGeo
  605.             }),
  606.             new GeometryInstance({
  607.                 geometry : wallGeo
  608.             })
  609.         ]);
  611.         return {
  612.             boundingSphere : boundingSphere,
  613.             attributes : geo[0].attributes,
  614.             indices : geo[0].indices
  615.         };
  616.     }
  618.     var scratchEnuToFixedMatrix = new Matrix4();
  619.     var scratchFixedToEnuMatrix = new Matrix4();
  620.     var scratchRotationMatrix = new Matrix3();
  621.     var scratchRectanglePoints = [new Cartesian3(), new Cartesian3(), new Cartesian3(), new Cartesian3()];
  622.     var scratchCartographicPoints = [new Cartographic(), new Cartographic(), new Cartographic(), new Cartographic()];
  624.     function computeRectangle(center, ellipsoid, semiMajorAxis, semiMinorAxis, rotation) {
  625.         Transforms.eastNorthUpToFixedFrame(center, ellipsoid, scratchEnuToFixedMatrix);
  626.         Matrix4.inverseTransformation(scratchEnuToFixedMatrix, scratchFixedToEnuMatrix);
  628.         // Find the 4 extreme points of the ellipse in ENU
  629.         for (var i = 0; i < 4; ++i) {
  630.             Cartesian3.clone(Cartesian3.ZERO, scratchRectanglePoints[i]);
  631.         }
  632.         scratchRectanglePoints[0].x += semiMajorAxis;
  633.         scratchRectanglePoints[1].x -= semiMajorAxis;
  634.         scratchRectanglePoints[2].y += semiMinorAxis;
  635.         scratchRectanglePoints[3].y -= semiMinorAxis;
  637.         Matrix3.fromRotationZ(rotation, scratchRotationMatrix);
  638.         for (i = 0; i < 4; ++i) {
  639.             // Apply the rotation
  640.             Matrix3.multiplyByVector(scratchRotationMatrix, scratchRectanglePoints[i], scratchRectanglePoints[i]);
  642.             // Convert back to fixed and then to cartographic
  643.             Matrix4.multiplyByPoint(scratchEnuToFixedMatrix, scratchRectanglePoints[i], scratchRectanglePoints[i]);
  644.             ellipsoid.cartesianToCartographic(scratchRectanglePoints[i], scratchCartographicPoints[i]);
  645.         }
  647.         return Rectangle.fromCartographicArray(scratchCartographicPoints);
  648.     }
  650.     /**
  651.      * A description of an ellipse on an ellipsoid. Ellipse geometry can be rendered with both {@link Primitive} and {@link GroundPrimitive}.
  652.      *
  653.      * @alias EllipseGeometry
  654.      * @constructor
  655.      *
  656.      * @param {Object} options Object with the following properties:
  657.      * @param {Cartesian3} options.center The ellipse's center point in the fixed frame.
  658.      * @param {Number} options.semiMajorAxis The length of the ellipse's semi-major axis in meters.
  659.      * @param {Number} options.semiMinorAxis The length of the ellipse's semi-minor axis in meters.
  660.      * @param {Ellipsoid} [options.ellipsoid=Ellipsoid.WGS84] The ellipsoid the ellipse will be on.
  661.      * @param {Number} [options.height=0.0] The distance in meters between the ellipse and the ellipsoid surface.
  662.      * @param {Number} [options.extrudedHeight] The distance in meters between the ellipse's extruded face and the ellipsoid surface.
  663.      * @param {Number} [options.rotation=0.0] The angle of rotation counter-clockwise from north.
  664.      * @param {Number} [options.stRotation=0.0] The rotation of the texture coordinates counter-clockwise from north.
  665.      * @param {Number} [options.granularity=CesiumMath.RADIANS_PER_DEGREE] The angular distance between points on the ellipse in radians.
  666.      * @param {VertexFormat} [options.vertexFormat=VertexFormat.DEFAULT] The vertex attributes to be computed.
  667.      *
  668.      * @exception {DeveloperError} semiMajorAxis and semiMinorAxis must be greater than zero.
  669.      * @exception {DeveloperError} semiMajorAxis must be greater than or equal to the semiMinorAxis.
  670.      * @exception {DeveloperError} granularity must be greater than zero.
  671.      *
  672.      *
  673.      * @example
  674.      * // Create an ellipse.
  675.      * var ellipse = new Cesium.EllipseGeometry({
  676.      *   center : Cesium.Cartesian3.fromDegrees(-75.59777, 40.03883),
  677.      *   semiMajorAxis : 500000.0,
  678.      *   semiMinorAxis : 300000.0,
  679.      *   rotation : Cesium.Math.toRadians(60.0)
  680.      * });
  681.      * var geometry = Cesium.EllipseGeometry.createGeometry(ellipse);
  682.      *
  683.      * @see EllipseGeometry.createGeometry
  684.      */
  685.     function EllipseGeometry(options) {
  686.         options = defaultValue(options, defaultValue.EMPTY_OBJECT);
  688.         var center = options.center;
  689.         var ellipsoid = defaultValue(options.ellipsoid, Ellipsoid.WGS84);
  690.         var semiMajorAxis = options.semiMajorAxis;
  691.         var semiMinorAxis = options.semiMinorAxis;
  692.         var granularity = defaultValue(options.granularity, CesiumMath.RADIANS_PER_DEGREE);
  693.         var height = defaultValue(options.height, 0.0);
  694.         var extrudedHeight = options.extrudedHeight;
  695.         var extrude = (defined(extrudedHeight) && Math.abs(height - extrudedHeight) > 1.0);
  696.         var vertexFormat = defaultValue(options.vertexFormat, VertexFormat.DEFAULT);
  698.         //>>includeStart('debug', pragmas.debug);
  699.         if (!defined(center)) {
  700.             throw new DeveloperError('center is required.');
  701.         }
  702.         if (!defined(semiMajorAxis)) {
  703.             throw new DeveloperError('semiMajorAxis is required.');
  704.         }
  705.         if (!defined(semiMinorAxis)) {
  706.             throw new DeveloperError('semiMinorAxis is required.');
  707.         }
  708.         if (semiMajorAxis < semiMinorAxis) {
  709.             throw new DeveloperError('semiMajorAxis must be greater than or equal to the semiMinorAxis.');
  710.         }
  711.         if (granularity <= 0.0) {
  712.             throw new DeveloperError('granularity must be greater than zero.');
  713.         }
  714.         //>>includeEnd('debug');
  716.         this._center = Cartesian3.clone(center);
  717.         this._semiMajorAxis = semiMajorAxis;
  718.         this._semiMinorAxis = semiMinorAxis;
  719.         this._ellipsoid = Ellipsoid.clone(ellipsoid);
  720.         this._rotation = defaultValue(options.rotation, 0.0);
  721.         this._stRotation = defaultValue(options.stRotation, 0.0);
  722.         this._height = height;
  723.         this._granularity = granularity;
  724.         this._vertexFormat = VertexFormat.clone(vertexFormat);
  725.         this._extrudedHeight = defaultValue(extrudedHeight, height);
  726.         this._extrude = extrude;
  727.         this._workerName = 'createEllipseGeometry';
  729.         this._rectangle = computeRectangle(this._center, this._ellipsoid, semiMajorAxis, semiMinorAxis, this._rotation);
  730.     }
  732.     /**
  733.      * The number of elements used to pack the object into an array.
  734.      * @type {Number}
  735.      */
  736.     EllipseGeometry.packedLength = Cartesian3.packedLength + Ellipsoid.packedLength + VertexFormat.packedLength + Rectangle.packedLength + 8;
  738.     /**
  739.      * Stores the provided instance into the provided array.
  740.      *
  741.      * @param {EllipseGeometry} value The value to pack.
  742.      * @param {Number[]} array The array to pack into.
  743.      * @param {Number} [startingIndex=0] The index into the array at which to start packing the elements.
  744.      *
  745.      * @returns {Number[]} The array that was packed into
  746.      */
  747.     EllipseGeometry.pack = function(value, array, startingIndex) {
  748.         //>>includeStart('debug', pragmas.debug);
  749.         if (!defined(value)) {
  750.             throw new DeveloperError('value is required');
  751.         }
  752.         if (!defined(array)) {
  753.             throw new DeveloperError('array is required');
  754.         }
  755.         //>>includeEnd('debug');
  757.         startingIndex = defaultValue(startingIndex, 0);
  759.         Cartesian3.pack(value._center, array, startingIndex);
  760.         startingIndex += Cartesian3.packedLength;
  762.         Ellipsoid.pack(value._ellipsoid, array, startingIndex);
  763.         startingIndex += Ellipsoid.packedLength;
  765.         VertexFormat.pack(value._vertexFormat, array, startingIndex);
  766.         startingIndex += VertexFormat.packedLength;
  768.         Rectangle.pack(value._rectangle, array, startingIndex);
  769.         startingIndex += Rectangle.packedLength;
  771.         array[startingIndex++] = value._semiMajorAxis;
  772.         array[startingIndex++] = value._semiMinorAxis;
  773.         array[startingIndex++] = value._rotation;
  774.         array[startingIndex++] = value._stRotation;
  775.         array[startingIndex++] = value._height;
  776.         array[startingIndex++] = value._granularity;
  777.         array[startingIndex++] = value._extrudedHeight;
  778.         array[startingIndex] = value._extrude ? 1.0 : 0.0;
  780.         return array;
  781.     };
  783.     var scratchCenter = new Cartesian3();
  784.     var scratchEllipsoid = new Ellipsoid();
  785.     var scratchVertexFormat = new VertexFormat();
  786.     var scratchRectangle = new Rectangle();
  787.     var scratchOptions = {
  788.         center : scratchCenter,
  789.         ellipsoid : scratchEllipsoid,
  790.         vertexFormat : scratchVertexFormat,
  791.         semiMajorAxis : undefined,
  792.         semiMinorAxis : undefined,
  793.         rotation : undefined,
  794.         stRotation : undefined,
  795.         height : undefined,
  796.         granularity : undefined,
  797.         extrudedHeight : undefined
  798.     };
  800.     /**
  801.      * Retrieves an instance from a packed array.
  802.      *
  803.      * @param {Number[]} array The packed array.
  804.      * @param {Number} [startingIndex=0] The starting index of the element to be unpacked.
  805.      * @param {EllipseGeometry} [result] The object into which to store the result.
  806.      * @returns {EllipseGeometry} The modified result parameter or a new EllipseGeometry instance if one was not provided.
  807.      */
  808.     EllipseGeometry.unpack = function(array, startingIndex, result) {
  809.         //>>includeStart('debug', pragmas.debug);
  810.         if (!defined(array)) {
  811.             throw new DeveloperError('array is required');
  812.         }
  813.         //>>includeEnd('debug');
  815.         startingIndex = defaultValue(startingIndex, 0);
  817.         var center = Cartesian3.unpack(array, startingIndex, scratchCenter);
  818.         startingIndex += Cartesian3.packedLength;
  820.         var ellipsoid = Ellipsoid.unpack(array, startingIndex, scratchEllipsoid);
  821.         startingIndex += Ellipsoid.packedLength;
  823.         var vertexFormat = VertexFormat.unpack(array, startingIndex, scratchVertexFormat);
  824.         startingIndex += VertexFormat.packedLength;
  826.         var rectangle = Rectangle.unpack(array, startingIndex, scratchRectangle);
  827.         startingIndex += Rectangle.packedLength;
  829.         var semiMajorAxis = array[startingIndex++];
  830.         var semiMinorAxis = array[startingIndex++];
  831.         var rotation = array[startingIndex++];
  832.         var stRotation = array[startingIndex++];
  833.         var height = array[startingIndex++];
  834.         var granularity = array[startingIndex++];
  835.         var extrudedHeight = array[startingIndex++];
  836.         var extrude = array[startingIndex] === 1.0;
  838.         if (!defined(result)) {
  839.             scratchOptions.height = height;
  840.             scratchOptions.extrudedHeight = extrudedHeight;
  841.             scratchOptions.granularity = granularity;
  842.             scratchOptions.stRotation = stRotation;
  843.             scratchOptions.rotation = rotation;
  844.             scratchOptions.semiMajorAxis = semiMajorAxis;
  845.             scratchOptions.semiMinorAxis = semiMinorAxis;
  846.             return new EllipseGeometry(scratchOptions);
  847.         }
  849.         result._center = Cartesian3.clone(center, result._center);
  850.         result._ellipsoid = Ellipsoid.clone(ellipsoid, result._ellipsoid);
  851.         result._vertexFormat = VertexFormat.clone(vertexFormat, result._vertexFormat);
  852.         result._semiMajorAxis = semiMajorAxis;
  853.         result._semiMinorAxis = semiMinorAxis;
  854.         result._rotation = rotation;
  855.         result._stRotation = stRotation;
  856.         result._height = height;
  857.         result._granularity = granularity;
  858.         result._extrudedHeight = extrudedHeight;
  859.         result._extrude = extrude;
  860.         result._rectangle = Rectangle.clone(rectangle);
  862.         return result;
  863.     };
  865.     /**
  866.      * Computes the geometric representation of a ellipse on an ellipsoid, including its vertices, indices, and a bounding sphere.
  867.      *
  868.      * @param {EllipseGeometry} ellipseGeometry A description of the ellipse.
  869.      * @returns {Geometry|undefined} The computed vertices and indices.
  870.      */
  871.     EllipseGeometry.createGeometry = function(ellipseGeometry) {
  872.         if ((ellipseGeometry._semiMajorAxis <= 0.0) || (ellipseGeometry._semiMinorAxis <= 0.0)) {
  873.             return;
  874.         }
  876.         ellipseGeometry._center = ellipseGeometry._ellipsoid.scaleToGeodeticSurface(ellipseGeometry._center, ellipseGeometry._center);
  877.         var options = {
  878.             center : ellipseGeometry._center,
  879.             semiMajorAxis : ellipseGeometry._semiMajorAxis,
  880.             semiMinorAxis : ellipseGeometry._semiMinorAxis,
  881.             ellipsoid : ellipseGeometry._ellipsoid,
  882.             rotation : ellipseGeometry._rotation,
  883.             height : ellipseGeometry._height,
  884.             extrudedHeight : ellipseGeometry._extrudedHeight,
  885.             granularity : ellipseGeometry._granularity,
  886.             vertexFormat : ellipseGeometry._vertexFormat,
  887.             stRotation : ellipseGeometry._stRotation
  888.         };
  889.         var geometry;
  890.         if (ellipseGeometry._extrude) {
  891.             options.extrudedHeight = Math.min(ellipseGeometry._extrudedHeight, ellipseGeometry._height);
  892.             options.height = Math.max(ellipseGeometry._extrudedHeight, ellipseGeometry._height);
  893.             geometry = computeExtrudedEllipse(options);
  894.         } else {
  895.             geometry = computeEllipse(options);
  896.         }
  898.         return new Geometry({
  899.             attributes : geometry.attributes,
  900.             indices : geometry.indices,
  901.             primitiveType : PrimitiveType.TRIANGLES,
  902.             boundingSphere : geometry.boundingSphere
  903.         });
  904.     };
  906.     /**
  907.      * @private
  908.      */
  909.     EllipseGeometry.createShadowVolume = function(ellipseGeometry, minHeightFunc, maxHeightFunc) {
  910.         var granularity = ellipseGeometry._granularity;
  911.         var ellipsoid = ellipseGeometry._ellipsoid;
  913.         var minHeight = minHeightFunc(granularity, ellipsoid);
  914.         var maxHeight = maxHeightFunc(granularity, ellipsoid);
  916.         return new EllipseGeometry({
  917.             center : ellipseGeometry._center,
  918.             semiMajorAxis : ellipseGeometry._semiMajorAxis,
  919.             semiMinorAxis : ellipseGeometry._semiMinorAxis,
  920.             ellipsoid : ellipsoid,
  921.             rotation : ellipseGeometry._rotation,
  922.             stRotation : ellipseGeometry._stRotation,
  923.             granularity : granularity,
  924.             extrudedHeight : minHeight,
  925.             height : maxHeight,
  926.             vertexFormat : VertexFormat.POSITION_ONLY
  927.         });
  928.     };
  930.     defineProperties(EllipseGeometry.prototype, {
  931.         /**
  932.          * @private
  933.          */
  934.         rectangle : {
  935.             get : function() {
  936.                 return this._rectangle;
  937.             }
  938.         }
  939.     });
  941.     return EllipseGeometry;
  942. });