== RED one (redcode) ==

This adds redcode (the file format of RED one, R3D) support to blender.
Seems to work fine with the footage I found on the web, but keep in
mind, that because of the unoptimized nature of libopenjpeg, frame
decoding isn't that fast.

It is also a rather challenging task, to make 4k-float-footage realtime :)

1 files changed, 129 insertions, 0 deletions

diff --git a/extern/libredcode/debayer.c b/extern/libredcode/debayer.c
new file mode 100644
index 00000000000..f7f22e1cc54
--- /dev/null
+++ b/extern/libredcode/debayer.c
@@ -0,0 +1,129 @@
+#include "debayer.h"
+
+/* pretty simple but astonishingly very effective "debayer" function 
+ */
+
+void redcode_ycbcr2rgb_fullscale(
+	int ** planes, int width, int height, float * out)
+{
+	int x,y;
+	int pix_max = 4096;
+	int mask = pix_max - 1;
+	float Kb = 0.0722;
+	float Kr = 0.2126;
+	float *o;
+
+	for (y = 0; y < height; y++) {
+		for (x = 0; x < width; x++) {
+			int i = x + y*width;
+			int i_p = (y > 0) ? i-width : i;
+			int i_n = (y < (height-1)) ? i + width : i;
+			float y1n = planes[0][i_n] & mask;
+			float y1  = planes[0][i] & mask;
+			float cb  = (planes[1][i] & mask)   - pix_max/2;
+			float cr  = (planes[2][i] & mask)   - pix_max/2;
+			float y2  = (planes[3][i] & mask);
+			float y2p = (planes[3][i_p] & mask);
+
+			float b_ = cb * (1.0 - Kb)/(pix_max/2);
+			float r_ = cr * (1.0 - Kr)/(pix_max/2);
+			float g_ = (- Kr * r_ - Kb * b_)/(1.0 - Kr - Kb);
+		
+			float y_[4] = {y1 / pix_max, 
+				       (y2 + y2p)/2 / pix_max, 
+				       (y1 + y1n)/2 / pix_max, 
+				       y2 / pix_max};
+
+			int j;
+			int yc = 0;
+
+			o = out + (2*height-1-2*y)*2*4*width 
+				+ x*2*4;
+
+			for (j = 0; j < 8; j += 4) {
+				o[j+0] = r_ + y_[yc];
+				o[j+1] = g_ + y_[yc];
+				o[j+2] = b_ + y_[yc];
+				o[j+3] = 1.0;
+				yc++;
+			}
+			
+			o = out + (2*height-1-2*y)*2*4*width 
+				+ x*2*4 - 2*4*width;
+
+			for (j = 0; j < 8; j += 4) {
+				o[j+0] = r_ + y_[yc];
+				o[j+1] = g_ + y_[yc];
+				o[j+2] = b_ + y_[yc];
+				o[j+3] = 1.0;
+				yc++;
+			}
+		}
+	}
+}
+
+void redcode_ycbcr2rgb_halfscale(
+	int ** planes, int width, int height, float * out)
+{
+	int x,y;
+	int pix_max = 4096;
+	int mask = pix_max - 1;
+	float Kb = 0.0722;
+	float Kr = 0.2126;
+
+	for (y = 0; y < height; y++) {
+		float *o = out + width * (height - y - 1);
+		for (x = 0; x < width; x++) {
+			int i = y*height + x;
+			float y1  = (planes[0][i] & mask);
+			float cb  = (planes[1][i] & mask)  - pix_max/2;
+			float cr  = (planes[2][i] & mask)  - pix_max/2;
+			float y2  = (planes[3][i] & mask);
+
+			float b_ = cb * (1.0 - Kb)/(pix_max/2);
+			float r_ = cr * (1.0 - Kr)/(pix_max/2);
+			float g_ = (- Kr * r_ - Kb * b_)/(1.0 - Kr - Kb);
+			
+			float y = (y1 + y2)/2 / pix_max;
+
+			*o++ = r_ + y;
+			*o++ = g_ + y;
+			*o++ = b_ + y;
+			*o++ = 1.0;
+		}
+	}
+}
+
+
+void redcode_ycbcr2rgb_quarterscale(
+	int ** planes, int width, int height, float * out)
+{
+	int x,y;
+	int pix_max = 4096;
+	int mask = pix_max - 1;
+	float Kb = 0.0722;
+	float Kr = 0.2126;
+
+	for (y = 0; y < height; y += 2) {
+		float *o = out + (width/2) * (height/2 - y/2 - 1);
+		for (x = 0; x < width; x += 2) {
+			int i = y * width + x;
+			float y1  = planes[0][i] & mask;
+			float cb  = (planes[1][i] & mask)  - pix_max/2;
+			float cr  = (planes[2][i] & mask)  - pix_max/2;
+			float y2  = planes[3][i] & mask;
+
+			float b_ = cb * (1.0 - Kb)/(pix_max/2);
+			float r_ = cr * (1.0 - Kr)/(pix_max/2);
+			float g_ = (- Kr * r_ - Kb * b_)/(1.0 - Kr - Kb);
+			
+			float y = (y1 + y2)/2 / pix_max;
+			
+			*o++ = r_ + y;
+			*o++ = g_ + y;
+			*o++ = b_ + y;
+			*o++ = 1.0;
+		}
+	}
+}
+