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

Implement insertion of conics by specifying 3 or 5 points.

This commit is contained in:
Alex Tsui 2012-06-29 04:12:55 +00:00
parent 1a5c7a49dc
commit ebb6cf73fe
3 changed files with 171 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
Arrangement_on_surface_2/demo/Arrangement_on_surface_2/ArrangementDemoWindow.cpp
View File

@ -107,11 +107,9 @@ setupUi( )
this->conicTypeGroup = new QActionGroup( this );
this->conicTypeGroup->addAction( this->ui->actionConicSegment );
this->conicTypeGroup->addAction( this->ui->actionConicCircle );
#if 0
this->conicTypeGroup->addAction( this->ui->actionConicEllipse );
this->conicTypeGroup->addAction( this->ui->actionConicThreePoint );
this->conicTypeGroup->addAction( this->ui->actionConicFivePoint );
#endif
}
void
@ -341,6 +339,18 @@ updateConicType( QAction* newType )
{
curveInputCallback->setConicType( ConicCurveInputCallback::CONIC_CIRCLE );
}
else if ( newType == this->ui->actionConicEllipse )
{
curveInputCallback->setConicType( ConicCurveInputCallback::CONIC_ELLIPSE );
}
else if ( newType == this->ui->actionConicThreePoint )
{
curveInputCallback->setConicType( ConicCurveInputCallback::CONIC_THREE_POINT );
}
else if ( newType == this->ui->actionConicFivePoint )
{
curveInputCallback->setConicType( ConicCurveInputCallback::CONIC_FIVE_POINT );
}
}
else
{

39
Arrangement_on_surface_2/demo/Arrangement_on_surface_2/ArrangementDemoWindow.ui
View File

@ -114,6 +114,9 @@
<addaction name="separator"/>
<addaction name="actionConicSegment"/>
<addaction name="actionConicCircle"/>
<addaction name="actionConicEllipse"/>
<addaction name="actionConicThreePoint"/>
<addaction name="actionConicFivePoint"/>
</widget>
<action name="actionQuit">
<property name="text">
@ -334,6 +337,42 @@
<string>Circle</string>
</property>
</action>
<action name="actionConicEllipse">
<property name="checkable">
<bool>true</bool>
</property>
<property name="icon">
<iconset resource="ArrangementDemoWindow.qrc">
<normaloff>:/icons/demo_conic_ellipse.xpm</normaloff>:/icons/demo_conic_ellipse.xpm</iconset>
</property>
<property name="text">
<string>Ellipse</string>
</property>
</action>
<action name="actionConicThreePoint">
<property name="checkable">
<bool>true</bool>
</property>
<property name="icon">
<iconset resource="ArrangementDemoWindow.qrc">
<normaloff>:/icons/demo_conic_3points.xpm</normaloff>:/icons/demo_conic_3points.xpm</iconset>
</property>
<property name="text">
<string>Three Point Conic</string>
</property>
</action>
<action name="actionConicFivePoint">
<property name="checkable">
<bool>true</bool>
</property>
<property name="icon">
<iconset resource="ArrangementDemoWindow.qrc">
<normaloff>:/icons/demo_conic_5points.xpm</normaloff>:/icons/demo_conic_5points.xpm</iconset>
</property>
<property name="text">
<string>Five Point Conic</string>
</property>
</action>
</widget>
<resources>
<include location="ArrangementDemoWindow.qrc"/>

124
Arrangement_on_surface_2/demo/Arrangement_on_surface_2/GraphicsViewCurveInput.h
View File

@ -11,6 +11,7 @@
#include <iostream>
#include "Callback.h"
#include "ISnappable.h"
#include <CGAL/CORE_algebraic_number_traits.h>
namespace CGAL {
namespace Qt {
@ -228,6 +229,7 @@ public:
typedef AlgKernel Kernel;
//typedef typename Kernel::Point_2 Point_2;
typedef typename Kernel::Segment_2 Segment_2;
typedef typename RatKernel::FT Rat_FT;
typedef typename RatKernel::Point_2 Rat_point_2;
typedef typename RatKernel::Segment_2 Rat_segment_2;
typedef typename RatKernel::Circle_2 Rat_circle_2;
@ -244,7 +246,8 @@ public:
GraphicsViewCurveInputBase( parent ),
construct_x_monotone_curve_2( this->traits.construct_x_monotone_curve_2_object( ) ),
conicType( CONIC_SEGMENT ),
circleItem( NULL )
circleItem( NULL ),
ellipseItem( NULL )
{ }
void setConicType( ConicType conicType_ )
@ -281,6 +284,15 @@ protected:
double d = radius * sqrt( 2.0 );
this->circleItem->setRect( p1.x( ) - radius, p1.y( ) - radius, 2*radius, 2*radius );
}
if ( this->ellipseItem != NULL )
{
Point_2 p1 = this->points.back( );
Point_2 p2 = this->snapPoint( event );
CGAL::Bbox_2 bb = p1.bbox( ) + p2.bbox( );
double w = bb.xmax( ) - bb.xmin( );
double h = bb.ymax( ) - bb.ymin( );
this->ellipseItem->setRect( bb.xmin( ), bb.ymin( ), w, h );
}
}
void mousePressEvent( QGraphicsSceneMouseEvent* event )
@ -307,8 +319,17 @@ protected:
QPointF pt = this->convert( clickedPoint );
if ( this->scene != NULL )
{
QGraphicsEllipseItem* circle = this->scene->addEllipse( pt.x( ), pt.y( ), 0, 0 );
this->circleItem = circle;
QGraphicsEllipseItem* ellipse = this->scene->addEllipse( pt.x( ), pt.y( ), 0, 0 );
this->circleItem = ellipse;
}
}
else if ( this->conicType == CONIC_ELLIPSE )
{
QPointF pt = this->convert( clickedPoint );
if ( this->scene != NULL )
{
QGraphicsEllipseItem* ellipse = this->scene->addEllipse( pt.x( ), pt.y( ), 0, 0 );
this->ellipseItem = ellipse;
}
}
}
@ -316,7 +337,6 @@ protected:
{
if ( this->conicType == CONIC_SEGMENT )
{
for ( int i = 0; i < this->polylineGuide.size( ); ++i )
{
if ( this->scene != NULL )
@ -357,6 +377,101 @@ protected:
this->points.clear( );
emit generate( CGAL::make_object( res ) );
}
else if ( this->conicType == CONIC_ELLIPSE )
{
if ( this->scene != NULL )
{
this->scene->removeItem( this->ellipseItem );
}
this->ellipseItem = NULL;
CGAL::Bbox_2 bb = this->points[ 0 ].bbox( ) + this->points[ 1 ].bbox( );
double x1 = CGAL::to_double( bb.xmin( ) );
double y1 = CGAL::to_double( bb.ymin( ) );
double x2 = CGAL::to_double( bb.xmax( ) );
double y2 = CGAL::to_double( bb.ymax( ) );
double sq_rad = CGAL::square(x2 - x1) + CGAL::square(y2 - y1);
Rat_FT a = CORE::abs( Rat_FT(x1) - Rat_FT(x2) )/2;
Rat_FT b = CORE::abs( Rat_FT(y1) - Rat_FT(y2) )/2;
Rat_FT a_sq = a*a;
Rat_FT b_sq = b*b;
Rat_FT x0 = (x2 + x1)/2;
Rat_FT y0 = (y2 + y1)/2;
Rat_FT r = b_sq;
Rat_FT s = a_sq;
Rat_FT t = 0;
Rat_FT u = -2*x0*b_sq;
Rat_FT v = -2*y0*a_sq;
Rat_FT ww = x0*x0*b_sq + y0*y0*a_sq - a_sq*b_sq;
Curve_2 res = Curve_2( r, s, t, u, v, ww );
this->points.clear( );
emit generate( CGAL::make_object( res ) );
}
else if ( this->conicType == CONIC_THREE_POINT )
{
if ( this->points.size( ) == 3 )
{
QPointF qp1 = this->convert( this->points[ 0 ] );
QPointF qp2 = this->convert( this->points[ 1 ] );
QPointF qp3 = this->convert( this->points[ 2 ] );
Rat_point_2 p1 = Rat_point_2( qp1.x( ), qp1.y( ) );
Rat_point_2 p2 = Rat_point_2( qp2.x( ), qp2.y( ) );
Rat_point_2 p3 = Rat_point_2( qp3.x( ), qp3.y( ) );
RatKernel ker;
if ( ! ker.collinear_2_object()( p1, p2, p3 ) )
{
Curve_2 res( p1, p2, p3 );
emit generate( CGAL::make_object( res ) );
}
else
{
std::cout << "Oops, points don't specify a valid conic. Try again!" << std::endl;
}
// TODO: make a valid curve and insert it
this->points.clear( );
}
}
else if ( this->conicType == CONIC_FIVE_POINT )
{
if ( this->points.size( ) == 5 )
{
QPointF qp1 = this->convert( this->points[ 0 ] );
QPointF qp2 = this->convert( this->points[ 1 ] );
QPointF qp3 = this->convert( this->points[ 2 ] );
QPointF qp4 = this->convert( this->points[ 3 ] );
QPointF qp5 = this->convert( this->points[ 4 ] );
Rat_point_2 p1 = Rat_point_2( qp1.x( ), qp1.y( ) );
Rat_point_2 p2 = Rat_point_2( qp2.x( ), qp2.y( ) );
Rat_point_2 p3 = Rat_point_2( qp3.x( ), qp3.y( ) );
Rat_point_2 p4 = Rat_point_2( qp4.x( ), qp4.y( ) );
Rat_point_2 p5 = Rat_point_2( qp5.x( ), qp5.y( ) );
try
{
Curve_2 res( p1, p2, p3, p4, p5 );
if ( res.is_valid( ) )
{
std::cout << "Meow, valid 5 point conic!" << std::endl;
emit generate( CGAL::make_object( res ) );
}
else
{
std::cout << "Oops, points don't specify a valid conic. Try again!" << std::endl;
}
this->points.clear( );
}
catch (...)
{
std::cout << "Oops, points don't specify a valid conic. Try again!" << std::endl;
this->points.clear( );
}
}
}
}
}
@ -371,6 +486,7 @@ protected:
std::vector< Point_2 > points;
std::vector< QGraphicsLineItem* > polylineGuide;
QGraphicsEllipseItem* circleItem;
QGraphicsEllipseItem* ellipseItem;
Traits traits;
Construct_x_monotone_curve_2 construct_x_monotone_curve_2;