dgeoutils 2.4.48 → 2.5.0

package/dist/cjs/DCircle.d.ts

@@ -0,0 +1,18 @@
+import { DPoint } from './DPoint';
+import { DPolygon } from './DPolygon';
+export declare class DCircle {
+    center: DPoint;
+    r: number;
+    constructor(center?: DPoint, r?: number);
+    toString(): string;
+    getValue(): {
+        center: DPoint;
+        r: number;
+    };
+    clone(): DCircle;
+    findPoints(c: DCircle): DPoint[] | number;
+    equal({ center, r }: DCircle): boolean;
+    findPolygonInside(pointCount?: number, startAngle?: number, stopAngle?: number): DPolygon;
+    findPolygonInsideOnSphere(pointCount?: number, startAngle?: number, stopAngle?: number): DPolygon;
+    private sphereOffset;
+}

package/dist/cjs/DCircle.js

@@ -0,0 +1,102 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.DCircle = void 0;
+var DPoint_1 = require("./DPoint");
+var DPolygon_1 = require("./DPolygon");
+var utils_1 = require("./utils");
+var DCircle = (function () {
+    function DCircle(center, r) {
+        if (center === void 0) { center = DPoint_1.DPoint.zero(); }
+        if (r === void 0) { r = 0; }
+        this.center = center;
+        this.r = r;
+    }
+    DCircle.prototype.toString = function () {
+        return "(".concat(this.center.toString(), ", ").concat(this.r, ")");
+    };
+    DCircle.prototype.getValue = function () {
+        return { center: this.center, r: this.r };
+    };
+    DCircle.prototype.clone = function () {
+        return new DCircle(this.center, this.r);
+    };
+    DCircle.prototype.findPoints = function (c) {
+        if (this.equal(c)) {
+            return Infinity;
+        }
+        var _a = this, _b = _a.center, x0 = _b.x, y0 = _b.y, r0 = _a.r;
+        var _c = c.center, x1 = _c.x, y1 = _c.y, r1 = c.r;
+        var r02 = r0 * r0;
+        var d = this.center.distance(c.center);
+        var a = (r02 - r1 * r1 + d * d) / (2 * d);
+        var h = Math.sqrt(r02 - a * a);
+        var ad = a / d;
+        var dy = y1 - y0;
+        var dx = x1 - x0;
+        var hd = h / d;
+        var x2 = x0 + ad * (x1 - x0);
+        var y2 = y0 + ad * (y1 - y0);
+        var x31 = x2 + hd * dy;
+        var y31 = y2 - hd * dx;
+        var x32 = x2 - hd * dy;
+        var y32 = y2 + hd * dx;
+        var res = [];
+        if (!isNaN(x31) && !isNaN(y31)) {
+            res.push(new DPoint_1.DPoint(x31, y31));
+        }
+        if (!isNaN(x32) && !isNaN(y32) && !(x31 === x32 && y31 === y32)) {
+            res.push(new DPoint_1.DPoint(x32, y32));
+        }
+        return res;
+    };
+    DCircle.prototype.equal = function (_a) {
+        var center = _a.center, r = _a.r;
+        return this.center.equal(center) && this.r === r;
+    };
+    DCircle.prototype.findPolygonInside = function (pointCount, startAngle, stopAngle) {
+        if (pointCount === void 0) { pointCount = 64; }
+        if (startAngle === void 0) { startAngle = 0; }
+        if (stopAngle === void 0) { stopAngle = 2 * Math.PI; }
+        var step = 2 * Math.PI / pointCount;
+        var points = [];
+        var angle = startAngle;
+        while (angle < stopAngle - step) {
+            points.push(new DPoint_1.DPoint(this.r).scale(Math.cos(angle), Math.sin(angle))
+                .move(this.center));
+            angle += step;
+        }
+        var x = this.r * Math.cos(stopAngle) + this.center.x;
+        var y = this.r * Math.sin(stopAngle) + this.center.y;
+        points.push(new DPoint_1.DPoint(x, y));
+        return new DPolygon_1.DPolygon(points);
+    };
+    DCircle.prototype.findPolygonInsideOnSphere = function (pointCount, startAngle, stopAngle) {
+        if (pointCount === void 0) { pointCount = 64; }
+        if (startAngle === void 0) { startAngle = 0; }
+        if (stopAngle === void 0) { stopAngle = 2 * Math.PI; }
+        (0, utils_1.checkFunction)('findPolygonInsideOnSphere')
+            .checkArgument('center')
+            .shouldBeDegree(this.center);
+        var step = 2 * Math.PI / pointCount;
+        var points = [];
+        var angle = startAngle;
+        while (angle < stopAngle - step) {
+            points.push(this.sphereOffset(angle));
+            angle += step;
+        }
+        points.push(this.sphereOffset(stopAngle));
+        return new DPolygon_1.DPolygon(points);
+    };
+    DCircle.prototype.sphereOffset = function (bearing, earthRadius) {
+        if (earthRadius === void 0) { earthRadius = DPoint_1.EARTH_RADIUS_IN_METERS; }
+        var _a = this.center.clone().degreeToRadians(), lon1 = _a.x, lat1 = _a.y;
+        var dByR = this.r / earthRadius;
+        var lat = Math.asin(Math.sin(lat1) * Math.cos(dByR) +
+            Math.cos(lat1) * Math.sin(dByR) * Math.cos(bearing));
+        var lon = lon1 +
+            Math.atan2(Math.sin(bearing) * Math.sin(dByR) * Math.cos(lat1), Math.cos(dByR) - Math.sin(lat1) * Math.sin(lat));
+        return new DPoint_1.DPoint(lon, lat).radiansToDegrees();
+    };
+    return DCircle;
+}());
+exports.DCircle = DCircle;

