@turf/directional-mean 7.2.0 → 7.3.0

package/dist/cjs/index.cjs

@@ -192,9 +192,9 @@ function getMeanLineString(centroidOfLine, angle, lenOfLine, isPlanar) {
     return [_invariant.getCoord.call(void 0, begin), _invariant.getCoord.call(void 0, end)];
   }
 }
-var turf_directional_mean_default = directionalMean;
+var index_default = directionalMean;
 
 
 
-exports.default = turf_directional_mean_default; exports.directionalMean = directionalMean;
+exports.default = index_default; exports.directionalMean = directionalMean;
 //# sourceMappingURL=index.cjs.map

package/dist/cjs/index.cjs.map

@@ -1 +1 @@
-{"version":3,"sources":["/home/runner/work/turf/turf/packages/turf-directional-mean/dist/cjs/index.cjs","../../index.ts"],"names":[],"mappings":"AAAA;ACCA,wCAAwB;AACxB,0CAAyB;AACzB,gDAA4B;AAC5B,wCAAqD;AACrD,4CAAyB;AACzB,sCAAuB;AACvB,kCAAwD;AA+CxD,SAAS,eAAA,CACP,KAAA,EACA,QAAA,EAGI,CAAC,CAAA,EACgB;AA5DvB,EAAA,IAAA,EAAA;AA6DE,EAAA,MAAM,SAAA,EAAW,CAAC,CAAC,OAAA,CAAQ,MAAA;AAC3B,EAAA,MAAM,UAAA,EAAA,CAAqB,GAAA,EAAA,OAAA,CAAQ,OAAA,EAAA,GAAR,KAAA,EAAA,GAAA,EAAmB,KAAA;AAC9C,EAAA,IAAI,SAAA,EAAW,CAAA;AACf,EAAA,IAAI,SAAA,EAAW,CAAA;AACf,EAAA,IAAI,aAAA,EAAe,CAAA;AACnB,EAAA,IAAI,SAAA,EAAW,CAAA;AACf,EAAA,MAAM,aAAA,EAAsC,CAAC,CAAA;AAE7C,EAAA,GAAA,CAAI,SAAA,EAAW;AACb,IAAA,+BAAA,KAAY,EAAO,CAAC,cAAA,EAAA,GAAwB;AAE1C,MAAA,MAAM,CAAC,IAAA,EAAM,IAAI,EAAA,EAAsB,YAAA;AAAA,QACrC,cAAA,CAAe,QAAA,CAAS,WAAA;AAAA,QACxB;AAAA,MACF,CAAA;AACA,MAAA,MAAM,UAAA,EAAY,qBAAA,CAAsB,cAAA,EAAgB,QAAQ,CAAA;AAChE,MAAA,GAAA,CAAI,KAAA,CAAM,IAAI,EAAA,GAAK,KAAA,CAAM,IAAI,CAAA,EAAG;AAC9B,QAAA,MAAA;AAAA,MACF,EAAA,KAAO;AACL,QAAA,SAAA,GAAY,IAAA;AACZ,QAAA,SAAA,GAAY,IAAA;AACZ,QAAA,aAAA,GAAgB,CAAA;AAChB,QAAA,SAAA,GAAY,SAAA;AACZ,QAAA,YAAA,CAAa,IAAA,CAAK,gCAAA,cAAuB,CAAC,CAAA;AAAA,MAC5C;AAAA,IACF,CAAC,CAAA;AAAA,EAEH,EAAA,KAAO;AAEL,IAAA,+BAAA,KAAY,EAAO,CAAC,cAAA,EAAA,GAAwC;AAC1D,MAAA,GAAA,CAAI,cAAA,CAAe,QAAA,CAAS,KAAA,IAAS,YAAA,EAAc;AACjD,QAAA,MAAM,IAAI,KAAA,CAAM,mCAAmC,CAAA;AAAA,MACrD;AACA,MAAA,MAAM,CAAC,IAAA,EAAM,IAAI,EAAA,EAAsB,YAAA;AAAA,QACrC,cAAA,CAAe,QAAA,CAAS,WAAA;AAAA,QACxB;AAAA,MACF,CAAA;AACA,MAAA,MAAM,UAAA,EAAY,qBAAA,CAAsB,cAAA,EAAgB,QAAQ,CAAA;AAChE,MAAA,GAAA,CAAI,KAAA,CAAM,IAAI,EAAA,GAAK,KAAA,CAAM,IAAI,CAAA,EAAG;AAC9B,QAAA,MAAA;AAAA,MACF,EAAA,KAAO;AACL,QAAA,SAAA,GAAY,IAAA;AACZ,QAAA,SAAA,GAAY,IAAA;AACZ,QAAA,aAAA,GAAgB,CAAA;AAChB,QAAA,SAAA,GAAY,SAAA;AACZ,QAAA,YAAA,CAAa,IAAA,CAAK,gCAAA,cAAuB,CAAC,CAAA;AAAA,MAC5C;AAAA,IACF,CAAC,CAAA;AAAA,EACH;AAEA,EAAA,MAAM,eAAA,EAAyB,mBAAA,CAAoB,QAAA,EAAU,QAAQ,CAAA;AACrE,EAAA,MAAM,aAAA,EAAuB,kBAAA,CAAmB,cAAc,CAAA;AAC9D,EAAA,MAAM,iBAAA,EAAmB,mBAAA;AAAA,IACvB,QAAA;AAAA,IACA,QAAA;AAAA,IACA;AAAA,EACF,CAAA;AACA,EAAA,MAAM,cAAA,EAAgB,SAAA,EAAW,YAAA;AACjC,EAAA,MAAM,gBAAA,EAAkB,gCAAA,wCAAS,YAA8B,CAAC,CAAA;AAChE,EAAA,MAAM,CAAC,QAAA,EAAU,QAAQ,EAAA,EAAc,iCAAA,eAAwB,CAAA;AAC/D,EAAA,IAAI,cAAA;AACJ,EAAA,GAAA,CAAI,QAAA,EAAU;AACZ,IAAA,eAAA,EAAiB,iBAAA;AAAA,MACf,CAAC,QAAA,EAAU,QAAQ,CAAA;AAAA,MACnB,cAAA;AAAA,MACA,aAAA;AAAA,MACA;AAAA,IACF,CAAA;AAAA,EACF,EAAA,KAAO;AACL,IAAA,eAAA,EAAiB,iBAAA;AAAA,MACf,CAAC,QAAA,EAAU,QAAQ,CAAA;AAAA,MACnB,YAAA;AAAA,MACA,aAAA;AAAA,MACA;AAAA,IACF,CAAA;AAAA,EACF;AAEA,EAAA,OAAO,iCAAA,cAAW,EAAgB;AAAA,IAChC,aAAA;AAAA,IACA,QAAA;AAAA,IACA,QAAA;AAAA,IACA,YAAA;AAAA,IACA,cAAA;AAAA,IACA,gBAAA;AAAA,IACA;AAAA,EACF,CAAC,CAAA;AACH;AAQA,SAAS,iBAAA,CAAkB,MAAA,EAAoB;AAC7C,EAAA,MAAM,CAAC,EAAA,EAAI,EAAE,EAAA,EAAc,MAAA,CAAO,CAAC,CAAA;AACnC,EAAA,MAAM,CAAC,EAAA,EAAI,EAAE,EAAA,EAAc,MAAA,CAAO,CAAC,CAAA;AACnC,EAAA,MAAM,GAAA,EAAa,GAAA,EAAK,EAAA;AACxB,EAAA,MAAM,GAAA,EAAa,GAAA,EAAK,EAAA;AACxB,EAAA,OAAO,IAAA,CAAK,IAAA,CAAK,IAAA,CAAK,GAAA,CAAI,EAAA,EAAI,CAAC,EAAA,EAAI,IAAA,CAAK,GAAA,CAAI,EAAA,EAAI,CAAC,CAAC,CAAA;AACpD;AASA,SAAS,qBAAA,CAAsB,IAAA,EAA2B,QAAA,EAAmB;AAC3E,EAAA,GAAA,CAAI,QAAA,EAAU;AACZ,IAAA,OAAO,iCAAA;AAAA,MACL,IAAA;AAAA,MACA,CAAC,aAAA,EAAwB,OAAA,EAAA,GAA0C;AACjE,QAAA,MAAM,OAAA,EAAS,OAAA,CAAS,QAAA,CAAS,WAAA;AACjC,QAAA,OAAO,cAAA,EAAiB,iBAAA,CAAkB,MAAM,CAAA;AAAA,MAClD,CAAA;AAAA,MACA;AAAA,IACF,CAAA;AAAA,EACF,EAAA,KAAO;AACL,IAAA,OAAO,4BAAA,IAAO,EAAM;AAAA,MAClB,KAAA,EAAO;AAAA,IACT,CAAC,CAAA;AAAA,EACH;AACF;AASA,SAAS,kBAAA,CAAmB,KAAA,EAAuB;AACjD,EAAA,IAAI,OAAA,EAAS,GAAA,EAAK,KAAA;AAClB,EAAA,GAAA,CAAI,OAAA,EAAS,GAAA,EAAK;AAChB,IAAA,OAAA,GAAU,GAAA;AAAA,EACZ;AACA,EAAA,OAAO,MAAA;AACT;AAQA,SAAS,YAAA,CACP,WAAA,EACA,QAAA,EACkB;AAClB,EAAA,MAAM,WAAA,EAAuB,WAAA,CAAY,CAAC,CAAA;AAC1C,EAAA,MAAM,SAAA,EAAqB,WAAA,CAAY,WAAA,CAAY,OAAA,EAAS,CAAC,CAAA;AAC7D,EAAA,GAAA,CAAI,QAAA,EAAU;AACZ,IAAA,MAAM,CAAC,EAAA,EAAI,EAAE,EAAA,EAAc,UAAA;AAC3B,IAAA,MAAM,CAAC,EAAA,EAAI,EAAE,EAAA,EAAc,QAAA;AAC3B,IAAA,MAAM,GAAA,EAAa,GAAA,EAAK,EAAA;AACxB,IAAA,MAAM,GAAA,EAAa,GAAA,EAAK,EAAA;AACxB,IAAA,MAAM,EAAA,EAAI,IAAA,CAAK,IAAA,CAAK,IAAA,CAAK,GAAA,CAAI,EAAA,EAAI,CAAC,EAAA,EAAI,IAAA,CAAK,GAAA,CAAI,EAAA,EAAI,CAAC,CAAC,CAAA;AACrD,IAAA,GAAA,CAAI,EAAA,EAAI,IAAA,EAAa;AACnB,MAAA,OAAO,CAAC,GAAA,EAAK,GAAG,CAAA;AAAA,IAClB;AACA,IAAA,MAAM,KAAA,EAAO,GAAA,EAAK,CAAA;AAClB,IAAA,MAAM,KAAA,EAAO,GAAA,EAAK,CAAA;AAClB,IAAA,OAAO,CAAC,IAAA,EAAM,IAAI,CAAA;AAAA,EACpB,EAAA,KAAO;AACL,IAAA,MAAM,MAAA,EAAQ,kBAAA,CAAmB,8BAAA,UAAQ,EAAY,QAAQ,CAAC,CAAA;AAC9D,IAAA,MAAM,OAAA,EAAU,MAAA,EAAQ,IAAA,CAAK,GAAA,EAAM,GAAA;AACnC,IAAA,OAAO,CAAC,IAAA,CAAK,GAAA,CAAI,MAAM,CAAA,EAAG,IAAA,CAAK,GAAA,CAAI,MAAM,CAAC,CAAA;AAAA,EAC5C;AACF;AAEA,SAAS,mBAAA,CAAoB,IAAA,EAAc,IAAA,EAAsB;AAC/D,EAAA,IAAI,MAAA,EAAQ,CAAA;AACZ,EAAA,GAAA,CAAI,IAAA,CAAK,GAAA,CAAI,IAAI,EAAA,EAAI,IAAA,EAAa;AAChC,IAAA,MAAA,EAAQ,EAAA;AAAA,EACV,EAAA,KAAO;AACL,IAAA,MAAA,EAAS,IAAA,CAAK,KAAA,CAAM,IAAA,EAAM,IAAI,EAAA,EAAI,IAAA,EAAO,IAAA,CAAK,EAAA;AAAA,EAChD;AACA,EAAA,GAAA,CAAI,KAAA,GAAQ,CAAA,EAAG;AACb,IAAA,GAAA,CAAI,KAAA,EAAO,CAAA,EAAG;AACZ,MAAA,MAAA,GAAS,GAAA;AAAA,IACX;AAAA,EACF,EAAA,KAAO;AACL,IAAA,GAAA,CAAI,KAAA,EAAO,CAAA,EAAG;AACZ,MAAA,MAAA,GAAS,GAAA;AAAA,IACX;AAAA,EACF;AACA,EAAA,OAAO,KAAA;AACT;AAEA,SAAS,mBAAA,CAAoB,IAAA,EAAc,IAAA,EAAc,GAAA,EAAa;AACpE,EAAA,GAAA,CAAI,IAAA,IAAQ,CAAA,EAAG;AACb,IAAA,MAAM,IAAI,KAAA,CAAM,qDAAqD,CAAA;AAAA,EACvE;AACA,EAAA,OAAO,EAAA,EAAI,IAAA,CAAK,IAAA,CAAK,IAAA,CAAK,GAAA,CAAI,IAAA,EAAM,CAAC,EAAA,EAAI,IAAA,CAAK,GAAA,CAAI,IAAA,EAAM,CAAC,CAAC,EAAA,EAAI,GAAA;AAChE;AAEA,SAAS,iBAAA,CACP,cAAA,EACA,KAAA,EACA,SAAA,EACA,QAAA,EACA;AACA,EAAA,GAAA,CAAI,QAAA,EAAU;AACZ,IAAA,MAAM,CAAC,QAAA,EAAU,QAAQ,EAAA,EAAc,cAAA;AACvC,IAAA,IAAI,MAAA;AACJ,IAAA,IAAI,MAAA;AACJ,IAAA,IAAI,IAAA;AACJ,IAAA,IAAI,IAAA;AACJ,IAAA,MAAM,EAAA,EAAa,MAAA,EAAQ,IAAA,CAAK,GAAA,EAAM,GAAA;AACtC,IAAA,MAAM,IAAA,EAAc,IAAA,CAAK,GAAA,CAAI,CAAC,CAAA;AAC9B,IAAA,MAAM,IAAA,EAAc,IAAA,CAAK,GAAA,CAAI,CAAC,CAAA;AAC9B,IAAA,OAAA,EAAS,SAAA,EAAY,UAAA,EAAY,EAAA,EAAK,GAAA;AACtC,IAAA,OAAA,EAAS,SAAA,EAAY,UAAA,EAAY,EAAA,EAAK,GAAA;AACtC,IAAA,KAAA,EAAO,SAAA,EAAY,UAAA,EAAY,EAAA,EAAK,GAAA;AACpC,IAAA,KAAA,EAAO,SAAA,EAAY,UAAA,EAAY,EAAA,EAAK,GAAA;AACpC,IAAA,OAAO;AAAA,MACL,CAAC,MAAA,EAAQ,MAAM,CAAA;AAAA,MACf,CAAC,IAAA,EAAM,IAAI;AAAA,IACb,CAAA;AAAA,EACF,EAAA,KAAO;AACL,IAAA,MAAM,IAAA,EAAM,sCAAA,4BAAY,cAAoB,CAAA,EAAG,UAAA,EAAY,CAAA,EAAG,KAAA,EAAO;AAAA,MACnE,KAAA,EAAO;AAAA,IACT,CAAC,CAAA;AACD,IAAA,MAAM,MAAA,EAAQ,sCAAA,4BAAY,cAAoB,CAAA,EAAG,CAAC,UAAA,EAAY,CAAA,EAAG,KAAA,EAAO;AAAA,MACtE,KAAA,EAAO;AAAA,IACT,CAAC,CAAA;AACD,IAAA,OAAO,CAAC,iCAAA,KAAc,CAAA,EAAG,iCAAA,GAAY,CAAC,CAAA;AAAA,EACxC;AACF;AAGA,IAAO,8BAAA,EAAQ,eAAA;ADnGf;AACE;AACA;AACF,2FAAC","file":"/home/runner/work/turf/turf/packages/turf-directional-mean/dist/cjs/index.cjs","sourcesContent":[null,"import { Feature, FeatureCollection, LineString, Point } from \"geojson\";\nimport { bearing } from \"@turf/bearing\";\nimport { centroid } from \"@turf/centroid\";\nimport { destination } from \"@turf/destination\";\nimport { featureCollection, lineString, point } from \"@turf/helpers\";\nimport { getCoord } from \"@turf/invariant\";\nimport { length } from \"@turf/length\";\nimport { featureEach, segmentEach, segmentReduce } from \"@turf/meta\";\n\ninterface DirectionalMeanLine extends Feature<LineString> {\n  properties: {\n    cartesianAngle: number;\n    bearingAngle: number;\n    circularVariance: number;\n    averageX: number;\n    averageY: number;\n    averageLength: number;\n    countOfLines: number;\n    [key: string]: any;\n  };\n}\n\n/**\n * @typedef {Object} DirectionalMeanLine\n * @property {number} cartesianAngle the mean angle of all lines. (measure from due earth counterclockwise).\n * @property {number} bearingAngle the mean angle of all lines. (bearing).\n * @property {number} circularVariance the extent to which features all point in the same direction.\n *  the value ranges 0-1, the bigger the value, the more variation in directions between lines.\n * @property {number} averageX the centroid of all lines.\n * @property {number} averageY the centroid of all line.\n * @property {number} averageLength the average length of line.\n * @property {number} countOfLines the count of features.\n */\n\n/**\n * This module calculate the average angle of a set of lines, measuring the trend of it.\n * It can be used in both project coordinate system and geography coordinate system.\n * It can handle segments of line or the whole line.\n *\n * @function\n * @param {FeatureCollection<LineString>} lines\n * @param {object} [options={}]\n * @param {boolean} [options.planar=true] whether the spatial reference system is projected or geographical.\n * @param {boolean} [options.segment=false] whether treat a LineString as a whole or a set of segments.\n * @returns {DirectionalMeanLine} Directional Mean Line\n * @example\n *\n * var lines = turf.lineStrings([\n *   [[110, 45], [120, 50]],\n *   [[100, 50], [115, 55]],\n * ])\n * var directionalMeanLine = turf.directionalMean(lines);\n * // => directionalMeanLine\n */\nfunction directionalMean(\n  lines: FeatureCollection<LineString>,\n  options: {\n    planar?: boolean;\n    segment?: boolean;\n  } = {}\n): DirectionalMeanLine {\n  const isPlanar = !!options.planar; // you can't use options.planar || true here.\n  const isSegment: boolean = options.segment ?? false;\n  let sigmaSin = 0;\n  let sigmaCos = 0;\n  let countOfLines = 0;\n  let sumOfLen = 0;\n  const centroidList: Array<Feature<Point>> = [];\n\n  if (isSegment) {\n    segmentEach(lines, (currentSegment: any) => {\n      // todo fix turf-meta's declaration file\n      const [sin1, cos1]: [number, number] = getCosAndSin(\n        currentSegment.geometry.coordinates,\n        isPlanar\n      );\n      const lenOfLine = getLengthOfLineString(currentSegment, isPlanar);\n      if (isNaN(sin1) || isNaN(cos1)) {\n        return;\n      } else {\n        sigmaSin += sin1;\n        sigmaCos += cos1;\n        countOfLines += 1;\n        sumOfLen += lenOfLine;\n        centroidList.push(centroid(currentSegment));\n      }\n    });\n    // planar and segment\n  } else {\n    // planar and non-segment\n    featureEach(lines, (currentFeature: Feature<LineString>) => {\n      if (currentFeature.geometry.type !== \"LineString\") {\n        throw new Error(\"shold to support MultiLineString?\");\n      }\n      const [sin1, cos1]: [number, number] = getCosAndSin(\n        currentFeature.geometry.coordinates,\n        isPlanar\n      );\n      const lenOfLine = getLengthOfLineString(currentFeature, isPlanar);\n      if (isNaN(sin1) || isNaN(cos1)) {\n        return;\n      } else {\n        sigmaSin += sin1;\n        sigmaCos += cos1;\n        countOfLines += 1;\n        sumOfLen += lenOfLine;\n        centroidList.push(centroid(currentFeature));\n      }\n    });\n  }\n\n  const cartesianAngle: number = getAngleBySinAndCos(sigmaSin, sigmaCos);\n  const bearingAngle: number = bearingToCartesian(cartesianAngle);\n  const circularVariance = getCircularVariance(\n    sigmaSin,\n    sigmaCos,\n    countOfLines\n  );\n  const averageLength = sumOfLen / countOfLines;\n  const centroidOfLines = centroid(featureCollection(centroidList));\n  const [averageX, averageY]: number[] = getCoord(centroidOfLines);\n  let meanLinestring;\n  if (isPlanar) {\n    meanLinestring = getMeanLineString(\n      [averageX, averageY],\n      cartesianAngle,\n      averageLength,\n      isPlanar\n    );\n  } else {\n    meanLinestring = getMeanLineString(\n      [averageX, averageY],\n      bearingAngle,\n      averageLength,\n      isPlanar\n    );\n  }\n\n  return lineString(meanLinestring, {\n    averageLength,\n    averageX,\n    averageY,\n    bearingAngle,\n    cartesianAngle,\n    circularVariance,\n    countOfLines,\n  });\n}\n\n/**\n * get euclidean distance between two points.\n * @private\n * @name euclideanDistance\n * @param coords\n */\nfunction euclideanDistance(coords: number[][]) {\n  const [x0, y0]: number[] = coords[0];\n  const [x1, y1]: number[] = coords[1];\n  const dx: number = x1 - x0;\n  const dy: number = y1 - y0;\n  return Math.sqrt(Math.pow(dx, 2) + Math.pow(dy, 2));\n}\n\n/**\n * get the length of a LineString, both in projected or geographical coordinate system.\n * @private\n * @name getLengthOfLineString\n * @param {Feature<LineString>} line\n * @param {boolean} isPlanar\n */\nfunction getLengthOfLineString(line: Feature<LineString>, isPlanar: boolean) {\n  if (isPlanar) {\n    return segmentReduce<number>(\n      line,\n      (previousValue?: number, segment?: Feature<LineString>): number => {\n        const coords = segment!.geometry.coordinates; // the signatrue of segmentReduce has problem ?\n        return previousValue! + euclideanDistance(coords);\n      },\n      0\n    );\n  } else {\n    return length(line, {\n      units: \"meters\",\n    });\n  }\n}\n\n/**\n * bearing to xy(from due earth counterclockwise 0-180)\n * convert between two forms\n * @private\n * @name bearingToCartesian\n * @param angle\n */\nfunction bearingToCartesian(angle: number): number {\n  let result = 90 - angle;\n  if (result > 180) {\n    result -= 360;\n  }\n  return result;\n}\n\n/**\n * @private\n * @name getCosAndSin\n * @param {Array<Array<number>>} coordinates\n * @returns {Array<number>} [cos, sin]\n */\nfunction getCosAndSin(\n  coordinates: number[][],\n  isPlanar: boolean\n): [number, number] {\n  const beginPoint: number[] = coordinates[0];\n  const endPoint: number[] = coordinates[coordinates.length - 1];\n  if (isPlanar) {\n    const [x0, y0]: number[] = beginPoint;\n    const [x1, y1]: number[] = endPoint;\n    const dx: number = x1 - x0;\n    const dy: number = y1 - y0;\n    const h = Math.sqrt(Math.pow(dx, 2) + Math.pow(dy, 2));\n    if (h < 0.000000001) {\n      return [NaN, NaN];\n    }\n    const sin1 = dy / h;\n    const cos1 = dx / h;\n    return [sin1, cos1];\n  } else {\n    const angle = bearingToCartesian(bearing(beginPoint, endPoint));\n    const radian = (angle * Math.PI) / 180;\n    return [Math.sin(radian), Math.cos(radian)];\n  }\n}\n\nfunction getAngleBySinAndCos(sin1: number, cos1: number): number {\n  let angle = 0;\n  if (Math.abs(cos1) < 0.000000001) {\n    angle = 90;\n  } else {\n    angle = (Math.atan2(sin1, cos1) * 180) / Math.PI;\n  }\n  if (sin1 >= 0) {\n    if (cos1 < 0) {\n      angle += 180;\n    }\n  } else {\n    if (cos1 < 0) {\n      angle -= 180;\n    }\n  }\n  return angle;\n}\n\nfunction getCircularVariance(sin1: number, cos1: number, len: number) {\n  if (len === 0) {\n    throw new Error(\"the size of the features set must be greater than 0\");\n  }\n  return 1 - Math.sqrt(Math.pow(sin1, 2) + Math.pow(cos1, 2)) / len;\n}\n\nfunction getMeanLineString(\n  centroidOfLine: number[],\n  angle: number,\n  lenOfLine: number,\n  isPlanar: boolean\n) {\n  if (isPlanar) {\n    const [averageX, averageY]: number[] = centroidOfLine;\n    let beginX: number;\n    let beginY: number;\n    let endX: number;\n    let endY: number;\n    const r: number = (angle * Math.PI) / 180;\n    const sin: number = Math.sin(r);\n    const cos: number = Math.cos(r);\n    beginX = averageX - (lenOfLine / 2) * cos;\n    beginY = averageY - (lenOfLine / 2) * sin;\n    endX = averageX + (lenOfLine / 2) * cos;\n    endY = averageY + (lenOfLine / 2) * sin;\n    return [\n      [beginX, beginY],\n      [endX, endY],\n    ];\n  } else {\n    const end = destination(point(centroidOfLine), lenOfLine / 2, angle, {\n      units: \"meters\",\n    });\n    const begin = destination(point(centroidOfLine), -lenOfLine / 2, angle, {\n      units: \"meters\",\n    });\n    return [getCoord(begin), getCoord(end)];\n  }\n}\n\nexport { directionalMean, DirectionalMeanLine };\nexport default directionalMean;\n"]}
+{"version":3,"sources":["/home/runner/work/turf/turf/packages/turf-directional-mean/dist/cjs/index.cjs","../../index.ts"],"names":[],"mappings":"AAAA;ACCA,wCAAwB;AACxB,0CAAyB;AACzB,gDAA4B;AAC5B,wCAAqD;AACrD,4CAAyB;AACzB,sCAAuB;AACvB,kCAAwD;AA+CxD,SAAS,eAAA,CACP,KAAA,EACA,QAAA,EAGI,CAAC,CAAA,EACgB;AA5DvB,EAAA,IAAA,EAAA;AA6DE,EAAA,MAAM,SAAA,EAAW,CAAC,CAAC,OAAA,CAAQ,MAAA;AAC3B,EAAA,MAAM,UAAA,EAAA,CAAqB,GAAA,EAAA,OAAA,CAAQ,OAAA,EAAA,GAAR,KAAA,EAAA,GAAA,EAAmB,KAAA;AAC9C,EAAA,IAAI,SAAA,EAAW,CAAA;AACf,EAAA,IAAI,SAAA,EAAW,CAAA;AACf,EAAA,IAAI,aAAA,EAAe,CAAA;AACnB,EAAA,IAAI,SAAA,EAAW,CAAA;AACf,EAAA,MAAM,aAAA,EAAsC,CAAC,CAAA;AAE7C,EAAA,GAAA,CAAI,SAAA,EAAW;AACb,IAAA,+BAAA,KAAY,EAAO,CAAC,cAAA,EAAA,GAAwB;AAE1C,MAAA,MAAM,CAAC,IAAA,EAAM,IAAI,EAAA,EAAsB,YAAA;AAAA,QACrC,cAAA,CAAe,QAAA,CAAS,WAAA;AAAA,QACxB;AAAA,MACF,CAAA;AACA,MAAA,MAAM,UAAA,EAAY,qBAAA,CAAsB,cAAA,EAAgB,QAAQ,CAAA;AAChE,MAAA,GAAA,CAAI,KAAA,CAAM,IAAI,EAAA,GAAK,KAAA,CAAM,IAAI,CAAA,EAAG;AAC9B,QAAA,MAAA;AAAA,MACF,EAAA,KAAO;AACL,QAAA,SAAA,GAAY,IAAA;AACZ,QAAA,SAAA,GAAY,IAAA;AACZ,QAAA,aAAA,GAAgB,CAAA;AAChB,QAAA,SAAA,GAAY,SAAA;AACZ,QAAA,YAAA,CAAa,IAAA,CAAK,gCAAA,cAAuB,CAAC,CAAA;AAAA,MAC5C;AAAA,IACF,CAAC,CAAA;AAAA,EAEH,EAAA,KAAO;AAEL,IAAA,+BAAA,KAAY,EAAO,CAAC,cAAA,EAAA,GAAwC;AAC1D,MAAA,GAAA,CAAI,cAAA,CAAe,QAAA,CAAS,KAAA,IAAS,YAAA,EAAc;AACjD,QAAA,MAAM,IAAI,KAAA,CAAM,mCAAmC,CAAA;AAAA,MACrD;AACA,MAAA,MAAM,CAAC,IAAA,EAAM,IAAI,EAAA,EAAsB,YAAA;AAAA,QACrC,cAAA,CAAe,QAAA,CAAS,WAAA;AAAA,QACxB;AAAA,MACF,CAAA;AACA,MAAA,MAAM,UAAA,EAAY,qBAAA,CAAsB,cAAA,EAAgB,QAAQ,CAAA;AAChE,MAAA,GAAA,CAAI,KAAA,CAAM,IAAI,EAAA,GAAK,KAAA,CAAM,IAAI,CAAA,EAAG;AAC9B,QAAA,MAAA;AAAA,MACF,EAAA,KAAO;AACL,QAAA,SAAA,GAAY,IAAA;AACZ,QAAA,SAAA,GAAY,IAAA;AACZ,QAAA,aAAA,GAAgB,CAAA;AAChB,QAAA,SAAA,GAAY,SAAA;AACZ,QAAA,YAAA,CAAa,IAAA,CAAK,gCAAA,cAAuB,CAAC,CAAA;AAAA,MAC5C;AAAA,IACF,CAAC,CAAA;AAAA,EACH;AAEA,EAAA,MAAM,eAAA,EAAyB,mBAAA,CAAoB,QAAA,EAAU,QAAQ,CAAA;AACrE,EAAA,MAAM,aAAA,EAAuB,kBAAA,CAAmB,cAAc,CAAA;AAC9D,EAAA,MAAM,iBAAA,EAAmB,mBAAA;AAAA,IACvB,QAAA;AAAA,IACA,QAAA;AAAA,IACA;AAAA,EACF,CAAA;AACA,EAAA,MAAM,cAAA,EAAgB,SAAA,EAAW,YAAA;AACjC,EAAA,MAAM,gBAAA,EAAkB,gCAAA,wCAAS,YAA8B,CAAC,CAAA;AAChE,EAAA,MAAM,CAAC,QAAA,EAAU,QAAQ,EAAA,EAAc,iCAAA,eAAwB,CAAA;AAC/D,EAAA,IAAI,cAAA;AACJ,EAAA,GAAA,CAAI,QAAA,EAAU;AACZ,IAAA,eAAA,EAAiB,iBAAA;AAAA,MACf,CAAC,QAAA,EAAU,QAAQ,CAAA;AAAA,MACnB,cAAA;AAAA,MACA,aAAA;AAAA,MACA;AAAA,IACF,CAAA;AAAA,EACF,EAAA,KAAO;AACL,IAAA,eAAA,EAAiB,iBAAA;AAAA,MACf,CAAC,QAAA,EAAU,QAAQ,CAAA;AAAA,MACnB,YAAA;AAAA,MACA,aAAA;AAAA,MACA;AAAA,IACF,CAAA;AAAA,EACF;AAEA,EAAA,OAAO,iCAAA,cAAW,EAAgB;AAAA,IAChC,aAAA;AAAA,IACA,QAAA;AAAA,IACA,QAAA;AAAA,IACA,YAAA;AAAA,IACA,cAAA;AAAA,IACA,gBAAA;AAAA,IACA;AAAA,EACF,CAAC,CAAA;AACH;AAQA,SAAS,iBAAA,CAAkB,MAAA,EAAoB;AAC7C,EAAA,MAAM,CAAC,EAAA,EAAI,EAAE,EAAA,EAAc,MAAA,CAAO,CAAC,CAAA;AACnC,EAAA,MAAM,CAAC,EAAA,EAAI,EAAE,EAAA,EAAc,MAAA,CAAO,CAAC,CAAA;AACnC,EAAA,MAAM,GAAA,EAAa,GAAA,EAAK,EAAA;AACxB,EAAA,MAAM,GAAA,EAAa,GAAA,EAAK,EAAA;AACxB,EAAA,OAAO,IAAA,CAAK,IAAA,CAAK,IAAA,CAAK,GAAA,CAAI,EAAA,EAAI,CAAC,EAAA,EAAI,IAAA,CAAK,GAAA,CAAI,EAAA,EAAI,CAAC,CAAC,CAAA;AACpD;AASA,SAAS,qBAAA,CAAsB,IAAA,EAA2B,QAAA,EAAmB;AAC3E,EAAA,GAAA,CAAI,QAAA,EAAU;AACZ,IAAA,OAAO,iCAAA;AAAA,MACL,IAAA;AAAA,MACA,CAAC,aAAA,EAAwB,OAAA,EAAA,GAA0C;AACjE,QAAA,MAAM,OAAA,EAAS,OAAA,CAAS,QAAA,CAAS,WAAA;AACjC,QAAA,OAAO,cAAA,EAAiB,iBAAA,CAAkB,MAAM,CAAA;AAAA,MAClD,CAAA;AAAA,MACA;AAAA,IACF,CAAA;AAAA,EACF,EAAA,KAAO;AACL,IAAA,OAAO,4BAAA,IAAO,EAAM;AAAA,MAClB,KAAA,EAAO;AAAA,IACT,CAAC,CAAA;AAAA,EACH;AACF;AASA,SAAS,kBAAA,CAAmB,KAAA,EAAuB;AACjD,EAAA,IAAI,OAAA,EAAS,GAAA,EAAK,KAAA;AAClB,EAAA,GAAA,CAAI,OAAA,EAAS,GAAA,EAAK;AAChB,IAAA,OAAA,GAAU,GAAA;AAAA,EACZ;AACA,EAAA,OAAO,MAAA;AACT;AAQA,SAAS,YAAA,CACP,WAAA,EACA,QAAA,EACkB;AAClB,EAAA,MAAM,WAAA,EAAuB,WAAA,CAAY,CAAC,CAAA;AAC1C,EAAA,MAAM,SAAA,EAAqB,WAAA,CAAY,WAAA,CAAY,OAAA,EAAS,CAAC,CAAA;AAC7D,EAAA,GAAA,CAAI,QAAA,EAAU;AACZ,IAAA,MAAM,CAAC,EAAA,EAAI,EAAE,EAAA,EAAc,UAAA;AAC3B,IAAA,MAAM,CAAC,EAAA,EAAI,EAAE,EAAA,EAAc,QAAA;AAC3B,IAAA,MAAM,GAAA,EAAa,GAAA,EAAK,EAAA;AACxB,IAAA,MAAM,GAAA,EAAa,GAAA,EAAK,EAAA;AACxB,IAAA,MAAM,EAAA,EAAI,IAAA,CAAK,IAAA,CAAK,IAAA,CAAK,GAAA,CAAI,EAAA,EAAI,CAAC,EAAA,EAAI,IAAA,CAAK,GAAA,CAAI,EAAA,EAAI,CAAC,CAAC,CAAA;AACrD,IAAA,GAAA,CAAI,EAAA,EAAI,IAAA,EAAa;AACnB,MAAA,OAAO,CAAC,GAAA,EAAK,GAAG,CAAA;AAAA,IAClB;AACA,IAAA,MAAM,KAAA,EAAO,GAAA,EAAK,CAAA;AAClB,IAAA,MAAM,KAAA,EAAO,GAAA,EAAK,CAAA;AAClB,IAAA,OAAO,CAAC,IAAA,EAAM,IAAI,CAAA;AAAA,EACpB,EAAA,KAAO;AACL,IAAA,MAAM,MAAA,EAAQ,kBAAA,CAAmB,8BAAA,UAAQ,EAAY,QAAQ,CAAC,CAAA;AAC9D,IAAA,MAAM,OAAA,EAAU,MAAA,EAAQ,IAAA,CAAK,GAAA,EAAM,GAAA;AACnC,IAAA,OAAO,CAAC,IAAA,CAAK,GAAA,CAAI,MAAM,CAAA,EAAG,IAAA,CAAK,GAAA,CAAI,MAAM,CAAC,CAAA;AAAA,EAC5C;AACF;AAEA,SAAS,mBAAA,CAAoB,IAAA,EAAc,IAAA,EAAsB;AAC/D,EAAA,IAAI,MAAA,EAAQ,CAAA;AACZ,EAAA,GAAA,CAAI,IAAA,CAAK,GAAA,CAAI,IAAI,EAAA,EAAI,IAAA,EAAa;AAChC,IAAA,MAAA,EAAQ,EAAA;AAAA,EACV,EAAA,KAAO;AACL,IAAA,MAAA,EAAS,IAAA,CAAK,KAAA,CAAM,IAAA,EAAM,IAAI,EAAA,EAAI,IAAA,EAAO,IAAA,CAAK,EAAA;AAAA,EAChD;AACA,EAAA,GAAA,CAAI,KAAA,GAAQ,CAAA,EAAG;AACb,IAAA,GAAA,CAAI,KAAA,EAAO,CAAA,EAAG;AACZ,MAAA,MAAA,GAAS,GAAA;AAAA,IACX;AAAA,EACF,EAAA,KAAO;AACL,IAAA,GAAA,CAAI,KAAA,EAAO,CAAA,EAAG;AACZ,MAAA,MAAA,GAAS,GAAA;AAAA,IACX;AAAA,EACF;AACA,EAAA,OAAO,KAAA;AACT;AAEA,SAAS,mBAAA,CAAoB,IAAA,EAAc,IAAA,EAAc,GAAA,EAAa;AACpE,EAAA,GAAA,CAAI,IAAA,IAAQ,CAAA,EAAG;AACb,IAAA,MAAM,IAAI,KAAA,CAAM,qDAAqD,CAAA;AAAA,EACvE;AACA,EAAA,OAAO,EAAA,EAAI,IAAA,CAAK,IAAA,CAAK,IAAA,CAAK,GAAA,CAAI,IAAA,EAAM,CAAC,EAAA,EAAI,IAAA,CAAK,GAAA,CAAI,IAAA,EAAM,CAAC,CAAC,EAAA,EAAI,GAAA;AAChE;AAEA,SAAS,iBAAA,CACP,cAAA,EACA,KAAA,EACA,SAAA,EACA,QAAA,EACA;AACA,EAAA,GAAA,CAAI,QAAA,EAAU;AACZ,IAAA,MAAM,CAAC,QAAA,EAAU,QAAQ,EAAA,EAAc,cAAA;AACvC,IAAA,IAAI,MAAA;AACJ,IAAA,IAAI,MAAA;AACJ,IAAA,IAAI,IAAA;AACJ,IAAA,IAAI,IAAA;AACJ,IAAA,MAAM,EAAA,EAAa,MAAA,EAAQ,IAAA,CAAK,GAAA,EAAM,GAAA;AACtC,IAAA,MAAM,IAAA,EAAc,IAAA,CAAK,GAAA,CAAI,CAAC,CAAA;AAC9B,IAAA,MAAM,IAAA,EAAc,IAAA,CAAK,GAAA,CAAI,CAAC,CAAA;AAC9B,IAAA,OAAA,EAAS,SAAA,EAAY,UAAA,EAAY,EAAA,EAAK,GAAA;AACtC,IAAA,OAAA,EAAS,SAAA,EAAY,UAAA,EAAY,EAAA,EAAK,GAAA;AACtC,IAAA,KAAA,EAAO,SAAA,EAAY,UAAA,EAAY,EAAA,EAAK,GAAA;AACpC,IAAA,KAAA,EAAO,SAAA,EAAY,UAAA,EAAY,EAAA,EAAK,GAAA;AACpC,IAAA,OAAO;AAAA,MACL,CAAC,MAAA,EAAQ,MAAM,CAAA;AAAA,MACf,CAAC,IAAA,EAAM,IAAI;AAAA,IACb,CAAA;AAAA,EACF,EAAA,KAAO;AACL,IAAA,MAAM,IAAA,EAAM,sCAAA,4BAAY,cAAoB,CAAA,EAAG,UAAA,EAAY,CAAA,EAAG,KAAA,EAAO;AAAA,MACnE,KAAA,EAAO;AAAA,IACT,CAAC,CAAA;AACD,IAAA,MAAM,MAAA,EAAQ,sCAAA,4BAAY,cAAoB,CAAA,EAAG,CAAC,UAAA,EAAY,CAAA,EAAG,KAAA,EAAO;AAAA,MACtE,KAAA,EAAO;AAAA,IACT,CAAC,CAAA;AACD,IAAA,OAAO,CAAC,iCAAA,KAAc,CAAA,EAAG,iCAAA,GAAY,CAAC,CAAA;AAAA,EACxC;AACF;AAGA,IAAO,cAAA,EAAQ,eAAA;ADnGf;AACE;AACA;AACF,2EAAC","file":"/home/runner/work/turf/turf/packages/turf-directional-mean/dist/cjs/index.cjs","sourcesContent":[null,"import { Feature, FeatureCollection, LineString, Point } from \"geojson\";\nimport { bearing } from \"@turf/bearing\";\nimport { centroid } from \"@turf/centroid\";\nimport { destination } from \"@turf/destination\";\nimport { featureCollection, lineString, point } from \"@turf/helpers\";\nimport { getCoord } from \"@turf/invariant\";\nimport { length } from \"@turf/length\";\nimport { featureEach, segmentEach, segmentReduce } from \"@turf/meta\";\n\ninterface DirectionalMeanLine extends Feature<LineString> {\n  properties: {\n    cartesianAngle: number;\n    bearingAngle: number;\n    circularVariance: number;\n    averageX: number;\n    averageY: number;\n    averageLength: number;\n    countOfLines: number;\n    [key: string]: any;\n  };\n}\n\n/**\n * @typedef {Object} DirectionalMeanLine\n * @property {number} cartesianAngle the mean angle of all lines. (measure from due earth counterclockwise).\n * @property {number} bearingAngle the mean angle of all lines. (bearing).\n * @property {number} circularVariance the extent to which features all point in the same direction.\n *  the value ranges 0-1, the bigger the value, the more variation in directions between lines.\n * @property {number} averageX the centroid of all lines.\n * @property {number} averageY the centroid of all line.\n * @property {number} averageLength the average length of line.\n * @property {number} countOfLines the count of features.\n */\n\n/**\n * This module calculate the average angle of a set of lines, measuring the trend of it.\n * It can be used in both project coordinate system and geography coordinate system.\n * It can handle segments of line or the whole line.\n *\n * @function\n * @param {FeatureCollection<LineString>} lines\n * @param {object} [options={}]\n * @param {boolean} [options.planar=true] whether the spatial reference system is projected or geographical.\n * @param {boolean} [options.segment=false] whether treat a LineString as a whole or a set of segments.\n * @returns {DirectionalMeanLine} Directional Mean Line\n * @example\n *\n * var lines = turf.lineStrings([\n *   [[110, 45], [120, 50]],\n *   [[100, 50], [115, 55]],\n * ])\n * var directionalMeanLine = turf.directionalMean(lines);\n * // => directionalMeanLine\n */\nfunction directionalMean(\n  lines: FeatureCollection<LineString>,\n  options: {\n    planar?: boolean;\n    segment?: boolean;\n  } = {}\n): DirectionalMeanLine {\n  const isPlanar = !!options.planar; // you can't use options.planar || true here.\n  const isSegment: boolean = options.segment ?? false;\n  let sigmaSin = 0;\n  let sigmaCos = 0;\n  let countOfLines = 0;\n  let sumOfLen = 0;\n  const centroidList: Array<Feature<Point>> = [];\n\n  if (isSegment) {\n    segmentEach(lines, (currentSegment: any) => {\n      // todo fix turf-meta's declaration file\n      const [sin1, cos1]: [number, number] = getCosAndSin(\n        currentSegment.geometry.coordinates,\n        isPlanar\n      );\n      const lenOfLine = getLengthOfLineString(currentSegment, isPlanar);\n      if (isNaN(sin1) || isNaN(cos1)) {\n        return;\n      } else {\n        sigmaSin += sin1;\n        sigmaCos += cos1;\n        countOfLines += 1;\n        sumOfLen += lenOfLine;\n        centroidList.push(centroid(currentSegment));\n      }\n    });\n    // planar and segment\n  } else {\n    // planar and non-segment\n    featureEach(lines, (currentFeature: Feature<LineString>) => {\n      if (currentFeature.geometry.type !== \"LineString\") {\n        throw new Error(\"shold to support MultiLineString?\");\n      }\n      const [sin1, cos1]: [number, number] = getCosAndSin(\n        currentFeature.geometry.coordinates,\n        isPlanar\n      );\n      const lenOfLine = getLengthOfLineString(currentFeature, isPlanar);\n      if (isNaN(sin1) || isNaN(cos1)) {\n        return;\n      } else {\n        sigmaSin += sin1;\n        sigmaCos += cos1;\n        countOfLines += 1;\n        sumOfLen += lenOfLine;\n        centroidList.push(centroid(currentFeature));\n      }\n    });\n  }\n\n  const cartesianAngle: number = getAngleBySinAndCos(sigmaSin, sigmaCos);\n  const bearingAngle: number = bearingToCartesian(cartesianAngle);\n  const circularVariance = getCircularVariance(\n    sigmaSin,\n    sigmaCos,\n    countOfLines\n  );\n  const averageLength = sumOfLen / countOfLines;\n  const centroidOfLines = centroid(featureCollection(centroidList));\n  const [averageX, averageY]: number[] = getCoord(centroidOfLines);\n  let meanLinestring;\n  if (isPlanar) {\n    meanLinestring = getMeanLineString(\n      [averageX, averageY],\n      cartesianAngle,\n      averageLength,\n      isPlanar\n    );\n  } else {\n    meanLinestring = getMeanLineString(\n      [averageX, averageY],\n      bearingAngle,\n      averageLength,\n      isPlanar\n    );\n  }\n\n  return lineString(meanLinestring, {\n    averageLength,\n    averageX,\n    averageY,\n    bearingAngle,\n    cartesianAngle,\n    circularVariance,\n    countOfLines,\n  });\n}\n\n/**\n * get euclidean distance between two points.\n * @private\n * @name euclideanDistance\n * @param coords\n */\nfunction euclideanDistance(coords: number[][]) {\n  const [x0, y0]: number[] = coords[0];\n  const [x1, y1]: number[] = coords[1];\n  const dx: number = x1 - x0;\n  const dy: number = y1 - y0;\n  return Math.sqrt(Math.pow(dx, 2) + Math.pow(dy, 2));\n}\n\n/**\n * get the length of a LineString, both in projected or geographical coordinate system.\n * @private\n * @name getLengthOfLineString\n * @param {Feature<LineString>} line\n * @param {boolean} isPlanar\n */\nfunction getLengthOfLineString(line: Feature<LineString>, isPlanar: boolean) {\n  if (isPlanar) {\n    return segmentReduce<number>(\n      line,\n      (previousValue?: number, segment?: Feature<LineString>): number => {\n        const coords = segment!.geometry.coordinates; // the signatrue of segmentReduce has problem ?\n        return previousValue! + euclideanDistance(coords);\n      },\n      0\n    );\n  } else {\n    return length(line, {\n      units: \"meters\",\n    });\n  }\n}\n\n/**\n * bearing to xy(from due earth counterclockwise 0-180)\n * convert between two forms\n * @private\n * @name bearingToCartesian\n * @param angle\n */\nfunction bearingToCartesian(angle: number): number {\n  let result = 90 - angle;\n  if (result > 180) {\n    result -= 360;\n  }\n  return result;\n}\n\n/**\n * @private\n * @name getCosAndSin\n * @param {Array<Array<number>>} coordinates\n * @returns {Array<number>} [cos, sin]\n */\nfunction getCosAndSin(\n  coordinates: number[][],\n  isPlanar: boolean\n): [number, number] {\n  const beginPoint: number[] = coordinates[0];\n  const endPoint: number[] = coordinates[coordinates.length - 1];\n  if (isPlanar) {\n    const [x0, y0]: number[] = beginPoint;\n    const [x1, y1]: number[] = endPoint;\n    const dx: number = x1 - x0;\n    const dy: number = y1 - y0;\n    const h = Math.sqrt(Math.pow(dx, 2) + Math.pow(dy, 2));\n    if (h < 0.000000001) {\n      return [NaN, NaN];\n    }\n    const sin1 = dy / h;\n    const cos1 = dx / h;\n    return [sin1, cos1];\n  } else {\n    const angle = bearingToCartesian(bearing(beginPoint, endPoint));\n    const radian = (angle * Math.PI) / 180;\n    return [Math.sin(radian), Math.cos(radian)];\n  }\n}\n\nfunction getAngleBySinAndCos(sin1: number, cos1: number): number {\n  let angle = 0;\n  if (Math.abs(cos1) < 0.000000001) {\n    angle = 90;\n  } else {\n    angle = (Math.atan2(sin1, cos1) * 180) / Math.PI;\n  }\n  if (sin1 >= 0) {\n    if (cos1 < 0) {\n      angle += 180;\n    }\n  } else {\n    if (cos1 < 0) {\n      angle -= 180;\n    }\n  }\n  return angle;\n}\n\nfunction getCircularVariance(sin1: number, cos1: number, len: number) {\n  if (len === 0) {\n    throw new Error(\"the size of the features set must be greater than 0\");\n  }\n  return 1 - Math.sqrt(Math.pow(sin1, 2) + Math.pow(cos1, 2)) / len;\n}\n\nfunction getMeanLineString(\n  centroidOfLine: number[],\n  angle: number,\n  lenOfLine: number,\n  isPlanar: boolean\n) {\n  if (isPlanar) {\n    const [averageX, averageY]: number[] = centroidOfLine;\n    let beginX: number;\n    let beginY: number;\n    let endX: number;\n    let endY: number;\n    const r: number = (angle * Math.PI) / 180;\n    const sin: number = Math.sin(r);\n    const cos: number = Math.cos(r);\n    beginX = averageX - (lenOfLine / 2) * cos;\n    beginY = averageY - (lenOfLine / 2) * sin;\n    endX = averageX + (lenOfLine / 2) * cos;\n    endY = averageY + (lenOfLine / 2) * sin;\n    return [\n      [beginX, beginY],\n      [endX, endY],\n    ];\n  } else {\n    const end = destination(point(centroidOfLine), lenOfLine / 2, angle, {\n      units: \"meters\",\n    });\n    const begin = destination(point(centroidOfLine), -lenOfLine / 2, angle, {\n      units: \"meters\",\n    });\n    return [getCoord(begin), getCoord(end)];\n  }\n}\n\nexport { directionalMean, DirectionalMeanLine };\nexport default directionalMean;\n"]}

