libs/canvas/utils.cc

   1 /*
   2     Copyright (C) 2011-2013 Paul Davis
   3     Author: Carl Hetherington <cth@carlh.net>
   4
   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.
   9
  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.
  14
  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.
  18 */
  19
  20 #include <algorithm>
  21 #include <cmath>
  22 #include <stdint.h>
  23 #include <cairomm/context.h>
  24 #include "canvas/utils.h"
  25
  26 using std::max;
  27 using std::min;
  28
  29 void
  30 ArdourCanvas::color_to_hsv (Color color, double& h, double& s, double& v)
  31 {
  32         double r, g, b, a;
  33         double cmax;
  34         double cmin;
  35         double delta;
  36
  37         color_to_rgba (color, r, g, b, a);
  38
  39         if (r > g) {
  40                 cmax = max (r, b);
  41         } else {
  42                 cmax = max (g, b);
  43         }
  44
  45         if (r < g) {
  46                 cmin = min (r, b);
  47         } else {
  48                 cmin = min (g, b);
  49         }
  50
  51         v = cmax;
  52
  53         delta = cmax - cmin;
  54
  55         if (cmax == 0) {
  56                 // r = g = b == 0 ... v is undefined, s = 0
  57                 s = 0.0;
  58                 h = -1.0;
  59         }
  60
  61         if (delta != 0.0) {
  62                 if (cmax == r) {
  63                         h = fmod ((g - b)/delta, 6.0);
  64                 } else if (cmax == g) {
  65                         h = ((b - r)/delta) + 2;
  66                 } else {
  67                         h = ((r - g)/delta) + 4;
  68                 }
  69
  70                 h *= 60.0;
  71         }
  72
  73         if (delta == 0 || cmax == 0) {
  74                 s = 0;
  75         } else {
  76                 s = delta / cmax;
  77         }
  78
  79 }
  80
  81 ArdourCanvas::Color
  82 ArdourCanvas::hsv_to_color (double h, double s, double v, double a)
  83 {
  84         s = min (1.0, max (0.0, s));
  85         v = min (1.0, max (0.0, v));
  86
  87         if (s == 0) {
  88                 // achromatic (grey)
  89                 return rgba_to_color (v, v, v, a);
  90         }
  91
  92         h = min (360.0, max (0.0, h));
  93
  94         double c = v * s;
  95         double x = c * (1.0 - fabs(fmod(h / 60.0, 2) - 1.0));
  96         double m = v - c;
  97
  98         if (h >= 0.0 && h < 60.0) {
  99                 return rgba_to_color (c + m, x + m, m, a);
 100         } else if (h >= 60.0 && h < 120.0) {
 101                 return rgba_to_color (x + m, c + m, m, a);
 102         } else if (h >= 120.0 && h < 180.0) {
 103                 return rgba_to_color (m, c + m, x + m, a);
 104         } else if (h >= 180.0 && h < 240.0) {
 105                 return rgba_to_color (m, x + m, c + m, a);
 106         } else if (h >= 240.0 && h < 300.0) {
 107                 return rgba_to_color (x + m, m, c + m, a);
 108         } else if (h >= 300.0 && h < 360.0) {
 109                 return rgba_to_color (c + m, m, x + m, a);
 110         }
 111         return rgba_to_color (m, m, m, a);
 112 }
 113
 114 void
 115 ArdourCanvas::color_to_rgba (Color color, double& r, double& g, double& b, double& a)
 116 {
 117         r = ((color >> 24) & 0xff) / 255.0;
 118         g = ((color >> 16) & 0xff) / 255.0;
 119         b = ((color >>  8) & 0xff) / 255.0;
 120         a = ((color >>  0) & 0xff) / 255.0;
 121 }
 122
 123 ArdourCanvas::Color
 124 ArdourCanvas::rgba_to_color (double r, double g, double b, double a)
 125 {
 126         /* clamp to [0 .. 1] range */
 127
 128         r = min (1.0, max (0.0, r));
 129         g = min (1.0, max (0.0, g));
 130         b = min (1.0, max (0.0, b));
 131         a = min (1.0, max (0.0, a));
 132
 133         /* convert to [0..255] range */
 134
 135         unsigned int rc, gc, bc, ac;
 136         rc = rint (r * 255.0);
 137         gc = rint (g * 255.0);
 138         bc = rint (b * 255.0);
 139         ac = rint (a * 255.0);
 140
 141         /* build-an-integer */
 142
 143         return (rc << 24) | (gc << 16) | (bc << 8) | ac;
 144 }
 145
 146 void
 147 ArdourCanvas::set_source_rgba (Cairo::RefPtr<Cairo::Context> context, Color color)
 148 {
 149         context->set_source_rgba (
 150                 ((color >> 24) & 0xff) / 255.0,
 151                 ((color >> 16) & 0xff) / 255.0,
 152                 ((color >>  8) & 0xff) / 255.0,
 153                 ((color >>  0) & 0xff) / 255.0
 154                 );
 155 }
 156
 157 ArdourCanvas::Distance
 158 ArdourCanvas::distance_to_segment_squared (Duple const & p, Duple const & p1, Duple const & p2, double& t, Duple& at)
 159 {
 160         static const double kMinSegmentLenSquared = 0.00000001;  // adjust to suit.  If you use float, you'll probably want something like 0.000001f
 161         static const double kEpsilon = 1.0E-14;  // adjust to suit.  If you use floats, you'll probably want something like 1E-7f
 162         double dx = p2.x - p1.x;
 163         double dy = p2.y - p1.y;
 164         double dp1x = p.x - p1.x;
 165         double dp1y = p.y - p1.y;
 166         const double segLenSquared = (dx * dx) + (dy * dy);
 167
 168         if (segLenSquared >= -kMinSegmentLenSquared && segLenSquared <= kMinSegmentLenSquared) {
 169                 // segment is a point.
 170                 at = p1;
 171                 t = 0.0;
 172                 return ((dp1x * dp1x) + (dp1y * dp1y));
 173         }
 174
 175
 176         // Project a line from p to the segment [p1,p2].  By considering the line
 177         // extending the segment, parameterized as p1 + (t * (p2 - p1)),
 178         // we find projection of point p onto the line.
 179         // It falls where t = [(p - p1) . (p2 - p1)] / |p2 - p1|^2
 180
 181         t = ((dp1x * dx) + (dp1y * dy)) / segLenSquared;
 182
 183         if (t < kEpsilon) {
 184                 // intersects at or to the "left" of first segment vertex (p1.x, p1.y).  If t is approximately 0.0, then
 185                 // intersection is at p1.  If t is less than that, then there is no intersection (i.e. p is not within
 186                 // the 'bounds' of the segment)
 187                 if (t > -kEpsilon) {
 188                         // intersects at 1st segment vertex
 189                         t = 0.0;
 190                 }
 191                 // set our 'intersection' point to p1.
 192                 at = p1;
 193                 // Note: If you wanted the ACTUAL intersection point of where the projected lines would intersect if
 194                 // we were doing PointLineDistanceSquared, then qx would be (p1.x + (t * dx)) and qy would be (p1.y + (t * dy)).
 195
 196         } else if (t > (1.0 - kEpsilon)) {
 197                 // intersects at or to the "right" of second segment vertex (p2.x, p2.y).  If t is approximately 1.0, then
 198                 // intersection is at p2.  If t is greater than that, then there is no intersection (i.e. p is not within
 199                 // the 'bounds' of the segment)
 200                 if (t < (1.0 + kEpsilon)) {
 201                         // intersects at 2nd segment vertex
 202                         t = 1.0;
 203                 }
 204                 // set our 'intersection' point to p2.
 205                 at = p2;
 206                 // Note: If you wanted the ACTUAL intersection point of where the projected lines would intersect if
 207                 // we were doing PointLineDistanceSquared, then qx would be (p1.x + (t * dx)) and qy would be (p1.y + (t * dy)).
 208         } else {
 209                 // The projection of the point to the point on the segment that is perpendicular succeeded and the point
 210                 // is 'within' the bounds of the segment.  Set the intersection point as that projected point.
 211                 at = Duple (p1.x + (t * dx), p1.y + (t * dy));
 212         }
 213
 214         // return the squared distance from p to the intersection point.  Note that we return the squared distance
 215         // as an optimization because many times you just need to compare relative distances and the squared values
 216         // works fine for that.  If you want the ACTUAL distance, just take the square root of this value.
 217         double dpqx = p.x - at.x;
 218         double dpqy = p.y - at.y;
 219
 220         return ((dpqx * dpqx) + (dpqy * dpqy));
 221 }
 222