package/dist/cjs/DLine.d.ts

@@ -0,0 +1,39 @@
+import { DPoint } from './DPoint';
+import { DCircle } from './DCircle';
+import { DPolygon } from './DPolygon';
+export declare class DLine {
+    a: number;
+    b: number;
+    c: number;
+    begin: DPoint;
+    end: DPoint;
+    constructor(a: number, b: number, c: number, begin?: DPoint, end?: DPoint);
+    clone(): DLine;
+    findPerpendicular(p: DPoint): DLine;
+    perpendicularDistance(p: DPoint): number;
+    intersection(l: DLine, d?: number, includeOnly?: boolean): DPoint | null;
+    intersectionWithCircle(circle: DCircle): DPoint | [DPoint, DPoint] | null;
+    inRange(p: DPoint, d?: number): boolean;
+    insideRange(p: DPoint, d?: number): boolean;
+    get center(): DPoint;
+    get minX(): number;
+    get minY(): number;
+    get maxX(): number;
+    get maxY(): number;
+    toString(): string;
+    getValue(): number[];
+    x(p: DPoint): DPoint;
+    y(p: DPoint): DPoint;
+    findPoint(l: DLine): DPoint | null;
+    get isParallel(): boolean;
+    get isParallelY(): boolean;
+    get isParallelX(): boolean;
+    get points(): [DPoint, DPoint];
+    getFi(): number;
+    toWKT(): string;
+    movePoint(point: DPoint, distance: number): DPoint;
+    movePoint(points: DPoint[], distances: number[]): DPoint[];
+    findFi({ a, b }: DLine, delta?: number): number;
+    vectorProduct({ a, b, c }: DLine): DLine;
+    bresenhamsLine(): DPolygon;
+}

package/dist/cjs/DLine.js

