基本信息

---

下载

   友情提示：（无需注册或充值，赞助后即可获取资源下载链接）

---

     嘿，亲！知识可是无价之宝呢，但咱这精心整理的资料也耗费了不少心血呀。小小地破费一下，绝对物超所值哦！如有下载和支付问题，请联系我们QQ(微信同号)：78630559

本次赞助数额为： 2 元　

支付宝支付

微信支付

   源码介绍

---

using System;
using System.Collections.Generic;
using System.ComponentModel;
using System.Data;
using System.Drawing;
using System.Linq;
using System.Text;
using System.Threading.Tasks;
using System.Windows.Forms;
using System.Drawing.Drawing2D;
using System.Drawing.Imaging;

namespace Ransac
{
    public partial class FormMain : Form
    {
        //成员
        Bitmap bmpSample;
        List<Point> points;

        public FormMain()
        {
            InitializeComponent();
        }

        /// <summary>
        /// 窗体加载事件
        /// </summary>
        /// <param name="sender"></param>
        /// <param name="e"></param>
        private void FormMain_Load(object sender, EventArgs e)
        {
            btnReset_Click(null, null);
        }

        /// <summary>
        /// 重置样本点
        /// </summary>
        /// <param name="sender"></param>
        /// <param name="e"></param>
        private void btnReset_Click(object sender, EventArgs e)
        {
            lblFormula.Text = "";
            txtRandomPoints.Text = "";
            pbSample.Image = null;
            bmpSample = new Bitmap(pbSample.Width, pbSample.Height);
        }

        /// <summary>
        /// 得到样本点
        /// </summary>
        /// <param name="sender"></param>
        /// <param name="e"></param>
        private void pbSample_Click(object sender, EventArgs e)
        {
            MouseEventArgs me=(MouseEventArgs)e;
            txtRandomPoints.Text  = string.Format("({0},{1}),", me.X, me.Y);
            DrawPoint(new Point(me.X, me.Y));
        }

        /// <summary>
        /// 在图例上画点
        /// </summary>
        /// <param name="p"></param>
        private void DrawPoint(Point p)
        {
            Graphics g = Graphics.FromImage(bmpSample);
            Pen pen = new Pen(Color.Red, 2);
            g.DrawEllipse(pen, p.X - 1, p.Y - 1, 2, 2);
            g.Flush();
            pbSample.Image = bmpSample;
        }

        /// <summary>
        /// 在图例上画直线
        /// </summary>
        /// <param name="line"></param>
        private void DrawLine(Line line)
        {
            Graphics g = Graphics.FromImage(bmpSample);
            Pen pen = new Pen(Color.Blue, 1);
            double x1 = 0;
            double y1 = line.GetY(x1);
            double x2 = pbSample.Width;
            double y2 = line.GetY(x2);
            g.DrawLine(pen, (int)x1, (int)y1, (int)x2, (int)y2);
            g.Flush();
            pbSample.Image = bmpSample;
        }

        /// <summary>
        /// 在图例上画圆
        /// </summary>
        /// <param name="circle"></param>
        private void DrawCircle(Circle circle)
        {
            Graphics g = Graphics.FromImage(bmpSample);
            Pen pen = new Pen(Color.Blue, 1);
            int x = (int)(circle.A - circle.R);
            int y = (int)(circle.B - circle.R);
            g.DrawEllipse(pen, x, y, (int)circle.R * 2, (int)circle.R * 2);
            g.Flush();
            pbSample.Image = bmpSample;
        }

        /// <summary>
        /// 得到样本点列表
        /// </summary>
        /// <returns></returns>
        private List<Point> GetSamplePoints()
        {
            List<Point> points=new List<Point>();
            string samplePoints = txtRandomPoints.Text.Replace("(", "").Replace(")", "");
            string[] arrPoints=samplePoints.Split(new char[] { ',' }, StringSplitOptions.RemoveEmptyEntries);
            int x, y;
            for (int idx = 0; idx < arrPoints.Length;idx =2 )
            {
                if (int.TryParse(arrPoints[idx], out x) && int.TryParse(arrPoints[idx   1], out y))
                    points.Add(new Point(x, y));
            }
            return points;
        }

