@webviz/subsurface-viewer 1.15.7 → 1.15.9

package/dist/layers/wells/wellsLayer.js

@@ -1,17 +1,20 @@
-import { interpolateNumberArray } from "d3";
 import _, { isEmpty, isEqual } from "lodash";
-import { distance, dot, subtract } from "mathjs";
 import { CompositeLayer } from "@deck.gl/core";
 import { PathStyleExtension } from "@deck.gl/extensions";
-import { GeoJsonLayer, PathLayer } from "@deck.gl/layers";
+import { GeoJsonLayer } from "@deck.gl/layers";
 import { GL } from "@luma.gl/constants";
-import { getColors, rgbValues } from "@emerson-eps/color-tables/";
-import { createPropertyData, getLayersById, isDrawingEnabled, } from "../utils/layerTools";
+import { createPropertyData, getFromAccessor, getLayersById, isDrawingEnabled, } from "../utils/layerTools";
 import { SectionViewport } from "../../viewports";
 import { WellLabelLayer } from "./layers/wellLabelLayer";
 import { abscissaTransform } from "./utils/abscissaTransform";
 import { GetBoundingBox, checkWells, coarsenWells, invertPath, splineRefine, } from "./utils/spline";
-import { getColor, getTrajectory } from "./utils/trajectory";
+import { getColor, getMd, getTrajectory, getTvd } from "./utils/trajectory";
+import { LogCurveLayer } from "./layers/logCurveLayer";
+import { getWellMds, getWellObjectByName, getWellObjectsByName, } from "./utils/wells";
+import { getLegendData, getLogProperty } from "./utils/log";
+import { getSize, LINE, POINT, DEFAULT_LINE_WIDTH } from "./utils/features";
+import { scaleArray } from "../../utils/arrays";
+const DEFAULT_DASH = [5, 5];
 export var SubLayerId;
 (function (SubLayerId) {
     SubLayerId["COLORS"] = "colors";
@@ -47,27 +50,6 @@ const defaultProps = {
     wellNameColor: [0, 0, 0, 255],
     wellNameSize: 10,
 };
-const LINE = "line";
-const POINT = "point";
-const DEFAULT_POINT_SIZE = 8;
-const DEFAULT_LINE_WIDTH = 5;
-const DEFAULT_DASH = [5, 5];
-export function getSize(type, accessor, offset = 0) {
-    if (typeof accessor == "function") {
-        return (object) => {
-            return (accessor(object) + offset);
-        };
-    }
-    if (accessor == 0)
-        return 0;
-    if (accessor > 0)
-        return accessor + offset;
-    if (type == LINE)
-        return DEFAULT_LINE_WIDTH + offset;
-    if (type == POINT)
-        return DEFAULT_POINT_SIZE + offset;
-    return 0;
-}
 export default class WellsLayer extends CompositeLayer {
     recomputeDataState() {
         const { data, refine, ZIncreasingDownwards, section } = this.props;
@@ -117,23 +99,28 @@ export default class WellsLayer extends CompositeLayer {
         }
     }
     onClick(info) {
-        var _a;
         // Make selection only when drawing is disabled
         if (isDrawingEnabled(this.context.layerManager)) {
             return false;
         }
         else {
             this.context.userData.setEditedData({
-                selectedWell: (_a = info.object) === null || _a === void 0 ? void 0 : _a.properties.name,
+                selectedWell: info.wellName,
             });
             return false; // do not return true to allow DeckGL props.onClick to be called
         }
     }
-    setSelection(well, _selection) {
-        if (this.internalState) {
+    setSelection(well, selection) {
+        if (!this.internalState)
+            return;
+        if (!well || !selection) {
+            this.setState({
+                logDomainSelection: undefined,
+            });
+        }
+        else {
             this.setState({
-                well: well,
-                selection: _selection,
+                logDomainSelection: { well, selection },
             });
         }
     }
@@ -151,6 +138,7 @@ export default class WellsLayer extends CompositeLayer {
         }
         return data;
     }
+    // ? Unused ?
     getLegendData(value) {
         return getLegendData(value, "", this.props.logName, this.props.logColor);
     }
@@ -160,20 +148,19 @@ export default class WellsLayer extends CompositeLayer {
         });
     }
     getLogLayer() {
-        var _a;
+        var _a, _b;
         const sub_layers = (_a = this.internalState) === null || _a === void 0 ? void 0 : _a.subLayers;
-        const log_layer = getLayersById(sub_layers, "wells-layer-log_curve");
-        return log_layer === null || log_layer === void 0 ? void 0 : log_layer[0];
+        const log_layers = getLayersById(sub_layers, "wells-layer-log_curve");
+        return (_b = log_layers === null || log_layers === void 0 ? void 0 : log_layers[0]) === null || _b === void 0 ? void 0 : _b.getCurveLayer();
     }
     getSelectionLayer() {
-        var _a;
+        var _a, _b;
         const sub_layers = (_a = this.internalState) === null || _a === void 0 ? void 0 : _a.subLayers;
-        const log_layer = getLayersById(sub_layers, "wells-layer-selection");
-        return log_layer === null || log_layer === void 0 ? void 0 : log_layer[0];
+        const log_layer = getLayersById(sub_layers, "wells-layer-log_curve");
+        return (_b = log_layer === null || log_layer === void 0 ? void 0 : log_layer[0]) === null || _b === void 0 ? void 0 : _b.getSelectionLayer();
     }
     getLogCurveData() {
-        const log_layer = this.getLogLayer();
-        return log_layer === null || log_layer === void 0 ? void 0 : log_layer.props.data;
+        return this.state.logData;
     }
     setupLegend() {
         const data = this.getLogCurveData();
@@ -264,92 +251,45 @@ export default class WellsLayer extends CompositeLayer {
         const highlightLayer = new GeoJsonLayer(highlightLayerProps);
         // Highlight the multi selected wells.
         const highlightMultiWellsLayer = new GeoJsonLayer(highlightMultiWellsLayerProps);
-        const logLayer = new PathLayer(this.getSubLayerProps({
+        const logLayer = new LogCurveLayer(Object.assign({}, this.getSubLayerProps({
             id: SubLayerId.LOG_CURVE,
-            data: this.props.logData,
-            positionFormat,
+            visible: this.props.logCurves && !fastDrawing,
             pickable: true,
-            widthScale: 10,
-            widthMinPixels: 1,
-            miterLimit: 100,
-            getPath: (d) => {
-                var _a;
-                return getLogPath(data.features, d, this.props.logrunName, (_a = this.props.lineStyle) === null || _a === void 0 ? void 0 : _a.color);
-            },
-            getColor: (d) => getLogColor(d, this.props.logrunName, this.props.logName, this.props.logColor, this.context.userData
-                .colorTables, this.props.colorMappingFunction, this.props.isLog),
-            getWidth: (d) => this.props.logRadius ||
-                getLogWidth(d, this.props.logrunName, this.props.logName),
-            updateTriggers: {
-                getColor: [
-                    this.props.logrunName,
-                    this.props.logName,
-                    this.props.logColor,
-                    this.context.userData
-                        .colorTables,
-                    this.props.isLog,
-                ],
-                getWidth: [
-                    this.props.logrunName,
-                    this.props.logName,
-                    this.props.logRadius,
-                ],
-                getPath: [positionFormat],
-            },
-            onDataLoad: (value) => {
-                this.setLegend(value);
-            },
             parameters,
-            visible: this.props.logCurves && !fastDrawing,
-        }));
-        const selectionLayer = new PathLayer(this.getSubLayerProps({
-            id: SubLayerId.SELECTION,
-            data: this.props.logData,
             positionFormat,
-            pickable: false,
-            widthScale: 10,
-            widthMinPixels: 1,
-            miterLimit: 100,
-            getPath: (d) => {
-                var _a;
-                return getLogPath1(data.features, d, this.state.well, this.state.selection, this.props.logrunName, (_a = this.props.lineStyle) === null || _a === void 0 ? void 0 : _a.color);
+            data: this.props.logData,
+            trajectoryData: data,
+            domainSelection: this.state.logDomainSelection,
+            logWidth: this.props.logRadius || "value",
+            activeLog: {
+                logRun: this.props.logrunName,
+                curveName: this.props.logName,
             },
-            getColor: (d) => getLogColor1(data.features, d, this.state.well, this.state.selection, this.props.logrunName),
-            getWidth: (d) => this.props.logRadius * 1.5 ||
-                getLogWidth(d, this.props.logrunName, this.props.logName),
-            updateTriggers: {
-                getColor: [
-                    this.props.logrunName,
-                    this.state.well,
-                    this.state.selection,
-                ],
-                getWidth: [
-                    this.props.logrunName,
-                    this.props.logName,
-                    this.props.logRadius,
-                ],
-                getPath: [
-                    positionFormat,
-                    this.props.logrunName,
-                    this.state.well,
-                    this.state.selection,
-                ],
+            colorScale: {
+                name: this.props.logColor,
+                isLogarithmic: this.props.isLog,
+                colorTables: this.context.userData
+                    .colorTables,
+                mappingFunction: this.props.colorMappingFunction,
             },
-            onDataLoad: (value) => {
+            getWellMds: getWellMds,
+            getTrajectoryPath: (d) => { var _a; return getTrajectory(d, (_a = this.props.lineStyle) === null || _a === void 0 ? void 0 : _a.color); },
+            onLogDataComputed: (value) => {
+                this.setState({ logData: value });
                 this.setLegend(value);
             },
-            parameters,
-            visible: this.props.logCurves && !fastDrawing,
-        }));
+            updateTriggers: {
+                getTrajectoryPath: [this.props.lineStyle],
+            },
+        })));
         const namesLayer = this.createWellLabelLayer(data.features);
         const layers = [
             outlineLayer,
-            logLayer,
             colorsLayer,
             highlightLayer,
             highlightMultiWellsLayer,
-            selectionLayer,
             namesLayer,
+            logLayer,
         ];
         if (fastDrawing) {
             layers.push(fastLayer);
@@ -357,7 +297,7 @@ export default class WellsLayer extends CompositeLayer {
         return layers;
     }
     getPickingInfo({ info }) {
-        var _a, _b, _c, _d, _e, _f;
+        var _a, _b, _c, _d, _e, _f, _g, _h, _j, _k;
         const noLog = {
             properties: [],
             logName: "",
@@ -371,6 +311,7 @@ export default class WellsLayer extends CompositeLayer {
         if (typeof coordinate[2] !== "undefined") {
             coordinate[2] /= Math.max(0.001, zScale);
         }
+        let wellName = undefined;
         let md_property = null;
         let tvd_property = null;
         let log_property = null;
@@ -378,6 +319,7 @@ export default class WellsLayer extends CompositeLayer {
         if (((_d = info.sourceLayer) === null || _d === void 0 ? void 0 : _d.id) === `${this.props.id}-${SubLayerId.LOG_CURVE}`) {
             // The user is hovering a well log entry
             const logPick = info;
+            wellName = (_f = (_e = logPick.object) === null || _e === void 0 ? void 0 : _e.header) === null || _f === void 0 ? void 0 : _f.well;
             md_property = getLogProperty(coordinate, features, logPick.object, this.props.logrunName, "MD");
             tvd_property = getLogProperty(coordinate, features, logPick.object, this.props.logrunName, "TVD");
             log_property = getLogProperty(coordinate, features, info.object, this.props.logrunName, this.props.logName);
@@ -385,8 +327,9 @@ export default class WellsLayer extends CompositeLayer {
         else {
             // User is hovering a wellbore path
             const wellpickInfo = info;
-            md_property = getMdProperty(coordinate, wellpickInfo.object, (_e = this.props.lineStyle) === null || _e === void 0 ? void 0 : _e.color, info.featureType);
-            tvd_property = getTvdProperty(coordinate, info.object, (_f = this.props.lineStyle) === null || _f === void 0 ? void 0 : _f.color, info.featureType);
+            wellName = (_h = (_g = wellpickInfo.object) === null || _g === void 0 ? void 0 : _g.properties) === null || _h === void 0 ? void 0 : _h.name;
+            md_property = getMdProperty(coordinate, wellpickInfo.object, (_j = this.props.lineStyle) === null || _j === void 0 ? void 0 : _j.color, info.featureType);
+            tvd_property = getTvdProperty(coordinate, info.object, (_k = this.props.lineStyle) === null || _k === void 0 ? void 0 : _k.color, info.featureType);
         }
         // Patch for inverting tvd readout to fix issue #830,
         // should make proper fix when handling z increase direction - issue #842
@@ -398,7 +341,7 @@ export default class WellsLayer extends CompositeLayer {
             layer_properties.push(inverted_tvd_property);
         if (log_property)
             layer_properties.push(log_property);
-        return Object.assign(Object.assign({}, info), { properties: layer_properties, logName: (log_property === null || log_property === void 0 ? void 0 : log_property.name) || "" });
+        return Object.assign(Object.assign({}, info), { properties: layer_properties, logName: log_property === null || log_property === void 0 ? void 0 : log_property.name, wellName: wellName });
     }
 }
 WellsLayer.layerName = "WellsLayer";
@@ -410,392 +353,27 @@ function onDataLoad(data, context) {
         context.layer.props.reportBoundingBox({ layerBoundingBox: bbox });
     }
 }
-function multiply(pair, factor) {
-    return [pair[0] * factor, pair[1] * factor];
-}
-function getDashFactor(accessor, width_accessor, offset = 0) {
-    return (object, objectInfo) => {
-        let width = DEFAULT_LINE_WIDTH;
-        if (typeof width_accessor == "function") {
-            width = width_accessor(object);
-        }
-        else if (width_accessor) {
-            width = width_accessor;
-        }
+// Create a dash accessor function that scales
+function getDashFactor(accessor, widthAccessor, offset = 0) {
+    return (object, ctx) => {
+        var _a;
+        const dashValue = getFromAccessor(accessor, object, ctx);
+        const width = (_a = getFromAccessor(widthAccessor, object, ctx)) !== null && _a !== void 0 ? _a : DEFAULT_LINE_WIDTH;
         const factor = width / (width + offset);
-        let dash = [0, 0];
-        if (typeof accessor == "function") {
-            dash = accessor(object, objectInfo);
+        let dashArray = [0, 0];
+        if (dashValue && Array.isArray(dashValue)) {
+            dashArray = dashValue;
         }
-        else if (accessor)
-            dash = accessor;
-        else if (accessor)
-            dash = DEFAULT_DASH;
-        if (dash.length == 2) {
-            return multiply(dash, factor);
-        }
-        else {
-            return multiply(DEFAULT_DASH, factor);
+        else if (dashValue) {
+            dashArray = DEFAULT_DASH;
         }
+        return scaleArray(dashArray, factor);
     };
 }
 function dataIsReady(layerData) {
     // Deck.gl always shows prop.data as `[]` while external data is being loaded
     return !!layerData && !isEmpty(layerData);
 }
-function getColumn(data, col) {
-    const column = [];
-    for (let i = 0; i < data.length; i++) {
-        column.push(data[i][col]);
-    }
-    return column;
-}
-function getLogMd(d, logrun_name) {
-    if (!isSelectedLogRun(d, logrun_name))
-        return [];
-    const names_md = ["DEPTH", "DEPT", "MD", "TDEP", "MD_RKB"]; // aliases for MD
-    const log_id = getLogIndexByNames(d, names_md);
-    return log_id >= 0 ? getColumn(d.data, log_id) : [];
-}
-export function getLogValues(d, logrun_name, log_name) {
-    if (!isSelectedLogRun(d, logrun_name))
-        return [];
-    const log_id = getLogIndexByName(d, log_name);
-    return log_id >= 0 ? getColumn(d.data, log_id) : [];
-}
-export function getLogInfo(d, logrun_name, log_name) {
-    if (!isSelectedLogRun(d, logrun_name))
-        return undefined;
-    const log_id = getLogIndexByName(d, log_name);
-    return d.curves[log_id];
-}
-function getDiscreteLogMetadata(d, log_name) {
-    return d === null || d === void 0 ? void 0 : d.metadata_discrete[log_name];
-}
-function isSelectedLogRun(d, logrun_name) {
-    return d.header.name.toLowerCase() === logrun_name.toLowerCase();
-}
-function getWellObjectByName(wells_data, name) {
-    return wells_data === null || wells_data === void 0 ? void 0 : wells_data.find((item) => { var _a, _b; return ((_b = (_a = item.properties) === null || _a === void 0 ? void 0 : _a.name) === null || _b === void 0 ? void 0 : _b.toLowerCase()) === (name === null || name === void 0 ? void 0 : name.toLowerCase()); });
-}
-function getWellObjectsByName(wells_data, name) {
-    const res = [];
-    for (let i = 0; i < (name === null || name === void 0 ? void 0 : name.length); i++) {
-        wells_data === null || wells_data === void 0 ? void 0 : wells_data.find((item) => {
-            var _a, _b, _c;
-            if (((_b = (_a = item.properties) === null || _a === void 0 ? void 0 : _a.name) === null || _b === void 0 ? void 0 : _b.toLowerCase()) === ((_c = name[i]) === null || _c === void 0 ? void 0 : _c.toLowerCase()))
-                res.push(item);
-        });
-    }
-    return res;
-}
-function getPointGeometry(well_object) {
-    const geometries = well_object.geometry.geometries;
-    return geometries.find((item) => item.type === "Point");
-}
-// Return well head position from Point Geometry
-function getWellHeadPosition(well_object) {
-    var _a, _b;
-    return (_b = (_a = getPointGeometry(well_object)) === null || _a === void 0 ? void 0 : _a.coordinates) !== null && _b !== void 0 ? _b : [-1, -1, -1];
-}
-function getWellMds(well_object) {
-    return well_object.properties.md[0];
-}
-function getNeighboringMdIndices(mds, md) {
-    const idx = mds.findIndex((x) => x >= md);
-    return idx === 0 ? [idx, idx + 1] : [idx - 1, idx];
-}
-function getPositionByMD(well_xyz, well_mds, md) {
-    const [l_idx, h_idx] = getNeighboringMdIndices(well_mds, md);
-    const md_low = well_mds[l_idx];
-    const md_normalized = (md - md_low) / (well_mds[h_idx] - md_low);
-    return interpolateNumberArray(well_xyz[l_idx], well_xyz[h_idx])(md_normalized);
-}
-function getLogPath(wells_data, d, logrun_name, trajectory_line_color) {
-    const well_object = getWellObjectByName(wells_data, d.header.well);
-    if (!well_object)
-        return [];
-    const well_xyz = getTrajectory(well_object, trajectory_line_color);
-    const well_mds = getWellMds(well_object);
-    if (well_xyz == undefined ||
-        well_mds == undefined ||
-        well_xyz.length == 0 ||
-        well_mds.length == 0)
-        return [];
-    const log_xyz = [];
-    const log_mds = getLogMd(d, logrun_name);
-    log_mds.forEach((md) => {
-        const xyz = getPositionByMD(well_xyz, well_mds, md);
-        log_xyz.push(xyz);
-    });
-    return log_xyz;
-}
-function getLogIndexByName(d, log_name) {
-    const name = log_name.toLowerCase();
-    return d.curves.findIndex((item) => item.name.toLowerCase() === name);
-}
-function getLogIndexByNames(d, names) {
-    for (const name of names) {
-        const index = getLogIndexByName(d, name);
-        if (index >= 0)
-            return index;
-    }
-    return -1;
-}
-function getLogColor(d, logrun_name, log_name, logColor, colorTables, colorMappingFunction, isLog) {
-    var _a;
-    const log_data = getLogValues(d, logrun_name, log_name);
-    const log_info = getLogInfo(d, logrun_name, log_name);
-    if (log_data.length == 0 || log_info == undefined)
-        return [];
-    const log_color = [];
-    if (log_info.description == "continuous") {
-        const min = Math.min(...log_data);
-        const max = Math.max(...log_data);
-        const max_delta = max - min;
-        log_data.forEach((value) => {
-            const adjustedVal = (value - min) / max_delta;
-            const rgb = colorMappingFunction
-                ? colorMappingFunction(adjustedVal)
-                : rgbValues(adjustedVal, logColor, colorTables, isLog);
-            if (rgb) {
-                log_color.push([rgb[0], rgb[1], rgb[2]]);
-            }
-            else {
-                log_color.push([0, 0, 0, 0]); // push transparent for null/undefined log values
-            }
-        });
-    }
-    else {
-        // well log data set for ex : H1: Array(2)0: (4) [255, 26, 202, 255] 1: 13
-        const log_attributes = (_a = getDiscreteLogMetadata(d, log_name)) === null || _a === void 0 ? void 0 : _a.objects;
-        const logLength = Object.keys(log_attributes).length;
-        const attributesObject = {};
-        const categorical = true;
-        Object.keys(log_attributes).forEach((key) => {
-            // get the point from log_attributes
-            const point = log_attributes[key][1];
-            const categoricalMin = 0;
-            const categoricalMax = logLength - 1;
-            let rgb;
-            if (colorMappingFunction) {
-                rgb = colorMappingFunction(point, categorical, categoricalMin, categoricalMax);
-            }
-            else {
-                // if color-map function is not defined
-                const arrayOfColors = getColors(logColor, colorTables, point);
-                if (!arrayOfColors.length)
-                    console.error(`Empty or missed '${logColor}' color table`);
-                else {
-                    rgb = arrayOfColors;
-                }
-            }
-            if (rgb) {
-                if (rgb.length === 3) {
-                    attributesObject[key] = [[rgb[0], rgb[1], rgb[2]], point];
-                }
-                else {
-                    // ? What is the point of this? Why do we offset the index in this case, isn't the fourth value the opacity?
-                    // (@anders2303)
-                    attributesObject[key] = [[rgb[1], rgb[2], rgb[3]], point];
-                }
-            }
-        });
-        log_data.forEach((log_value) => {
-            var _a;
-            const dl_attrs = (_a = Object.entries(attributesObject).find(([, value]) => value[1] == log_value)) === null || _a === void 0 ? void 0 : _a[1];
-            if (dl_attrs)
-                log_color.push(dl_attrs[0]);
-            else
-                log_color.push([0, 0, 0, 0]); // use transparent for undefined/null log values
-        });
-    }
-    return log_color;
-}
-function getLogPath1(wells_data, d, selectedWell, selection, logrun_name, trajectory_line_color) {
-    if (!selection || selectedWell !== d.header.well)
-        return [];
-    const well_object = getWellObjectByName(wells_data, d.header.well);
-    if (!well_object)
-        return [];
-    const well_xyz = getTrajectory(well_object, trajectory_line_color);
-    const well_mds = getWellMds(well_object);
-    if (well_xyz == undefined ||
-        well_mds == undefined ||
-        well_xyz.length == 0 ||
-        well_mds.length == 0)
-        return [];
-    const log_mds = getLogMd(d, logrun_name);
-    if (!log_mds)
-        return [];
-    const log_xyz = [];
-    let md0 = selection[0];
-    if (md0 !== undefined) {
-        let md1 = selection[1];
-        if (md1 == md0)
-            md1 = undefined;
-        const mdFirst = well_mds[0];
-        const mdLast = well_mds[well_mds.length - 1];
-        if (md1 !== undefined) {
-            if (md0 > md1) {
-                const tmp = md0;
-                md0 = md1;
-                md1 = tmp;
-            }
-        }
-        const delta = 2;
-        if (md0 - delta > mdFirst) {
-            let xyz = getPositionByMD(well_xyz, well_mds, md0 - delta);
-            log_xyz.push(xyz);
-            xyz = getPositionByMD(well_xyz, well_mds, md0);
-            log_xyz.push(xyz);
-        }
-        if (md1 !== undefined) {
-            const _md1 = md1;
-            let index = 0;
-            well_mds.forEach((md) => {
-                if (md0 <= md && md <= _md1) {
-                    const xyz = well_xyz[index];
-                    log_xyz.push(xyz);
-                }
-                index++;
-            });
-            if (_md1 + delta < mdLast) {
-                let xyz = getPositionByMD(well_xyz, well_mds, _md1);
-                log_xyz.push(xyz);
-                xyz = getPositionByMD(well_xyz, well_mds, _md1 + delta);
-                log_xyz.push(xyz);
-            }
-        }
-    }
-    return log_xyz;
-}
-function getLogColor1(wells_data, d, selectedWell, selection, logrun_name) {
-    if (!selection || selectedWell !== d.header.well)
-        return [];
-    const well_object = getWellObjectByName(wells_data, d.header.well);
-    if (!well_object)
-        return [];
-    const well_mds = getWellMds(well_object);
-    const log_mds = getLogMd(d, logrun_name);
-    if (!log_mds || log_mds.length === 0)
-        return [];
-    const log_color = [];
-    let md0 = selection[0];
-    if (md0 !== undefined) {
-        const mdFirst = well_mds[0];
-        const mdLast = well_mds[well_mds.length - 1];
-        let md1 = selection[1];
-        if (md1 == md0)
-            md1 = undefined;
-        let swap = false;
-        if (md1 !== undefined) {
-            if (md0 > md1) {
-                const tmp = md0;
-                md0 = md1;
-                md1 = tmp;
-                swap = true;
-            }
-        }
-        const delta = 2;
-        if (md0 - delta > mdFirst)
-            log_color.push(swap ? [0, 255, 0, 128] : [255, 0, 0, 128]);
-        if (md1 !== undefined) {
-            const _md1 = md1;
-            log_color.push([128, 128, 128, 128]);
-            well_mds.forEach((md) => {
-                if (md0 <= md && md <= _md1) {
-                    log_color.push([128, 128, 128, 128]);
-                }
-            });
-            if (_md1 + delta < mdLast)
-                log_color.push(swap ? [255, 0, 0, 128] : [0, 255, 0, 128]);
-        }
-    }
-    return log_color;
-}
-function getLogWidth(d, logrun_name, log_name) {
-    return getLogValues(d, logrun_name, log_name);
-}
-function squared_distance(a, b) {
-    const dx = a[0] - b[0];
-    const dy = a[1] - b[1];
-    return dx * dx + dy * dy;
-}
-// Return distance between line segment vw and point p
-function distToSegmentSquared(v, w, p) {
-    const l2 = squared_distance(v, w);
-    if (l2 == 0)
-        return squared_distance(p, v);
-    let t = ((p[0] - v[0]) * (w[0] - v[0]) + (p[1] - v[1]) * (w[1] - v[1])) / l2;
-    t = Math.max(0, Math.min(1, t));
-    return squared_distance(p, [
-        v[0] + t * (w[0] - v[0]),
-        v[1] + t * (w[1] - v[1]),
-    ]);
-}
-function isPointAwayFromLineEnd(point, line) {
-    const ab = subtract(line[1], line[0]);
-    const cb = subtract(line[1], point);
-    const dotProduct = dot(ab, cb);
-    // If the dot product is negative, the point has moved past the end of the line
-    return dotProduct < 0;
-}
-// Interpolates point closest to the coords on trajectory
-function interpolateDataOnTrajectory(coord, data, trajectory) {
-    // if number of data points in less than 1 or
-    // length of data and trajectory are different we cannot interpolate.
-    if (data.length <= 1 || data.length != trajectory.length)
-        return -1;
-    // Identify closest well path leg to coord.
-    const segment_index = getSegmentIndex(coord, trajectory);
-    const index0 = segment_index;
-    const index1 = index0 + 1;
-    // Get the nearest data.
-    const data0 = data[index0];
-    const data1 = data[index1];
-    // Get the nearest survey points.
-    const survey0 = trajectory[index0];
-    const survey1 = trajectory[index1];
-    // To avoid interpolating longer than the actual wellbore path we ignore the coordinate if it's moved beyond the last line
-    if (index1 === trajectory.length - 1 &&
-        isPointAwayFromLineEnd(coord, [survey0, survey1])) {
-        coord = survey1;
-    }
-    const dv = distance(survey0, survey1);
-    if (dv === 0) {
-        return -1;
-    }
-    // Calculate the scalar projection onto segment.
-    const v0 = subtract(coord, survey0);
-    const v1 = subtract(survey1, survey0);
-    // scalar_projection in interval [0,1]
-    const scalar_projection = dot(v0, v1) / (dv * dv);
-    // Interpolate data.
-    return data0 * (1.0 - scalar_projection) + data1 * scalar_projection;
-}
-function getMd(coord, feature, accessor) {
-    var _a, _b;
-    if (!((_b = (_a = feature.properties) === null || _a === void 0 ? void 0 : _a["md"]) === null || _b === void 0 ? void 0 : _b[0]) || !feature.geometry)
-        return null;
-    const measured_depths = feature.properties.md[0];
-    const trajectory3D = getTrajectory(feature, accessor);
-    if (trajectory3D == undefined)
-        return null;
-    let trajectory;
-    // In 2D view coord is of type Point2D and in 3D view it is Point3D,
-    // so use appropriate trajectory for interpolation
-    if (coord.length == 2) {
-        const trajectory2D = trajectory3D.map((v) => {
-            return v.slice(0, 2);
-        });
-        trajectory = trajectory2D;
-    }
-    else {
-        trajectory = trajectory3D;
-    }
-    return interpolateDataOnTrajectory(coord, measured_depths, trajectory);
-}
 function getMdProperty(coord, feature, accessor, featureType) {
     var _a, _b;
     if (featureType === "points") {
@@ -808,31 +386,6 @@ function getMdProperty(coord, feature, accessor, featureType) {
     }
     return null;
 }
-function getTvd(coord, feature, accessor) {
-    var _a;
-    const trajectory3D = getTrajectory(feature, accessor);
-    // if trajectory is not found or if it has a data single point then get tvd from well head
-    if (trajectory3D == undefined || (trajectory3D === null || trajectory3D === void 0 ? void 0 : trajectory3D.length) <= 1) {
-        const wellhead_xyz = getWellHeadPosition(feature);
-        return (_a = wellhead_xyz === null || wellhead_xyz === void 0 ? void 0 : wellhead_xyz[2]) !== null && _a !== void 0 ? _a : null;
-    }
-    let trajectory;
-    // For 2D view coord is Point2D and for 3D view it is Point3D
-    if (coord.length == 2) {
-        const trajectory2D = trajectory3D === null || trajectory3D === void 0 ? void 0 : trajectory3D.map((v) => {
-            return v.slice(0, 2);
-        });
-        trajectory = trajectory2D;
-    }
-    else {
-        trajectory = trajectory3D;
-    }
-    const tvds = trajectory3D.map((v) => {
-        return v[2];
-    });
-    // TVD goes downwards, so it's reversed
-    return interpolateDataOnTrajectory(coord, tvds, trajectory);
-}
 function getTvdProperty(coord, feature, accessor, featureType) {
     var _a, _b;
     if (featureType === "points") {
@@ -845,80 +398,4 @@ function getTvdProperty(coord, feature, accessor, featureType) {
     }
     return null;
 }
-// Identify closest path leg to coord.
-function getSegmentIndex(coord, path) {
-    let min_d = Number.MAX_VALUE;
-    let segment_index = 0;
-    for (let i = 0; i < (path === null || path === void 0 ? void 0 : path.length) - 1; i++) {
-        const d = distToSegmentSquared(path[i], path[i + 1], coord);
-        if (d > min_d)
-            continue;
-        segment_index = i;
-        min_d = d;
-    }
-    return segment_index;
-}
-// Returns segment index of discrete logs
-function getLogSegmentIndex(coord, wells_data, log_data, logrun_name) {
-    const trajectory = getLogPath(wells_data, log_data, logrun_name);
-    return getSegmentIndex(coord, trajectory);
-}
-function getLogProperty(coord, wells_data, log_data, logrun_name, log_name) {
-    var _a, _b, _c;
-    if (!log_data.data)
-        return null;
-    const segment_index = getLogSegmentIndex(coord, wells_data, log_data, logrun_name);
-    let log_value = getLogValues(log_data, logrun_name, log_name)[segment_index];
-    let dl_attrs = undefined;
-    const dl_metadata = (_a = getDiscreteLogMetadata(log_data, log_name)) === null || _a === void 0 ? void 0 : _a.objects;
-    if (dl_metadata) {
-        dl_attrs = Object.entries(dl_metadata).find(([, value]) => value[1] == log_value);
-    }
-    const log = (_b = getLogInfo(log_data, logrun_name, log_name)) === null || _b === void 0 ? void 0 : _b.name;
-    const prop_name = log + " " + log_data.header.well;
-    log_value = dl_attrs ? dl_attrs[0] + " (" + log_value + ")" : log_value;
-    if (log_value) {
-        const well_object = getWellObjectByName(wells_data, log_data.header.well);
-        return createPropertyData(prop_name, log_value, (_c = well_object === null || well_object === void 0 ? void 0 : well_object.properties) === null || _c === void 0 ? void 0 : _c["color"]);
-    }
-    else
-        return null;
-}
-// Return data required to build well layer legend
-function getLegendData(logs, wellName, logName, logColor) {
-    if (!logs)
-        return null;
-    const log = wellName
-        ? logs.find((log) => log.header.well == wellName)
-        : logs[0];
-    const logInfo = !log
-        ? undefined
-        : getLogInfo(log, log.header.name, logName);
-    const title = "Wells / " + logName;
-    if (log && (logInfo === null || logInfo === void 0 ? void 0 : logInfo.description) == "discrete") {
-        const meta = log["metadata_discrete"];
-        const metadataDiscrete = meta[logName].objects;
-        return {
-            title: title,
-            colorName: logColor,
-            discrete: true,
-            metadata: metadataDiscrete,
-        };
-    }
-    else {
-        const minArray = [];
-        const maxArray = [];
-        logs.forEach(function (log) {
-            const logValues = getLogValues(log, log.header.name, logName);
-            minArray.push(Math.min(...logValues));
-            maxArray.push(Math.max(...logValues));
-        });
-        return {
-            title: title,
-            colorName: logColor,
-            discrete: false,
-            valueRange: [Math.min(...minArray), Math.max(...maxArray)],
-        };
-    }
-}
 //# sourceMappingURL=wellsLayer.js.map