+++ /dev/null
-/*
- * Cinelerra :: Blue Banana - color modification plugin for Cinelerra-CV
- * Copyright (C) 2012-2013 Monty <monty@xiph.org>
- *
- * This program is free software; you can redistribute it and/or modify
- * it under the terms of the GNU General Public License as published by
- * the Free Software Foundation; either version 2 of the License, or
- * (at your option) any later version.
- *
- * This program is distributed in the hope that it will be useful,
- * but WITHOUT ANY WARRANTY; without even the implied warranty of
- * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
- * GNU General Public License for more details.
- *
- * You should have received a copy of the GNU General Public License
- * along with this program; if not, write to the Free Software
- * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
- *
- */
-
-static inline void rgb8_to_RGB(unsigned char *row, float *R, float *G, float *B, int w){
- while(w--){
- *R++ = *row++*.0039215686f;
- *G++ = *row++*.0039215686f;
- *B++ = *row++*.0039215686f;
- }
-}
-
-static inline void rgba8_to_RGBA(unsigned char *row, float *R, float *G, float *B, float *A, int w){
- while(w--){
- *R++ = *row++*.0039215686f;
- *G++ = *row++*.0039215686f;
- *B++ = *row++*.0039215686f;
- *A++ = *row++*.0039215686f;
- }
-}
-
-static inline void rgbF_to_RGB(float *row, float *R, float *G, float *B, int w){
- while(w--){
- *R++ = *row++;
- *G++ = *row++;
- *B++ = *row++;
- }
-}
-
-static inline void rgbaF_to_RGBA(float *row, float *R, float *G, float *B, float *A, int w){
- while(w--){
- *R++ = *row++;
- *G++ = *row++;
- *B++ = *row++;
- *A++ = *row++;
- }
-}
-
-
-// Full swing and bt601 are chosen below to conform to existing Cinelerra convention.
-// It's not correct, but it's consistent.
-
-#if 1
-// Full swing
-#define Y_SWING 255.f
-#define Y_SHIFT 0
-#define C_SWING 255.f
-#define C_SHIFT 128
-#define CLAMP_LO 0
-#define CLAMP_HI 255
-#else
-// Studio swing
-#define Y_SWING 219.f
-#define Y_SHIFT 16
-#define C_SWING 224.f
-#define C_SHIFT 128
-#define CLAMP_LO 1
-#define CLAMP_HI 254
-#endif
-
-#if 1
-// bt601
-#define Kr .299f
-#define Kb .114f
-#else
-// bt709
-#define Kr .2126f
-#define Kb .0722f
-#endif
-#define Kg (1.f-Kr-Kb)
-
-#define Br (.5f*-Kr/(1.f-Kb))
-#define Bg (.5f*-Kg/(1.f-Kb))
-#define Bb (.5f)
-
-#define Rr (.5f)
-#define Rg (.5f*-Kg/(1.f-Kr))
-#define Rb (.5f*-Kb/(1.f-Kr))
-
-#define Rv (2.f*(1.f-Kr))
-#define Gu (2.f*(1.f-Kb)*Kb/Kg)
-#define Gv (2.f*(1.f-Kr)*Kr/Kg)
-#define Bu (2.f*(1.f-Kb))
-
-static inline void yuv8_to_RGB(unsigned char *row, float *R, float *G, float *B, int w){
- while(w--){
- float y = (*row++-Y_SHIFT) * (1.f/Y_SWING);
- float u = (*row++-C_SHIFT);
- float v = (*row++-C_SHIFT);
- *R++ = y + (Rv/C_SWING)*v;
- *G++ = y - (Gu/C_SWING)*u - (Gv/C_SWING)*v;
- *B++ = y + (Bu/C_SWING)*u;
- }
-}
-
-static inline void yuva8_to_RGBA(unsigned char *row, float *R, float *G, float *B, float *A, int w){
- while(w--){
- float y = (*row++-Y_SHIFT) * (1.f/Y_SWING);
- float u = (*row++-C_SHIFT);
- float v = (*row++-C_SHIFT);
- *R++ = y + (Rv/C_SWING)*v;
- *G++ = y - (Gu/C_SWING)*u - (Gv/C_SWING)*v;
- *B++ = y + (Bu/C_SWING)*u;
- *A++ = *row++*(1.f/255.f);
- }
-}
-
-static inline void RGB_to_yuv8(float *R, float *G, float *B, float *S, float F, unsigned char *row, int w, int bpp){
- if(F>SELECT_THRESH){
- if(S){
- if(F<1.-SELECT_THRESH){
- while(w--){
- float s = *S++*F;
- float y = (Y_SHIFT + (Kr*Y_SWING)**R + (Kg*Y_SWING)**G + (Kb*Y_SWING)**B - row[0])*s + row[0] +.5f;
- float u = (C_SHIFT + (Br*C_SWING)**R + (Bg*C_SWING)**G + (Bb*C_SWING)**B - row[1])*s + row[1] +.5f;
- float v = (C_SHIFT + (Rr*C_SWING)**R++ + (Rg*C_SWING)**G++ + (Rb*C_SWING)**B++ - row[2])*s + row[2] +.5f;
- row[0] = CLAMP(y,CLAMP_LO,CLAMP_HI);
- row[1] = CLAMP(u,CLAMP_LO,CLAMP_HI);
- row[2] = CLAMP(v,CLAMP_LO,CLAMP_HI);
- row+=bpp;
- }
- }else{
- while(w--){
- float y = (Y_SHIFT + (Kr*Y_SWING)**R + (Kg*Y_SWING)**G + (Kb*Y_SWING)**B - row[0])**S + row[0] +.5f;
- float u = (C_SHIFT + (Br*C_SWING)**R + (Bg*C_SWING)**G + (Bb*C_SWING)**B - row[1])**S + row[1] +.5f;
- float v = (C_SHIFT + (Rr*C_SWING)**R++ + (Rg*C_SWING)**G++ + (Rb*C_SWING)**B++ - row[2])**S++ + row[2] +.5f;
- row[0] = CLAMP(y,CLAMP_LO,CLAMP_HI);
- row[1] = CLAMP(u,CLAMP_LO,CLAMP_HI);
- row[2] = CLAMP(v,CLAMP_LO,CLAMP_HI);
- row+=bpp;
- }
- }
- }else{
- if(F<1.-SELECT_THRESH){
- while(w--){
- float y = (Y_SHIFT + (Kr*Y_SWING)**R + (Kg*Y_SWING)**G + (Kb*Y_SWING)**B - row[0])*F + row[0] +.5f;
- float u = (C_SHIFT + (Br*C_SWING)**R + (Bg*C_SWING)**G + (Bb*C_SWING)**B - row[1])*F + row[1] +.5f;
- float v = (C_SHIFT + (Rr*C_SWING)**R++ + (Rg*C_SWING)**G++ + (Rb*C_SWING)**B++ - row[2])*F + row[2] +.5f;
- row[0] = CLAMP(y,CLAMP_LO,CLAMP_HI);
- row[1] = CLAMP(u,CLAMP_LO,CLAMP_HI);
- row[2] = CLAMP(v,CLAMP_LO,CLAMP_HI);
- row+=bpp;
- }
- }else{
- while(w--){
- float y = (Y_SHIFT+.5f) + (Kr*Y_SWING)**R + (Kg*Y_SWING)**G + (Kb*Y_SWING)**B ;
- float u = (C_SHIFT+.5f) + (Br*C_SWING)**R + (Bg*C_SWING)**G + (Bb*C_SWING)**B ;
- float v = (C_SHIFT+.5f) + (Rr*C_SWING)**R++ + (Rg*C_SWING)**G++ + (Rb*C_SWING)**B++;
- row[0] = CLAMP(y,CLAMP_LO,CLAMP_HI);
- row[1] = CLAMP(u,CLAMP_LO,CLAMP_HI);
- row[2] = CLAMP(v,CLAMP_LO,CLAMP_HI);
- row+=bpp;
- }
- }
- }
- }
-}
-
-static inline void RGB_to_rgb8(float *R, float *G, float *B, float *S, float F, unsigned char *row, int w, int bpp){
- if(F>SELECT_THRESH){
- if(S){
- if(F<1.-SELECT_THRESH){
- while(w--){
- float s = *S++*F;
- float r = (*R++*255.f-row[0])* s +row[0] +.5f;
- float g = (*G++*255.f-row[1])* s +row[1] +.5f;
- float b = (*B++*255.f-row[2])* s +row[2] +.5f;
- row[0] = CLAMP(r,0,255);
- row[1] = CLAMP(g,0,255);
- row[2] = CLAMP(b,0,255);
- row+=bpp;
- }
-
-
- }else{
- while(w--){
- float r = (*R++*255.f-row[0])* *S +row[0] +.5f;
- float g = (*G++*255.f-row[1])* *S +row[1] +.5f;
- float b = (*B++*255.f-row[2])* *S++ +row[2] +.5f;
- row[0] = CLAMP(r,0,255);
- row[1] = CLAMP(g,0,255);
- row[2] = CLAMP(b,0,255);
- row+=bpp;
- }
- }
- }else{
- if(F<1.-SELECT_THRESH){
- while(w--){
- float r = (*R++*255.f-row[0])* F +row[0] +.5f;
- float g = (*G++*255.f-row[1])* F +row[1] +.5f;
- float b = (*B++*255.f-row[2])* F +row[2] +.5f;
- row[0] = CLAMP(r,0,255);
- row[1] = CLAMP(g,0,255);
- row[2] = CLAMP(b,0,255);
- row+=bpp;
- }
- }else{
- while(w--){
- float r = *R++*255.f +.5f;
- float g = *G++*255.f +.5f;
- float b = *B++*255.f +.5f;
- row[0] = CLAMP(r,0,255);
- row[1] = CLAMP(g,0,255);
- row[2] = CLAMP(b,0,255);
- row+=bpp;
- }
- }
- }
- }
-}
-
-static inline void RGB_to_rgbF(float *R, float *G, float *B, float *S, float F, float *row, int w, int bpp){
- if(F>SELECT_THRESH){
- if(S){
- if(F<1.-SELECT_THRESH){
- while(w--){
- float s = *S++*F;
- row[0] = (*R++-row[0])* s +row[0];
- row[1] = (*G++-row[1])* s +row[1];
- row[2] = (*B++-row[2])* s +row[2];
- row+=bpp;
- }
- }else{
- while(w--){
- row[0] = (*R++-row[0])* *S +row[0];
- row[1] = (*G++-row[1])* *S +row[1];
- row[2] = (*B++-row[2])* *S++ +row[2];
- row+=bpp;
- }
- }
- }else{
- if(F<1.-SELECT_THRESH){
- while(w--){
- row[0] = (*R++-row[0])* F +row[0];
- row[1] = (*G++-row[1])* F +row[1];
- row[2] = (*B++-row[2])* F +row[2];
- row+=bpp;
- }
- }else{
- while(w--){
- row[0] = *R++;
- row[1] = *G++;
- row[2] = *B++;
- row+=bpp;
- }
- }
- }
- }
-}
-
-static inline void unmask_rgba8(unsigned char *row,int w){
- while(w--){
- row[3] = 255;
- row+=4;
- }
-}
-
-static inline void unmask_rgbaF(float *row, int w){
- while(w--){
- row[3] = 1.f;
- row+=4;
- }
-}
-
-static inline void unmask_yuva8(unsigned char *row, int w){
- while(w--){
- row[3] = 255;
- row+=4;
- }
-}
-
-// A modified HSV colorspace: saturation is calculated using |V| + a
-// fixed bias to shift the absolute blackpoint well below the in-use
-// color range. This eliminates the singularity and desaturates
-// values near the functional blackpoint (zero). The modifications
-// also allows us to use HSV with footroom without clipping chroma
-// when V<=0.
-
-// H is constrained to the range 0 to 1.
-// S is strictly >=0 but the upper range is technically unbounded.
-// It is reasonable to clip to the reference range [0:1]
-// V is unbounded and may extend considerably below 0 and above 1 depending on foot and head room
-
-/* Simple mods of the original Cinelerra conversion routines */
-static inline void RGB_to_HSpV(float R, float G, float B, float &H, float &Sp, float &V){
- if(R<G){
- if(B<R){
- V = G;
- H = 2 + (B-R) / (G-B);
- Sp = (G-B)/(fabs(V)+HSpV_SATURATION_BIAS)*HSpV_SATURATION_SCALE;
- }else{
- if(B<G){
- V = G;
- H = 2 + (B-R) / (G-R);
- Sp = (G-R)/(fabs(V)+HSpV_SATURATION_BIAS)*HSpV_SATURATION_SCALE;
- }else{
- V = B;
- H = 4 + (R-G) / (B-R);
- Sp = (B-R)/(fabs(V)+HSpV_SATURATION_BIAS)*HSpV_SATURATION_SCALE;
- }
- }
- }else{
- if(B<G){
- V = R;
- H = (G-B) / (R-B);
- Sp = (R-B)/(fabs(V)+HSpV_SATURATION_BIAS)*HSpV_SATURATION_SCALE;
- }else{
- if(B<R){
- V = R;
- H = 6 + (G-B) / (R-G);
- Sp = (R-G)/(fabs(V)+HSpV_SATURATION_BIAS)*HSpV_SATURATION_SCALE;
- }else{
- V = B;
- H = 4 + (R-G) / (B-G+.0001);
- Sp = (B-G)/(fabs(V)+HSpV_SATURATION_BIAS)*HSpV_SATURATION_SCALE;
- }
- }
- }
-}
-
-static inline void HSpV_to_RGB(float H, float Sp, float V, float &R, float &G, float &B){
- float vp = (fabs(V)+HSpV_SATURATION_BIAS)*Sp*HSpV_SATURATION_ISCALE;
- int i = (int)H;
- switch(i){
- default:
- R = V;
- G = V + (H-i-1)*vp;
- B = V - vp;
- return;
- case 1:
- R = V - (H-i)*vp;
- G = V;
- B = V - vp;
- return;
- case 2:
- R = V - vp;
- G = V;
- B = V + (H-i-1)*vp;
- return;
- case 3:
- R = V - vp;
- G = V - (H-i)*vp;
- B = V;
- return;
- case 4:
- R = V + (H-i-1)*vp;
- G = V - vp;
- B = V;
- return;
- case 5:
- R = V;
- G = V - vp;
- B = V - (H-i)*vp;
- return;
- }
-}
-
-static inline void HSpV_correct_RGB(float H, float Sp, float &V, float &R, float &G, float &B){
- float vp = Sp<0.f?0.f:(fabs(V)+HSpV_SATURATION_BIAS)*Sp*HSpV_SATURATION_ISCALE;
- int i = (int)H;
- switch(i){
- default:
- R = V = R*R_TO_Y + (G-(H-i-1)*vp)*G_TO_Y + (B+vp)*B_TO_Y;
- G = V + (H-i-1)*vp;
- B = V - vp;
- return;
- case 1:
- G = V = (R+(H-i)*vp)*R_TO_Y + G*G_TO_Y + (B+vp)*B_TO_Y;
- R = V - (H-i)*vp;
- B = V - vp;
- return;
- case 2:
- G = V = (R+vp)*R_TO_Y + G*G_TO_Y + (B-(H-i-1)*vp)*B_TO_Y;
- R = V - vp;
- B = V + (H-i-1)*vp;
- return;
- case 3:
- B = V = (R+vp)*R_TO_Y + (G+(H-i)*vp)*G_TO_Y + B*B_TO_Y;
- R = V - vp;
- G = V - (H-i)*vp;
- return;
- case 4:
- B = V = (R-(H-i-1)*vp)*R_TO_Y + (G+vp)*G_TO_Y + B*B_TO_Y;
- R = V + (H-i-1)*vp;
- G = V - vp;
- return;
- case 5:
- R = V = R*R_TO_Y + (G+vp)*G_TO_Y + (B+(H-i)*vp)*B_TO_Y;
- G = V - vp;
- B = V - (H-i)*vp;
- return;
- }
-}