@loaders.gl/gis 4.3.1 → 4.4.0-alpha.1

package/dist/lib/geometry-converters/wkb/convert-wkb-to-binary-geometry.js

@@ -0,0 +1,164 @@
+// loaders.gl
+// SPDX-License-Identifier: MIT
+// Copyright (c) vis.gl contributors
+import { concatenateBinaryPointGeometries, concatenateBinaryLineGeometries, concatenateBinaryPolygonGeometries } from "../../binary-geometry-api/concat-binary-geometry.js";
+import { concatTypedArrays } from "../../utils/concat-typed-arrays.js";
+import { WKBGeometryType } from "./helpers/wkb-types.js";
+import { parseWKBHeader } from "./helpers/parse-wkb-header.js";
+export function convertWKBToBinaryGeometry(arrayBuffer) {
+    const dataView = new DataView(arrayBuffer);
+    const wkbHeader = parseWKBHeader(dataView);
+    const { geometryType, dimensions, littleEndian } = wkbHeader;
+    const offset = wkbHeader.byteOffset;
+    switch (geometryType) {
+        case WKBGeometryType.Point:
+            const point = parsePoint(dataView, offset, dimensions, littleEndian);
+            return point.geometry;
+        case WKBGeometryType.LineString:
+            const line = parseLineString(dataView, offset, dimensions, littleEndian);
+            return line.geometry;
+        case WKBGeometryType.Polygon:
+            const polygon = parsePolygon(dataView, offset, dimensions, littleEndian);
+            return polygon.geometry;
+        case WKBGeometryType.MultiPoint:
+            const multiPoint = parseMultiPoint(dataView, offset, dimensions, littleEndian);
+            multiPoint.type = 'Point';
+            return multiPoint;
+        case WKBGeometryType.MultiLineString:
+            const multiLine = parseMultiLineString(dataView, offset, dimensions, littleEndian);
+            multiLine.type = 'LineString';
+            return multiLine;
+        case WKBGeometryType.MultiPolygon:
+            const multiPolygon = parseMultiPolygon(dataView, offset, dimensions, littleEndian);
+            multiPolygon.type = 'Polygon';
+            return multiPolygon;
+        // case WKBGeometryType.GeometryCollection:
+        // TODO: handle GeometryCollections
+        // return parseGeometryCollection(dataView, offset, dimensions, littleEndian);
+        default:
+            throw new Error(`WKB: Unsupported geometry type: ${geometryType}`);
+    }
+}
+// Primitives; parse point and linear ring
+function parsePoint(dataView, offset, dimension, littleEndian) {
+    const positions = new Float64Array(dimension);
+    for (let i = 0; i < dimension; i++) {
+        positions[i] = dataView.getFloat64(offset, littleEndian);
+        offset += 8;
+    }
+    return {
+        geometry: { type: 'Point', positions: { value: positions, size: dimension } },
+        offset
+    };
+}
+function parseLineString(dataView, offset, dimension, littleEndian) {
+    const nPoints = dataView.getUint32(offset, littleEndian);
+    offset += 4;
+    // Instantiate array
+    const positions = new Float64Array(nPoints * dimension);
+    for (let i = 0; i < nPoints * dimension; i++) {
+        positions[i] = dataView.getFloat64(offset, littleEndian);
+        offset += 8;
+    }
+    const pathIndices = [0];
+    if (nPoints > 0) {
+        pathIndices.push(nPoints);
+    }
+    return {
+        geometry: {
+            type: 'LineString',
+            positions: { value: positions, size: dimension },
+            pathIndices: { value: new Uint32Array(pathIndices), size: 1 }
+        },
+        offset
+    };
+}
+// https://stackoverflow.com/a/55261098
+const cumulativeSum = (sum) => (value) => (sum += value);
+function parsePolygon(dataView, offset, dimension, littleEndian) {
+    const nRings = dataView.getUint32(offset, littleEndian);
+    offset += 4;
+    const rings = [];
+    for (let i = 0; i < nRings; i++) {
+        const parsed = parseLineString(dataView, offset, dimension, littleEndian);
+        const { positions } = parsed.geometry;
+        offset = parsed.offset;
+        rings.push(positions.value);
+    }
+    const concatenatedPositions = new Float64Array(concatTypedArrays(rings).buffer);
+    const polygonIndices = [0];
+    if (concatenatedPositions.length > 0) {
+        polygonIndices.push(concatenatedPositions.length / dimension);
+    }
+    const primitivePolygonIndices = rings.map((l) => l.length / dimension).map(cumulativeSum(0));
+    primitivePolygonIndices.unshift(0);
+    return {
+        geometry: {
+            type: 'Polygon',
+            positions: { value: concatenatedPositions, size: dimension },
+            polygonIndices: {
+                value: new Uint32Array(polygonIndices),
+                size: 1
+            },
+            primitivePolygonIndices: { value: new Uint32Array(primitivePolygonIndices), size: 1 }
+        },
+        offset
+    };
+}
+function parseMultiPoint(dataView, offset, dimension, littleEndian) {
+    const nPoints = dataView.getUint32(offset, littleEndian);
+    offset += 4;
+    const binaryPointGeometries = [];
+    for (let i = 0; i < nPoints; i++) {
+        // Byte order for point
+        const littleEndianPoint = dataView.getUint8(offset) === 1;
+        offset++;
+        // Assert point type
+        if (dataView.getUint32(offset, littleEndianPoint) % 1000 !== 1) {
+            throw new Error('WKB: Inner geometries of MultiPoint not of type Point');
+        }
+        offset += 4;
+        const parsed = parsePoint(dataView, offset, dimension, littleEndianPoint);
+        offset = parsed.offset;
+        binaryPointGeometries.push(parsed.geometry);
+    }
+    return concatenateBinaryPointGeometries(binaryPointGeometries, dimension);
+}
+function parseMultiLineString(dataView, offset, dimension, littleEndian) {
+    const nLines = dataView.getUint32(offset, littleEndian);
+    offset += 4;
+    const binaryLineGeometries = [];
+    for (let i = 0; i < nLines; i++) {
+        // Byte order for line
+        const littleEndianLine = dataView.getUint8(offset) === 1;
+        offset++;
+        // Assert type LineString
+        if (dataView.getUint32(offset, littleEndianLine) % 1000 !== 2) {
+            throw new Error('WKB: Inner geometries of MultiLineString not of type LineString');
+        }
+        offset += 4;
+        const parsed = parseLineString(dataView, offset, dimension, littleEndianLine);
+        offset = parsed.offset;
+        binaryLineGeometries.push(parsed.geometry);
+    }
+    return concatenateBinaryLineGeometries(binaryLineGeometries, dimension);
+}
+function parseMultiPolygon(dataView, offset, dimension, littleEndian) {
+    const nPolygons = dataView.getUint32(offset, littleEndian);
+    offset += 4;
+    const binaryPolygonGeometries = [];
+    for (let i = 0; i < nPolygons; i++) {
+        // Byte order for polygon
+        const littleEndianPolygon = dataView.getUint8(offset) === 1;
+        offset++;
+        // Assert type Polygon
+        if (dataView.getUint32(offset, littleEndianPolygon) % 1000 !== 3) {
+            throw new Error('WKB: Inner geometries of MultiPolygon not of type Polygon');
+        }
+        offset += 4;
+        const parsed = parsePolygon(dataView, offset, dimension, littleEndianPolygon);
+        offset = parsed.offset;
+        binaryPolygonGeometries.push(parsed.geometry);
+    }
+    return concatenateBinaryPolygonGeometries(binaryPolygonGeometries, dimension);
+}

