]> git.sesse.net Git - mlt/blob - src/modules/avformat/filter_avcolour_space.c
Fix regression initializing coefficients.
[mlt] / src / modules / avformat / filter_avcolour_space.c
1 /*
2  * filter_avcolour_space.c -- Colour space filter
3  * Copyright (C) 2004-2005 Ushodaya Enterprises Limited
4  * Author: Charles Yates <charles.yates@pandora.be>
5  *
6  * This library is free software; you can redistribute it and/or
7  * modify it under the terms of the GNU Lesser General Public
8  * License as published by the Free Software Foundation; either
9  * version 2.1 of the License, or (at your option) any later version.
10  *
11  * This library 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 GNU
14  * Lesser General Public License for more details.
15  *
16  * You should have received a copy of the GNU Lesser General Public
17  * License along with this library; if not, write to the Free Software
18  * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
19  */
20
21 #include <framework/mlt_filter.h>
22 #include <framework/mlt_frame.h>
23 #include <framework/mlt_log.h>
24 #include <framework/mlt_profile.h>
25
26 // ffmpeg Header files
27 #include <libavformat/avformat.h>
28 #ifdef SWSCALE
29 #include <libswscale/swscale.h>
30 #endif
31
32 #if LIBAVUTIL_VERSION_INT < (50<<16)
33 #define PIX_FMT_YUYV422 PIX_FMT_YUV422
34 #endif
35
36 #include <stdio.h>
37 #include <stdlib.h>
38
39 #if 0 // This test might come in handy elsewhere someday.
40 static int is_big_endian( )
41 {
42         union { int i; char c[ 4 ]; } big_endian_test;
43         big_endian_test.i = 1;
44
45         return big_endian_test.c[ 0 ] != 1;
46 }
47 #endif
48
49 static int convert_mlt_to_av_cs( mlt_image_format format )
50 {
51         int value = 0;
52
53         switch( format )
54         {
55                 case mlt_image_rgb24:
56                         value = PIX_FMT_RGB24;
57                         break;
58                 case mlt_image_rgb24a:
59                 case mlt_image_opengl:
60                         value = PIX_FMT_RGBA;
61                         break;
62                 case mlt_image_yuv422:
63                         value = PIX_FMT_YUYV422;
64                         break;
65                 case mlt_image_yuv420p:
66                         value = PIX_FMT_YUV420P;
67                         break;
68                 case mlt_image_none:
69                         mlt_log_error( NULL, "[filter avcolor_space] Invalid format\n" );
70                         break;
71         }
72
73         return value;
74 }
75
76 static void set_luma_transfer( struct SwsContext *context, int colorspace, int use_full_range )
77 {
78 #if defined(SWSCALE) && (LIBSWSCALE_VERSION_INT >= ((0<<16)+(7<<8)+2))
79         int *coefficients;
80         const int *new_coefficients;
81         int full_range;
82         int brightness, contrast, saturation;
83
84         if ( sws_getColorspaceDetails( context, &coefficients, &full_range, &coefficients, &full_range,
85                         &brightness, &contrast, &saturation ) != -1 )
86         {
87                 // Don't change these from defaults unless explicitly told to.
88                 if ( use_full_range >= 0 )
89                         full_range = use_full_range;
90                 switch ( colorspace )
91                 {
92                 case 170:
93                 case 470:
94                 case 601:
95                 case 624:
96                         new_coefficients = sws_getCoefficients( SWS_CS_ITU601 );
97                         break;
98                 case 240:
99                         new_coefficients = sws_getCoefficients( SWS_CS_SMPTE240M );
100                         break;
101                 case 709:
102                         new_coefficients = sws_getCoefficients( SWS_CS_ITU709 );
103                         break;
104                 default:
105                         new_coefficients = coefficients;
106                         break;
107                 }
108                 sws_setColorspaceDetails( context, new_coefficients, full_range, new_coefficients, full_range,
109                         brightness, contrast, saturation );
110         }
111 #endif
112 }
113
114 static void av_convert_image( uint8_t *out, uint8_t *in, int out_fmt, int in_fmt,
115         int width, int height, int colorspace, int use_full_range )
116 {
117         AVPicture input;
118         AVPicture output;
119 #ifdef SWSCALE
120         int flags = SWS_BILINEAR | SWS_ACCURATE_RND;
121
122         if ( out_fmt == PIX_FMT_YUYV422 )
123                 flags |= SWS_FULL_CHR_H_INP;
124         else
125                 flags |= SWS_FULL_CHR_H_INT;
126 #ifdef USE_MMX
127         flags |= SWS_CPU_CAPS_MMX;
128 #endif
129 #ifdef USE_SSE
130         flags |= SWS_CPU_CAPS_MMX2;
131 #endif
132 #endif /* SWSCALE */
133
134         avpicture_fill( &input, in, in_fmt, width, height );
135         avpicture_fill( &output, out, out_fmt, width, height );
136 #ifdef SWSCALE
137         struct SwsContext *context = sws_getContext( width, height, in_fmt,
138                 width, height, out_fmt, flags, NULL, NULL, NULL);
139         set_luma_transfer( context, colorspace, use_full_range );
140         sws_scale( context, (const uint8_t* const*) input.data, input.linesize, 0, height,
141                 output.data, output.linesize);
142         sws_freeContext( context );
143 #else
144         img_convert( &output, out_fmt, &input, in_fmt, width, height );
145 #endif
146 }
147
148 /** Do it :-).
149 */
150
151 static int convert_image( mlt_frame frame, uint8_t **image, mlt_image_format *format, mlt_image_format output_format )
152 {
153         mlt_properties properties = MLT_FRAME_PROPERTIES( frame );
154         int width = mlt_properties_get_int( properties, "width" );
155         int height = mlt_properties_get_int( properties, "height" );
156         int error = 0;
157
158         if ( *format != output_format )
159         {
160                 int colorspace = mlt_properties_get_int( properties, "colorspace" );
161                 int force_full_luma = -1;
162                 
163                 mlt_log_debug( NULL, "[filter avcolor_space] %s -> %s @ %dx%d space %d\n",
164                         mlt_image_format_name( *format ), mlt_image_format_name( output_format ),
165                         width, height, colorspace );
166
167                 int in_fmt = convert_mlt_to_av_cs( *format );
168                 int out_fmt = convert_mlt_to_av_cs( output_format );
169                 int size = avpicture_get_size( out_fmt, width, height );
170                 uint8_t *output = mlt_pool_alloc( size );
171
172                 if ( *format == mlt_image_rgb24a || *format == mlt_image_opengl )
173                 {
174                         register int len = width * height;
175                         uint8_t *alpha = mlt_pool_alloc( len );
176
177                         if ( alpha )
178                         {
179                                 // Extract the alpha mask from the RGBA image using Duff's Device
180                                 register uint8_t *s = *image + 3; // start on the alpha component
181                                 register uint8_t *d = alpha;
182                                 register int n = ( len + 7 ) / 8;
183
184                                 switch ( len % 8 )
185                                 {
186                                         case 0: do { *d++ = *s; s += 4;
187                                         case 7:          *d++ = *s; s += 4;
188                                         case 6:          *d++ = *s; s += 4;
189                                         case 5:          *d++ = *s; s += 4;
190                                         case 4:          *d++ = *s; s += 4;
191                                         case 3:          *d++ = *s; s += 4;
192                                         case 2:          *d++ = *s; s += 4;
193                                         case 1:          *d++ = *s; s += 4;
194                                                         }
195                                                         while ( --n > 0 );
196                                 }
197                                 mlt_frame_set_alpha( frame, alpha, len, mlt_pool_release );
198                                 frame->get_alpha_mask = NULL;
199                         }
200                 }
201
202                 // Update the output
203                 if ( *format == mlt_image_yuv422 && mlt_properties_get( properties, "force_full_luma" )
204                      && ( output_format == mlt_image_rgb24 || output_format == mlt_image_rgb24a ) )
205                 {
206                         // By removing the frame property we only permit the luma to skip scaling once.
207                         // Thereafter, we let swscale scale the luma range as it pleases since it seems
208                         // we do not have control over the RGB to YUV conversion.
209                         force_full_luma = mlt_properties_get_int( properties, "force_full_luma" );                      
210                         mlt_properties_set( properties, "force_full_luma", NULL );
211                 }
212                 av_convert_image( output, *image, out_fmt, in_fmt, width, height, colorspace, force_full_luma );
213                 *image = output;
214                 *format = output_format;
215                 mlt_frame_set_image( frame, output, size, mlt_pool_release );
216                 mlt_properties_set_int( properties, "format", output_format );
217
218                 if ( output_format == mlt_image_rgb24a || output_format == mlt_image_opengl )
219                 {
220                         register int len = width * height;
221                         int alpha_size = 0;
222                         uint8_t *alpha = mlt_frame_get_alpha_mask( frame );
223                         mlt_properties_get_data( properties, "alpha", &alpha_size );
224
225                         if ( alpha && alpha_size >= len )
226                         {
227                                 // Merge the alpha mask from into the RGBA image using Duff's Device
228                                 register uint8_t *s = alpha;
229                                 register uint8_t *d = *image + 3; // start on the alpha component
230                                 register int n = ( len + 7 ) / 8;
231
232                                 switch ( len % 8 )
233                                 {
234                                         case 0: do { *d = *s++; d += 4;
235                                         case 7:          *d = *s++; d += 4;
236                                         case 6:          *d = *s++; d += 4;
237                                         case 5:          *d = *s++; d += 4;
238                                         case 4:          *d = *s++; d += 4;
239                                         case 3:          *d = *s++; d += 4;
240                                         case 2:          *d = *s++; d += 4;
241                                         case 1:          *d = *s++; d += 4;
242                                                         }
243                                                         while ( --n > 0 );
244                                 }
245                         }
246                 }
247         }
248         return error;
249 }
250
251 /* TODO: The below is not working because swscale does not have
252  * adjustable coefficients yet for RGB->YUV */
253 #if 0
254 static int get_image( mlt_frame frame, uint8_t **image, mlt_image_format *format, int *width, int *height, int writable )
255 {
256         int error = 0;
257         mlt_profile profile = (mlt_profile) mlt_frame_pop_get_image( frame );
258         mlt_properties properties = MLT_FRAME_PROPERTIES(frame);
259         mlt_image_format format_from = *format;
260         mlt_image_format format_to = mlt_image_rgb24;
261         
262         error = mlt_frame_get_image( frame, image, format, width, height, writable );
263         
264         int frame_colorspace = mlt_properties_get_int( properties, "colorspace" );
265         
266         if ( !error && *format == mlt_image_yuv422 && profile->colorspace > 0 &&
267              frame_colorspace > 0 && frame_colorspace != profile->colorspace )
268         {
269                 mlt_log_debug( NULL, "[filter avcolor_space] colorspace %d -> %d\n",
270                         frame_colorspace, profile->colorspace );
271                 
272                 // Convert to RGB using frame's colorspace
273                 error = convert_image( frame, image, &format_from, format_to );
274
275                 // Convert to YUV using profile's colorspace
276                 if ( !error )
277                 {
278                         *image = mlt_properties_get_data( properties, "image", NULL );
279                         format_from = mlt_image_rgb24;
280                         format_to = *format;
281                         mlt_properties_set_int( properties, "colorspace", profile->colorspace );
282                         error = convert_image( frame, image, &format_from, format_to );
283                         *image = mlt_properties_get_data( properties, "image", NULL );
284                 }
285         }
286         
287         return error;
288 }
289 #endif
290
291 /** Filter processing.
292 */
293
294 static mlt_frame filter_process( mlt_filter filter, mlt_frame frame )
295 {
296         // Set a default colorspace on the frame if not yet set by the producer.
297         // The producer may still change it during get_image.
298         // This way we do not have to modify each producer to set a valid colorspace.
299         mlt_properties properties = MLT_FRAME_PROPERTIES(frame);
300         if ( mlt_properties_get_int( properties, "colorspace" ) <= 0 )
301                 mlt_properties_set_int( properties, "colorspace", mlt_service_profile( MLT_FILTER_SERVICE(filter) )->colorspace );
302
303         frame->convert_image = convert_image;
304     
305 //      Not working yet - see comment for get_image() above.
306 //      mlt_frame_push_service( frame, mlt_service_profile( MLT_FILTER_SERVICE( filter ) ) );
307 //      mlt_frame_push_get_image( frame, get_image );
308
309         return frame;
310 }
311
312 /** Constructor for the filter.
313 */
314
315 mlt_filter filter_avcolour_space_init( void *arg )
316 {
317 #ifdef SWSCALE
318 #if (LIBSWSCALE_VERSION_INT >= ((0<<16)+(7<<8)+2))
319         // Test to see if swscale accepts the arg as resolution
320         if ( arg )
321         {
322                 int width = (int) arg;
323                 struct SwsContext *context = sws_getContext( width, width, PIX_FMT_RGB32, 64, 64, PIX_FMT_RGB32, SWS_BILINEAR, NULL, NULL, NULL);
324                 if ( context )
325                         sws_freeContext( context );
326                 else
327                         return NULL;
328         }               
329 #else
330         return NULL;
331 #endif
332 #endif
333         mlt_filter filter = mlt_filter_new( );
334         if ( filter != NULL )
335                 filter->process = filter_process;
336         return filter;
337 }
338