如何使用turf.centerOfMass计算多多边形FeatureCollection的质心？

问题场景 #

构造包含多个Polygon要素的FeatureCollection后，调用turf.centroid()和turf.center()均正常返回结果，但调用turf.centerOfMass()时抛出错误：Uncaught TypeError TypeError: Qt is not a constructor，仅传入单个Polygon要素时centerOfMass可正常运行。

可能原因 #

1. 使用的Turf.js版本存在Bug，部分旧版本（如v6.x早期版本）对FeatureCollection类型的输入处理逻辑有缺陷，导致批量处理多要素时触发异常。
2. 尽管官方文档标注centerOfMass支持FeatureCollection，但部分版本底层实现未正确兼容该类型，需先将多Polygon转换为MultiPolygon再处理。

解决方案 #

方案1：升级/调整Turf.js版本 #

部分旧版本的Turf.js存在centerOfMass处理FeatureCollection的已知Bug，建议：

• 升级至最新稳定版Turf.js（如v7.x系列）
• 若无法升级，可尝试降级到v6.5.0版本，该版本对多要素处理的兼容性已验证。

方案2：转换为MultiPolygon后计算 #

若暂时无法调整版本，可先将FeatureCollection中的多个Polygon合并为MultiPolygon要素，再传入centerOfMass：

// 将FeatureCollection的所有多边形坐标提取并组合为MultiPolygon
const multiPolygon = turf.multiPolygon(
  featCollection.features.map(feature => feature.geometry.coordinates)
);
// 计算整体质心
const centreofMass = turf.centerOfMass(multiPolygon);

方案3：手动计算加权质心 #

如果上述方案均不适用，可通过遍历单个要素质心并按面积加权的方式，手动计算整体质心：

let totalArea = 0;
let weightedX = 0;
let weightedY = 0;

featCollection.features.forEach(feature => {
  const area = turf.area(feature);
  const singleCenter = turf.centerOfMass(feature);
  const [x, y] = singleCenter.geometry.coordinates;
  
  totalArea += area;
  weightedX += x * area;
  weightedY += y * area;
});

// 生成最终质心Feature
const overallCenterOfMass = {
  type: 'Feature',
  properties: {},
  geometry: {
    type: 'Point',
    coordinates: [weightedX / totalArea, weightedY / totalArea]
  }
};

额外检查 #

确保每个Polygon要素符合GeoJSON规范：坐标数组首尾点一致（闭合）、坐标顺序为[经度, 纬度]，避免因单个要素格式问题导致批量处理异常。

内容的提问来源于stack exchange，提问作者tj26