songsenand
/
cgal
mirror of https://github.com/CGAL/cgal
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

*** empty log message ***

This commit is contained in:
Matthias Bäsken 2003-01-03 15:29:54 +00:00
parent d4e650847c
commit b0ffdcec2e
5 changed files with 66 additions and 62 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
Packages/GeoWin/changes.txt
View File

@ -1,3 +1,6 @@
1.6.4 (3 Jan 2003)
- arrangement and planar map demos changed
1.6.3 (27 Nov 2002)
- IRIX' make is picky about spaces.

37
Packages/GeoWin/demo/GeoWin/CGAL_arrangement_2.C
View File

@ -53,13 +53,18 @@ int main(int argc, char *argv[])
#include <CGAL/Cartesian.h>
#include <CGAL/Arrangement_2.h>
#include <CGAL/Arr_2_default_dcel.h>
#include <CGAL/Arr_leda_segment_exact_traits.h>
#include <CGAL/Pm_walk_along_line_point_location.h>
#include <CGAL/Arr_segment_traits_2.h>
#include <CGAL/leda_rational.h>
typedef CGAL::Arr_leda_segment_exact_traits Traits;
#if defined(LEDA_NAMESPACE)
using namespace leda;
#endif
typedef leda_rational NT;
typedef CGAL::Cartesian<NT> K;
typedef CGAL::Arr_segment_exact_traits<K> Traits;
typedef Traits::Point IPoint;
typedef CGAL::Cartesian<double> K;
typedef K::Point_2 Point;
typedef K::Segment_2 Segment;
typedef Traits::X_curve X_curve;
@ -68,16 +73,12 @@ typedef CGAL::Arr_base_node<X_curve> Base_node;
typedef CGAL::Arr_2_default_dcel<Traits> Dcel;
typedef CGAL::Arrangement_2<Dcel,Traits,Base_node > Arr_2;
typedef CGAL::Arrangement_2<Dcel,Traits,Base_node >::Curve_iterator Curve_iterator;
typedef CGAL::Arrangement_2<Dcel,Traits,Base_node >::Vertex_iterator Vertex_iterator;
typedef CGAL::Arrangement_2<Dcel,Traits,Base_node >::Vertex Vertex;
typedef Arr_2::Curve_iterator Curve_iterator;
typedef Arr_2::Vertex_iterator Vertex_iterator;
typedef Arr_2::Vertex Vertex;
#include <CGAL/geowin_support.h>
#if defined(LEDA_NAMESPACE)
using namespace leda;
#endif
Arr_2 arr;
class geo_locate : public geowin_update<std::list<Point>,std::list<Segment> >
@ -108,9 +109,7 @@ public:
if (f->does_outer_ccb_exist()) {
Arr_2::Ccb_halfedge_circulator cc=f->outer_ccb();
do {
leda_segment seg = (cc->curve()).to_float();
Segment cseg(Point(seg.xcoord1(),seg.ycoord1()), Point(seg.xcoord2(),seg.ycoord2()) );
Segment cseg = cc->curve();
Sl.push_back(cseg);
} while (++cc != f->outer_ccb());
@ -120,9 +119,7 @@ public:
for (;hit!=eit; ++hit) {
Arr_2::Ccb_halfedge_circulator cc=*hit;
do {
leda_segment seg = (cc->curve()).to_float();
Segment cseg(Point(seg.xcoord1(),seg.ycoord1()), Point(seg.xcoord2(),seg.ycoord2()) );
Segment cseg = cc->curve();
Sl.push_back(cseg);
} while (++cc != *hit);
}
@ -138,8 +135,7 @@ public:
bool insert(const Segment& seg)
{
leda_rat_segment ls(convert_to_leda(seg));
arr.insert(X_curve(ls.start(), ls.end() ));
arr.insert(X_curve(seg));
//std::cout << "insert:" << seg << "\n";
return true;
}
@ -167,8 +163,7 @@ public:
for(;it1 != L.end(); ++it1){
const Segment& s1= *it1;
leda_rat_segment ls(convert_to_leda(s1));
arr.insert(X_curve(ls.start(), ls.end() ));
arr.insert(s1);
}
}
};

82
Packages/GeoWin/demo/GeoWin/CGAL_planar_map_2.C
View File

