| 123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104 | import { isObject, isNumber, degreesToRadians, polygon } from '@turf/helpers';import rhumbDestination from '@turf/rhumb-destination';import transformRotate from '@turf/transform-rotate';import { getCoord } from '@turf/invariant';/** * Takes a {@link Point} and calculates the ellipse polygon given two semi-axes expressed in variable units and steps for precision. * * @param {Coord} center center point * @param {number} xSemiAxis semi (major) axis of the ellipse along the x-axis * @param {number} ySemiAxis semi (minor) axis of the ellipse along the y-axis * @param {Object} [options={}] Optional parameters * @param {number} [options.angle=0] angle of rotation in decimal degrees, positive clockwise * @param {Coord} [options.pivot='origin'] point around which the rotation will be performed * @param {number} [options.steps=64] number of steps * @param {string} [options.units='kilometers'] unit of measurement for axes * @param {Object} [options.properties={}] properties * @returns {Feature<Polygon>} ellipse polygon * @example * var center = [-75, 40]; * var xSemiAxis = 5; * var ySemiAxis = 2; * var ellipse = turf.ellipse(center, xSemiAxis, ySemiAxis); * * //addToMap * var addToMap = [turf.point(center), ellipse] */function ellipse(center, xSemiAxis, ySemiAxis, options) {  // Optional params  options = options || {};  var steps = options.steps || 64;  var units = options.units || "kilometers";  var angle = options.angle || 0;  var pivot = options.pivot || center;  var properties = options.properties || center.properties || {};  // validation  if (!center) throw new Error("center is required");  if (!xSemiAxis) throw new Error("xSemiAxis is required");  if (!ySemiAxis) throw new Error("ySemiAxis is required");  if (!isObject(options)) throw new Error("options must be an object");  if (!isNumber(steps)) throw new Error("steps must be a number");  if (!isNumber(angle)) throw new Error("angle must be a number");  var centerCoords = getCoord(center);  if (units === "degrees") {    var angleRad = degreesToRadians(angle);  } else {    xSemiAxis = rhumbDestination(center, xSemiAxis, 90, { units: units });    ySemiAxis = rhumbDestination(center, ySemiAxis, 0, { units: units });    xSemiAxis = getCoord(xSemiAxis)[0] - centerCoords[0];    ySemiAxis = getCoord(ySemiAxis)[1] - centerCoords[1];  }  var coordinates = [];  for (var i = 0; i < steps; i += 1) {    var stepAngle = (i * -360) / steps;    var x =      (xSemiAxis * ySemiAxis) /      Math.sqrt(        Math.pow(ySemiAxis, 2) +          Math.pow(xSemiAxis, 2) * Math.pow(getTanDeg(stepAngle), 2)      );    var y =      (xSemiAxis * ySemiAxis) /      Math.sqrt(        Math.pow(xSemiAxis, 2) +          Math.pow(ySemiAxis, 2) / Math.pow(getTanDeg(stepAngle), 2)      );    if (stepAngle < -90 && stepAngle >= -270) x = -x;    if (stepAngle < -180 && stepAngle >= -360) y = -y;    if (units === "degrees") {      var newx = x * Math.cos(angleRad) + y * Math.sin(angleRad);      var newy = y * Math.cos(angleRad) - x * Math.sin(angleRad);      x = newx;      y = newy;    }    coordinates.push([x + centerCoords[0], y + centerCoords[1]]);  }  coordinates.push(coordinates[0]);  if (units === "degrees") {    return polygon([coordinates], properties);  } else {    return transformRotate(polygon([coordinates], properties), angle, {      pivot: pivot,    });  }}/** * Get Tan Degrees * * @private * @param {number} deg Degrees * @returns {number} Tan Degrees */function getTanDeg(deg) {  var rad = (deg * Math.PI) / 180;  return Math.tan(rad);}export default ellipse;
 |