openskidata-format 9.0.0 → 10.0.0

package/dist/ElevationProfile.d.ts

@@ -6,10 +6,12 @@ import { LineString } from 'geojson';
  * except for the final height which corresponds to the LineString endpoint and may have a different spacing.
  * Height values can be mapped to geographical coordinates using Turf.js: `turf.lineChunk(geometry, resolution, {units: 'meters'})`.
  * @property {number} resolution - The horizontal sampling distance in meters between consecutive height measurements.
+ * @property {number} targetResolution - The horizontal sampling distance in meters that was requested for the elevation profile. The actual resolution may be less to ensure an integer number of segments along the line. This is used for reconstructing the profile geometry from the feature geometry.
  */
 export type ElevationProfile = {
     heights: number[];
     resolution: number;
+    targetResolution: number;
 };
 /**
  * Elevation data that can be computed on demand fom a Run or Lift feature.
@@ -32,11 +34,33 @@ type PitchData = {
     pitchCalculationResolutionInMeters: number;
 };
 export declare function getElevationData(profileGeometry: GeoJSON.LineString): ElevationData;
-export declare function getPitchData(profileGeometry: GeoJSON.LineString, resolutionInMeters?: number): PitchData;
+export declare function getPitchData(profileGeometry: GeoJSON.LineString, minResolutionInMeters?: number): PitchData;
 export declare function getProfileGeometry(geometry: GeoJSON.LineString, elevationProfile: ElevationProfile): GeoJSON.LineString;
 /**
- * Determines the points to use for an elevation profile.
+ * Determines the points to use for an elevation profile, given a minimum horizontal resolution. The resolution is selected to have evenly spaced points along the line.
  */
-export declare function extractPointsForElevationProfile(geometry: LineString, resolution: number): GeoJSON.LineString;
+export declare function extractPointsForElevationProfile(geometry: LineString, minResolutionInMeters: number): {
+    resolutionInMeters: number;
+    geometry: GeoJSON.LineString;
+};
 export declare function getAscentAndDescent(profileGeometry: GeoJSON.LineString): AscentDescentData;
+/**
+ * Turf's lineChunk has several issues:
+ * - it doesn't allow for a custom resolution that divides the line evenly, which is important for consistent pitch calculations
+ *   change: calculate a resolution that divides the line into an integer number of segments, with a maximum of the provided resolution
+ * - it includes elevations from the original geometry on newly created points which are wrong (https://github.com/Turfjs/turf/issues/3007)
+ *   fix: strip Z coordinates from non-original points
+ * - it can produce a last point that is very close but not exactly the same as the original endpoint
+ *   fix: set the last point to be exactly the original endpoint
+ * - it can produce two points very close to the end of line due to floating point precision issues
+ *   fix: drop the second to last point if it's very close to the end
+ *
+ * @param geometry
+ * @param minResolutionInMeters
+ * @returns
+ */
+export declare function lineChunkPatched(geometry: LineString, minResolutionInMeters: number): {
+    resolutionInMeters: number;
+    geometry: GeoJSON.LineString[];
+};
 export {};

package/dist/ElevationProfile.js

@@ -5,6 +5,7 @@ exports.getPitchData = getPitchData;
 exports.getProfileGeometry = getProfileGeometry;
 exports.extractPointsForElevationProfile = extractPointsForElevationProfile;
 exports.getAscentAndDescent = getAscentAndDescent;
+exports.lineChunkPatched = lineChunkPatched;
 const tslib_1 = require("tslib");
 const distance_1 = tslib_1.__importDefault(require("@turf/distance"));
 const length_1 = tslib_1.__importDefault(require("@turf/length"));
@@ -16,7 +17,7 @@ function getElevationData(profileGeometry) {
         profileGeometry,
     };
 }
