landxml 0.4.2 → 0.5.1

package/CHANGELOG.md

@@ -1,5 +1,22 @@
 # landxml
 
+## 0.5.1
+
+### Patch Changes
+
+- 2d30d02: Fix: Edge case generating contour lines that are shared between multiple faces
+
+## 0.5.0
+
+### Minor Changes
+
+- df95e64: Added ability to generate surface outline geojson
+
+### Patch Changes
+
+- df95e64: LandXML Face with with "i" attribute equal to "1" will now be correctly ignored
+- df95e64: Improved contour generation algorithm
+
 ## 0.4.2
 
 ### Patch Changes

package/dist/index.d.mts

@@ -39,7 +39,7 @@ declare const toGlb: (landXmlString: string, center?: "auto" | "origin" | [
  * @returns {string} surfaceContours[].wktString - WKT string of LandXML coordinate system projection
  * @returns {Object} surfaceContours[].geojson - Geojson feature collection of contour lines
  */
-declare const toGeojsonContours: (landXmlString: string, contourInterval?: number, surfaceId?: string | number) => Promise<{
+declare const toGeojsonContours: (landXmlString: string, contourInterval?: number, generateOutline?: boolean, surfaceId?: string | number) => Promise<{
     name: string;
     description: string;
     sourceFile: string;

package/dist/index.d.ts

@@ -39,7 +39,7 @@ declare const toGlb: (landXmlString: string, center?: "auto" | "origin" | [
  * @returns {string} surfaceContours[].wktString - WKT string of LandXML coordinate system projection
  * @returns {Object} surfaceContours[].geojson - Geojson feature collection of contour lines
  */
-declare const toGeojsonContours: (landXmlString: string, contourInterval?: number, surfaceId?: string | number) => Promise<{
+declare const toGeojsonContours: (landXmlString: string, contourInterval?: number, generateOutline?: boolean, surfaceId?: string | number) => Promise<{
     name: string;
     description: string;
     sourceFile: string;

package/dist/index.js

@@ -164,14 +164,16 @@ var parseXML = (xmlString) => __async(void 0, null, function* () {
       points.push([x, y, z]);
       pointIdMap[attr.id] = points.length - 1;
     }, []);
-    faces = surface.Definition.Faces.F.map((face) => {
-      const { content } = face;
+    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 [a, b, c];
-    });
+      return faceList.concat([[a, b, c]]);
+    }, []);
     return {
       sourceFile: LandXML.Project.attr.name,
       timeStamp: LandXML.Application.attr.timeStamp,
@@ -216,15 +218,25 @@ 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 && vertex2[2] >= z || vertex1[2] >= z && vertex2[2] <= z) {
+    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 linesToPolyLines = (lineSegments) => {
@@ -325,29 +337,80 @@ var getContours = (data, interval = 2) => __async(void 0, null, function* () {
     [Infinity, -Infinity]
   );
   const elevations = contourElevations(minElevation, maxElevation, interval);
-  const elevationPolylines = elevations.map((e) => ({
-    elevation: e,
-    polylines: linesToPolyLines(
-      triangles.reduce((prev, curr) => {
-        const line = contourLineOnFace(curr, e);
-        if (line)
-          prev.push(line);
-        return prev;
-      }, [])
-    )
-  }));
+  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
+    };
+  });
   return constructGeojson(elevationPolylines);
 });
 var get_contours_default = getContours;
 
+// src/private/get-outline.ts
+var getOutline = (surface) => {
+  const triangleVertexIdEdgePairs = [];
+  let i = -1;
+  let pairs = [];
+  surface.surfaceDefinition.faces.forEach((f) => {
+    pairs = [];
+    [
+      [f[0], f[1]],
+      [f[1], f[2]],
+      [f[2], f[0]]
+    ].forEach(([a, b]) => {
+      if (a < b) {
+        pairs.push(`${a};${b}`);
+      } else {
+        pairs.push(`${b};${a}`);
+      }
+    });
+    pairs.forEach((pair) => {
+      i = triangleVertexIdEdgePairs.indexOf(pair);
+      if (i >= 0) {
+        triangleVertexIdEdgePairs.splice(i, 1);
+      } else {
+        triangleVertexIdEdgePairs.push(pair);
+      }
+    });
+  });
+  const edges = [];
+  triangleVertexIdEdgePairs.map((pair) => {
+    const [v1, v2] = pair.split(";").map((v) => {
+      var _a;
+      return (_a = surface.surfaceDefinition.points[parseInt(v, 10)]) == null ? void 0 : _a.slice(0, 2);
+    });
+    if (v1 && v2)
+      edges.push([v1, v2]);
+  });
+  return constructGeojson([{ elevation: 0, polylines: linesToPolyLines(edges) }]);
+};
+var get_outline_default = getOutline;
+
 // src/public/to-geojson-contours.ts
-var toGeojsonContours = (landXmlString, contourInterval = 2, surfaceId = -1) => __async(void 0, null, function* () {
+var toGeojsonContours = (landXmlString, contourInterval = 2, generateOutline = true, surfaceId = -1) => __async(void 0, null, function* () {
   let requestedParsedSurfaces = filter_by_surfaceId_default(yield parse_xml_default(landXmlString), surfaceId);
   const contours = yield Promise.all(
     requestedParsedSurfaces.map(
       (surface) => new Promise((resolve, reject) => __async(void 0, null, function* () {
         try {
           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"]);
           resolve(__spreadProps(__spreadValues({}, rest), {
             geojson

package/dist/index.mjs

@@ -129,14 +129,16 @@ var parseXML = (xmlString) => __async(void 0, null, function* () {
       points.push([x, y, z]);
       pointIdMap[attr.id] = points.length - 1;
     }, []);
-    faces = surface.Definition.Faces.F.map((face) => {
-      const { content } = face;
+    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 [a, b, c];
-    });
+      return faceList.concat([[a, b, c]]);
+    }, []);
     return {
       sourceFile: LandXML.Project.attr.name,
       timeStamp: LandXML.Application.attr.timeStamp,
@@ -181,15 +183,25 @@ 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 && vertex2[2] >= z || vertex1[2] >= z && vertex2[2] <= z) {
+    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 linesToPolyLines = (lineSegments) => {
@@ -290,29 +302,80 @@ var getContours = (data, interval = 2) => __async(void 0, null, function* () {
     [Infinity, -Infinity]
   );
   const elevations = contourElevations(minElevation, maxElevation, interval);
-  const elevationPolylines = elevations.map((e) => ({
-    elevation: e,
-    polylines: linesToPolyLines(
-      triangles.reduce((prev, curr) => {
-        const line = contourLineOnFace(curr, e);
-        if (line)
-          prev.push(line);
-        return prev;
-      }, [])
-    )
-  }));
+  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
+    };
+  });
   return constructGeojson(elevationPolylines);
 });
 var get_contours_default = getContours;
 
+// src/private/get-outline.ts
+var getOutline = (surface) => {
+  const triangleVertexIdEdgePairs = [];
+  let i = -1;
+  let pairs = [];
+  surface.surfaceDefinition.faces.forEach((f) => {
+    pairs = [];
+    [
+      [f[0], f[1]],
+      [f[1], f[2]],
+      [f[2], f[0]]
+    ].forEach(([a, b]) => {
+      if (a < b) {
+        pairs.push(`${a};${b}`);
+      } else {
+        pairs.push(`${b};${a}`);
+      }
+    });
+    pairs.forEach((pair) => {
+      i = triangleVertexIdEdgePairs.indexOf(pair);
+      if (i >= 0) {
+        triangleVertexIdEdgePairs.splice(i, 1);
+      } else {
+        triangleVertexIdEdgePairs.push(pair);
+      }
+    });
+  });
+  const edges = [];
+  triangleVertexIdEdgePairs.map((pair) => {
+    const [v1, v2] = pair.split(";").map((v) => {
+      var _a;
+      return (_a = surface.surfaceDefinition.points[parseInt(v, 10)]) == null ? void 0 : _a.slice(0, 2);
+    });
+    if (v1 && v2)
+      edges.push([v1, v2]);
+  });
+  return constructGeojson([{ elevation: 0, polylines: linesToPolyLines(edges) }]);
+};
+var get_outline_default = getOutline;
+
 // src/public/to-geojson-contours.ts
-var toGeojsonContours = (landXmlString, contourInterval = 2, surfaceId = -1) => __async(void 0, null, function* () {
+var toGeojsonContours = (landXmlString, contourInterval = 2, generateOutline = true, surfaceId = -1) => __async(void 0, null, function* () {
   let requestedParsedSurfaces = filter_by_surfaceId_default(yield parse_xml_default(landXmlString), surfaceId);
   const contours = yield Promise.all(
     requestedParsedSurfaces.map(
       (surface) => new Promise((resolve, reject) => __async(void 0, null, function* () {
         try {
           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"]);
           resolve(__spreadProps(__spreadValues({}, rest), {
             geojson

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "landxml",
-  "version": "0.4.2",
+  "version": "0.5.1",
   "description": "Parse LandXML surfaces on the modern web.",
   "main": "dist/index.js",
   "module": "dist/index.mjs",