landxml 0.5.2 → 0.6.1

package/CHANGELOG.md

@@ -1,5 +1,17 @@
 # landxml
 
+## 0.6.1
+
+### Patch Changes
+
+- 6c547e2: deploy options fix
+
+## 0.6.0
+
+### Minor Changes
+
+- 330e3e8: Introduced web workers for some of the processing heavy workflows
+
 ## 0.5.2
 
 ### Patch Changes

package/dist/index.js

@@ -145,47 +145,109 @@ var downloadGlb = (glbData, fileName) => {
 var download_glb_default = downloadGlb;
 
 // src/private/parse-xml.ts
+var import_easy_web_worker = require("easy-web-worker");
 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 == null ? void 0 : 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 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;
 
@@ -217,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 = [];
@@ -337,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;
@@ -408,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"]);

package/dist/index.mjs

@@ -110,47 +110,109 @@ var downloadGlb = (glbData, fileName) => {
 var download_glb_default = downloadGlb;
 
 // src/private/parse-xml.ts
+import { createEasyWebWorker } from "easy-web-worker";
 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 == null ? void 0 : 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 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;
 
@@ -182,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 = [];
@@ -302,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;
@@ -373,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"]);

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "landxml",
-  "version": "0.5.2",
+  "version": "0.6.1",
   "description": "Parse LandXML surfaces on the modern web.",
   "main": "dist/index.js",
   "module": "dist/index.mjs",
@@ -12,8 +12,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 +30,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"
   }
 }