git.sesse.net Git - mlt/blob - src/modules/core/filter_resize.c

   1 /*
   2  * filter_resize.c -- resizing filter
   3  * Copyright (C) 2003-2004 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
  24 #include <stdio.h>
  25 #include <string.h>
  26 #include <stdlib.h>
  27 #include <math.h>
  28
  29 /** Swapbytes inline.
  30 */
  31
  32 static inline void swap_bytes( uint8_t *upper, uint8_t *lower )
  33 {
  34         uint8_t t = *lower;
  35         *lower = *upper;
  36         *upper = t;
  37 }
  38
  39 static uint8_t *resize_alpha( uint8_t *input, int owidth, int oheight, int iwidth, int iheight, uint8_t alpha_value )
  40 {
  41         uint8_t *output = NULL;
  42
  43         if ( input != NULL && ( iwidth != owidth || iheight != oheight ) && ( owidth > 6 && oheight > 6 ) )
  44         {
  45                 uint8_t *out_line;
  46                 int offset_x = ( owidth - iwidth ) / 2;
  47                 int offset_y = ( oheight - iheight ) / 2;
  48                 int iused = iwidth;
  49
  50                 output = mlt_pool_alloc( owidth * oheight );
  51                 memset( output, alpha_value, owidth * oheight );
  52
  53                 offset_x -= offset_x % 2;
  54
  55                 out_line = output + offset_y * owidth;
  56                 out_line += offset_x;
  57
  58                 // Loop for the entirety of our output height.
  59                 while ( iheight -- )
  60                 {
  61                         // We're in the input range for this row.
  62                         memcpy( out_line, input, iused );
  63
  64                         // Move to next input line
  65                         input += iwidth;
  66
  67                         // Move to next output line
  68                         out_line += owidth;
  69                 }
  70         }
  71
  72         return output;
  73 }
  74
  75 static void resize_image( uint8_t *output, int owidth, int oheight, uint8_t *input, int iwidth, int iheight, int bpp )
  76 {
  77         // Calculate strides
  78         int istride = iwidth * bpp;
  79         int ostride = owidth * bpp;
  80         int offset_x = ( owidth - iwidth ) / 2 * bpp;
  81         int offset_y = ( oheight - iheight ) / 2;
  82         uint8_t *in_line = input;
  83         uint8_t *out_line;
  84         int size = owidth * oheight;
  85         uint8_t *p = output;
  86
  87         // Optimisation point
  88         if ( output == NULL || input == NULL || ( owidth <= 6 || oheight <= 6 || iwidth <= 6 || oheight <= 6 ) )
  89         {
  90                 return;
  91         }
  92         else if ( iwidth == owidth && iheight == oheight )
  93         {
  94                 memcpy( output, input, iheight * istride );
  95                 return;
  96         }
  97
  98         if ( bpp == 2 )
  99         {
 100                 while( size -- )
 101                 {
 102                         *p ++ = 16;
 103                         *p ++ = 128;
 104                 }
 105                 offset_x -= offset_x % 4;
 106         }
 107         else
 108         {
 109                 size *= bpp;
 110                 while ( size-- )
 111                         *p ++ = 0;
 112         }
 113
 114         out_line = output + offset_y * ostride;
 115         out_line += offset_x;
 116
 117         // Loop for the entirety of our output height.
 118         while ( iheight -- )
 119         {
 120                 // We're in the input range for this row.
 121                 memcpy( out_line, in_line, istride );
 122
 123                 // Move to next input line
 124                 in_line += istride;
 125
 126                 // Move to next output line
 127                 out_line += ostride;
 128         }
 129 }
 130
 131 /** A padding function for frames - this does not rescale, but simply
 132         resizes.
 133 */
 134
 135 static uint8_t *frame_resize_image( mlt_frame this, int owidth, int oheight, int bpp )
 136 {
 137         // Get properties
 138         mlt_properties properties = MLT_FRAME_PROPERTIES( this );
 139
 140         // Get the input image, width and height
 141         uint8_t *input = mlt_properties_get_data( properties, "image", NULL );
 142         uint8_t *alpha = mlt_frame_get_alpha_mask( this );
 143         int alpha_size = 0;
 144         mlt_properties_get_data( properties, "alpha", &alpha_size );
 145
 146         int iwidth = mlt_properties_get_int( properties, "width" );
 147         int iheight = mlt_properties_get_int( properties, "height" );
 148
 149         // If width and height are correct, don't do anything
 150         if ( iwidth < owidth || iheight < oheight )
 151         {
 152                 uint8_t alpha_value = mlt_properties_get_int( properties, "resize_alpha" );
 153
 154                 // Create the output image
 155                 uint8_t *output = mlt_pool_alloc( owidth * ( oheight + 1 ) * bpp );
 156
 157                 // Call the generic resize
 158                 resize_image( output, owidth, oheight, input, iwidth, iheight, bpp );
 159
 160                 // Now update the frame
 161                 mlt_properties_set_data( properties, "image", output, owidth * ( oheight + 1 ) * bpp, ( mlt_destructor )mlt_pool_release, NULL );
 162                 mlt_properties_set_int( properties, "width", owidth );
 163                 mlt_properties_set_int( properties, "height", oheight );
 164
 165                 // We should resize the alpha too
 166                 if ( alpha && alpha_size >= iwidth * iheight )
 167                 {
 168                         alpha = resize_alpha( alpha, owidth, oheight, iwidth, iheight, alpha_value );
 169                         if ( alpha )
 170                         {
 171                                 mlt_properties_set_data( properties, "alpha", alpha, owidth * oheight, ( mlt_destructor )mlt_pool_release, NULL );
 172                                 this->get_alpha_mask = NULL;
 173                         }
 174                 }
 175
 176                 // Return the output
 177                 return output;
 178         }
 179         // No change, return input
 180         return input;
 181 }
 182
 183 /** Do it :-).
 184 */
 185
 186 static int filter_get_image( mlt_frame this, uint8_t **image, mlt_image_format *format, int *width, int *height, int writable )
 187 {
 188         int error = 0;
 189
 190         // Get the properties from the frame
 191         mlt_properties properties = MLT_FRAME_PROPERTIES( this );
 192
 193         // Pop the top of stack now
 194         mlt_filter filter = mlt_frame_pop_service( this );
 195
 196         // Retrieve the aspect ratio
 197         double aspect_ratio = mlt_deque_pop_back_double( MLT_FRAME_IMAGE_STACK( this ) );
 198
 199         // Correct Width/height if necessary
 200         if ( *width == 0 || *height == 0 )
 201         {
 202                 *width = mlt_properties_get_int( properties, "normalised_width" );
 203                 *height = mlt_properties_get_int( properties, "normalised_height" );
 204         }
 205
 206         // Assign requested width/height from our subordinate
 207         int owidth = *width;
 208         int oheight = *height;
 209
 210         // Check for the special case - no aspect ratio means no problem :-)
 211         if ( aspect_ratio == 0.0 )
 212                 aspect_ratio = mlt_properties_get_double( properties, "consumer_aspect_ratio" );
 213
 214         // Reset the aspect ratio
 215         mlt_properties_set_double( properties, "aspect_ratio", aspect_ratio );
 216
 217         // Hmmm...
 218         char *rescale = mlt_properties_get( properties, "rescale.interp" );
 219         if ( rescale != NULL && !strcmp( rescale, "none" ) )
 220                 return mlt_frame_get_image( this, image, format, width, height, writable );
 221
 222         if ( mlt_properties_get_int( properties, "distort" ) == 0 )
 223         {
 224                 // Normalise the input and out display aspect
 225                 int normalised_width = mlt_properties_get_int( properties, "normalised_width" );
 226                 int normalised_height = mlt_properties_get_int( properties, "normalised_height" );
 227                 int real_width = mlt_properties_get_int( properties, "real_width" );
 228                 int real_height = mlt_properties_get_int( properties, "real_height" );
 229                 if ( real_width == 0 )
 230                         real_width = mlt_properties_get_int( properties, "width" );
 231                 if ( real_height == 0 )
 232                         real_height = mlt_properties_get_int( properties, "height" );
 233                 double input_ar = aspect_ratio * real_width / real_height;
 234                 double output_ar = mlt_properties_get_double( properties, "consumer_aspect_ratio" ) * owidth / oheight;
 235
 236 //              fprintf( stderr, "real %dx%d normalised %dx%d output %dx%d sar %f in-dar %f out-dar %f\n",
 237 //              real_width, real_height, normalised_width, normalised_height, owidth, oheight, aspect_ratio, input_ar, output_ar);
 238
 239                 // Optimised for the input_ar > output_ar case (e.g. widescreen on standard)
 240                 int scaled_width = rint( ( input_ar * normalised_width ) / output_ar );
 241                 int scaled_height = normalised_height;
 242
 243                 // Now ensure that our images fit in the output frame
 244                 if ( scaled_width > normalised_width )
 245                 {
 246                         scaled_width = normalised_width;
 247                         scaled_height = rint( ( output_ar * normalised_height ) / input_ar );
 248                 }
 249
 250                 // Now calculate the actual image size that we want
 251                 owidth = rint( scaled_width * owidth / normalised_width );
 252                 oheight = rint( scaled_height * oheight / normalised_height );
 253
 254                 // Tell frame we have conformed the aspect to the consumer
 255                 mlt_frame_set_aspect_ratio( this, mlt_properties_get_double( properties, "consumer_aspect_ratio" ) );
 256         }
 257
 258         mlt_properties_set_int( properties, "distort", 0 );
 259
 260         // Now pass on the calculations down the line
 261         mlt_properties_set_int( properties, "resize_width", *width );
 262         mlt_properties_set_int( properties, "resize_height", *height );
 263
 264         // Now get the image
 265         if ( *format == mlt_image_yuv422 )
 266                 owidth -= owidth % 2;
 267         error = mlt_frame_get_image( this, image, format, &owidth, &oheight, writable );
 268
 269         if ( error == 0 && *image )
 270         {
 271                 // Get the requested scale operation
 272                 char *op = mlt_properties_get( MLT_FILTER_PROPERTIES( filter ), "scale" );
 273                 int bpp;
 274
 275                 switch ( *format )
 276                 {
 277                         case mlt_image_yuv422:
 278                                 bpp = 2;
 279                                 break;
 280                         case mlt_image_rgb24:
 281                                 bpp = 3;
 282                                 break;
 283                         case mlt_image_rgb24a:
 284                         case mlt_image_opengl:
 285                                 bpp = 4;
 286                                 break;
 287                         default:
 288                                 // XXX: we only know how to resize packed formats
 289                                 return 1;
 290                 }
 291
 292                 // Provides a manual override for misreported field order
 293                 if ( mlt_properties_get( properties, "meta.top_field_first" ) )
 294                         mlt_properties_set_int( properties, "top_field_first", mlt_properties_get_int( properties, "meta.top_field_first" ) );
 295
 296                 // Correct field order if needed
 297                 if ( mlt_properties_get_int( properties, "top_field_first" ) == 1 &&
 298                      mlt_properties_get( properties, "progressive" ) &&
 299                      mlt_properties_get_int( properties, "progressive" ) == 0 )
 300                 {
 301                         // Get the input image, width and height
 302                         int size = owidth * oheight * bpp;
 303                         uint8_t *new_image = mlt_pool_alloc( size );
 304                         mlt_properties_set_data( properties, "image", new_image, size, ( mlt_destructor )mlt_pool_release, NULL );
 305                         uint8_t *ptr = new_image + owidth * bpp;
 306                         memcpy( new_image, *image, owidth * bpp );
 307                         memcpy( ptr, *image, size - owidth * bpp );
 308                         *image = new_image;
 309
 310                         // Set the normalised field order
 311                         mlt_properties_set_int( properties, "top_field_first", 0 );
 312                         mlt_properties_set_int( properties, "meta.top_field_first", 0 );
 313                 }
 314
 315                 if ( !strcmp( op, "affine" ) )
 316                 {
 317                         // TODO: Determine where this is needed and find a different way
 318                         // *image = mlt_frame_rescale_image( this, *width, *height, bpp );
 319                 }
 320                 else if ( strcmp( op, "none" ) != 0 )
 321                 {
 322                         *image = frame_resize_image( this, *width, *height, bpp );
 323                 }
 324                 else
 325                 {
 326                         *width = owidth;
 327                         *height = oheight;
 328                 }
 329         }
 330
 331         return error;
 332 }
 333
 334 /** Filter processing.
 335 */
 336
 337 static mlt_frame filter_process( mlt_filter this, mlt_frame frame )
 338 {
 339         // Store the aspect ratio reported by the source
 340         mlt_deque_push_back_double( MLT_FRAME_IMAGE_STACK( frame ), mlt_frame_get_aspect_ratio( frame ) );
 341
 342         // Push this on to the service stack
 343         mlt_frame_push_service( frame, this );
 344
 345         // Push the get_image method on to the stack
 346         mlt_frame_push_get_image( frame, filter_get_image );
 347
 348         return frame;
 349 }
 350
 351 /** Constructor for the filter.
 352 */
 353
 354 mlt_filter filter_resize_init( mlt_profile profile, mlt_service_type type, const char *id, char *arg )
 355 {
 356         mlt_filter this = calloc( sizeof( struct mlt_filter_s ), 1 );
 357         if ( mlt_filter_init( this, this ) == 0 )
 358         {
 359                 this->process = filter_process;
 360                 mlt_properties_set( MLT_FILTER_PROPERTIES( this ), "scale", arg == NULL ? "off" : arg );
 361         }
 362         return this;
 363 }