c# MathNet.Numerics 圆拟合使用案例

下载地址

https://www.nuget.org/packages/MathNet.Numerics/4.15.0#supportedframeworks-body-tab

原理

从标准圆方程(x-c1)2+(y-c2)2=r2中进行方程变换得到2xc1+2yc2+(r2−c12−c22)=x2+y2,其中,我们c3替换常量值r2−c12−c22,即:r2−c12−c22=c3。由此,我们得到2xc1+2yc2+c3=x2+y2,将点集带入,方程就只剩三个未知数`c1,c2 和 c3。

案例

csharp 复制代码
using MathNet.Numerics;
using MathNet.Numerics.LinearAlgebra.Double;
using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;

namespace MathNetTest
{
    class Program
    {
        static void Main(string[] args)
        {
            FitCircle();
        }

        private static void FitCircle()
        {
            var count = 5000;
            var step = 2 * Math.PI / 100;
            var rd = new Random();
            //参照圆
            var x0 = 204.1;
            var y0 = 213.1;
            var r0 = 98.4;
            //噪音绝对差
            var diff = (int)(r0 * 0.1);
            var x = new double[count];
            var y = new double[count];
            //输出点集
            for (int i = 0; i < count; i++)
            {
                //circle
                x[i] = x0 + Math.Cos(i * step) * r0;
                y[i] = y0 + Math.Sin(i * step) * r0;
                //noise
                x[i] += Math.Cos(rd.Next() % 2 * Math.PI) * rd.Next(diff);
                y[i] += Math.Cos(rd.Next() % 2 * Math.PI) * rd.Next(diff);
            }

            var d = LinearAlgebraFit(x, y);


        }


        public static Tuple<double,double,double> LinearAlgebraFit(double[] X, double[] Y)
        {
            if (X.Length < 3)
            {
                return null;
            }
            var count = X.Length;
            var a = new double[count, 3];
            var c = new double[count, 1];
            for (int i = 0; i < count; i++)
            {
                //matrix
                a[i, 0] = 2 * X[i];
                a[i, 1] = 2 * Y[i];
                a[i, 2] = 1;
                c[i, 0] = X[i] * X[i] + Y[i] * Y[i];
            }
            var A = DenseMatrix.OfArray(a);
            var C = DenseMatrix.OfArray(c);
            //A*B=C
            var B = A.Solve(C);
            double c1 = B.At(0, 0),
                c2 = B.At(1, 0),
                r = Math.Sqrt(B.At(2, 0) + c1 * c1 + c2 * c2);
            Tuple<double, double, double> XYR = new Tuple<double, double, double>(c1, c2, r);
            return XYR;
        }
}
相关推荐
香蕉可乐荷包蛋1 分钟前
Python学习之路(十三)-常用函数的使用,及优化
开发语言·python·学习
惜.己9 分钟前
使用python的读取xml文件,简单的处理成元组数组
xml·开发语言·python·测试工具
apihz32 分钟前
域名WHOIS信息查询免费API使用指南
android·开发语言·数据库·网络协议·tcp/ip
coding随想1 小时前
掌控网页的魔法之书:JavaScript DOM的奇幻之旅
开发语言·javascript·ecmascript
爱吃烤鸡翅的酸菜鱼1 小时前
IDEA高效开发:Database Navigator插件安装与核心使用指南
java·开发语言·数据库·编辑器·intellij-idea·database
神仙别闹1 小时前
基于C#+SQL Server实现(Web)学生选课管理系统
前端·数据库·c#
心情好的小球藻2 小时前
Python应用进阶DAY9--类型注解Type Hinting
开发语言·python
惜.己2 小时前
使用python读取json数据,简单的处理成元组数组
开发语言·python·测试工具·json
Y4090012 小时前
C语言转Java语言,相同与相异之处
java·c语言·开发语言·笔记
向宇it2 小时前
【unity组件介绍】URP Decal Projector贴花投影器,将特定材质(贴花)投影到场景中的其他对象上。
游戏·3d·unity·c#·游戏引擎·材质