package/dist/lib/geometry-converters/wkb/convert-wkb-to-geometry.d.ts

@@ -0,0 +1,11 @@
+import type { Geometry } from '@loaders.gl/schema';
+export type convertWKBOptions = {
+    /** Does the GeoJSON input have Z values? */
+    hasZ?: boolean;
+    /** Does the GeoJSON input have M values? */
+    hasM?: boolean;
+    /** Spatial reference for input GeoJSON */
+    srid?: any;
+};
+export declare function convertWKBToGeometry(arrayBuffer: ArrayBuffer): Geometry;
+//# sourceMappingURL=convert-wkb-to-geometry.d.ts.map

package/dist/lib/geometry-converters/wkb/convert-wkb-to-geometry.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"convert-wkb-to-geometry.d.ts","sourceRoot":"","sources":["../../../../src/lib/geometry-converters/wkb/convert-wkb-to-geometry.ts"],"names":[],"mappings":"AAIA,OAAO,KAAK,EAAC,QAAQ,EAAC,MAAM,oBAAoB,CAAC;AAIjD,MAAM,MAAM,iBAAiB,GAAG;IAC9B,4CAA4C;IAC5C,IAAI,CAAC,EAAE,OAAO,CAAC;IAEf,4CAA4C;IAC5C,IAAI,CAAC,EAAE,OAAO,CAAC;IAEf,0CAA0C;IAC1C,IAAI,CAAC,EAAE,GAAG,CAAC;CACZ,CAAC;AAEF,wBAAgB,oBAAoB,CAAC,WAAW,EAAE,WAAW,GAAG,QAAQ,CAGvE"}

