.NET TopologySuite多边形扩展实现方法技术问询

嘿，我刚好在项目里用.NET TopologySuite（NTS）处理过类似的多边形扩展需求，这个问题核心是用**缓冲（Buffer）**方法，但要特别注意坐标系统的坑——不然你按米传参数，结果会完全不对！下面给你一步步拆解实现方式：

首先得明确：NTS的Buffer方法是基于当前几何对象的坐标单位计算的。如果你的多边形是WGS84经纬度（EPSG:4326），直接传米数会出错（因为经纬度的单位是度，1度的实际距离在不同纬度差很多）。所以第一步必须处理坐标转换。

步骤1：把经纬度坐标转成米制投影坐标系 #

如果原始多边形是WGS84，先转成米为单位的投影坐标系，比如Web墨卡托（EPSG:3857，多数场景够用）或者对应区域的UTM坐标系：

using NetTopologySuite.Geometries;
using NetTopologySuite.CoordinateSystems.Transformations;
using ProjNet.CoordinateSystems;

// 定义原始坐标系（WGS84经纬度）
var wgs84 = GeographicCoordinateSystem.WGS84;
// 定义目标米制坐标系（Web墨卡托）
var webMercator = ProjectedCoordinateSystem.WebMercator;
// 创建坐标转换器
var transform = CoordinateTransformationFactory.CreateFromCoordinateSystems(wgs84, webMercator);

// 转换原始多边形到米制坐标系
var transformedPolygon = (Polygon)transform.MathTransform.Transform(originalPolygon);

步骤2：执行多边形扩展（缓冲） #

坐标转换完成后，直接调用Buffer方法，传入你要扩展的距离（单位和投影坐标系一致，比如米）：

// 扩展1千米（1000米）
var expandedPolygon = (Polygon)transformedPolygon.Buffer(1000);
// 如果需要缩小多边形，传负数即可，比如缩小100米
// var shrunkPolygon = (Polygon)transformedPolygon.Buffer(-100);

步骤3：转换回原始坐标系统（可选） #

如果最终结果需要回到WGS84经纬度，做一次反向转换就行：

var reverseTransform = CoordinateTransformationFactory.CreateFromCoordinateSystems(webMercator, wgs84);
var finalPolygon = (Polygon)reverseTransform.MathTransform.Transform(expandedPolygon);

额外优化：控制缓冲的样式和精度 #

Buffer方法有重载可以自定义缓冲的细节，比如端帽形状、圆弧平滑度：

var bufferParams = new BufferParameters
{
    EndCapStyle = EndCapStyle.Round, // 圆形端帽（默认），也可选Flat/Square
    JoinStyle = JoinStyle.Round,
    QuadrantSegments = 16 // 分段数越高，缓冲的圆弧越平滑，默认是8
};
var expandedPolygon = (Polygon)transformedPolygon.Buffer(1000, bufferParams);

另外，缓冲后偶尔会出现无效几何（比如自相交），可以用下面的代码修复：

if (!expandedPolygon.IsValid)
{
    expandedPolygon = (Polygon)expandedPolygon.MakeValid();
}

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