@flatten-js/core

<!DOCTYPE html> <html lang="en"> <head> <meta charset="utf-8"> <meta name="viewport" content="width=device-width,initial-scale=1"> <title>algorithms/distance.js - Documentation</title> <script src="scripts/prettify/prettify.js"></script> <script src="scripts/prettify/lang-css.js"></script>  <link type="text/css" rel="stylesheet" href="https://code.ionicframework.com/ionicons/2.0.1/css/ionicons.min.css"> <link type="text/css" rel="stylesheet" href="styles/prettify-tomorrow.css"> <link type="text/css" rel="stylesheet" href="styles/jsdoc-default.css"> </head> <body> <input type="checkbox" id="nav-trigger" class="nav-trigger" /> <label for="nav-trigger" class="navicon-button x"> <div class="navicon"></div> </label> <label for="nav-trigger" class="overlay"></label> <nav> <li class="nav-link nav-home-link"><a href="index.html">Home</a></li><li class="nav-heading">Classes</li><li class="nav-heading">C<a href="Arc.html">Arc</a></li><li class="nav-item">F<a href="Arc.html#clone">clone</a></li><li class="nav-item">F<a href="Arc.html#contains">contains</a></li><li class="nav-item">F<a href="Arc.html#split">split</a></li><li class="nav-item">F<a href="Arc.html#middle">middle</a></li><li class="nav-item">F<a href="Arc.html#chordHeight">chordHeight</a></li><li class="nav-item">F<a href="Arc.html#intersect">intersect</a></li><li class="nav-item">F<a href="Arc.html#distanceTo">distanceTo</a></li><li class="nav-item">F<a href="Arc.html#breakToFunctional">breakToFunctional</a></li><li class="nav-item">F<a href="Arc.html#tangentInStart">tangentInStart</a></li><li class="nav-item">F<a href="Arc.html#tangentInEnd">tangentInEnd</a></li><li class="nav-item">F<a href="Arc.html#reverse">reverse</a></li><li class="nav-item">F<a href="Arc.html#translate">translate</a></li><li class="nav-item">F<a href="Arc.html#rotate">rotate</a></li><li class="nav-item">F<a href="Arc.html#scale">scale</a></li><li class="nav-item">F<a href="Arc.html#transform">transform</a></li><li class="nav-item">F<a href="Arc.html#sortPoints">sortPoints</a></li><li class="nav-item">F<a href="Arc.html#toJSON">toJSON</a></li><li class="nav-item">F<a href="Arc.html#svg">svg</a></li><li class="nav-heading">C<a href="Box.html">Box</a></li><li class="nav-item">F<a href="Box.html#clone">clone</a></li><li class="nav-item">F<a href="Box.html#not_intersect">not_intersect</a></li><li class="nav-item">F<a href="Box.html#intersect">intersect</a></li><li class="nav-item">F<a href="Box.html#merge">merge</a></li><li class="nav-item">F<a href="Box.html#less_than">less_than</a></li><li class="nav-item">F<a href="Box.html#equal_to">equal_to</a></li><li class="nav-item">F<a href="Box.html#set">set</a></li><li class="nav-item">F<a href="Box.html#toPoints">toPoints</a></li><li class="nav-item">F<a href="Box.html#toSegments">toSegments</a></li><li class="nav-item">F<a href="Box.html#svg">svg</a></li><li class="nav-heading">C<a href="Circle.html">Circle</a></li><li class="nav-item">F<a href="Circle.html#clone">clone</a></li><li class="nav-item">F<a href="Circle.html#contains">contains</a></li><li class="nav-item">F<a href="Circle.html#toArc">toArc</a></li><li class="nav-item">F<a href="Circle.html#intersect">intersect</a></li><li class="nav-item">F<a href="Circle.html#distanceTo">distanceTo</a></li><li class="nav-item">F<a href="Circle.html#toJSON">toJSON</a></li><li class="nav-item">F<a href="Circle.html#svg">svg</a></li><li class="nav-heading">C<a href="Edge.html">Edge</a></li><li class="nav-item">F<a href="Edge.html#middle">middle</a></li><li class="nav-item">F<a href="Edge.html#contains">contains</a></li><li class="nav-item">F<a href="Edge.html#setInclusion">setInclusion</a></li><li class="nav-item">F<a href="Edge.html#setOverlap">setOverlap</a></li><li class="nav-heading">C<a href="Face.html">Face</a></li><li class="nav-item">F<a href="Face.html#append">append</a></li><li class="nav-item">F<a href="Face.html#insert">insert</a></li><li class="nav-item">F<a href="Face.html#remove">remove</a></li><li class="nav-item">F<a href="Face.html#reverse">reverse</a></li><li class="nav-item">F<a href="Face.html#setArcLength">setArcLength</a></li><li class="nav-item">F<a href="Face.html#area">area</a></li><li class="nav-item">F<a href="Face.html#signedArea">signedArea</a></li><li class="nav-item">F<a href="Face.html#orientation">orientation</a></li><li class="nav-item">F<a href="Face.html#isSimple">isSimple</a></li><li class="nav-item">F<a href="Face.html#findEdgeByPoint">findEdgeByPoint</a></li><li class="nav-item">F<a href="Face.html#toPolygon">toPolygon</a></li><li class="nav-item">F<a href="Face.html#svg">svg</a></li><li class="nav-heading">C<a href="Line.html">Line</a></li><li class="nav-item">F<a href="Line.html#clone">clone</a></li><li class="nav-item">F<a href="Line.html#parallelTo">parallelTo</a></li><li class="nav-item">F<a href="Line.html#incidentTo">incidentTo</a></li><li class="nav-item">F<a href="Line.html#contains">contains</a></li><li class="nav-item">F<a href="Line.html#coord">coord</a></li><li class="nav-item">F<a href="Line.html#intersect">intersect</a></li><li class="nav-item">F<a href="Line.html#distanceTo">distanceTo</a></li><li class="nav-item">F<a href="Line.html#split">split</a></li><li class="nav-item">F<a href="Line.html#sortPoints">sortPoints</a></li><li class="nav-item">F<a href="Line.html#toJSON">toJSON</a></li><li class="nav-item">F<a href="Line.html#svg">svg</a></li><li class="nav-heading">C<a href="Matrix.html">Matrix</a></li><li class="nav-item">F<a href="Matrix.html#clone">clone</a></li><li class="nav-item">F<a href="Matrix.html#transform">transform</a></li><li class="nav-item">F<a href="Matrix.html#multiply">multiply</a></li><li class="nav-item">F<a href="Matrix.html#translate">translate</a></li><li class="nav-item">F<a href="Matrix.html#rotate">rotate</a></li><li class="nav-item">F<a href="Matrix.html#scale">scale</a></li><li class="nav-item">F<a href="Matrix.html#equalTo">equalTo</a></li><li class="nav-heading">C<a href="Multiline.html">Multiline</a></li><li class="nav-item">F<a href="Multiline.html#clone">clone</a></li><li class="nav-item">F<a href="Multiline.html#addVertex">addVertex</a></li><li class="nav-item">F<a href="Multiline.html#split">split</a></li><li class="nav-item">F<a href="Multiline.html#findEdgeByPoint">findEdgeByPoint</a></li><li class="nav-item">F<a href="Multiline.html#translate">translate</a></li><li class="nav-item">F<a href="Multiline.html#rotate">rotate</a></li><li class="nav-item">F<a href="Multiline.html#transform">transform</a></li><li class="nav-item">F<a href="Multiline.html#toShapes">toShapes</a></li><li class="nav-item">F<a href="Multiline.html#toJSON">toJSON</a></li><li class="nav-item">F<a href="Multiline.html#svg">svg</a></li><li class="nav-heading">C<a href="Point.html">Point</a></li><li class="nav-item">F<a href="Point.html#clone">clone</a></li><li class="nav-item">F<a href="Point.html#equalTo">equalTo</a></li><li class="nav-item">F<a href="Point.html#lessThan">lessThan</a></li><li class="nav-item">F<a href="Point.html#rotate">rotate</a></li><li class="nav-item">F<a href="Point.html#translate">translate</a></li><li class="nav-item">F<a href="Point.html#transform">transform</a></li><li class="nav-item">F<a href="Point.html#projectionOn">projectionOn</a></li><li class="nav-item">F<a href="Point.html#leftTo">leftTo</a></li><li class="nav-item">F<a href="Point.html#distanceTo">distanceTo</a></li><li class="nav-item">F<a href="Point.html#on">on</a></li><li class="nav-item">F<a href="Point.html#toJSON">toJSON</a></li><li class="nav-item">F<a href="Point.html#svg">svg</a></li><li class="nav-heading">C<a href="Polygon.html">Polygon</a></li><li class="nav-item">F<a href="Polygon.html#clone">clone</a></li><li class="nav-item">F<a href="Polygon.html#isEmpty">isEmpty</a></li><li class="nav-item">F<a href="Polygon.html#isValid">isValid</a></li><li class="nav-item">F<a href="Polygon.html#area">area</a></li><li class="nav-item">F<a href="Polygon.html#addFace">addFace</a></li><li class="nav-item">F<a href="Polygon.html#deleteFace">deleteFace</a></li><li class="nav-item">F<a href="Polygon.html#removeChain">removeChain</a></li><li class="nav-item">F<a href="Polygon.html#addVertex">addVertex</a></li><li class="nav-item">F<a href="Polygon.html#cut">cut</a></li><li class="nav-item">F<a href="Polygon.html#cutFace">cutFace</a></li><li class="nav-item">F<a href="Polygon.html#findEdgeByPoint">findEdgeByPoint</a></li><li class="nav-item">F<a href="Polygon.html#splitToIslands">splitToIslands</a></li><li class="nav-item">F<a href="Polygon.html#reverse">reverse</a></li><li class="nav-item">F<a href="Polygon.html#contains">contains</a></li><li class="nav-item">F<a href="Polygon.html#distanceTo">distanceTo</a></li><li class="nav-item">F<a href="Polygon.html#intersect">intersect</a></li><li class="nav-item">F<a href="Polygon.html#translate">translate</a></li><li class="nav-item">F<a href="Polygon.html#rotate">rotate</a></li><li class="nav-item">F<a href="Polygon.html#transform">transform</a></li><li class="nav-item">F<a href="Polygon.html#toJSON">toJSON</a></li><li class="nav-item">F<a href="Polygon.html#toArray">toArray</a></li><li class="nav-item">F<a href="Polygon.html#svg">svg</a></li><li class="nav-heading">C<a href="Ray.html">Ray</a></li><li class="nav-item">F<a href="Ray.html#clone">clone</a></li><li class="nav-item">F<a href="Ray.html#contains">contains</a></li><li class="nav-item">F<a href="Ray.html#split">split</a></li><li class="nav-item">F<a href="Ray.html#intersect">intersect</a></li><li class="nav-item">F<a href="Ray.html#svg">svg</a></li><li class="nav-heading">C<a href="Segment.html">Segment</a></li><li class="nav-item">F<a href="Segment.html#clone">clone</a></li><li class="nav-item">F<a href="Segment.html#equalTo">equalTo</a></li><li class="nav-item">F<a href="Segment.html#contains">contains</a></li><li class="nav-item">F<a href="Segment.html#intersect">intersect</a></li><li class="nav-item">F<a href="Segment.html#distanceTo">distanceTo</a></li><li class="nav-item">F<a href="Segment.html#tangentInStart">tangentInStart</a></li><li class="nav-item">F<a href="Segment.html#tangentInEnd">tangentInEnd</a></li><li class="nav-item">F<a href="Segment.html#reverse">reverse</a></li><li class="nav-item">F<a href="Segment.html#split">split</a></li><li class="nav-item">F<a href="Segment.html#middle">middle</a></li><li class="nav-item">F<a href="Segment.html#translate">translate</a></li><li class="nav-item">F<a href="Segment.html#rotate">rotate</a></li><li class="nav-item">F<a href="Segment.html#transform">transform</a></li><li class="nav-item">F<a href="Segment.html#isZeroLength">isZeroLength</a></li><li class="nav-item">F<a href="Segment.html#sortPoints">sortPoints</a></li><li class="nav-item">F<a href="Segment.html#toJSON">toJSON</a></li><li class="nav-item">F<a href="Segment.html#svg">svg</a></li><li class="nav-heading">C<a href="Vector.html">Vector</a></li><li class="nav-item">F<a href="Vector.html#clone">clone</a></li><li class="nav-item">F<a href="Vector.html#equalTo">equalTo</a></li><li class="nav-item">F<a href="Vector.html#multiply">multiply</a></li><li class="nav-item">F<a href="Vector.html#dot">dot</a></li><li class="nav-item">F<a href="Vector.html#cross">cross</a></li><li class="nav-item">F<a href="Vector.html#normalize">normalize</a></li><li class="nav-item">F<a href="Vector.html#rotate">rotate</a></li><li class="nav-item">F<a href="Vector.html#rotate90CCW">rotate90CCW</a></li><li class="nav-item">F<a href="Vector.html#rotate90CW">rotate90CW</a></li><li class="nav-item">F<a href="Vector.html#invert">invert</a></li><li class="nav-item">F<a href="Vector.html#add">add</a></li><li class="nav-item">F<a href="Vector.html#subtract">subtract</a></li><li class="nav-item">F<a href="Vector.html#angleTo">angleTo</a></li><li class="nav-item">F<a href="Vector.html#projectionOn">projectionOn</a></li><li class="nav-item">F<a href="Vector.html#toJSON">toJSON</a></li><li class="nav-heading">C<a href="CircularLinkedList.html">CircularLinkedList</a></li><li class="nav-item">F<a href="CircularLinkedList.html#append">append</a></li><li class="nav-item">F<a href="CircularLinkedList.html#insert">insert</a></li><li class="nav-item">F<a href="CircularLinkedList.html#remove">remove</a></li><li class="nav-heading">C<a href="DE9IM.html">DE9IM</a></li><li class="nav-item">F<a href="DE9IM.html#toString">toString</a></li><li class="nav-heading">C<a href="LinkedList.html">LinkedList</a></li><li class="nav-item">F<a href="LinkedList.html#toArray">toArray</a></li><li class="nav-item">F<a href="LinkedList.html#append">append</a></li><li class="nav-item">F<a href="LinkedList.html#insert">insert</a></li><li class="nav-item">F<a href="LinkedList.html#remove">remove</a></li><li class="nav-item">F<a href="LinkedList.html#isEmpty">isEmpty</a></li><li class="nav-item">F<a href="LinkedList.html#.testInfiniteLoop">testInfiniteLoop</a></li><li class="nav-heading">C<a href="PlanarSet.html">PlanarSet</a></li><li class="nav-item">F<a href="PlanarSet.html#add">add</a></li><li class="nav-item">F<a href="PlanarSet.html#delete">delete</a></li><li class="nav-item">F<a href="PlanarSet.html#clear">clear</a></li><li class="nav-item">F<a href="PlanarSet.html#search">search</a></li><li class="nav-item">F<a href="PlanarSet.html#hit">hit</a></li><li class="nav-item">F<a href="PlanarSet.html#svg">svg</a></li><li class="nav-heading">C<a href="Errors.html">Errors</a></li><li class="nav-heading">Modules</li><li class="nav-heading">M<a href="module-BooleanOperations.html">BooleanOperations</a></li><li class="nav-item">F<a href="module-BooleanOperations.html#.unify">unify</a></li><li class="nav-item">F<a href="module-BooleanOperations.html#.subtract">subtract</a></li><li class="nav-item">F<a href="module-BooleanOperations.html#.intersect">intersect</a></li><li class="nav-item">F<a href="module-BooleanOperations.html#.innerClip">innerClip</a></li><li class="nav-item">F<a href="module-BooleanOperations.html#.outerClip">outerClip</a></li><li class="nav-item">F<a href="module-BooleanOperations.html#.calculateIntersections">calculateIntersections</a></li><li class="nav-heading">M<a href="module-RayShoot.html">RayShoot</a></li><li class="nav-item">F<a href="module-RayShoot.html#.ray_shoot">ray_shoot</a></li><li class="nav-heading">M<a href="module-Relation.html">Relation</a></li><li class="nav-item">F<a href="module-Relation.html#.equal">equal</a></li><li class="nav-item">F<a href="module-Relation.html#.intersect">intersect</a></li><li class="nav-item">F<a href="module-Relation.html#.touch">touch</a></li><li class="nav-item">F<a href="module-Relation.html#.disjoint">disjoint</a></li><li class="nav-item">F<a href="module-Relation.html#.inside">inside</a></li><li class="nav-item">F<a href="module-Relation.html#.covered">covered</a></li><li class="nav-item">F<a href="module-Relation.html#.contain">contain</a></li><li class="nav-item">F<a href="module-Relation.html#.cover">cover</a></li><li class="nav-item">F<a href="module-Relation.html#.relate">relate</a></li><li class="nav-heading"><a href="global.html">Globals</a></li><li class="nav-item">F<a href="global.html#inverse">inverse</a></li><li class="nav-item">M<a href="global.html#arc">arc</a></li><li class="nav-item">M<a href="global.html#box">box</a></li><li class="nav-item">M<a href="global.html#circle">circle</a></li><li class="nav-item">M<a href="global.html#line">line</a></li><li class="nav-item">M<a href="global.html#matrix">matrix</a></li><li class="nav-item">M<a href="global.html#multiline">multiline</a></li><li class="nav-item">M<a href="global.html#point">point</a></li><li class="nav-item">M<a href="global.html#polygon">polygon</a></li><li class="nav-item">M<a href="global.html#segment">segment</a></li><li class="nav-item">M<a href="global.html#vector">vector</a></li><li class="nav-item">M<a href="global.html#CCW">CCW</a></li><li class="nav-item">M<a href="global.html#CW">CW</a></li><li class="nav-item">M<a href="global.html#ORIENTATION">ORIENTATION</a></li><li class="nav-item">M<a href="global.html#DP_TOL">DP_TOL</a></li><li class="nav-item">F<a href="global.html#setTolerance">setTolerance</a></li><li class="nav-item">F<a href="global.html#getTolerance">getTolerance</a></li><li class="nav-item">F<a href="global.html#EQ_0">EQ_0</a></li><li class="nav-item">F<a href="global.html#EQ">EQ</a></li><li class="nav-item">F<a href="global.html#GT">GT</a></li><li class="nav-item">F<a href="global.html#GE">GE</a></li><li class="nav-item">F<a href="global.html#LT">LT</a></li><li class="nav-item">F<a href="global.html#LE">LE</a></li> </nav> <div id="main"> <h1 class="page-title">algorithms/distance.js</h1> <section> <article> <pre class="prettyprint source linenums"><code>"use strict"; // let IntervalTree = require('flatten-interval-tree'); import IntervalTree from '@flatten-js/interval-tree'; import * as Intersection from '../algorithms/intersection'; import Flatten from '../flatten'; export class Distance { /** * Calculate distance and shortest segment between points * @param pt1 * @param pt2 * @returns {Number | Segment} - distance and shortest segment */ static point2point(pt1, pt2) { return pt1.distanceTo(pt2); } /** * Calculate distance and shortest segment between point and line * @param pt * @param line * @returns {Number | Segment} - distance and shortest segment */ static point2line(pt, line) { let closest_point = pt.projectionOn(line); let vec = new Flatten.Vector(pt, closest_point); return [vec.length, new Flatten.Segment(pt, closest_point)]; } /** * Calculate distance and shortest segment between point and circle * @param pt * @param circle * @returns {Number | Segment} - distance and shortest segment */ static point2circle(pt, circle) { let [dist2center, shortest_dist] = pt.distanceTo(circle.center); if (Flatten.Utils.EQ_0(dist2center)) { return [circle.r, new Flatten.Segment(pt, circle.toArc().start)]; } else { let dist = Math.abs(dist2center - circle.r); let v = new Flatten.Vector(circle.pc, pt).normalize().multiply(circle.r); let closest_point = circle.pc.translate(v); return [dist, new Flatten.Segment(pt, closest_point)]; } } /** * Calculate distance and shortest segment between point and segment * @param pt * @param segment * @returns {Number | Segment} - distance and shortest segment */ static point2segment(pt, segment) { /* Degenerated case of zero-length segment */ if (segment.start.equalTo(segment.end)) { return Distance.point2point(pt, segment.start); } let v_seg = new Flatten.Vector(segment.start, segment.end); let v_ps2pt = new Flatten.Vector(segment.start, pt); let v_pe2pt = new Flatten.Vector(segment.end, pt); let start_sp = v_seg.dot(v_ps2pt); /* dot product v_seg * v_ps2pt */ let end_sp = -v_seg.dot(v_pe2pt); /* minus dot product v_seg * v_pe2pt */ let dist; let closest_point; if (Flatten.Utils.GE(start_sp, 0) && Flatten.Utils.GE(end_sp, 0)) { /* point inside segment scope */ let v_unit = segment.tangentInStart(); // new Flatten.Vector(v_seg.x / this.length, v_seg.y / this.length); /* unit vector ||v_unit|| = 1 */ dist = Math.abs(v_unit.cross(v_ps2pt)); /* dist = abs(v_unit x v_ps2pt) */ closest_point = segment.start.translate(v_unit.multiply(v_unit.dot(v_ps2pt))); return [dist, new Flatten.Segment(pt, closest_point)]; } else if (start_sp < 0) { /* point is out of scope closer to ps */ return pt.distanceTo(segment.start); } else { /* point is out of scope closer to pe */ return pt.distanceTo(segment.end); } }; /** * Calculate distance and shortest segment between point and arc * @param pt * @param arc * @returns {Number | Segment} - distance and shortest segment */ static point2arc(pt, arc) { let circle = new Flatten.Circle(arc.pc, arc.r); let dist_and_segment = []; let dist, shortest_segment; [dist, shortest_segment] = Distance.point2circle(pt, circle); if (shortest_segment.end.on(arc)) { dist_and_segment.push(Distance.point2circle(pt, circle)); } dist_and_segment.push(Distance.point2point(pt, arc.start)); dist_and_segment.push(Distance.point2point(pt, arc.end)); Distance.sort(dist_and_segment); return dist_and_segment[0]; } /** * Calculate distance and shortest segment between segment and line * @param seg * @param line * @returns {Number | Segment} */ static segment2line(seg, line) { let ip = seg.intersect(line); if (ip.length > 0) { return [0, new Flatten.Segment(ip[0], ip[0])]; // distance = 0, closest point is the first point } let dist_and_segment = []; dist_and_segment.push(Distance.point2line(seg.start, line)); dist_and_segment.push(Distance.point2line(seg.end, line)); Distance.sort(dist_and_segment); return dist_and_segment[0]; } /** * Calculate distance and shortest segment between two segments * @param seg1 * @param seg2 * @returns {Number | Segment} - distance and shortest segment */ static segment2segment(seg1, seg2) { let ip = Intersection.intersectSegment2Segment(seg1, seg2); if (ip.length > 0) { return [0, new Flatten.Segment(ip[0], ip[0])]; // distance = 0, closest point is the first point } // Seg1 and seg2 not intersected let dist_and_segment = []; let dist_tmp, shortest_segment_tmp [dist_tmp, shortest_segment_tmp] = Distance.point2segment(seg2.start, seg1); dist_and_segment.push([dist_tmp, shortest_segment_tmp.reverse()]); [dist_tmp, shortest_segment_tmp] = Distance.point2segment(seg2.end, seg1); dist_and_segment.push([dist_tmp, shortest_segment_tmp.reverse()]); dist_and_segment.push(Distance.point2segment(seg1.start, seg2)); dist_and_segment.push(Distance.point2segment(seg1.end, seg2)); Distance.sort(dist_and_segment); return dist_and_segment[0]; } /** * Calculate distance and shortest segment between segment and circle * @param seg * @param circle * @returns {Number | Segment} - distance and shortest segment */ static segment2circle(seg, circle) { /* Case 1 Segment and circle intersected. Return the first point and zero distance */ let ip = seg.intersect(circle); if (ip.length > 0) { return [0, new Flatten.Segment(ip[0], ip[0])]; } // No intersection between segment and circle /* Case 2. Distance to projection of center point to line bigger than radius * And projection point belong to segment * Then measure again distance from projection to circle and return it */ let line = new Flatten.Line(seg.ps, seg.pe); let [dist, shortest_segment] = Distance.point2line(circle.center, line); if (Flatten.Utils.GE(dist, circle.r) && shortest_segment.end.on(seg)) { return Distance.point2circle(shortest_segment.end, circle); } /* Case 3. Otherwise closest point is one of the end points of the segment */ else { let [dist_from_start, shortest_segment_from_start] = Distance.point2circle(seg.start, circle); let [dist_from_end, shortest_segment_from_end] = Distance.point2circle(seg.end, circle); return Flatten.Utils.LT(dist_from_start, dist_from_end) ? [dist_from_start, shortest_segment_from_start] : [dist_from_end, shortest_segment_from_end]; } } /** * Calculate distance and shortest segment between segment and arc * @param seg * @param arc * @returns {Number | Segment} - distance and shortest segment */ static segment2arc(seg, arc) { /* Case 1 Segment and arc intersected. Return the first point and zero distance */ let ip = seg.intersect(arc); if (ip.length > 0) { return [0, new Flatten.Segment(ip[0], ip[0])]; } // No intersection between segment and arc let line = new Flatten.Line(seg.ps, seg.pe); let circle = new Flatten.Circle(arc.pc, arc.r); /* Case 2. Distance to projection of center point to line bigger than radius AND * projection point belongs to segment AND * distance from projection point to circle belongs to arc => * return this distance from projection to circle */ let [dist_from_center, shortest_segment_from_center] = Distance.point2line(circle.center, line); if (Flatten.Utils.GE(dist_from_center, circle.r) && shortest_segment_from_center.end.on(seg)) { let [dist_from_projection, shortest_segment_from_projection] = Distance.point2circle(shortest_segment_from_center.end, circle); if (shortest_segment_from_projection.end.on(arc)) { return [dist_from_projection, shortest_segment_from_projection]; } } /* Case 3. Otherwise closest point is one of the end points of the segment */ let dist_and_segment = []; dist_and_segment.push(Distance.point2arc(seg.start, arc)); dist_and_segment.push(Distance.point2arc(seg.end, arc)); let dist_tmp, segment_tmp; [dist_tmp, segment_tmp] = Distance.point2segment(arc.start, seg); dist_and_segment.push([dist_tmp, segment_tmp.reverse()]); [dist_tmp, segment_tmp] = Distance.point2segment(arc.end, seg); dist_and_segment.push([dist_tmp, segment_tmp.reverse()]); Distance.sort(dist_and_segment); return dist_and_segment[0]; } /** * Calculate distance and shortest segment between two circles * @param circle1 * @param circle2 * @returns {Number | Segment} - distance and shortest segment */ static circle2circle(circle1, circle2) { let ip = circle1.intersect(circle2); if (ip.length > 0) { return [0, new Flatten.Segment(ip[0], ip[0])]; } // Case 1. Concentric circles. Convert to arcs and take distance between two arc starts if (circle1.center.equalTo(circle2.center)) { let arc1 = circle1.toArc(); let arc2 = circle2.toArc(); return Distance.point2point(arc1.start, arc2.start); } else { // Case 2. Not concentric circles let line = new Flatten.Line(circle1.center, circle2.center); let ip1 = line.intersect(circle1); let ip2 = line.intersect(circle2); let dist_and_segment = []; dist_and_segment.push(Distance.point2point(ip1[0], ip2[0])); dist_and_segment.push(Distance.point2point(ip1[0], ip2[1])); dist_and_segment.push(Distance.point2point(ip1[1], ip2[0])); dist_and_segment.push(Distance.point2point(ip1[1], ip2[1])); Distance.sort(dist_and_segment); return dist_and_segment[0]; } } /** * Calculate distance and shortest segment between two circles * @param circle * @param line * @returns {Number | Segment} - distance and shortest segment */ static circle2line(circle, line) { let ip = circle.intersect(line); if (ip.length > 0) { return [0, new Flatten.Segment(ip[0], ip[0])]; } let [dist_from_center, shortest_segment_from_center] = Distance.point2line(circle.center, line); let [dist, shortest_segment] = Distance.point2circle(shortest_segment_from_center.end, circle); shortest_segment = shortest_segment.reverse(); return [dist, shortest_segment]; } /** * Calculate distance and shortest segment between arc and line * @param arc * @param line * @returns {Number | Segment} - distance and shortest segment */ static arc2line(arc, line) { /* Case 1 Line and arc intersected. Return the first point and zero distance */ let ip = line.intersect(arc); if (ip.length > 0) { return [0, new Flatten.Segment(ip[0], ip[0])]; } let circle = new Flatten.Circle(arc.center, arc.r); /* Case 2. Distance to projection of center point to line bigger than radius AND * projection point belongs to segment AND * distance from projection point to circle belongs to arc => * return this distance from projection to circle */ let [dist_from_center, shortest_segment_from_center] = Distance.point2line(circle.center, line); if (Flatten.Utils.GE(dist_from_center, circle.r)) { let [dist_from_projection, shortest_segment_from_projection] = Distance.point2circle(shortest_segment_from_center.end, circle); if (shortest_segment_from_projection.end.on(arc)) { return [dist_from_projection, shortest_segment_from_projection]; } } else { let dist_and_segment = []; dist_and_segment.push(Distance.point2line(arc.start, line)); dist_and_segment.push(Distance.point2line(arc.end, line)); Distance.sort(dist_and_segment); return dist_and_segment[0]; } } /** * Calculate distance and shortest segment between arc and circle * @param arc * @param circle2 * @returns {Number | Segment} - distance and shortest segment */ static arc2circle(arc, circle2) { let ip = arc.intersect(circle2); if (ip.length > 0) { return [0, new Flatten.Segment(ip[0], ip[0])]; } let circle1 = new Flatten.Circle(arc.center, arc.r); let [dist, shortest_segment] = Distance.circle2circle(circle1, circle2); if (shortest_segment.start.on(arc)) { return [dist, shortest_segment]; } else { let dist_and_segment = []; dist_and_segment.push(Distance.point2circle(arc.start, circl