package/dist/lib/geometry-converters/wkb/convert-wkb-to-geometry.js

@@ -0,0 +1,9 @@
+// loaders.gl
+// SPDX-License-Identifier: MIT
+// Copyright (c) vis.gl contributors
+import { convertBinaryGeometryToGeometry } from "../convert-binary-geometry-to-geojson.js";
+import { convertWKBToBinaryGeometry } from "./convert-wkb-to-binary-geometry.js";
+export function convertWKBToGeometry(arrayBuffer) {
+    const binaryGeometry = convertWKBToBinaryGeometry(arrayBuffer);
+    return convertBinaryGeometryToGeometry(binaryGeometry);
+}

package/dist/lib/geometry-converters/wkb/convert-wkt-to-geometry.d.ts

@@ -0,0 +1,16 @@
+import type { Geometry } from '@loaders.gl/schema';
+export type ParseWKTOptions = {
+    wkt?: {
+        /** Whether to add any CRS, if found, as undocumented CRS property on the return geometry */
+        crs?: boolean;
+    };
+};
+/**
+ * Parse WKT and return GeoJSON.
+ * @param input A WKT geometry string
+ * @return A GeoJSON geometry object
+ *
+ * @note We only support the "geojson" subset of the OGC simple features standard
+ **/
+export declare function convertWKTToGeometry(input: string, options?: ParseWKTOptions): Geometry | null;
+//# sourceMappingURL=convert-wkt-to-geometry.d.ts.map

package/dist/lib/geometry-converters/wkb/convert-wkt-to-geometry.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"convert-wkt-to-geometry.d.ts","sourceRoot":"","sources":["../../../../src/lib/geometry-converters/wkb/convert-wkt-to-geometry.ts"],"names":[],"mappings":"AAKA,OAAO,KAAK,EAAC,QAAQ,EAAC,MAAM,oBAAoB,CAAC;AASjD,MAAM,MAAM,eAAe,GAAG;IAC5B,GAAG,CAAC,EAAE;QACJ,4FAA4F;QAC5F,GAAG,CAAC,EAAE,OAAO,CAAC;KACf,CAAC;CACH,CAAC;AASF;;;;;;IAMI;AACJ,wBAAgB,oBAAoB,CAAC,KAAK,EAAE,MAAM,EAAE,OAAO,CAAC,EAAE,eAAe,GAAG,QAAQ,GAAG,IAAI,CAU9F"}

package/dist/lib/geometry-converters/wkb/convert-wkt-to-geometry.js