@@ -0,0 +1,343 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.DLine = void 0;
+var DPoint_1 = require("./DPoint");
+var utils_1 = require("./utils");
+var DNumbers_1 = require("./DNumbers");
+var DPolygon_1 = require("./DPolygon");
+var DLine = (function () {
+    function DLine(a, b, c, begin, end) {
+        if (begin === void 0) { begin = DPoint_1.DPoint.zero(); }
+        if (end === void 0) { end = DPoint_1.DPoint.zero(); }
+        this.a = a;
+        this.b = b;
+        this.c = c;
+        this.begin = begin;
+        this.end = end;
+    }
+    DLine.prototype.clone = function () {
+        return new DLine(this.a, this.b, this.c, this.begin.clone(), this.end.clone());
+    };
+    DLine.prototype.findPerpendicular = function (p) {
+        (0, utils_1.checkFunction)('findPerpendicular')
+            .checkArgument('this.begin')
+            .shouldBeMeters(this.begin)
+            .checkArgument('this.end')
+            .shouldBeMeters(this.end)
+            .checkArgument('p')
+            .shouldBeMeters(p);
+        return new DLine(-this.b, this.a, this.b * p.x - this.a * p.y);
+    };
+    DLine.prototype.perpendicularDistance = function (p) {
+        (0, utils_1.checkFunction)('perpendicularDistance')
+            .checkArgument('this.begin')
+            .shouldBeMeters(this.begin)
+            .checkArgument('this.end')
+            .shouldBeMeters(this.end)
+            .checkArgument('p')
+            .shouldBeMeters(p);
+        var perpendicularLine = this.findPerpendicular(p);
+        var targetPoint = perpendicularLine.findPoint(this);
+        return targetPoint.distance(p);
+    };
+    DLine.prototype.intersection = function (l, d, includeOnly) {
+        if (d === void 0) { d = 0; }
+        if (includeOnly === void 0) { includeOnly = false; }
+        var p = this.findPoint(l);
+        if (p) {
+            if (includeOnly) {
+                return this.insideRange(p, d) && l.insideRange(p, d) ? p : null;
+            }
+            return this.inRange(p, d) && l.inRange(p, d) ? p : null;
+        }
+        return null;
+    };
+    DLine.prototype.intersectionWithCircle = function (circle) {
+        var center = circle.center, r = circle.r;
+        var per = this.findPerpendicular(center);
+        var t = this.findPoint(per);
+        var distance = t.distance(center);
+        if (this.begin.equal(center)) {
+            distance = 0;
+        }
+        if (this.end.equal(center)) {
+            distance = 0;
+        }
+        if (distance < r) {
+            var _a = this, a = _a.a, b = _a.b, c = _a.c;
+            if (this.isParallel) {
+                var ct = center.distance(t);
+                var move = Math.sqrt(r * r - ct * ct);
+                if (this.isParallelY) {
+                    t.x = this.begin.x;
+                    var r1_1 = t.clone().move(0, -move);
+                    var r2_1 = t.clone().move(0, move);
+                    return [r1_1, r2_1];
+                }
+                if (this.isParallelX) {
+                    t.y = this.begin.y;
+                    var r1_2 = t.clone().move(move, 0);
+                    var r2_2 = t.clone().move(-move, 0);
+                    return [r1_2, r2_2];
+                }
+            }
+            if (this.begin.like(center)) {
+                var p = this.begin.clone();
+                return [this.movePoint(p, r), this.movePoint(p, -r)];
+            }
+            if (this.end.like(center)) {
+                var p = this.end.clone();
+                return [this.movePoint(p, r), this.movePoint(p, -r)];
+            }
+            var s = a * a + b * b;
+            var d = r * r - c * c / s;
+            var mult = Math.sqrt(d / s);
+            var r1 = t.clone().move(b * mult, -a * mult);
+            var r2 = t.clone().move(-b * mult, a * mult);
+            return [r1, r2];
+        }
+        if (distance === r) {
+            return t;
+        }
+        return null;
+    };
+    DLine.prototype.inRange = function (p, d) {
+        if (d === void 0) { d = 0; }
+        var _a = this, minX = _a.minX, minY = _a.minY, maxX = _a.maxX, maxY = _a.maxY;
+        var x = p.x, y = p.y;
+        var isInX = x >= minX - d && x <= maxX + d;
+        var isInY = y >= minY - d && y <= maxY + d;
+        return isInX && isInY;
+    };
+    DLine.prototype.insideRange = function (p, d) {
+        if (d === void 0) { d = 0; }
+        var _a = this, begin = _a.begin, end = _a.end;
+        return this.inRange(p, d) && !begin.like(p, 0.00001) && !end.like(p, 0.00001);
+    };
+    Object.defineProperty(DLine.prototype, "center", {
+        get: function () {
+            return this.begin
+                .clone()
+                .setIfLessThan(this.end)
+                .move(this.end
+                .clone()
+                .move(this.begin
+                .clone()
+                .minus())
+                .abs()
+                .minus()
+                .divide(2));
+        },
+        enumerable: false,
+        configurable: true
+    });
+    Object.defineProperty(DLine.prototype, "minX", {
+        get: function () {
+            return Math.min(this.begin.x, this.end.x);
+        },
+        enumerable: false,
+        configurable: true
+    });
+    Object.defineProperty(DLine.prototype, "minY", {
+        get: function () {
+            return Math.min(this.begin.y, this.end.y);
+        },
+        enumerable: false,
+        configurable: true
+    });
+    Object.defineProperty(DLine.prototype, "maxX", {
+        get: function () {
+            return Math.max(this.begin.x, this.end.x);
+        },
+        enumerable: false,
+        configurable: true
+    });
+    Object.defineProperty(DLine.prototype, "maxY", {
+        get: function () {
+            return Math.max(this.begin.y, this.end.y);
+        },
+        enumerable: false,
+        configurable: true
+    });
+    DLine.prototype.toString = function () {
+        return "(".concat(this.a, ", ").concat(this.b, ", ").concat(this.c, ")");
+    };
+    DLine.prototype.getValue = function () {
+        return [this.a, this.b, this.c];
+    };
+    DLine.prototype.x = function (p) {
+        if (this.isParallelY) {
+            return new DPoint_1.DPoint(-this.c / this.a, p.y);
+        }
+        if (this.isParallelX) {
+            return new DPoint_1.DPoint(p.x, -this.c / this.b);
+        }
+        return new DPoint_1.DPoint(-this.b / this.a * p.y - this.c / this.a, p.y);
+    };
+    DLine.prototype.y = function (p) {
+        if (this.isParallelY) {
+            return new DPoint_1.DPoint(-this.c / this.a, p.y);
+        }
+        if (this.isParallelX) {
+            return new DPoint_1.DPoint(p.x, -this.c / this.b);
+        }
+        return new DPoint_1.DPoint(p.x, -this.a / this.b * p.x - this.c / this.b);
+    };
+    DLine.prototype.findPoint = function (l) {
+        if (this.isParallelY && l.isParallelY) {
+            return null;
+        }
+        if (this.isParallelX && l.isParallelX) {
+            return null;
+        }
+        if (this.isParallelX && l.isParallelY) {
+            return new DPoint_1.DPoint(-l.c / l.a, -this.c / this.b);
+        }
+        if (this.isParallelY && l.isParallelX) {
+            return new DPoint_1.DPoint(-this.c / this.a, -l.c / l.b);
+        }
+        if (this.isParallelY) {
+            var x = -this.c / this.a;
+            return l.y(new DPoint_1.DPoint(x));
+        }
+        if (this.isParallelX) {
+            var y = -this.c / this.b;
+            return l.x(new DPoint_1.DPoint(0, y));
+        }
+        if (l.isParallelY) {
+            var x = -l.c / l.a;
+            return this.y(new DPoint_1.DPoint(x));
+        }
+        if (l.isParallelX) {
+            var y = -l.c / l.b;
+            return this.x(new DPoint_1.DPoint(0, y));
+        }
+        var res = this.y(new DPoint_1.DPoint((l.c / l.b - this.c / this.b) / (this.a / this.b - l.a / l.b)));
+        if (!isFinite(res.x) && !isFinite(res.y)) {
+            return null;
+        }
+        return res;
+    };
+    Object.defineProperty(DLine.prototype, "isParallel", {
+        get: function () {
+            return this.isParallelX || this.isParallelY;
+        },
+        enumerable: false,
+        configurable: true
+    });
+    Object.defineProperty(DLine.prototype, "isParallelY", {
+        get: function () {
+            return Math.abs(this.b) < 0.001;
+        },
+        enumerable: false,
+        configurable: true
+    });
+    Object.defineProperty(DLine.prototype, "isParallelX", {
+        get: function () {
+            return Math.abs(this.a) < 0.001;
+        },
+        enumerable: false,
+        configurable: true
+    });
+    Object.defineProperty(DLine.prototype, "points", {
+        get: function () {
+            return [this.begin, this.end];
+        },
+        enumerable: false,
+        configurable: true
+    });
+    DLine.prototype.getFi = function () {
+        (0, utils_1.checkFunction)('getFi')
+            .checkArgument('this.begin')
+            .shouldBeMeters(this.begin)
+            .checkArgument('this.end')
+            .shouldBeMeters(this.end);
+        var _a = this.end.clone().move(this.begin.clone().minus()), x = _a.x, y = _a.y;
+        var v = Math.atan2(y, x) - Math.PI;
+        if (v > 0) {
+            v = Math.PI - v;
+        }
+        return (Math.PI - v) % (Math.PI * 2);
+    };
+    DLine.prototype.toWKT = function () {
+        var _a = this, _b = _a.begin, x1 = _b.x, y1 = _b.y, _c = _a.end, x2 = _c.x, y2 = _c.y;
+        return "LINESTRING (".concat(x1, " ").concat(y1, ", ").concat(x2, " ").concat(y2, ")");
+    };
+    DLine.prototype.movePoint = function (i, k) {
+        var isArray = Array.isArray(i);
+        var p = isArray ? i : [i];
+        var d = (isArray ? k : [k]);
+        var fi = this.findFi(new DLine(1, 0, 0));
+        var td = this.x(new DPoint_1.DPoint(1, 1)).distance(this.x(new DPoint_1.DPoint(2, 2))) / 2;
+        var sinCos = new DPoint_1.DPoint(Math.sin(fi), Math.cos(fi));
+        var dt = sinCos.clone().scale(td);
+        var p1T = p[0].clone().move(dt.clone().minus());
+        var p2T = p[0].clone().move(dt);
+        var res = [];
+        if (DNumbers_1.DNumbers.like(this.y(p1T).y, p1T.y) || DNumbers_1.DNumbers.like(this.y(p2T).y, p2T.y)) {
+            res = p.map(function (t, index) { return t.clone()
+                .move(sinCos.clone().scale(d[index])); });
+        }
+        else {
+            res = p.map(function (t, index) { return t.clone()
+                .move(sinCos.clone().scale(d[index])
+                .setX(function (_a) {
+                var x = _a.x;
+                return -x;
+            })); });
+        }
+        return isArray ? res : res[0];
+    };
+    DLine.prototype.findFi = function (_a, delta) {
+        var a = _a.a, b = _a.b;
+        if (delta === void 0) { delta = 1.0001; }
+        var _b = this, q = _b.a, w = _b.b;
+        var val = (q * a + w * b) / (Math.sqrt(q * q + w * w) * Math.sqrt(a * a + b * b));
+        if (val > 1 && val < delta) {
+            val = 1;
+        }
+        else if (val < -1 && val > -delta) {
+            val = -1;
+        }
+        return Math.acos(val);
+    };
+    DLine.prototype.vectorProduct = function (_a) {
+        var a = _a.a, b = _a.b, c = _a.c;
+        var _b = this, q = _b.a, w = _b.b, e = _b.c;
+        return new DLine(w * c - e * b, e * a - q * c, q * b - w * a);
+    };
+    DLine.prototype.bresenhamsLine = function () {
+        var _a = this.begin, x0 = _a.x, y0 = _a.y;
+        var _b = this.end, x1 = _b.x, y1 = _b.y;
+        var dx = Math.abs(x1 - x0);
+        var sx = x0 < x1 ? 1 : -1;
+        var dy = -Math.abs(y1 - y0);
+        var sy = y0 < y1 ? 1 : -1;
+        var error = dx + dy;
+        var res = new DPolygon_1.DPolygon();
+        while (true) {
+            res.push(new DPoint_1.DPoint(x0, y0));
+            if (x0 === x1 && y0 === y1) {
+                break;
+            }
+            var e2 = 2 * error;
+            if (e2 >= dy) {
+                if (x0 === x1) {
+                    break;
+                }
+                error += dy;
+                x0 += sx;
+            }
+            if (e2 <= dx) {
+                if (y0 === y1) {
+                    break;
+                }
+                error += dx;
+                y0 += sy;
+            }
+        }
+        return res;
+    };
+    return DLine;
+}());
+exports.DLine = DLine;

