/*global define*/
define(['Core/defaultValue', 'Core/defined', 'Core/BoundingSphere', 'Core/Cartesian3', 'Core/ComponentDatatype', 'Core/DeveloperError', 'Core/Ellipsoid', 'Core/EllipsoidTangentPlane', 'Core/Geometry', 'Core/GeometryAttribute', 'Core/GeometryAttributes', 'Core/GeometryInstance', 'Core/GeometryPipeline', 'Core/IndexDatatype', 'Core/Math', 'Core/PolygonGeometryLibrary', 'Core/PolygonPipeline', 'Core/PrimitiveType', 'Core/Queue', 'Core/WindingOrder'], function(
defaultValue,
defined,
BoundingSphere,
Cartesian3,
ComponentDatatype,
DeveloperError,
Ellipsoid,
EllipsoidTangentPlane,
Geometry,
GeometryAttribute,
GeometryAttributes,
GeometryInstance,
GeometryPipeline,
IndexDatatype,
CesiumMath,
PolygonGeometryLibrary,
PolygonPipeline,
PrimitiveType,
Queue,
WindingOrder) {
"use strict";
var createGeometryFromPositionsPositions = [];
function createGeometryFromPositions(ellipsoid, positions, granularity, perPositionHeight) {
var cleanedPositions = PolygonPipeline.removeDuplicates(positions);
//>>includeStart('debug', pragmas.debug);
if (cleanedPositions.length < 3) {
throw new DeveloperError('Duplicate positions result in not enough positions to form a polygon.');
}
//>>includeEnd('debug');
var tangentPlane = EllipsoidTangentPlane.fromPoints(cleanedPositions, ellipsoid);
var positions2D = tangentPlane.projectPointsOntoPlane(cleanedPositions, createGeometryFromPositionsPositions);
var originalWindingOrder = PolygonPipeline.computeWindingOrder2D(positions2D);
if (originalWindingOrder === WindingOrder.CLOCKWISE) {
positions2D.reverse();
cleanedPositions.reverse();
}
var subdividedPositions = [];
var length = cleanedPositions.length;
var i;
if (!perPositionHeight) {
for (i = 0; i < length; i++) {
subdividedPositions = subdividedPositions.concat(PolygonGeometryLibrary.subdivideLine(cleanedPositions[i], cleanedPositions[(i + 1) % length], granularity));
}
} else {
for (i = 0; i < length; i++) {
var p0 = cleanedPositions[i];
var p1 = cleanedPositions[(i + 1) % length];
subdividedPositions.push(p0.x, p0.y, p0.z, p1.x, p1.y, p1.z);
}
}
length = subdividedPositions.length / 3;
var indicesSize = length * 2;
var indices = IndexDatatype.createTypedArray(subdividedPositions.length / 3, indicesSize);
var index = 0;
for (i = 0; i < length - 1; i++) {
indices[index++] = i;
indices[index++] = i + 1;
}
indices[index++] = length - 1;
indices[index++] = 0;
return new GeometryInstance({
geometry : new Geometry({
attributes : new GeometryAttributes({
position : new GeometryAttribute({
componentDatatype : ComponentDatatype.DOUBLE,
componentsPerAttribute : 3,
values : new Float64Array(subdividedPositions)
})
}),
indices : indices,
primitiveType : PrimitiveType.LINES
})
});
}
var scratchPosition = new Cartesian3();
var scratchNormal = new Cartesian3();
var scratchBoundingSphere = new BoundingSphere();
function createGeometryFromPositionsExtruded(ellipsoid, positions, granularity, perPositionHeight) {
var cleanedPositions = PolygonPipeline.removeDuplicates(positions);
//>>includeStart('debug', pragmas.debug);
if (cleanedPositions.length < 3) {
throw new DeveloperError('Duplicate positions result in not enough positions to form a polygon.');
}
//>>includeEnd('debug');
var tangentPlane = EllipsoidTangentPlane.fromPoints(cleanedPositions, ellipsoid);
var positions2D = tangentPlane.projectPointsOntoPlane(cleanedPositions, createGeometryFromPositionsPositions);
var originalWindingOrder = PolygonPipeline.computeWindingOrder2D(positions2D);
if (originalWindingOrder === WindingOrder.CLOCKWISE) {
positions2D.reverse();
cleanedPositions.reverse();
}
var subdividedPositions = [];
var length = cleanedPositions.length;
var i;
var corners = new Array(length);
corners[0] = 0;
if (!perPositionHeight) {
for (i = 0; i < length - 1; i++) {
subdividedPositions = subdividedPositions.concat(PolygonGeometryLibrary.subdivideLine(cleanedPositions[i], cleanedPositions[i + 1], granularity));
corners[i + 1] = subdividedPositions.length / 3;
}
subdividedPositions = subdividedPositions.concat(PolygonGeometryLibrary.subdivideLine(cleanedPositions[length - 1], cleanedPositions[0], granularity));
} else {
var p0;
var p1;
for (i = 0; i < length - 1; i++) {
p0 = cleanedPositions[i];
p1 = cleanedPositions[(i + 1) % length];
subdividedPositions.push(p0.x, p0.y, p0.z, p1.x, p1.y, p1.z);
corners[i + 1] = subdividedPositions.length / 3;
}
p0 = cleanedPositions[length - 1];
p1 = cleanedPositions[0];
subdividedPositions.push(p0.x, p0.y, p0.z, p1.x, p1.y, p1.z);
}
length = subdividedPositions.length / 3;
var indicesSize = ((length * 2) + corners.length) * 2;
var indices = IndexDatatype.createTypedArray(subdividedPositions.length / 3, indicesSize);
var index = 0;
for (i = 0; i < length - 1; i++) {
indices[index++] = i;
indices[index++] = i + 1;
indices[index++] = i + length;
indices[index++] = i + 1 + length;
}
indices[index++] = length - 1;
indices[index++] = 0;
indices[index++] = length + length - 1;
indices[index++] = length;
for (i = 0; i < corners.length; i++) {
var corner = corners[i];
indices[index++] = corner;
indices[index++] = corner + length;
}
subdividedPositions = subdividedPositions.concat(subdividedPositions);
return new GeometryInstance({
geometry : new Geometry({
attributes : new GeometryAttributes({
position : new GeometryAttribute({
componentDatatype : ComponentDatatype.DOUBLE,
componentsPerAttribute : 3,
values : new Float64Array(subdividedPositions)
})
}),
indices : indices,
primitiveType : PrimitiveType.LINES
})
});
}
/**
* A description of the outline of a polygon on the ellipsoid. The polygon is defined by a polygon hierarchy.
*
* @alias PolygonOutlineGeometry
* @constructor
*
* @param {Object} options.polygonHierarchy A polygon hierarchy that can include holes.
* @param {Number} [options.height=0.0] The height of the polygon.
* @param {Number} [options.extrudedHeight] The height of the polygon.
* @param {VertexFormat} [options.vertexFormat=VertexFormat.DEFAULT] The vertex attributes to be computed.
* @param {Ellipsoid} [options.ellipsoid=Ellipsoid.WGS84] The ellipsoid to be used as a reference.
* @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 {Boolean} [options.perPositionHeight=false] Use the height of options.positions for each position instead of using options.height to determine the height.
*
* @exception {DeveloperError} polygonHierarchy is required.
*
* @see PolygonOutlineGeometry#createGeometry
* @see PolygonOutlineGeometry#fromPositions
*
* @example
* // 1. create a polygon outline from points
* var polygon = new Cesium.PolygonOutlineGeometry({
* polygonHierarchy : {
* positions : ellipsoid.cartographicArrayToCartesianArray([
* Cesium.Cartographic.fromDegrees(-72.0, 40.0),
* Cesium.Cartographic.fromDegrees(-70.0, 35.0),
* Cesium.Cartographic.fromDegrees(-75.0, 30.0),
* Cesium.Cartographic.fromDegrees(-70.0, 30.0),
* Cesium.Cartographic.fromDegrees(-68.0, 40.0)
* ])
* }
* });
* var geometry = Cesium.PolygonOutlineGeometry.createGeometry(polygon);
*
* // 2. create a nested polygon with holes outline
* var polygonWithHole = new Cesium.PolygonOutlineGeometry({
* polygonHierarchy : {
* positions : ellipsoid.cartographicArrayToCartesianArray([
* Cesium.Cartographic.fromDegrees(-109.0, 30.0),
* Cesium.Cartographic.fromDegrees(-95.0, 30.0),
* Cesium.Cartographic.fromDegrees(-95.0, 40.0),
* Cesium.Cartographic.fromDegrees(-109.0, 40.0)
* ]),
* holes : [{
* positions : ellipsoid.cartographicArrayToCartesianArray([
* Cesium.Cartographic.fromDegrees(-107.0, 31.0),
* Cesium.Cartographic.fromDegrees(-107.0, 39.0),
* Cesium.Cartographic.fromDegrees(-97.0, 39.0),
* Cesium.Cartographic.fromDegrees(-97.0, 31.0)
* ]),
* holes : [{
* positions : ellipsoid.cartographicArrayToCartesianArray([
* Cesium.Cartographic.fromDegrees(-105.0, 33.0),
* Cesium.Cartographic.fromDegrees(-99.0, 33.0),
* Cesium.Cartographic.fromDegrees(-99.0, 37.0),
* Cesium.Cartographic.fromDegrees(-105.0, 37.0)
* ]),
* holes : [{
* positions : ellipsoid.cartographicArrayToCartesianArray([
* Cesium.Cartographic.fromDegrees(-103.0, 34.0),
* Cesium.Cartographic.fromDegrees(-101.0, 34.0),
* Cesium.Cartographic.fromDegrees(-101.0, 36.0),
* Cesium.Cartographic.fromDegrees(-103.0, 36.0)
* ])
* }]
* }]
* }]
* }
* });
* var geometry = PolygonOutlineGeometry.createGeometry(polygonWithHole);
*
* // 3. create extruded polygon outline
* var extrudedPolygon = new PolygonOutlineGeometry({
* positions : ellipsoid.cartographicArrayToCartesianArray([
* Cesium.Cartographic.fromDegrees(-72.0, 40.0),
* Cesium.Cartographic.fromDegrees(-70.0, 35.0),
* Cesium.Cartographic.fromDegrees(-75.0, 30.0),
* Cesium.Cartographic.fromDegrees(-70.0, 30.0),
* Cesium.Cartographic.fromDegrees(-68.0, 40.0)
* ]),
* extrudedHeight: 300000
* });
* var geometry = PolygonOutlineGeometry.createGeometry(extrudedPolygon);
*
*/
var PolygonOutlineGeometry = function(options) {
options = defaultValue(options, defaultValue.EMPTY_OBJECT);
var ellipsoid = defaultValue(options.ellipsoid, Ellipsoid.WGS84);
var granularity = defaultValue(options.granularity, CesiumMath.RADIANS_PER_DEGREE);
var height = defaultValue(options.height, 0.0);
var perPositionHeight = defaultValue(options.perPositionHeight, false);
var extrudedHeight = options.extrudedHeight;
var extrude = (defined(extrudedHeight) && (!CesiumMath.equalsEpsilon(height, extrudedHeight, CesiumMath.EPSILON6) || perPositionHeight));
if (extrude) {
var h = extrudedHeight;
extrudedHeight = Math.min(h, height);
height = Math.max(h, height);
}
var polygonHierarchy = options.polygonHierarchy;
//>>includeStart('debug', pragmas.debug);
if (!defined(polygonHierarchy)) {
throw new DeveloperError('options.polygonHierarchy is required.');
}
//>>includeEnd('debug');
this._ellipsoid = ellipsoid;
this._granularity = granularity;
this._height = height;
this._extrudedHeight = extrudedHeight;
this._extrude = extrude;
this._polygonHierarchy = polygonHierarchy;
this._perPositionHeight = perPositionHeight;
this._workerName = 'createPolygonOutlineGeometry';
};
/**
* A description of a polygon outline from an array of positions.
*
* @memberof PolygonOutlineGeometry
*
* @param {Array} options.positions An array of positions that defined the corner points of the polygon.
* @param {Number} [options.height=0.0] The height of the polygon.
* @param {Number} [options.extrudedHeight] The height of the polygon extrusion.
* @param {Ellipsoid} [options.ellipsoid=Ellipsoid.WGS84] The ellipsoid to be used as a reference.
* @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 {Boolean} [options.perPositionHeight=false] Use the height of options.positions for each position instead of using options.height to determine the height.
*
* @exception {DeveloperError} options.positions is required.
*
* @see PolygonGeometry#createGeometry
*
* @example
* // create a polygon from points
* var polygon = Cesium.PolygonOutlineGeometry.fromPositions({
* positions : ellipsoid.cartographicArrayToCartesianArray([
* Cesium.Cartographic.fromDegrees(-72.0, 40.0),
* Cesium.Cartographic.fromDegrees(-70.0, 35.0),
* Cesium.Cartographic.fromDegrees(-75.0, 30.0),
* Cesium.Cartographic.fromDegrees(-70.0, 30.0),
* Cesium.Cartographic.fromDegrees(-68.0, 40.0)
* ])
* });
* var geometry = Cesium.PolygonOutlineGeometry.createGeometry(polygon);
*/
PolygonOutlineGeometry.fromPositions = function(options) {
options = defaultValue(options, defaultValue.EMPTY_OBJECT);
//>>includeStart('debug', pragmas.debug);
if (!defined(options.positions)) {
throw new DeveloperError('options.positions is required.');
}
//>>includeEnd('debug');
var newOptions = {
polygonHierarchy : {
positions : options.positions
},
height : options.height,
extrudedHeight : options.extrudedHeight,
ellipsoid : options.ellipsoid,
granularity : options.granularity,
perPositionHeight : options.perPositionHeight
};
return new PolygonOutlineGeometry(newOptions);
};
/**
* Computes the geometric representation of a polygon outline, including its vertices, indices, and a bounding sphere.
* @memberof PolygonOutlineGeometry
*
* @param {PolygonOutlineGeometry} polygonGeometry A description of the polygon outline.
* @returns {Geometry} The computed vertices and indices.
*
* @exception {DeveloperError} At least three positions are required.
* @exception {DeveloperError} Duplicate positions result in not enough positions to form a polygon.
*/
PolygonOutlineGeometry.createGeometry = function(polygonGeometry) {
var ellipsoid = polygonGeometry._ellipsoid;
var granularity = polygonGeometry._granularity;
var height = polygonGeometry._height;
var extrudedHeight = polygonGeometry._extrudedHeight;
var extrude = polygonGeometry._extrude;
var polygonHierarchy = polygonGeometry._polygonHierarchy;
var perPositionHeight = polygonGeometry._perPositionHeight;
var boundingSphere;
var outerPositions;
// create from a polygon hierarchy
// Algorithm adapted from http://www.geometrictools.com/Documentation/TriangulationByEarClipping.pdf
var polygons = [];
var queue = new Queue();
queue.enqueue(polygonHierarchy);
var i;
while (queue.length !== 0) {
var outerNode = queue.dequeue();
var outerRing = outerNode.positions;
if (outerRing.length < 3) {
throw new DeveloperError('At least three positions are required.');
}
var numChildren = outerNode.holes ? outerNode.holes.length : 0;
// The outer polygon contains inner polygons
for (i = 0; i < numChildren; i++) {
var hole = outerNode.holes[i];
polygons.push(hole.positions);
var numGrandchildren = 0;
if (defined(hole.holes)) {
numGrandchildren = hole.holes.length;
}
for ( var j = 0; j < numGrandchildren; j++) {
queue.enqueue(hole.holes[j]);
}
}
polygons.push(outerRing);
}
outerPositions = polygons[0];
// The bounding volume is just around the boundary points, so there could be cases for
// contrived polygons on contrived ellipsoids - very oblate ones - where the bounding
// volume doesn't cover the polygon.
boundingSphere = BoundingSphere.fromPoints(outerPositions);
var geometry;
var geometries = [];
if (extrude) {
for (i = 0; i < polygons.length; i++) {
geometry = createGeometryFromPositionsExtruded(ellipsoid, polygons[i], granularity, perPositionHeight);
if (defined(geometry)) {
geometry.geometry = PolygonGeometryLibrary.scaleToGeodeticHeightExtruded(geometry.geometry, height, extrudedHeight, ellipsoid, perPositionHeight);
geometries.push(geometry);
}
}
} else {
for (i = 0; i < polygons.length; i++) {
geometry = createGeometryFromPositions(ellipsoid, polygons[i], granularity, perPositionHeight);
if (defined(geometry)) {
geometry.geometry = PolygonPipeline.scaleToGeodeticHeight(geometry.geometry, height, ellipsoid, !perPositionHeight);
geometries.push(geometry);
}
}
}
geometry = GeometryPipeline.combine(geometries);
var center = boundingSphere.center;
scratchNormal = ellipsoid.geodeticSurfaceNormal(center, scratchNormal);
scratchPosition = Cartesian3.multiplyByScalar(scratchNormal, height, scratchPosition);
center = Cartesian3.add(center, scratchPosition, center);
if (extrude) {
scratchBoundingSphere = BoundingSphere.clone(boundingSphere, scratchBoundingSphere);
center = scratchBoundingSphere.center;
scratchPosition = Cartesian3.multiplyByScalar(scratchNormal, extrudedHeight, scratchPosition);
center = Cartesian3.add(center, scratchPosition, center);
boundingSphere = BoundingSphere.union(boundingSphere, scratchBoundingSphere, boundingSphere);
}
return new Geometry({
attributes : geometry.attributes,
indices : geometry.indices,
primitiveType : geometry.primitiveType,
boundingSphere : boundingSphere
});
};
return PolygonOutlineGeometry;
});