使用Turf.js缓冲OpenLayers几何图形时出现投影与EPSG不匹配导致的位移及几何拆分问题

看起来你遇到的问题确实是投影转换和坐标系不匹配导致的，尤其是在Turf.js和OpenLayers之间切换坐标系时，处理步骤出了点小问题，才会导致缓冲后的图形位移、甚至拆分成多个几何。我来帮你梳理下问题所在，再给你修正后的解决方案。

问题根源分析 #

你的代码里有个关键误区：Turf.js的buffer方法在使用kilometers这类球面单位时，必须基于WGS84（EPSG:4326）坐标系来计算。而你现在先把OpenLayers的3857几何转成4326，又用turf.toMercator转回3857再缓冲，这就犯了两个核心错误：

• 在平面投影（3857）下用球面单位（公里）缓冲，会因为投影变形（高纬度地区尤其明显）导致计算结果严重失真；
• 来回多次投影转换，容易引入坐标精度误差，甚至触发几何拓扑异常，导致图形被拆分。

修正后的解决方案 #

正确的思路是尽量减少投影转换次数，在Turf原生支持的WGS84坐标系下完成缓冲计算，再转回OpenLayers的3857投影即可：

this.#modifyEvtNullifierFeature = evt.feature;
const modifyEvtNullifierAsGeom = this.#modifyEvtNullifierFeature.getGeometry();
const geojsonFormat = new GeoJSON();

// 1. 把OpenLayers的EPSG:3857几何转换为Turf支持的EPSG:4326 GeoJSON对象
const wgs84Geometry = geojsonFormat.writeGeometryObject(modifyEvtNullifierAsGeom, {
    featureProjection: 'EPSG:3857', // OpenLayers当前使用的投影
    dataProjection: 'EPSG:4326'     // 转换为WGS84坐标系
});

// 2. 直接在WGS84坐标系下执行缓冲计算，公里单位完全适配
const bufferedWgs84 = turf.buffer(wgs84Geometry, 60, {units: 'kilometers'});

// 3. 将缓冲后的WGS84 GeoJSON转换回OpenLayers的EPSG:3857特征
const asOLFeature = geojsonFormat.readFeature(bufferedWgs84, {
    dataProjection: 'EPSG:4326',     // 源数据的坐标系
    featureProjection: 'EPSG:3857'   // 目标坐标系（OpenLayers使用的投影）
});

this.#modifyEvtNullifierFeature.set(ModifyEvtNullifierFeaturesConstants.CONST_IS_MODIFY_EVT_NULLIFIER_FEATURE.description, true);
this.#vectorSource.addFeature(asOLFeature);

额外说明 #

Web墨卡托（EPSG:3857）是平面投影，越靠近南北极坐标拉伸越严重，在这种坐标系下用球面距离单位（如公里）做空间计算，结果必然会出现位移、变形甚至几何拆分的问题。而Turf.js的空间分析方法默认是基于球面坐标系（WGS84）设计的，所以在WGS84下完成缓冲是最稳妥的选择。

另外，尽量避免多次来回转换坐标系，每一次转换都可能带来微小的精度损失，累积起来就可能引发几何异常。

备注：内容来源于stack exchange，提问作者Amr