2 Copyright (C) 2011-2013 Paul Davis
3 Author: Carl Hetherington <cth@carlh.net>
5 This program is free software; you can redistribute it and/or modify
6 it under the terms of the GNU General Public License as published by
7 the Free Software Foundation; either version 2 of the License, or
8 (at your option) any later version.
10 This program is distributed in the hope that it will be useful,
11 but WITHOUT ANY WARRANTY; without even the implied warranty of
12 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
13 GNU General Public License for more details.
15 You should have received a copy of the GNU General Public License
16 along with this program; if not, write to the Free Software
17 Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
23 #include <cairomm/context.h>
25 #include "canvas/utils.h"
28 using namespace ArdourCanvas;
31 ArdourCanvas::set_source_rgba (Cairo::RefPtr<Cairo::Context> context, Color color)
33 context->set_source_rgba (
34 ((color >> 24) & 0xff) / 255.0,
35 ((color >> 16) & 0xff) / 255.0,
36 ((color >> 8) & 0xff) / 255.0,
37 ((color >> 0) & 0xff) / 255.0
42 ArdourCanvas::set_source_rgb_a (Cairo::RefPtr<Cairo::Context> context, Color color, float alpha)
44 context->set_source_rgba (
45 ((color >> 24) & 0xff) / 255.0,
46 ((color >> 16) & 0xff) / 255.0,
47 ((color >> 8) & 0xff) / 255.0,
53 ArdourCanvas::set_source_rgba (cairo_t *cr, Color color)
55 cairo_set_source_rgba ( cr,
56 ((color >> 24) & 0xff) / 255.0,
57 ((color >> 16) & 0xff) / 255.0,
58 ((color >> 8) & 0xff) / 255.0,
59 ((color >> 0) & 0xff) / 255.0
64 ArdourCanvas::set_source_rgb_a (cairo_t *cr, Color color, float alpha)
66 cairo_set_source_rgba ( cr,
67 ((color >> 24) & 0xff) / 255.0,
68 ((color >> 16) & 0xff) / 255.0,
69 ((color >> 8) & 0xff) / 255.0,
74 ArdourCanvas::Distance
75 ArdourCanvas::distance_to_segment_squared (Duple const & p, Duple const & p1, Duple const & p2, double& t, Duple& at)
77 static const double kMinSegmentLenSquared = 0.00000001; // adjust to suit. If you use float, you'll probably want something like 0.000001f
78 static const double kEpsilon = 1.0E-14; // adjust to suit. If you use floats, you'll probably want something like 1E-7f
79 double dx = p2.x - p1.x;
80 double dy = p2.y - p1.y;
81 double dp1x = p.x - p1.x;
82 double dp1y = p.y - p1.y;
83 const double segLenSquared = (dx * dx) + (dy * dy);
85 if (segLenSquared >= -kMinSegmentLenSquared && segLenSquared <= kMinSegmentLenSquared) {
86 // segment is a point.
89 return ((dp1x * dp1x) + (dp1y * dp1y));
93 // Project a line from p to the segment [p1,p2]. By considering the line
94 // extending the segment, parameterized as p1 + (t * (p2 - p1)),
95 // we find projection of point p onto the line.
96 // It falls where t = [(p - p1) . (p2 - p1)] / |p2 - p1|^2
98 t = ((dp1x * dx) + (dp1y * dy)) / segLenSquared;
101 // intersects at or to the "left" of first segment vertex (p1.x, p1.y). If t is approximately 0.0, then
102 // intersection is at p1. If t is less than that, then there is no intersection (i.e. p is not within
103 // the 'bounds' of the segment)
105 // intersects at 1st segment vertex
108 // set our 'intersection' point to p1.
110 // Note: If you wanted the ACTUAL intersection point of where the projected lines would intersect if
111 // we were doing PointLineDistanceSquared, then qx would be (p1.x + (t * dx)) and qy would be (p1.y + (t * dy)).
113 } else if (t > (1.0 - kEpsilon)) {
114 // intersects at or to the "right" of second segment vertex (p2.x, p2.y). If t is approximately 1.0, then
115 // intersection is at p2. If t is greater than that, then there is no intersection (i.e. p is not within
116 // the 'bounds' of the segment)
117 if (t < (1.0 + kEpsilon)) {
118 // intersects at 2nd segment vertex
121 // set our 'intersection' point to p2.
123 // Note: If you wanted the ACTUAL intersection point of where the projected lines would intersect if
124 // we were doing PointLineDistanceSquared, then qx would be (p1.x + (t * dx)) and qy would be (p1.y + (t * dy)).
126 // The projection of the point to the point on the segment that is perpendicular succeeded and the point
127 // is 'within' the bounds of the segment. Set the intersection point as that projected point.
128 at = Duple (p1.x + (t * dx), p1.y + (t * dy));
131 // return the squared distance from p to the intersection point. Note that we return the squared distance
132 // as an optimization because many times you just need to compare relative distances and the squared values
133 // works fine for that. If you want the ACTUAL distance, just take the square root of this value.
134 double dpqx = p.x - at.x;
135 double dpqy = p.y - at.y;
137 return ((dpqx * dpqx) + (dpqy * dpqy));