4 * Copyright (C) 2008 Adam Williams <broadcast at earthling dot net>
6 * This program is free software; you can redistribute it and/or modify
7 * it under the terms of the GNU General Public License as published by
8 * the Free Software Foundation; either version 2 of the License, or
9 * (at your option) any later version.
11 * This program is distributed in the hope that it will be useful,
12 * but WITHOUT ANY WARRANTY; without even the implied warranty of
13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14 * GNU General Public License for more details.
16 * You should have received a copy of the GNU General Public License
17 * along with this program; if not, write to the Free Software
18 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
27 #define SQR(x) ((x) * (x))
29 #define CLIP(x, y, z) ((x) < (y) ? (y) : ((x) > (z) ? (z) : (x)))
31 #define RECLIP(x, y, z) ((x) = ((x) < (y) ? (y) : ((x) > (z) ? (z) : (x))))
33 #define CLAMP(x, y, z) ((x) = ((x) < (y) ? (y) : ((x) > (z) ? (z) : (x))))
35 #define MAX(x, y) ((x) > (y) ? (x) : (y))
37 #define MIN(x, y) ((x) < (y) ? (x) : (y))
39 #define EQUIV(x, y) (fabs((x) - (y)) < 0.001)
41 #define DISTANCE(x1, y1, x2, y2) \
42 (sqrt(((x2) - (x1)) * ((x2) - (x1)) + ((y2) - (y1)) * ((y2) - (y1))))
43 #define TO_RAD(x) ((x) * 2 * M_PI / 360)
44 #define TO_DEG(x) ((x) * 360 / 2 / M_PI)
46 static inline int bclip(int &iv, int imn, int imx) {
47 return iv < imn ? imn : iv > imx ? imx : iv;
49 static inline float bclip(float &fv, float fmn, float fmx) {
50 return fv < fmn ? fmn : fv > fmx ? fmx : fv;
52 static inline double bclip(double &dv, double dmn, double dmx) {
53 return dv < dmn ? dmn : dv > dmx ? dmx : dv;
55 static inline void bclamp(int &iv, int imn, int imx) {
56 if( iv < imn ) iv = imn; else if( iv > imx ) iv = imx;
58 static inline void bclamp(float &fv, float fmn, float fmx) {
59 if( fv < fmn ) fv = fmn; else if( fv > fmx ) fv = fmx;
61 static inline void bclamp(double &dv, double dmn, double dmx) {
62 if( dv < dmn ) dv = dmn; else if( dv > dmx ) dv = dmx;
64 static inline void bc_rgb2yuv(int r, int g, int b, int &y, int &u, int &v, int max=255)
65 { //bt601, jpeg, unclipped
66 double mx = max, rr = r/mx, gg = g/mx, bb = b/mx;
67 y = (int)(( 0.29900*rr + 0.58700*gg + 0.11400*bb) * mx + 0.5);
68 u = (int)((-0.16874*rr - 0.33126*gg + 0.50000*bb + 0.5) * mx + 0.5);
69 v = (int)(( 0.50000*rr - 0.41869*gg - 0.08131*bb + 0.5) * mx + 0.5);
71 static inline void bc_yuv2rgb(int y, int u, int v, int &r, int &g, int &b, int max=255)
73 int ofs = (max + 1) / 2;
74 double mx = max, yy = y/mx, uu = (u-ofs)/mx, vv = (v-ofs)/mx;
75 r = (int)((yy + 1.40200*vv) * mx + 0.5);
76 g = (int)((yy - 0.34414*uu - 0.71414*vv) * mx + 0.5);
77 b = (int)((yy + 1.77200*uu) * mx + 0.5);