@@ -0,0 +1,258 @@
+// loaders.gl
+// SPDX-License-Identifier: MIT
+// Copyright (c) vis.gl contributors
+// Fork of https://github.com/mapbox/wellknown under ISC license (MIT/BSD-2-clause equivalent)
+/* eslint-disable */
+// @ts-nocheck
+const numberRegexp = /[-+]?([0-9]*\.[0-9]+|[0-9]+)([eE][-+]?[0-9]+)?/;
+// Matches sequences like '100 100' or '100 100 100'.
+const tuples = new RegExp('^' + numberRegexp.source + '(\\s' + numberRegexp.source + '){1,}');
+/**
+ * Parse WKT and return GeoJSON.
+ * @param input A WKT geometry string
+ * @return A GeoJSON geometry object
+ *
+ * @note We only support the "geojson" subset of the OGC simple features standard
+ **/
+export function convertWKTToGeometry(input, options) {
+    const parts = input.split(';');
+    let _ = parts.pop();
+    const srid = (parts.shift() || '').split('=').pop();
+    const state = { parts, _, i: 0 };
+    const geometry = parseGeometry(state);
+    return options?.wkt?.crs ? addCRS(geometry, srid) : geometry;
+}
+function parseGeometry(state) {
+    return (parsePoint(state) ||
+        parseLineString(state) ||
+        parsePolygon(state) ||
+        parseMultiPoint(state) ||
+        parseMultiLineString(state) ||
+        parseMultiPolygon(state) ||
+        parseGeometryCollection(state));
+}
+/** Adds a coordinate reference system as an undocumented  */
+function addCRS(obj, srid) {
+    if (obj && srid?.match(/\d+/)) {
+        const crs = {
+            type: 'name',
+            properties: {
+                name: 'urn:ogc:def:crs:EPSG::' + srid
+            }
+        };
+        // @ts-expect-error we assign an undocumented property on the geometry
+        obj.crs = crs;
+    }
+    return obj;
+}
+// GEOMETRIES
+function parsePoint(state) {
+    if (!$(/^(POINT(\sz)?)/i, state)) {
+        return null;
+    }
+    white(state);
+    if (!$(/^(\()/, state)) {
+        return null;
+    }
+    const c = coords(state);
+    if (!c) {
+        return null;
+    }
+    white(state);
+    if (!$(/^(\))/, state)) {
+        return null;
+    }
+    return {
+        type: 'Point',
+        coordinates: c[0]
+    };
+}
+function parseMultiPoint(state) {
+    if (!$(/^(MULTIPOINT)/i, state)) {
+        return null;
+    }
+    white(state);
+    const newCoordsFormat = state._?.substring(state._?.indexOf('(') + 1, state._.length - 1)
+        .replace(/\(/g, '')
+        .replace(/\)/g, '');
+    state._ = 'MULTIPOINT (' + newCoordsFormat + ')';
+    const c = multicoords(state);
+    if (!c) {
+        return null;
+    }
+    white(state);
+    return {
+        type: 'MultiPoint',
+        coordinates: c
+    };
+}
+function parseLineString(state) {
+    if (!$(/^(LINESTRING(\sz)?)/i, state)) {
+        return null;
+    }
+    white(state);
+    if (!$(/^(\()/, state)) {
+        return null;
+    }
+    const c = coords(state);
+    if (!c) {
+        return null;
+    }
+    if (!$(/^(\))/, state)) {
+        return null;
+    }
+    return {
+        type: 'LineString',
+        coordinates: c
+    };
+}
+function parseMultiLineString(state) {
+    if (!$(/^(MULTILINESTRING)/i, state))
+        return null;
+    white(state);
+    const c = multicoords(state);
+    if (!c) {
+        return null;
+    }
+    white(state);
+    return {
+        // @ts-ignore
+        type: 'MultiLineString',
+        // @ts-expect-error
+        coordinates: c
+    };
+}
+function parsePolygon(state) {
+    if (!$(/^(POLYGON(\sz)?)/i, state)) {
+        return null;
+    }
+    white(state);
+    const c = multicoords(state);
+    if (!c) {
+        return null;
+    }
+    return {
+        // @ts-ignore
+        type: 'Polygon',
+        // @ts-expect-error
+        coordinates: c
+    };
+}
+function parseMultiPolygon(state) {
+    if (!$(/^(MULTIPOLYGON)/i, state)) {
+        return null;
+    }
+    white(state);
+    const c = multicoords(state);
+    if (!c) {
+        return null;
+    }
+    return {
+        type: 'MultiPolygon',
+        // @ts-expect-error
+        coordinates: c
+    };
+}
+function parseGeometryCollection(state) {
+    const geometries = [];
+    let geometry;
+    if (!$(/^(GEOMETRYCOLLECTION)/i, state)) {
+        return null;
+    }
+    white(state);
+    if (!$(/^(\()/, state)) {
+        return null;
+    }
+    while ((geometry = parseGeometry(state))) {
+        geometries.push(geometry);
+        white(state);
+        $(/^(,)/, state);
+        white(state);
+    }
+    if (!$(/^(\))/, state)) {
+        return null;
+    }
+    return {
+        type: 'GeometryCollection',
+        geometries: geometries
+    };
+}
+// COORDINATES
+function multicoords(state) {
+    white(state);
+    let depth = 0;
+    const rings = [];
+    const stack = [rings];
+    let pointer = rings;
+    let elem;
+    while ((elem = $(/^(\()/, state) || $(/^(\))/, state) || $(/^(,)/, state) || $(tuples, state))) {
+        if (elem === '(') {
+            stack.push(pointer);
+            pointer = [];
+            stack[stack.length - 1].push(pointer);
+            depth++;
+        }
+        else if (elem === ')') {
+            // For the case: Polygon(), ...
+            if (pointer.length === 0)
+                return null;
+            // @ts-ignore
+            pointer = stack.pop();
+            // the stack was empty, input was malformed
+            if (!pointer)
+                return null;
+            depth--;
+            if (depth === 0)
+                break;
+        }
+        else if (elem === ',') {
+            pointer = [];
+            stack[stack.length - 1].push(pointer);
+        }
+        else if (!elem.split(/\s/g).some(isNaN)) {
+            Array.prototype.push.apply(pointer, elem.split(/\s/g).map(parseFloat));
+        }
+        else {
+            return null;
+        }
+        white(state);
+    }
+    if (depth !== 0)
+        return null;
+    return rings;
+}
+function coords(state) {
+    const list = [];
+    let item;
+    let pt;
+    while ((pt = $(tuples, state) || $(/^(,)/, state))) {
+        if (pt === ',') {
+            list.push(item);
+            item = [];
+        }
+        else if (!pt.split(/\s/g).some(isNaN)) {
+            if (!item)
+                item = [];
+            Array.prototype.push.apply(item, pt.split(/\s/g).map(parseFloat));
+        }
+        white(state);
+    }
+    if (item)
+        list.push(item);
+    else
+        return null;
+    return list.length ? list : null;
+}
+// HELPERS
+function $(regexp, state) {
+    const match = state._?.substring(state.i).match(regexp);
+    if (!match)
+        return null;
+    else {
+        state.i += match[0].length;
+        return match[0];
+    }
+}
+function white(state) {
+    $(/^\s*/, state);
+}

