research.m1stereo.tv/gitweb/ - melted/blob - src/framework/mlt_frame.c

   1 /*
   2  * mlt_frame.c -- interface for all frame classes
   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 "config.h"
  22 #include "mlt_frame.h"
  23 #include "mlt_producer.h"
  24 #include "mlt_factory.h"
  25 #include <stdio.h>
  26 #include <stdlib.h>
  27 #include <string.h>
  28 #include <math.h>
  29
  30 /** Constructor for a frame.
  31 */
  32
  33 mlt_frame mlt_frame_init( )
  34 {
  35         // Allocate a frame
  36         mlt_frame this = calloc( sizeof( struct mlt_frame_s ), 1 );
  37
  38         if ( this != NULL )
  39         {
  40                 // Get the normalisation
  41                 char *normalisation = mlt_environment( "MLT_NORMALISATION" );
  42
  43                 // Initialise the properties
  44                 mlt_properties properties = &this->parent;
  45                 mlt_properties_init( properties, this );
  46
  47                 // Set default properties on the frame
  48                 mlt_properties_set_position( properties, "_position", 0.0 );
  49                 mlt_properties_set_data( properties, "image", NULL, 0, NULL, NULL );
  50
  51                 if ( normalisation == NULL || strcmp( normalisation, "NTSC" ) )
  52                 {
  53                         mlt_properties_set_int( properties, "width", 720 );
  54                         mlt_properties_set_int( properties, "height", 576 );
  55                         mlt_properties_set_int( properties, "normalised_width", 720 );
  56                         mlt_properties_set_int( properties, "normalised_height", 576 );
  57                         mlt_properties_set_double( properties, "aspect_ratio", 59.0/54.0 );
  58                 }
  59                 else
  60                 {
  61                         mlt_properties_set_int( properties, "width", 720 );
  62                         mlt_properties_set_int( properties, "height", 480 );
  63                         mlt_properties_set_int( properties, "normalised_width", 720 );
  64                         mlt_properties_set_int( properties, "normalised_height", 480 );
  65                         mlt_properties_set_double( properties, "aspect_ratio", 10.0/11.0 );
  66                 }
  67
  68                 mlt_properties_set_data( properties, "audio", NULL, 0, NULL, NULL );
  69                 mlt_properties_set_data( properties, "alpha", NULL, 0, NULL, NULL );
  70
  71                 // Construct stacks for frames and methods
  72                 this->stack_image = mlt_deque_init( );
  73                 this->stack_audio = mlt_deque_init( );
  74                 this->stack_service = mlt_deque_init( );
  75         }
  76
  77         return this;
  78 }
  79
  80 /** Fetch the frames properties.
  81 */
  82
  83 mlt_properties mlt_frame_properties( mlt_frame this )
  84 {
  85         return this != NULL ? &this->parent : NULL;
  86 }
  87
  88 /** Check if we have a way to derive something other than a test card.
  89 */
  90
  91 int mlt_frame_is_test_card( mlt_frame this )
  92 {
  93         return mlt_deque_count( this->stack_image ) == 0 || mlt_properties_get_int( MLT_FRAME_PROPERTIES( this ), "test_image" );
  94 }
  95
  96 /** Check if we have a way to derive something other than test audio.
  97 */
  98
  99 int mlt_frame_is_test_audio( mlt_frame this )
 100 {
 101         return mlt_deque_count( this->stack_audio ) == 0 || mlt_properties_get_int( MLT_FRAME_PROPERTIES( this ), "test_audio" );
 102 }
 103
 104 /** Get the aspect ratio of the frame.
 105 */
 106
 107 double mlt_frame_get_aspect_ratio( mlt_frame this )
 108 {
 109         return mlt_properties_get_double( MLT_FRAME_PROPERTIES( this ), "aspect_ratio" );
 110 }
 111
 112 /** Set the aspect ratio of the frame.
 113 */
 114
 115 int mlt_frame_set_aspect_ratio( mlt_frame this, double value )
 116 {
 117         return mlt_properties_set_double( MLT_FRAME_PROPERTIES( this ), "aspect_ratio", value );
 118 }
 119
 120 /** Get the position of this frame.
 121 */
 122
 123 mlt_position mlt_frame_get_position( mlt_frame this )
 124 {
 125         int pos = mlt_properties_get_position( MLT_FRAME_PROPERTIES( this ), "_position" );
 126         return pos < 0 ? 0 : pos;
 127 }
 128
 129 /** Set the position of this frame.
 130 */
 131
 132 int mlt_frame_set_position( mlt_frame this, mlt_position value )
 133 {
 134         return mlt_properties_set_position( MLT_FRAME_PROPERTIES( this ), "_position", value );
 135 }
 136
 137 /** Stack a get_image callback.
 138 */
 139
 140 int mlt_frame_push_get_image( mlt_frame this, mlt_get_image get_image )
 141 {
 142         return mlt_deque_push_back( this->stack_image, get_image );
 143 }
 144
 145 /** Pop a get_image callback.
 146 */
 147
 148 mlt_get_image mlt_frame_pop_get_image( mlt_frame this )
 149 {
 150         return mlt_deque_pop_back( this->stack_image );
 151 }
 152
 153 /** Push a frame.
 154 */
 155
 156 int mlt_frame_push_frame( mlt_frame this, mlt_frame that )
 157 {
 158         return mlt_deque_push_back( this->stack_image, that );
 159 }
 160
 161 /** Pop a frame.
 162 */
 163
 164 mlt_frame mlt_frame_pop_frame( mlt_frame this )
 165 {
 166         return mlt_deque_pop_back( this->stack_image );
 167 }
 168
 169 /** Push a service.
 170 */
 171
 172 int mlt_frame_push_service( mlt_frame this, void *that )
 173 {
 174         return mlt_deque_push_back( this->stack_image, that );
 175 }
 176
 177 /** Pop a service.
 178 */
 179
 180 void *mlt_frame_pop_service( mlt_frame this )
 181 {
 182         return mlt_deque_pop_back( this->stack_image );
 183 }
 184
 185 /** Push a service.
 186 */
 187
 188 int mlt_frame_push_service_int( mlt_frame this, int that )
 189 {
 190         return mlt_deque_push_back_int( this->stack_image, that );
 191 }
 192
 193 /** Pop a service.
 194 */
 195
 196 int mlt_frame_pop_service_int( mlt_frame this )
 197 {
 198         return mlt_deque_pop_back_int( this->stack_image );
 199 }
 200
 201 /** Push an audio item on the stack.
 202 */
 203
 204 int mlt_frame_push_audio( mlt_frame this, void *that )
 205 {
 206         return mlt_deque_push_back( this->stack_audio, that );
 207 }
 208
 209 /** Pop an audio item from the stack
 210 */
 211
 212 void *mlt_frame_pop_audio( mlt_frame this )
 213 {
 214         return mlt_deque_pop_back( this->stack_audio );
 215 }
 216
 217 /** Return the service stack
 218 */
 219
 220 mlt_deque mlt_frame_service_stack( mlt_frame this )
 221 {
 222         return this->stack_service;
 223 }
 224
 225 /** Replace image stack with the information provided.
 226
 227         This might prove to be unreliable and restrictive - the idea is that a transition
 228         which normally uses two images may decide to only use the b frame (ie: in the case
 229         of a composite where the b frame completely obscures the a frame).
 230
 231         The image must be writable and the destructor for the image itself must be taken
 232         care of on another frame and that frame cannot have a replace applied to it...
 233         Further it assumes that no alpha mask is in use.
 234
 235         For these reasons, it can only be used in a specific situation - when you have
 236         multiple tracks each with their own transition and these transitions are applied
 237         in a strictly reversed order (ie: highest numbered [lowest track] is processed
 238         first).
 239
 240         More reliable approach - the cases should be detected during the process phase
 241         and the upper tracks should simply not be invited to stack...
 242 */
 243
 244 void mlt_frame_replace_image( mlt_frame this, uint8_t *image, mlt_image_format format, int width, int height )
 245 {
 246         // Remove all items from the stack
 247         while( mlt_deque_pop_back( this->stack_image ) ) ;
 248
 249         // Update the information
 250         mlt_properties_set_data( MLT_FRAME_PROPERTIES( this ), "image", image, 0, NULL, NULL );
 251         mlt_properties_set_int( MLT_FRAME_PROPERTIES( this ), "width", width );
 252         mlt_properties_set_int( MLT_FRAME_PROPERTIES( this ), "height", height );
 253         mlt_properties_set_int( MLT_FRAME_PROPERTIES( this ), "format", format );
 254         this->get_alpha_mask = NULL;
 255 }
 256
 257 /** Get the image associated to the frame.
 258 */
 259
 260 int mlt_frame_get_image( mlt_frame this, uint8_t **buffer, mlt_image_format *format, int *width, int *height, int writable )
 261 {
 262         mlt_properties properties = MLT_FRAME_PROPERTIES( this );
 263         mlt_get_image get_image = mlt_frame_pop_get_image( this );
 264         mlt_producer producer = mlt_properties_get_data( properties, "test_card_producer", NULL );
 265         int error = 0;
 266
 267         if ( get_image != NULL )
 268         {
 269                 mlt_properties_set_int( properties, "image_count", mlt_properties_get_int( properties, "image_count" ) - 1 );
 270                 mlt_position position = mlt_frame_get_position( this );
 271                 error = get_image( this, buffer, format, width, height, writable );
 272                 mlt_properties_set_int( properties, "width", *width );
 273                 mlt_properties_set_int( properties, "height", *height );
 274                 mlt_properties_set_int( properties, "format", *format );
 275                 mlt_frame_set_position( this, position );
 276         }
 277         else if ( mlt_properties_get_data( properties, "image", NULL ) != NULL )
 278         {
 279                 *format = mlt_properties_get_int( properties, "format" );
 280                 *buffer = mlt_properties_get_data( properties, "image", NULL );
 281                 *width = mlt_properties_get_int( properties, "width" );
 282                 *height = mlt_properties_get_int( properties, "height" );
 283         }
 284         else if ( producer != NULL )
 285         {
 286                 mlt_frame test_frame = NULL;
 287                 mlt_service_get_frame( MLT_PRODUCER_SERVICE( producer ), &test_frame, 0 );
 288                 if ( test_frame != NULL )
 289                 {
 290                         mlt_properties test_properties = MLT_FRAME_PROPERTIES( test_frame );
 291                         mlt_properties_set_double( test_properties, "consumer_aspect_ratio", mlt_properties_get_double( properties, "consumer_aspect_ratio" ) );
 292                         mlt_properties_set( test_properties, "rescale.interp", mlt_properties_get( properties, "rescale.interp" ) );
 293                         mlt_frame_get_image( test_frame, buffer, format, width, height, writable );
 294                         mlt_properties_set_data( properties, "test_card_frame", test_frame, 0, ( mlt_destructor )mlt_frame_close, NULL );
 295                         mlt_properties_set_data( properties, "image", *buffer, *width * *height * 2, NULL, NULL );
 296                         mlt_properties_set_int( properties, "width", *width );
 297                         mlt_properties_set_int( properties, "height", *height );
 298                         mlt_properties_set_int( properties, "format", *format );
 299                         mlt_properties_set_double( properties, "aspect_ratio", mlt_frame_get_aspect_ratio( test_frame ) );
 300                 }
 301                 else
 302                 {
 303                         mlt_properties_set_data( properties, "test_card_producer", NULL, 0, NULL, NULL );
 304                         mlt_frame_get_image( this, buffer, format, width, height, writable );
 305                 }
 306         }
 307         else
 308         {
 309                 register uint8_t *p;
 310                 register uint8_t *q;
 311                 int size = 0;
 312
 313                 *width = *width == 0 ? 720 : *width;
 314                 *height = *height == 0 ? 576 : *height;
 315                 size = *width * *height;
 316
 317                 mlt_properties_set_int( properties, "format", *format );
 318                 mlt_properties_set_int( properties, "width", *width );
 319                 mlt_properties_set_int( properties, "height", *height );
 320                 mlt_properties_set_int( properties, "aspect_ratio", 0 );
 321
 322                 switch( *format )
 323                 {
 324                         case mlt_image_none:
 325                                 size = 0;
 326                                 *buffer = NULL;
 327                                 break;
 328                         case mlt_image_rgb24:
 329                                 size *= 3;
 330                                 size += *width * 3;
 331                                 *buffer = mlt_pool_alloc( size );
 332                                 if ( *buffer )
 333                                         memset( *buffer, 255, size );
 334                                 break;
 335                         case mlt_image_rgb24a:
 336                         case mlt_image_opengl:
 337                                 size *= 4;
 338                                 size += *width * 4;
 339                                 *buffer = mlt_pool_alloc( size );
 340                                 if ( *buffer )
 341                                         memset( *buffer, 255, size );
 342                                 break;
 343                         case mlt_image_yuv422:
 344                                 size *= 2;
 345                                 size += *width * 2;
 346                                 *buffer = mlt_pool_alloc( size );
 347                                 p = *buffer;
 348                                 q = p + size;
 349                                 while ( p != NULL && p != q )
 350                                 {
 351                                         *p ++ = 235;
 352                                         *p ++ = 128;
 353                                 }
 354                                 break;
 355                         case mlt_image_yuv420p:
 356                                 size = size * 3 / 2;
 357                                 *buffer = mlt_pool_alloc( size );
 358                                 if ( *buffer )
 359                                         memset( *buffer, 255, size );
 360                                 break;
 361                 }
 362
 363                 mlt_properties_set_data( properties, "image", *buffer, size, ( mlt_destructor )mlt_pool_release, NULL );
 364                 mlt_properties_set_int( properties, "test_image", 1 );
 365         }
 366
 367         mlt_properties_set_int( properties, "scaled_width", *width );
 368         mlt_properties_set_int( properties, "scaled_height", *height );
 369
 370         return error;
 371 }
 372
 373 uint8_t *mlt_frame_get_alpha_mask( mlt_frame this )
 374 {
 375         uint8_t *alpha = NULL;
 376         if ( this != NULL )
 377         {
 378                 if ( this->get_alpha_mask != NULL )
 379                         alpha = this->get_alpha_mask( this );
 380                 if ( alpha == NULL )
 381                         alpha = mlt_properties_get_data( &this->parent, "alpha", NULL );
 382                 if ( alpha == NULL )
 383                 {
 384                         int size = mlt_properties_get_int( &this->parent, "scaled_width" ) * mlt_properties_get_int( &this->parent, "scaled_height" );
 385                         alpha = mlt_pool_alloc( size );
 386                         memset( alpha, 255, size );
 387                         mlt_properties_set_data( &this->parent, "alpha", alpha, size, mlt_pool_release, NULL );
 388                 }
 389         }
 390         return alpha;
 391 }
 392
 393 int mlt_frame_get_audio( mlt_frame this, int16_t **buffer, mlt_audio_format *format, int *frequency, int *channels, int *samples )
 394 {
 395         mlt_get_audio get_audio = mlt_frame_pop_audio( this );
 396         mlt_properties properties = MLT_FRAME_PROPERTIES( this );
 397         int hide = mlt_properties_get_int( properties, "test_audio" );
 398
 399         if ( hide == 0 && get_audio != NULL )
 400         {
 401                 mlt_position position = mlt_frame_get_position( this );
 402                 get_audio( this, buffer, format, frequency, channels, samples );
 403                 mlt_frame_set_position( this, position );
 404         }
 405         else if ( mlt_properties_get_data( properties, "audio", NULL ) )
 406         {
 407                 *buffer = mlt_properties_get_data( properties, "audio", NULL );
 408                 *frequency = mlt_properties_get_int( properties, "audio_frequency" );
 409                 *channels = mlt_properties_get_int( properties, "audio_channels" );
 410                 *samples = mlt_properties_get_int( properties, "audio_samples" );
 411         }
 412         else
 413         {
 414                 int size = 0;
 415                 *samples = *samples <= 0 ? 1920 : *samples;
 416                 *channels = *channels <= 0 ? 2 : *channels;
 417                 *frequency = *frequency <= 0 ? 48000 : *frequency;
 418                 size = *samples * *channels * sizeof( int16_t );
 419                 *buffer = mlt_pool_alloc( size );
 420                 if ( *buffer != NULL )
 421                         memset( *buffer, 0, size );
 422                 mlt_properties_set_data( properties, "audio", *buffer, size, ( mlt_destructor )mlt_pool_release, NULL );
 423                 mlt_properties_set_int( properties, "test_audio", 1 );
 424         }
 425
 426         mlt_properties_set_int( properties, "audio_frequency", *frequency );
 427         mlt_properties_set_int( properties, "audio_channels", *channels );
 428         mlt_properties_set_int( properties, "audio_samples", *samples );
 429
 430         if ( mlt_properties_get( properties, "meta.volume" ) )
 431         {
 432                 double value = mlt_properties_get_double( properties, "meta.volume" );
 433
 434                 if ( value == 0.0 )
 435                 {
 436                         memset( *buffer, 0, *samples * *channels * 2 );
 437                 }
 438                 else if ( value != 1.0 )
 439                 {
 440                         int total = *samples * *channels;
 441                         int16_t *p = *buffer;
 442                         while ( total -- )
 443                         {
 444                                 *p = *p * value;
 445                                 p ++;
 446                         }
 447                 }
 448
 449                 mlt_properties_set( properties, "meta.volume", NULL );
 450         }
 451
 452         return 0;
 453 }
 454
 455 unsigned char *mlt_frame_get_waveform( mlt_frame this, int w, int h )
 456 {
 457         int16_t *pcm = NULL;
 458         mlt_properties properties = MLT_FRAME_PROPERTIES( this );
 459         mlt_audio_format format = mlt_audio_pcm;
 460         int frequency = 32000; // lower frequency available?
 461         int channels = 2;
 462         double fps = mlt_properties_get_double( properties, "fps" );
 463         int samples = mlt_sample_calculator( fps, frequency, mlt_frame_get_position( this ) );
 464
 465         // Get the pcm data
 466         mlt_frame_get_audio( this, &pcm, &format, &frequency, &channels, &samples );
 467
 468         // Make an 8-bit buffer large enough to hold rendering
 469         int size = w * h;
 470         unsigned char *bitmap = ( unsigned char* )mlt_pool_alloc( size );
 471         if ( bitmap != NULL )
 472                 memset( bitmap, 0, size );
 473         mlt_properties_set_data( properties, "waveform", bitmap, size, ( mlt_destructor )mlt_pool_release, NULL );
 474
 475         // Render vertical lines
 476         int16_t *ubound = pcm + samples * channels;
 477         int skip = samples / w - 1;
 478         int i, j, k;
 479
 480         // Iterate sample stream and along x coordinate
 481         for ( i = 0; i < w && pcm < ubound; i++ )
 482         {
 483                 // pcm data has channels interleaved
 484                 for ( j = 0; j < channels; j++ )
 485                 {
 486                         // Determine sample's magnitude from 2s complement;
 487                         int pcm_magnitude = *pcm < 0 ? ~(*pcm) + 1 : *pcm;
 488                         // The height of a line is the ratio of the magnitude multiplied by
 489                         // half the vertical resolution
 490                         int height = ( int )( ( double )( pcm_magnitude ) / 32768 * h / 2 );
 491                         // Determine the starting y coordinate - left channel above center,
 492                         // right channel below - currently assumes 2 channels
 493                         int displacement = ( h / 2 ) - ( 1 - j ) * height;
 494                         // Position buffer pointer using y coordinate, stride, and x coordinate
 495                         unsigned char *p = &bitmap[ i + displacement * w ];
 496
 497                         // Draw vertical line
 498                         for ( k = 0; k < height; k++ )
 499                                 p[ w * k ] = 0xFF;
 500
 501                         pcm++;
 502                 }
 503                 pcm += skip * channels;
 504         }
 505
 506         return bitmap;
 507 }
 508
 509 mlt_producer mlt_frame_get_original_producer( mlt_frame this )
 510 {
 511         if ( this != NULL )
 512                 return mlt_properties_get_data( MLT_FRAME_PROPERTIES( this ), "_producer", NULL );
 513         return NULL;
 514 }
 515
 516 void mlt_frame_close( mlt_frame this )
 517 {
 518         if ( this != NULL && mlt_properties_dec_ref( MLT_FRAME_PROPERTIES( this ) ) <= 0 )
 519         {
 520                 mlt_deque_close( this->stack_image );
 521                 mlt_deque_close( this->stack_audio );
 522                 while( mlt_deque_peek_back( this->stack_service ) )
 523                         mlt_service_close( mlt_deque_pop_back( this->stack_service ) );
 524                 mlt_deque_close( this->stack_service );
 525                 mlt_properties_close( &this->parent );
 526                 free( this );
 527         }
 528 }
 529
 530 /***** convenience functions *****/
 531
 532 int mlt_convert_yuv422_to_rgb24a( uint8_t *yuv, uint8_t *rgba, unsigned int total )
 533 {
 534         int ret = 0;
 535         int yy, uu, vv, ug_plus_vg, ub, vr;
 536         int r,g,b;
 537         total /= 2;
 538         while (total--)
 539         {
 540                 yy = yuv[0] << 8;
 541                 uu = yuv[1] - 128;
 542                 vv = yuv[3] - 128;
 543                 ug_plus_vg = uu * 88 + vv * 183;
 544                 ub = uu * 454;
 545                 vr = vv * 359;
 546                 r = (yy + vr) >> 8;
 547                 g = (yy - ug_plus_vg) >> 8;
 548                 b = (yy + ub) >> 8;
 549                 rgba[0] = r < 0 ? 0 : (r > 255 ? 255 : (unsigned char)r);
 550                 rgba[1] = g < 0 ? 0 : (g > 255 ? 255 : (unsigned char)g);
 551                 rgba[2] = b < 0 ? 0 : (b > 255 ? 255 : (unsigned char)b);
 552                 rgba[3] = 255;
 553                 yy = yuv[2] << 8;
 554                 r = (yy + vr) >> 8;
 555                 g = (yy - ug_plus_vg) >> 8;
 556                 b = (yy + ub) >> 8;
 557                 rgba[4] = r < 0 ? 0 : (r > 255 ? 255 : (unsigned char)r);
 558                 rgba[5] = g < 0 ? 0 : (g > 255 ? 255 : (unsigned char)g);
 559                 rgba[6] = b < 0 ? 0 : (b > 255 ? 255 : (unsigned char)b);
 560                 rgba[7] = 255;
 561                 yuv += 4;
 562                 rgba += 8;
 563         }
 564         return ret;
 565 }
 566
 567 int mlt_convert_rgb24a_to_yuv422( uint8_t *rgba, int width, int height, int stride, uint8_t *yuv, uint8_t *alpha )
 568 {
 569         int ret = 0;
 570         register int y0, y1, u0, u1, v0, v1;
 571         register int r, g, b;
 572         register uint8_t *d = yuv;
 573         register int i, j;
 574
 575         for ( i = 0; i < height; i++ )
 576         {
 577                 register uint8_t *s = rgba + ( stride * i );
 578                 for ( j = 0; j < ( width / 2 ); j++ )
 579                 {
 580                         r = *s++;
 581                         g = *s++;
 582                         b = *s++;
 583                         *alpha++ = *s++;
 584                         RGB2YUV (r, g, b, y0, u0 , v0);
 585                         r = *s++;
 586                         g = *s++;
 587                         b = *s++;
 588                         *alpha++ = *s++;
 589                         RGB2YUV (r, g, b, y1, u1 , v1);
 590                         *d++ = y0;
 591                         *d++ = (u0+u1) >> 1;
 592                         *d++ = y1;
 593                         *d++ = (v0+v1) >> 1;
 594                 }
 595                 if ( width % 2 )
 596                 {
 597                         r = *s++;
 598                         g = *s++;
 599                         b = *s++;
 600                         *alpha++ = *s++;
 601                         RGB2YUV (r, g, b, y0, u0 , v0);
 602                         *d++ = y0;
 603                         *d++ = u0;
 604                 }
 605         }
 606         return ret;
 607 }
 608
 609 int mlt_convert_rgb24_to_yuv422( uint8_t *rgb, int width, int height, int stride, uint8_t *yuv )
 610 {
 611         int ret = 0;
 612         register int y0, y1, u0, u1, v0, v1;
 613         register int r, g, b;
 614         register uint8_t *d = yuv;
 615         register int i, j;
 616
 617         for ( i = 0; i < height; i++ )
 618         {
 619                 register uint8_t *s = rgb + ( stride * i );
 620                 for ( j = 0; j < ( width / 2 ); j++ )
 621                 {
 622                         r = *s++;
 623                         g = *s++;
 624                         b = *s++;
 625                         RGB2YUV (r, g, b, y0, u0 , v0);
 626                         r = *s++;
 627                         g = *s++;
 628                         b = *s++;
 629                         RGB2YUV (r, g, b, y1, u1 , v1);
 630                         *d++ = y0;
 631                         *d++ = (u0+u1) >> 1;
 632                         *d++ = y1;
 633                         *d++ = (v0+v1) >> 1;
 634                 }
 635                 if ( width % 2 )
 636                 {
 637                         r = *s++;
 638                         g = *s++;
 639                         b = *s++;
 640                         RGB2YUV (r, g, b, y0, u0 , v0);
 641                         *d++ = y0;
 642                         *d++ = u0;
 643                 }
 644         }
 645         return ret;
 646 }
 647
 648 int mlt_convert_bgr24a_to_yuv422( uint8_t *rgba, int width, int height, int stride, uint8_t *yuv, uint8_t *alpha )
 649 {
 650         int ret = 0;
 651         register int y0, y1, u0, u1, v0, v1;
 652         register int r, g, b;
 653         register uint8_t *d = yuv;
 654         register int i, j;
 655
 656         for ( i = 0; i < height; i++ )
 657         {
 658                 register uint8_t *s = rgba + ( stride * i );
 659                 for ( j = 0; j < ( width / 2 ); j++ )
 660                 {
 661                         b = *s++;
 662                         g = *s++;
 663                         r = *s++;
 664                         *alpha++ = *s++;
 665                         RGB2YUV (r, g, b, y0, u0 , v0);
 666                         b = *s++;
 667                         g = *s++;
 668                         r = *s++;
 669                         *alpha++ = *s++;
 670                         RGB2YUV (r, g, b, y1, u1 , v1);
 671                         *d++ = y0;
 672                         *d++ = (u0+u1) >> 1;
 673                         *d++ = y1;
 674                         *d++ = (v0+v1) >> 1;
 675                 }
 676                 if ( width % 2 )
 677                 {
 678                         b = *s++;
 679                         g = *s++;
 680                         r = *s++;
 681                         *alpha++ = *s++;
 682                         RGB2YUV (r, g, b, y0, u0 , v0);
 683                         *d++ = y0;
 684                         *d++ = u0;
 685                 }
 686         }
 687         return ret;
 688 }
 689
 690 int mlt_convert_bgr24_to_yuv422( uint8_t *rgb, int width, int height, int stride, uint8_t *yuv )
 691 {
 692         int ret = 0;
 693         register int y0, y1, u0, u1, v0, v1;
 694         register int r, g, b;
 695         register uint8_t *d = yuv;
 696         register int i, j;
 697
 698         for ( i = 0; i < height; i++ )
 699         {
 700                 register uint8_t *s = rgb + ( stride * i );
 701                 for ( j = 0; j < ( width / 2 ); j++ )
 702                 {
 703                         b = *s++;
 704                         g = *s++;
 705                         r = *s++;
 706                         RGB2YUV (r, g, b, y0, u0 , v0);
 707                         b = *s++;
 708                         g = *s++;
 709                         r = *s++;
 710                         RGB2YUV (r, g, b, y1, u1 , v1);
 711                         *d++ = y0;
 712                         *d++ = (u0+u1) >> 1;
 713                         *d++ = y1;
 714                         *d++ = (v0+v1) >> 1;
 715                 }
 716                 if ( width % 2 )
 717                 {
 718                         b = *s++;
 719                         g = *s++;
 720                         r = *s++;
 721                         RGB2YUV (r, g, b, y0, u0 , v0);
 722                         *d++ = y0;
 723                         *d++ = u0;
 724                 }
 725         }
 726         return ret;
 727 }
 728
 729 int mlt_convert_argb_to_yuv422( uint8_t *rgba, int width, int height, int stride, uint8_t *yuv, uint8_t *alpha )
 730 {
 731         int ret = 0;
 732         register int y0, y1, u0, u1, v0, v1;
 733         register int r, g, b;
 734         register uint8_t *d = yuv;
 735         register int i, j;
 736
 737         for ( i = 0; i < height; i++ )
 738         {
 739                 register uint8_t *s = rgba + ( stride * i );
 740                 for ( j = 0; j < ( width / 2 ); j++ )
 741                 {
 742                         *alpha++ = *s++;
 743                         r = *s++;
 744                         g = *s++;
 745                         b = *s++;
 746                         RGB2YUV (r, g, b, y0, u0 , v0);
 747                         *alpha++ = *s++;
 748                         r = *s++;
 749                         g = *s++;
 750                         b = *s++;
 751                         RGB2YUV (r, g, b, y1, u1 , v1);
 752                         *d++ = y0;
 753                         *d++ = (u0+u1) >> 1;
 754                         *d++ = y1;
 755                         *d++ = (v0+v1) >> 1;
 756                 }
 757                 if ( width % 2 )
 758                 {
 759                         *alpha++ = *s++;
 760                         r = *s++;
 761                         g = *s++;
 762                         b = *s++;
 763                         RGB2YUV (r, g, b, y0, u0 , v0);
 764                         *d++ = y0;
 765                         *d++ = u0;
 766                 }
 767         }
 768         return ret;
 769 }
 770
 771 int mlt_convert_yuv420p_to_yuv422( uint8_t *yuv420p, int width, int height, int stride, uint8_t *yuv )
 772 {
 773         int ret = 0;
 774         register int i, j;
 775
 776         int half = width >> 1;
 777
 778         uint8_t *Y = yuv420p;
 779         uint8_t *U = Y + width * height;
 780         uint8_t *V = U + width * height / 4;
 781
 782         register uint8_t *d = yuv;
 783
 784         for ( i = 0; i < height; i++ )
 785         {
 786                 register uint8_t *u = U + ( i / 2 ) * ( half );
 787                 register uint8_t *v = V + ( i / 2 ) * ( half );
 788
 789                 for ( j = 0; j < half; j++ )
 790                 {
 791                         *d ++ = *Y ++;
 792                         *d ++ = *u ++;
 793                         *d ++ = *Y ++;
 794                         *d ++ = *v ++;
 795                 }
 796         }
 797         return ret;
 798 }
 799
 800 uint8_t *mlt_resize_alpha( uint8_t *input, int owidth, int oheight, int iwidth, int iheight, uint8_t alpha_value )
 801 {
 802         uint8_t *output = NULL;
 803
 804         if ( input != NULL && ( iwidth != owidth || iheight != oheight ) && ( owidth > 6 && oheight > 6 ) )
 805         {
 806                 uint8_t *out_line;
 807                 int offset_x = ( owidth - iwidth ) / 2;
 808                 int offset_y = ( oheight - iheight ) / 2;
 809                 int iused = iwidth;
 810
 811                 output = mlt_pool_alloc( owidth * oheight );
 812                 memset( output, alpha_value, owidth * oheight );
 813
 814                 offset_x -= offset_x % 2;
 815
 816                 out_line = output + offset_y * owidth;
 817                 out_line += offset_x;
 818
 819                 // Loop for the entirety of our output height.
 820                 while ( iheight -- )
 821                 {
 822                         // We're in the input range for this row.
 823                         memcpy( out_line, input, iused );
 824
 825                         // Move to next input line
 826                         input += iwidth;
 827
 828                 // Move to next output line
 829                 out_line += owidth;
 830                 }
 831         }
 832
 833         return output;
 834 }
 835
 836 void mlt_resize_yuv422( uint8_t *output, int owidth, int oheight, uint8_t *input, int iwidth, int iheight )
 837 {
 838         // Calculate strides
 839         int istride = iwidth * 2;
 840         int ostride = owidth * 2;
 841         int offset_x = ( owidth - iwidth );
 842         int offset_y = ( oheight - iheight ) / 2;
 843         uint8_t *in_line = input;
 844         uint8_t *out_line;
 845         int size = owidth * oheight;
 846         uint8_t *p = output;
 847
 848         // Optimisation point
 849         if ( output == NULL || input == NULL || ( owidth <= 6 || oheight <= 6 || iwidth <= 6 || oheight <= 6 ) )
 850         {
 851                 return;
 852         }
 853         else if ( iwidth == owidth && iheight == oheight )
 854         {
 855                 memcpy( output, input, iheight * istride );
 856                 return;
 857         }
 858
 859         while( size -- )
 860         {
 861                 *p ++ = 16;
 862                 *p ++ = 128;
 863         }
 864
 865         offset_x -= offset_x % 4;
 866
 867         out_line = output + offset_y * ostride;
 868         out_line += offset_x;
 869
 870         // Loop for the entirety of our output height.
 871         while ( iheight -- )
 872         {
 873                 // We're in the input range for this row.
 874                 memcpy( out_line, in_line, iwidth * 2 );
 875
 876                 // Move to next input line
 877                 in_line += istride;
 878
 879                 // Move to next output line
 880                 out_line += ostride;
 881         }
 882 }
 883
 884 /** A resizing function for yuv422 frames - this does not rescale, but simply
 885         resizes. It assumes yuv422 images available on the frame so use with care.
 886 */
 887
 888 uint8_t *mlt_frame_resize_yuv422( mlt_frame this, int owidth, int oheight )
 889 {
 890         // Get properties
 891         mlt_properties properties = MLT_FRAME_PROPERTIES( this );
 892
 893         // Get the input image, width and height
 894         uint8_t *input = mlt_properties_get_data( properties, "image", NULL );
 895         uint8_t *alpha = mlt_frame_get_alpha_mask( this );
 896
 897         int iwidth = mlt_properties_get_int( properties, "width" );
 898         int iheight = mlt_properties_get_int( properties, "height" );
 899
 900         // If width and height are correct, don't do anything
 901         if ( iwidth != owidth || iheight != oheight )
 902         {
 903                 uint8_t alpha_value = mlt_properties_get_int( properties, "resize_alpha" );
 904
 905                 // Create the output image
 906                 uint8_t *output = mlt_pool_alloc( owidth * ( oheight + 1 ) * 2 );
 907
 908                 // Call the generic resize
 909                 mlt_resize_yuv422( output, owidth, oheight, input, iwidth, iheight );
 910
 911                 // Now update the frame
 912                 mlt_properties_set_data( properties, "image", output, owidth * ( oheight + 1 ) * 2, ( mlt_destructor )mlt_pool_release, NULL );
 913                 mlt_properties_set_int( properties, "width", owidth );
 914                 mlt_properties_set_int( properties, "height", oheight );
 915
 916                 // We should resize the alpha too
 917                 alpha = mlt_resize_alpha( alpha, owidth, oheight, iwidth, iheight, alpha_value );
 918                 if ( alpha != NULL )
 919                 {
 920                         mlt_properties_set_data( properties, "alpha", alpha, owidth * oheight, ( mlt_destructor )mlt_pool_release, NULL );
 921                         this->get_alpha_mask = NULL;
 922                 }
 923
 924                 // Return the output
 925                 return output;
 926         }
 927         // No change, return input
 928         return input;
 929 }
 930
 931 /** A rescaling function for yuv422 frames - low quality, and provided for testing
 932         only. It assumes yuv422 images available on the frame so use with care.
 933 */
 934
 935 uint8_t *mlt_frame_rescale_yuv422( mlt_frame this, int owidth, int oheight )
 936 {
 937         // Get properties
 938         mlt_properties properties = MLT_FRAME_PROPERTIES( this );
 939
 940         // Get the input image, width and height
 941         uint8_t *input = mlt_properties_get_data( properties, "image", NULL );
 942         int iwidth = mlt_properties_get_int( properties, "width" );
 943         int iheight = mlt_properties_get_int( properties, "height" );
 944
 945         // If width and height are correct, don't do anything
 946         if ( iwidth != owidth || iheight != oheight )
 947         {
 948                 // Create the output image
 949                 uint8_t *output = mlt_pool_alloc( owidth * ( oheight + 1 ) * 2 );
 950
 951                 // Calculate strides
 952                 int istride = iwidth * 2;
 953                 int ostride = owidth * 2;
 954
 955                 iwidth = iwidth - ( iwidth % 4 );
 956
 957                 // Derived coordinates
 958                 int dy, dx;
 959
 960         // Calculate ranges
 961         int out_x_range = owidth / 2;
 962         int out_y_range = oheight / 2;
 963         int in_x_range = iwidth / 2;
 964         int in_y_range = iheight / 2;
 965
 966         // Output pointers
 967         register uint8_t *out_line = output;
 968         register uint8_t *out_ptr;
 969
 970         // Calculate a middle pointer
 971         uint8_t *in_middle = input + istride * in_y_range + in_x_range * 2;
 972         uint8_t *in_line;
 973
 974                 // Generate the affine transform scaling values
 975                 register int scale_width = ( iwidth << 16 ) / owidth;
 976                 register int scale_height = ( iheight << 16 ) / oheight;
 977                 register int base = 0;
 978
 979                 int outer = out_x_range * scale_width;
 980                 int bottom = out_y_range * scale_height;
 981
 982         // Loop for the entirety of our output height.
 983         for ( dy = - bottom; dy < bottom; dy += scale_height )
 984         {
 985                 // Start at the beginning of the line
 986                 out_ptr = out_line;
 987
 988                 // Pointer to the middle of the input line
 989                 in_line = in_middle + ( dy >> 16 ) * istride;
 990
 991                 // Loop for the entirety of our output row.
 992                 for ( dx = - outer; dx < outer; dx += scale_width )
 993                 {
 994                                 base = dx >> 15;
 995                                 base &= 0xfffffffe;
 996                                 *out_ptr ++ = *( in_line + base );
 997                                 base &= 0xfffffffc;
 998                                 *out_ptr ++ = *( in_line + base + 1 );
 999                                 dx += scale_width;
1000                                 base = dx >> 15;
1001                                 base &= 0xfffffffe;
1002                                 *out_ptr ++ = *( in_line + base );
1003                                 base &= 0xfffffffc;
1004                                 *out_ptr ++ = *( in_line + base + 3 );
1005                 }
1006
1007                 // Move to next output line
1008                 out_line += ostride;
1009         }
1010
1011                 // Now update the frame
1012                 mlt_properties_set_data( properties, "image", output, owidth * ( oheight + 1 ) * 2, ( mlt_destructor )mlt_pool_release, NULL );
1013                 mlt_properties_set_int( properties, "width", owidth );
1014                 mlt_properties_set_int( properties, "height", oheight );
1015
1016                 // Return the output
1017                 return output;
1018         }
1019
1020         // No change, return input
1021         return input;
1022 }
1023
1024 int mlt_frame_mix_audio( mlt_frame this, mlt_frame that, float weight_start, float weight_end, int16_t **buffer, mlt_audio_format *format, int *frequency, int *channels, int *samples )
1025 {
1026         int ret = 0;
1027         int16_t *src, *dest;
1028         int frequency_src = *frequency, frequency_dest = *frequency;
1029         int channels_src = *channels, channels_dest = *channels;
1030         int samples_src = *samples, samples_dest = *samples;
1031         int i, j;
1032         double d = 0, s = 0;
1033
1034         mlt_frame_get_audio( that, &src, format, &frequency_src, &channels_src, &samples_src );
1035         mlt_frame_get_audio( this, &dest, format, &frequency_dest, &channels_dest, &samples_dest );
1036
1037         int silent = mlt_properties_get_int( MLT_FRAME_PROPERTIES( this ), "silent_audio" );
1038         mlt_properties_set_int( MLT_FRAME_PROPERTIES( this ), "silent_audio", 0 );
1039         if ( silent )
1040                 memset( dest, 0, samples_dest * channels_dest * sizeof( int16_t ) );
1041
1042         silent = mlt_properties_get_int( MLT_FRAME_PROPERTIES( that ), "silent_audio" );
1043         mlt_properties_set_int( MLT_FRAME_PROPERTIES( that ), "silent_audio", 0 );
1044         if ( silent )
1045                 memset( src, 0, samples_src * channels_src * sizeof( int16_t ) );
1046
1047         if ( channels_src > 6 )
1048                 channels_src = 0;
1049         if ( channels_dest > 6 )
1050                 channels_dest = 0;
1051         if ( samples_src > 4000 )
1052                 samples_src = 0;
1053         if ( samples_dest > 4000 )
1054                 samples_dest = 0;
1055
1056         // determine number of samples to process
1057         *samples = samples_src < samples_dest ? samples_src : samples_dest;
1058         *channels = channels_src < channels_dest ? channels_src : channels_dest;
1059         *buffer = dest;
1060         *frequency = frequency_dest;
1061
1062         // Compute a smooth ramp over start to end
1063         float weight = weight_start;
1064         float weight_step = ( weight_end - weight_start ) / *samples;
1065
1066         if ( src == dest )
1067         {
1068                 *samples = samples_src;
1069                 *channels = channels_src;
1070                 *buffer = src;
1071                 *frequency = frequency_src;
1072                 return ret;
1073         }
1074
1075         // Mixdown
1076         for ( i = 0; i < *samples; i++ )
1077         {
1078                 for ( j = 0; j < *channels; j++ )
1079                 {
1080                         if ( j < channels_dest )
1081                                 d = (double) dest[ i * channels_dest + j ];
1082                         if ( j < channels_src )
1083                                 s = (double) src[ i * channels_src + j ];
1084                         dest[ i * channels_dest + j ] = s * weight + d * ( 1.0 - weight );
1085                 }
1086                 weight += weight_step;
1087         }
1088
1089         return ret;
1090 }
1091
1092 // Replacement for broken mlt_frame_audio_mix - this filter uses an inline low pass filter
1093 // to allow mixing without volume hacking
1094 int mlt_frame_combine_audio( mlt_frame this, mlt_frame that, int16_t **buffer, mlt_audio_format *format, int *frequency, int *channels, int *samples )
1095 {
1096         int ret = 0;
1097         int16_t *src, *dest;
1098         int frequency_src = *frequency, frequency_dest = *frequency;
1099         int channels_src = *channels, channels_dest = *channels;
1100         int samples_src = *samples, samples_dest = *samples;
1101         int i, j;
1102         double vp[ 6 ];
1103         double b_weight = 1.0;
1104
1105         if ( mlt_properties_get_int( MLT_FRAME_PROPERTIES( this ), "meta.mixdown" ) )
1106                 b_weight = 1.0 - mlt_properties_get_double( MLT_FRAME_PROPERTIES( this ), "meta.volume" );
1107
1108         mlt_frame_get_audio( that, &src, format, &frequency_src, &channels_src, &samples_src );
1109         mlt_frame_get_audio( this, &dest, format, &frequency_dest, &channels_dest, &samples_dest );
1110
1111         int silent = mlt_properties_get_int( MLT_FRAME_PROPERTIES( this ), "silent_audio" );
1112         mlt_properties_set_int( MLT_FRAME_PROPERTIES( this ), "silent_audio", 0 );
1113         if ( silent )
1114                 memset( dest, 0, samples_dest * channels_dest * sizeof( int16_t ) );
1115
1116         silent = mlt_properties_get_int( MLT_FRAME_PROPERTIES( that ), "silent_audio" );
1117         mlt_properties_set_int( MLT_FRAME_PROPERTIES( that ), "silent_audio", 0 );
1118         if ( silent )
1119                 memset( src, 0, samples_src * channels_src * sizeof( int16_t ) );
1120
1121         if ( src == dest )
1122         {
1123                 *samples = samples_src;
1124                 *channels = channels_src;
1125                 *buffer = src;
1126                 *frequency = frequency_src;
1127                 return ret;
1128         }
1129
1130         // determine number of samples to process
1131         *samples = samples_src < samples_dest ? samples_src : samples_dest;
1132         *channels = channels_src < channels_dest ? channels_src : channels_dest;
1133         *buffer = dest;
1134         *frequency = frequency_dest;
1135
1136         for ( j = 0; j < *channels; j++ )
1137                 vp[ j ] = ( double )dest[ j ];
1138
1139         double Fc = 0.5;
1140         double B = exp(-2.0 * M_PI * Fc);
1141         double A = 1.0 - B;
1142         double v;
1143
1144         for ( i = 0; i < *samples; i++ )
1145         {
1146                 for ( j = 0; j < *channels; j++ )
1147                 {
1148                         v = ( double )( b_weight * dest[ i * channels_dest + j ] + src[ i * channels_src + j ] );
1149                         v = v < -32767 ? -32767 : v > 32768 ? 32768 : v;
1150                         vp[ j ] = dest[ i * channels_dest + j ] = ( int16_t )( v * A + vp[ j ] * B );
1151                 }
1152         }
1153
1154         return ret;
1155 }
1156
1157 /* Will this break when mlt_position is converted to double? -Zach */
1158 int mlt_sample_calculator( float fps, int frequency, int64_t position )
1159 {
1160         int samples = 0;
1161
1162         if ( ( int )( fps * 100 ) == 2997 )
1163         {
1164                 samples = frequency / 30;
1165
1166                 switch ( frequency )
1167                 {
1168                         case 48000:
1169                                 if ( position % 5 != 0 )
1170                                         samples += 2;
1171                                 break;
1172                         case 44100:
1173                                 if ( position % 300 == 0 )
1174                                         samples = 1471;
1175                                 else if ( position % 30 == 0 )
1176                                         samples = 1470;
1177                                 else if ( position % 2 == 0 )
1178                                         samples = 1472;
1179                                 else
1180                                         samples = 1471;
1181                                 break;
1182                         case 32000:
1183                                 if ( position % 30 == 0 )
1184                                         samples = 1068;
1185                                 else if ( position % 29 == 0 )
1186                                         samples = 1067;
1187                                 else if ( position % 4 == 2 )
1188                                         samples = 1067;
1189                                 else
1190                                         samples = 1068;
1191                                 break;
1192                         default:
1193                                 samples = 0;
1194                 }
1195         }
1196         else if ( fps != 0 )
1197         {
1198                 samples = frequency / fps;
1199         }
1200
1201         return samples;
1202 }
1203
1204 int64_t mlt_sample_calculator_to_now( float fps, int frequency, int64_t frame )
1205 {
1206         int64_t samples = 0;
1207
1208         // TODO: Correct rules for NTSC and drop the * 100 hack
1209         if ( ( int )( fps * 100 ) == 2997 )
1210         {
1211                 samples = ( ( double )( frame * frequency ) / 30 );
1212                 switch( frequency )
1213                 {
1214                         case 48000:
1215                                 samples += 2 * ( frame / 5 );
1216                                 break;
1217                         case 44100:
1218                                 samples += frame + ( frame / 2 ) - ( frame / 30 ) + ( frame / 300 );
1219                                 break;
1220                         case 32000:
1221                                 samples += ( 2 * frame ) - ( frame / 4 ) - ( frame / 29 );
1222                                 break;
1223                 }
1224         }
1225         else if ( fps != 0 )
1226         {
1227                 samples = ( ( frame * frequency ) / ( int )fps );
1228         }
1229
1230         return samples;
1231 }