package/dist/cjs/DNumbers.d.ts

@@ -0,0 +1,8 @@
+export declare class DNumbers {
+    static like(v: number, s: number, d?: number): boolean;
+    static likeZero(v: number): boolean;
+    static like2PI(v: number): boolean;
+    static likePI(v: number): boolean;
+    static rad2Deg(v: number): number;
+    static deg2Rad(v: number): number;
+}

package/dist/cjs/DNumbers.js

@@ -0,0 +1,30 @@
+"use strict";
+Object.defineProperty(exports, "__esModule", { value: true });
+exports.DNumbers = void 0;
+var DPoint_1 = require("./DPoint");
+var delta = 0.001;
+var DNumbers = (function () {
+    function DNumbers() {
+    }
+    DNumbers.like = function (v, s, d) {
+        if (d === void 0) { d = delta; }
+        return Math.abs(v - s) < d;
+    };
+    DNumbers.likeZero = function (v) {
+        return DNumbers.like(v, 0);
+    };
+    DNumbers.like2PI = function (v) {
+        return DNumbers.like(DNumbers.rad2Deg(v), DPoint_1.DOUBLE_PI_IN_DEGREE);
+    };
+    DNumbers.likePI = function (v) {
+        return DNumbers.like(DNumbers.rad2Deg(v), DPoint_1.PI_IN_DEGREE);
+    };
+    DNumbers.rad2Deg = function (v) {
+        return v * DPoint_1.DEGREE_TO_PI;
+    };
+    DNumbers.deg2Rad = function (v) {
+        return v * DPoint_1.PI_TO_DEGREE;
+    };
+    return DNumbers;
+}());
+exports.DNumbers = DNumbers;

package/dist/cjs/DPlane.d.ts

@@ -0,0 +1,25 @@
+import { DPoint } from './DPoint';
+import { DPolygon } from './DPolygon';
+export declare class DPlane {
+    a: number;
+    b: number;
+    c: number;
+    d: number;
+    p1: DPoint;
+    p2: DPoint;
+    p3: DPoint;
+    constructor(a: number, b: number, c: number, d: number, p1?: DPoint, p2?: DPoint, p3?: DPoint);
+    static find(p1: DPoint, p2: DPoint, p3: DPoint): DPlane;
+    x(p: DPoint): DPoint;
+    x(p: DPolygon): DPolygon;
+    y(p: DPoint): DPoint;
+    y(p: DPolygon): DPolygon;
+    z(p: DPoint): DPoint;
+    z(p: DPolygon): DPolygon;
+    clone(): DPlane;
+    distance(p: DPoint): number;
+    distance(p: DPlane): number;
+    equal(p: DPlane): boolean;
+    same(p: DPlane): boolean;
+    parallel(p: DPlane): boolean;
+}