void draw(cairo_t *cr,
std::ostringstream *notify,
int width, int height, bool save, std::ostringstream *timer_stream) {
if (first_time)
{
first_time = false;
sliders[0].geometry(Point(50, 50), 100);
}
size_t const num_points = static_cast<size_t>(sliders[0].value());
D2<SBasis> B1 = b_handle.asBezier();
Piecewise<D2<SBasis> >B;
B.concat(Piecewise<D2<SBasis> >(B1));
// testing fuse_nearby_ends
std::vector< Piecewise<D2<SBasis> > > pieces;
pieces = fuse_nearby_ends(split_at_discontinuities(B),9);
Piecewise<D2<SBasis> > C;
for (unsigned i=0; i<pieces.size(); i++){
C.concat(pieces[i]);
}
// testing fuse_nearby_ends
cairo_set_line_width (cr, 2.);
cairo_set_source_rgba (cr, 0., 0.5, 0., 1);
//cairo_d2_sb(cr, B1);
//cairo_pw_d2_sb(cr, C);
//cairo_pw_d2_sb(cr, B);
cairo_stroke(cr);
Timer tm;
Timer::Time als_time = tm.lap();
cairo_set_source_rgba (cr, 0., 0., 0.9, 1);
//dot_plot(cr,uniform_B);
cairo_stroke(cr);
std::cout << B[0] << std::endl;
Geom::Affine translation;
Geom::Path original_path;
//original_bezier.append(B[0]);
//original_bezier.appendNew<CubicBezier> (B[0]);
CubicBezier original_bezier(b_handle.pts);
original_path.append(original_bezier);
std::vector<double> initial_t;
std::vector<Geom::Point> curve_points;
if (randomize_times) {
std::uniform_real_distribution<double> dist_t(0,1);
for (size_t ii = 0; ii < num_points; ++ii) {
double const t = dist_t(generator);
initial_t.push_back(t);
}
std::sort(initial_t.begin(), initial_t.end());
double const min = initial_t.front();
double const max = initial_t.back();
for (auto& t : initial_t) {
t = (t-min)/(max-min);
}
for (auto const t : initial_t) {
curve_points.push_back(original_bezier.pointAt(t));
}
}
else {
for (size_t ii = 0; ii < num_points; ++ii) {
double const t = static_cast<double>(ii) / (num_points-1);
Geom::Point const p = original_bezier.pointAt(t);
initial_t.push_back(t);
curve_points.push_back(p);
}
}
cairo_set_source_rgba (cr, 0., 0., .9, 1);
cairo_path(cr, original_path);
draw_text(cr, original_path.initialPoint(), "original curve and old fit");
Geom::CubicBezier fitted_new;
Geom::CubicBezier fitted_new_a;
Geom::Point very_old_version_raw[4];
bezier_fit_cubic(very_old_version_raw, curve_points.data(), curve_points.size(), 0.);
Geom::CubicBezier very_old_bezier(
very_old_version_raw[0],
very_old_version_raw[1],
very_old_version_raw[2],
very_old_version_raw[3]
);
Geom::Path very_old_version_path;
very_old_version_path.append(very_old_bezier);
cairo_set_source_rgba (cr, .7, .7, 0., 1);
cairo_stroke(cr);
cairo_path(cr, very_old_version_path);
cairo_set_source_rgba (cr, 0., 0., .9, 1);
cairo_stroke(cr);
cross_plot(cr, curve_points);
if(1) {
Geom::CubicBezier combination(very_old_bezier);
tm.ask_for_timeslice();
tm.start();
auto new_result_ig_a = experiment::fit_bezier(combination, curve_points);
als_time = tm.lap();
*notify << "Bezier fit a, old algorithm as initial guess, time = " << als_time << std::endl
<< "Worst residual: " << new_result_ig_a.first << " at t=" << new_result_ig_a.second << std::endl;
Geom::Path combination_path;
translation.setTranslation(Geom::Point(300,300));
combination_path.append(combination.transformed(translation));
cairo_set_source_rgba (cr, .0, .0, .9, 1);
cross_plot(cr, curve_points, translation.translation());
cairo_path(cr, combination_path);
draw_text(cr, combination_path.initialPoint(), "old fit as i.g.");
}
{
tm.ask_for_timeslice();
tm.start();
auto new_result = fit_bezier(fitted_new, curve_points);
als_time = tm.lap();
*notify << "Bezier fit, time = " << als_time << std::endl
<< "Worst residual: " << new_result.first << " at t=" << new_result.second << std::endl;
Geom::Path fitted_new_path;
translation.setTranslation(Geom::Point(300,0));
fitted_new_path.append(fitted_new.transformed(translation));
cairo_set_source_rgba (cr, .0, .9, .0, 1);
cross_plot(cr, curve_points, translation.translation());
cairo_path(cr, fitted_new_path);
draw_text(cr, fitted_new_path.initialPoint(), "new fit");
}
{
tm.ask_for_timeslice();
tm.start();
auto new_result_a = experiment::fit_bezier(fitted_new_a, curve_points);
als_time = tm.lap();
*notify << "Bezier fit a, time = " << als_time << std::endl
<< "Worst residual: " << new_result_a.first << " at t=" << new_result_a.second << std::endl;
Geom::Path fitted_new_a_path;
translation.setTranslation(Geom::Point(0,300));
fitted_new_a_path.append(fitted_new_a.transformed(translation));
cairo_set_source_rgba (cr, .9, .0, .0, 1);
cross_plot(cr, curve_points, translation.translation());
cairo_path(cr, fitted_new_a_path);
draw_text(cr, fitted_new_a_path.initialPoint(), "new fit (a)");
}
Geom::CubicBezier fixed_times_bezier;
{
tm.ask_for_timeslice();
tm.start();
auto fixed_times_result = experiment::fit_bezier_fixed_times(fixed_times_bezier, curve_points);
als_time = tm.lap();
*notify << "Bezier fit a (fixed times), time = " << als_time << std::endl
<< "Worst residual: " << fixed_times_result.first << " at t=" << fixed_times_result.second << std::endl;
Geom::Path fixed_times_path;
translation.setTranslation(Geom::Point(600,300));
fixed_times_path.append(fixed_times_bezier.transformed(translation));
cairo_set_source_rgba (cr, .9, .0, .0, 1);
cross_plot(cr, curve_points, translation.translation());
cairo_path(cr, fixed_times_path);
draw_text(cr, fixed_times_path.initialPoint(), "fixed t fit (a)");
}
Geom::CubicBezier fixed_times_ig_bezier = fixed_times_bezier;
{
tm.ask_for_timeslice();
tm.start();
auto fixed_times_ig_result = experiment::fit_bezier(fixed_times_ig_bezier, curve_points);
als_time = tm.lap();
*notify << "Bezier fit a (with fixed times as i.g.), time = " << als_time << std::endl
<< "Worst residual: " << fixed_times_ig_result.first << " at t=" << fixed_times_ig_result.second << std::endl;
Geom::Path fixed_times_path;
translation.setTranslation(Geom::Point(900,300));
fixed_times_path.append(fixed_times_ig_bezier.transformed(translation));
cairo_set_source_rgba (cr, .9, .0, .0, 1);
cross_plot(cr, curve_points, translation.translation());
cairo_path(cr, fixed_times_path);
draw_text(cr, fixed_times_path.initialPoint(), "new (a) with fixed t as i.g.");
}
Geom::CubicBezier icp_bezier;
{
tm.ask_for_timeslice();
tm.start();
auto icp_result = experiment::fit_bezier_icp(icp_bezier, curve_points);
als_time = tm.lap();
*notify << "Bezier fit icp, time = " << als_time << std::endl
<< "Worst residual: " << icp_result.first << " at t=" << icp_result.second << std::endl;
Geom::Path icp_path;
translation.setTranslation(Geom::Point(600,600));
icp_path.append(icp_bezier.transformed(translation));
cairo_set_source_rgba (cr, .9, .0, .0, 1);
cross_plot(cr, curve_points, translation.translation());
cairo_path(cr, icp_path);
draw_text(cr, icp_path.initialPoint(), "icp fit");
}
Geom::CubicBezier icp_ig_bezier(icp_bezier);
{
tm.ask_for_timeslice();
tm.start();
auto icp_ig_result = experiment::fit_bezier(icp_ig_bezier, curve_points);
als_time = tm.lap();
*notify << "Bezier fit with icp i.g., time = " << als_time << std::endl
<< "Worst residual: " << icp_ig_result.first << " at t=" << icp_ig_result.second << std::endl;
Geom::Path icp_ig_path;
translation.setTranslation(Geom::Point(900,600));
icp_ig_path.append(icp_ig_bezier.transformed(translation));
cairo_set_source_rgba (cr, .9, .0, .0, 1);
cross_plot(cr, curve_points, translation.translation());
cairo_path(cr, icp_ig_path);
draw_text(cr, icp_ig_path.initialPoint(), "bezier fit with icp as i.g.");
}
std::cout << "original: " << write_svg_path(original_path) << std::endl;
Geom::CubicBezier initial_guess(
curve_points.front(), curve_points.front(),
curve_points.back(), curve_points.back()
);
{
experiment::get_initial_guess(initial_guess, curve_points);
Geom::Path initial_guess_path;
translation.setTranslation(Geom::Point(600,0));
initial_guess_path.append(initial_guess.transformed(translation));
cairo_set_source_rgba (cr, .8, .0, .8, 1);
cross_plot(cr, curve_points, translation.translation());
cairo_path(cr, initial_guess_path);
draw_text(cr, initial_guess_path.initialPoint(), "initial guess");
}
cairo_stroke(cr);
Toy::draw(cr, notify, width, height, save,timer_stream);
}
Geom::OptRect SPShape::bbox(Geom::Affine const &transform, SPItem::BBoxType bboxtype) const {
Geom::OptRect bbox;
if (!this->_curve) {
return bbox;
}
bbox = bounds_exact_transformed(this->_curve->get_pathvector(), transform);
if (!bbox) {
return bbox;
}
if (bboxtype == SPItem::VISUAL_BBOX) {
// convert the stroke to a path and calculate that path's geometric bbox
if (!this->style->stroke.isNone()) {
Geom::PathVector *pathv = item_outline(this, true); // calculate bbox_only
if (pathv) {
bbox |= bounds_exact_transformed(*pathv, transform);
delete pathv;
}
}
// Union with bboxes of the markers, if any
if ( this->hasMarkers() && !this->_curve->get_pathvector().empty() ) {
/** \todo make code prettier! */
Geom::PathVector const & pathv = this->_curve->get_pathvector();
// START marker
for (unsigned i = 0; i < 2; i++) { // SP_MARKER_LOC and SP_MARKER_LOC_START
if ( this->_marker[i] ) {
SPItem* marker_item = sp_item_first_item_child( _marker[i] );
if (marker_item) {
Geom::Affine tr(sp_shape_marker_get_transform_at_start(pathv.begin()->front()));
if (_marker[i]->orient_mode == MARKER_ORIENT_AUTO_START_REVERSE) {
// Reverse start marker if necessary
tr = Geom::Rotate::from_degrees( 180.0 ) * tr;
} else if (_marker[i]->orient_mode == MARKER_ORIENT_ANGLE) {
Geom::Point transl = tr.translation();
tr = Geom::Rotate::from_degrees(_marker[i]->orient) * Geom::Translate(transl);
}
if (_marker[i]->markerUnits == SP_MARKER_UNITS_STROKEWIDTH) {
tr = Geom::Scale(this->style->stroke_width.computed) * tr;
}
// total marker transform
tr = marker_item->transform * _marker[i]->c2p * tr * transform;
// get bbox of the marker with that transform
bbox |= marker_item->visualBounds(tr);
}
}
}
// MID marker
for (unsigned i = 0; i < 3; i += 2) { // SP_MARKER_LOC and SP_MARKER_LOC_MID
if ( !this->_marker[i] ) {
continue;
}
SPMarker* marker = _marker[i];
SPItem* marker_item = sp_item_first_item_child( marker );
if ( !marker_item ) {
continue;
}
for(Geom::PathVector::const_iterator path_it = pathv.begin(); path_it != pathv.end(); ++path_it) {
// START position
if ( path_it != pathv.begin()
&& ! ((path_it == (pathv.end()-1)) && (path_it->size_default() == 0)) ) // if this is the last path and it is a moveto-only, there is no mid marker there
{
Geom::Affine tr(sp_shape_marker_get_transform_at_start(path_it->front()));
if (marker->orient_mode == MARKER_ORIENT_ANGLE) {
Geom::Point transl = tr.translation();
tr = Geom::Rotate::from_degrees(marker->orient) * Geom::Translate(transl);
}
if (marker->markerUnits == SP_MARKER_UNITS_STROKEWIDTH) {
tr = Geom::Scale(this->style->stroke_width.computed) * tr;
}
tr = marker_item->transform * marker->c2p * tr * transform;
bbox |= marker_item->visualBounds(tr);
}
// MID position
if ( path_it->size_default() > 1) {
Geom::Path::const_iterator curve_it1 = path_it->begin(); // incoming curve
Geom::Path::const_iterator curve_it2 = ++(path_it->begin()); // outgoing curve
while (curve_it2 != path_it->end_default())
{
/* Put marker between curve_it1 and curve_it2.
* Loop to end_default (so including closing segment), because when a path is closed,
* there should be a midpoint marker between last segment and closing straight line segment */
SPMarker* marker = _marker[i];
SPItem* marker_item = sp_item_first_item_child( marker );
if (marker_item) {
Geom::Affine tr(sp_shape_marker_get_transform(*curve_it1, *curve_it2));
if (marker->orient_mode == MARKER_ORIENT_ANGLE) {
Geom::Point transl = tr.translation();
tr = Geom::Rotate::from_degrees(marker->orient) * Geom::Translate(transl);
}
if (marker->markerUnits == SP_MARKER_UNITS_STROKEWIDTH) {
tr = Geom::Scale(this->style->stroke_width.computed) * tr;
}
tr = marker_item->transform * marker->c2p * tr * transform;
bbox |= marker_item->visualBounds(tr);
}
++curve_it1;
++curve_it2;
}
}
// END position
if ( path_it != (pathv.end()-1) && !path_it->empty()) {
Geom::Curve const &lastcurve = path_it->back_default();
Geom::Affine tr = sp_shape_marker_get_transform_at_end(lastcurve);
if (marker->orient_mode == MARKER_ORIENT_ANGLE) {
Geom::Point transl = tr.translation();
tr = Geom::Rotate::from_degrees(marker->orient) * Geom::Translate(transl);
}
if (marker->markerUnits == SP_MARKER_UNITS_STROKEWIDTH) {
tr = Geom::Scale(this->style->stroke_width.computed) * tr;
}
tr = marker_item->transform * marker->c2p * tr * transform;
bbox |= marker_item->visualBounds(tr);
}
}
}
// END marker
for (unsigned i = 0; i < 4; i += 3) { // SP_MARKER_LOC and SP_MARKER_LOC_END
if ( _marker[i] ) {
SPMarker* marker = _marker[i];
SPItem* marker_item = sp_item_first_item_child( marker );
if (marker_item) {
/* Get reference to last curve in the path.
* For moveto-only path, this returns the "closing line segment". */
Geom::Path const &path_last = pathv.back();
unsigned int index = path_last.size_default();
if (index > 0) {
index--;
}
Geom::Curve const &lastcurve = path_last[index];
Geom::Affine tr = sp_shape_marker_get_transform_at_end(lastcurve);
if (marker->orient_mode == MARKER_ORIENT_ANGLE) {
Geom::Point transl = tr.translation();
tr = Geom::Rotate::from_degrees(marker->orient) * Geom::Translate(transl);
}
if (marker->markerUnits == SP_MARKER_UNITS_STROKEWIDTH) {
tr = Geom::Scale(this->style->stroke_width.computed) * tr;
}
// total marker transform
tr = marker_item->transform * marker->c2p * tr * transform;
// get bbox of the marker with that transform
bbox |= marker_item->visualBounds(tr);
}
}
}
}
}
return bbox;
}
