landxml 0.5.1 → 0.6.0

package/CHANGELOG.md

@@ -1,5 +1,17 @@
 # landxml
 
+## 0.6.0
+
+### Minor Changes
+
+- 330e3e8: Introduced web workers for some of the processing heavy workflows
+
+## 0.5.2
+
+### Patch Changes
+
+- 6cafe85: Bugfix parsing landxml's without extra triangulated faces
+
 ## 0.5.1
 
 ### Patch Changes

package/dist/{index.mjs → index.cjs}

@@ -1,7 +1,12 @@
+"use strict";
+var __create = Object.create;
 var __defProp = Object.defineProperty;
 var __defProps = Object.defineProperties;
+var __getOwnPropDesc = Object.getOwnPropertyDescriptor;
 var __getOwnPropDescs = Object.getOwnPropertyDescriptors;
+var __getOwnPropNames = Object.getOwnPropertyNames;
 var __getOwnPropSymbols = Object.getOwnPropertySymbols;
+var __getProtoOf = Object.getPrototypeOf;
 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;
@@ -29,6 +34,27 @@ var __objRest = (source, exclude) => {
     }
   return target;
 };
+var __export = (target, all) => {
+  for (var name in all)
+    __defProp(target, name, { get: all[name], enumerable: true });
+};
+var __copyProps = (to, from, except, desc) => {
+  if (from && typeof from === "object" || typeof from === "function") {
+    for (let key of __getOwnPropNames(from))
+      if (!__hasOwnProp.call(to, key) && key !== except)
+        __defProp(to, key, { get: () => from[key], enumerable: !(desc = __getOwnPropDesc(from, key)) || desc.enumerable });
+  }
+  return to;
+};
+var __toESM = (mod, isNodeMode, target) => (target = mod != null ? __create(__getProtoOf(mod)) : {}, __copyProps(
+  // If the importer is in node compatibility mode or this is not an ESM
+  // file that has been converted to a CommonJS file using a Babel-
+  // compatible transform (i.e. "__esModule" has not been set), then set
+  // "default" to the CommonJS "module.exports" for node compatibility.
+  isNodeMode || !mod || !mod.__esModule ? __defProp(target, "default", { value: mod, enumerable: true }) : target,
+  mod
+));
+var __toCommonJS = (mod) => __copyProps(__defProp({}, "__esModule", { value: true }), mod);
 var __async = (__this, __arguments, generator) => {
   return new Promise((resolve, reject) => {
     var fulfilled = (value) => {
@@ -50,6 +76,15 @@ var __async = (__this, __arguments, generator) => {
   });
 };
 
+// src/index.ts
+var src_exports = {};
+__export(src_exports, {
+  reprojectGeoJson: () => reproject_geojson_default,
+  toGeojsonContours: () => to_geojson_contours_default,
+  toGlb: () => to_glb_default
+});
+module.exports = __toCommonJS(src_exports);
+
 // src/private/filter-by-surfaceId.ts
 var filterBySurfaceId = (parsedSurfaces, surfaceId) => {
   let filtered = [...parsedSurfaces];
@@ -68,7 +103,7 @@ var filterBySurfaceId = (parsedSurfaces, surfaceId) => {
 var filter_by_surfaceId_default = filterBySurfaceId;
 
 // src/private/get-glb.ts
-import { Document, WebIO } from "@gltf-transform/core";
+var import_core = require("@gltf-transform/core");
 var findXYAxisMedians = (vertices) => {
   var _a, _b;
   vertices = vertices.slice().filter(Boolean);
@@ -83,7 +118,7 @@ var getGlb = (data, customCenter) => __async(void 0, null, function* () {
     return [x - center[0], z, -(y - center[1])];
   }).reduce((prev, curr) => prev.concat(curr), []);
   const triangles = data.surfaceDefinition.faces.reduce((prev, curr) => prev.concat(curr), []);
-  const doc = new Document();
+  const doc = new import_core.Document();
   const buffer = doc.createBuffer();
   const position = doc.createAccessor().setType("VEC3").setArray(new Float32Array(vertices)).setBuffer(buffer);
   const indices = doc.createAccessor().setType("SCALAR").setArray(new Uint32Array(triangles)).setBuffer(buffer);
@@ -91,7 +126,7 @@ var getGlb = (data, customCenter) => __async(void 0, null, function* () {
   const mesh = doc.createMesh().addPrimitive(prim);
   const node = doc.createNode().setMesh(mesh);
   const scene = doc.createScene().addChild(node);
-  const glb = yield new WebIO().writeBinary(doc);
+  const glb = yield new import_core.WebIO().writeBinary(doc);
   return { glb, center };
 });
 var get_glb_default = getGlb;
@@ -110,47 +145,109 @@ var downloadGlb = (glbData, fileName) => {
 var download_glb_default = downloadGlb;
 
 // src/private/parse-xml.ts
-import convert from "xml-js";
-var parseXML = (xmlString) => __async(void 0, null, function* () {
-  const { LandXML } = convert.xml2js(xmlString, {
-    compact: true,
-    attributesKey: "attr",
-    textKey: "content"
-  });
-  const landXML_surfaces = Array.isArray(LandXML.Surfaces.Surface) ? LandXML.Surfaces.Surface : [LandXML.Surfaces.Surface];
-  return landXML_surfaces.map((surface) => {
-    var _a, _b;
-    let points = [];
-    let faces = [];
-    let pointIdMap = {};
-    surface.Definition.Pnts.P.forEach((pt) => {
-      const { attr, content } = pt;
-      const [y, x, z] = content.split(" ").map((v) => parseFloat(v));
-      points.push([x, y, z]);
-      pointIdMap[attr.id] = points.length - 1;
-    }, []);
-    faces = surface.Definition.Faces.F.reduce((faceList, face) => {
-      const { attr, content } = face;
-      if (attr.i === "1")
-        return faceList;
-      const [a, b, c] = content.split(" ").map((v) => pointIdMap[v]);
-      if ([a, b, c].filter((v) => typeof v === "undefined").length > 0) {
-        throw `Invalid LandXML. A face is referencing a point that doesn't exist. Face is referencing points: ${content}`;
-      }
-      return faceList.concat([[a, b, c]]);
-    }, []);
-    return {
-      sourceFile: LandXML.Project.attr.name,
-      timeStamp: LandXML.Application.attr.timeStamp,
-      name: surface.attr.name,
-      description: surface.attr.desc,
-      wktString: (_b = (_a = LandXML.CoordinateSystem) == null ? void 0 : _a.attr) == null ? void 0 : _b.ogcWktCode,
-      surfaceDefinition: {
-        points,
-        faces
+var import_easy_web_worker = require("easy-web-worker");
+var import_xml_js = __toESM(require("xml-js"), 1);
+var surfaceDefWorker = (0, import_easy_web_worker.createEasyWebWorker)(
+  ({ onMessage }) => {
+    onMessage((message) => {
+      const { isPoint, arr, idMap } = message.payload;
+      if (isPoint) {
+        message.resolve(
+          arr.map((pt) => [pt.attr.id, pt.content.split(" ").map(Number)]).map((v) => [v[0], [v[1][1], v[1][0], v[1][2]]])
+        );
+      } else {
+        message.resolve(
+          arr.flatMap((f) => {
+            var _a;
+            if (typeof f === "string")
+              return [f.split(" ").map((id) => idMap == null ? void 0 : idMap.indexOf(id))];
+            if (((_a = f == null ? void 0 : f.attr) == null ? void 0 : _a.i) === "1")
+              return [];
+            return [f.content.split(" ").map((id) => idMap == null ? void 0 : idMap.indexOf(id))];
+          })
+        );
       }
-    };
-  });
+    });
+  },
+  { maxWorkers: 10 }
+);
+var parseXML = (xmlString) => __async(void 0, null, function* () {
+  return new Promise((resolve, reject) => __async(void 0, null, function* () {
+    var _a, _b, _c, _d, _e;
+    const parsed = import_xml_js.default.xml2js(xmlString, {
+      compact: true,
+      attributesKey: "attr",
+      textKey: "content"
+    });
+    if (typeof ((_b = (_a = parsed.LandXML) == null ? void 0 : _a.Surfaces) == null ? void 0 : _b.Surface) === "undefined") {
+      throw new Error("LandXML doesn't contain any surfaces");
+      return;
+    }
+    if (!Array.isArray(parsed.LandXML.Surfaces.Surface)) {
+      parsed.LandXML.Surfaces.Surface = [parsed.LandXML.Surfaces.Surface];
+    }
+    let sourceFile = parsed.LandXML.Project.attr.name || "Undefined source";
+    let timeStamp = parsed.LandXML.Application.attr.timeStamp || "";
+    let wktString = ((_e = (_d = (_c = parsed.LandXML) == null ? void 0 : _c.CoordinateSystem) == null ? void 0 : _d.attr) == null ? void 0 : _e.ogcWktCode) || void 0;
+    const surfaces = parsed.LandXML.Surfaces.Surface.map(
+      (surface) => __async(void 0, null, function* () {
+        return new Promise((resolve2, reject2) => __async(void 0, null, function* () {
+          const { name, desc } = surface.attr;
+          const Pnts = surface.Definition.Pnts.P;
+          const Faces = surface.Definition.Faces.F;
+          let ptsIdArray = [];
+          let faces = [];
+          if (Pnts.length > 1e4) {
+            const sliceIndexes = [...Array(20).keys()].map((i) => Math.round(Pnts.length / 20 * i));
+            ptsIdArray = (yield Promise.all(
+              sliceIndexes.map(
+                (v, i, a) => new Promise((resolve3, reject3) => __async(void 0, null, function* () {
+                  const pts = yield surfaceDefWorker.send({
+                    isPoint: true,
+                    arr: Pnts.slice(a[i], a[i + 1] || Pnts.length)
+                  });
+                  resolve3(pts);
+                }))
+              )
+            )).reduce((prev, curr) => [...prev, ...curr], []);
+          } else {
+            ptsIdArray = yield surfaceDefWorker.send({ arr: Pnts, isPoint: true });
+          }
+          const points = ptsIdArray.map((v) => v[1]);
+          const pointsIdMap = ptsIdArray.map((v) => v[0]);
+          if (Faces.length > 1e4) {
+            const sliceIndexes = [...Array(20).keys()].map((i) => Math.round(Faces.length / 20 * i));
+            faces = (yield Promise.all(
+              sliceIndexes.map(
+                (v, i, a) => new Promise((resolve3, reject3) => __async(void 0, null, function* () {
+                  const fcs = yield surfaceDefWorker.send({
+                    isPoint: false,
+                    arr: Faces.slice(a[i], a[i + 1] || Faces.length),
+                    idMap: pointsIdMap
+                  });
+                  resolve3(fcs);
+                }))
+              )
+            )).reduce((prev, curr) => [...prev, ...curr], []);
+          } else {
+            faces = yield surfaceDefWorker.send({ arr: Faces, isPoint: false, idMap: pointsIdMap });
+          }
+          resolve2({
+            sourceFile,
+            timeStamp: timeStamp || "",
+            name,
+            description: desc || "",
+            wktString,
+            surfaceDefinition: {
+              points,
+              faces
+            }
+          });
+        }));
+      })
+    );
+    resolve(yield Promise.all(surfaces));
+  }));
 });
 var parse_xml_default = parseXML;
 
@@ -182,32 +279,104 @@ var toGlb = (landXmlString, center = "auto", surfaceId = -1) => __async(void 0,
 var to_glb_default = toGlb;
 
 // src/private/get-contours.ts
-var contourLineOnFace = (face, z) => {
-  let vertsAtElevation = 0;
-  let line = [];
-  for (let i = 0; i < face.length; i++) {
-    let vertex1 = face[i];
-    let vertex2 = face[(i + 1) % face.length];
-    if (vertex1[2] === z)
-      vertsAtElevation++;
-    if ((vertex1[2] <= z && vertex2[2] >= z || vertex1[2] >= z && vertex2[2] <= z) && !Number.isNaN((z - vertex1[2]) / (vertex2[2] - vertex1[2]))) {
-      let t = (z - vertex1[2]) / (vertex2[2] - vertex1[2]);
-      line.push([vertex1[0] + t * (vertex2[0] - vertex1[0]), vertex1[1] + t * (vertex2[1] - vertex1[1])]);
-    }
-  }
-  if (vertsAtElevation >= 2 && face.map((f) => f[2]).reduce((a, b) => a + b) > z * face.length)
-    return void 0;
-  if (line.length === 2 && line[0][0] === line[1][0] && line[0][1] === line[1][1])
-    return void 0;
-  if (line.length > 2) {
-    line = [...new Set(line.map((v) => JSON.stringify(v)))].map((s) => JSON.parse(s));
-  }
-  return line.length > 0 ? line : void 0;
-};
+var import_easy_web_worker2 = require("easy-web-worker");
+var contoursWorker = (0, import_easy_web_worker2.createEasyWebWorker)(
+  ({ onMessage }) => {
+    const contourLineOnFace = (face, z) => {
+      let vertsAtElevation = 0;
+      let line = [];
+      for (let i = 0; i < face.length; i++) {
+        let vertex1 = face[i];
+        let vertex2 = face[(i + 1) % face.length];
+        if (vertex1[2] === z)
+          vertsAtElevation++;
+        if ((vertex1[2] <= z && vertex2[2] >= z || vertex1[2] >= z && vertex2[2] <= z) && !Number.isNaN((z - vertex1[2]) / (vertex2[2] - vertex1[2]))) {
+          let t = (z - vertex1[2]) / (vertex2[2] - vertex1[2]);
+          line.push([vertex1[0] + t * (vertex2[0] - vertex1[0]), vertex1[1] + t * (vertex2[1] - vertex1[1])]);
+        }
+      }
+      if (vertsAtElevation >= 2 && face.map((f) => f[2]).reduce((a, b) => a + b) > z * face.length)
+        return void 0;
+      if (line.length === 2 && line[0][0] === line[1][0] && line[0][1] === line[1][1])
+        return void 0;
+      if (line.length > 2) {
+        line = [...new Set(line.map((v) => JSON.stringify(v)))].map((s) => JSON.parse(s));
+      }
+      return line.length > 0 ? line : void 0;
+    };
+    const linesToPolyLines3 = (lineSegments) => {
+      var _a, _b;
+      if (!Array.isArray(lineSegments) || (lineSegments == null ? void 0 : lineSegments.length) === 0) {
+        return [];
+      }
+      const segmentsMapIndexes = {};
+      const polylines = [];
+      const parsedSegmentIndexes = [];
+      const lineSegmentStrings = lineSegments.map((v) => v.map((c) => c.join(",")));
+      lineSegmentStrings.forEach(([start, end], i) => {
+        segmentsMapIndexes[start] = segmentsMapIndexes[start] ? [...segmentsMapIndexes[start] || [], i] : [i];
+        segmentsMapIndexes[end] = segmentsMapIndexes[end] ? [...segmentsMapIndexes[end] || [], i] : [i];
+      });
+      for (let i = 0; i < lineSegmentStrings.length; i++) {
+        if (parsedSegmentIndexes.includes(i))
+          continue;
+        parsedSegmentIndexes.push(i);
+        let [start, end] = lineSegmentStrings[i];
+        let polyline = [start, end];
+        while (start && segmentsMapIndexes[start]) {
+          const nextLineIndex = (_a = segmentsMapIndexes[start]) == null ? void 0 : _a.find(
+            (lineIndex) => !parsedSegmentIndexes.includes(lineIndex)
+          );
+          if (nextLineIndex) {
+            parsedSegmentIndexes.push(nextLineIndex);
+            const nextLineSegment = lineSegmentStrings[nextLineIndex];
+            const nextLineSegmentPointIndex = nextLineSegment[0] === start ? 1 : 0;
+            const newPoint = nextLineSegment[nextLineSegmentPointIndex];
+            polyline.unshift(newPoint);
+            start = newPoint;
+          } else {
+            start = null;
+          }
+        }
+        while (end && segmentsMapIndexes[end]) {
+          const nextLineIndex = (_b = segmentsMapIndexes[end]) == null ? void 0 : _b.find(
+            (lineIndex) => !parsedSegmentIndexes.includes(lineIndex)
+          );
+          if (nextLineIndex) {
+            parsedSegmentIndexes.push(nextLineIndex);
+            const nextLineSegment = lineSegmentStrings[nextLineIndex];
+            const nextLineSegmentPointIndex = nextLineSegment[0] === end ? 1 : 0;
+            const newPoint = nextLineSegment[nextLineSegmentPointIndex];
+            polyline.push(newPoint);
+            end = newPoint;
+          } else {
+            end = null;
+          }
+        }
+        polylines.push(polyline.map((coord) => coord.split(",").map((v) => parseFloat(v))));
+      }
+      return polylines;
+    };
+    onMessage((message) => {
+      const { triangles, elevation } = message.payload;
+      const linesAtElevationE = triangles.reduce((prev, curr) => {
+        const line = contourLineOnFace(curr, elevation);
+        if (line)
+          prev.push(line);
+        return prev;
+      }, []);
+      message.resolve({
+        elevation,
+        polylines: linesToPolyLines3(linesAtElevationE)
+      });
+    });
+  },
+  { maxWorkers: 10 }
+);
 var linesToPolyLines = (lineSegments) => {
   var _a, _b;
-  if (!Array.isArray(lineSegments) || lineSegments.length === 0) {
-    throw new Error("Invalid input: Please provide a non-empty array of line segments.");
+  if (!Array.isArray(lineSegments) || (lineSegments == null ? void 0 : lineSegments.length) === 0) {
+    return [];
   }
   const segmentsMapIndexes = {};
   const polylines = [];
@@ -302,23 +471,7 @@ var getContours = (data, interval = 2) => __async(void 0, null, function* () {
     [Infinity, -Infinity]
   );
   const elevations = contourElevations(minElevation, maxElevation, interval);
-  const elevationPolylines = elevations.map((e) => {
-    const linesAtElevationE = triangles.reduce((prev, curr) => {
-      const line = contourLineOnFace(curr, e);
-      if (line)
-        prev.push(line);
-      return prev;
-    }, []);
-    const polylinesAtElevationE = linesToPolyLines(linesAtElevationE);
-    if (e === 442) {
-      console.log("linesAtElevationE", JSON.stringify(linesAtElevationE));
-      console.log("polylinesAtElevationE", JSON.stringify(polylinesAtElevationE));
-    }
-    return {
-      elevation: e,
-      polylines: polylinesAtElevationE
-    };
-  });
+  const elevationPolylines = yield Promise.all(elevations.map((elevation) => contoursWorker.send({ triangles, elevation })));
   return constructGeojson(elevationPolylines);
 });
 var get_contours_default = getContours;
@@ -373,7 +526,6 @@ var toGeojsonContours = (landXmlString, contourInterval = 2, generateOutline = t
           const geojson = yield get_contours_default(surface, contourInterval);
           if (generateOutline) {
             const outlineGeojson = get_outline_default(surface);
-            console.log(outlineGeojson.features);
             geojson.features = [...geojson.features, ...outlineGeojson.features];
           }
           const _a = surface, { surfaceDefinition } = _a, rest = __objRest(_a, ["surfaceDefinition"]);
@@ -391,7 +543,7 @@ var toGeojsonContours = (landXmlString, contourInterval = 2, generateOutline = t
 var to_geojson_contours_default = toGeojsonContours;
 
 // src/public/reproject-geojson.ts
-import proj4 from "proj4";
+var import_proj4 = __toESM(require("proj4"), 1);
 var reprojectGeoJson = (geojson, sourceProjection, targetProjection = "WGS84", keepOriginalGeometryAsFeatureProperty = true) => {
   const transformCoordinates = (coordinates, sourceProjection2, targetProjection2) => {
     if (Array.isArray(coordinates[0])) {
@@ -399,7 +551,7 @@ var reprojectGeoJson = (geojson, sourceProjection, targetProjection = "WGS84", k
         (subCoordinates) => transformCoordinates(subCoordinates, sourceProjection2, targetProjection2)
       );
     } else {
-      coordinates = proj4(sourceProjection2, targetProjection2, coordinates);
+      coordinates = (0, import_proj4.default)(sourceProjection2, targetProjection2, coordinates);
     }
     return coordinates;
   };
@@ -424,8 +576,9 @@ var reprojectGeoJson = (geojson, sourceProjection, targetProjection = "WGS84", k
   return geojson;
 };
 var reproject_geojson_default = reprojectGeoJson;
-export {
-  reproject_geojson_default as reprojectGeoJson,
-  to_geojson_contours_default as toGeojsonContours,
-  to_glb_default as toGlb
-};
+// Annotate the CommonJS export names for ESM import in node:
+0 && (module.exports = {
+  reprojectGeoJson,
+  toGeojsonContours,
+  toGlb
+});

package/dist/index.js

@@ -1,12 +1,7 @@
-"use strict";
-var __create = Object.create;
 var __defProp = Object.defineProperty;
 var __defProps = Object.defineProperties;
-var __getOwnPropDesc = Object.getOwnPropertyDescriptor;
 var __getOwnPropDescs = Object.getOwnPropertyDescriptors;
-var __getOwnPropNames = Object.getOwnPropertyNames;
 var __getOwnPropSymbols = Object.getOwnPropertySymbols;
-var __getProtoOf = Object.getPrototypeOf;
 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;
@@ -34,27 +29,6 @@ var __objRest = (source, exclude) => {
     }
   return target;
 };
-var __export = (target, all) => {
-  for (var name in all)
-    __defProp(target, name, { get: all[name], enumerable: true });
-};
-var __copyProps = (to, from, except, desc) => {
-  if (from && typeof from === "object" || typeof from === "function") {
-    for (let key of __getOwnPropNames(from))
-      if (!__hasOwnProp.call(to, key) && key !== except)
-        __defProp(to, key, { get: () => from[key], enumerable: !(desc = __getOwnPropDesc(from, key)) || desc.enumerable });
-  }
-  return to;
-};
-var __toESM = (mod, isNodeMode, target) => (target = mod != null ? __create(__getProtoOf(mod)) : {}, __copyProps(
-  // If the importer is in node compatibility mode or this is not an ESM
-  // file that has been converted to a CommonJS file using a Babel-
-  // compatible transform (i.e. "__esModule" has not been set), then set
-  // "default" to the CommonJS "module.exports" for node compatibility.
-  isNodeMode || !mod || !mod.__esModule ? __defProp(target, "default", { value: mod, enumerable: true }) : target,
-  mod
-));
-var __toCommonJS = (mod) => __copyProps(__defProp({}, "__esModule", { value: true }), mod);
 var __async = (__this, __arguments, generator) => {
   return new Promise((resolve, reject) => {
     var fulfilled = (value) => {
@@ -76,15 +50,6 @@ var __async = (__this, __arguments, generator) => {
   });
 };
 
-// src/index.ts
-var src_exports = {};
-__export(src_exports, {
-  reprojectGeoJson: () => reproject_geojson_default,
-  toGeojsonContours: () => to_geojson_contours_default,
-  toGlb: () => to_glb_default
-});
-module.exports = __toCommonJS(src_exports);
-
 // src/private/filter-by-surfaceId.ts
 var filterBySurfaceId = (parsedSurfaces, surfaceId) => {
   let filtered = [...parsedSurfaces];
@@ -103,7 +68,7 @@ var filterBySurfaceId = (parsedSurfaces, surfaceId) => {
 var filter_by_surfaceId_default = filterBySurfaceId;
 
 // src/private/get-glb.ts
-var import_core = require("@gltf-transform/core");
+import { Document, WebIO } from "@gltf-transform/core";
 var findXYAxisMedians = (vertices) => {
   var _a, _b;
   vertices = vertices.slice().filter(Boolean);
@@ -118,7 +83,7 @@ var getGlb = (data, customCenter) => __async(void 0, null, function* () {
     return [x - center[0], z, -(y - center[1])];
   }).reduce((prev, curr) => prev.concat(curr), []);
   const triangles = data.surfaceDefinition.faces.reduce((prev, curr) => prev.concat(curr), []);
-  const doc = new import_core.Document();
+  const doc = new Document();
   const buffer = doc.createBuffer();
   const position = doc.createAccessor().setType("VEC3").setArray(new Float32Array(vertices)).setBuffer(buffer);
   const indices = doc.createAccessor().setType("SCALAR").setArray(new Uint32Array(triangles)).setBuffer(buffer);
@@ -126,7 +91,7 @@ var getGlb = (data, customCenter) => __async(void 0, null, function* () {
   const mesh = doc.createMesh().addPrimitive(prim);
   const node = doc.createNode().setMesh(mesh);
   const scene = doc.createScene().addChild(node);
-  const glb = yield new import_core.WebIO().writeBinary(doc);
+  const glb = yield new WebIO().writeBinary(doc);
   return { glb, center };
 });
 var get_glb_default = getGlb;
@@ -145,47 +110,109 @@ var downloadGlb = (glbData, fileName) => {
 var download_glb_default = downloadGlb;
 
 // src/private/parse-xml.ts
-var import_xml_js = __toESM(require("xml-js"));
-var parseXML = (xmlString) => __async(void 0, null, function* () {
-  const { LandXML } = import_xml_js.default.xml2js(xmlString, {
-    compact: true,
-    attributesKey: "attr",
-    textKey: "content"
-  });
-  const landXML_surfaces = Array.isArray(LandXML.Surfaces.Surface) ? LandXML.Surfaces.Surface : [LandXML.Surfaces.Surface];
-  return landXML_surfaces.map((surface) => {
-    var _a, _b;
-    let points = [];
-    let faces = [];
-    let pointIdMap = {};
-    surface.Definition.Pnts.P.forEach((pt) => {
-      const { attr, content } = pt;
-      const [y, x, z] = content.split(" ").map((v) => parseFloat(v));
-      points.push([x, y, z]);
-      pointIdMap[attr.id] = points.length - 1;
-    }, []);
-    faces = surface.Definition.Faces.F.reduce((faceList, face) => {
-      const { attr, content } = face;
-      if (attr.i === "1")
-        return faceList;
-      const [a, b, c] = content.split(" ").map((v) => pointIdMap[v]);
-      if ([a, b, c].filter((v) => typeof v === "undefined").length > 0) {
-        throw `Invalid LandXML. A face is referencing a point that doesn't exist. Face is referencing points: ${content}`;
-      }
-      return faceList.concat([[a, b, c]]);
-    }, []);
-    return {
-      sourceFile: LandXML.Project.attr.name,
-      timeStamp: LandXML.Application.attr.timeStamp,
-      name: surface.attr.name,
-      description: surface.attr.desc,
-      wktString: (_b = (_a = LandXML.CoordinateSystem) == null ? void 0 : _a.attr) == null ? void 0 : _b.ogcWktCode,
-      surfaceDefinition: {
-        points,
-        faces
+import { createEasyWebWorker } from "easy-web-worker";
+import convert from "xml-js";
+var surfaceDefWorker = createEasyWebWorker(
+  ({ onMessage }) => {
+    onMessage((message) => {
+      const { isPoint, arr, idMap } = message.payload;
+      if (isPoint) {
+        message.resolve(
+          arr.map((pt) => [pt.attr.id, pt.content.split(" ").map(Number)]).map((v) => [v[0], [v[1][1], v[1][0], v[1][2]]])
+        );
+      } else {
+        message.resolve(
+          arr.flatMap((f) => {
+            var _a;
+            if (typeof f === "string")
+              return [f.split(" ").map((id) => idMap == null ? void 0 : idMap.indexOf(id))];
+            if (((_a = f == null ? void 0 : f.attr) == null ? void 0 : _a.i) === "1")
+              return [];
+            return [f.content.split(" ").map((id) => idMap == null ? void 0 : idMap.indexOf(id))];
+          })
+        );
       }
-    };
-  });
+    });
+  },
+  { maxWorkers: 10 }
+);
+var parseXML = (xmlString) => __async(void 0, null, function* () {
+  return new Promise((resolve, reject) => __async(void 0, null, function* () {
+    var _a, _b, _c, _d, _e;
+    const parsed = convert.xml2js(xmlString, {
+      compact: true,
+      attributesKey: "attr",
+      textKey: "content"
+    });
+    if (typeof ((_b = (_a = parsed.LandXML) == null ? void 0 : _a.Surfaces) == null ? void 0 : _b.Surface) === "undefined") {
+      throw new Error("LandXML doesn't contain any surfaces");
+      return;
+    }
+    if (!Array.isArray(parsed.LandXML.Surfaces.Surface)) {
+      parsed.LandXML.Surfaces.Surface = [parsed.LandXML.Surfaces.Surface];
+    }
+    let sourceFile = parsed.LandXML.Project.attr.name || "Undefined source";
+    let timeStamp = parsed.LandXML.Application.attr.timeStamp || "";
+    let wktString = ((_e = (_d = (_c = parsed.LandXML) == null ? void 0 : _c.CoordinateSystem) == null ? void 0 : _d.attr) == null ? void 0 : _e.ogcWktCode) || void 0;
+    const surfaces = parsed.LandXML.Surfaces.Surface.map(
+      (surface) => __async(void 0, null, function* () {
+        return new Promise((resolve2, reject2) => __async(void 0, null, function* () {
+          const { name, desc } = surface.attr;
+          const Pnts = surface.Definition.Pnts.P;
+          const Faces = surface.Definition.Faces.F;
+          let ptsIdArray = [];
+          let faces = [];
+          if (Pnts.length > 1e4) {
+            const sliceIndexes = [...Array(20).keys()].map((i) => Math.round(Pnts.length / 20 * i));
+            ptsIdArray = (yield Promise.all(
+              sliceIndexes.map(
+                (v, i, a) => new Promise((resolve3, reject3) => __async(void 0, null, function* () {
+                  const pts = yield surfaceDefWorker.send({
+                    isPoint: true,
+                    arr: Pnts.slice(a[i], a[i + 1] || Pnts.length)
+                  });
+                  resolve3(pts);
+                }))
+              )
+            )).reduce((prev, curr) => [...prev, ...curr], []);
+          } else {
+            ptsIdArray = yield surfaceDefWorker.send({ arr: Pnts, isPoint: true });
+          }
+          const points = ptsIdArray.map((v) => v[1]);
+          const pointsIdMap = ptsIdArray.map((v) => v[0]);
+          if (Faces.length > 1e4) {
+            const sliceIndexes = [...Array(20).keys()].map((i) => Math.round(Faces.length / 20 * i));
+            faces = (yield Promise.all(
+              sliceIndexes.map(
+                (v, i, a) => new Promise((resolve3, reject3) => __async(void 0, null, function* () {
+                  const fcs = yield surfaceDefWorker.send({
+                    isPoint: false,
+                    arr: Faces.slice(a[i], a[i + 1] || Faces.length),
+                    idMap: pointsIdMap
+                  });
+                  resolve3(fcs);
+                }))
+              )
+            )).reduce((prev, curr) => [...prev, ...curr], []);
+          } else {
+            faces = yield surfaceDefWorker.send({ arr: Faces, isPoint: false, idMap: pointsIdMap });
+          }
+          resolve2({
+            sourceFile,
+            timeStamp: timeStamp || "",
+            name,
+            description: desc || "",
+            wktString,
+            surfaceDefinition: {
+              points,
+              faces
+            }
+          });
+        }));
+      })
+    );
+    resolve(yield Promise.all(surfaces));
+  }));
 });
 var parse_xml_default = parseXML;
 
@@ -217,32 +244,104 @@ var toGlb = (landXmlString, center = "auto", surfaceId = -1) => __async(void 0,
 var to_glb_default = toGlb;
 
 // src/private/get-contours.ts
-var contourLineOnFace = (face, z) => {
-  let vertsAtElevation = 0;
-  let line = [];
-  for (let i = 0; i < face.length; i++) {
-    let vertex1 = face[i];
-    let vertex2 = face[(i + 1) % face.length];
-    if (vertex1[2] === z)
-      vertsAtElevation++;
-    if ((vertex1[2] <= z && vertex2[2] >= z || vertex1[2] >= z && vertex2[2] <= z) && !Number.isNaN((z - vertex1[2]) / (vertex2[2] - vertex1[2]))) {
-      let t = (z - vertex1[2]) / (vertex2[2] - vertex1[2]);
-      line.push([vertex1[0] + t * (vertex2[0] - vertex1[0]), vertex1[1] + t * (vertex2[1] - vertex1[1])]);
-    }
-  }
-  if (vertsAtElevation >= 2 && face.map((f) => f[2]).reduce((a, b) => a + b) > z * face.length)
-    return void 0;
-  if (line.length === 2 && line[0][0] === line[1][0] && line[0][1] === line[1][1])
-    return void 0;
-  if (line.length > 2) {
-    line = [...new Set(line.map((v) => JSON.stringify(v)))].map((s) => JSON.parse(s));
-  }
-  return line.length > 0 ? line : void 0;
-};
+import { createEasyWebWorker as createEasyWebWorker2 } from "easy-web-worker";
+var contoursWorker = createEasyWebWorker2(
+  ({ onMessage }) => {
+    const contourLineOnFace = (face, z) => {
+      let vertsAtElevation = 0;
+      let line = [];
+      for (let i = 0; i < face.length; i++) {
+        let vertex1 = face[i];
+        let vertex2 = face[(i + 1) % face.length];
+        if (vertex1[2] === z)
+          vertsAtElevation++;
+        if ((vertex1[2] <= z && vertex2[2] >= z || vertex1[2] >= z && vertex2[2] <= z) && !Number.isNaN((z - vertex1[2]) / (vertex2[2] - vertex1[2]))) {
+          let t = (z - vertex1[2]) / (vertex2[2] - vertex1[2]);
+          line.push([vertex1[0] + t * (vertex2[0] - vertex1[0]), vertex1[1] + t * (vertex2[1] - vertex1[1])]);
+        }
+      }
+      if (vertsAtElevation >= 2 && face.map((f) => f[2]).reduce((a, b) => a + b) > z * face.length)
+        return void 0;
+      if (line.length === 2 && line[0][0] === line[1][0] && line[0][1] === line[1][1])
+        return void 0;
+      if (line.length > 2) {
+        line = [...new Set(line.map((v) => JSON.stringify(v)))].map((s) => JSON.parse(s));
+      }
+      return line.length > 0 ? line : void 0;
+    };
+    const linesToPolyLines3 = (lineSegments) => {
+      var _a, _b;
+      if (!Array.isArray(lineSegments) || (lineSegments == null ? void 0 : lineSegments.length) === 0) {
+        return [];
+      }
+      const segmentsMapIndexes = {};
+      const polylines = [];
+      const parsedSegmentIndexes = [];
+      const lineSegmentStrings = lineSegments.map((v) => v.map((c) => c.join(",")));
+      lineSegmentStrings.forEach(([start, end], i) => {
+        segmentsMapIndexes[start] = segmentsMapIndexes[start] ? [...segmentsMapIndexes[start] || [], i] : [i];
+        segmentsMapIndexes[end] = segmentsMapIndexes[end] ? [...segmentsMapIndexes[end] || [], i] : [i];
+      });
+      for (let i = 0; i < lineSegmentStrings.length; i++) {
+        if (parsedSegmentIndexes.includes(i))
+          continue;
+        parsedSegmentIndexes.push(i);
+        let [start, end] = lineSegmentStrings[i];
+        let polyline = [start, end];
+        while (start && segmentsMapIndexes[start]) {
+          const nextLineIndex = (_a = segmentsMapIndexes[start]) == null ? void 0 : _a.find(
+            (lineIndex) => !parsedSegmentIndexes.includes(lineIndex)
+          );
+          if (nextLineIndex) {
+            parsedSegmentIndexes.push(nextLineIndex);
+            const nextLineSegment = lineSegmentStrings[nextLineIndex];
+            const nextLineSegmentPointIndex = nextLineSegment[0] === start ? 1 : 0;
+            const newPoint = nextLineSegment[nextLineSegmentPointIndex];
+            polyline.unshift(newPoint);
+            start = newPoint;
+          } else {
+            start = null;
+          }
+        }
+        while (end && segmentsMapIndexes[end]) {
+          const nextLineIndex = (_b = segmentsMapIndexes[end]) == null ? void 0 : _b.find(
+            (lineIndex) => !parsedSegmentIndexes.includes(lineIndex)
+          );
+          if (nextLineIndex) {
+            parsedSegmentIndexes.push(nextLineIndex);
+            const nextLineSegment = lineSegmentStrings[nextLineIndex];
+            const nextLineSegmentPointIndex = nextLineSegment[0] === end ? 1 : 0;
+            const newPoint = nextLineSegment[nextLineSegmentPointIndex];
+            polyline.push(newPoint);
+            end = newPoint;
+          } else {
+            end = null;
+          }
+        }
+        polylines.push(polyline.map((coord) => coord.split(",").map((v) => parseFloat(v))));
+      }
+      return polylines;
+    };
+    onMessage((message) => {
+      const { triangles, elevation } = message.payload;
+      const linesAtElevationE = triangles.reduce((prev, curr) => {
+        const line = contourLineOnFace(curr, elevation);
+        if (line)
+          prev.push(line);
+        return prev;
+      }, []);
+      message.resolve({
+        elevation,
+        polylines: linesToPolyLines3(linesAtElevationE)
+      });
+    });
+  },
+  { maxWorkers: 10 }
+);
 var linesToPolyLines = (lineSegments) => {
   var _a, _b;
-  if (!Array.isArray(lineSegments) || lineSegments.length === 0) {
-    throw new Error("Invalid input: Please provide a non-empty array of line segments.");
+  if (!Array.isArray(lineSegments) || (lineSegments == null ? void 0 : lineSegments.length) === 0) {
+    return [];
   }
   const segmentsMapIndexes = {};
   const polylines = [];
@@ -337,23 +436,7 @@ var getContours = (data, interval = 2) => __async(void 0, null, function* () {
     [Infinity, -Infinity]
   );
   const elevations = contourElevations(minElevation, maxElevation, interval);
-  const elevationPolylines = elevations.map((e) => {
-    const linesAtElevationE = triangles.reduce((prev, curr) => {
-      const line = contourLineOnFace(curr, e);
-      if (line)
-        prev.push(line);
-      return prev;
-    }, []);
-    const polylinesAtElevationE = linesToPolyLines(linesAtElevationE);
-    if (e === 442) {
-      console.log("linesAtElevationE", JSON.stringify(linesAtElevationE));
-      console.log("polylinesAtElevationE", JSON.stringify(polylinesAtElevationE));
-    }
-    return {
-      elevation: e,
-      polylines: polylinesAtElevationE
-    };
-  });
+  const elevationPolylines = yield Promise.all(elevations.map((elevation) => contoursWorker.send({ triangles, elevation })));
   return constructGeojson(elevationPolylines);
 });
 var get_contours_default = getContours;
@@ -408,7 +491,6 @@ var toGeojsonContours = (landXmlString, contourInterval = 2, generateOutline = t
           const geojson = yield get_contours_default(surface, contourInterval);
           if (generateOutline) {
             const outlineGeojson = get_outline_default(surface);
-            console.log(outlineGeojson.features);
             geojson.features = [...geojson.features, ...outlineGeojson.features];
           }
           const _a = surface, { surfaceDefinition } = _a, rest = __objRest(_a, ["surfaceDefinition"]);
@@ -426,7 +508,7 @@ var toGeojsonContours = (landXmlString, contourInterval = 2, generateOutline = t
 var to_geojson_contours_default = toGeojsonContours;
 
 // src/public/reproject-geojson.ts
-var import_proj4 = __toESM(require("proj4"));
+import proj4 from "proj4";
 var reprojectGeoJson = (geojson, sourceProjection, targetProjection = "WGS84", keepOriginalGeometryAsFeatureProperty = true) => {
   const transformCoordinates = (coordinates, sourceProjection2, targetProjection2) => {
     if (Array.isArray(coordinates[0])) {
@@ -434,7 +516,7 @@ var reprojectGeoJson = (geojson, sourceProjection, targetProjection = "WGS84", k
         (subCoordinates) => transformCoordinates(subCoordinates, sourceProjection2, targetProjection2)
       );
     } else {
-      coordinates = (0, import_proj4.default)(sourceProjection2, targetProjection2, coordinates);
+      coordinates = proj4(sourceProjection2, targetProjection2, coordinates);
     }
     return coordinates;
   };
@@ -459,9 +541,8 @@ var reprojectGeoJson = (geojson, sourceProjection, targetProjection = "WGS84", k
   return geojson;
 };
 var reproject_geojson_default = reprojectGeoJson;
-// Annotate the CommonJS export names for ESM import in node:
-0 && (module.exports = {
-  reprojectGeoJson,
-  toGeojsonContours,
-  toGlb
-});
+export {
+  reproject_geojson_default as reprojectGeoJson,
+  to_geojson_contours_default as toGeojsonContours,
+  to_glb_default as toGlb
+};

package/package.json

@@ -1,6 +1,7 @@
 {
   "name": "landxml",
-  "version": "0.5.1",
+  "type": "module",
+  "version": "0.6.0",
   "description": "Parse LandXML surfaces on the modern web.",
   "main": "dist/index.js",
   "module": "dist/index.mjs",
@@ -12,8 +13,8 @@
     "url": "https://github.com/abrman/landxml"
   },
   "scripts": {
-    "dev": "vitest",
-    "test": "vitest run",
+    "dev": "vitest --config ./vitest.config.ts --slow-test-threshold=0",
+    "test": "vitest --config ./vitest.config.ts run",
     "build": "tsup src/index.ts --format cjs,esm --dts",
     "lint": "tsc",
     "ci": "npm run lint && npm run test && npm run build",
@@ -30,14 +31,20 @@
     "@changesets/cli": "^2.26.2",
     "@types/geojson": "^7946.0.13",
     "@types/proj4": "^2.5.5",
+    "@types/sax": "^1.2.7",
     "@types/xml2json": "^0.11.6",
+    "jsdom": "^24.0.0",
     "tsup": "^8.0.0",
     "typescript": "^5.2.2",
-    "vitest": "^0.34.6"
+    "vitest": "^1.6.0"
   },
   "dependencies": {
     "@gltf-transform/core": "^3.9.0",
+    "@vitest/web-worker": "^1.6.0",
+    "easy-web-worker": "^6.2.0",
     "proj4": "^2.9.2",
-    "xml-js": "^1.6.11"
+    "sax": "^1.3.0",
+    "xml-js": "^1.6.11",
+    "xml2json": "^0.12.0"
   }
 }