开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 数值算法/人工智能 > 查看源码
line_arrange.c
资源名称：leda.tar.gz　[点击查看]
上传用户：gzelex
上传日期：2007-01-07
资源大小：707k
文件大小：3k
源码类别：

数值算法/人工智能

开发平台：

MultiPlatform

line_arrange.c：源码内容
  1. #include <LEDA/subdivision.h>
  2. #include <LEDA/plane_alg.h>
  3. #include <LEDA/window.h>
  8. void  show_face(subdivision<segment>& G, face f, window& W)
  9. { list<node> L = G.adj_nodes(f);
  10.   list<point> P;
  11.   node v;
  12.   forall(v,L) P.append(G[v]);
  13.   W.draw_filled_polygon(P,red);
  14. }
  17. main()
  18. {
  20.   window W;
  22.   W.init(0,1000,0);
  24.   W.message("            Arrangement of line segments                ");
  25.   W.message("                                                        ");
  26.   W.message("This program computes an arrangement of straight lines. ");
  27.   W.message("Use the left button to insert a sequence of lines and ");
  28.   W.message("terminate the input by clicking the right button. Now");
  29.   W.message("you can input query points with the left button. For ");
  30.   W.message("each point p the face of the arrangment containing p ");
  31.   W.message("is computed and displayed. Terminate the program by  "); 
  32.   W.message("pressing the right mouse key.                        ");
  33.   W.message("                                                     ");
  34.   W.message("Click any button to start.                           ");
  35.   W.read_mouse();
  36.   W.clear();
  40.   list<segment> seglist;
  42.   // build up a frame
  44.   double x0 = W.xmin();
  45.   double x1 = W.xmax();
  46.   double y0 = W.ymin();
  47.   double y1 = W.ymax();
  49.   segment b(x0,y0,x1,y0);
  50.   segment t(x0,y1,x1,y1);
  51.   segment l(x0,y0,x0,y1);
  52.   segment r(x1,y0,x1,y1);
  54.   seglist.append(b);
  55.   seglist.append(t);
  56.   seglist.append(l);
  57.   seglist.append(r);
  59.   line L;
  61.   while ( W >> L ) 
  62.   {  W << L;
  63.      if (L.vertical()) 
  64.         seglist.append(segment(L.x_proj(y0),y0,L.x_proj(y1),y1));
  65.      else 
  66.         seglist.append(segment(x0,L.y_proj(x0),x1,L.y_proj(x1)));
  67.    }
  69.   W.message("Computing subdivision");
  70.   newline;
  72.   GRAPH<point,segment> G;
  74.   SWEEP_SEGMENTS(seglist,G);
  77.   // insert reverse edges
  79.   edge e;
  80.   list<edge> E = G.all_edges();
  81.   forall(e,E) G.new_edge(target(e),source(e),G[e]);
  83.   // sort edges counter-clockwise
  85.   edge_array<double> angle(G);
  87.   forall_edges(e,G)
  88.   { segment s(G[source(e)],G[target(e)]);
  89.     angle[e] = s.angle();
  90.    }
  92.   G.sort_edges(angle);
  94.   subdivision<segment> S(G);
  96.   W.clear(); 
  97.   W.set_line_width(2);
  99.   forall_edges(e,G)
  100.   W.draw_segment(G[source(e)], G[target(e)],blue);
  103.   // Locate points
  105.   W.set_mode(xor_mode);
  107.   point p;
  108.   face  f=nil;
  110.   W.message("Give query points !");
  112.   while (W >> p)
  113.   { if (f) show_face(S,f,W);
  114.     f = S.locate_point(p);
  115.     show_face(S,f,W);
  116.    }
  120. /*
  121.   // show strips
  123.   segment s;
  124.   slist<segment> SL;
  126.   while (W >> p)
  127.   { forall(s,SL) W.draw_segment(s);
  128.     SL.clear();
  129.     SL = S.compute_strip(p);
  130.     W.set_line_width(1);
  131.     W.draw_vline(p.xcoord());
  132.     W.set_line_width(5);
  133.     forall(s,SL) W.draw_segment(s);
  134.    }
  135. */;
  137.   return 0;
  138. }