package/dist/lib/geometry-converters/wkb/helpers/parse-wkb-header.d.ts

@@ -0,0 +1,43 @@
+import type { WKBGeometryType, WKBHeader } from "./wkb-types.js";
+/**
+ * Check if a string is WKT.
+ * @param input A potential WKT geometry string
+ * @return `true` if input appears to be a WKT geometry string, `false` otherwise
+
+ * @note We only support the "geojson" subset of the OGC simple features standard
+ * @todo Does not handle leading spaces which appear to be permitted per the spec:
+ * "A WKT string contains no white space outside of double quotes.
+ * However padding with white space to improve human readability is permitted;
+ * the examples of WKT that are included in this document have
+ * spaces and line feeds inserted to improve clarity. Any padding is stripped out or ignored by parsers."
+ */
+export declare function isWKT(input: string | ArrayBuffer): boolean;
+/**
+ * Get the geometry type of a WKT string.
+ * @param input A potential WKT geometry string
+ * @return `true` if input appears to be a WKT geometry string, `false` otherwise
+
+ * @note We only support the "geojson" subset of the OGC simple features standard
+ * @todo Does not handle leading spaces which appear to be permitted per the spec:
+ * "A WKT string contains no white space outside of double quotes.
+ * However padding with white space to improve human readability is permitted;
+ * the examples of WKT that are included in this document have
+ * spaces and line feeds inserted to improve clarity. Any padding is stripped out or ignored by parsers."
+ */
+export declare function getWKTGeometryType(input: string | ArrayBuffer): WKBGeometryType | null;
+/**
+ * Check if an array buffer might be a TWKB array buffer
+ * @param arrayBuffer The array buffer to check
+ * @returns false if this is definitely not a TWKB array buffer, true if it might be a TWKB array buffer
+ */
+export declare function isTWKB(arrayBuffer: ArrayBuffer): boolean;
+/** Sanity checks that first to 5-9 bytes could represent a supported WKB dialect header */
+export declare function isWKB(arrayBuffer: ArrayBuffer): boolean;
+/**
+ * Parses header and provides a byteOffset to start of geometry data
+ * @param dataView
+ * @param target optionally supply a WKBHeader object to avoid creating a new object for every call
+ * @returns a header object describing the WKB data
+ */
+export declare function parseWKBHeader(dataView: DataView, target?: WKBHeader): WKBHeader;
+//# sourceMappingURL=parse-wkb-header.d.ts.map

