码迷,mamicode.com
首页 > 其他好文 > 详细

几何算法:点集合构造简单多边形

时间:2014-05-29 04:07:56      阅读:365      评论:0      收藏:0      [点我收藏+]

标签:style   c   class   blog   code   tar   

问题:给定平面中n个点所组成的集合,将它们连接起来形成一条简单的封闭路径。所谓简单路径,是指边与边无交叉。

如下图所示10个点组成的简单轮廓:

bubuko.com,布布扣

思路:取x坐标最大的点A(如果最大x坐标的点不止一个,则取Y坐标最小的点),依次计算A点与其余各点的连线与水平线之间夹角的正切值,然后按照正切值排序,依次连接排序后的各点即组成一个简单图形。

原理:其它所有点都在A点的左侧,所有夹角的范围为-Pi/2~Pi/2,单调递增函数。

bubuko.com,布布扣

举一个例子如下:

bubuko.com,布布扣

各点坐标与A点的角度斜率如下(已经排序好):

x:426.192518536091,y:30.5668629242884,slope:-2.21036105157629
x:132.904271903869,y:111.805767306036,slope:0.0233827696146631
x:209.153583263584,y:158.396180071121,slope:0.216615047225945
x:51.2625493860163,y:271.425922467106,slope:0.409713066051227
x:172.80558813494,y:320.363658168522,slope:0.754116336162768
x:174.841647802313,y:361.474091434606,slope:0.903935084923323
x:262.993097888768,y:306.679940091763,slope:1.03059799172764
x:405.520514378101,y:212.478244240618,slope:2.00680658499766
x:410.405247491042,y:324.597360433357,slope:4.49064367657446
x:459.491329337233,y:104.169257382941,slope:1.79769313486232E+308

其中A点为:x:459.491329337233,y:104.169257382941,slope:1.79769313486232E+308

下面给出具体算法(C#实现):

几何点定义,实现IComparable<T>接口,按照正切值排序要用到:

bubuko.com,布布扣
  public struct GeometryPoint : IComparable<GeometryPoint>
    {        
        public GeometryPoint(double x, double y, double slope = double.NaN)
        {
            this.x = x;
            this.y = y;
            this.slope = slope;
        }
        private double x;
        public double X
        {
            get { return x; }
            set { x = value; }
        }
        private double y;
        public double Y
        {
            get { return y; }
            set { y = value; }
        }
        private double slope;
        public double SLOPE
        {
            get { return slope; }
            set { slope = value; }
        }
      
        public int CompareTo(GeometryPoint p)
        {
            if (this.slope < p.slope)
            {
                return -1;
            }
            else if (this.slope > p.slope)
            {
                return 1;
            }
            else
            {
                if (this.x == p.x && this.SLOPE == p.SLOPE && this.SLOPE == double.MaxValue)
                {
                    if (this.y == p.y)
                    {
                        return 0;
                    }
                    else if (this.y < p.y)
                    {
                        return 1;
                    }
                    else//(this.y > p.y)
                    {
                        return -1;
                    }
                }
                return 0;
            }
        }
        public override string ToString()
        {
            return string.Format("x:{0},y:{1},slope:{2}", x, y, slope);
        }
    }
GeometryPoint 定义

简单封闭图形定义,并定义初始化简单封闭图形的方法,该方法随机产生多边形的顶点:

bubuko.com,布布扣
 public class SimplePolygon
    {
        private GeometryPoint[] geometrypoints;

        public GeometryPoint[] GeometryPoints
        {
            get { return geometrypoints; }
            set { geometrypoints = value; }
        }


        public SimplePolygon()
        {
        }
        public void Initialize(int size, double minX, double maxX, double minY, double maxY)
        {
            if (size <= 0) throw new ArgumentOutOfRangeException();
            geometrypoints = new GeometryPoint[size];
            Random rnd = new Random(DateTime.Now.Millisecond);
            double xRange = maxX - minX;
            double yRange = maxY - minY;
            int MaxXPointIndex = 0;//选取x坐标最大的点
            for (int i = 0; i < size; i++)
            {
                GeometryPoint gp = new GeometryPoint(minX + xRange * rnd.NextDouble(), minY + yRange * rnd.NextDouble());
                geometrypoints[i] = gp;
                if (geometrypoints[MaxXPointIndex].X < gp.X)////选取x坐标最大的点
                {
                    MaxXPointIndex = i;
                }
                else if (geometrypoints[MaxXPointIndex].X < gp.X && geometrypoints[MaxXPointIndex].Y > gp.Y)//选取x坐标最大的点,如果最大x坐标点有多个,去y最小者
                {
                    MaxXPointIndex = i;
                }
            }
            //计算斜率
            for (int i = 0; i < size; i++)
            {
                if (i == MaxXPointIndex)
                {
                    geometrypoints[MaxXPointIndex].SLOPE = double.MaxValue;
                }
                else
                {
                    if (geometrypoints[i].X == geometrypoints[MaxXPointIndex].X)//与最大x坐标的x相同的点,因为x坐标之差为零,所以取SLOPE最大值
                    {
                        geometrypoints[i].SLOPE = double.MaxValue;
                    }
                    else//计算斜率,注意正切函数在-0.5Pi和0.5Pi之间是单调递增的
                    {
                        geometrypoints[i].SLOPE = (geometrypoints[i].Y - geometrypoints[MaxXPointIndex].Y) / (geometrypoints[MaxXPointIndex].X - geometrypoints[i].X);
                    }
                }
            }
            //按照斜率slope排序,取稳定排序方法的堆排序。
            HeapSort<GeometryPoint> heapsort = new HeapSort<GeometryPoint>();
            heapsort.Sort(this.geometrypoints,0,size-1);
        }
    }
SimplePolygon定义

控制台程序调用方法,按照连线顺序打印顶点:

bubuko.com,布布扣
 class Program
    {
        static void Main(string[] args)
        {
            SimplePolygon sp = new SimplePolygon();
            sp.Initialize(10, -50, 50, -50, 50);
            for (int i = 0; i < sp.GeometryPoints.Length; i++)
            {
                Console.WriteLine(sp.GeometryPoints[i]);
            }
            Console.ReadKey();
        }
    }
SimplePolygon 控制台调用方法

如果用界面绘图,应用WPF几何绘图可实现如下效果,红线为计算正切值的示例连线,绿色线为生成的简单多边形:

关于坐标系与绘图的方法,请参照另一篇文章“轮廓算法”。

bubuko.com,布布扣

完毕。

作者:Andy Zeng

欢迎任何形式的转载,但请务必注明出处。

http://www.cnblogs.com/andyzeng/p/3754005.html

几何算法:点集合构造简单多边形,布布扣,bubuko.com

几何算法:点集合构造简单多边形

标签:style   c   class   blog   code   tar   

原文地址:http://www.cnblogs.com/andyzeng/p/3754005.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!