package/dist/cjs/index.d.cts

@@ -1,4 +1,4 @@
-import { Feature, LineString, FeatureCollection } from 'geojson';
+import { FeatureCollection, LineString, Feature } from 'geojson';
 
 interface DirectionalMeanLine extends Feature<LineString> {
     properties: {

package/dist/esm/index.d.ts

@@ -1,4 +1,4 @@
-import { Feature, LineString, FeatureCollection } from 'geojson';
+import { FeatureCollection, LineString, Feature } from 'geojson';
 
 interface DirectionalMeanLine extends Feature<LineString> {
     properties: {

package/dist/esm/index.js

@@ -192,9 +192,9 @@ function getMeanLineString(centroidOfLine, angle, lenOfLine, isPlanar) {
     return [getCoord(begin), getCoord(end)];
   }
 }
-var turf_directional_mean_default = directionalMean;
+var index_default = directionalMean;
 export {
-  turf_directional_mean_default as default,
+  index_default as default,
   directionalMean
 };
 //# sourceMappingURL=index.js.map

package/dist/esm/index.js.map

@@ -1 +1 @@
-{"version":3,"sources":["../../index.ts"],"sourcesContent":["import { Feature, FeatureCollection, LineString, Point } from \"geojson\";\nimport { bearing } from \"@turf/bearing\";\nimport { centroid } from \"@turf/centroid\";\nimport { destination } from \"@turf/destination\";\nimport { featureCollection, lineString, point } from \"@turf/helpers\";\nimport { getCoord } from \"@turf/invariant\";\nimport { length } from \"@turf/length\";\nimport { featureEach, segmentEach, segmentReduce } from \"@turf/meta\";\n\ninterface DirectionalMeanLine extends Feature<LineString> {\n  properties: {\n    cartesianAngle: number;\n    bearingAngle: number;\n    circularVariance: number;\n    averageX: number;\n    averageY: number;\n    averageLength: number;\n    countOfLines: number;\n    [key: string]: any;\n  };\n}\n\n/**\n * @typedef {Object} DirectionalMeanLine\n * @property {number} cartesianAngle the mean angle of all lines. (measure from due earth counterclockwise).\n * @property {number} bearingAngle the mean angle of all lines. (bearing).\n * @property {number} circularVariance the extent to which features all point in the same direction.\n *  the value ranges 0-1, the bigger the value, the more variation in directions between lines.\n * @property {number} averageX the centroid of all lines.\n * @property {number} averageY the centroid of all line.\n * @property {number} averageLength the average length of line.\n * @property {number} countOfLines the count of features.\n */\n\n/**\n * This module calculate the average angle of a set of lines, measuring the trend of it.\n * It can be used in both project coordinate system and geography coordinate system.\n * It can handle segments of line or the whole line.\n *\n * @function\n * @param {FeatureCollection<LineString>} lines\n * @param {object} [options={}]\n * @param {boolean} [options.planar=true] whether the spatial reference system is projected or geographical.\n * @param {boolean} [options.segment=false] whether treat a LineString as a whole or a set of segments.\n * @returns {DirectionalMeanLine} Directional Mean Line\n * @example\n *\n * var lines = turf.lineStrings([\n *   [[110, 45], [120, 50]],\n *   [[100, 50], [115, 55]],\n * ])\n * var directionalMeanLine = turf.directionalMean(lines);\n * // => directionalMeanLine\n */\nfunction directionalMean(\n  lines: FeatureCollection<LineString>,\n  options: {\n    planar?: boolean;\n    segment?: boolean;\n  } = {}\n): DirectionalMeanLine {\n  const isPlanar = !!options.planar; // you can't use options.planar || true here.\n  const isSegment: boolean = options.segment ?? false;\n  let sigmaSin = 0;\n  let sigmaCos = 0;\n  let countOfLines = 0;\n  let sumOfLen = 0;\n  const centroidList: Array<Feature<Point>> = [];\n\n  if (isSegment) {\n    segmentEach(lines, (currentSegment: any) => {\n      // todo fix turf-meta's declaration file\n      const [sin1, cos1]: [number, number] = getCosAndSin(\n        currentSegment.geometry.coordinates,\n        isPlanar\n      );\n      const lenOfLine = getLengthOfLineString(currentSegment, isPlanar);\n      if (isNaN(sin1) || isNaN(cos1)) {\n        return;\n      } else {\n        sigmaSin += sin1;\n        sigmaCos += cos1;\n        countOfLines += 1;\n        sumOfLen += lenOfLine;\n        centroidList.push(centroid(currentSegment));\n      }\n    });\n    // planar and segment\n  } else {\n    // planar and non-segment\n    featureEach(lines, (currentFeature: Feature<LineString>) => {\n      if (currentFeature.geometry.type !== \"LineString\") {\n        throw new Error(\"shold to support MultiLineString?\");\n      }\n      const [sin1, cos1]: [number, number] = getCosAndSin(\n        currentFeature.geometry.coordinates,\n        isPlanar\n      );\n      const lenOfLine = getLengthOfLineString(currentFeature, isPlanar);\n      if (isNaN(sin1) || isNaN(cos1)) {\n        return;\n      } else {\n        sigmaSin += sin1;\n        sigmaCos += cos1;\n        countOfLines += 1;\n        sumOfLen += lenOfLine;\n        centroidList.push(centroid(currentFeature));\n      }\n    });\n  }\n\n  const cartesianAngle: number = getAngleBySinAndCos(sigmaSin, sigmaCos);\n  const bearingAngle: number = bearingToCartesian(cartesianAngle);\n  const circularVariance = getCircularVariance(\n    sigmaSin,\n    sigmaCos,\n    countOfLines\n  );\n  const averageLength = sumOfLen / countOfLines;\n  const centroidOfLines = centroid(featureCollection(centroidList));\n  const [averageX, averageY]: number[] = getCoord(centroidOfLines);\n  let meanLinestring;\n  if (isPlanar) {\n    meanLinestring = getMeanLineString(\n      [averageX, averageY],\n      cartesianAngle,\n      averageLength,\n      isPlanar\n    );\n  } else {\n    meanLinestring = getMeanLineString(\n      [averageX, averageY],\n      bearingAngle,\n      averageLength,\n      isPlanar\n    );\n  }\n\n  return lineString(meanLinestring, {\n    averageLength,\n    averageX,\n    averageY,\n    bearingAngle,\n    cartesianAngle,\n    circularVariance,\n    countOfLines,\n  });\n}\n\n/**\n * get euclidean distance between two points.\n * @private\n * @name euclideanDistance\n * @param coords\n */\nfunction euclideanDistance(coords: number[][]) {\n  const [x0, y0]: number[] = coords[0];\n  const [x1, y1]: number[] = coords[1];\n  const dx: number = x1 - x0;\n  const dy: number = y1 - y0;\n  return Math.sqrt(Math.pow(dx, 2) + Math.pow(dy, 2));\n}\n\n/**\n * get the length of a LineString, both in projected or geographical coordinate system.\n * @private\n * @name getLengthOfLineString\n * @param {Feature<LineString>} line\n * @param {boolean} isPlanar\n */\nfunction getLengthOfLineString(line: Feature<LineString>, isPlanar: boolean) {\n  if (isPlanar) {\n    return segmentReduce<number>(\n      line,\n      (previousValue?: number, segment?: Feature<LineString>): number => {\n        const coords = segment!.geometry.coordinates; // the signatrue of segmentReduce has problem ?\n        return previousValue! + euclideanDistance(coords);\n      },\n      0\n    );\n  } else {\n    return length(line, {\n      units: \"meters\",\n    });\n  }\n}\n\n/**\n * bearing to xy(from due earth counterclockwise 0-180)\n * convert between two forms\n * @private\n * @name bearingToCartesian\n * @param angle\n */\nfunction bearingToCartesian(angle: number): number {\n  let result = 90 - angle;\n  if (result > 180) {\n    result -= 360;\n  }\n  return result;\n}\n\n/**\n * @private\n * @name getCosAndSin\n * @param {Array<Array<number>>} coordinates\n * @returns {Array<number>} [cos, sin]\n */\nfunction getCosAndSin(\n  coordinates: number[][],\n  isPlanar: boolean\n): [number, number] {\n  const beginPoint: number[] = coordinates[0];\n  const endPoint: number[] = coordinates[coordinates.length - 1];\n  if (isPlanar) {\n    const [x0, y0]: number[] = beginPoint;\n    const [x1, y1]: number[] = endPoint;\n    const dx: number = x1 - x0;\n    const dy: number = y1 - y0;\n    const h = Math.sqrt(Math.pow(dx, 2) + Math.pow(dy, 2));\n    if (h < 0.000000001) {\n      return [NaN, NaN];\n    }\n    const sin1 = dy / h;\n    const cos1 = dx / h;\n    return [sin1, cos1];\n  } else {\n    const angle = bearingToCartesian(bearing(beginPoint, endPoint));\n    const radian = (angle * Math.PI) / 180;\n    return [Math.sin(radian), Math.cos(radian)];\n  }\n}\n\nfunction getAngleBySinAndCos(sin1: number, cos1: number): number {\n  let angle = 0;\n  if (Math.abs(cos1) < 0.000000001) {\n    angle = 90;\n  } else {\n    angle = (Math.atan2(sin1, cos1) * 180) / Math.PI;\n  }\n  if (sin1 >= 0) {\n    if (cos1 < 0) {\n      angle += 180;\n    }\n  } else {\n    if (cos1 < 0) {\n      angle -= 180;\n    }\n  }\n  return angle;\n}\n\nfunction getCircularVariance(sin1: number, cos1: number, len: number) {\n  if (len === 0) {\n    throw new Error(\"the size of the features set must be greater than 0\");\n  }\n  return 1 - Math.sqrt(Math.pow(sin1, 2) + Math.pow(cos1, 2)) / len;\n}\n\nfunction getMeanLineString(\n  centroidOfLine: number[],\n  angle: number,\n  lenOfLine: number,\n  isPlanar: boolean\n) {\n  if (isPlanar) {\n    const [averageX, averageY]: number[] = centroidOfLine;\n    let beginX: number;\n    let beginY: number;\n    let endX: number;\n    let endY: number;\n    const r: number = (angle * Math.PI) / 180;\n    const sin: number = Math.sin(r);\n    const cos: number = Math.cos(r);\n    beginX = averageX - (lenOfLine / 2) * cos;\n    beginY = averageY - (lenOfLine / 2) * sin;\n    endX = averageX + (lenOfLine / 2) * cos;\n    endY = averageY + (lenOfLine / 2) * sin;\n    return [\n      [beginX, beginY],\n      [endX, endY],\n    ];\n  } else {\n    const end = destination(point(centroidOfLine), lenOfLine / 2, angle, {\n      units: \"meters\",\n    });\n    const begin = destination(point(centroidOfLine), -lenOfLine / 2, angle, {\n      units: \"meters\",\n    });\n    return [getCoord(begin), getCoord(end)];\n  }\n}\n\nexport { directionalMean, DirectionalMeanLine };\nexport default directionalMean;\n"],"mappings":";AACA,SAAS,eAAe;AACxB,SAAS,gBAAgB;AACzB,SAAS,mBAAmB;AAC5B,SAAS,mBAAmB,YAAY,aAAa;AACrD,SAAS,gBAAgB;AACzB,SAAS,cAAc;AACvB,SAAS,aAAa,aAAa,qBAAqB;AA+CxD,SAAS,gBACP,OACA,UAGI,CAAC,GACgB;AA5DvB;AA6DE,QAAM,WAAW,CAAC,CAAC,QAAQ;AAC3B,QAAM,aAAqB,aAAQ,YAAR,YAAmB;AAC9C,MAAI,WAAW;AACf,MAAI,WAAW;AACf,MAAI,eAAe;AACnB,MAAI,WAAW;AACf,QAAM,eAAsC,CAAC;AAE7C,MAAI,WAAW;AACb,gBAAY,OAAO,CAAC,mBAAwB;AAE1C,YAAM,CAAC,MAAM,IAAI,IAAsB;AAAA,QACrC,eAAe,SAAS;AAAA,QACxB;AAAA,MACF;AACA,YAAM,YAAY,sBAAsB,gBAAgB,QAAQ;AAChE,UAAI,MAAM,IAAI,KAAK,MAAM,IAAI,GAAG;AAC9B;AAAA,MACF,OAAO;AACL,oBAAY;AACZ,oBAAY;AACZ,wBAAgB;AAChB,oBAAY;AACZ,qBAAa,KAAK,SAAS,cAAc,CAAC;AAAA,MAC5C;AAAA,IACF,CAAC;AAAA,EAEH,OAAO;AAEL,gBAAY,OAAO,CAAC,mBAAwC;AAC1D,UAAI,eAAe,SAAS,SAAS,cAAc;AACjD,cAAM,IAAI,MAAM,mCAAmC;AAAA,MACrD;AACA,YAAM,CAAC,MAAM,IAAI,IAAsB;AAAA,QACrC,eAAe,SAAS;AAAA,QACxB;AAAA,MACF;AACA,YAAM,YAAY,sBAAsB,gBAAgB,QAAQ;AAChE,UAAI,MAAM,IAAI,KAAK,MAAM,IAAI,GAAG;AAC9B;AAAA,MACF,OAAO;AACL,oBAAY;AACZ,oBAAY;AACZ,wBAAgB;AAChB,oBAAY;AACZ,qBAAa,KAAK,SAAS,cAAc,CAAC;AAAA,MAC5C;AAAA,IACF,CAAC;AAAA,EACH;AAEA,QAAM,iBAAyB,oBAAoB,UAAU,QAAQ;AACrE,QAAM,eAAuB,mBAAmB,cAAc;AAC9D,QAAM,mBAAmB;AAAA,IACvB;AAAA,IACA;AAAA,IACA;AAAA,EACF;AACA,QAAM,gBAAgB,WAAW;AACjC,QAAM,kBAAkB,SAAS,kBAAkB,YAAY,CAAC;AAChE,QAAM,CAAC,UAAU,QAAQ,IAAc,SAAS,eAAe;AAC/D,MAAI;AACJ,MAAI,UAAU;AACZ,qBAAiB;AAAA,MACf,CAAC,UAAU,QAAQ;AAAA,MACnB;AAAA,MACA;AAAA,MACA;AAAA,IACF;AAAA,EACF,OAAO;AACL,qBAAiB;AAAA,MACf,CAAC,UAAU,QAAQ;AAAA,MACnB;AAAA,MACA;AAAA,MACA;AAAA,IACF;AAAA,EACF;AAEA,SAAO,WAAW,gBAAgB;AAAA,IAChC;AAAA,IACA;AAAA,IACA;AAAA,IACA;AAAA,IACA;AAAA,IACA;AAAA,IACA;AAAA,EACF,CAAC;AACH;AAQA,SAAS,kBAAkB,QAAoB;AAC7C,QAAM,CAAC,IAAI,EAAE,IAAc,OAAO,CAAC;AACnC,QAAM,CAAC,IAAI,EAAE,IAAc,OAAO,CAAC;AACnC,QAAM,KAAa,KAAK;AACxB,QAAM,KAAa,KAAK;AACxB,SAAO,KAAK,KAAK,KAAK,IAAI,IAAI,CAAC,IAAI,KAAK,IAAI,IAAI,CAAC,CAAC;AACpD;AASA,SAAS,sBAAsB,MAA2B,UAAmB;AAC3E,MAAI,UAAU;AACZ,WAAO;AAAA,MACL;AAAA,MACA,CAAC,eAAwB,YAA0C;AACjE,cAAM,SAAS,QAAS,SAAS;AACjC,eAAO,gBAAiB,kBAAkB,MAAM;AAAA,MAClD;AAAA,MACA;AAAA,IACF;AAAA,EACF,OAAO;AACL,WAAO,OAAO,MAAM;AAAA,MAClB,OAAO;AAAA,IACT,CAAC;AAAA,EACH;AACF;AASA,SAAS,mBAAmB,OAAuB;AACjD,MAAI,SAAS,KAAK;AAClB,MAAI,SAAS,KAAK;AAChB,cAAU;AAAA,EACZ;AACA,SAAO;AACT;AAQA,SAAS,aACP,aACA,UACkB;AAClB,QAAM,aAAuB,YAAY,CAAC;AAC1C,QAAM,WAAqB,YAAY,YAAY,SAAS,CAAC;AAC7D,MAAI,UAAU;AACZ,UAAM,CAAC,IAAI,EAAE,IAAc;AAC3B,UAAM,CAAC,IAAI,EAAE,IAAc;AAC3B,UAAM,KAAa,KAAK;AACxB,UAAM,KAAa,KAAK;AACxB,UAAM,IAAI,KAAK,KAAK,KAAK,IAAI,IAAI,CAAC,IAAI,KAAK,IAAI,IAAI,CAAC,CAAC;AACrD,QAAI,IAAI,MAAa;AACnB,aAAO,CAAC,KAAK,GAAG;AAAA,IAClB;AACA,UAAM,OAAO,KAAK;AAClB,UAAM,OAAO,KAAK;AAClB,WAAO,CAAC,MAAM,IAAI;AAAA,EACpB,OAAO;AACL,UAAM,QAAQ,mBAAmB,QAAQ,YAAY,QAAQ,CAAC;AAC9D,UAAM,SAAU,QAAQ,KAAK,KAAM;AACnC,WAAO,CAAC,KAAK,IAAI,MAAM,GAAG,KAAK,IAAI,MAAM,CAAC;AAAA,EAC5C;AACF;AAEA,SAAS,oBAAoB,MAAc,MAAsB;AAC/D,MAAI,QAAQ;AACZ,MAAI,KAAK,IAAI,IAAI,IAAI,MAAa;AAChC,YAAQ;AAAA,EACV,OAAO;AACL,YAAS,KAAK,MAAM,MAAM,IAAI,IAAI,MAAO,KAAK;AAAA,EAChD;AACA,MAAI,QAAQ,GAAG;AACb,QAAI,OAAO,GAAG;AACZ,eAAS;AAAA,IACX;AAAA,EACF,OAAO;AACL,QAAI,OAAO,GAAG;AACZ,eAAS;AAAA,IACX;AAAA,EACF;AACA,SAAO;AACT;AAEA,SAAS,oBAAoB,MAAc,MAAc,KAAa;AACpE,MAAI,QAAQ,GAAG;AACb,UAAM,IAAI,MAAM,qDAAqD;AAAA,EACvE;AACA,SAAO,IAAI,KAAK,KAAK,KAAK,IAAI,MAAM,CAAC,IAAI,KAAK,IAAI,MAAM,CAAC,CAAC,IAAI;AAChE;AAEA,SAAS,kBACP,gBACA,OACA,WACA,UACA;AACA,MAAI,UAAU;AACZ,UAAM,CAAC,UAAU,QAAQ,IAAc;AACvC,QAAI;AACJ,QAAI;AACJ,QAAI;AACJ,QAAI;AACJ,UAAM,IAAa,QAAQ,KAAK,KAAM;AACtC,UAAM,MAAc,KAAK,IAAI,CAAC;AAC9B,UAAM,MAAc,KAAK,IAAI,CAAC;AAC9B,aAAS,WAAY,YAAY,IAAK;AACtC,aAAS,WAAY,YAAY,IAAK;AACtC,WAAO,WAAY,YAAY,IAAK;AACpC,WAAO,WAAY,YAAY,IAAK;AACpC,WAAO;AAAA,MACL,CAAC,QAAQ,MAAM;AAAA,MACf,CAAC,MAAM,IAAI;AAAA,IACb;AAAA,EACF,OAAO;AACL,UAAM,MAAM,YAAY,MAAM,cAAc,GAAG,YAAY,GAAG,OAAO;AAAA,MACnE,OAAO;AAAA,IACT,CAAC;AACD,UAAM,QAAQ,YAAY,MAAM,cAAc,GAAG,CAAC,YAAY,GAAG,OAAO;AAAA,MACtE,OAAO;AAAA,IACT,CAAC;AACD,WAAO,CAAC,SAAS,KAAK,GAAG,SAAS,GAAG,CAAC;AAAA,EACxC;AACF;AAGA,IAAO,gCAAQ;","names":[]}
+{"version":3,"sources":["../../index.ts"],"sourcesContent":["import { Feature, FeatureCollection, LineString, Point } from \"geojson\";\nimport { bearing } from \"@turf/bearing\";\nimport { centroid } from \"@turf/centroid\";\nimport { destination } from \"@turf/destination\";\nimport { featureCollection, lineString, point } from \"@turf/helpers\";\nimport { getCoord } from \"@turf/invariant\";\nimport { length } from \"@turf/length\";\nimport { featureEach, segmentEach, segmentReduce } from \"@turf/meta\";\n\ninterface DirectionalMeanLine extends Feature<LineString> {\n  properties: {\n    cartesianAngle: number;\n    bearingAngle: number;\n    circularVariance: number;\n    averageX: number;\n    averageY: number;\n    averageLength: number;\n    countOfLines: number;\n    [key: string]: any;\n  };\n}\n\n/**\n * @typedef {Object} DirectionalMeanLine\n * @property {number} cartesianAngle the mean angle of all lines. (measure from due earth counterclockwise).\n * @property {number} bearingAngle the mean angle of all lines. (bearing).\n * @property {number} circularVariance the extent to which features all point in the same direction.\n *  the value ranges 0-1, the bigger the value, the more variation in directions between lines.\n * @property {number} averageX the centroid of all lines.\n * @property {number} averageY the centroid of all line.\n * @property {number} averageLength the average length of line.\n * @property {number} countOfLines the count of features.\n */\n\n/**\n * This module calculate the average angle of a set of lines, measuring the trend of it.\n * It can be used in both project coordinate system and geography coordinate system.\n * It can handle segments of line or the whole line.\n *\n * @function\n * @param {FeatureCollection<LineString>} lines\n * @param {object} [options={}]\n * @param {boolean} [options.planar=true] whether the spatial reference system is projected or geographical.\n * @param {boolean} [options.segment=false] whether treat a LineString as a whole or a set of segments.\n * @returns {DirectionalMeanLine} Directional Mean Line\n * @example\n *\n * var lines = turf.lineStrings([\n *   [[110, 45], [120, 50]],\n *   [[100, 50], [115, 55]],\n * ])\n * var directionalMeanLine = turf.directionalMean(lines);\n * // => directionalMeanLine\n */\nfunction directionalMean(\n  lines: FeatureCollection<LineString>,\n  options: {\n    planar?: boolean;\n    segment?: boolean;\n  } = {}\n): DirectionalMeanLine {\n  const isPlanar = !!options.planar; // you can't use options.planar || true here.\n  const isSegment: boolean = options.segment ?? false;\n  let sigmaSin = 0;\n  let sigmaCos = 0;\n  let countOfLines = 0;\n  let sumOfLen = 0;\n  const centroidList: Array<Feature<Point>> = [];\n\n  if (isSegment) {\n    segmentEach(lines, (currentSegment: any) => {\n      // todo fix turf-meta's declaration file\n      const [sin1, cos1]: [number, number] = getCosAndSin(\n        currentSegment.geometry.coordinates,\n        isPlanar\n      );\n      const lenOfLine = getLengthOfLineString(currentSegment, isPlanar);\n      if (isNaN(sin1) || isNaN(cos1)) {\n        return;\n      } else {\n        sigmaSin += sin1;\n        sigmaCos += cos1;\n        countOfLines += 1;\n        sumOfLen += lenOfLine;\n        centroidList.push(centroid(currentSegment));\n      }\n    });\n    // planar and segment\n  } else {\n    // planar and non-segment\n    featureEach(lines, (currentFeature: Feature<LineString>) => {\n      if (currentFeature.geometry.type !== \"LineString\") {\n        throw new Error(\"shold to support MultiLineString?\");\n      }\n      const [sin1, cos1]: [number, number] = getCosAndSin(\n        currentFeature.geometry.coordinates,\n        isPlanar\n      );\n      const lenOfLine = getLengthOfLineString(currentFeature, isPlanar);\n      if (isNaN(sin1) || isNaN(cos1)) {\n        return;\n      } else {\n        sigmaSin += sin1;\n        sigmaCos += cos1;\n        countOfLines += 1;\n        sumOfLen += lenOfLine;\n        centroidList.push(centroid(currentFeature));\n      }\n    });\n  }\n\n  const cartesianAngle: number = getAngleBySinAndCos(sigmaSin, sigmaCos);\n  const bearingAngle: number = bearingToCartesian(cartesianAngle);\n  const circularVariance = getCircularVariance(\n    sigmaSin,\n    sigmaCos,\n    countOfLines\n  );\n  const averageLength = sumOfLen / countOfLines;\n  const centroidOfLines = centroid(featureCollection(centroidList));\n  const [averageX, averageY]: number[] = getCoord(centroidOfLines);\n  let meanLinestring;\n  if (isPlanar) {\n    meanLinestring = getMeanLineString(\n      [averageX, averageY],\n      cartesianAngle,\n      averageLength,\n      isPlanar\n    );\n  } else {\n    meanLinestring = getMeanLineString(\n      [averageX, averageY],\n      bearingAngle,\n      averageLength,\n      isPlanar\n    );\n  }\n\n  return lineString(meanLinestring, {\n    averageLength,\n    averageX,\n    averageY,\n    bearingAngle,\n    cartesianAngle,\n    circularVariance,\n    countOfLines,\n  });\n}\n\n/**\n * get euclidean distance between two points.\n * @private\n * @name euclideanDistance\n * @param coords\n */\nfunction euclideanDistance(coords: number[][]) {\n  const [x0, y0]: number[] = coords[0];\n  const [x1, y1]: number[] = coords[1];\n  const dx: number = x1 - x0;\n  const dy: number = y1 - y0;\n  return Math.sqrt(Math.pow(dx, 2) + Math.pow(dy, 2));\n}\n\n/**\n * get the length of a LineString, both in projected or geographical coordinate system.\n * @private\n * @name getLengthOfLineString\n * @param {Feature<LineString>} line\n * @param {boolean} isPlanar\n */\nfunction getLengthOfLineString(line: Feature<LineString>, isPlanar: boolean) {\n  if (isPlanar) {\n    return segmentReduce<number>(\n      line,\n      (previousValue?: number, segment?: Feature<LineString>): number => {\n        const coords = segment!.geometry.coordinates; // the signatrue of segmentReduce has problem ?\n        return previousValue! + euclideanDistance(coords);\n      },\n      0\n    );\n  } else {\n    return length(line, {\n      units: \"meters\",\n    });\n  }\n}\n\n/**\n * bearing to xy(from due earth counterclockwise 0-180)\n * convert between two forms\n * @private\n * @name bearingToCartesian\n * @param angle\n */\nfunction bearingToCartesian(angle: number): number {\n  let result = 90 - angle;\n  if (result > 180) {\n    result -= 360;\n  }\n  return result;\n}\n\n/**\n * @private\n * @name getCosAndSin\n * @param {Array<Array<number>>} coordinates\n * @returns {Array<number>} [cos, sin]\n */\nfunction getCosAndSin(\n  coordinates: number[][],\n  isPlanar: boolean\n): [number, number] {\n  const beginPoint: number[] = coordinates[0];\n  const endPoint: number[] = coordinates[coordinates.length - 1];\n  if (isPlanar) {\n    const [x0, y0]: number[] = beginPoint;\n    const [x1, y1]: number[] = endPoint;\n    const dx: number = x1 - x0;\n    const dy: number = y1 - y0;\n    const h = Math.sqrt(Math.pow(dx, 2) + Math.pow(dy, 2));\n    if (h < 0.000000001) {\n      return [NaN, NaN];\n    }\n    const sin1 = dy / h;\n    const cos1 = dx / h;\n    return [sin1, cos1];\n  } else {\n    const angle = bearingToCartesian(bearing(beginPoint, endPoint));\n    const radian = (angle * Math.PI) / 180;\n    return [Math.sin(radian), Math.cos(radian)];\n  }\n}\n\nfunction getAngleBySinAndCos(sin1: number, cos1: number): number {\n  let angle = 0;\n  if (Math.abs(cos1) < 0.000000001) {\n    angle = 90;\n  } else {\n    angle = (Math.atan2(sin1, cos1) * 180) / Math.PI;\n  }\n  if (sin1 >= 0) {\n    if (cos1 < 0) {\n      angle += 180;\n    }\n  } else {\n    if (cos1 < 0) {\n      angle -= 180;\n    }\n  }\n  return angle;\n}\n\nfunction getCircularVariance(sin1: number, cos1: number, len: number) {\n  if (len === 0) {\n    throw new Error(\"the size of the features set must be greater than 0\");\n  }\n  return 1 - Math.sqrt(Math.pow(sin1, 2) + Math.pow(cos1, 2)) / len;\n}\n\nfunction getMeanLineString(\n  centroidOfLine: number[],\n  angle: number,\n  lenOfLine: number,\n  isPlanar: boolean\n) {\n  if (isPlanar) {\n    const [averageX, averageY]: number[] = centroidOfLine;\n    let beginX: number;\n    let beginY: number;\n    let endX: number;\n    let endY: number;\n    const r: number = (angle * Math.PI) / 180;\n    const sin: number = Math.sin(r);\n    const cos: number = Math.cos(r);\n    beginX = averageX - (lenOfLine / 2) * cos;\n    beginY = averageY - (lenOfLine / 2) * sin;\n    endX = averageX + (lenOfLine / 2) * cos;\n    endY = averageY + (lenOfLine / 2) * sin;\n    return [\n      [beginX, beginY],\n      [endX, endY],\n    ];\n  } else {\n    const end = destination(point(centroidOfLine), lenOfLine / 2, angle, {\n      units: \"meters\",\n    });\n    const begin = destination(point(centroidOfLine), -lenOfLine / 2, angle, {\n      units: \"meters\",\n    });\n    return [getCoord(begin), getCoord(end)];\n  }\n}\n\nexport { directionalMean, DirectionalMeanLine };\nexport default directionalMean;\n"],"mappings":";AACA,SAAS,eAAe;AACxB,SAAS,gBAAgB;AACzB,SAAS,mBAAmB;AAC5B,SAAS,mBAAmB,YAAY,aAAa;AACrD,SAAS,gBAAgB;AACzB,SAAS,cAAc;AACvB,SAAS,aAAa,aAAa,qBAAqB;AA+CxD,SAAS,gBACP,OACA,UAGI,CAAC,GACgB;AA5DvB;AA6DE,QAAM,WAAW,CAAC,CAAC,QAAQ;AAC3B,QAAM,aAAqB,aAAQ,YAAR,YAAmB;AAC9C,MAAI,WAAW;AACf,MAAI,WAAW;AACf,MAAI,eAAe;AACnB,MAAI,WAAW;AACf,QAAM,eAAsC,CAAC;AAE7C,MAAI,WAAW;AACb,gBAAY,OAAO,CAAC,mBAAwB;AAE1C,YAAM,CAAC,MAAM,IAAI,IAAsB;AAAA,QACrC,eAAe,SAAS;AAAA,QACxB;AAAA,MACF;AACA,YAAM,YAAY,sBAAsB,gBAAgB,QAAQ;AAChE,UAAI,MAAM,IAAI,KAAK,MAAM,IAAI,GAAG;AAC9B;AAAA,MACF,OAAO;AACL,oBAAY;AACZ,oBAAY;AACZ,wBAAgB;AAChB,oBAAY;AACZ,qBAAa,KAAK,SAAS,cAAc,CAAC;AAAA,MAC5C;AAAA,IACF,CAAC;AAAA,EAEH,OAAO;AAEL,gBAAY,OAAO,CAAC,mBAAwC;AAC1D,UAAI,eAAe,SAAS,SAAS,cAAc;AACjD,cAAM,IAAI,MAAM,mCAAmC;AAAA,MACrD;AACA,YAAM,CAAC,MAAM,IAAI,IAAsB;AAAA,QACrC,eAAe,SAAS;AAAA,QACxB;AAAA,MACF;AACA,YAAM,YAAY,sBAAsB,gBAAgB,QAAQ;AAChE,UAAI,MAAM,IAAI,KAAK,MAAM,IAAI,GAAG;AAC9B;AAAA,MACF,OAAO;AACL,oBAAY;AACZ,oBAAY;AACZ,wBAAgB;AAChB,oBAAY;AACZ,qBAAa,KAAK,SAAS,cAAc,CAAC;AAAA,MAC5C;AAAA,IACF,CAAC;AAAA,EACH;AAEA,QAAM,iBAAyB,oBAAoB,UAAU,QAAQ;AACrE,QAAM,eAAuB,mBAAmB,cAAc;AAC9D,QAAM,mBAAmB;AAAA,IACvB;AAAA,IACA;AAAA,IACA;AAAA,EACF;AACA,QAAM,gBAAgB,WAAW;AACjC,QAAM,kBAAkB,SAAS,kBAAkB,YAAY,CAAC;AAChE,QAAM,CAAC,UAAU,QAAQ,IAAc,SAAS,eAAe;AAC/D,MAAI;AACJ,MAAI,UAAU;AACZ,qBAAiB;AAAA,MACf,CAAC,UAAU,QAAQ;AAAA,MACnB;AAAA,MACA;AAAA,MACA;AAAA,IACF;AAAA,EACF,OAAO;AACL,qBAAiB;AAAA,MACf,CAAC,UAAU,QAAQ;AAAA,MACnB;AAAA,MACA;AAAA,MACA;AAAA,IACF;AAAA,EACF;AAEA,SAAO,WAAW,gBAAgB;AAAA,IAChC;AAAA,IACA;AAAA,IACA;AAAA,IACA;AAAA,IACA;AAAA,IACA;AAAA,IACA;AAAA,EACF,CAAC;AACH;AAQA,SAAS,kBAAkB,QAAoB;AAC7C,QAAM,CAAC,IAAI,EAAE,IAAc,OAAO,CAAC;AACnC,QAAM,CAAC,IAAI,EAAE,IAAc,OAAO,CAAC;AACnC,QAAM,KAAa,KAAK;AACxB,QAAM,KAAa,KAAK;AACxB,SAAO,KAAK,KAAK,KAAK,IAAI,IAAI,CAAC,IAAI,KAAK,IAAI,IAAI,CAAC,CAAC;AACpD;AASA,SAAS,sBAAsB,MAA2B,UAAmB;AAC3E,MAAI,UAAU;AACZ,WAAO;AAAA,MACL;AAAA,MACA,CAAC,eAAwB,YAA0C;AACjE,cAAM,SAAS,QAAS,SAAS;AACjC,eAAO,gBAAiB,kBAAkB,MAAM;AAAA,MAClD;AAAA,MACA;AAAA,IACF;AAAA,EACF,OAAO;AACL,WAAO,OAAO,MAAM;AAAA,MAClB,OAAO;AAAA,IACT,CAAC;AAAA,EACH;AACF;AASA,SAAS,mBAAmB,OAAuB;AACjD,MAAI,SAAS,KAAK;AAClB,MAAI,SAAS,KAAK;AAChB,cAAU;AAAA,EACZ;AACA,SAAO;AACT;AAQA,SAAS,aACP,aACA,UACkB;AAClB,QAAM,aAAuB,YAAY,CAAC;AAC1C,QAAM,WAAqB,YAAY,YAAY,SAAS,CAAC;AAC7D,MAAI,UAAU;AACZ,UAAM,CAAC,IAAI,EAAE,IAAc;AAC3B,UAAM,CAAC,IAAI,EAAE,IAAc;AAC3B,UAAM,KAAa,KAAK;AACxB,UAAM,KAAa,KAAK;AACxB,UAAM,IAAI,KAAK,KAAK,KAAK,IAAI,IAAI,CAAC,IAAI,KAAK,IAAI,IAAI,CAAC,CAAC;AACrD,QAAI,IAAI,MAAa;AACnB,aAAO,CAAC,KAAK,GAAG;AAAA,IAClB;AACA,UAAM,OAAO,KAAK;AAClB,UAAM,OAAO,KAAK;AAClB,WAAO,CAAC,MAAM,IAAI;AAAA,EACpB,OAAO;AACL,UAAM,QAAQ,mBAAmB,QAAQ,YAAY,QAAQ,CAAC;AAC9D,UAAM,SAAU,QAAQ,KAAK,KAAM;AACnC,WAAO,CAAC,KAAK,IAAI,MAAM,GAAG,KAAK,IAAI,MAAM,CAAC;AAAA,EAC5C;AACF;AAEA,SAAS,oBAAoB,MAAc,MAAsB;AAC/D,MAAI,QAAQ;AACZ,MAAI,KAAK,IAAI,IAAI,IAAI,MAAa;AAChC,YAAQ;AAAA,EACV,OAAO;AACL,YAAS,KAAK,MAAM,MAAM,IAAI,IAAI,MAAO,KAAK;AAAA,EAChD;AACA,MAAI,QAAQ,GAAG;AACb,QAAI,OAAO,GAAG;AACZ,eAAS;AAAA,IACX;AAAA,EACF,OAAO;AACL,QAAI,OAAO,GAAG;AACZ,eAAS;AAAA,IACX;AAAA,EACF;AACA,SAAO;AACT;AAEA,SAAS,oBAAoB,MAAc,MAAc,KAAa;AACpE,MAAI,QAAQ,GAAG;AACb,UAAM,IAAI,MAAM,qDAAqD;AAAA,EACvE;AACA,SAAO,IAAI,KAAK,KAAK,KAAK,IAAI,MAAM,CAAC,IAAI,KAAK,IAAI,MAAM,CAAC,CAAC,IAAI;AAChE;AAEA,SAAS,kBACP,gBACA,OACA,WACA,UACA;AACA,MAAI,UAAU;AACZ,UAAM,CAAC,UAAU,QAAQ,IAAc;AACvC,QAAI;AACJ,QAAI;AACJ,QAAI;AACJ,QAAI;AACJ,UAAM,IAAa,QAAQ,KAAK,KAAM;AACtC,UAAM,MAAc,KAAK,IAAI,CAAC;AAC9B,UAAM,MAAc,KAAK,IAAI,CAAC;AAC9B,aAAS,WAAY,YAAY,IAAK;AACtC,aAAS,WAAY,YAAY,IAAK;AACtC,WAAO,WAAY,YAAY,IAAK;AACpC,WAAO,WAAY,YAAY,IAAK;AACpC,WAAO;AAAA,MACL,CAAC,QAAQ,MAAM;AAAA,MACf,CAAC,MAAM,IAAI;AAAA,IACb;AAAA,EACF,OAAO;AACL,UAAM,MAAM,YAAY,MAAM,cAAc,GAAG,YAAY,GAAG,OAAO;AAAA,MACnE,OAAO;AAAA,IACT,CAAC;AACD,UAAM,QAAQ,YAAY,MAAM,cAAc,GAAG,CAAC,YAAY,GAAG,OAAO;AAAA,MACtE,OAAO;AAAA,IACT,CAAC;AACD,WAAO,CAAC,SAAS,KAAK,GAAG,SAAS,GAAG,CAAC;AAAA,EACxC;AACF;AAGA,IAAO,gBAAQ;","names":[]}

package/package.json

@@ -1,7 +1,7 @@
 {
   "name": "@turf/directional-mean",
-  "version": "7.2.0",
-  "description": "turf directional-mean module",
+  "version": "7.3.0",
+  "description": "Calculates the average angle of a set of lines, measuring the trend of it.",
   "author": "Turf Authors",
   "contributors": [
     "Haoming Zhuang <@zhuang-hao-ming>"
@@ -53,26 +53,26 @@
   },
   "devDependencies": {
     "@types/benchmark": "^2.1.5",
-    "@types/tape": "^4.13.4",
+    "@types/tape": "^5.8.1",
     "benchmark": "^2.1.4",
     "load-json-file": "^7.0.1",
     "npm-run-all": "^4.1.5",
     "tape": "^5.9.0",
-    "tsup": "^8.3.5",
-    "tsx": "^4.19.2",
-    "typescript": "^5.5.4",
-    "write-json-file": "^5.0.0"
+    "tsup": "^8.4.0",
+    "tsx": "^4.19.4",
+    "typescript": "^5.8.3",
+    "write-json-file": "^6.0.0"
   },
   "dependencies": {
-    "@turf/bearing": "^7.2.0",
-    "@turf/centroid": "^7.2.0",
-    "@turf/destination": "^7.2.0",
-    "@turf/helpers": "^7.2.0",
-    "@turf/invariant": "^7.2.0",
-    "@turf/length": "^7.2.0",
-    "@turf/meta": "^7.2.0",
+    "@turf/bearing": "7.3.0",
+    "@turf/centroid": "7.3.0",
+    "@turf/destination": "7.3.0",
+    "@turf/helpers": "7.3.0",
+    "@turf/invariant": "7.3.0",
+    "@turf/length": "7.3.0",
+    "@turf/meta": "7.3.0",
     "@types/geojson": "^7946.0.10",
     "tslib": "^2.8.1"
   },
-  "gitHead": "7b0f0374c4668cd569f8904c71e2ae7d941be867"
+  "gitHead": "9f58a103e8f9a587ab640307ed03ba5233913ddd"
 }