2 * filter_deinterlace.c -- deinterlace filter
3 * Copyright (C) 2003-2004 Ushodaya Enterprises Limited
4 * Author: Charles Yates <charles.yates@pandora.be>
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 <framework/mlt_filter.h>
22 #include <framework/mlt_log.h>
23 #include <framework/mlt_producer.h>
24 #include <framework/mlt_events.h>
25 #include "deinterlace.h"
28 #include <framework/mlt_frame.h>
33 #define YADIF_MODE_TEMPORAL_SPATIAL (0)
34 #define YADIF_MODE_TEMPORAL (2)
36 static yadif_filter *init_yadif( int width, int height )
38 yadif_filter *yadif = mlt_pool_alloc( sizeof( *yadif ) );
40 yadif->cpu = 0; // Pure C
42 yadif->cpu |= AVS_CPU_INTEGER_SSE;
45 yadif->cpu |= AVS_CPU_SSE2;
47 // Create intermediate planar planes
48 yadif->yheight = height;
49 yadif->ywidth = width;
50 yadif->uvwidth = yadif->ywidth / 2;
51 yadif->ypitch = ( yadif->ywidth + 15 ) / 16 * 16;
52 yadif->uvpitch = ( yadif->uvwidth + 15 ) / 16 * 16;
53 yadif->ysrc = (unsigned char *) mlt_pool_alloc( yadif->yheight * yadif->ypitch );
54 yadif->usrc = (unsigned char *) mlt_pool_alloc( yadif->yheight * yadif->uvpitch);
55 yadif->vsrc = (unsigned char *) mlt_pool_alloc( yadif->yheight * yadif->uvpitch );
56 yadif->yprev = (unsigned char *) mlt_pool_alloc( yadif->yheight * yadif->ypitch );
57 yadif->uprev = (unsigned char *) mlt_pool_alloc( yadif->yheight * yadif->uvpitch );
58 yadif->vprev = (unsigned char *) mlt_pool_alloc( yadif->yheight * yadif->uvpitch );
59 yadif->ynext = (unsigned char *) mlt_pool_alloc( yadif->yheight * yadif->ypitch );
60 yadif->unext = (unsigned char *) mlt_pool_alloc( yadif->yheight * yadif->uvpitch );
61 yadif->vnext = (unsigned char *) mlt_pool_alloc( yadif->yheight * yadif->uvpitch );
62 yadif->ydest = (unsigned char *) mlt_pool_alloc( yadif->yheight * yadif->ypitch );
63 yadif->udest = (unsigned char *) mlt_pool_alloc( yadif->yheight * yadif->uvpitch );
64 yadif->vdest = (unsigned char *) mlt_pool_alloc( yadif->yheight * yadif->uvpitch );
69 static void close_yadif(yadif_filter *yadif)
71 mlt_pool_release( yadif->ysrc );
72 mlt_pool_release( yadif->usrc );
73 mlt_pool_release( yadif->vsrc );
74 mlt_pool_release( yadif->yprev );
75 mlt_pool_release( yadif->uprev );
76 mlt_pool_release( yadif->vprev );
77 mlt_pool_release( yadif->ynext );
78 mlt_pool_release( yadif->unext );
79 mlt_pool_release( yadif->vnext );
80 mlt_pool_release( yadif->ydest );
81 mlt_pool_release( yadif->udest );
82 mlt_pool_release( yadif->vdest );
83 mlt_pool_release( yadif );
85 #if defined(__GNUC__) && !defined(PIC)
86 // Set SSSE3 bit to cpu
92 "mov %%ecx, %%edx \n\t"\
96 "and $511, %%ebx \n\t"\
97 "or %%edx, %%ebx \n\t"\
98 : "=b"(yadif->cpu) : "p"(yadif->cpu) : "%eax", "%ecx", "%edx");
102 static int deinterlace_yadif( mlt_frame frame, mlt_filter filter, uint8_t **image, mlt_image_format *format, int *width, int *height, int mode )
104 mlt_properties properties = MLT_FRAME_PROPERTIES( frame );
105 mlt_frame previous_frame = mlt_properties_get_data( properties, "previous frame", NULL );
106 uint8_t* previous_image = NULL;
107 int previous_width = *width;
108 int previous_height = *height;
109 mlt_frame next_frame = mlt_properties_get_data( properties, "next frame", NULL );
110 uint8_t* next_image = NULL;
111 int next_width = *width;
112 int next_height = *height;
114 mlt_log_debug( MLT_FILTER_SERVICE(filter), "previous %d current %d next %d\n",
115 previous_frame? mlt_frame_get_position(previous_frame) : -1,
116 mlt_frame_get_position(frame),
117 next_frame? mlt_frame_get_position(next_frame) : -1);
119 if ( !previous_frame || !next_frame )
122 // Get the preceding frame's image
123 int error = mlt_frame_get_image( previous_frame, &previous_image, format, &previous_width, &previous_height, 0 );
124 int progressive = mlt_properties_get_int( MLT_FRAME_PROPERTIES( previous_frame ), "progressive" );
126 // Check that we aren't already progressive
127 if ( !error && previous_image && !progressive )
129 // OK, now we know we have work to do and can request the image in our format
130 *format = mlt_image_yuv422;
131 error = mlt_frame_get_image( previous_frame, &previous_image, format, &previous_width, &previous_height, 0 );
133 // Get the current frame's image
134 *format = mlt_image_yuv422;
135 error = mlt_frame_get_image( frame, image, format, width, height, 0 );
137 if ( !error && *image && *format == mlt_image_yuv422 )
139 // Get the following frame's image
140 error = mlt_frame_get_image( next_frame, &next_image, format, &next_width, &next_height, 0 );
142 if ( !error && next_image && *format == mlt_image_yuv422 )
144 yadif_filter *yadif = init_yadif( *width, *height );
147 const int order = mlt_properties_get_int( properties, "top_field_first" );
148 const int pitch = *width << 1;
149 const int parity = 0;
151 // Convert packed to planar
152 YUY2ToPlanes( *image, pitch, *width, *height, yadif->ysrc,
153 yadif->ypitch, yadif->usrc, yadif->vsrc, yadif->uvpitch, yadif->cpu );
154 YUY2ToPlanes( previous_image, pitch, *width, *height, yadif->yprev,
155 yadif->ypitch, yadif->uprev, yadif->vprev, yadif->uvpitch, yadif->cpu );
156 YUY2ToPlanes( next_image, pitch, *width, *height, yadif->ynext,
157 yadif->ypitch, yadif->unext, yadif->vnext, yadif->uvpitch, yadif->cpu );
159 // Deinterlace each plane
160 filter_plane( mode, yadif->ydest, yadif->ypitch, yadif->yprev, yadif->ysrc,
161 yadif->ynext, yadif->ypitch, *width, *height, parity, order, yadif->cpu);
162 filter_plane( mode, yadif->udest, yadif->uvpitch,yadif->uprev, yadif->usrc,
163 yadif->unext, yadif->uvpitch, *width >> 1, *height, parity, order, yadif->cpu);
164 filter_plane( mode, yadif->vdest, yadif->uvpitch, yadif->vprev, yadif->vsrc,
165 yadif->vnext, yadif->uvpitch, *width >> 1, *height, parity, order, yadif->cpu);
167 // Convert planar to packed
168 YUY2FromPlanes( *image, pitch, *width, *height, yadif->ydest,
169 yadif->ypitch, yadif->udest, yadif->vdest, yadif->uvpitch, yadif->cpu);
171 close_yadif( yadif );
178 // Get the current frame's image
179 error = mlt_frame_get_image( frame, image, format, width, height, 0 );
187 static int filter_get_image( mlt_frame frame, uint8_t **image, mlt_image_format *format, int *width, int *height, int writable )
190 mlt_filter filter = mlt_frame_pop_service( frame );
191 mlt_properties properties = MLT_FRAME_PROPERTIES( frame );
192 int deinterlace = mlt_properties_get_int( properties, "consumer_deinterlace" );
193 // The progressive var should only represent the frame's original state and not the state
194 // as modified by this filter!
195 int progressive = mlt_properties_get_int( properties, "progressive" );
196 // At this point - before image was requested - (progressive == 0) cannot be trusted because
197 // some producers (avformat) do not yet know.
201 // Determine deinterlace method
202 char *method_str = mlt_properties_get( MLT_FILTER_PROPERTIES( filter ), "method" );
203 int method = DEINTERLACE_NONE;
204 char *frame_method_str = mlt_properties_get( properties, "deinterlace_method" );
206 if ( frame_method_str )
207 method_str = frame_method_str;
209 if ( !method_str || strcmp( method_str, "yadif" ) == 0 )
210 method = DEINTERLACE_YADIF;
211 else if ( strcmp( method_str, "yadif-nospatial" ) == 0 )
212 method = DEINTERLACE_YADIF_NOSPATIAL;
213 else if ( strcmp( method_str, "onefield" ) == 0 )
214 method = DEINTERLACE_ONEFIELD;
215 else if ( strcmp( method_str, "linearblend" ) == 0 )
216 method = DEINTERLACE_LINEARBLEND;
217 else if ( strcmp( method_str, "bob" ) == 0 )
218 method = DEINTERLACE_BOB;
219 else if ( strcmp( method_str, "weave" ) == 0 )
220 method = DEINTERLACE_BOB;
221 else if ( strcmp( method_str, "greedy" ) == 0 )
222 method = DEINTERLACE_GREEDY;
224 if ( method == DEINTERLACE_YADIF )
226 error = deinterlace_yadif( frame, filter, image, format, width, height, YADIF_MODE_TEMPORAL_SPATIAL );
228 else if ( method == DEINTERLACE_YADIF_NOSPATIAL )
230 error = deinterlace_yadif( frame, filter, image, format, width, height, YADIF_MODE_TEMPORAL );
232 if ( error || ( method > DEINTERLACE_NONE && method < DEINTERLACE_YADIF ) )
234 // Signal that we no longer need previous and next frames
235 mlt_service service = mlt_properties_get_data( MLT_FILTER_PROPERTIES(filter), "service", NULL );
236 mlt_properties_set_int( MLT_SERVICE_PROPERTIES(service), "_need_previous_next", 0 );
239 method = DEINTERLACE_ONEFIELD;
241 // Get the current frame's image
242 error = mlt_frame_get_image( frame, image, format, width, height, writable );
243 progressive = mlt_properties_get_int( properties, "progressive" );
244 if ( !error && !progressive )
246 // OK, now we know we have work to do and can request the image in our format
247 error = frame->convert_image( frame, image, format, mlt_image_yuv422 );
249 // Check that we aren't already progressive
250 if ( !error && *image && *format == mlt_image_yuv422 )
252 // Deinterlace the image using one of the Xine deinterlacers
253 int image_size = *width * *height * 2;
254 uint8_t *new_image = mlt_pool_alloc( image_size );
256 deinterlace_yuv( new_image, image, *width * 2, *height, method );
257 mlt_frame_set_image( frame, new_image, image_size, mlt_pool_release );
262 else if ( method == DEINTERLACE_NONE )
264 error = mlt_frame_get_image( frame, image, format, width, height, writable );
267 // update progressive flag after having obtained image
268 progressive = mlt_properties_get_int( properties, "progressive" );
270 mlt_log_debug( MLT_FILTER_SERVICE( filter ), "error %d deint %d prog %d fmt %s method %s\n",
271 error, deinterlace, progressive, mlt_image_format_name( *format ), method_str ? method_str : "yadif" );
275 // Make sure that others know the frame is deinterlaced
276 mlt_properties_set_int( properties, "progressive", 1 );
282 error = mlt_frame_get_image( frame, image, format, width, height, writable );
285 if ( !deinterlace || progressive )
287 // Signal that we no longer need previous and next frames
288 mlt_service service = mlt_properties_get_data( MLT_FILTER_PROPERTIES(filter), "service", NULL );
290 mlt_properties_set_int( MLT_SERVICE_PROPERTIES(service), "_need_previous_next", 0 );
296 /** Deinterlace filter processing - this should be lazy evaluation here...
299 static mlt_frame deinterlace_process( mlt_filter filter, mlt_frame frame )
301 // Push filter on to the service stack
302 mlt_frame_push_service( frame, filter );
304 // Push the get_image method on to the stack
305 mlt_frame_push_get_image( frame, filter_get_image );
310 static void on_service_changed( mlt_service owner, mlt_service filter )
312 mlt_service service = mlt_properties_get_data( MLT_SERVICE_PROPERTIES(filter), "service", NULL );
313 mlt_properties_set_int( MLT_SERVICE_PROPERTIES(service), "_need_previous_next", 1 );
316 /** Constructor for the filter.
319 mlt_filter filter_deinterlace_init( mlt_profile profile, mlt_service_type type, const char *id, char *arg )
321 mlt_filter filter = mlt_filter_new( );
322 if ( filter != NULL )
324 filter->process = deinterlace_process;
325 mlt_properties_set( MLT_FILTER_PROPERTIES( filter ), "method", arg );
326 mlt_events_listen( MLT_FILTER_PROPERTIES( filter ), filter, "service-changed", (mlt_listener) on_service_changed );