node-red-contrib-tak-registration
Version:
A Node-RED node to register to TAK and to help wrap files as datapackages to send to TAK
162 lines (154 loc) • 6.06 kB
JavaScript
Object.defineProperty(exports, "__esModule", {value: true});var __defProp = Object.defineProperty;
var __defProps = Object.defineProperties;
var __getOwnPropDescs = Object.getOwnPropertyDescriptors;
var __getOwnPropSymbols = Object.getOwnPropertySymbols;
var __hasOwnProp = Object.prototype.hasOwnProperty;
var __propIsEnum = Object.prototype.propertyIsEnumerable;
var __defNormalProp = (obj, key, value) => key in obj ? __defProp(obj, key, { enumerable: true, configurable: true, writable: true, value }) : obj[key] = value;
var __spreadValues = (a, b) => {
for (var prop in b || (b = {}))
if (__hasOwnProp.call(b, prop))
__defNormalProp(a, prop, b[prop]);
if (__getOwnPropSymbols)
for (var prop of __getOwnPropSymbols(b)) {
if (__propIsEnum.call(b, prop))
__defNormalProp(a, prop, b[prop]);
}
return a;
};
var __spreadProps = (a, b) => __defProps(a, __getOwnPropDescs(b));
// index.ts
var _distance = require('@turf/distance');
var _meta = require('@turf/meta');
var _helpers = require('@turf/helpers');
var _invariant = require('@turf/invariant');
function nearestPointOnLine(lines, pt, options = {}) {
if (!lines || !pt) {
throw new Error("lines and pt are required arguments");
}
const ptPos = _invariant.getCoord.call(void 0, pt);
let closestPt = _helpers.point.call(void 0, [Infinity, Infinity], {
dist: Infinity,
index: -1,
multiFeatureIndex: -1,
location: -1
});
let length = 0;
_meta.flattenEach.call(void 0,
lines,
function(line, _featureIndex, multiFeatureIndex) {
const coords = _invariant.getCoords.call(void 0, line);
for (let i = 0; i < coords.length - 1; i++) {
const start = _helpers.point.call(void 0, coords[i]);
start.properties.dist = _distance.distance.call(void 0, pt, start, options);
const startPos = _invariant.getCoord.call(void 0, start);
const stop = _helpers.point.call(void 0, coords[i + 1]);
stop.properties.dist = _distance.distance.call(void 0, pt, stop, options);
const stopPos = _invariant.getCoord.call(void 0, stop);
const sectionLength = _distance.distance.call(void 0, start, stop, options);
let intersectPos;
let wasEnd;
if (startPos[0] === ptPos[0] && startPos[1] === ptPos[1]) {
[intersectPos, , wasEnd] = [startPos, void 0, false];
} else if (stopPos[0] === ptPos[0] && stopPos[1] === ptPos[1]) {
[intersectPos, , wasEnd] = [stopPos, void 0, true];
} else {
[intersectPos, , wasEnd] = nearestPointOnSegment(
start.geometry.coordinates,
stop.geometry.coordinates,
_invariant.getCoord.call(void 0, pt)
);
}
let intersectPt;
if (intersectPos) {
intersectPt = _helpers.point.call(void 0, intersectPos, {
dist: _distance.distance.call(void 0, pt, intersectPos, options),
multiFeatureIndex,
location: length + _distance.distance.call(void 0, start, intersectPos, options)
});
}
if (intersectPt && intersectPt.properties.dist < closestPt.properties.dist) {
closestPt = __spreadProps(__spreadValues({}, intersectPt), {
properties: __spreadProps(__spreadValues({}, intersectPt.properties), {
// Legacy behaviour where index progresses to next segment # if we
// went with the end point this iteration.
index: wasEnd ? i + 1 : i
})
});
}
length += sectionLength;
}
}
);
return closestPt;
}
function dot(v1, v2) {
const [v1x, v1y, v1z] = v1;
const [v2x, v2y, v2z] = v2;
return v1x * v2x + v1y * v2y + v1z * v2z;
}
function cross(v1, v2) {
const [v1x, v1y, v1z] = v1;
const [v2x, v2y, v2z] = v2;
return [v1y * v2z - v1z * v2y, v1z * v2x - v1x * v2z, v1x * v2y - v1y * v2x];
}
function magnitude(v) {
return Math.sqrt(Math.pow(v[0], 2) + Math.pow(v[1], 2) + Math.pow(v[2], 2));
}
function angle(v1, v2) {
const theta = dot(v1, v2) / (magnitude(v1) * magnitude(v2));
return Math.acos(Math.min(Math.max(theta, -1), 1));
}
function lngLatToVector(a) {
const lat = _helpers.degreesToRadians.call(void 0, a[1]);
const lng = _helpers.degreesToRadians.call(void 0, a[0]);
return [
Math.cos(lat) * Math.cos(lng),
Math.cos(lat) * Math.sin(lng),
Math.sin(lat)
];
}
function vectorToLngLat(v) {
const [x, y, z] = v;
const lat = _helpers.radiansToDegrees.call(void 0, Math.asin(z));
const lng = _helpers.radiansToDegrees.call(void 0, Math.atan2(y, x));
return [lng, lat];
}
function nearestPointOnSegment(posA, posB, posC) {
const A = lngLatToVector(posA);
const B = lngLatToVector(posB);
const C = lngLatToVector(posC);
const [Cx, Cy, Cz] = C;
const [D, E, F] = cross(A, B);
const a = E * Cz - F * Cy;
const b = F * Cx - D * Cz;
const c = D * Cy - E * Cx;
const f = c * E - b * F;
const g = a * F - c * D;
const h = b * D - a * E;
const t = 1 / Math.sqrt(Math.pow(f, 2) + Math.pow(g, 2) + Math.pow(h, 2));
const I1 = [f * t, g * t, h * t];
const I2 = [-1 * f * t, -1 * g * t, -1 * h * t];
const angleAB = angle(A, B);
const angleAI1 = angle(A, I1);
const angleBI1 = angle(B, I1);
const angleAI2 = angle(A, I2);
const angleBI2 = angle(B, I2);
let I;
if (angleAI1 < angleAI2 && angleAI1 < angleBI2 || angleBI1 < angleAI2 && angleBI1 < angleBI2) {
I = I1;
} else {
I = I2;
}
if (angle(A, I) > angleAB || angle(B, I) > angleAB) {
if (_distance.distance.call(void 0, vectorToLngLat(I), vectorToLngLat(A)) <= _distance.distance.call(void 0, vectorToLngLat(I), vectorToLngLat(B))) {
return [vectorToLngLat(A), true, false];
} else {
return [vectorToLngLat(B), false, true];
}
}
return [vectorToLngLat(I), false, false];
}
var turf_nearest_point_on_line_default = nearestPointOnLine;
exports.default = turf_nearest_point_on_line_default; exports.nearestPointOnLine = nearestPointOnLine;
//# sourceMappingURL=index.cjs.map
;