        /// <summary>
        /// 尝试获取直线
        /// </summary>
        /// <param name="sender"></param>
        /// <param name="e"></param>
        private void btnGetLine_Click(object sender, EventArgs e)
        {
            //用RANSAC方法获取最佳直线
            points = GetSamplePoints();
            Line bestLine = null;           //最佳直线
            double bestInliersCount = 0;    //最佳模型的局内点数目
            Random random = new Random();
            for (int idx = 0; idx < nudIterCount.Value; idx  )
            {
                int idx1, idx2;
                GetRandomInliersPoints(random, out idx1, out idx2);
                int inliersCount = 2;
                Line line = new Line(points[idx1], points[idx2]);
                for (int i = 0; i < points.Count; i  )
                {
                    if (i != idx1 && i != idx2)
                    {
                        if (line.GetDistance(points[i]) <= (double)nudMinDistance.Value)
                            inliersCount  ;
                    }
                }
                if (inliersCount >= nudMinPointCount.Value)
                {
                    if (inliersCount > bestInliersCount)
                    {
                        bestLine = line;
                        bestInliersCount = inliersCount;
                    }
                }
            }
            //显示最佳直线
            if (bestLine != null)
            {
                lblFormula.Text = string.Format("方程:{0}\r\nA:{1}\r\nB:{2}\r\nC:{3}\r\n局内点数目:{4}",
                    bestLine.ToString(), bestLine.A, bestLine.B, bestLine.C, bestInliersCount);
                DrawLine(bestLine);
            }
            else
                lblFormula.Text = "没有获取到最佳直线。";
        }

        /// <summary>
        /// 随机获取两个局内点所在的索引
        /// </summary>
        /// <param name="idx1"></param>
        /// <param name="idx2"></param>
        private void GetRandomInliersPoints(Random random, out int idx1, out int idx2)
        {
            while (true)
            {
                idx1 = random.Next(points.Count);
                idx2 = random.Next(points.Count);
                if (idx1 != idx2)
                    break;
            }
        }

        /// <summary>
        /// 随机获取三个局内点所在的索引
        /// </summary>
        /// <param name="random"></param>
        /// <param name="idx1"></param>
        /// <param name="idx2"></param>
        /// <param name="idx3"></param>
        private void GetRandomInliersPoints(Random random, out int idx1, out int idx2, out int idx3)
        {
            while (true)
            {
                idx1 = random.Next(points.Count);
                idx2 = random.Next(points.Count);
                idx3 = random.Next(points.Count);
                if (idx1 != idx2 && idx2 != idx3 && idx1 != idx3)
                    break;
            }
        }

        /// <summary>
        /// 尝试获取圆
        /// </summary>
        /// <param name="sender"></param>
        /// <param name="e"></param>
        private void btnGetCircle_Click(object sender, EventArgs e)
        {
            //用RANSAC方法获取最佳直线
            points = GetSamplePoints();
            Circle bestCircle = null;       //最佳圆
            double bestInliersCount = 0;    //最佳模型的局内点数目
            Random random = new Random();
            for (int idx = 0; idx < nudIterCount.Value; idx  )
            {
                int idx1, idx2, idx3;
                GetRandomInliersPoints(random, out idx1, out idx2, out idx3);
                int inliersCount = 3;
                Circle circle;
                try
                {
                    circle = new Circle(points[idx1], points[idx2], points[idx3]);
                }
                catch
                {
                    continue;
                }
                for (int i = 0; i < points.Count; i  )
                {
                    if (i != idx1 && i != idx2 && i!=idx3)
                    {
                        if (circle.GetDistance(points[i]) <= (double)nudMinDistance.Value)
                            inliersCount  ;
                    }
                }
                if (inliersCount >= nudMinPointCount.Value)
                {
                    if (inliersCount > bestInliersCount)
                    {
                        bestCircle = circle;
                        bestInliersCount = inliersCount;
                    }
                }
            }
            //显示最佳圆
            if (bestCircle != null)
            {
                lblFormula.Text = string.Format("方程:{0}\r\nA:{1}\r\nB:{2}\r\nR:{3}\r\n局内点数目:{4}",
                    bestCircle.ToString(), bestCircle.A, bestCircle.B, bestCircle.R, bestInliersCount);
                DrawCircle(bestCircle);
            }
            else
                lblFormula.Text = "没有获取到最佳圆。";
        }
    }
}