2 * transition_luma.c -- a generic dissolve/wipe processor
3 * Copyright (C) 2003-2004 Ushodaya Enterprises Limited
4 * Author: Dan Dennedy <dan@dennedy.org>
6 * This program is free software; you can redistribute it and/or modify
7 * it under the terms of the GNU General Public License as published by
8 * the Free Software Foundation; either version 2 of the License, or
9 * (at your option) any later version.
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.
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., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
21 #include "transition_luma.h"
22 #include <framework/mlt_frame.h>
34 struct mlt_transition_s parent
;
43 // forward declarations
44 static void transition_close( mlt_transition parent
);
47 // image processing functions
49 static inline double smoothstep( double edge1
, double edge2
, double a
)
57 a
= ( a
- edge1
) / ( edge2
- edge1
);
59 return ( a
* a
* ( 3 - 2 * a
) );
64 \param field_order -1 = progressive, 0 = lower field first, 1 = top field first
66 static void luma_composite( mlt_frame
this, mlt_frame b_frame
, int luma_width
, int luma_height
,
67 double *luma_bitmap
, double pos
, double frame_delta
, double softness
, int field_order
)
69 int width_src
, height_src
;
70 int width_dest
, height_dest
;
71 mlt_image_format format_src
, format_dest
;
72 uint8_t *p_src
, *p_dest
;
79 format_src
= mlt_image_yuv422
;
80 format_dest
= mlt_image_yuv422
;
82 mlt_frame_get_image( this, &p_dest
, &format_dest
, &width_dest
, &height_dest
, 1 /* writable */ );
83 mlt_frame_get_image( b_frame
, &p_src
, &format_src
, &width_src
, &height_src
, 0 /* writable */ );
85 stride_src
= width_src
* 2;
86 stride_dest
= width_dest
* 2;
88 // composite using luma map
89 for ( field
= 0; field
< ( field_order
< 0 ?
1 : 2 ); ++field
)
91 // Offset the position based on which field we're looking at ...
92 double field_pos
= pos
+ ( ( field_order
== 0 ?
1 - field
: field
) * frame_delta
* 0.5 );
94 // adjust the position for the softness level
95 field_pos
*= ( 1.0 + softness
);
97 for ( i
= field
; i
< height_src
; i
+= ( field_order
< 0 ?
1 : 2 ) )
99 uint8_t *p
= &p_src
[ i
* stride_src
];
100 uint8_t *q
= &p_dest
[ i
* stride_dest
];
101 uint8_t *o
= &p_dest
[ i
* stride_dest
];
102 double *l
= &luma_bitmap
[ i
* luma_width
];
104 for ( j
= 0; j
< width_src
; j
++ )
109 double value
= smoothstep( weight
, weight
+ softness
, field_pos
);
111 *o
++ = (uint8_t)( y
* value
+ *q
++ * ( 1 - value
) );
112 *o
++ = (uint8_t)( uv
* value
+ *q
++ * ( 1 - value
) );
121 static int transition_get_image( mlt_frame
this, uint8_t **image
, mlt_image_format
*format
, int *width
, int *height
, int writable
)
123 // Get the properties of the a frame
124 mlt_properties a_props
= mlt_frame_properties( this );
126 // Get the b frame from the stack
127 mlt_frame b_frame
= mlt_frame_pop_frame( this );
129 // Get the properties of the b frame
130 mlt_properties b_props
= mlt_frame_properties( b_frame
);
132 // Arbitrary composite defaults
133 static double previous_mix
= 0;
137 double *luma_bitmap
= NULL
;
138 double luma_softness
= 0;
140 int top_field_first
= 0;
142 // mix is the offset time value in the duration of the transition
143 // - also used as the mixing level for a dissolve
144 if ( mlt_properties_get( b_props
, "image.mix" ) != NULL
)
145 mix
= mlt_properties_get_double( b_props
, "image.mix" );
147 // (mix - previous_mix) is the animation delta, if backwards reset previous
148 if ( mix
< previous_mix
)
151 // Get the interlace and field properties of the frame
152 if ( mlt_properties_get( b_props
, "progressive" ) != NULL
)
153 progressive
= mlt_properties_get_int( b_props
, "progressive" );
154 if ( mlt_properties_get( b_props
, "top_field_first" ) != NULL
)
155 top_field_first
= mlt_properties_get_int( b_props
, "top_field_first" );
157 // Get the luma map parameters
158 if ( mlt_properties_get( b_props
, "luma.width" ) != NULL
)
159 luma_width
= mlt_properties_get_int( b_props
, "luma.width" );
160 if ( mlt_properties_get( b_props
, "luma.height" ) != NULL
)
161 luma_height
= mlt_properties_get_int( b_props
, "luma.height" );
162 if ( mlt_properties_get( b_props
, "luma.softness" ) != NULL
)
163 luma_softness
= mlt_properties_get_double( b_props
, "luma.softness" );
164 luma_bitmap
= (double*) mlt_properties_get_data( b_props
, "luma.bitmap", NULL
);
166 if ( luma_width
> 0 && luma_height
> 0 && luma_bitmap
!= NULL
)
167 // Composite the frames using a luma map
168 luma_composite( this, b_frame
, luma_width
, luma_height
, luma_bitmap
, mix
, mix
- previous_mix
,
169 luma_softness
, progressive
> 0 ?
-1 : top_field_first
);
171 // Dissolve the frames using the time offset for mix value
172 mlt_frame_composite_yuv( this, b_frame
, 0, 0, mix
);
174 // Extract the a_frame image info
175 *width
= mlt_properties_get_int( a_props
, "width" );
176 *height
= mlt_properties_get_int( a_props
, "height" );
177 *image
= mlt_properties_get_data( a_props
, "image", NULL
);
184 /** Load the luma map from PGM stream.
187 static void luma_read_pgm( FILE *f
, double **map
, int *width
, int *height
)
189 uint8_t *data
= NULL
;
200 // get the magic code
201 if ( fgets( line
, 127, f
) == NULL
)
203 if ( line
[0] != 'P' || line
[1] != '5' )
206 // skip white space and see if a new line must be fetched
207 for ( i
= 2; i
< 127 && line
[i
] != '\0' && isspace( line
[i
] ); i
++ );
208 if ( line
[i
] == '\0' && fgets( line
, 127, f
) == NULL
)
211 // get the dimensions
212 if ( line
[0] == 'P' )
213 i
= sscanf( line
, "P5 %d %d %d", width
, height
, &maxval
);
215 i
= sscanf( line
, "%d %d %d", width
, height
, &maxval
);
217 // get the height value, if not yet
220 if ( fgets( line
, 127, f
) == NULL
)
222 i
= sscanf( line
, "%d", height
);
229 // get the maximum gray value, if not yet
232 if ( fgets( line
, 127, f
) == NULL
)
234 i
= sscanf( line
, "%d", &maxval
);
239 // determine if this is one or two bytes per pixel
240 bpp
= maxval
> 255 ?
2 : 1;
241 // allocate temporary storage for the raw data
242 data
= malloc( *width
* *height
* bpp
);
247 if ( fread( data
, *width
* *height
* bpp
, 1, f
) != 1 )
250 // allocate the luma bitmap
251 *map
= p
= (double*) malloc( *width
* *height
* sizeof( double ) );
255 // proces the raw data into the luma bitmap
256 for ( i
= 0; i
< *width
* *height
* bpp
; i
+= bpp
)
259 *p
++ = (double) data
[ i
] / (double) maxval
;
261 *p
++ = (double) ( ( data
[ i
] << 8 ) + data
[ i
+1 ] ) / (double) maxval
;
272 /** Luma transition processing.
275 static mlt_frame
transition_process( mlt_transition transition
, mlt_frame a_frame
, mlt_frame b_frame
)
277 transition_luma
*this = (transition_luma
*) transition
->child
;
279 // Get the properties of the transition
280 mlt_properties properties
= mlt_transition_properties( transition
);
282 // Get the properties of the b frame
283 mlt_properties b_props
= mlt_frame_properties( b_frame
);
285 // If the filename property changed, reload the map
286 char *luma_file
= mlt_properties_get( properties
, "filename" );
287 if ( luma_file
!= NULL
&& ( this->filename
== NULL
|| ( this->filename
&& strcmp( luma_file
, this->filename
) ) ) )
289 int width
= mlt_properties_get_int( b_props
, "width" );
290 int height
= mlt_properties_get_int( b_props
, "height" );
294 command
[ 511 ] = '\0';
295 if ( this->filename
)
296 free( this->filename
);
297 this->filename
= strdup( luma_file
);
298 snprintf( command
, 511, "anytopnm %s | pnmscale -width %d -height %d", luma_file
, width
, height
);
299 //pipe = popen( command, "r" );
300 pipe
= fopen( luma_file
, "r" );
304 free( this->bitmap
);
305 luma_read_pgm( pipe
, &this->bitmap
, &this->width
, &this->height
);
312 // Determine the time position of this frame in the transition duration
313 mlt_position in
= mlt_transition_get_in( transition
);
314 mlt_position out
= mlt_transition_get_out( transition
);
315 mlt_position time
= mlt_frame_get_position( b_frame
);
316 double pos
= ( double )( time
- in
) / ( double )( out
- in
+ 1 );
318 // Set the b frame properties
319 mlt_properties_set_double( b_props
, "image.mix", pos
);
320 mlt_properties_set_int( b_props
, "luma.width", this->width
);
321 mlt_properties_set_int( b_props
, "luma.height", this->height
);
322 mlt_properties_set_data( b_props
, "luma.bitmap", this->bitmap
, 0, NULL
, NULL
);
323 if ( mlt_properties_get( properties
, "softness" ) != NULL
)
324 mlt_properties_set_double( b_props
, "luma.softness", mlt_properties_get_double( properties
, "softness" ) );
326 mlt_frame_push_get_image( a_frame
, transition_get_image
);
327 mlt_frame_push_frame( a_frame
, b_frame
);
332 /** Constructor for the filter.
335 mlt_transition
transition_luma_init( char *lumafile
)
337 transition_luma
*this = calloc( sizeof( transition_luma
), 1 );
340 mlt_transition transition
= &this->parent
;
341 mlt_transition_init( transition
, this );
342 transition
->process
= transition_process
;
343 transition
->close
= transition_close
;
345 if ( lumafile
!= NULL
)
346 mlt_properties_set( mlt_transition_properties( transition
), "filename", lumafile
);
348 return &this->parent
;
353 /** Close the transition.
356 static void transition_close( mlt_transition parent
)
358 transition_luma
*this = (transition_luma
*) parent
->child
;
361 free( this->bitmap
);
363 if ( this->filename
)
364 free( this->filename
);