X-Git-Url: http://git.cinelerra-gg.org/git/?a=blobdiff_plain;f=cinelerra-5.1%2Fplugins%2Fbluebanana%2Fbluebananacolor.c;fp=cinelerra-5.1%2Fplugins%2Fbluebanana%2Fbluebananacolor.c;h=0de216de661580925b992bb469c7782188d7af13;hb=30bdb85eb33a8ee7ba675038a86c6be59c43d7bd;hp=0000000000000000000000000000000000000000;hpb=52fcc46226f9df46f9ce9d0566dc568455a7db0b;p=goodguy%2Fhistory.git

diff --git a/cinelerra-5.1/plugins/bluebanana/bluebananacolor.c b/cinelerra-5.1/plugins/bluebanana/bluebananacolor.c
new file mode 100644
index 00000000..0de216de
--- /dev/null
+++ b/cinelerra-5.1/plugins/bluebanana/bluebananacolor.c
@@ -0,0 +1,409 @@
+/*
+ * 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;
+  }
+}