Planar Subdivision in CSharp: Difference between revisions

From EMGU
Jump to navigation Jump to search
No edit summary
No edit summary
Line 1: Line 1:
'''This example requires [[Version_History#Emgu.CV-1.5.0.0|Emgu CV 1.5.0.0]]'''
[[image:PlanarSubdivisionExample.png]]
[[image:PlanarSubdivisionExample.png]]


Line 8: Line 10:
using Emgu.CV;
using Emgu.CV;
using Emgu.CV.UI;
using Emgu.CV.UI;
using Emgu.CV.Structure;


namespace PlannarSubdivision
namespace PlanarSubdivisionExample
{
{
   static class Program
   static class Program
Line 29: Line 32:


         #region create random points in the range of [0, maxValue]
         #region create random points in the range of [0, maxValue]
         Point2D<float>[] pts = new Point2D<float>[20];
         PointF[] pts = new PointF[20];
         Random r = new Random((int)(DateTime.Now.Ticks & 0x0000ffff));
         Random r = new Random((int)(DateTime.Now.Ticks & 0x0000ffff));
         for (int i = 0; i < pts.Length; i++)
         for (int i = 0; i < pts.Length; i++)
             pts[i] = new Point2D<float>((float)r.NextDouble() * maxValue, (float)r.NextDouble() * maxValue);
             pts[i] = new PointF((float)r.NextDouble() * maxValue, (float)r.NextDouble() * maxValue);
         #endregion
         #endregion


         List<Triangle2D<float>> delaunayTriangles;
         Triangle2DF[] delaunayTriangles;
         List<VoronoiFacet> voronoiFacets;
         VoronoiFacet[] voronoiFacets;
         using (PlanarSubdivision subdivision = new PlanarSubdivision(pts))
         using (PlanarSubdivision subdivision = new PlanarSubdivision(pts))
         {
         {
Line 52: Line 55:
         foreach (VoronoiFacet facet in voronoiFacets)
         foreach (VoronoiFacet facet in voronoiFacets)
         {
         {
             MCvPoint[] points = Array.ConvertAll<Point2D<float>, MCvPoint>(facet.Vertices, delegate(Point2D<float> p) { return p.MCvPoint; });
             Point[] points = Array.ConvertAll<PointF, Point>(facet.Vertices, Point.Round);


             //Draw the facet in color
             //Draw the facet in color
Line 64: Line 67:


             //draw the points associated with each facet in red
             //draw the points associated with each facet in red
             img.Draw(new Circle<float>(facet.Point, 5), new Bgr(Color.Red), 0);
             img.Draw(new CircleF(facet.Point, 5.0f), new Bgr(Color.Red), 0);
         }
         }


         //Draw the Delaunay triangulation
         //Draw the Delaunay triangulation
         foreach (Triangle2D<float> triangles in delaunayTriangles)
         foreach (Triangle2DF triangles in delaunayTriangles)
         {
         {
             img.Draw(triangles, new Bgr(Color.White), 1);
             img.Draw(triangles, new Bgr(Color.White), 1);

Revision as of 03:44, 25 February 2009

This example requires Emgu CV 1.5.0.0

using System;
using System.Collections.Generic;
using System.Windows.Forms;
using System.Drawing;
using Emgu.CV;
using Emgu.CV.UI;
using Emgu.CV.Structure;

namespace PlanarSubdivisionExample
{
   static class Program
   {
      /// <summary>
      /// The main entry point for the application.
      /// </summary>
      [STAThread]
      static void Main()
      {
         Application.EnableVisualStyles();
         Application.SetCompatibleTextRenderingDefault(false);
         Run();
      }

      static void Run()
      {
         float maxValue = 600;

         #region create random points in the range of [0, maxValue]
         PointF[] pts = new PointF[20];
         Random r = new Random((int)(DateTime.Now.Ticks & 0x0000ffff));
         for (int i = 0; i < pts.Length; i++)
            pts[i] = new PointF((float)r.NextDouble() * maxValue, (float)r.NextDouble() * maxValue);
         #endregion

         Triangle2DF[] delaunayTriangles;
         VoronoiFacet[] voronoiFacets;
         using (PlanarSubdivision subdivision = new PlanarSubdivision(pts))
         {
            //Obtain the delaunay's triangulation from the set of points;
            delaunayTriangles = subdivision.GetDelaunayTriangles();

            //Obtain the voronoi facets from the set of points
            voronoiFacets = subdivision.GetVoronoiFacets();
         }

         //create an image for display purpose
         Image<Bgr, Byte> img = new Image<Bgr, byte>((int)maxValue, (int) maxValue);

         //Draw the voronoi Facets
         foreach (VoronoiFacet facet in voronoiFacets)
         {
            Point[] points = Array.ConvertAll<PointF, Point>(facet.Vertices, Point.Round);

            //Draw the facet in color
            img.FillConvexPoly(
                points,
                new Bgr(r.NextDouble() * 120, r.NextDouble() * 120, r.NextDouble() * 120)
                );

            //highlight the edge of the facet in black
            img.DrawPolyline(points, true, new Bgr(Color.Black), 2);

            //draw the points associated with each facet in red
            img.Draw(new CircleF(facet.Point, 5.0f), new Bgr(Color.Red), 0);
         }

         //Draw the Delaunay triangulation
         foreach (Triangle2DF triangles in delaunayTriangles)
         {
            img.Draw(triangles, new Bgr(Color.White), 1);
         }

         //display the image
         ImageViewer.Show(img, "Plannar Subdivision");
      }
   }
}