package/dist/lib/geometry-converters/wkb/helpers/parse-wkb-header.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"parse-wkb-header.d.ts","sourceRoot":"","sources":["../../../../../src/lib/geometry-converters/wkb/helpers/parse-wkb-header.ts"],"names":[],"mappings":"AAIA,OAAO,KAAK,EAAC,eAAe,EAAE,SAAS,EAAC,uBAAoB;AAU5D;;;;;;;;;;;GAWG;AACH,wBAAgB,KAAK,CAAC,KAAK,EAAE,MAAM,GAAG,WAAW,GAAG,OAAO,CAE1D;AAED;;;;;;;;;;;GAWG;AACH,wBAAgB,kBAAkB,CAAC,KAAK,EAAE,MAAM,GAAG,WAAW,GAAG,eAAe,GAAG,IAAI,CAUtF;AAED;;;;GAIG;AACH,wBAAgB,MAAM,CAAC,WAAW,EAAE,WAAW,GAAG,OAAO,CAaxD;AAED,2FAA2F;AAC3F,wBAAgB,KAAK,CAAC,WAAW,EAAE,WAAW,GAAG,OAAO,CAkDvD;AAED;;;;;GAKG;AAEH,wBAAgB,cAAc,CAAC,QAAQ,EAAE,QAAQ,EAAE,MAAM,CAAC,EAAE,SAAS,GAAG,SAAS,CA8EhF"}

package/dist/lib/geometry-converters/wkb/helpers/parse-wkb-header.js

