2 * filter_vignette.c -- vignette filter
3 * Copyright (c) 2007 Marco Gittler <g.marco@freenet.de>
5 * This library is free software; you can redistribute it and/or
6 * modify it under the terms of the GNU Lesser General Public
7 * License as published by the Free Software Foundation; either
8 * version 2.1 of the License, or (at your option) any later version.
10 * This library is distributed in the hope that it will be useful,
11 * but WITHOUT ANY WARRANTY; without even the implied warranty of
12 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13 * Lesser General Public License for more details.
15 * You should have received a copy of the GNU Lesser General Public
16 * License along with this library; if not, write to the Free Software
17 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
20 #include <framework/mlt_filter.h>
21 #include <framework/mlt_frame.h>
22 #include <framework/mlt_geometry.h>
27 #define MIN(a,b) (a<b?a:b)
28 #define MAX(a,b) (a<b?b:a)
30 #define SIGMOD_STEPS 1000
31 #define POSITION_VALUE(p,s,e) (s+((double)(e-s)*p ))
32 //static double pow2[SIGMOD_STEPS];
34 static int filter_get_image( mlt_frame this, uint8_t **image, mlt_image_format *format, int *width, int *height, int writable )
37 mlt_filter filter = mlt_frame_pop_service( this );
38 *format = mlt_image_yuv422;
39 int error = mlt_frame_get_image( this, image, format, width, height, 1 );
41 if ( error == 0 && *image )
43 mlt_position in = mlt_filter_get_in( filter );
44 //mlt_position out = mlt_filter_get_out( filter );
45 mlt_position time = mlt_frame_get_position( this );
47 mlt_geometry geom=mlt_geometry_init();
48 struct mlt_geometry_item_s item;
49 float smooth, radius, cx, cy, opac;
50 char *val=mlt_properties_get(MLT_FILTER_PROPERTIES( filter ), "geometry" );
51 mlt_geometry_parse(geom,val,-1,-1,-1);
52 mlt_geometry_fetch(geom,&item,time-in);
58 mlt_geometry_close(geom);
60 int video_width = *width;
61 int video_height = *height;
66 double max_opac=opac/100.0;
68 for (y=0;y<video_height;y++){
69 int h2_pow2=pow(y-h2,2.0);
70 for (x=0;x<video_width;x++){
71 uint8_t *pix=(*image+y*video_width*2+x*2);
72 int dx=sqrt(h2_pow2+pow(x-w2,2.0));
74 if (radius-smooth>dx){ //center, make not darker
77 else if (radius+smooth<=dx){//max dark after smooth area
80 //double sigx=5.0-10.0*(double)(dx-radius+smooth)/(2.0*smooth);//smooth >10 inner area, <-10 in dark area
81 //delta=pow2[((int)((sigx+10.0)*SIGMOD_STEPS/20.0))];//sigmoidal
82 delta = ((double)(radius+smooth-dx)/(2.0*smooth));//linear
84 delta=MAX(max_opac,delta);
85 *pix=(double)(*pix)*delta;
86 *(pix+1)=((double)(*(pix+1)-127.0)*delta)+127.0;
89 // short a, short b, short c, short d
91 // b: +=blue, -=yellow
99 static mlt_frame filter_process( mlt_filter this, mlt_frame frame )
102 mlt_frame_push_service( frame, this );
103 mlt_frame_push_get_image( frame, filter_get_image );
108 mlt_filter filter_vignette_init( mlt_profile profile, mlt_service_type type, const char *id, char *arg )
110 mlt_filter this = mlt_filter_new( );
115 for (i=-SIGMOD_STEPS/2;i<SIGMOD_STEPS/2;i++){
116 pow2[i+SIGMOD_STEPS/2]=1.0/(1.0+pow(2.0,-((double)i)/((double)SIGMOD_STEPS/20.0)));
120 this->process = filter_process;
121 mlt_properties_set( MLT_FILTER_PROPERTIES( this ), "geometry", "80:50%:50%:50%:0" );
122 //mlt_properties_set( MLT_FILTER_PROPERTIES( this ), "end", "" );