git.sesse.net Git - mlt/blob - src/modules/oldfilm/filter_vignette.c

   1 /*
   2  * filter_vignette.c -- vignette filter
   3  * Copyright (c) 2007 Marco Gittler <g.marco@freenet.de>
   4  *
   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.
   9  *
  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.
  14  *
  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
  18  */
  19
  20 #include <framework/mlt_filter.h>
  21 #include <framework/mlt_frame.h>
  22 #include <framework/mlt_geometry.h>
  23
  24 #include <stdio.h>
  25 #include <stdlib.h>
  26 #include <math.h>
  27 #define MIN(a,b) (a<b?a:b)
  28 #define MAX(a,b) (a<b?b:a)
  29
  30 #define SIGMOD_STEPS 1000
  31 #define POSITION_VALUE(p,s,e) (s+((double)(e-s)*p ))
  32 //static double pow2[SIGMOD_STEPS];
  33
  34 static int filter_get_image( mlt_frame this, uint8_t **image, mlt_image_format *format, int *width, int *height, int writable )
  35 {
  36
  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 );
  40
  41         if ( error == 0 && *image )
  42         {
  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 );
  46
  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);
  53                 smooth=item.x;
  54                 radius=item.y;
  55                 cx=item.w;
  56                 cy=item.h;
  57                 opac=item.mix;
  58                 mlt_geometry_close(geom);
  59
  60                 int video_width = *width;
  61                 int video_height = *height;
  62
  63                 int x,y;
  64                 int w2=cx,h2=cy;
  65                 double delta=1.0;
  66                 double max_opac=opac/100.0;
  67
  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));
  73
  74                                 if (radius-smooth>dx){  //center, make not darker
  75                                         continue;
  76                                 }
  77                                 else if (radius+smooth<=dx){//max dark after smooth area
  78                                         delta=0.0;
  79                                 }else{
  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
  83                                 }
  84                                 delta=MAX(max_opac,delta);
  85                                 *pix=(double)(*pix)*delta;
  86                                 *(pix+1)=((double)(*(pix+1)-127.0)*delta)+127.0;
  87                         }
  88                 }
  89                 // short a, short b, short c, short d
  90                 // a= gray val pix 1
  91                 // b: +=blue, -=yellow
  92                 // c: =gray pix 2
  93                 // d: +=red,-=green
  94         }
  95
  96         return error;
  97 }
  98
  99 static mlt_frame filter_process( mlt_filter this, mlt_frame frame )
 100 {
 101
 102         mlt_frame_push_service( frame, this );
 103         mlt_frame_push_get_image( frame, filter_get_image );
 104         return frame;
 105 }
 106
 107
 108 mlt_filter filter_vignette_init( mlt_profile profile, mlt_service_type type, const char *id, char *arg )
 109 {
 110         mlt_filter this = mlt_filter_new( );
 111         //int i=0;
 112         if ( this != NULL )
 113         {
 114                 /*
 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)));
 117                 }
 118                 */
 119
 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", "" );
 123
 124         }
 125         return this;
 126 }
 127
 128