Optimisations (part 0), pixel v percentage, reworked aspect ratio calcs, ante/post...
[melted] / src / modules / core / filter_obscure.c
1 /*
2 * filter_obscure.c -- obscure filter
3 * Copyright (C) 2003-2004 Ushodaya Enterprises Limited
4 * Author: Charles Yates <charles.yates@pandora.be>
5 *
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.
10 *
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.
15 *
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.
19 */
20
21 #include "filter_obscure.h"
22
23 #include <framework/mlt_frame.h>
24 #include <framework/mlt_deque.h>
25
26 #include <stdio.h>
27 #include <stdlib.h>
28
29 /** Geometry struct.
30 */
31
32 struct geometry_s
33 {
34 int nw;
35 int nh;
36 float x;
37 float y;
38 float w;
39 float h;
40 int mask_w;
41 int mask_h;
42 };
43
44 /** Parse a value from a geometry string.
45 */
46
47 static float parse_value( char **ptr, int normalisation, char delim, float defaults )
48 {
49 float value = defaults;
50
51 if ( *ptr != NULL && **ptr != '\0' )
52 {
53 char *end = NULL;
54 value = strtod( *ptr, &end );
55 if ( end != NULL )
56 {
57 if ( *end == '%' )
58 value = ( value / 100.0 ) * normalisation;
59 while ( *end == delim || *end == '%' )
60 end ++;
61 }
62 *ptr = end;
63 }
64
65 return value;
66 }
67
68 /** Parse a geometry property string.
69 */
70
71 static void geometry_parse( struct geometry_s *geometry, struct geometry_s *defaults, char *property, int nw, int nh )
72 {
73 // Assign normalised width and height
74 geometry->nw = nw;
75 geometry->nh = nh;
76
77 // Assign from defaults if available
78 if ( defaults != NULL )
79 {
80 geometry->x = defaults->x;
81 geometry->y = defaults->y;
82 geometry->w = defaults->w;
83 geometry->h = defaults->h;
84 geometry->mask_w = defaults->mask_w;
85 geometry->mask_h = defaults->mask_h;
86 }
87 else
88 {
89 geometry->mask_w = 20;
90 geometry->mask_h = 20;
91 }
92
93 // Parse the geomtry string
94 if ( property != NULL )
95 {
96 char *ptr = property;
97 geometry->x = parse_value( &ptr, nw, ',', geometry->x );
98 geometry->y = parse_value( &ptr, nh, ':', geometry->y );
99 geometry->w = parse_value( &ptr, nw, 'x', geometry->w );
100 geometry->h = parse_value( &ptr, nh, ':', geometry->h );
101 geometry->mask_w = parse_value( &ptr, nw, 'x', geometry->mask_w );
102 geometry->mask_h = parse_value( &ptr, nh, ' ', geometry->mask_h );
103 }
104 }
105
106 /** A Timism but not as clean ;-).
107 */
108
109 static float lerp( float value, float lower, float upper )
110 {
111 if ( value < lower )
112 return lower;
113 else if ( value > upper )
114 return upper;
115 return value;
116 }
117
118 /** Calculate real geometry.
119 */
120
121 static void geometry_calculate( struct geometry_s *output, struct geometry_s *in, struct geometry_s *out, float position, int ow, int oh )
122 {
123 // Calculate this frames geometry
124 output->x = lerp( ( in->x + ( out->x - in->x ) * position ) / ( float )out->nw * ow, 0, ow );
125 output->y = lerp( ( in->y + ( out->y - in->y ) * position ) / ( float )out->nh * oh, 0, oh );
126 output->w = lerp( ( in->w + ( out->w - in->w ) * position ) / ( float )out->nw * ow, 0, ow - output->x );
127 output->h = lerp( ( in->h + ( out->h - in->h ) * position ) / ( float )out->nh * oh, 0, oh - output->y );
128 output->mask_w = in->mask_w + ( out->mask_w - in->mask_w ) * position;
129 output->mask_h = in->mask_h + ( out->mask_h - in->mask_h ) * position;
130 }
131
132 /** Calculate the position for this frame.
133 */
134
135 static float position_calculate( mlt_filter this, mlt_frame frame )
136 {
137 // Get the in and out position
138 mlt_position in = mlt_filter_get_in( this );
139 mlt_position out = mlt_filter_get_out( this );
140
141 // Get the position of the frame
142 mlt_position position = mlt_frame_get_position( frame );
143
144 // Now do the calcs
145 return ( float )( position - in ) / ( float )( out - in + 1 );
146 }
147
148 /** The averaging function...
149 */
150
151 void obscure_average( uint8_t *start, int width, int height, int stride )
152 {
153 int y;
154 int x;
155 int Y = ( *start + *( start + 2 ) ) / 2;
156 int U = *( start + 1 );
157 int V = *( start + 3 );
158
159 for ( y = 0; y < height; y ++ )
160 {
161 uint8_t *p = start + y * stride;
162 for ( x = 0; x < width / 2; x ++ )
163 {
164 Y = ( Y + *p ++ ) / 2;
165 U = ( U + *p ++ ) / 2;
166 Y = ( Y + *p ++ ) / 2;
167 V = ( V + *p ++ ) / 2;
168 }
169 }
170
171 for ( y = 0; y < height; y ++ )
172 {
173 uint8_t *p = start + y * stride;
174 for ( x = 0; x < width / 2; x ++ )
175 {
176 *p ++ = Y;
177 *p ++ = U;
178 *p ++ = Y;
179 *p ++ = V;
180 }
181 }
182 }
183
184
185 /** The obscurer rendering function...
186 */
187
188 static void obscure_render( uint8_t *image, int width, int height, struct geometry_s result )
189 {
190 int area_x = result.x;
191 int area_y = result.y;
192 int area_w = result.w;
193 int area_h = result.h;
194
195 int mw = result.mask_w;
196 int mh = result.mask_h;
197 int w;
198 int h;
199
200 uint8_t *p = image + area_y * width * 2 + area_x * 2;
201
202 for ( w = 0; w < area_w; w += mw )
203 {
204 for ( h = 0; h < area_h; h += mh )
205 {
206 int aw = w + mw > area_w ? mw - ( w + mw - area_w ) : mw;
207 int ah = h + mh > area_h ? mh - ( h + mh - area_h ) : mh;
208 if ( aw > 1 && ah > 1 );
209 obscure_average( p + h * width * 2 + w * 2, aw, ah, width * 2 );
210 }
211 }
212 }
213
214 /** Do it :-).
215 */
216
217 static int filter_get_image( mlt_frame frame, uint8_t **image, mlt_image_format *format, int *width, int *height, int writable )
218 {
219 // Get the frame properties
220 mlt_properties frame_properties = mlt_frame_properties( frame );
221
222 // Fetch the obscure stack for this frame
223 mlt_deque deque = mlt_properties_get_data( frame_properties, "filter_obscure", NULL );
224
225 // Pop the top of stack now
226 mlt_filter this = mlt_deque_pop_back( deque );
227
228 // Get the image from the frame
229 if ( mlt_frame_get_image( frame, image, format, width, height, 1 ) == 0 )
230 {
231 if ( this != NULL )
232 {
233 // Get the filter properties
234 mlt_properties properties = mlt_filter_properties( this );
235
236 // Obtain the normalised width and height from the frame
237 int normalised_width = mlt_properties_get_int( frame_properties, "normalised_width" );
238 int normalised_height = mlt_properties_get_int( frame_properties, "normalised_height" );
239
240 // Structures for geometry
241 struct geometry_s result;
242 struct geometry_s start;
243 struct geometry_s end;
244
245 // Calculate the position
246 float position = position_calculate( this, frame );
247
248 // Now parse the geometries
249 geometry_parse( &start, NULL, mlt_properties_get( properties, "start" ), normalised_width, normalised_height );
250 geometry_parse( &end, &start, mlt_properties_get( properties, "end" ), normalised_width, normalised_height );
251
252 // Do the calculation
253 geometry_calculate( &result, &start, &end, position, *width, *height );
254
255 // Now actually render it
256 obscure_render( *image, *width, *height, result );
257 }
258 }
259
260 return 0;
261 }
262
263 /** Filter processing.
264 */
265
266 static mlt_frame filter_process( mlt_filter this, mlt_frame frame )
267 {
268 // Get the frame properties
269 mlt_properties frame_properties = mlt_frame_properties( frame );
270
271 // Fetch the obscure stack for this frame
272 mlt_deque deque = mlt_properties_get_data( frame_properties, "filter_obscure", NULL );
273
274 // Create stack if necessary
275 if ( deque == NULL )
276 {
277 // Create the deque
278 deque = mlt_deque_init( );
279
280 // Assign to the frame
281 mlt_properties_set_data( frame_properties, "filter_obscure", deque, 0, ( mlt_destructor )mlt_deque_close, NULL );
282 }
283
284 // Push this on to the obscure stack
285 mlt_deque_push_back( deque, this );
286
287 // Push the get image call
288 mlt_frame_push_get_image( frame, filter_get_image );
289
290 return frame;
291 }
292
293 /** Constructor for the filter.
294 */
295
296 mlt_filter filter_obscure_init( void *arg )
297 {
298 mlt_filter this = calloc( 1, sizeof( struct mlt_filter_s ) );
299 mlt_properties properties = mlt_filter_properties( this );
300 mlt_filter_init( this, NULL );
301 this->process = filter_process;
302 mlt_properties_set( properties, "start", arg != NULL ? arg : "0%,0%:100%x100%" );
303 mlt_properties_set( properties, "end", "" );
304 return this;
305 }
306