Change compositor color balance to match the sequencer exactly, for this to work...
[blender.git] / source / blender / nodes / intern / CMP_nodes / CMP_colorbalance.c
1 /**
2  * 
3  *
4  * ***** BEGIN GPL LICENSE BLOCK *****
5  *
6  * This program is free software; you can redistribute it and/or
7  * modify it under the terms of the GNU General Public License
8  * as published by the Free Software Foundation; either version 2
9  * of the License, or (at your option) any later version. 
10  *
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.
15  *
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 Foundation,
18  * Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
19  *
20  * The Original Code is Copyright (C) 2006 Blender Foundation.
21  * All rights reserved.
22  *
23  * The Original Code is: all of this file.
24  *
25  * Contributor(s): Matt Ebb.
26  *
27  * ***** END GPL LICENSE BLOCK *****
28  */
29
30
31 #include "../CMP_util.h"
32
33
34 /* ******************* Color Balance ********************************* */
35 static bNodeSocketType cmp_node_colorbalance_in[]={
36         {SOCK_VALUE, 1, "Fac",  1.0f, 0.0f, 0.0f, 1.0f, 0.0f, 1.0f},
37         {SOCK_RGBA,1,"Image", 0.8f, 0.8f, 0.8f, 1.0f, 0.0f, 1.0f},
38         {-1,0,""}
39 };
40
41 static bNodeSocketType cmp_node_colorbalance_out[]={
42         {SOCK_RGBA,0,"Image", 0.0f, 0.0f, 1.0f, 1.0f, -1.0f, 1.0f},
43         {-1,0,""}
44 };
45
46 DO_INLINE float colorbalance_cdl(float in, float offset, float power, float slope)
47 {
48         float x = in * slope + offset;
49         
50         /* prevent NaN */
51         if (x < 0.f) x = 0.f;
52         
53         //powf(in * slope + offset, power)
54         return powf(x, 1.f/power);
55 }
56
57 /* note: lift_lgg is just 2-lift */
58 DO_INLINE float colorbalance_lgg(float in, float lift_lgg, float gamma, float gain)
59 {
60         /* 1:1 match with the sequencer with linear/srgb conversions, the conversion isnt pretty
61          * but best keep it this way, sice testing for durian shows a similar calculation
62          * without lin/srgb conversions gives bad results (over-saturated shadows) with colors
63          * slightly below 1.0. some correction can be done but it ends up looking bad for shadows or lighter tones - campbell */
64         float x= (((linearrgb_to_srgb(in) - 1.0f) * lift_lgg) + 1.0f) * gain;
65
66         /* prevent NaN */
67         if (x < 0.f) x = 0.f;
68
69         return powf(srgb_to_linearrgb(x), gamma);
70 }
71
72 static void do_colorbalance_cdl(bNode *node, float* out, float *in)
73 {
74         NodeColorBalance *n= (NodeColorBalance *)node->storage;
75         
76         out[0] = colorbalance_cdl(in[0], n->lift[0], n->gamma[0], n->gain[0]);
77         out[1] = colorbalance_cdl(in[1], n->lift[1], n->gamma[1], n->gain[1]);
78         out[2] = colorbalance_cdl(in[2], n->lift[2], n->gamma[2], n->gain[2]);
79         out[3] = in[3];
80 }
81
82 static void do_colorbalance_cdl_fac(bNode *node, float* out, float *in, float *fac)
83 {
84         NodeColorBalance *n= (NodeColorBalance *)node->storage;
85         const float mfac= 1.0f - *fac;
86         
87         out[0] = mfac*in[0] + *fac * colorbalance_cdl(in[0], n->lift[0], n->gamma[0], n->gain[0]);
88         out[1] = mfac*in[1] + *fac * colorbalance_cdl(in[1], n->lift[1], n->gamma[1], n->gain[1]);
89         out[2] = mfac*in[2] + *fac * colorbalance_cdl(in[2], n->lift[2], n->gamma[2], n->gain[2]);
90         out[3] = in[3];
91 }
92
93 static void do_colorbalance_lgg(bNode *node, float* out, float *in)
94 {
95         NodeColorBalance *n= (NodeColorBalance *)node->storage;
96
97         out[0] = colorbalance_lgg(in[0], n->lift_lgg[0], n->gamma[0], n->gain[0]);
98         out[1] = colorbalance_lgg(in[1], n->lift_lgg[1], n->gamma[1], n->gain[1]);
99         out[2] = colorbalance_lgg(in[2], n->lift_lgg[2], n->gamma[2], n->gain[2]);
100         out[3] = in[3];
101 }
102
103 static void do_colorbalance_lgg_fac(bNode *node, float* out, float *in, float *fac)
104 {
105         NodeColorBalance *n= (NodeColorBalance *)node->storage;
106         const float mfac= 1.0f - *fac;
107
108         out[0] = mfac*in[0] + *fac * colorbalance_lgg(in[0], n->lift_lgg[0], n->gamma[0], n->gain[0]);
109         out[1] = mfac*in[1] + *fac * colorbalance_lgg(in[1], n->lift_lgg[1], n->gamma[1], n->gain[1]);
110         out[2] = mfac*in[2] + *fac * colorbalance_lgg(in[2], n->lift_lgg[2], n->gamma[2], n->gain[2]);
111         out[3] = in[3];
112 }
113
114 static void node_composit_exec_colorbalance(void *data, bNode *node, bNodeStack **in, bNodeStack **out)
115 {
116         CompBuf *cbuf= in[1]->data;
117         CompBuf *stackbuf;
118         
119         /* stack order input:  fac, image */
120         /* stack order output: image */
121         if(out[0]->hasoutput==0) return;
122         
123         if(in[0]->vec[0] == 0.f && in[0]->data == NULL) {
124                 out[0]->data = pass_on_compbuf(cbuf);
125                 return;
126         }
127
128         {
129                 NodeColorBalance *n= (NodeColorBalance *)node->storage;
130                 int c;
131
132                 for (c = 0; c < 3; c++) {
133                         n->lift_lgg[c] = 2.0f - n->lift[c];
134                 }
135         }
136
137         if (cbuf) {
138                 stackbuf= alloc_compbuf(cbuf->x, cbuf->y, CB_RGBA, 1); /* create output based on image input */
139                         
140                 if (node->custom1 == 0) {
141                         /* lift gamma gain */
142                         if ((in[0]->data==NULL) && (in[0]->vec[0] >= 1.f)) {
143                                 composit1_pixel_processor(node, stackbuf, in[1]->data, in[1]->vec, do_colorbalance_lgg, CB_RGBA);
144                         }
145                         else {
146                                 composit2_pixel_processor(node, stackbuf, in[1]->data, in[1]->vec, in[0]->data, in[0]->vec, do_colorbalance_lgg_fac, CB_RGBA, CB_VAL);
147                         }
148                 } else {
149                         /* offset/power/slope : ASC-CDL */
150                         if ((in[0]->data==NULL) && (in[0]->vec[0] >= 1.f)) {
151                                 composit1_pixel_processor(node, stackbuf, in[1]->data, in[1]->vec, do_colorbalance_cdl, CB_RGBA);
152                         }
153                         else {
154                                 composit2_pixel_processor(node, stackbuf, in[1]->data, in[1]->vec, in[0]->data, in[0]->vec, do_colorbalance_cdl_fac, CB_RGBA, CB_VAL);
155                         }
156                         
157                 }
158
159                 out[0]->data=stackbuf;
160         }
161 }
162
163 static void node_composit_init_colorbalance(bNode *node)
164 {
165         NodeColorBalance *n= node->storage= MEM_callocN(sizeof(NodeColorBalance), "node colorbalance");
166
167         n->lift[0] = n->lift[1] = n->lift[2] = 1.0f;
168         n->gamma[0] = n->gamma[1] = n->gamma[2] = 1.0f;
169         n->gain[0] = n->gain[1] = n->gain[2] = 1.0f;
170 }
171
172 bNodeType cmp_node_colorbalance={
173         /* *next,*prev */       NULL, NULL,
174         /* type code   */       CMP_NODE_COLORBALANCE,
175         /* name        */       "Color Balance",
176         /* width+range */       400, 200, 400,
177         /* class+opts  */       NODE_CLASS_OP_COLOR, NODE_OPTIONS,
178         /* input sock  */       cmp_node_colorbalance_in,
179         /* output sock */       cmp_node_colorbalance_out,
180         /* storage     */       "NodeColorBalance",
181         /* execfunc    */       node_composit_exec_colorbalance,
182         /* butfunc     */       NULL,
183         /* initfunc    */       node_composit_init_colorbalance,
184         /* freestoragefunc    */        node_free_standard_storage,
185         /* copystoragefunc    */        node_copy_standard_storage,
186         /* id          */       NULL
187 };
188