@ -51,40 +51,37 @@ int main(int argc, char *argv[])
#else
#include <CGAL/Cartesian.h>
#include <CGAL/Pm_segment_epsilon_traits.h>
#include <CGAL/Pm_segment_traits_2.h>
#include <CGAL/Pm_default_dcel.h>
#include <CGAL/Planar_map_2.h>
#include <CGAL/convex_hull_2.h>
#include <CGAL/Polygon_2_algorithms.h>
#include <CGAL/Pm_with_intersections.h>
#include<list>
#include <CGAL/geowin_support.h>
#include <CGAL/distance_predicates_2.h>
#include <CGAL/leda_rational.h>
#if defined(LEDA_NAMESPACE)
using namespace leda;
#endif
typedef double number_type;
//typedef double number_type;
typedef leda_rational number_type;
typedef CGAL::Cartesian<number_type> K;
typedef CGAL::Pm_segment_epsilon_traits<K> pmtraits;
typedef CGAL::Pm_segment_traits_2<K> pmtraits;
typedef K::Point_2 Point;
typedef K::Segment_2 Segment;
typedef CGAL::Polygon_traits_2<K> PTraits;
typedef CGAL::Polygon_2<PTraits,std::list<Point> > Polygon;
typedef pmtraits::X_curve Curve;
typedef CGAL::Pm_default_dcel<pmtraits> pmdcel;
typedef CGAL::Planar_map_2<pmdcel, pmtraits> planar_map;
typedef planar_map::Locate_type loc_type;
typedef planar_map::Halfedge_handle halfedge_handle;
typedef CGAL::Planar_map_2<pmdcel, pmtraits> plmap;
typedef plmap::Locate_type loc_type;
typedef plmap::Halfedge_handle halfedge_handle;
std::list<Point> L;
std::list<Point> shoot;
std::list<Segment> Lseg;
std::list<Point> Loc;
bool snap_to_points(Segment& seg)
{
if (L.size() < 2) return false;
@ -122,7 +119,8 @@ bool segment_start_change(GeoWin& gw, const Segment& segold)
bool point_start_change(GeoWin& gw, const Point& p)
{ return false; }
CGAL::Planar_map_2<pmdcel, pmtraits>* pm;
//plmap pm; - seems not to work
plmap* pm;
class geo_shoot : public geowin_redraw, public geowin_update<std::list<Point>, std::list<Point> >
{
@ -145,13 +143,14 @@ public:
//draw shooting arrow ...
Point src = *pit;
double x1,y1,x2,y2,m,n;
x1 = (*it).source().x(); x2 = (*it).target().x();
y1 = (*it).source().y(); y2 = (*it).target().y();
x1 = CGAL::to_double ((*it).source().x()); x2 = CGAL::to_double ((*it).target().x());
y1 = CGAL::to_double ((*it).source().y()); y2 = CGAL::to_double ((*it).target().y());
m = (y2-y1)/(x2-x1);
n = y2-m*x2;
W.set_color(leda_red);
W.draw_arrow((*pit).x(),(*pit).y(),(*pit).x(),m*((*pit).x())+n);
W.draw_arrow(CGAL::to_double((*pit).x()),CGAL::to_double((*pit).y()),
CGAL::to_double((*pit).x()),CGAL::to_double(m*((*pit).x())+n));
W.set_color(leda_blue);
}
@ -169,6 +168,10 @@ public:
halfedge_handle he = pm->vertical_ray_shoot(*it,lt,true);
switch(lt){
case 1: { //Vertex
Point pt1 = (*he).source()->point();
Point pt2 = (*he).target()->point();
PSEG.push_back(Segment(pt1,pt2));
PT.push_back(*it);
break;
}
case 2: { //Edge
@ -208,11 +211,14 @@ public:
std::list<Point>::const_iterator rt;
W.set_color(leda_red);
for (rt=red.begin(); rt != red.end(); rt++) W.draw_disc((*rt).x(),(*rt).y(),5);
for (rt=red.begin(); rt != red.end(); rt++)
W.draw_disc(CGAL::to_double((*rt).x()),CGAL::to_double((*rt).y()),5);
W.set_color(leda_green);
for (rt=green.begin(); rt != green.end(); rt++) W.draw_disc((*rt).x(),(*rt).y(),5);
for (rt=green.begin(); rt != green.end(); rt++)
W.draw_disc(CGAL::to_double((*rt).x()),CGAL::to_double((*rt).y()),5);
W.set_color(leda_white);
for (rt=white.begin(); rt != white.end(); rt++) W.draw_disc((*rt).x(),(*rt).y(),5);
for (rt=white.begin(); rt != white.end(); rt++)
W.draw_disc(CGAL::to_double((*rt).x()),CGAL::to_double((*rt).y()),5);
}
virtual bool write_postscript(ps_file& PS, leda_color c1, leda_color c2)
@ -225,14 +231,17 @@ public:
if (pol.is_simple()) PS << convert_to_leda(pol);
}
// output the lists for location ...
std::list<CGALPoint>::const_iterator rt;
std::list<Point>::const_iterator rt;
PS.set_color(leda_red);
for (rt=red.begin(); rt != red.end(); rt++) PS.draw_disc((*rt).x(),(*rt).y(),5);
for (rt=red.begin(); rt != red.end(); rt++)
PS.draw_disc(CGAL::to_double((*rt).x()),CGAL::to_double((*rt).y()),5);
PS.set_color(leda_green);
for (rt=green.begin(); rt != green.end(); rt++) PS.draw_disc((*rt).x(),(*rt).y(),5);
for (rt=green.begin(); rt != green.end(); rt++)
PS.draw_disc(CGAL::to_double((*rt).x()),CGAL::to_double((*rt).y()),5);
PS.set_color(leda_white);
for (rt=white.begin(); rt != white.end(); rt++) PS.draw_disc((*rt).x(),(*rt).y(),5);
for (rt=white.begin(); rt != white.end(); rt++)
PS.draw_disc(CGAL::to_double((*rt).x()),CGAL::to_double((*rt).y()),5);
return false;
}
@ -292,40 +301,37 @@ public:
}
};
int main()
{
geowin_init_default_type((std::list<Point>*)0, leda_string("CGALPointList"));
geowin_init_default_type((std::list<Segment>*)0, leda_string("CGALSegmentList"));
CGAL::Planar_map_2<pmdcel, pmtraits> pmap;
plmap pmap;
pm = &pmap;
GeoWin GW("CGAL - Planar map demo");
GW.add_help_text(leda_string("CGAL_planar_map_2"));
GW.message("Attention - the input segments must not intersect");
geo_scene point_scene = GW.new_scene(L);
GeoEditScene<std::list<Point> >* PTR1 =
(GeoEditScene<std::list<Point> >*)point_scene;
GW.set_start_change_handler(PTR1, point_start_change);
GeoEditScene<std::list<Point> >* point_scene = GW.new_scene(L);
GW.set_name(point_scene,"Start/end points of line segments for the planar map");
GW.set_start_change_handler(point_scene, point_start_change);
// vertical ray shooting scene ...
geo_scene shoot_scene = GW.new_scene(shoot);
GW.set_name(shoot_scene,"Input for vertical ray shooting");
GW.set_color(shoot_scene, leda_blue);
GW.set_point_style(shoot_scene, leda_disc_point);
geo_scene segment_scene = GW.new_scene(Lseg);
GeoEditScene<std::list<Segment> >* PTR2 =
(GeoEditScene<std::list<Segment> >*)segment_scene;
GW.set_pre_add_change_handler(PTR2,segment_add_changer);
GW.set_start_change_handler(PTR2,segment_start_change);
GeoEditScene<std::list<Segment> >* segment_scene = GW.new_scene(Lseg);
GW.set_name(segment_scene,"Input segments of the planar map");
GW.set_pre_add_change_handler(segment_scene,segment_add_changer);
GW.set_start_change_handler(segment_scene,segment_start_change);
// second point scene for locating ...
geo_scene ploc_scene = GW.new_scene(Loc);
geo_scene ploc_scene = GW.new_scene(Loc);
GW.set_name(ploc_scene,"Input for point location");
GW.set_color(ploc_scene,leda_black);
geo_plmap PM;

4
Packages/GeoWin/demo/GeoWin/CGAL_planar_map_2.hlp
View File

@ -3,9 +3,9 @@ CGAL Planar map demo
This demo program lets you construct a planar map.
Construct the planar map using the scenes
- CGALPointList and
- Start/end points of line segments for the planar map and
- CGALSegmentList.
- Input segments of the planar map
The program will show you the faces of the planar map in
different colors.

2
Packages/GeoWin/version
View File

@ -1 +1 @@
1.6.3 (27 Nov 2002)
1.6.4 (3 Jan 2003)