From 24fd9b1360d877a8db8b3ef472c7651ec2cf71a2 Mon Sep 17 00:00:00 2001
From: lilo_booter <lilo_booter@d19143bc-622f-0410-bfdd-b5b2a6649095>
Date: Fri, 24 Dec 2004 09:54:24 +0000
Subject: [PATCH] Luma generation and use

git-svn-id: https://mlt.svn.sourceforge.net/svnroot/mlt/trunk/mlt@571 d19143bc-622f-0410-bfdd-b5b2a6649095
---
 src/modules/core/transition_composite.c |   8 +
 src/modules/core/transition_luma.c      |  10 +-
 src/modules/lumas/Makefile              |  22 ++
 src/modules/lumas/configure             |   0
 src/modules/lumas/create_lumas          |  30 ++
 src/modules/lumas/luma.c                | 390 ++++++++++++++++++++++++
 6 files changed, 459 insertions(+), 1 deletion(-)
 create mode 100644 src/modules/lumas/Makefile
 create mode 100755 src/modules/lumas/configure
 create mode 100755 src/modules/lumas/create_lumas
 create mode 100644 src/modules/lumas/luma.c

diff --git a/src/modules/core/transition_composite.c b/src/modules/core/transition_composite.c
index 5bc6aff5..6cc1055b 100644
--- a/src/modules/core/transition_composite.c
+++ b/src/modules/core/transition_composite.c
@@ -554,6 +554,14 @@ static uint16_t* get_luma( mlt_properties properties, int width, int height )
 	// If the filename property changed, reload the map
 	char *resource = mlt_properties_get( properties, "luma" );
 