@@ -0,0 +1,174 @@
+// loaders.gl
+// SPDX-License-Identifier: MIT
+// Copyright (c) vis.gl contributors
+import { EWKB_FLAG_Z, EWKB_FLAG_M, EWKB_FLAG_SRID, MAX_SRID, WKT_MAGIC_STRINGS, WKT_MAGIC_BYTES } from "./wkb-types.js";
+/**
+ * Check if a string is WKT.
+ * @param input A potential WKT geometry string
+ * @return `true` if input appears to be a WKT geometry string, `false` otherwise
+
+ * @note We only support the "geojson" subset of the OGC simple features standard
+ * @todo Does not handle leading spaces which appear to be permitted per the spec:
+ * "A WKT string contains no white space outside of double quotes.
+ * However padding with white space to improve human readability is permitted;
+ * the examples of WKT that are included in this document have
+ * spaces and line feeds inserted to improve clarity. Any padding is stripped out or ignored by parsers."
+ */
+export function isWKT(input) {
+    return getWKTGeometryType(input) !== null;
+}
+/**
+ * Get the geometry type of a WKT string.
+ * @param input A potential WKT geometry string
+ * @return `true` if input appears to be a WKT geometry string, `false` otherwise
+
+ * @note We only support the "geojson" subset of the OGC simple features standard
+ * @todo Does not handle leading spaces which appear to be permitted per the spec:
+ * "A WKT string contains no white space outside of double quotes.
+ * However padding with white space to improve human readability is permitted;
+ * the examples of WKT that are included in this document have
+ * spaces and line feeds inserted to improve clarity. Any padding is stripped out or ignored by parsers."
+ */
+export function getWKTGeometryType(input) {
+    if (typeof input === 'string') {
+        const index = WKT_MAGIC_STRINGS.findIndex((magicString) => input.startsWith(magicString));
+        return index >= 0 ? (index + 1) : null;
+    }
+    const inputArray = new Uint8Array(input);
+    const index = WKT_MAGIC_BYTES.findIndex((magicBytes) => magicBytes.every((value, index) => value === inputArray[index]));
+    return index >= 0 ? (index + 1) : null;
+}
+/**
+ * Check if an array buffer might be a TWKB array buffer
+ * @param arrayBuffer The array buffer to check
+ * @returns false if this is definitely not a TWKB array buffer, true if it might be a TWKB array buffer
+ */
+export function isTWKB(arrayBuffer) {
+    const dataView = new DataView(arrayBuffer);
+    const byteOffset = 0;
+    const type = dataView.getUint8(byteOffset);
+    const geometryType = type & 0x0f;
+    // Only geometry types 1 to 7 (point to geometry collection are currently defined)
+    if (geometryType < 1 || geometryType > 7) {
+        return false;
+    }
+    return true;
+}
+/** Sanity checks that first to 5-9 bytes could represent a supported WKB dialect header */
+export function isWKB(arrayBuffer) {
+    const dataView = new DataView(arrayBuffer);
+    let byteOffset = 0;
+    const endianness = dataView.getUint8(byteOffset);
+    byteOffset += 1;
+    // Check valid endianness (only 0 or 1 are allowed)
+    if (endianness > 1) {
+        return false;
+    }
+    const littleEndian = endianness === 1;
+    const geometry = dataView.getUint32(byteOffset, littleEndian);
+    byteOffset += 4;
+    // check valid geometry type (we don't support extension geometries)
+    const geometryType = geometry & 0x07;
+    if (geometryType === 0 || geometryType > 7) {
+        return false;
+    }
+    const geometryFlags = geometry - geometryType;
+    // Accept iso-wkb flags
+    if (geometryFlags === 0 ||
+        geometryFlags === 1000 ||
+        geometryFlags === 2000 ||
+        geometryFlags === 3000) {
+        return true;
+    }
+    // Accept ewkb flags but reject otherwise
+    if ((geometryFlags & ~(EWKB_FLAG_Z | EWKB_FLAG_M | EWKB_FLAG_SRID)) !== 0) {
+        return false;
+    }
+    if (geometryFlags & EWKB_FLAG_SRID) {
+        const srid = dataView.getUint32(byteOffset, littleEndian);
+        byteOffset += 4;
+        if (srid > MAX_SRID) {
+            return false;
+        }
+    }
+    return true;
+}
+/**
+ * Parses header and provides a byteOffset to start of geometry data
+ * @param dataView
+ * @param target optionally supply a WKBHeader object to avoid creating a new object for every call
+ * @returns a header object describing the WKB data
+ */
+// eslint-disable-next-line max-statements, complexity
+export function parseWKBHeader(dataView, target) {
+    const wkbHeader = Object.assign(target || {}, {
+        type: 'wkb',
+        variant: 'wkb',
+        geometryType: 1,
+        dimensions: 2,
+        coordinates: 'xy',
+        littleEndian: true,
+        byteOffset: 0
+    });
+    if (isWKT(dataView.buffer)) {
+        // TODO - WKB header could include WKT type
+        throw new Error('WKB: Cannot parse WKT data');
+    }
+    // Check endianness of data
+    wkbHeader.littleEndian = dataView.getUint8(wkbHeader.byteOffset) === 1;
+    wkbHeader.byteOffset++;
+    // 4-digit code representing dimension and type of geometry
+    const geometryCode = dataView.getUint32(wkbHeader.byteOffset, wkbHeader.littleEndian);
+    wkbHeader.byteOffset += 4;
+    wkbHeader.geometryType = (geometryCode & 0x7);
+    // Check if iso-wkb variant: iso-wkb adds 1000, 2000 or 3000 to the geometry code
+    const isoType = (geometryCode - wkbHeader.geometryType) / 1000;
+    switch (isoType) {
+        case 0:
+            break;
+        case 1:
+            wkbHeader.variant = 'iso-wkb';
+            wkbHeader.dimensions = 3;
+            wkbHeader.coordinates = 'xyz';
+            break;
+        case 2:
+            wkbHeader.variant = 'iso-wkb';
+            wkbHeader.dimensions = 3;
+            wkbHeader.coordinates = 'xym';
+            break;
+        case 3:
+            wkbHeader.variant = 'iso-wkb';
+            wkbHeader.dimensions = 4;
+            wkbHeader.coordinates = 'xyzm';
+            break;
+        default:
+            throw new Error(`WKB: Unsupported iso-wkb type: ${isoType}`);
+    }
+    // Check if EWKB variant. Uses bitmasks for Z&M dimensions as well as optional SRID field
+    const ewkbZ = geometryCode & EWKB_FLAG_Z;
+    const ewkbM = geometryCode & EWKB_FLAG_M;
+    const ewkbSRID = geometryCode & EWKB_FLAG_SRID;
+    if (ewkbZ && ewkbM) {
+        wkbHeader.variant = 'ewkb';
+        wkbHeader.dimensions = 4;
+        wkbHeader.coordinates = 'xyzm';
+    }
+    else if (ewkbZ) {
+        wkbHeader.variant = 'ewkb';
+        wkbHeader.dimensions = 3;
+        wkbHeader.coordinates = 'xyz';
+    }
+    else if (ewkbM) {
+        wkbHeader.variant = 'ewkb';
+        wkbHeader.dimensions = 3;
+        wkbHeader.coordinates = 'xym';
+    }
+    // If SRID present read four more bytes
+    if (ewkbSRID) {
+        wkbHeader.variant = 'ewkb';
+        // 4-digit code representing dimension and type of geometry
+        wkbHeader.srid = dataView.getUint32(wkbHeader.byteOffset, wkbHeader.littleEndian);
+        wkbHeader.byteOffset += 4;
+    }
+    return wkbHeader;
+}