-function getPitchData(profileGeometry, resolutionInMeters = 25) {
+function getPitchData(profileGeometry, minResolutionInMeters = 25) {
     const coordinates = profileGeometry.coordinates;
     if (coordinates[0].length < 3) {
         throw 'Elevation data is required for slope analysis';
@@ -47,7 +48,7 @@ function getPitchData(profileGeometry, resolutionInMeters = 25) {
     // If the line is too short relative to the resolution, pitch calculations are unreliable
     // due to elevation data resolution. A 1m elevation change over 1m distance would be a 45 degree slope,
     // which can easily happen due to elevation data precision issues.
-    if (totalLength < resolutionInMeters / 2) {
+    if (totalLength < minResolutionInMeters / 2) {
         return {
             averagePitchInPercent: null,
             maxPitchInPercent: null,
@@ -56,18 +57,8 @@ function getPitchData(profileGeometry, resolutionInMeters = 25) {
             pitchCalculationResolutionInMeters: totalLength,
         };
     }
-    // Calculate a resolution that divides the distance evenly, with the constraint that it must be <= resolutionInMeters
-    const numSegments = Math.ceil(totalLength / resolutionInMeters);
-    const actualResolution = totalLength / numSegments;
     // Chunk the line into segments using the calculated resolution for max pitch calculation
-    const chunkedGeometries = (0, line_chunk_1.lineChunk)(profileGeometry, actualResolution, {
-        units: 'meters',
-    }).features.map((feature) => feature.geometry);
-    // turf 7.3+ adds Z coordinates to newly created points using segment-end elevation,
-    // but we need distance-based interpolation. Strip Z from non-original vertices so
-    // interpolateElevation can add them correctly.
-    // https://github.com/Turfjs/turf/issues/3007
-    stripNonOriginalElevations(chunkedGeometries, coordinates);
+    const { geometry: chunkedGeometries, resolutionInMeters: actualResolutionInMeters, } = lineChunkPatched(profileGeometry, minResolutionInMeters);
     // Mutates the geometries in place to add elevation data
     interpolateElevation(chunkedGeometries);
     // Calculate max pitch over the evenly-divided chunks
@@ -90,14 +81,18 @@ function getPitchData(profileGeometry, resolutionInMeters = 25) {
         maxPitchInPercent: maxPitchValue,
         inclinedLengthInMeters: totalInclinedLength,
         overallPitchInPercent: overallElevationChange / totalLength,
-        pitchCalculationResolutionInMeters: actualResolution,
+        pitchCalculationResolutionInMeters: actualResolutionInMeters,
     };
 }
 function getProfileGeometry(geometry, elevationProfile) {
-    const profileLine = extractPointsForElevationProfile(geometry, elevationProfile.resolution);
+    const { geometry: profileLine, resolutionInMeters: actualResolution } = extractPointsForElevationProfile(geometry, elevationProfile.targetResolution);
     if (profileLine.coordinates.length !== elevationProfile.heights.length) {
         throw `Mismatch of points & elevation profile`;
     }
+    if (Math.abs(actualResolution - elevationProfile.resolution) >
+        elevationProfile.resolution * 0.001) {
+        throw `Resolution mismatch between profile geometry (${actualResolution}) and elevation profile (${elevationProfile.resolution})`;
+    }
     for (let i = 0; i < profileLine.coordinates.length; i++) {
         const point = profileLine.coordinates[i];
         const height = elevationProfile.heights[i];
@@ -106,13 +101,10 @@ function getProfileGeometry(geometry, elevationProfile) {
     return profileLine;
 }
 /**
- * Determines the points to use for an elevation profile.
+ * Determines the points to use for an elevation profile, given a minimum horizontal resolution. The resolution is selected to have evenly spaced points along the line.
  */
-function extractPointsForElevationProfile(geometry, resolution) {
-    // Important: Z coordinates in the output of lineChunk may be incorrect, ignore them (https://github.com/Turfjs/turf/issues/3007)
-    const lineChunks = (0, line_chunk_1.lineChunk)(geometry, resolution, {
-        units: 'meters',
-    }).features.map((feature) => feature.geometry);
+function extractPointsForElevationProfile(geometry, minResolutionInMeters) {
+    const { geometry: lineChunks, resolutionInMeters } = lineChunkPatched(geometry, minResolutionInMeters);
     const points = [];
     for (let subline of lineChunks) {
         const point = subline.coordinates[0];
@@ -127,8 +119,11 @@ function extractPointsForElevationProfile(geometry, resolution) {
         }
     }
     return {
-        type: 'LineString',
-        coordinates: points,
+        resolutionInMeters,
+        geometry: {
+            type: 'LineString',
+            coordinates: points,
+        },
     };
 }
 function getAscentAndDescent(profileGeometry) {
@@ -169,6 +164,53 @@ function getAscentAndDescent(profileGeometry) {
         verticalInMeters: result.maxElevation - result.minElevation,
     };
 }
+/**
+ * Turf's lineChunk has several issues:
+ * - it doesn't allow for a custom resolution that divides the line evenly, which is important for consistent pitch calculations
+ *   change: calculate a resolution that divides the line into an integer number of segments, with a maximum of the provided resolution
+ * - it includes elevations from the original geometry on newly created points which are wrong (https://github.com/Turfjs/turf/issues/3007)
+ *   fix: strip Z coordinates from non-original points
+ * - it can produce a last point that is very close but not exactly the same as the original endpoint
+ *   fix: set the last point to be exactly the original endpoint
+ * - it can produce two points very close to the end of line due to floating point precision issues
+ *   fix: drop the second to last point if it's very close to the end
+ *
+ * @param geometry
+ * @param minResolutionInMeters
+ * @returns
+ */
+function lineChunkPatched(geometry, minResolutionInMeters) {
+    const resolutionInMeters = horizontalResolutionFitting(geometry, minResolutionInMeters);
+    const lineChunks = (0, line_chunk_1.lineChunk)(geometry, resolutionInMeters, {
+        units: 'meters',
+    }).features.map((feature) => feature.geometry);
+    // Check the last segment has the expected length
+    const lastSegment = lineChunks[lineChunks.length - 1];
+    const lastSegmentLength = (0, length_1.default)({ type: 'Feature', geometry: lastSegment, properties: {} }, { units: 'meters' });
+    if (lastSegmentLength < resolutionInMeters * 0.01) {
+        // If the last segment is very short, which can happen due to floating point precision issues, drop it to avoid issues with elevation interpolation
+        lineChunks.pop();
+    }
+    // Set the last point from the original coordinates, turf line chunk can sometimes produce a last point that
+    // is very close but not exactly the same as the original endpoint, which can cause issues with elevation interpolation
+    const lastOriginalPoint = geometry.coordinates[geometry.coordinates.length - 1];
+    const lastChunkCoords = lineChunks[lineChunks.length - 1].coordinates;
+    lastChunkCoords[lastChunkCoords.length - 1] = lastOriginalPoint;
+    // turf 7.3+ adds Z coordinates to newly created points using segment-end elevation,
+    // but we need distance-based interpolation. Strip Z from non-original vertices so
+    // interpolateElevation can add them correctly.
+    // https://github.com/Turfjs/turf/issues/3007
+    stripNonOriginalElevations(lineChunks, geometry.coordinates);
+    return { resolutionInMeters, geometry: lineChunks };
+}
+function horizontalResolutionFitting(geometry, minResolution) {
+    const totalLength = (0, length_1.default)({ type: 'Feature', geometry, properties: {} }, { units: 'meters' });
+    if (totalLength === 0) {
+        return 0;
+    }
+    const numSegments = Math.ceil(totalLength / minResolution);
+    return totalLength / numSegments;
+}
 /**
  * Strip Z coordinates from points that were created by lineChunk (not original vertices).
  * Original vertices are identified by matching lon/lat coordinates.
@@ -224,7 +266,7 @@ function interpolateElevation(geometries) {
     // Find elevation reference points
     const referencePoints = allPoints.filter((p) => p.hasElevation);
     if (referencePoints.length <= 1) {
-        throw 'At least two points with elevation data are required';
+        throw `At least two points with elevation data are required ${JSON.stringify(allPoints)}`;
     }
     // Process segments between each pair of reference points
     for (let i = 0; i < referencePoints.length - 1; i++) {

package/dist/ElevationProfile.js.map

@@ -1 +1 @@
-{"version":3,"file":"ElevationProfile.js","sourceRoot":"","sources":["../src/ElevationProfile.ts"],"names":[],"mappings":";;AA0CA,4CAQC;AAED,oCAoGC;AAED,gDAmBC;AAKD,4EA2BC;AAED,kDAiDC;;AAhQD,sEAAqC;AACrC,kEAAiC;AACjC,iDAA4C;AAwC5C,SAAgB,gBAAgB,CAC9B,eAAmC;IAEnC,OAAO;QACL,GAAG,mBAAmB,CAAC,eAAe,CAAC;QACvC,GAAG,YAAY,CAAC,eAAe,CAAC;QAChC,eAAe;KAChB,CAAA;AACH,CAAC;AAED,SAAgB,YAAY,CAC1B,eAAmC,EACnC,qBAA6B,EAAE;IAE/B,MAAM,WAAW,GAAG,eAAe,CAAC,WAAW,CAAA;IAC/C,IAAI,WAAW,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QAC9B,MAAM,+CAA+C,CAAA;IACvD,CAAC;IAED,+DAA+D;IAC/D,IAAI,WAAW,GAAG,CAAC,CAAA;IACnB,IAAI,mBAAmB,GAAG,CAAC,CAAA;IAC3B,IAAI,oBAAoB,GAAG,CAAC,CAAA;IAE5B,2DAA2D;IAC3D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,WAAW,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAChD,MAAM,MAAM,GAAG,WAAW,CAAC,CAAC,CAAC,CAAA;QAC7B,MAAM,KAAK,GAAG,WAAW,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;QAChC,MAAM,SAAS,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAA;QACtC,MAAM,QAAQ,GAAuB;YACnC,IAAI,EAAE,YAAY;YAClB,WAAW,EAAE,CAAC,MAAM,EAAE,KAAK,CAAC;SAC7B,CAAA;QACD,MAAM,OAAO,GAAoB;YAC/B,IAAI,EAAE,SAAS;YACf,QAAQ;YACR,UAAU,EAAE,EAAE;SACf,CAAA;QACD,MAAM,cAAc,GAAG,IAAA,gBAAM,EAAC,OAAO,EAAE,EAAE,KAAK,EAAE,QAAQ,EAAE,CAAC,CAAA;QAE3D,WAAW,IAAI,cAAc,CAAA;QAC7B,mBAAmB,IAAI,IAAI,CAAC,IAAI,CAC9B,IAAI,CAAC,GAAG,CAAC,cAAc,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,EAAE,CAAC,CAAC,CACrD,CAAA;QACD,oBAAoB,IAAI,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAA;IAC7C,CAAC;IAED,yFAAyF;IACzF,uGAAuG;IACvG,kEAAkE;IAClE,IAAI,WAAW,GAAG,kBAAkB,GAAG,CAAC,EAAE,CAAC;QACzC,OAAO;YACL,qBAAqB,EAAE,IAAI;YAC3B,iBAAiB,EAAE,IAAI;YACvB,sBAAsB,EAAE,mBAAmB;YAC3C,qBAAqB,EAAE,IAAI;YAC3B,kCAAkC,EAAE,WAAW;SAChD,CAAA;IACH,CAAC;IAED,qHAAqH;IACrH,MAAM,WAAW,GAAG,IAAI,CAAC,IAAI,CAAC,WAAW,GAAG,kBAAkB,CAAC,CAAA;IAC/D,MAAM,gBAAgB,GAAG,WAAW,GAAG,WAAW,CAAA;IAElD,yFAAyF;IACzF,MAAM,iBAAiB,GAAG,IAAA,sBAAS,EAAC,eAAe,EAAE,gBAAgB,EAAE;QACrE,KAAK,EAAE,QAAQ;KAChB,CAAC,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,OAAO,EAAE,EAAE,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAA;IAE9C,oFAAoF;IACpF,kFAAkF;IAClF,+CAA+C;IAC/C,6CAA6C;IAC7C,0BAA0B,CAAC,iBAAiB,EAAE,WAAW,CAAC,CAAA;IAE1D,wDAAwD;IACxD,oBAAoB,CAAC,iBAAiB,CAAC,CAAA;IAEvC,qDAAqD;IACrD,IAAI,aAAa,GAAG,CAAC,CAAA;IACrB,KAAK,MAAM,KAAK,IAAI,iBAAiB,EAAE,CAAC;QACtC,MAAM,WAAW,GAAG,KAAK,CAAC,WAAW,CAAA;QAErC,MAAM,UAAU,GAAG,WAAW,CAAC,CAAC,CAAC,CAAA;QACjC,MAAM,QAAQ,GAAG,WAAW,CAAC,WAAW,CAAC,MAAM,GAAG,CAAC,CAAC,CAAA;QAEpD,MAAM,eAAe,GAAG,QAAQ,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAA;QACnD,MAAM,mBAAmB,GAAG,IAAA,gBAAM,EAChC,EAAE,IAAI,EAAE,SAAS,EAAE,QAAQ,EAAE,KAAK,EAAE,UAAU,EAAE,EAAE,EAAE,EACpD,EAAE,KAAK,EAAE,QAAQ,EAAE,CACpB,CAAA;QAED,MAAM,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,eAAe,GAAG,mBAAmB,CAAC,CAAA;QAClE,IAAI,UAAU,GAAG,aAAa,EAAE,CAAC;YAC/B,aAAa,GAAG,UAAU,CAAA;QAC5B,CAAC;IACH,CAAC;IAED,MAAM,YAAY,GAAG,oBAAoB,GAAG,WAAW,CAAA;IACvD,MAAM,sBAAsB,GAAG,IAAI,CAAC,GAAG,CACrC,WAAW,CAAC,WAAW,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAC3D,CAAA;IAED,OAAO;QACL,qBAAqB,EAAE,YAAY;QACnC,iBAAiB,EAAE,aAAa;QAChC,sBAAsB,EAAE,mBAAmB;QAC3C,qBAAqB,EAAE,sBAAsB,GAAG,WAAW;QAC3D,kCAAkC,EAAE,gBAAgB;KACrD,CAAA;AACH,CAAC;AAED,SAAgB,kBAAkB,CAChC,QAA4B,EAC5B,gBAAkC;IAElC,MAAM,WAAW,GAAG,gCAAgC,CAClD,QAAQ,EACR,gBAAgB,CAAC,UAAU,CAC5B,CAAA;IACD,IAAI,WAAW,CAAC,WAAW,CAAC,MAAM,KAAK,gBAAgB,CAAC,OAAO,CAAC,MAAM,EAAE,CAAC;QACvE,MAAM,wCAAwC,CAAA;IAChD,CAAC;IAED,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,WAAW,CAAC,WAAW,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACxD,MAAM,KAAK,GAAG,WAAW,CAAC,WAAW,CAAC,CAAC,CAAC,CAAA;QACxC,MAAM,MAAM,GAAG,gBAAgB,CAAC,OAAO,CAAC,CAAC,CAAC,CAAA;QAC1C,KAAK,CAAC,IAAI,CAAC,MAAM,CAAC,CAAA;IACpB,CAAC;IAED,OAAO,WAAW,CAAA;AACpB,CAAC;AAED;;GAEG;AACH,SAAgB,gCAAgC,CAC9C,QAAoB,EACpB,UAAkB;IAElB,iIAAiI;IACjI,MAAM,UAAU,GAAG,IAAA,sBAAS,EAAC,QAAQ,EAAE,UAAU,EAAE;QACjD,KAAK,EAAE,QAAQ;KAChB,CAAC,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,OAAO,EAAE,EAAE,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAA;IAE9C,MAAM,MAAM,GAAuB,EAAE,CAAA;IACrC,KAAK,IAAI,OAAO,IAAI,UAAU,EAAE,CAAC;QAC/B,MAAM,KAAK,GAAG,OAAO,CAAC,WAAW,CAAC,CAAC,CAAC,CAAA;QACpC,MAAM,CAAC,IAAI,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;IACnC,CAAC;IACD,IAAI,UAAU,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QAC1B,MAAM,QAAQ,GAAG,UAAU,CAAC,UAAU,CAAC,MAAM,GAAG,CAAC,CAAC,CAAA;QAClD,MAAM,MAAM,GAAG,QAAQ,CAAC,WAAW,CAAA;QACnC,IAAI,MAAM,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;YACtB,MAAM,KAAK,GAAG,MAAM,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,CAAA;YACvC,MAAM,CAAC,IAAI,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;QACnC,CAAC;IACH,CAAC;IAED,OAAO;QACL,IAAI,EAAE,YAAY;QAClB,WAAW,EAAE,MAAM;KACpB,CAAA;AACH,CAAC;AAED,SAAgB,mBAAmB,CACjC,eAAmC;IAEnC,MAAM,WAAW,GAAG,eAAe,CAAC,WAAW,CAAA;IAC/C,IAAI,WAAW,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QAC7B,MAAM,iEAAiE,CAAA;IACzE,CAAC;IACD,IAAI,WAAW,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QAC9B,MAAM,wDAAwD,CAAA;IAChE,CAAC;IACD,MAAM,gBAAgB,GAAG,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;IAC1C,MAAM,MAAM,GAAG,WAAW;SACvB,GAAG,CAAC,CAAC,KAAK,EAAE,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;SACxB,MAAM,CACL,CAAC,WAAW,EAAE,gBAAgB,EAAE,EAAE;QAChC,MAAM,MAAM,GAAG,gBAAgB,GAAG,WAAW,CAAC,aAAa,CAAA;QAC3D,IAAI,MAAM,GAAG,CAAC,EAAE,CAAC;YACf,WAAW,CAAC,MAAM,IAAI,MAAM,CAAA;QAC9B,CAAC;aAAM,CAAC;YACN,WAAW,CAAC,OAAO,IAAI,MAAM,CAAA;QAC/B,CAAC;QACD,WAAW,CAAC,aAAa,GAAG,gBAAgB,CAAA;QAC5C,WAAW,CAAC,YAAY,GAAG,IAAI,CAAC,GAAG,CACjC,gBAAgB,EAChB,WAAW,CAAC,YAAY,CACzB,CAAA;QACD,WAAW,CAAC,YAAY,GAAG,IAAI,CAAC,GAAG,CACjC,gBAAgB,EAChB,WAAW,CAAC,YAAY,CACzB,CAAA;QAED,OAAO,WAAW,CAAA;IACpB,CAAC,EACD;QACE,MAAM,EAAE,CAAC;QACT,OAAO,EAAE,CAAC;QACV,YAAY,EAAE,gBAAgB;QAC9B,YAAY,EAAE,gBAAgB;QAC9B,aAAa,EAAE,gBAAgB;KAChC,CACF,CAAA;IAEH,OAAO;QACL,cAAc,EAAE,MAAM,CAAC,MAAM;QAC7B,eAAe,EAAE,MAAM,CAAC,OAAO;QAC/B,oBAAoB,EAAE,MAAM,CAAC,YAAY;QACzC,oBAAoB,EAAE,MAAM,CAAC,YAAY;QACzC,gBAAgB,EAAE,MAAM,CAAC,YAAY,GAAG,MAAM,CAAC,YAAY;KAC5D,CAAA;AACH,CAAC;AAED;;;;GAIG;AACH,SAAS,0BAA0B,CACjC,UAAgC,EAChC,mBAAuC;IAEvC,MAAM,WAAW,GAAG,IAAI,GAAG,CAAC,mBAAmB,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAA;IAE9E,KAAK,MAAM,QAAQ,IAAI,UAAU,EAAE,CAAC;QAClC,KAAK,MAAM,KAAK,IAAI,QAAQ,CAAC,WAAW,EAAE,CAAC;YACzC,MAAM,GAAG,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC,CAAC,CAAC,EAAE,CAAA;YACrC,IAAI,CAAC,WAAW,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,KAAK,CAAC,MAAM,IAAI,CAAC,EAAE,CAAC;gBAC/C,KAAK,CAAC,MAAM,GAAG,CAAC,CAAA;YAClB,CAAC;QACH,CAAC;IACH,CAAC;AACH,CAAC;AAED;;;;;;GAMG;AACH,SAAS,oBAAoB,CAAC,UAAgC;IAC5D,IAAI,UAAU,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QAC5B,OAAM;IACR,CAAC;IAED,oEAAoE;IACpE,MAAM,SAAS,GAKT,EAAE,CAAA;IAER,8DAA8D;IAC9D,IAAI,UAAU,GAAG,CAAC,CAAA;IAClB,IAAI,aAAa,GAAG,CAAC,CAAA;IACrB,IAAI,aAAa,GAA4B,IAAI,CAAA;IAEjD,KAAK,MAAM,QAAQ,IAAI,UAAU,EAAE,CAAC;QAClC,KAAK,MAAM,KAAK,IAAI,QAAQ,CAAC,WAAW,EAAE,CAAC;YACzC,yCAAyC;YACzC,IAAI,aAAa,EAAE,CAAC;gBAClB,MAAM,aAAa,GAAG,IAAA,kBAAQ,EAAC,aAAa,EAAE,KAAK,EAAE;oBACnD,KAAK,EAAE,QAAQ;iBAChB,CAAC,CAAA;gBACF,aAAa,IAAI,aAAa,CAAA;YAChC,CAAC;YAED,MAAM,YAAY,GAAG,KAAK,CAAC,MAAM,IAAI,CAAC,CAAA;YACtC,SAAS,CAAC,IAAI,CAAC;gBACb,KAAK;gBACL,YAAY;gBACZ,KAAK,EAAE,UAAU,EAAE;gBACnB,iBAAiB,EAAE,aAAa;aACjC,CAAC,CAAA;YAEF,aAAa,GAAG,KAAK,CAAA;QACvB,CAAC;IACH,CAAC;IAED,kCAAkC;IAClC,MAAM,eAAe,GAAG,SAAS,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,YAAY,CAAC,CAAA;IAE/D,IAAI,eAAe,CAAC,MAAM,IAAI,CAAC,EAAE,CAAC;QAChC,MAAM,sDAAsD,CAAA;IAC9D,CAAC;IAED,yDAAyD;IACzD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,eAAe,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QACpD,MAAM,UAAU,GAAG,eAAe,CAAC,CAAC,CAAC,CAAA;QACrC,MAAM,QAAQ,GAAG,eAAe,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;QAEvC,6BAA6B;QAC7B,IAAI,QAAQ,CAAC,KAAK,IAAI,UAAU,CAAC,KAAK,EAAE,CAAC;YACvC,MAAM,sCAAsC,CAAA;QAC9C,CAAC;QAED,MAAM,cAAc,GAAG,UAAU,CAAC,KAAK,CAAC,CAAC,CAAC,CAAA;QAC1C,MAAM,YAAY,GAAG,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAA;QACtC,MAAM,cAAc,GAAG,YAAY,GAAG,cAAc,CAAA;QAEpD,6DAA6D;QAC7D,MAAM,mBAAmB,GACvB,QAAQ,CAAC,iBAAiB,GAAG,UAAU,CAAC,iBAAiB,CAAA;QAE3D,uDAAuD;QACvD,KAAK,IAAI,CAAC,GAAG,UAAU,CAAC,KAAK,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3D,MAAM,SAAS,GAAG,SAAS,CAAC,CAAC,CAAC,CAAA;YAE9B,2DAA2D;YAC3D,IAAI,SAAS,CAAC,YAAY,EAAE,CAAC;gBAC3B,MAAM,6DAA6D,CAAA;YACrE,CAAC;YACD,qEAAqE;YACrE,MAAM,iBAAiB,GACrB,SAAS,CAAC,iBAAiB,GAAG,UAAU,CAAC,iBAAiB,CAAA;YAE5D,0FAA0F;YAC1F,IAAI,qBAAqB,GAAG,cAAc,CAAA;YAC1C,IAAI,mBAAmB,GAAG,CAAC,EAAE,CAAC;gBAC5B,MAAM,mBAAmB,GAAG,iBAAiB,GAAG,mBAAmB,CAAA;gBACnE,qBAAqB;oBACnB,cAAc,GAAG,cAAc,GAAG,mBAAmB,CAAA;YACzD,CAAC;YAED,wEAAwE;YACxE,SAAS,CAAC,KAAK,CAAC,IAAI,CAAC,qBAAqB,CAAC,CAAA;QAC7C,CAAC;IACH,CAAC;AACH,CAAC"}
+{"version":3,"file":"ElevationProfile.js","sourceRoot":"","sources":["../src/ElevationProfile.ts"],"names":[],"mappings":";;AA4CA,4CAQC;AAED,oCA2FC;AAED,gDA2BC;AAKD,4EA8BC;AAED,kDAiDC;AAiBD,4CAsCC;;AA3TD,sEAAqC;AACrC,kEAAiC;AACjC,iDAA4C;AA0C5C,SAAgB,gBAAgB,CAC9B,eAAmC;IAEnC,OAAO;QACL,GAAG,mBAAmB,CAAC,eAAe,CAAC;QACvC,GAAG,YAAY,CAAC,eAAe,CAAC;QAChC,eAAe;KAChB,CAAA;AACH,CAAC;AAED,SAAgB,YAAY,CAC1B,eAAmC,EACnC,wBAAgC,EAAE;IAElC,MAAM,WAAW,GAAG,eAAe,CAAC,WAAW,CAAA;IAC/C,IAAI,WAAW,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QAC9B,MAAM,+CAA+C,CAAA;IACvD,CAAC;IAED,+DAA+D;IAC/D,IAAI,WAAW,GAAG,CAAC,CAAA;IACnB,IAAI,mBAAmB,GAAG,CAAC,CAAA;IAC3B,IAAI,oBAAoB,GAAG,CAAC,CAAA;IAE5B,2DAA2D;IAC3D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,WAAW,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QAChD,MAAM,MAAM,GAAG,WAAW,CAAC,CAAC,CAAC,CAAA;QAC7B,MAAM,KAAK,GAAG,WAAW,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;QAChC,MAAM,SAAS,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAA;QACtC,MAAM,QAAQ,GAAuB;YACnC,IAAI,EAAE,YAAY;YAClB,WAAW,EAAE,CAAC,MAAM,EAAE,KAAK,CAAC;SAC7B,CAAA;QACD,MAAM,OAAO,GAAoB;YAC/B,IAAI,EAAE,SAAS;YACf,QAAQ;YACR,UAAU,EAAE,EAAE;SACf,CAAA;QACD,MAAM,cAAc,GAAG,IAAA,gBAAM,EAAC,OAAO,EAAE,EAAE,KAAK,EAAE,QAAQ,EAAE,CAAC,CAAA;QAE3D,WAAW,IAAI,cAAc,CAAA;QAC7B,mBAAmB,IAAI,IAAI,CAAC,IAAI,CAC9B,IAAI,CAAC,GAAG,CAAC,cAAc,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,EAAE,CAAC,CAAC,CACrD,CAAA;QACD,oBAAoB,IAAI,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAA;IAC7C,CAAC;IAED,yFAAyF;IACzF,uGAAuG;IACvG,kEAAkE;IAClE,IAAI,WAAW,GAAG,qBAAqB,GAAG,CAAC,EAAE,CAAC;QAC5C,OAAO;YACL,qBAAqB,EAAE,IAAI;YAC3B,iBAAiB,EAAE,IAAI;YACvB,sBAAsB,EAAE,mBAAmB;YAC3C,qBAAqB,EAAE,IAAI;YAC3B,kCAAkC,EAAE,WAAW;SAChD,CAAA;IACH,CAAC;IAED,yFAAyF;IACzF,MAAM,EACJ,QAAQ,EAAE,iBAAiB,EAC3B,kBAAkB,EAAE,wBAAwB,GAC7C,GAAG,gBAAgB,CAAC,eAAe,EAAE,qBAAqB,CAAC,CAAA;IAE5D,wDAAwD;IACxD,oBAAoB,CAAC,iBAAiB,CAAC,CAAA;IAEvC,qDAAqD;IACrD,IAAI,aAAa,GAAG,CAAC,CAAA;IACrB,KAAK,MAAM,KAAK,IAAI,iBAAiB,EAAE,CAAC;QACtC,MAAM,WAAW,GAAG,KAAK,CAAC,WAAW,CAAA;QAErC,MAAM,UAAU,GAAG,WAAW,CAAC,CAAC,CAAC,CAAA;QACjC,MAAM,QAAQ,GAAG,WAAW,CAAC,WAAW,CAAC,MAAM,GAAG,CAAC,CAAC,CAAA;QAEpD,MAAM,eAAe,GAAG,QAAQ,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,CAAC,CAAC,CAAA;QACnD,MAAM,mBAAmB,GAAG,IAAA,gBAAM,EAChC,EAAE,IAAI,EAAE,SAAS,EAAE,QAAQ,EAAE,KAAK,EAAE,UAAU,EAAE,EAAE,EAAE,EACpD,EAAE,KAAK,EAAE,QAAQ,EAAE,CACpB,CAAA;QAED,MAAM,UAAU,GAAG,IAAI,CAAC,GAAG,CAAC,eAAe,GAAG,mBAAmB,CAAC,CAAA;QAClE,IAAI,UAAU,GAAG,aAAa,EAAE,CAAC;YAC/B,aAAa,GAAG,UAAU,CAAA;QAC5B,CAAC;IACH,CAAC;IAED,MAAM,YAAY,GAAG,oBAAoB,GAAG,WAAW,CAAA;IACvD,MAAM,sBAAsB,GAAG,IAAI,CAAC,GAAG,CACrC,WAAW,CAAC,WAAW,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAC3D,CAAA;IAED,OAAO;QACL,qBAAqB,EAAE,YAAY;QACnC,iBAAiB,EAAE,aAAa;QAChC,sBAAsB,EAAE,mBAAmB;QAC3C,qBAAqB,EAAE,sBAAsB,GAAG,WAAW;QAC3D,kCAAkC,EAAE,wBAAwB;KAC7D,CAAA;AACH,CAAC;AAED,SAAgB,kBAAkB,CAChC,QAA4B,EAC5B,gBAAkC;IAElC,MAAM,EAAE,QAAQ,EAAE,WAAW,EAAE,kBAAkB,EAAE,gBAAgB,EAAE,GACnE,gCAAgC,CAC9B,QAAQ,EACR,gBAAgB,CAAC,gBAAgB,CAClC,CAAA;IACH,IAAI,WAAW,CAAC,WAAW,CAAC,MAAM,KAAK,gBAAgB,CAAC,OAAO,CAAC,MAAM,EAAE,CAAC;QACvE,MAAM,wCAAwC,CAAA;IAChD,CAAC;IAED,IACE,IAAI,CAAC,GAAG,CAAC,gBAAgB,GAAG,gBAAgB,CAAC,UAAU,CAAC;QACxD,gBAAgB,CAAC,UAAU,GAAG,KAAK,EACnC,CAAC;QACD,MAAM,iDAAiD,gBAAgB,4BAA4B,gBAAgB,CAAC,UAAU,GAAG,CAAA;IACnI,CAAC;IAED,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,WAAW,CAAC,WAAW,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;QACxD,MAAM,KAAK,GAAG,WAAW,CAAC,WAAW,CAAC,CAAC,CAAC,CAAA;QACxC,MAAM,MAAM,GAAG,gBAAgB,CAAC,OAAO,CAAC,CAAC,CAAC,CAAA;QAC1C,KAAK,CAAC,IAAI,CAAC,MAAM,CAAC,CAAA;IACpB,CAAC;IAED,OAAO,WAAW,CAAA;AACpB,CAAC;AAED;;GAEG;AACH,SAAgB,gCAAgC,CAC9C,QAAoB,EACpB,qBAA6B;IAE7B,MAAM,EAAE,QAAQ,EAAE,UAAU,EAAE,kBAAkB,EAAE,GAAG,gBAAgB,CACnE,QAAQ,EACR,qBAAqB,CACtB,CAAA;IAED,MAAM,MAAM,GAAuB,EAAE,CAAA;IACrC,KAAK,IAAI,OAAO,IAAI,UAAU,EAAE,CAAC;QAC/B,MAAM,KAAK,GAAG,OAAO,CAAC,WAAW,CAAC,CAAC,CAAC,CAAA;QACpC,MAAM,CAAC,IAAI,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;IACnC,CAAC;IACD,IAAI,UAAU,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QAC1B,MAAM,QAAQ,GAAG,UAAU,CAAC,UAAU,CAAC,MAAM,GAAG,CAAC,CAAC,CAAA;QAClD,MAAM,MAAM,GAAG,QAAQ,CAAC,WAAW,CAAA;QACnC,IAAI,MAAM,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;YACtB,MAAM,KAAK,GAAG,MAAM,CAAC,MAAM,CAAC,MAAM,GAAG,CAAC,CAAC,CAAA;YACvC,MAAM,CAAC,IAAI,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;QACnC,CAAC;IACH,CAAC;IAED,OAAO;QACL,kBAAkB;QAClB,QAAQ,EAAE;YACR,IAAI,EAAE,YAAY;YAClB,WAAW,EAAE,MAAM;SACpB;KACF,CAAA;AACH,CAAC;AAED,SAAgB,mBAAmB,CACjC,eAAmC;IAEnC,MAAM,WAAW,GAAG,eAAe,CAAC,WAAW,CAAA;IAC/C,IAAI,WAAW,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QAC7B,MAAM,iEAAiE,CAAA;IACzE,CAAC;IACD,IAAI,WAAW,CAAC,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QAC9B,MAAM,wDAAwD,CAAA;IAChE,CAAC;IACD,MAAM,gBAAgB,GAAG,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;IAC1C,MAAM,MAAM,GAAG,WAAW;SACvB,GAAG,CAAC,CAAC,KAAK,EAAE,EAAE,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;SACxB,MAAM,CACL,CAAC,WAAW,EAAE,gBAAgB,EAAE,EAAE;QAChC,MAAM,MAAM,GAAG,gBAAgB,GAAG,WAAW,CAAC,aAAa,CAAA;QAC3D,IAAI,MAAM,GAAG,CAAC,EAAE,CAAC;YACf,WAAW,CAAC,MAAM,IAAI,MAAM,CAAA;QAC9B,CAAC;aAAM,CAAC;YACN,WAAW,CAAC,OAAO,IAAI,MAAM,CAAA;QAC/B,CAAC;QACD,WAAW,CAAC,aAAa,GAAG,gBAAgB,CAAA;QAC5C,WAAW,CAAC,YAAY,GAAG,IAAI,CAAC,GAAG,CACjC,gBAAgB,EAChB,WAAW,CAAC,YAAY,CACzB,CAAA;QACD,WAAW,CAAC,YAAY,GAAG,IAAI,CAAC,GAAG,CACjC,gBAAgB,EAChB,WAAW,CAAC,YAAY,CACzB,CAAA;QAED,OAAO,WAAW,CAAA;IACpB,CAAC,EACD;QACE,MAAM,EAAE,CAAC;QACT,OAAO,EAAE,CAAC;QACV,YAAY,EAAE,gBAAgB;QAC9B,YAAY,EAAE,gBAAgB;QAC9B,aAAa,EAAE,gBAAgB;KAChC,CACF,CAAA;IAEH,OAAO;QACL,cAAc,EAAE,MAAM,CAAC,MAAM;QAC7B,eAAe,EAAE,MAAM,CAAC,OAAO;QAC/B,oBAAoB,EAAE,MAAM,CAAC,YAAY;QACzC,oBAAoB,EAAE,MAAM,CAAC,YAAY;QACzC,gBAAgB,EAAE,MAAM,CAAC,YAAY,GAAG,MAAM,CAAC,YAAY;KAC5D,CAAA;AACH,CAAC;AAED;;;;;;;;;;;;;;GAcG;AACH,SAAgB,gBAAgB,CAC9B,QAAoB,EACpB,qBAA6B;IAE7B,MAAM,kBAAkB,GAAG,2BAA2B,CACpD,QAAQ,EACR,qBAAqB,CACtB,CAAA;IAED,MAAM,UAAU,GAAG,IAAA,sBAAS,EAAC,QAAQ,EAAE,kBAAkB,EAAE;QACzD,KAAK,EAAE,QAAQ;KAChB,CAAC,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,OAAO,EAAE,EAAE,CAAC,OAAO,CAAC,QAAQ,CAAC,CAAA;IAE9C,iDAAiD;IACjD,MAAM,WAAW,GAAG,UAAU,CAAC,UAAU,CAAC,MAAM,GAAG,CAAC,CAAC,CAAA;IACrD,MAAM,iBAAiB,GAAG,IAAA,gBAAM,EAC9B,EAAE,IAAI,EAAE,SAAS,EAAE,QAAQ,EAAE,WAAW,EAAE,UAAU,EAAE,EAAE,EAAE,EAC1D,EAAE,KAAK,EAAE,QAAQ,EAAE,CACpB,CAAA;IACD,IAAI,iBAAiB,GAAG,kBAAkB,GAAG,IAAI,EAAE,CAAC;QAClD,mJAAmJ;QACnJ,UAAU,CAAC,GAAG,EAAE,CAAA;IAClB,CAAC;IAED,4GAA4G;IAC5G,uHAAuH;IACvH,MAAM,iBAAiB,GACrB,QAAQ,CAAC,WAAW,CAAC,QAAQ,CAAC,WAAW,CAAC,MAAM,GAAG,CAAC,CAAC,CAAA;IACvD,MAAM,eAAe,GAAG,UAAU,CAAC,UAAU,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,WAAW,CAAA;IACrE,eAAe,CAAC,eAAe,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,iBAAiB,CAAA;IAE/D,oFAAoF;IACpF,kFAAkF;IAClF,+CAA+C;IAC/C,6CAA6C;IAC7C,0BAA0B,CAAC,UAAU,EAAE,QAAQ,CAAC,WAAW,CAAC,CAAA;IAE5D,OAAO,EAAE,kBAAkB,EAAE,QAAQ,EAAE,UAAU,EAAE,CAAA;AACrD,CAAC;AAED,SAAS,2BAA2B,CAClC,QAAoB,EACpB,aAAqB;IAErB,MAAM,WAAW,GAAG,IAAA,gBAAM,EACxB,EAAE,IAAI,EAAE,SAAS,EAAE,QAAQ,EAAE,UAAU,EAAE,EAAE,EAAE,EAC7C,EAAE,KAAK,EAAE,QAAQ,EAAE,CACpB,CAAA;IAED,IAAI,WAAW,KAAK,CAAC,EAAE,CAAC;QACtB,OAAO,CAAC,CAAA;IACV,CAAC;IAED,MAAM,WAAW,GAAG,IAAI,CAAC,IAAI,CAAC,WAAW,GAAG,aAAa,CAAC,CAAA;IAC1D,OAAO,WAAW,GAAG,WAAW,CAAA;AAClC,CAAC;AAED;;;;GAIG;AACH,SAAS,0BAA0B,CACjC,UAAgC,EAChC,mBAAuC;IAEvC,MAAM,WAAW,GAAG,IAAI,GAAG,CAAC,mBAAmB,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAA;IAE9E,KAAK,MAAM,QAAQ,IAAI,UAAU,EAAE,CAAC;QAClC,KAAK,MAAM,KAAK,IAAI,QAAQ,CAAC,WAAW,EAAE,CAAC;YACzC,MAAM,GAAG,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,IAAI,KAAK,CAAC,CAAC,CAAC,EAAE,CAAA;YACrC,IAAI,CAAC,WAAW,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,KAAK,CAAC,MAAM,IAAI,CAAC,EAAE,CAAC;gBAC/C,KAAK,CAAC,MAAM,GAAG,CAAC,CAAA;YAClB,CAAC;QACH,CAAC;IACH,CAAC;AACH,CAAC;AAED;;;;;;GAMG;AACH,SAAS,oBAAoB,CAAC,UAAgC;IAC5D,IAAI,UAAU,CAAC,MAAM,KAAK,CAAC,EAAE,CAAC;QAC5B,OAAM;IACR,CAAC;IAED,oEAAoE;IACpE,MAAM,SAAS,GAKT,EAAE,CAAA;IAER,8DAA8D;IAC9D,IAAI,UAAU,GAAG,CAAC,CAAA;IAClB,IAAI,aAAa,GAAG,CAAC,CAAA;IACrB,IAAI,aAAa,GAA4B,IAAI,CAAA;IAEjD,KAAK,MAAM,QAAQ,IAAI,UAAU,EAAE,CAAC;QAClC,KAAK,MAAM,KAAK,IAAI,QAAQ,CAAC,WAAW,EAAE,CAAC;YACzC,yCAAyC;YACzC,IAAI,aAAa,EAAE,CAAC;gBAClB,MAAM,aAAa,GAAG,IAAA,kBAAQ,EAAC,aAAa,EAAE,KAAK,EAAE;oBACnD,KAAK,EAAE,QAAQ;iBAChB,CAAC,CAAA;gBACF,aAAa,IAAI,aAAa,CAAA;YAChC,CAAC;YAED,MAAM,YAAY,GAAG,KAAK,CAAC,MAAM,IAAI,CAAC,CAAA;YACtC,SAAS,CAAC,IAAI,CAAC;gBACb,KAAK;gBACL,YAAY;gBACZ,KAAK,EAAE,UAAU,EAAE;gBACnB,iBAAiB,EAAE,aAAa;aACjC,CAAC,CAAA;YAEF,aAAa,GAAG,KAAK,CAAA;QACvB,CAAC;IACH,CAAC;IAED,kCAAkC;IAClC,MAAM,eAAe,GAAG,SAAS,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,YAAY,CAAC,CAAA;IAE/D,IAAI,eAAe,CAAC,MAAM,IAAI,CAAC,EAAE,CAAC;QAChC,MAAM,wDAAwD,IAAI,CAAC,SAAS,CAAC,SAAS,CAAC,EAAE,CAAA;IAC3F,CAAC;IAED,yDAAyD;IACzD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,eAAe,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;QACpD,MAAM,UAAU,GAAG,eAAe,CAAC,CAAC,CAAC,CAAA;QACrC,MAAM,QAAQ,GAAG,eAAe,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;QAEvC,6BAA6B;QAC7B,IAAI,QAAQ,CAAC,KAAK,IAAI,UAAU,CAAC,KAAK,EAAE,CAAC;YACvC,MAAM,sCAAsC,CAAA;QAC9C,CAAC;QAED,MAAM,cAAc,GAAG,UAAU,CAAC,KAAK,CAAC,CAAC,CAAC,CAAA;QAC1C,MAAM,YAAY,GAAG,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAA;QACtC,MAAM,cAAc,GAAG,YAAY,GAAG,cAAc,CAAA;QAEpD,6DAA6D;QAC7D,MAAM,mBAAmB,GACvB,QAAQ,CAAC,iBAAiB,GAAG,UAAU,CAAC,iBAAiB,CAAA;QAE3D,uDAAuD;QACvD,KAAK,IAAI,CAAC,GAAG,UAAU,CAAC,KAAK,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3D,MAAM,SAAS,GAAG,SAAS,CAAC,CAAC,CAAC,CAAA;YAE9B,2DAA2D;YAC3D,IAAI,SAAS,CAAC,YAAY,EAAE,CAAC;gBAC3B,MAAM,6DAA6D,CAAA;YACrE,CAAC;YACD,qEAAqE;YACrE,MAAM,iBAAiB,GACrB,SAAS,CAAC,iBAAiB,GAAG,UAAU,CAAC,iBAAiB,CAAA;YAE5D,0FAA0F;YAC1F,IAAI,qBAAqB,GAAG,cAAc,CAAA;YAC1C,IAAI,mBAAmB,GAAG,CAAC,EAAE,CAAC;gBAC5B,MAAM,mBAAmB,GAAG,iBAAiB,GAAG,mBAAmB,CAAA;gBACnE,qBAAqB;oBACnB,cAAc,GAAG,cAAc,GAAG,mBAAmB,CAAA;YACzD,CAAC;YAED,wEAAwE;YACxE,SAAS,CAAC,KAAK,CAAC,IAAI,CAAC,qBAAqB,CAAC,CAAA;QAC7C,CAAC;IACH,CAAC;AACH,CAAC"}

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "openskidata-format",
-  "version": "9.0.0",
+  "version": "10.0.0",
   "description": "Data format for OpenSkiMap.org",
   "main": "dist/index.js",
   "types": "dist/index.d.ts",

package/src/ElevationProfile.test.ts

@@ -1,16 +1,13 @@
+import length from '@turf/length'
 import {
   extractPointsForElevationProfile,
   getAscentAndDescent,
   getPitchData,
+  lineChunkPatched,
 } from './ElevationProfile'
 import { FeatureType } from './FeatureType'
 import { getLiftElevationData, LiftFeature, LiftType } from './Lift'
-import {
-  getRunElevationData,
-  RunDifficulty,
-  RunFeature,
-  RunUse,
-} from './Run'
+import { getRunElevationData, RunDifficulty, RunFeature, RunUse } from './Run'
 import { RunDifficultyConvention } from './RunDifficultyConvention'
 import { Status } from './Status'
 
@@ -75,10 +72,10 @@ describe('ElevationProfile', () => {
       const result = getPitchData(profileGeometry)
 
       // With the new resolution calculation, the actual resolution divides the distance evenly
-      // Total length is ~60.44m, with default resolution of 25m: numSegments = ceil(60.44/25) = 3, actualResolution = 60.44/3 ≈ 20.15m
-      expect(result.pitchCalculationResolutionInMeters).toBeCloseTo(20.15, 1)
+      // Total length is ~60.44m, with default resolution of 25m: numSegments = ceil(60.44/25) = 3, actualResolution = 60.44/3 ≈ 20.14m
+      expect(result.pitchCalculationResolutionInMeters).toBeCloseTo(20.14)
       expect(result.pitchCalculationResolutionInMeters).toBeLessThanOrEqual(25)
-      expect(result.maxPitchInPercent).toBeCloseTo(0.489, 2)
+      expect(result.maxPitchInPercent).toBeCloseTo(0.489)
       expect(result.averagePitchInPercent).toBeCloseTo(0.2482)
       expect(result.inclinedLengthInMeters).toBeCloseTo(63.0597)
     })
@@ -119,8 +116,8 @@ describe('ElevationProfile', () => {
 
       // With ~222m total length and 50m input resolution:
       // numSegments = ceil(222/50) = 5
-      // actualResolution = 222/5 = 44.4m
-      expect(result.pitchCalculationResolutionInMeters).toBeCloseTo(44.4, 0)
+      // actualResolution = 222/5 ≈ 44.48m
+      expect(result.pitchCalculationResolutionInMeters).toBeCloseTo(44.48)
     })
 
     it('uses total length as resolution when line is shorter than input resolution', () => {
@@ -135,8 +132,8 @@ describe('ElevationProfile', () => {
       const result = getPitchData(profileGeometry, 50)
 
       // Resolution should equal the total length since numSegments = ceil(11/50) = 1
-      // actualResolution = 11/1 = 11m
-      expect(result.pitchCalculationResolutionInMeters).toBeCloseTo(11.1, 0)
+      // actualResolution = 11/1 ≈ 11.12m
+      expect(result.pitchCalculationResolutionInMeters).toBeCloseTo(11.12)
       expect(result.pitchCalculationResolutionInMeters).toBeLessThanOrEqual(50)
     })
 
@@ -156,7 +153,7 @@ describe('ElevationProfile', () => {
       // With evenly divided chunks, no chunks should be skipped
       // The max pitch should be calculated from all segments
       expect(result.maxPitchInPercent).toBeGreaterThan(0)
-      expect(result.pitchCalculationResolutionInMeters).toBeCloseTo(20.15, 1)
+      expect(result.pitchCalculationResolutionInMeters).toBeCloseTo(20.14)
     })
 
     it('returns null pitch data for very short lines', () => {
@@ -172,7 +169,7 @@ describe('ElevationProfile', () => {
 
       // Lines shorter than half the resolution should return null pitch data
       // because elevation data resolution makes pitch calculations unreliable
-      expect(result.pitchCalculationResolutionInMeters).toBeCloseTo(1.11, 1)
+      expect(result.pitchCalculationResolutionInMeters).toBeCloseTo(1.11)
       expect(result.maxPitchInPercent).toBeNull()
       expect(result.averagePitchInPercent).toBeNull()
       expect(result.overallPitchInPercent).toBeNull()
@@ -191,7 +188,7 @@ describe('ElevationProfile', () => {
       const result = getPitchData(profileGeometry, 25)
 
       // Lines at or above half the resolution should calculate pitch data
-      expect(result.pitchCalculationResolutionInMeters).toBeCloseTo(16.7, 0)
+      expect(result.pitchCalculationResolutionInMeters).toBeCloseTo(16.68)
       expect(result.maxPitchInPercent).not.toBeNull()
       expect(result.averagePitchInPercent).not.toBeNull()
       expect(result.overallPitchInPercent).not.toBeNull()
@@ -211,13 +208,13 @@ describe('ElevationProfile', () => {
       // With 20m resolution, there should be 12 points
       const result = extractPointsForElevationProfile(geometry, 20)
 
-      expect(result.coordinates.length).toEqual(12)
-      expect(result.coordinates[0]).toEqual([
+      expect(result.geometry.coordinates.length).toEqual(12)
+      expect(result.geometry.coordinates[0]).toEqual([
         11.177452968770694, 47.312650638218656,
       ])
-      expect(result.coordinates[result.coordinates.length - 1]).toEqual([
-        11.175409464719593, 47.31138883724759,
-      ])
+      expect(
+        result.geometry.coordinates[result.geometry.coordinates.length - 1],
+      ).toEqual([11.175409464719593, 47.31138883724759])
     })
 
     it('handles points that already have elevation', () => {
@@ -232,9 +229,9 @@ describe('ElevationProfile', () => {
       const result = extractPointsForElevationProfile(geometry, 20)
 
       // Should still work properly, returning 2D points regardless of input dimension
-      expect(result.coordinates.length).toBeGreaterThan(1)
-      expect(result.coordinates[0].length).toBe(2) // Should always return 2D points
-      expect(result.coordinates[0]).toEqual([
+      expect(result.geometry.coordinates.length).toBeGreaterThan(1)
+      expect(result.geometry.coordinates[0].length).toBe(2) // Should always return 2D points
+      expect(result.geometry.coordinates[0]).toEqual([
         11.177452968770694, 47.312650638218656,
       ])
     })
@@ -331,7 +328,12 @@ describe('ElevationProfile', () => {
         type: 'Feature',
         geometry: {
           type: 'MultiLineString',
-          coordinates: [[[0, 0, 100], [1, 1, 200]]],
+          coordinates: [
+            [
+              [0, 0, 100],
+              [1, 1, 200],
+            ],
+          ],
         },
         properties: {
           type: FeatureType.Lift,
@@ -395,7 +397,8 @@ describe('ElevationProfile', () => {
               1574, 1572, 1571, 1571, 1572, 1572, 1573, 1574, 1573, 1572, 1572,
               1573, 1575, 1576, 1577, 1580, 1583, 1584, 1586,
             ],
-            resolution: 25,
+            resolution: 24.737575778032642,
+            targetResolution: 25,
           },
           sources: [],
           websites: [],
@@ -472,7 +475,255 @@ describe('ElevationProfile', () => {
       } as RunFeature
 
       const result = getRunElevationData(runFeature)
-      expect(result?.maxPitchInPercent).toBeCloseTo(0.12, 2)
+      expect(result?.maxPitchInPercent).toBeCloseTo(0.12)
+    })
+  })
+
+  describe('lineChunkPatched', () => {
+    it('chunks a line into evenly spaced segments', () => {
+      const geometry: GeoJSON.LineString = {
+        type: 'LineString',
+        coordinates: [
+          [0, 0, 100],
+          [0, 0.001, 110], // ~111m distance
+          [0, 0.002, 105], // ~111m distance
+        ],
+      }
+
+      const result = lineChunkPatched(geometry, 50)
+
+      // Total length ~222m, with 50m input resolution:
+      // numSegments = ceil(222/50) = 5, actualResolution = 222/5 ≈ 44.48m
+      expect(result.resolutionInMeters).toBeCloseTo(44.48)
+      expect(result.geometry.length).toBe(5)
+    })
+
+    it('drops very short last segment caused by floating point precision issues', () => {
+      // This geometry used to produce a 0-length last segment
+      const geometry: GeoJSON.LineString = {
+        type: 'LineString',
+        coordinates: [
+          [11.0796655, 47.4527256, 1692.2],
+          [11.079720619903608, 47.45294351099136, 1681.8],
+          [11.079741829533312, 47.4531647479081, 1668.6],
+          [11.079727503541202, 47.453383892073894, 1656.7],
+          [11.07965005139028, 47.453599318647775, 1644.3],
+          [11.079572598604779, 47.4538147451695, 1630.9],
+          [11.07942834590789, 47.45401183160125, 1618.7],
+          [11.079261771477254, 47.45420278959679, 1607.3],
+          [11.079114389973567, 47.454399266146474, 1594.7],
+          [11.0790705, 47.45461691834462, 1583.1],
+          [11.079181898893397, 47.45481861566543, 1576.6],
+          [11.079441578857713, 47.45494597233321, 1569.5],
+          [11.079759061257176, 47.454999538815876, 1562],
+          [11.080078284221905, 47.4550501436389, 1556.1],
+          [11.080397507801036, 47.45510074757592, 1552.5],
+          [11.080713989344984, 47.45515789593482, 1548.1],
+          [11.081008027508169, 47.45525121278351, 1543.2],
+          [11.081212874246143, 47.455416412781105, 1540.9],
+          [11.081238270576348, 47.455636585144894, 1535.5],
+          [11.081208499760809, 47.45585723881769, 1526.8],
+          [11.081175991164459, 47.456077847750485, 1510.4],
+          [11.081099444790018, 47.456292146692356, 1497.9],
+          [11.080996819007606, 47.45650270883295, 1484.3],
+          [11.080897169545302, 47.45671391641284, 1472.4],
+          [11.080799550285963, 47.45692556428702, 1456.4],
+          [11.080701930240755, 47.457137212078344, 1443.9],
+          [11.08059047871202, 47.4573452579998, 1434.7],
+          [11.080452900000004, 47.45754649999998, 1427.4],
+        ],
+      }
+
+      const result = lineChunkPatched(geometry, 25)
+
+      // Verify that all segments have correct length
+      for (const chunk of result.geometry) {
+        const chunkLength = length(
+          { type: 'Feature', geometry: chunk, properties: {} },
+          { units: 'meters' },
+        )
+        expect(chunkLength).toBeCloseTo(result.resolutionInMeters)
+      }
+    })
+
+    it('preserves the exact last point from the original geometry', () => {
+      // This geometry had a split last point that was different from the input last point
+      const geometry: GeoJSON.LineString = {
+        type: 'LineString',
+        coordinates: [
+          [11.1215994, 47.466431, 886.6],
+          [11.118092699999998, 47.4702057, 768.5],
+        ],
+      }
+
+      const result = lineChunkPatched(geometry, 25)
+
+      // The last point of the last chunk should be exactly the same as the original last point
+      const lastChunk = result.geometry[result.geometry.length - 1]
+      const lastPoint = lastChunk.coordinates[lastChunk.coordinates.length - 1]
+      const originalLastPoint =
+        geometry.coordinates[geometry.coordinates.length - 1]
+
+      expect(lastPoint[0]).toBe(originalLastPoint[0])
+      expect(lastPoint[1]).toBe(originalLastPoint[1])
+      expect(lastPoint[2]).toBe(originalLastPoint[2])
+    })
+
+    it('strips elevation from non-original points', () => {
+      const geometry: GeoJSON.LineString = {
+        type: 'LineString',
+        coordinates: [
+          [0, 0, 100],
+          [0, 0.001, 110], // ~111m distance
+          [0, 0.002, 120], // ~111m distance
+        ],
+      }
+
+      const result = lineChunkPatched(geometry, 50)
+
+      // Check that interpolated points (not from original geometry) don't have elevation
+      for (const chunk of result.geometry) {
+        for (const point of chunk.coordinates) {
+          const key = `${point[0]},${point[1]}`
+          const originalKey = geometry.coordinates.some(
+            (c) => `${c[0]},${c[1]}` === key,
+          )
+
+          if (!originalKey) {
+            // Non-original points should not have elevation (should be 2D)
+            expect(point.length).toBe(2)
+          } else {
+            // Original points should still have elevation
+            expect(point.length).toBe(3)
+          }
+        }
+      }
+    })
+
+    it('handles very short lines', () => {
+      const geometry: GeoJSON.LineString = {
+        type: 'LineString',
+        coordinates: [
+          [0, 0, 100],
+          [0, 0.00001, 101], // ~1.11m distance
+        ],
+      }
+
+      const result = lineChunkPatched(geometry, 25)
+
+      // Should create a single chunk for very short lines
+      expect(result.geometry.length).toBe(1)
+      expect(result.resolutionInMeters).toBeCloseTo(1.11)
+    })
+
+    it('handles lines with exact multiples of resolution', () => {
+      const geometry: GeoJSON.LineString = {
+        type: 'LineString',
+        coordinates: [
+          [0, 0, 100],
+          [0, 0.0009, 110],
+        ],
+      }
+
+      const result = lineChunkPatched(geometry, 50)
+
+      // The actual distance is calculated and divided evenly
+      // Resolution should be <= 50m
+      expect(result.resolutionInMeters).toBeLessThanOrEqual(50)
+      expect(result.geometry.length).toBeGreaterThanOrEqual(1)
+    })
+
+    it('returns first and last points matching original geometry', () => {
+      const geometry: GeoJSON.LineString = {
+        type: 'LineString',
+        coordinates: [
+          [11.177425600412874, 47.31265682344346, 100],
+          [11.177224899122194, 47.312533118812354, 110],
+          [11.176823496540862, 47.31229807921545, 105],
+        ],
+      }
+
+      const result = lineChunkPatched(geometry, 25)
+
+      // First point of first chunk should match original first point
+      const firstPoint = result.geometry[0].coordinates[0]
+      expect(firstPoint[0]).toBe(geometry.coordinates[0][0])
+      expect(firstPoint[1]).toBe(geometry.coordinates[0][1])
+      expect(firstPoint[2]).toBe(geometry.coordinates[0][2])
+
+      // Last point of last chunk should match original last point
+      const lastChunk = result.geometry[result.geometry.length - 1]
+      const lastPoint = lastChunk.coordinates[lastChunk.coordinates.length - 1]
+      const originalLastPoint =
+        geometry.coordinates[geometry.coordinates.length - 1]
+      expect(lastPoint[0]).toBe(originalLastPoint[0])
+      expect(lastPoint[1]).toBe(originalLastPoint[1])
+      expect(lastPoint[2]).toBe(originalLastPoint[2])
+    })
+
+    it('calculates resolution that is less than or equal to requested resolution', () => {
+      const geometry: GeoJSON.LineString = {
+        type: 'LineString',
+        coordinates: [
+          [0, 0, 100],
+          [0, 0.0015, 110], // ~166.7m distance
+        ],
+      }
+
+      const result = lineChunkPatched(geometry, 50)
+
+      // Resolution should never exceed the requested minimum resolution
+      expect(result.resolutionInMeters).toBeLessThanOrEqual(50)
+    })
+
+    it('handles zero-length segments in original geometry', () => {
+      const geometry: GeoJSON.LineString = {
+        type: 'LineString',
+        coordinates: [
+          [0, 0, 100],
+          [0, 0, 100], // Duplicate point (0 length segment)
+          [0, 0.001, 110],
+        ],
+      }
+
+      const result = lineChunkPatched(geometry, 25)
+
+      // Should still work and produce valid chunks
+      expect(result.geometry.length).toBeGreaterThan(0)
+      expect(result.resolutionInMeters).toBeGreaterThan(0)
+    })
+
+    it('preserves all original points that fall on chunk boundaries', () => {
+      const geometry: GeoJSON.LineString = {
+        type: 'LineString',
+        coordinates: [
+          [0, 0, 100],
+          [0, 0.0005, 105], // ~55.6m - should be close to a chunk boundary with 50m resolution
+          [0, 0.001, 110],
+        ],
+      }
+
+      const result = lineChunkPatched(geometry, 50)
+
+      // Count how many original points are preserved in the chunks
+      let preservedOriginalPoints = 0
+      for (const chunk of result.geometry) {
+        for (const point of chunk.coordinates) {
+          if (
+            geometry.coordinates.some(
+              (c) => c[0] === point[0] && c[1] === point[1],
+            )
+          ) {
+            preservedOriginalPoints++
+          }
+        }
+      }
+
+      // All original points should be present somewhere in the chunks
+      // (Note: points may appear multiple times if they're at chunk boundaries)
+      expect(preservedOriginalPoints).toBeGreaterThanOrEqual(
+        geometry.coordinates.length,
+      )
     })
   })
 })

package/src/ElevationProfile.ts

@@ -10,10 +10,12 @@ import { LineString } from 'geojson'
  * except for the final height which corresponds to the LineString endpoint and may have a different spacing.
  * Height values can be mapped to geographical coordinates using Turf.js: `turf.lineChunk(geometry, resolution, {units: 'meters'})`.
  * @property {number} resolution - The horizontal sampling distance in meters between consecutive height measurements.
+ * @property {number} targetResolution - The horizontal sampling distance in meters that was requested for the elevation profile. The actual resolution may be less to ensure an integer number of segments along the line. This is used for reconstructing the profile geometry from the feature geometry.
  */
 export type ElevationProfile = {
   heights: number[]
   resolution: number
+  targetResolution: number
 }
 
 /**
@@ -52,7 +54,7 @@ export function getElevationData(
 
 export function getPitchData(
   profileGeometry: GeoJSON.LineString,
-  resolutionInMeters: number = 25, // Default resolution of 25 meters for pitch calculation
+  minResolutionInMeters: number = 25, // Default resolution of 25 meters for pitch calculation
 ): PitchData {
   const coordinates = profileGeometry.coordinates
   if (coordinates[0].length < 3) {
@@ -90,7 +92,7 @@ export function getPitchData(
   // If the line is too short relative to the resolution, pitch calculations are unreliable
   // due to elevation data resolution. A 1m elevation change over 1m distance would be a 45 degree slope,
   // which can easily happen due to elevation data precision issues.
-  if (totalLength < resolutionInMeters / 2) {
+  if (totalLength < minResolutionInMeters / 2) {
     return {
       averagePitchInPercent: null,
       maxPitchInPercent: null,
@@ -100,20 +102,11 @@ export function getPitchData(
     }
   }
 
-  // Calculate a resolution that divides the distance evenly, with the constraint that it must be <= resolutionInMeters
-  const numSegments = Math.ceil(totalLength / resolutionInMeters)
-  const actualResolution = totalLength / numSegments
-
   // Chunk the line into segments using the calculated resolution for max pitch calculation
-  const chunkedGeometries = lineChunk(profileGeometry, actualResolution, {
-    units: 'meters',
-  }).features.map((feature) => feature.geometry)
-
-  // turf 7.3+ adds Z coordinates to newly created points using segment-end elevation,
-  // but we need distance-based interpolation. Strip Z from non-original vertices so
-  // interpolateElevation can add them correctly.
-  // https://github.com/Turfjs/turf/issues/3007
-  stripNonOriginalElevations(chunkedGeometries, coordinates)
+  const {
+    geometry: chunkedGeometries,
+    resolutionInMeters: actualResolutionInMeters,
+  } = lineChunkPatched(profileGeometry, minResolutionInMeters)
 
   // Mutates the geometries in place to add elevation data
   interpolateElevation(chunkedGeometries)
@@ -148,7 +141,7 @@ export function getPitchData(
     maxPitchInPercent: maxPitchValue,
     inclinedLengthInMeters: totalInclinedLength,
     overallPitchInPercent: overallElevationChange / totalLength,
-    pitchCalculationResolutionInMeters: actualResolution,
+    pitchCalculationResolutionInMeters: actualResolutionInMeters,
   }
 }
 
@@ -156,14 +149,22 @@ export function getProfileGeometry(
   geometry: GeoJSON.LineString,
   elevationProfile: ElevationProfile,
 ): GeoJSON.LineString {
-  const profileLine = extractPointsForElevationProfile(
-    geometry,
-    elevationProfile.resolution,
-  )
+  const { geometry: profileLine, resolutionInMeters: actualResolution } =
+    extractPointsForElevationProfile(
+      geometry,
+      elevationProfile.targetResolution,
+    )
   if (profileLine.coordinates.length !== elevationProfile.heights.length) {
     throw `Mismatch of points & elevation profile`
   }
 
+  if (
+    Math.abs(actualResolution - elevationProfile.resolution) >
+    elevationProfile.resolution * 0.001
+  ) {
+    throw `Resolution mismatch between profile geometry (${actualResolution}) and elevation profile (${elevationProfile.resolution})`
+  }
+
   for (let i = 0; i < profileLine.coordinates.length; i++) {
     const point = profileLine.coordinates[i]
     const height = elevationProfile.heights[i]
@@ -174,16 +175,16 @@ export function getProfileGeometry(
 }
 
 /**
- * Determines the points to use for an elevation profile.
+ * Determines the points to use for an elevation profile, given a minimum horizontal resolution. The resolution is selected to have evenly spaced points along the line.
  */
 export function extractPointsForElevationProfile(
   geometry: LineString,
-  resolution: number,
-): GeoJSON.LineString {
-  // Important: Z coordinates in the output of lineChunk may be incorrect, ignore them (https://github.com/Turfjs/turf/issues/3007)
-  const lineChunks = lineChunk(geometry, resolution, {
-    units: 'meters',
-  }).features.map((feature) => feature.geometry)
+  minResolutionInMeters: number,
+): { resolutionInMeters: number; geometry: GeoJSON.LineString } {
+  const { geometry: lineChunks, resolutionInMeters } = lineChunkPatched(
+    geometry,
+    minResolutionInMeters,
+  )
 
   const points: GeoJSON.Position[] = []
   for (let subline of lineChunks) {
@@ -200,8 +201,11 @@ export function extractPointsForElevationProfile(
   }
 
   return {
-    type: 'LineString',
-    coordinates: points,
+    resolutionInMeters,
+    geometry: {
+      type: 'LineString',
+      coordinates: points,
+    },
   }
 }
 
@@ -256,6 +260,78 @@ export function getAscentAndDescent(
   }
 }
 
+/**
+ * Turf's lineChunk has several issues:
+ * - it doesn't allow for a custom resolution that divides the line evenly, which is important for consistent pitch calculations
+ *   change: calculate a resolution that divides the line into an integer number of segments, with a maximum of the provided resolution
+ * - it includes elevations from the original geometry on newly created points which are wrong (https://github.com/Turfjs/turf/issues/3007)
+ *   fix: strip Z coordinates from non-original points
+ * - it can produce a last point that is very close but not exactly the same as the original endpoint
+ *   fix: set the last point to be exactly the original endpoint
+ * - it can produce two points very close to the end of line due to floating point precision issues
+ *   fix: drop the second to last point if it's very close to the end
+ *
+ * @param geometry
+ * @param minResolutionInMeters
+ * @returns
+ */
+export function lineChunkPatched(
+  geometry: LineString,
+  minResolutionInMeters: number,
+): { resolutionInMeters: number; geometry: GeoJSON.LineString[] } {
+  const resolutionInMeters = horizontalResolutionFitting(
+    geometry,
+    minResolutionInMeters,
+  )
+
+  const lineChunks = lineChunk(geometry, resolutionInMeters, {
+    units: 'meters',
+  }).features.map((feature) => feature.geometry)
+
+  // Check the last segment has the expected length
+  const lastSegment = lineChunks[lineChunks.length - 1]
+  const lastSegmentLength = length(
+    { type: 'Feature', geometry: lastSegment, properties: {} },
+    { units: 'meters' },
+  )
+  if (lastSegmentLength < resolutionInMeters * 0.01) {
+    // If the last segment is very short, which can happen due to floating point precision issues, drop it to avoid issues with elevation interpolation
+    lineChunks.pop()
+  }
+
+  // Set the last point from the original coordinates, turf line chunk can sometimes produce a last point that
+  // is very close but not exactly the same as the original endpoint, which can cause issues with elevation interpolation
+  const lastOriginalPoint =
+    geometry.coordinates[geometry.coordinates.length - 1]
+  const lastChunkCoords = lineChunks[lineChunks.length - 1].coordinates
+  lastChunkCoords[lastChunkCoords.length - 1] = lastOriginalPoint
+
+  // turf 7.3+ adds Z coordinates to newly created points using segment-end elevation,
+  // but we need distance-based interpolation. Strip Z from non-original vertices so
+  // interpolateElevation can add them correctly.
+  // https://github.com/Turfjs/turf/issues/3007
+  stripNonOriginalElevations(lineChunks, geometry.coordinates)
+
+  return { resolutionInMeters, geometry: lineChunks }
+}
+
+function horizontalResolutionFitting(
+  geometry: LineString,
+  minResolution: number,
+): number {
+  const totalLength = length(
+    { type: 'Feature', geometry, properties: {} },
+    { units: 'meters' },
+  )
+
+  if (totalLength === 0) {
+    return 0
+  }
+
+  const numSegments = Math.ceil(totalLength / minResolution)
+  return totalLength / numSegments
+}
+
 /**
  * Strip Z coordinates from points that were created by lineChunk (not original vertices).
  * Original vertices are identified by matching lon/lat coordinates.
@@ -328,7 +404,7 @@ function interpolateElevation(geometries: GeoJSON.LineString[]) {
   const referencePoints = allPoints.filter((p) => p.hasElevation)
 
   if (referencePoints.length <= 1) {
-    throw 'At least two points with elevation data are required'
+    throw `At least two points with elevation data are required ${JSON.stringify(allPoints)}`
   }
 
   // Process segments between each pair of reference points