+	char temp[ 512 ];
+
+	if ( strchr( resource, '%' ) )
+	{
+		sprintf( temp, "%s/lumas/%s/%s", mlt_factory_prefix( ), mlt_environment( "MLT_NORMALISATION" ), strchr( resource, '%' ) + 1 );
+		resource = temp;
+	}
+
 	if ( resource != NULL && ( luma_bitmap == NULL || luma_width != width || luma_height != height ) )
 	{
 		uint16_t *orig_bitmap = mlt_properties_get_data( properties, "_luma.orig_bitmap", NULL );
diff --git a/src/modules/core/transition_luma.c b/src/modules/core/transition_luma.c
index c241b3ec..248d8485 100644
--- a/src/modules/core/transition_luma.c
+++ b/src/modules/core/transition_luma.c
@@ -372,7 +372,15 @@ static int transition_get_image( mlt_frame a_frame, uint8_t **image, mlt_image_f
 
 	if ( luma_bitmap == NULL && resource != NULL )
 	{
-		char *extension = extension = strrchr( resource, '.' );
+		char temp[ 512 ];
+		char *extension = strrchr( resource, '.' );
+
+		if ( strchr( resource, '%' ) )
+		{
+			sprintf( temp, "%s/lumas/%s/%s", mlt_factory_prefix( ), mlt_environment( "MLT_NORMALISATION" ), strchr( resource, '%' ) + 1 );
+			resource = temp;
+			extension = strrchr( resource, '.' );
+		}
 
 		// See if it is a PGM
 		if ( extension != NULL && strcmp( extension, ".pgm" ) == 0 )
diff --git a/src/modules/lumas/Makefile b/src/modules/lumas/Makefile
new file mode 100644
index 00000000..94a9fa7b
--- /dev/null
+++ b/src/modules/lumas/Makefile
@@ -0,0 +1,22 @@
+include ../../../config.mak
+
+all:	luma create_lumas
+	@./create_lumas 
+
+luma:	luma.c
+
+create_lumas:
+
+depend:
+
+dist-clean:
+	rm -rf PAL NTSC luma 
+
+clean:
+	rm -f luma 
+
+install:	all
+	install -d $(prefix)/share/mlt/modules/lumas/PAL
+	install -d $(prefix)/share/mlt/modules/lumas/NTSC
+	install -m 644 PAL/* $(prefix)/share/mlt/modules/lumas/PAL
+	install -m 644 NTSC/* $(prefix)/share/mlt/modules/lumas/NTSC
diff --git a/src/modules/lumas/configure b/src/modules/lumas/configure
new file mode 100755
index 00000000..e69de29b
diff --git a/src/modules/lumas/create_lumas b/src/modules/lumas/create_lumas
new file mode 100755
index 00000000..dcc603c6
--- /dev/null
+++ b/src/modules/lumas/create_lumas
@@ -0,0 +1,30 @@
+#!/bin/sh
+
+[ \( -d PAL \) -a \( ! $0 -nt PAL/luma01.pgm \) ] && exit 0
+
+for i in PAL NTSC
+do
+	mkdir -p $i
+	[ "$i" == "PAL" ] && h=576 || h=480
+	./luma -h $h -bpp 16 > $i/luma01.pgm
+	./luma -h $h -bpp 16 -offset 32768 -dmirror 1 > $i/luma02.pgm
+	./luma -h $h -bpp 16 -bands 8 -rband 1 > $i/luma03.pgm
+	./luma -h $h -bpp 16 -hmirror 1 > $i/luma04.pgm
+	./luma -h $h -bpp 16 -bands 576 -vmirror 1 > $i/luma05.pgm
+	./luma -h $h -bpp 16 -bands 576 -vmirror 0 > $i/luma06.pgm
+	./luma -h $h -bpp 16 -bands 8 -rband 0 > $i/luma07.pgm
+	./luma -h $h -bpp 16 -offset 32768 -dmirror 1 -hmirror 1 > $i/luma08.pgm
+	./luma -h $h -bpp 16 -bands 8 -rband 0 > $i/luma09.pgm
+	./luma -h $h -bpp 16 -offset 32768 -dmirror 1 -flip 1 > $i/luma10.pgm
+	./luma -h $h -bpp 16 -offset 32768 -dmirror 1 -quart 1 > $i/luma11.pgm
+	./luma -h $h -bpp 16 -offset 32768 -dmirror 1 -quart 1 -flip 1 > $i/luma12.pgm
+	./luma -h $h -bpp 16 -offset 32768 -dmirror 1 -quart 1 -pflop 1 > $i/luma13.pgm
+	./luma -h $h -bpp 16 -offset 32768 -dmirror 1 -quart 1 -flip 1 -pflop 1 > $i/luma14.pgm
+	./luma -h $h -bpp 16 -type 1 > $i/luma15.pgm
+	./luma -h $h -bpp 16 -type 1 -bands 2 -rband 1 > $i/luma16.pgm
+	./luma -h $h -bpp 16 -type 2 > $i/luma17.pgm
+	./luma -h $h -bpp 16 -type 2 -quart 1 > $i/luma18.pgm
+	./luma -h $h -bpp 16 -type 2 -quart 1 -flip 1 > $i/luma19.pgm
+	./luma -h $h -bpp 16 -type 2 -quart 1 -bands 2 > $i/luma20.pgm
+	./luma -h $h -bpp 16 -type 3 > $i/luma21.pgm
+done
diff --git a/src/modules/lumas/luma.c b/src/modules/lumas/luma.c
new file mode 100644
index 00000000..76f927a1
--- /dev/null
+++ b/src/modules/lumas/luma.c
@@ -0,0 +1,390 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <stdint.h>
+#include <string.h>
+
+typedef struct 
+{
+	int type;
+	int w;
+	int h;
+	int bands;
+	int rband;
+	int vmirror;
+	int hmirror;
+	int dmirror;
+	int invert;
+	int offset;
+	int flip;
+	int flop;
+	int pflip;
+	int pflop;
+	int quart;
+}
+luma;
+
+void luma_init( luma *this )
+{
+	this->type = 0;
+	this->w = 720;
+	this->h = 576;
+	this->bands = 1;
+	this->rband = 0;
+	this->vmirror = 0;
+	this->hmirror = 0;
+	this->dmirror = 0;
+	this->invert = 0;
+	this->offset = 0;
+	this->flip = 0;
+	this->flop = 0;
+	this->quart = 0;
+	this->pflop = 0;
+	this->pflip = 0;
+}
+
+static inline int sqrti( int n )
+{
+    int p = 0;
+	int q = 1;
+	int r = n;
+	int h = 0;
+
+    while( q <= n ) 
+		q = 4 * q;
+
+    while( q != 1 )
+    {
+        q = q / 4;
+        h = p + q;
+        p = p / 2;
+        if ( r >= h )
+        {
+            p = p + q;
+            r = r - h;
+        } 
+    }
+
+    return p;
+}
+
+uint16_t *luma_render( luma *this )
+{
+	int i = 0;
+	int j = 0;
+	int k = 0;
+
+	if ( this->quart )
+	{
+		this->w *= 2;
+		this->h *= 2;
+	}
+
+	int max = ( 1 << 16 ) - 1;
+	uint16_t *image = malloc( this->w * this->h * sizeof( uint16_t ) );
+	uint16_t *end = image + this->w * this->h;
+	uint16_t *p = image;
+	uint16_t *r = image;
+	int lower = 0;
+	int lpb = this->h / this->bands;
+	int rpb = max / this->bands;
+	int direction = 1;
+
+	int half_w = this->w / 2;
+	int half_h = this->h / 2;
+
+	if ( !this->dmirror && ( this->hmirror || this->vmirror ) )
+		rpb *= 2;
+
+	for ( i = 0; i < this->bands; i ++ )
+	{
+		lower = i * rpb;
+		direction = 1;
+
+		if ( this->rband && i % 2 == 1 )
+		{
+			direction = -1;
+			lower += rpb;
+		}
+
+		switch( this->type )
+		{
+			case 1:
+				{
+					int length = sqrti( half_w * half_w + lpb * lpb / 4 );
+					int value;
+					int x = 0;
+					int y = 0;
+					for ( j = 0; j < lpb; j ++ )
+					{
+						y = j - lpb / 2;
+						for ( k = 0; k < this->w; k ++ )
+						{
+							x = k - half_w;
+							value = sqrti( x * x + y * y );
+							*p ++ = lower + ( direction * rpb * ( ( max * value ) / length ) / max ) + ( j * this->offset * 2 / lpb ) + ( j * this->offset / lpb );
+						}
+					}
+				}
+				break;
+
+			case 2:
+				{
+					for ( j = 0; j < lpb; j ++ )
+					{
+						int value = ( ( j * this->w ) / lpb ) - half_w;
+						if ( value > 0 )
+							value = - value;
+						for ( k = - half_w; k < value; k ++ )
+							*p ++ = lower + ( direction * rpb * ( ( max * abs( k ) ) / half_w ) / max );
+						for ( k = value; k < abs( value ); k ++ )
+							*p ++ = lower + ( direction * rpb * ( ( max * abs( value ) ) / half_w ) / max ) + ( j * this->offset * 2 / lpb ) + ( j * this->offset / lpb );
+						for ( k = abs( value ); k < half_w; k ++ )
+							*p ++ = lower + ( direction * rpb * ( ( max * abs( k ) ) / half_w ) / max );
+					}
+				}
+				break;
+
+			case 3:
+				{
+					int length;
+					for ( j = -half_h; j < half_h; j ++ )
+					{
+						if ( j < 0 )
+						{
+							for ( k = - half_w; k < half_w; k ++ )
+							{
+								length = sqrti( k * k + j * j );
+								*p ++ = ( max / 4 * k ) / ( length + 1 );
+							}
+						}
+						else
+						{
+							for ( k = half_w; k > - half_w; k -- )
+							{
+								length = sqrti( k * k + j * j );
+								*p ++ = ( max / 2 ) + ( max / 4 * k ) / ( length + 1 );
+							}
+						}
+					}
+				}
+				break;
+
+			default:
+				for ( j = 0; j < lpb; j ++ )
+					for ( k = 0; k < this->w; k ++ )
+						*p ++ = lower + ( direction * ( rpb * ( ( k * max ) / this->w ) / max ) ) + ( j * this->offset * 2 / lpb );
+				break;
+		}
+	}
+
+	if ( this->quart )
+	{
+		this->w /= 2;
+		this->h /= 2;
+		for ( i = 1; i < this->h; i ++ )
+		{
+			p = image + i * this->w;
+			r = image + i * 2 * this->w;
+			j = this->w;
+			while ( j -- > 0 )
+				*p ++ = *r ++;
+		}
+	}
+
+	if ( this->dmirror )
+	{
+		for ( i = 0; i < this->h; i ++ )
+		{
+			p = image + i * this->w;
+			r = end - i * this->w;
+			j = ( this->w * ( this->h - i ) ) / this->h;
+			while ( j -- )
+				*( -- r ) = *p ++;
+		}
+	}
+
+	if ( this->flip )
+	{
+		uint16_t t;
+		for ( i = 0; i < this->h; i ++ )
+		{
+			p = image + i * this->w;
+			r = p + this->w;
+			while( p != r )
+			{
+				t = *p;
+				*p ++ = *( -- r );
+				*r = t;
+			}
+		}
+	}
+
+	if ( this->flop )
+	{
+		uint16_t t;
+		r = end;
+		for ( i = 1; i < this->h / 2; i ++ )
+		{
+			p = image + i * this->w;
+			j = this->w;
+			while( j -- )
+			{
+				t = *( -- p );
+				*p = *( -- r );
+				*r = t;
+			}
+		}
+	}
+
+	if ( this->hmirror )
+	{
+		p = image;
+		while ( p < end )
+		{
+			r = p + this->w;
+			while ( p != r )
+				*( -- r ) = *p ++;
+			p += this->w / 2;
+		}
+	}
+
+	if ( this->vmirror )
+	{
+		p = image;
+		r = end;
+		while ( p != r )
+			*( -- r ) = *p ++;
+	}
+
+	if ( this->invert )
+	{
+		p = image;
+		r = image;
+		while ( p < end )
+			*p ++ = max - *r ++;
+	}
+
+	if ( this->pflip )
+	{
+		uint16_t t;
+		for ( i = 0; i < this->h; i ++ )
+		{
+			p = image + i * this->w;
+			r = p + this->w;
+			while( p != r )
+			{
+				t = *p;
+				*p ++ = *( -- r );
+				*r = t;
+			}
+		}
+	}
+
+	if ( this->pflop )
+	{
+		uint16_t t;
+		end = image + this->w * this->h;
+		r = end;
+		for ( i = 1; i < this->h / 2; i ++ )
+		{
+			p = image + i * this->w;
+			j = this->w;
+			while( j -- )
+			{
+				t = *( -- p );
+				*p = *( -- r );
+				*r = t;
+			}
+		}
+	}
+
+	return image;
+}
+
+int main( int argc, char **argv )
+{
+	int arg = 1;
+	int bpp = 8;
+
+	luma this;
+	uint16_t *image = NULL;
+
+	luma_init( &this );
+
+	for ( arg = 1; arg < argc - 1; arg ++ )
+	{
+		if ( !strcmp( argv[ arg ], "-bpp" ) )
+			bpp = atoi( argv[ ++ arg ] );
+		else if ( !strcmp( argv[ arg ], "-type" ) )
+			this.type = atoi( argv[ ++ arg ] );
+		else if ( !strcmp( argv[ arg ], "-w" ) )
+			this.w = atoi( argv[ ++ arg ] );
+		else if ( !strcmp( argv[ arg ], "-h" ) )
+			this.h = atoi( argv[ ++ arg ] );
+		else if ( !strcmp( argv[ arg ], "-bands" ) )
+			this.bands = atoi( argv[ ++ arg ] );
+		else if ( !strcmp( argv[ arg ], "-rband" ) )
+			this.rband = atoi( argv[ ++ arg ] );
+		else if ( !strcmp( argv[ arg ], "-hmirror" ) )
+			this.hmirror = atoi( argv[ ++ arg ] );
+		else if ( !strcmp( argv[ arg ], "-vmirror" ) )
+			this.vmirror = atoi( argv[ ++ arg ] );
+		else if ( !strcmp( argv[ arg ], "-dmirror" ) )
+			this.dmirror = atoi( argv[ ++ arg ] );
+		else if ( !strcmp( argv[ arg ], "-offset" ) )
+			this.offset = atoi( argv[ ++ arg ] );
+		else if ( !strcmp( argv[ arg ], "-invert" ) )
+			this.invert = atoi( argv[ ++ arg ] );
+		else if ( !strcmp( argv[ arg ], "-flip" ) )
+			this.flip = atoi( argv[ ++ arg ] );
+		else if ( !strcmp( argv[ arg ], "-flop" ) )
+			this.flop = atoi( argv[ ++ arg ] );
+		else if ( !strcmp( argv[ arg ], "-pflip" ) )
+			this.pflip = atoi( argv[ ++ arg ] );
+		else if ( !strcmp( argv[ arg ], "-pflop" ) )
+			this.pflop = atoi( argv[ ++ arg ] );
+		else if ( !strcmp( argv[ arg ], "-quart" ) )
+			this.quart = atoi( argv[ ++ arg ] );
+		else
+			fprintf( stderr, "ignoring %s\n", argv[ arg ] );
+	}
+
+	if ( bpp != 8 && bpp != 16 )
+	{
+		fprintf( stderr, "Invalid bpp %d\n", bpp );
+		return 1;
+	}
+
+	image = luma_render( &this );
+
+	if ( bpp == 16 )
+	{
+		uint16_t *end = image + this.w * this.h;
+		uint16_t *p = image;
+		uint8_t *q = ( uint8_t * )image;
+		while ( p < end )
+		{
+			*p ++ = ( *q << 8 ) + *( q + 1 );
+			q += 2;
+		}
+		printf( "P5\n" );
+		printf( "%d %d\n", this.w, this.h );
+		printf( "65535\n" );
+		fwrite( image, this.w * this.h * sizeof( uint16_t ), 1, stdout );
+	}
+	else
+	{
+		uint16_t *end = image + this.w * this.h;
+		uint16_t *p = image;
+		uint8_t *q = ( uint8_t * )image;
+		while ( p < end )
+			*q ++ = ( uint8_t )( *p ++ >> 8 );
+		printf( "P5\n" );
+		printf( "%d %d\n", this.w, this.h );
+		printf( "255\n" );
+		fwrite( image, this.w * this.h, 1, stdout );
+	}
+
+	return 0;
+}
+
-- 
2.39.2