【NetTopologySuite类库】计算多边形最大内接圆

项目版本与依赖库

项目用的.NET Framework4.7.2。

核心库：在NuGet中安装NetTopologySuite2.1及以上版本。（低版本不具有MaximumInscribedCircle这个类）

辅助库：在NuGet中安装NetTopologySuite.IO.Esri.Shapefile1.0.0版本，用于输出shapefile格式的图形。

Github上的NetTopologySuite源码地址

测试代码

   [TestClass]
    public class UnitTest1
    {
        [TestMethod]
        public void TestMethod1()
        {
            var wktReader = new WKTReader();
            var geo = wktReader.Read("polygon((0 0,0 30,40 0,0 0))");

            var maxCircleMaker = new MaximumInscribedCircle(geo, 0.1);
            var centerPoint = maxCircleMaker.GetCenter();
            var radiusPoint = maxCircleMaker.GetRadiusPoint();
            var radiusLine = maxCircleMaker.GetRadiusLine();

            var circleLine = centerPoint.Buffer(radiusLine.Length).Boundary;

            Save("./ATest/geo.shp", geo);
            Save("./ATest/centerPoint.shp", centerPoint);
            Save("./ATest/radiusPoint.shp", radiusPoint);
            Save("./ATest/radiusLine.shp", radiusLine);
            Save("./ATest/circle.shp", circleLine);
        }

        public static void Save(string shpPath, params Geometry[] geometries)
        {
            var fs = geometries.Select(g =>
            {
                return new Feature()
                {
                    Geometry = g,
                };
            });
            if (fs.Count() > 0)
            {
                var dir = Path.GetDirectoryName(shpPath);
                Directory.CreateDirectory(dir);
                Shapefile.WriteAllFeatures(fs, shpPath);
            }
        }
    }

其中new MaximumInscribedCircle(geo, 0.1);中指定了0.1的容差，发现如果容差为0，则长时间运行且不出结果。所以至少应指定一个大于0的容差。

效果图

MaximumInscribedCircle介绍

MaximumInscribedCircle出处

翻译：

为多边形构造最大内接圆，最大公差为指定公差。最大内接圆由区域内部距离区域边界最远的点以及该距离处的边界点确定。
在地理学的背景下，最大内接圆的中心被称为无法到达的极点。制图用例是确定一个合适的点来在多边形中放置地图标签。

最大内接圆的半径长度是衡量多边形“窄”程度的指标。它是负缓冲区变空的距离。

该类支持带孔的多边形和多多边形（即MultiPolygon）。

该实现在覆盖区域几何形状的正方形单元网格上使用逐次逼近技术。网格使用分支定界算法进行细化。点包含度和距离是通过使用空间索引以性能化的方式计算的。