@webviz/subsurface-viewer 0.0.1-alpha.1 → 0.0.2-alpha.2

package/src/layers/map/webworker.ts

@@ -1,479 +0,0 @@
-import { MeshType, MeshTypeLines } from "./privateMapLayer";
-import { Params } from "./mapLayer";
-
-type Vec = [number, number, number];
-
-export function makeFullMesh(e: { data: Params }): void {
-    const params = e.data;
-
-    // Keep
-    //const t0 = performance.now();
-
-    const meshData = params.meshData;
-    const propertiesData = params.propertiesData;
-    const isMesh = params.isMesh;
-    const frame = params.frame;
-    const smoothShading = params.smoothShading;
-
-    function getFloat32ArrayMinMax(data: Float32Array) {
-        let max = -99999999;
-        let min = 99999999;
-        for (let i = 0; i < data.length; i++) {
-            max = data[i] > max ? data[i] : max;
-            min = data[i] < min ? data[i] : min;
-        }
-        return [min, max];
-    }
-
-    function crossProduct(a: Vec, b: Vec): Vec {
-        const c = [
-            a[1] * b[2] - a[2] * b[1],
-            a[2] * b[0] - a[0] * b[2],
-            a[0] * b[1] - a[1] * b[0],
-        ];
-        return c as Vec;
-    }
-
-    function normalize(a: Vec): void {
-        const L = Math.sqrt(a[0] * a[0] + a[1] * a[1] + a[2] * a[2]);
-        a[0] /= L;
-        a[1] /= L;
-        a[2] /= L;
-    }
-
-    function calcNormal(
-        w: number,
-        h: number,
-        nx: number,
-        ny: number,
-        isMesh: boolean,
-        smoothShading: boolean,
-        meshData: Float32Array,
-        ox: number,
-        oy: number
-    ) {
-        if (!smoothShading) {
-            return [1, 1, 1];
-        }
-
-        if (!isMesh) {
-            return [0, 0, 1];
-        }
-
-        const i0 = h * nx + w;
-        const i1 = h * nx + (w - 1);
-        const i2 = (h + 1) * nx + w;
-        const i3 = h * nx + (w + 1);
-        const i4 = (h - 1) * nx + w;
-
-        const i0_act =                 !isNaN(meshData[i0]); // eslint-disable-line
-        const i1_act = (w - 1) >= 0 && !isNaN(meshData[i1]); // eslint-disable-line
-        const i2_act = (h + 1) < ny && !isNaN(meshData[i2]); // eslint-disable-line
-        const i3_act = (w + 1) < nx && !isNaN(meshData[i3]); // eslint-disable-line
-        const i4_act = (h - 1) >= 0 && !isNaN(meshData[i4]); // eslint-disable-line
-
-        const noNormal = [0, 0, 1]; // signals a normal could not be calculated.
-        if (!i0_act) {
-            return noNormal;
-        }
-
-        const hh = ny - 1 - h; // Note use hh for h for getting y values.
-        const p0 = [ox + w * dx,         oy + hh * dy,        i0_act ? -meshData[i0] : 0]; // eslint-disable-line
-        const p1 = [ ox + (w - 1) * dx,  oy + hh * dy,        i1_act ? -meshData[i1] : 0]; // eslint-disable-line
-        const p2 = [ ox + w * dx,        oy + (hh + 1) * dy,  i2_act ? -meshData[i2] : 0]; // eslint-disable-line
-        const p3 = [ ox + (w + 1) * dx,  oy + hh * dy,        i3_act ? -meshData[i3] : 0]; // eslint-disable-line
-        const p4 = [ ox + w * dx,        oy + (hh - 1) * dy,  i4_act ? -meshData[i4] : 0]; // eslint-disable-line
-
-        const v1 = [p1[0] - p0[0], p1[1] - p0[1], p1[2] - p0[2]] as Vec;
-        const v2 = [p2[0] - p0[0], p2[1] - p0[1], p2[2] - p0[2]] as Vec;
-        const v3 = [p3[0] - p0[0], p3[1] - p0[1], p3[2] - p0[2]] as Vec;
-        const v4 = [p4[0] - p0[0], p4[1] - p0[1], p4[2] - p0[2]] as Vec;
-
-        // Estimating a normal vector at p0:
-        // Take cross product of vectors v1, v2,
-        // Do this for all 4 quadrants.
-        // The resulting normal will be the mean of these four normals.
-        //        p2
-        //         |
-        //   p1 - p0 - p3
-        //         |
-        //        p4
-
-        const normals: Vec[] = [];
-        if (i1_act && i2_act) {
-            const normal = crossProduct(v2, v1);
-            normals.push(normal);
-        }
-
-        if (i2_act && i3_act) {
-            const normal = crossProduct(v3, v2);
-            normals.push(normal);
-        }
-
-        if (i3_act && i4_act) {
-            const normal = crossProduct(v4, v3);
-            normals.push(normal);
-        }
-
-        if (i4_act && i1_act) {
-            const normal = crossProduct(v1, v4);
-            normals.push(normal);
-        }
-
-        if (normals.length === 0) {
-            return noNormal;
-        }
-
-        const mean = normals[0];
-        for (let i = 1; i < normals.length; i++) {
-            mean[0] += normals[i][0];
-            mean[1] += normals[i][1];
-            mean[2] += normals[i][2];
-        }
-
-        normalize(mean);
-        return mean;
-    }
-
-    const meshZValueRange = getFloat32ArrayMinMax(meshData);
-    const propertyValueRange = getFloat32ArrayMinMax(propertiesData);
-
-    // Dimensions.
-    const ox = frame.origin[0];
-    const oy = frame.origin[1];
-
-    const dx = frame.increment[0];
-    const dy = frame.increment[1];
-
-    const nx = frame.count[0];
-    const ny = frame.count[1];
-
-    const propLength = propertiesData.length;
-    const isCellCenteredProperties = propLength === (nx - 1) * (ny - 1);
-
-    if (propLength !== (nx - 1) * (ny - 1) && propLength !== nx * ny) {
-        console.error(
-            "There should be as many property values as nodes (nx*ny) OR as many as cells (nx - 1) * (ny - 1)."
-        );
-    }
-
-    const positions: number[] = [];
-    const normals: number[] = [];
-    const indices: number[] = [];
-    const vertexProperties: number[] = [];
-    const vertexIndexs: number[] = [];
-    const line_positions: number[] = [];
-
-    // Note: Assumed layout of the incomming 2D array of data:
-    // First coloumn corresponds to lowest x value. Last column highest x value.
-    // First row corresponds to max y value. Last row corresponds to lowest y value.
-    // This must be taken into account when calculating vertex x,y values and texture coordinates.
-
-    if (!isCellCenteredProperties) {
-        // PROPERTIES IS SET INTERPOLATED OVER A CELL.
-        let i = 0;
-        // Loop over nodes.
-        for (let h = 0; h < ny; h++) {
-            for (let w = 0; w < nx; w++) {
-                const i0 = h * nx + w;
-
-                const x0 = ox + w * dx;
-                const y0 = oy + (ny - 1 - h) * dy; // See note above.
-                const z = isMesh ? -meshData[i0] : 0;
-
-                const propertyValue = propertiesData[i0];
-
-                positions.push(x0, y0, z);
-
-                const normal = calcNormal(w, h, nx, ny, isMesh, smoothShading, meshData, ox, oy); // eslint-disable-line
-                normals.push(normal[0], normal[1], normal[2]);
-
-                vertexProperties.push(propertyValue);
-                vertexIndexs.push(i++);
-            }
-        }
-
-        for (let h = 0; h < ny - 1; h++) {
-            for (let w = 0; w < nx - 1; w++) {
-                const i0 = h * nx + w;
-                const i1 = h * nx + (w + 1);
-                const i2 = (h + 1) * nx + (w + 1);
-                const i3 = (h + 1) * nx + w;
-
-                const i0_act = !isMesh || (!isNaN(meshData[i0]) && !isNaN(propertiesData[i0])); // eslint-disable-line
-                const i1_act = !isMesh || (!isNaN(meshData[i1]) && !isNaN(propertiesData[i1])); // eslint-disable-line
-                const i2_act = !isMesh || (!isNaN(meshData[i2]) && !isNaN(propertiesData[i2])); // eslint-disable-line
-                const i3_act = !isMesh || (!isNaN(meshData[i3]) && !isNaN(propertiesData[i3])); // eslint-disable-line
-
-                const hh = ny - h - 1; // See note above.
-
-                const x0 = ox + w * dx;
-                const y0 = oy + hh * dy;
-                const z0 = isMesh ? -meshData[i0] : 0;
-
-                const x1 = ox + (w + 1) * dx;
-                const y1 = oy + hh * dy;
-                const z1 = isMesh ? -meshData[i1] : 0;
-
-                const x2 = ox + (w + 1) * dx;
-                const y2 = oy + (hh - 1) * dy;
-                const z2 = isMesh ? -meshData[i2] : 0;
-
-                const x3 = ox + w * dx;
-                const y3 = oy + (hh - 1) * dy;
-                const z3 = isMesh ? -meshData[i3] : 0;
-
-                if (i1_act && i3_act) {
-                    // diagonal i1, i3
-                    if (i0_act) {
-                        indices.push(i1, i3, i0); // t1 - i0 provoking index.
-
-                        line_positions.push(x0, y0, z0);
-                        line_positions.push(x3, y3, z3);
-
-                        line_positions.push(x0, y0, z0);
-                        line_positions.push(x1, y1, z1);
-                    }
-
-                    if (i2_act) {
-                        indices.push(i1, i3, i2); // t2 - i2 provoking index.
-
-                        line_positions.push(x2, y2, z2);
-                        line_positions.push(x3, y3, z3);
-
-                        line_positions.push(x2, y2, z2);
-                        line_positions.push(x1, y1, z1);
-                    }
-
-                    // diagonal
-                    if ((i0_act && !i2_act) || (!i0_act && i2_act)) {
-                        line_positions.push(x1, y1, z1);
-                        line_positions.push(x3, y3, z3);
-                    }
-                } else if (i0_act && i2_act) {
-                    // diagonal i0, i2
-                    if (i1_act) {
-                        indices.push(i1, i2, i0); // t1 - i0 provoking index.
-                    }
-
-                    if (i3_act) {
-                        indices.push(i3, i0, i2); // t2 - i2 provoking index.
-                    }
-
-                    // diagonal
-                    if ((i3_act && !i1_act) || (!i3_act && i1_act)) {
-                        line_positions.push(x0, y0, z0);
-                        line_positions.push(x2, y2, z2);
-                    }
-                }
-            }
-        }
-    } else {
-        // PROPERTIES IS SET CONSTANT OVER A CELL.
-        let i_indices = 0;
-        let i_vertices = 0;
-        // Loop over cells.
-        for (let h = 0; h < ny - 1; h++) {
-            for (let w = 0; w < nx - 1; w++) {
-                const hh = ny - 1 - h; // See note above.
-
-                const i0 = h * nx + w;
-                const i1 = h * nx + (w + 1);
-                const i2 = (h + 1) * nx + (w + 1);
-                const i3 = (h + 1) * nx + w;
-
-                const normal0 = calcNormal(w, h, nx, ny, isMesh, smoothShading, meshData, ox, oy);         // eslint-disable-line
-                const normal1 = calcNormal(w + 1, h, nx, ny, isMesh, smoothShading, meshData, ox, oy);     // eslint-disable-line
-                const normal2 = calcNormal(w + 1, h + 1, nx, ny, isMesh, smoothShading, meshData, ox, oy); // eslint-disable-line
-                const normal3 = calcNormal(w, h + 1, nx, ny, isMesh, smoothShading, meshData, ox, oy);     // eslint-disable-line
-
-                const i0_act = !isMesh || !isNaN(meshData[i0]); // eslint-disable-line
-                const i1_act = !isMesh || !isNaN(meshData[i1]); // eslint-disable-line
-                const i2_act = !isMesh || !isNaN(meshData[i2]); // eslint-disable-line
-                const i3_act = !isMesh || !isNaN(meshData[i3]); // eslint-disable-line
-
-                const x0 = ox + w * dx;
-                const y0 = oy + hh * dy;
-                const z0 = isMesh ? -meshData[i0] : 0;
-
-                const x1 = ox + (w + 1) * dx;
-                const y1 = oy + hh * dy;
-                const z1 = isMesh ? -meshData[i1] : 0;
-
-                const x2 = ox + (w + 1) * dx;
-                const y2 = oy + (hh - 1) * dy; // Note hh - 1 here.
-                const z2 = isMesh ? -meshData[i2] : 0;
-
-                const x3 = ox + w * dx;
-                const y3 = oy + (hh - 1) * dy; // Note hh - 1 here.
-                const z3 = isMesh ? -meshData[i3] : 0;
-
-                const propertyIndex = h * (nx - 1) + w; // (nx - 1) -> the width of the property 2D array is one less than for the nodes in this case.
-                const propertyValue = propertiesData[propertyIndex];
-
-                if (isNaN(propertyValue)) {
-                    // Inactive cell, dont draw.
-                    continue;
-                }
-
-                if (i1_act && i3_act) {
-                    // diagonal i1, i3
-                    if (i0_act) {
-                        // t1 - i0 provoking index.
-                        positions.push(x1, y1, z1);
-                        positions.push(x3, y3, z3);
-                        positions.push(x0, y0, z0);
-
-                        normals.push(normal1[0], normal1[1], normal1[2]);
-                        normals.push(normal3[0], normal3[1], normal3[2]);
-                        normals.push(normal0[0], normal0[1], normal0[2]);
-
-                        vertexIndexs.push(
-                            i_vertices++,
-                            i_vertices++,
-                            i_vertices++
-                        );
-
-                        indices.push(i_indices++, i_indices++, i_indices++);
-                        vertexProperties.push(propertyValue);
-                        vertexProperties.push(propertyValue);
-                        vertexProperties.push(propertyValue);
-
-                        line_positions.push(x0, y0, z0);
-                        line_positions.push(x3, y3, z3);
-
-                        line_positions.push(x0, y0, z0);
-                        line_positions.push(x1, y1, z1);
-                    }
-
-                    if (i2_act) {
-                        // t2 - i2 provoking index.
-                        positions.push(x1, y1, z1);
-                        positions.push(x3, y3, z3);
-                        positions.push(x2, y2, z2);
-
-                        normals.push(normal1[0], normal1[1], normal1[2]);
-                        normals.push(normal3[0], normal3[1], normal3[2]);
-                        normals.push(normal2[0], normal2[1], normal2[2]);
-
-                        vertexIndexs.push(
-                            i_vertices++,
-                            i_vertices++,
-                            i_vertices++
-                        );
-
-                        indices.push(i_indices++, i_indices++, i_indices++);
-                        vertexProperties.push(propertyValue);
-                        vertexProperties.push(propertyValue);
-                        vertexProperties.push(propertyValue);
-
-                        line_positions.push(x2, y2, z2);
-                        line_positions.push(x3, y3, z3);
-
-                        line_positions.push(x2, y2, z2);
-                        line_positions.push(x1, y1, z1);
-                    }
-
-                    // diagonal
-                    if ((i0_act && !i2_act) || (!i0_act && i2_act)) {
-                        line_positions.push(x1, y1, z1);
-                        line_positions.push(x3, y3, z3);
-                    }
-                } else if (i0_act && i2_act) {
-                    // diagonal i0, i2
-                    if (i1_act) {
-                        // t1 - i0 provoking index.
-                        positions.push(x1, y1, z1);
-                        positions.push(x2, y2, z2);
-                        positions.push(x0, y0, z0);
-
-                        normals.push(normal1[0], normal1[1], normal1[2]);
-                        normals.push(normal2[0], normal2[1], normal2[2]);
-                        normals.push(normal0[0], normal0[1], normal0[2]);
-
-                        vertexIndexs.push(
-                            i_vertices++,
-                            i_vertices++,
-                            i_vertices++
-                        );
-
-                        indices.push(i_indices++, i_indices++, i_indices++);
-                        vertexProperties.push(propertyValue);
-                        vertexProperties.push(propertyValue);
-                        vertexProperties.push(propertyValue);
-
-                        line_positions.push(x1, y1, z1);
-                        line_positions.push(x0, y0, z0);
-
-                        line_positions.push(x1, y1, z1);
-                        line_positions.push(x2, y2, z2);
-                    }
-
-                    if (i3_act) {
-                        // t2 - i2 provoking index.
-                        positions.push(x0, y0, z0);
-                        positions.push(x3, y3, z3);
-                        positions.push(x2, y2, z2);
-
-                        normals.push(normal0[0], normal0[1], normal0[2]);
-                        normals.push(normal3[0], normal3[1], normal3[2]);
-                        normals.push(normal2[0], normal2[1], normal2[2]);
-
-                        vertexIndexs.push(
-                            i_vertices++,
-                            i_vertices++,
-                            i_vertices++
-                        );
-
-                        indices.push(i_indices++, i_indices++, i_indices++);
-                        vertexProperties.push(propertyValue);
-                        vertexProperties.push(propertyValue);
-                        vertexProperties.push(propertyValue);
-
-                        line_positions.push(x3, y3, z3);
-                        line_positions.push(x0, y0, z0);
-
-                        line_positions.push(x3, y3, z3);
-                        line_positions.push(x2, y2, z2);
-                    }
-
-                    // diagonal
-                    if ((i3_act && !i1_act) || (!i3_act && i1_act)) {
-                        line_positions.push(x0, y0, z0);
-                        line_positions.push(x2, y2, z2);
-                    }
-                }
-            }
-        }
-    }
-
-    const mesh: MeshType = {
-        drawMode: 4, // corresponds to GL.TRIANGLES,
-        attributes: {
-            positions: { value: new Float32Array(positions), size: 3 },
-            normals: { value: new Float32Array(normals), size: 3 },
-            properties: { value: new Float32Array(vertexProperties), size: 1 },
-            vertex_indexs: { value: new Int32Array(vertexIndexs), size: 1 },
-        },
-        vertexCount: indices.length,
-        indices: { value: new Uint32Array(indices), size: 1 },
-    };
-
-    const mesh_lines: MeshTypeLines = {
-        drawMode: 1, // corresponds to GL.LINES,
-        attributes: {
-            positions: { value: new Float32Array(line_positions), size: 3 },
-        },
-        vertexCount: line_positions.length / 3,
-    };
-
-    //const t1 = performance.now();
-    // Keep this.
-    //console.log(`Task makeMesh took ${(t1 - t0) * 0.001}  seconds.`);
-
-    // Note: typescript gives this error "error TS2554: Expected 2-3 arguments, but got 1."
-    // Disabling this for now as the second argument should be optional.
-    // eslint-disable-next-line @typescript-eslint/ban-ts-comment
-    // @ts-ignore
-    postMessage([mesh, mesh_lines, meshZValueRange, propertyValueRange]);
-}

package/src/layers/northarrow/northArrow.stories.tsx

@@ -1,108 +0,0 @@
-import React from "react";
-import SubsurfaceViewer from "../../SubsurfaceViewer";
-import { ComponentStory, ComponentMeta } from "@storybook/react";
-
-export default {
-    component: SubsurfaceViewer,
-    title: "SubsurfaceViewer / North Arrow Layer",
-} as ComponentMeta<typeof SubsurfaceViewer>;
-
-type NumberQuad = [number, number, number, number];
-
-const white = [255, 255, 255, 255];
-
-const axes = {
-    "@@type": "AxesLayer",
-    id: "axes-layer2",
-    bounds: [432150, 6475800, 0, 439400, 6481500, 100],
-};
-
-const north_arrow_layer = {
-    "@@type": "NorthArrow3DLayer",
-    id: "north-arrow-layer",
-};
-
-const defaultParameters = {
-    docs: {
-        inlineStories: false,
-        iframeHeight: 500,
-    },
-};
-
-const defaultArgs = {
-    id: "map",
-    bounds: [432150, 6475800, 439400, 6481500] as NumberQuad,
-    layers: [axes, north_arrow_layer],
-    views: {
-        layout: [1, 1] as [number, number],
-        viewports: [
-            {
-                id: "view_1",
-                show3D: true,
-            },
-        ],
-    },
-};
-
-export const NorthArrow3d: ComponentStory<typeof SubsurfaceViewer> = (args) => {
-    return <SubsurfaceViewer {...args} />;
-};
-
-NorthArrow3d.args = {
-    ...defaultArgs,
-    views: {
-        layout: [1, 1],
-        viewports: [
-            {
-                id: "view_1",
-                show3D: true,
-            },
-        ],
-    },
-};
-
-NorthArrow3d.parameters = {
-    docs: {
-        ...defaultParameters.docs,
-        description: {
-            story: "Example using north arrow in 3D.",
-        },
-    },
-};
-
-export const NorthArrow2dDarkMode: ComponentStory<typeof SubsurfaceViewer> = (
-    args
-) => {
-    return <SubsurfaceViewer {...args} />;
-};
-
-NorthArrow2dDarkMode.args = {
-    ...defaultArgs,
-    layers: [
-        { ...axes, labelColor: white, axisColor: white },
-        { ...north_arrow_layer, color: white },
-    ],
-    views: {
-        layout: [1, 1],
-        viewports: [
-            {
-                id: "view_1",
-                show3D: false,
-            },
-        ],
-    },
-    scale: {
-        visible: true,
-        cssStyle: { color: "white" },
-    },
-};
-
-NorthArrow2dDarkMode.parameters = {
-    docs: {
-        ...defaultParameters.docs,
-        description: {
-            story: "Example using north arrow in 2D Dark Mode.",
-        },
-    },
-    backgrounds: { default: "dark" },
-};