mxfdec.c 90.7 KB
Newer Older
1 2
/*
 * MXF demuxer.
3
 * Copyright (c) 2006 SmartJog S.A., Baptiste Coudurier <baptiste dot coudurier at smartjog dot com>
4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48
 *
 * This file is part of FFmpeg.
 *
 * FFmpeg is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Lesser General Public
 * License as published by the Free Software Foundation; either
 * version 2.1 of the License, or (at your option) any later version.
 *
 * FFmpeg is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * Lesser General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General Public
 * License along with FFmpeg; if not, write to the Free Software
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
 */

/*
 * References
 * SMPTE 336M KLV Data Encoding Protocol Using Key-Length-Value
 * SMPTE 377M MXF File Format Specifications
 * SMPTE 378M Operational Pattern 1a
 * SMPTE 379M MXF Generic Container
 * SMPTE 381M Mapping MPEG Streams into the MXF Generic Container
 * SMPTE 382M Mapping AES3 and Broadcast Wave Audio into the MXF Generic Container
 * SMPTE 383M Mapping DV-DIF Data to the MXF Generic Container
 *
 * Principle
 * Search for Track numbers which will identify essence element KLV packets.
 * Search for SourcePackage which define tracks which contains Track numbers.
 * Material Package contains tracks with reference to SourcePackage tracks.
 * Search for Descriptors (Picture, Sound) which contains codec info and parameters.
 * Assign Descriptors to correct Tracks.
 *
 * Metadata reading functions read Local Tags, get InstanceUID(0x3C0A) then add MetaDataSet to MXFContext.
 * Metadata parsing resolves Strong References to objects.
 *
 * Simple demuxer, only OP1A supported and some files might not work at all.
 * Only tracks with associated descriptors will be decoded. "Highly Desirable" SMPTE 377M D.1
 */

//#define DEBUG

#include "libavutil/aes.h"
49
#include "libavutil/avassert.h"
50
#include "libavutil/mathematics.h"
51
#include "libavcodec/bytestream.h"
52
#include "libavutil/timecode.h"
53
#include "avformat.h"
54
#include "internal.h"
55
#include "mxf.h"
56

Tomas Härdin's avatar
Tomas Härdin committed
57 58 59 60 61 62 63
typedef enum {
    Header,
    BodyPartition,
    Footer
} MXFPartitionType;

typedef enum {
64
    OP1a = 1,
Tomas Härdin's avatar
Tomas Härdin committed
65 66 67 68 69 70 71 72 73
    OP1b,
    OP1c,
    OP2a,
    OP2b,
    OP2c,
    OP3a,
    OP3b,
    OP3c,
    OPAtom,
74
    OPSONYOpt,  /* FATE sample, violates the spec in places */
Tomas Härdin's avatar
Tomas Härdin committed
75 76 77 78 79 80 81 82 83
} MXFOP;

typedef struct {
    int closed;
    int complete;
    MXFPartitionType type;
    uint64_t previous_partition;
    int index_sid;
    int body_sid;
84
    int64_t this_partition;
85
    int64_t essence_offset;         ///< absolute offset of essence
86
    int64_t essence_length;
87 88 89 90
    int32_t kag_size;
    int64_t header_byte_count;
    int64_t index_byte_count;
    int pack_length;
Tomas Härdin's avatar
Tomas Härdin committed
91 92
} MXFPartition;

93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117
typedef struct {
    UID uid;
    enum MXFMetadataSetType type;
    UID source_container_ul;
} MXFCryptoContext;

typedef struct {
    UID uid;
    enum MXFMetadataSetType type;
    UID source_package_uid;
    UID data_definition_ul;
    int64_t duration;
    int64_t start_position;
    int source_track_id;
} MXFStructuralComponent;

typedef struct {
    UID uid;
    enum MXFMetadataSetType type;
    UID data_definition_ul;
    UID *structural_components_refs;
    int structural_components_count;
    int64_t duration;
} MXFSequence;

118 119 120 121 122 123 124 125 126
typedef struct {
    UID uid;
    enum MXFMetadataSetType type;
    int drop_frame;
    int start_frame;
    struct AVRational rate;
    AVTimecode tc;
} MXFTimecodeComponent;

127 128 129 130 131 132 133 134
typedef struct {
    UID uid;
    enum MXFMetadataSetType type;
    MXFSequence *sequence; /* mandatory, and only one */
    UID sequence_ref;
    int track_id;
    uint8_t track_number[4];
    AVRational edit_rate;
135
    int intra_only;
136
    uint64_t sample_count;
137
    int64_t original_duration;  ///< duration before multiplying st->duration by SampleRate/EditRate
138 139 140 141 142 143 144 145 146 147
} MXFTrack;

typedef struct {
    UID uid;
    enum MXFMetadataSetType type;
    UID essence_container_ul;
    UID essence_codec_ul;
    AVRational sample_rate;
    AVRational aspect_ratio;
    int width;
148 149
    int height; /* Field height, not frame height */
    int frame_layout; /* See MXFFrameLayout enum */
150 151
    int channels;
    int bits_per_sample;
152 153 154
    unsigned int component_depth;
    unsigned int horiz_subsampling;
    unsigned int vert_subsampling;
155 156 157 158 159
    UID *sub_descriptors_refs;
    int sub_descriptors_count;
    int linked_track_id;
    uint8_t *extradata;
    int extradata_size;
160
    enum AVPixelFormat pix_fmt;
161 162
} MXFDescriptor;

163 164 165
typedef struct {
    UID uid;
    enum MXFMetadataSetType type;
166 167 168 169 170 171
    int edit_unit_byte_count;
    int index_sid;
    int body_sid;
    AVRational index_edit_rate;
    uint64_t index_start_position;
    uint64_t index_duration;
172
    int8_t *temporal_offset_entries;
173 174 175
    int *flag_entries;
    uint64_t *stream_offset_entries;
    int nb_index_entries;
176 177
} MXFIndexTableSegment;

178 179 180 181 182 183 184 185 186 187 188 189 190 191 192
typedef struct {
    UID uid;
    enum MXFMetadataSetType type;
    UID package_uid;
    UID *tracks_refs;
    int tracks_count;
    MXFDescriptor *descriptor; /* only one */
    UID descriptor_ref;
} MXFPackage;

typedef struct {
    UID uid;
    enum MXFMetadataSetType type;
} MXFMetadataSet;

193 194 195 196 197 198 199 200 201
/* decoded index table */
typedef struct {
    int index_sid;
    int body_sid;
    int nb_ptses;               /* number of PTSes or total duration of index */
    int64_t first_dts;          /* DTS = EditUnit + first_dts */
    int64_t *ptses;             /* maps EditUnit -> PTS */
    int nb_segments;
    MXFIndexTableSegment **segments;    /* sorted by IndexStartPosition */
202
    AVIndexEntry *fake_index;   /* used for calling ff_index_search_timestamp() */
203 204
} MXFIndexTable;

205
typedef struct {
Tomas Härdin's avatar
Tomas Härdin committed
206 207 208
    MXFPartition *partitions;
    unsigned partitions_count;
    MXFOP op;
209 210 211 212 213 214 215 216
    UID *packages_refs;
    int packages_count;
    MXFMetadataSet **metadata_sets;
    int metadata_sets_count;
    AVFormatContext *fc;
    struct AVAES *aesc;
    uint8_t *local_tags;
    int local_tags_count;
217
    uint64_t footer_partition;
218 219
    KLVPacket current_klv_data;
    int current_klv_index;
220 221 222 223
    int run_in;
    MXFPartition *current_partition;
    int parsing_backward;
    int64_t last_forward_tell;
224
    int last_forward_partition;
225
    int current_edit_unit;
226 227
    int nb_index_tables;
    MXFIndexTable *index_tables;
228
    int edit_units_per_packet;      ///< how many edit units to read at a time (PCM, OPAtom)
229 230 231 232 233 234 235
} MXFContext;

enum MXFWrappingScheme {
    Frame,
    Clip,
};

236 237
/* NOTE: klv_offset is not set (-1) for local keys */
typedef int MXFMetadataReadFunc(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset);
238

239 240
typedef struct {
    const UID key;
241
    MXFMetadataReadFunc *read;
242 243 244 245
    int ctx_size;
    enum MXFMetadataSetType type;
} MXFMetadataReadTableEntry;

246 247
static int mxf_read_close(AVFormatContext *s);

248 249 250
/* partial keys to match */
static const uint8_t mxf_header_partition_pack_key[]       = { 0x06,0x0e,0x2b,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x02 };
static const uint8_t mxf_essence_element_key[]             = { 0x06,0x0e,0x2b,0x34,0x01,0x02,0x01,0x01,0x0d,0x01,0x03,0x01 };
251
static const uint8_t mxf_avid_essence_element_key[]        = { 0x06,0x0e,0x2b,0x34,0x01,0x02,0x01,0x01,0x0e,0x04,0x03,0x01 };
252
static const uint8_t mxf_system_item_key[]                 = { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0D,0x01,0x03,0x01,0x04 };
253 254 255 256 257 258 259 260 261
static const uint8_t mxf_klv_key[]                         = { 0x06,0x0e,0x2b,0x34 };
/* complete keys to match */
static const uint8_t mxf_crypto_source_container_ul[]      = { 0x06,0x0e,0x2b,0x34,0x01,0x01,0x01,0x09,0x06,0x01,0x01,0x02,0x02,0x00,0x00,0x00 };
static const uint8_t mxf_encrypted_triplet_key[]           = { 0x06,0x0e,0x2b,0x34,0x02,0x04,0x01,0x07,0x0d,0x01,0x03,0x01,0x02,0x7e,0x01,0x00 };
static const uint8_t mxf_encrypted_essence_container[]     = { 0x06,0x0e,0x2b,0x34,0x04,0x01,0x01,0x07,0x0d,0x01,0x03,0x01,0x02,0x0b,0x01,0x00 };
static const uint8_t mxf_sony_mpeg4_extradata[]            = { 0x06,0x0e,0x2b,0x34,0x04,0x01,0x01,0x01,0x0e,0x06,0x06,0x02,0x02,0x01,0x00,0x00 };

#define IS_KLV_KEY(x, y) (!memcmp(x, y, sizeof(y)))

262
static int64_t klv_decode_ber_length(AVIOContext *pb)
263
{
264
    uint64_t size = avio_r8(pb);
265 266 267 268
    if (size & 0x80) { /* long form */
        int bytes_num = size & 0x7f;
        /* SMPTE 379M 5.3.4 guarantee that bytes_num must not exceed 8 bytes */
        if (bytes_num > 8)
269
            return AVERROR_INVALIDDATA;
270 271
        size = 0;
        while (bytes_num--)
272
            size = size << 8 | avio_r8(pb);
273 274 275 276
    }
    return size;
}

277
static int mxf_read_sync(AVIOContext *pb, const uint8_t *key, unsigned size)
278 279 280
{
    int i, b;
    for (i = 0; i < size && !url_feof(pb); i++) {
281
        b = avio_r8(pb);
282 283 284 285 286 287 288 289
        if (b == key[0])
            i = 0;
        else if (b != key[i])
            i = -1;
    }
    return i == size;
}

290
static int klv_read_packet(KLVPacket *klv, AVIOContext *pb)
291 292
{
    if (!mxf_read_sync(pb, mxf_klv_key, 4))
293
        return AVERROR_INVALIDDATA;
294
    klv->offset = avio_tell(pb) - 4;
295
    memcpy(klv->key, mxf_klv_key, 4);
296
    avio_read(pb, klv->key + 4, 12);
297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315
    klv->length = klv_decode_ber_length(pb);
    return klv->length == -1 ? -1 : 0;
}

static int mxf_get_stream_index(AVFormatContext *s, KLVPacket *klv)
{
    int i;

    for (i = 0; i < s->nb_streams; i++) {
        MXFTrack *track = s->streams[i]->priv_data;
        /* SMPTE 379M 7.3 */
        if (!memcmp(klv->key + sizeof(mxf_essence_element_key), track->track_number, sizeof(track->track_number)))
            return i;
    }
    /* return 0 if only one stream, for OP Atom files with 0 as track number */
    return s->nb_streams == 1 ? 0 : -1;
}

/* XXX: use AVBitStreamFilter */
316
static int mxf_get_d10_aes3_packet(AVIOContext *pb, AVStream *st, AVPacket *pkt, int64_t length)
317 318 319 320 321 322
{
    const uint8_t *buf_ptr, *end_ptr;
    uint8_t *data_ptr;
    int i;

    if (length > 61444) /* worst case PAL 1920 samples 8 channels */
323
        return AVERROR_INVALIDDATA;
324 325 326
    length = av_get_packet(pb, pkt, length);
    if (length < 0)
        return length;
327
    data_ptr = pkt->data;
328 329
    end_ptr = pkt->data + length;
    buf_ptr = pkt->data + 4; /* skip SMPTE 331M header */
330
    for (; buf_ptr + st->codec->channels*4 <= end_ptr; ) {
331 332
        for (i = 0; i < st->codec->channels; i++) {
            uint32_t sample = bytestream_get_le32(&buf_ptr);
333
            if (st->codec->bits_per_coded_sample == 24)
334 335 336 337 338 339
                bytestream_put_le24(&data_ptr, (sample >> 4) & 0xffffff);
            else
                bytestream_put_le16(&data_ptr, (sample >> 12) & 0xffff);
        }
        buf_ptr += 32 - st->codec->channels*4; // always 8 channels stored SMPTE 331M
    }
340
    av_shrink_packet(pkt, data_ptr - pkt->data);
341 342 343 344 345 346 347
    return 0;
}

static int mxf_decrypt_triplet(AVFormatContext *s, AVPacket *pkt, KLVPacket *klv)
{
    static const uint8_t checkv[16] = {0x43, 0x48, 0x55, 0x4b, 0x43, 0x48, 0x55, 0x4b, 0x43, 0x48, 0x55, 0x4b, 0x43, 0x48, 0x55, 0x4b};
    MXFContext *mxf = s->priv_data;
348
    AVIOContext *pb = s->pb;
349
    int64_t end = avio_tell(pb) + klv->length;
350
    int64_t size;
351 352 353 354 355 356 357
    uint64_t orig_size;
    uint64_t plaintext_size;
    uint8_t ivec[16];
    uint8_t tmpbuf[16];
    int index;

    if (!mxf->aesc && s->key && s->keylen == 16) {
358
        mxf->aesc = av_aes_alloc();
359
        if (!mxf->aesc)
360
            return AVERROR(ENOMEM);
361 362 363
        av_aes_init(mxf->aesc, s->key, 128, 1);
    }
    // crypto context
364
    avio_skip(pb, klv_decode_ber_length(pb));
365 366
    // plaintext offset
    klv_decode_ber_length(pb);
367
    plaintext_size = avio_rb64(pb);
368 369
    // source klv key
    klv_decode_ber_length(pb);
370
    avio_read(pb, klv->key, 16);
371
    if (!IS_KLV_KEY(klv, mxf_essence_element_key))
372
        return AVERROR_INVALIDDATA;
373 374
    index = mxf_get_stream_index(s, klv);
    if (index < 0)
375
        return AVERROR_INVALIDDATA;
376 377
    // source size
    klv_decode_ber_length(pb);
378
    orig_size = avio_rb64(pb);
379
    if (orig_size < plaintext_size)
380
        return AVERROR_INVALIDDATA;
381 382 383
    // enc. code
    size = klv_decode_ber_length(pb);
    if (size < 32 || size - 32 < orig_size)
384
        return AVERROR_INVALIDDATA;
385 386
    avio_read(pb, ivec, 16);
    avio_read(pb, tmpbuf, 16);
387 388 389 390 391
    if (mxf->aesc)
        av_aes_crypt(mxf->aesc, tmpbuf, tmpbuf, 1, ivec, 1);
    if (memcmp(tmpbuf, checkv, 16))
        av_log(s, AV_LOG_ERROR, "probably incorrect decryption key\n");
    size -= 32;
392 393 394 395 396
    size = av_get_packet(pb, pkt, size);
    if (size < 0)
        return size;
    else if (size < plaintext_size)
        return AVERROR_INVALIDDATA;
397 398 399 400
    size -= plaintext_size;
    if (mxf->aesc)
        av_aes_crypt(mxf->aesc, &pkt->data[plaintext_size],
                     &pkt->data[plaintext_size], size >> 4, ivec, 1);
401
    av_shrink_packet(pkt, orig_size);
402
    pkt->stream_index = index;
403
    avio_skip(pb, end - avio_tell(pb));
404 405 406
    return 0;
}

407
static int mxf_read_primer_pack(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
408
{
409
    MXFContext *mxf = arg;
410 411
    int item_num = avio_rb32(pb);
    int item_len = avio_rb32(pb);
412 413

    if (item_len != 18) {
414 415 416
        av_log_ask_for_sample(pb, "unsupported primer pack item length %d\n",
                              item_len);
        return AVERROR_PATCHWELCOME;
417
    }
418 419
    if (item_num > 65536) {
        av_log(mxf->fc, AV_LOG_ERROR, "item_num %d is too large\n", item_num);
420
        return AVERROR_INVALIDDATA;
421 422
    }
    mxf->local_tags = av_calloc(item_num, item_len);
423
    if (!mxf->local_tags)
424
        return AVERROR(ENOMEM);
425
    mxf->local_tags_count = item_num;
426
    avio_read(pb, mxf->local_tags, item_num*item_len);
427 428 429
    return 0;
}

430
static int mxf_read_partition_pack(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
Tomas Härdin's avatar
Tomas Härdin committed
431 432
{
    MXFContext *mxf = arg;
433
    MXFPartition *partition, *tmp_part;
Tomas Härdin's avatar
Tomas Härdin committed
434
    UID op;
435
    uint64_t footer_partition;
436
    uint32_t nb_essence_containers;
Tomas Härdin's avatar
Tomas Härdin committed
437 438 439 440

    if (mxf->partitions_count+1 >= UINT_MAX / sizeof(*mxf->partitions))
        return AVERROR(ENOMEM);

441 442
    tmp_part = av_realloc(mxf->partitions, (mxf->partitions_count + 1) * sizeof(*mxf->partitions));
    if (!tmp_part)
Tomas Härdin's avatar
Tomas Härdin committed
443
        return AVERROR(ENOMEM);
444
    mxf->partitions = tmp_part;
Tomas Härdin's avatar
Tomas Härdin committed
445

446 447 448 449 450 451 452 453 454 455 456 457 458 459
    if (mxf->parsing_backward) {
        /* insert the new partition pack in the middle
         * this makes the entries in mxf->partitions sorted by offset */
        memmove(&mxf->partitions[mxf->last_forward_partition+1],
                &mxf->partitions[mxf->last_forward_partition],
                (mxf->partitions_count - mxf->last_forward_partition)*sizeof(*mxf->partitions));
        partition = mxf->current_partition = &mxf->partitions[mxf->last_forward_partition];
    } else {
        mxf->last_forward_partition++;
        partition = mxf->current_partition = &mxf->partitions[mxf->partitions_count];
    }

    memset(partition, 0, sizeof(*partition));
    mxf->partitions_count++;
460
    partition->pack_length = avio_tell(pb) - klv_offset + size;
Tomas Härdin's avatar
Tomas Härdin committed
461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479

    switch(uid[13]) {
    case 2:
        partition->type = Header;
        break;
    case 3:
        partition->type = BodyPartition;
        break;
    case 4:
        partition->type = Footer;
        break;
    default:
        av_log(mxf->fc, AV_LOG_ERROR, "unknown partition type %i\n", uid[13]);
        return AVERROR_INVALIDDATA;
    }

    /* consider both footers to be closed (there is only Footer and CompleteFooter) */
    partition->closed = partition->type == Footer || !(uid[14] & 1);
    partition->complete = uid[14] > 2;
480 481
    avio_skip(pb, 4);
    partition->kag_size = avio_rb32(pb);
482
    partition->this_partition = avio_rb64(pb);
Tomas Härdin's avatar
Tomas Härdin committed
483
    partition->previous_partition = avio_rb64(pb);
484
    footer_partition = avio_rb64(pb);
485 486
    partition->header_byte_count = avio_rb64(pb);
    partition->index_byte_count = avio_rb64(pb);
Tomas Härdin's avatar
Tomas Härdin committed
487 488 489 490
    partition->index_sid = avio_rb32(pb);
    avio_skip(pb, 8);
    partition->body_sid = avio_rb32(pb);
    avio_read(pb, op, sizeof(UID));
491
    nb_essence_containers = avio_rb32(pb);
Tomas Härdin's avatar
Tomas Härdin committed
492

493 494 495
    /* some files don'thave FooterPartition set in every partition */
    if (footer_partition) {
        if (mxf->footer_partition && mxf->footer_partition != footer_partition) {
496 497
            av_log(mxf->fc, AV_LOG_ERROR,
                   "inconsistent FooterPartition value: %"PRIu64" != %"PRIu64"\n",
498 499 500 501 502 503
                   mxf->footer_partition, footer_partition);
        } else {
            mxf->footer_partition = footer_partition;
        }
    }

504 505 506 507
    av_dlog(mxf->fc,
            "PartitionPack: ThisPartition = 0x%"PRIX64
            ", PreviousPartition = 0x%"PRIX64", "
            "FooterPartition = 0x%"PRIX64", IndexSID = %i, BodySID = %i\n",
508
            partition->this_partition,
509
            partition->previous_partition, footer_partition,
Tomas Härdin's avatar
Tomas Härdin committed
510 511
            partition->index_sid, partition->body_sid);

512 513 514
    /* sanity check PreviousPartition if set */
    if (partition->previous_partition &&
        mxf->run_in + partition->previous_partition >= klv_offset) {
515 516
        av_log(mxf->fc, AV_LOG_ERROR,
               "PreviousPartition points to this partition or forward\n");
517 518 519
        return AVERROR_INVALIDDATA;
    }

Tomas Härdin's avatar
Tomas Härdin committed
520 521 522 523 524 525 526 527 528
    if      (op[12] == 1 && op[13] == 1) mxf->op = OP1a;
    else if (op[12] == 1 && op[13] == 2) mxf->op = OP1b;
    else if (op[12] == 1 && op[13] == 3) mxf->op = OP1c;
    else if (op[12] == 2 && op[13] == 1) mxf->op = OP2a;
    else if (op[12] == 2 && op[13] == 2) mxf->op = OP2b;
    else if (op[12] == 2 && op[13] == 3) mxf->op = OP2c;
    else if (op[12] == 3 && op[13] == 1) mxf->op = OP3a;
    else if (op[12] == 3 && op[13] == 2) mxf->op = OP3b;
    else if (op[12] == 3 && op[13] == 3) mxf->op = OP3c;
529
    else if (op[12] == 64&& op[13] == 1) mxf->op = OPSONYOpt;
530 531
    else if (op[12] == 0x10) {
        /* SMPTE 390m: "There shall be exactly one essence container"
532 533 534
         * The following block deals with files that violate this, namely:
         * 2011_DCPTEST_24FPS.V.mxf - two ECs, OP1a
         * abcdefghiv016f56415e.mxf - zero ECs, OPAtom, output by Avid AirSpeed */
535
        if (nb_essence_containers != 1) {
536 537
            MXFOP op = nb_essence_containers ? OP1a : OPAtom;

538 539
            /* only nag once */
            if (!mxf->op)
540 541
                av_log(mxf->fc, AV_LOG_WARNING, "\"OPAtom\" with %u ECs - assuming %s\n",
                       nb_essence_containers, op == OP1a ? "OP1a" : "OPAtom");
542

543
            mxf->op = op;
544 545 546
        } else
            mxf->op = OPAtom;
    } else {
547 548 549 550 551 552 553 554 555 556 557 558 559 560
        av_log(mxf->fc, AV_LOG_ERROR, "unknown operational pattern: %02xh %02xh - guessing OP1a\n", op[12], op[13]);
        mxf->op = OP1a;
    }

    if (partition->kag_size <= 0 || partition->kag_size > (1 << 20)) {
        av_log(mxf->fc, AV_LOG_WARNING, "invalid KAGSize %i - guessing ", partition->kag_size);

        if (mxf->op == OPSONYOpt)
            partition->kag_size = 512;
        else
            partition->kag_size = 1;

        av_log(mxf->fc, AV_LOG_WARNING, "%i\n", partition->kag_size);
    }
Tomas Härdin's avatar
Tomas Härdin committed
561 562 563 564

    return 0;
}

565 566
static int mxf_add_metadata_set(MXFContext *mxf, void *metadata_set)
{
567
    MXFMetadataSet **tmp;
568 569
    if (mxf->metadata_sets_count+1 >= UINT_MAX / sizeof(*mxf->metadata_sets))
        return AVERROR(ENOMEM);
570 571
    tmp = av_realloc(mxf->metadata_sets, (mxf->metadata_sets_count + 1) * sizeof(*mxf->metadata_sets));
    if (!tmp)
572
        return AVERROR(ENOMEM);
573
    mxf->metadata_sets = tmp;
574 575 576 577 578
    mxf->metadata_sets[mxf->metadata_sets_count] = metadata_set;
    mxf->metadata_sets_count++;
    return 0;
}

579
static int mxf_read_cryptographic_context(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
580
{
581
    MXFCryptoContext *cryptocontext = arg;
582
    if (size != 16)
583
        return AVERROR_INVALIDDATA;
584
    if (IS_KLV_KEY(uid, mxf_crypto_source_container_ul))
585
        avio_read(pb, cryptocontext->source_container_ul, 16);
586 587 588
    return 0;
}

589
static int mxf_read_content_storage(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
590
{
591
    MXFContext *mxf = arg;
592 593
    switch (tag) {
    case 0x1901:
594
        mxf->packages_count = avio_rb32(pb);
595
        mxf->packages_refs = av_calloc(mxf->packages_count, sizeof(UID));
596
        if (!mxf->packages_refs)
597
            return AVERROR(ENOMEM);
598
        avio_skip(pb, 4); /* useless size of objects, always 16 according to specs */
599
        avio_read(pb, (uint8_t *)mxf->packages_refs, mxf->packages_count * sizeof(UID));
600 601 602 603 604
        break;
    }
    return 0;
}

605
static int mxf_read_source_clip(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
606
{
607
    MXFStructuralComponent *source_clip = arg;
608 609
    switch(tag) {
    case 0x0202:
610
        source_clip->duration = avio_rb64(pb);
611 612
        break;
    case 0x1201:
613
        source_clip->start_position = avio_rb64(pb);
614 615 616
        break;
    case 0x1101:
        /* UMID, only get last 16 bytes */
617
        avio_skip(pb, 16);
618
        avio_read(pb, source_clip->source_package_uid, 16);
619 620
        break;
    case 0x1102:
621
        source_clip->source_track_id = avio_rb32(pb);
622 623 624 625 626
        break;
    }
    return 0;
}

627
static int mxf_read_material_package(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
628
{
629
    MXFPackage *package = arg;
630 631
    switch(tag) {
    case 0x4403:
632
        package->tracks_count = avio_rb32(pb);
633
        package->tracks_refs = av_calloc(package->tracks_count, sizeof(UID));
634
        if (!package->tracks_refs)
635
            return AVERROR(ENOMEM);
636
        avio_skip(pb, 4); /* useless size of objects, always 16 according to specs */
637
        avio_read(pb, (uint8_t *)package->tracks_refs, package->tracks_count * sizeof(UID));
638 639 640 641 642
        break;
    }
    return 0;
}

643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659
static int mxf_read_timecode_component(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
{
    MXFTimecodeComponent *mxf_timecode = arg;
    switch(tag) {
    case 0x1501:
        mxf_timecode->start_frame = avio_rb64(pb);
        break;
    case 0x1502:
        mxf_timecode->rate = (AVRational){avio_rb16(pb), 1};
        break;
    case 0x1503:
        mxf_timecode->drop_frame = avio_r8(pb);
        break;
    }
    return 0;
}

660
static int mxf_read_track(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
661
{
662
    MXFTrack *track = arg;
663 664
    switch(tag) {
    case 0x4801:
665
        track->track_id = avio_rb32(pb);
666 667
        break;
    case 0x4804:
668
        avio_read(pb, track->track_number, 4);
669 670
        break;
    case 0x4B01:
671
        track->edit_rate.num = avio_rb32(pb);
672
        track->edit_rate.den = avio_rb32(pb);
673 674
        break;
    case 0x4803:
675
        avio_read(pb, track->sequence_ref, 16);
676 677 678 679 680
        break;
    }
    return 0;
}

681
static int mxf_read_sequence(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
682
{
683
    MXFSequence *sequence = arg;
684 685
    switch(tag) {
    case 0x0202:
686
        sequence->duration = avio_rb64(pb);
687 688
        break;
    case 0x0201:
689
        avio_read(pb, sequence->data_definition_ul, 16);
690 691
        break;
    case 0x1001:
692
        sequence->structural_components_count = avio_rb32(pb);
693
        sequence->structural_components_refs = av_calloc(sequence->structural_components_count, sizeof(UID));
694
        if (!sequence->structural_components_refs)
695
            return AVERROR(ENOMEM);
696
        avio_skip(pb, 4); /* useless size of objects, always 16 according to specs */
697
        avio_read(pb, (uint8_t *)sequence->structural_components_refs, sequence->structural_components_count * sizeof(UID));
698 699 700 701 702
        break;
    }
    return 0;
}

703
static int mxf_read_source_package(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
704
{
705
    MXFPackage *package = arg;
706 707
    switch(tag) {
    case 0x4403:
708
        package->tracks_count = avio_rb32(pb);
709
        package->tracks_refs = av_calloc(package->tracks_count, sizeof(UID));
710
        if (!package->tracks_refs)
711
            return AVERROR(ENOMEM);
712
        avio_skip(pb, 4); /* useless size of objects, always 16 according to specs */
713
        avio_read(pb, (uint8_t *)package->tracks_refs, package->tracks_count * sizeof(UID));
714 715 716
        break;
    case 0x4401:
        /* UMID, only get last 16 bytes */
717
        avio_skip(pb, 16);
718
        avio_read(pb, package->package_uid, 16);
719 720
        break;
    case 0x4701:
721
        avio_read(pb, package->descriptor_ref, 16);
722 723 724 725 726
        break;
    }
    return 0;
}

727 728
static int mxf_read_index_entry_array(AVIOContext *pb, MXFIndexTableSegment *segment)
{
729
    int i, length;
730 731

    segment->nb_index_entries = avio_rb32(pb);
732

733 734
    length = avio_rb32(pb);

735 736
    if (!(segment->temporal_offset_entries=av_calloc(segment->nb_index_entries, sizeof(*segment->temporal_offset_entries))) ||
        !(segment->flag_entries          = av_calloc(segment->nb_index_entries, sizeof(*segment->flag_entries))) ||
737 738 739 740
        !(segment->stream_offset_entries = av_calloc(segment->nb_index_entries, sizeof(*segment->stream_offset_entries))))
        return AVERROR(ENOMEM);

    for (i = 0; i < segment->nb_index_entries; i++) {
741 742
        segment->temporal_offset_entries[i] = avio_r8(pb);
        avio_r8(pb);                                        /* KeyFrameOffset */
743 744
        segment->flag_entries[i] = avio_r8(pb);
        segment->stream_offset_entries[i] = avio_rb64(pb);
745
        avio_skip(pb, length - 11);
746 747 748 749
    }
    return 0;
}

750
static int mxf_read_index_table_segment(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
751
{
752
    MXFIndexTableSegment *segment = arg;
753
    switch(tag) {
754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782
    case 0x3F05:
        segment->edit_unit_byte_count = avio_rb32(pb);
        av_dlog(NULL, "EditUnitByteCount %d\n", segment->edit_unit_byte_count);
        break;
    case 0x3F06:
        segment->index_sid = avio_rb32(pb);
        av_dlog(NULL, "IndexSID %d\n", segment->index_sid);
        break;
    case 0x3F07:
        segment->body_sid = avio_rb32(pb);
        av_dlog(NULL, "BodySID %d\n", segment->body_sid);
        break;
    case 0x3F0A:
        av_dlog(NULL, "IndexEntryArray found\n");
        return mxf_read_index_entry_array(pb, segment);
    case 0x3F0B:
        segment->index_edit_rate.num = avio_rb32(pb);
        segment->index_edit_rate.den = avio_rb32(pb);
        av_dlog(NULL, "IndexEditRate %d/%d\n", segment->index_edit_rate.num,
                segment->index_edit_rate.den);
        break;
    case 0x3F0C:
        segment->index_start_position = avio_rb64(pb);
        av_dlog(NULL, "IndexStartPosition %"PRId64"\n", segment->index_start_position);
        break;
    case 0x3F0D:
        segment->index_duration = avio_rb64(pb);
        av_dlog(NULL, "IndexDuration %"PRId64"\n", segment->index_duration);
        break;
783 784 785 786
    }
    return 0;
}

787
static void mxf_read_pixel_layout(AVIOContext *pb, MXFDescriptor *descriptor)
788
{
789
    int code, value, ofs = 0;
Tomas Härdin's avatar
Tomas Härdin committed
790
    char layout[16] = {0}; /* not for printing, may end up not terminated on purpose */
791 792

    do {
793 794
        code = avio_r8(pb);
        value = avio_r8(pb);
795
        av_dlog(NULL, "pixel layout: code %#x\n", code);
796

Tomas Härdin's avatar
Tomas Härdin committed
797
        if (ofs <= 14) {
798 799
            layout[ofs++] = code;
            layout[ofs++] = value;
800 801
        } else
            break;  /* don't read byte by byte on sneaky files filled with lots of non-zeroes */
802
    } while (code != 0); /* SMPTE 377M E.2.46 */
803 804

    ff_mxf_decode_pixel_layout(layout, &descriptor->pix_fmt);
805 806
}

807
static int mxf_read_generic_descriptor(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
808
{
809
    MXFDescriptor *descriptor = arg;
810
    descriptor->pix_fmt = AV_PIX_FMT_NONE;
811 812
    switch(tag) {
    case 0x3F01:
813
        descriptor->sub_descriptors_count = avio_rb32(pb);
814
        descriptor->sub_descriptors_refs = av_calloc(descriptor->sub_descriptors_count, sizeof(UID));
815
        if (!descriptor->sub_descriptors_refs)
816
            return AVERROR(ENOMEM);
817
        avio_skip(pb, 4); /* useless size of objects, always 16 according to specs */
818
        avio_read(pb, (uint8_t *)descriptor->sub_descriptors_refs, descriptor->sub_descriptors_count * sizeof(UID));
819 820
        break;
    case 0x3004:
821
        avio_read(pb, descriptor->essence_container_ul, 16);
822 823
        break;
    case 0x3006:
824
        descriptor->linked_track_id = avio_rb32(pb);
825 826
        break;
    case 0x3201: /* PictureEssenceCoding */
827
        avio_read(pb, descriptor->essence_codec_ul, 16);
828 829
        break;
    case 0x3203:
830
        descriptor->width = avio_rb32(pb);
831 832
        break;
    case 0x3202:
833
        descriptor->height = avio_rb32(pb);
834
        break;
835 836 837
    case 0x320C:
        descriptor->frame_layout = avio_r8(pb);
        break;
838
    case 0x320E:
839 840
        descriptor->aspect_ratio.num = avio_rb32(pb);
        descriptor->aspect_ratio.den = avio_rb32(pb);
841
        break;
842 843 844 845 846 847 848 849 850
    case 0x3301:
        descriptor->component_depth = avio_rb32(pb);
        break;
    case 0x3302:
        descriptor->horiz_subsampling = avio_rb32(pb);
        break;
    case 0x3308:
        descriptor->vert_subsampling = avio_rb32(pb);
        break;
851
    case 0x3D03:
852 853
        descriptor->sample_rate.num = avio_rb32(pb);
        descriptor->sample_rate.den = avio_rb32(pb);
854 855
        break;
    case 0x3D06: /* SoundEssenceCompression */
856
        avio_read(pb, descriptor->essence_codec_ul, 16);
857 858
        break;
    case 0x3D07:
859
        descriptor->channels = avio_rb32(pb);
860 861
        break;
    case 0x3D01:
862
        descriptor->bits_per_sample = avio_rb32(pb);
863 864 865 866 867 868 869
        break;
    case 0x3401:
        mxf_read_pixel_layout(pb, descriptor);
        break;
    default:
        /* Private uid used by SONY C0023S01.mxf */
        if (IS_KLV_KEY(uid, mxf_sony_mpeg4_extradata)) {
870 871 872 873
            if (descriptor->extradata)
                av_log(NULL, AV_LOG_WARNING, "Duplicate sony_mpeg4_extradata\n");
            av_free(descriptor->extradata);
            descriptor->extradata_size = 0;
874
            descriptor->extradata = av_malloc(size);
875
            if (!descriptor->extradata)
876
                return AVERROR(ENOMEM);
877
            descriptor->extradata_size = size;
878
            avio_read(pb, descriptor->extradata, size);
879 880 881 882 883 884 885 886 887 888 889 890 891 892 893 894 895 896 897 898 899 900
        }
        break;
    }
    return 0;
}

/*
 * Match an uid independently of the version byte and up to len common bytes
 * Returns: boolean
 */
static int mxf_match_uid(const UID key, const UID uid, int len)
{
    int i;
    for (i = 0; i < len; i++) {
        if (i != 7 && key[i] != uid[i])
            return 0;
    }
    return 1;
}

static const MXFCodecUL *mxf_get_codec_ul(const MXFCodecUL *uls, UID *uid)
{
901
    while (uls->uid[0]) {
902 903 904 905 906 907 908 909 910 911 912 913 914 915 916 917 918 919 920 921 922 923
        if(mxf_match_uid(uls->uid, *uid, uls->matching_len))
            break;
        uls++;
    }
    return uls;
}

static void *mxf_resolve_strong_ref(MXFContext *mxf, UID *strong_ref, enum MXFMetadataSetType type)
{
    int i;

    if (!strong_ref)
        return NULL;
    for (i = 0; i < mxf->metadata_sets_count; i++) {
        if (!memcmp(*strong_ref, mxf->metadata_sets[i]->uid, 16) &&
            (type == AnyType || mxf->metadata_sets[i]->type == type)) {
            return mxf->metadata_sets[i];
        }
    }
    return NULL;
}

924
static const MXFCodecUL mxf_picture_essence_container_uls[] = {
925
    // video essence container uls
926 927 928 929
    { { 0x06,0x0E,0x2B,0x34,0x04,0x01,0x01,0x02,0x0D,0x01,0x03,0x01,0x02,0x04,0x60,0x01 }, 14, AV_CODEC_ID_MPEG2VIDEO }, /* MPEG-ES Frame wrapped */
    { { 0x06,0x0E,0x2B,0x34,0x04,0x01,0x01,0x01,0x0D,0x01,0x03,0x01,0x02,0x02,0x41,0x01 }, 14,    AV_CODEC_ID_DVVIDEO }, /* DV 625 25mbps */
    { { 0x06,0x0E,0x2B,0x34,0x04,0x01,0x01,0x01,0x0D,0x01,0x03,0x01,0x02,0x05,0x00,0x00 }, 14,   AV_CODEC_ID_RAWVIDEO }, /* Uncompressed Picture */
    { { 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00 },  0,      AV_CODEC_ID_NONE },
930
};
931 932 933

/* EC ULs for intra-only formats */
static const MXFCodecUL mxf_intra_only_essence_container_uls[] = {
934 935
    { { 0x06,0x0E,0x2B,0x34,0x04,0x01,0x01,0x01,0x0D,0x01,0x03,0x01,0x02,0x01,0x00,0x00 }, 14, AV_CODEC_ID_MPEG2VIDEO }, /* MXF-GC SMPTE D-10 Mappings */
    { { 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00 },  0,       AV_CODEC_ID_NONE },
936 937 938 939
};

/* intra-only PictureEssenceCoding ULs, where no corresponding EC UL exists */
static const MXFCodecUL mxf_intra_only_picture_essence_coding_uls[] = {
940 941
    { { 0x06,0x0E,0x2B,0x34,0x04,0x01,0x01,0x0A,0x04,0x01,0x02,0x02,0x01,0x32,0x00,0x00 }, 14,       AV_CODEC_ID_H264 }, /* H.264/MPEG-4 AVC Intra Profiles */
    { { 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00 },  0,       AV_CODEC_ID_NONE },
942 943
};

944
static const MXFCodecUL mxf_sound_essence_container_uls[] = {
945
    // sound essence container uls
946 947 948 949 950
    { { 0x06,0x0E,0x2B,0x34,0x04,0x01,0x01,0x01,0x0D,0x01,0x03,0x01,0x02,0x06,0x01,0x00 }, 14, AV_CODEC_ID_PCM_S16LE }, /* BWF Frame wrapped */
    { { 0x06,0x0E,0x2B,0x34,0x04,0x01,0x01,0x02,0x0D,0x01,0x03,0x01,0x02,0x04,0x40,0x01 }, 14,       AV_CODEC_ID_MP2 }, /* MPEG-ES Frame wrapped, 0x40 ??? stream id */
    { { 0x06,0x0E,0x2B,0x34,0x04,0x01,0x01,0x01,0x0D,0x01,0x03,0x01,0x02,0x01,0x01,0x01 }, 14, AV_CODEC_ID_PCM_S16LE }, /* D-10 Mapping 50Mbps PAL Extended Template */
    { { 0x06,0x0E,0x2B,0x34,0x01,0x01,0x01,0xFF,0x4B,0x46,0x41,0x41,0x00,0x0D,0x4D,0x4F }, 14, AV_CODEC_ID_PCM_S16LE }, /* 0001GL00.MXF.A1.mxf_opatom.mxf */
    { { 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00 },  0,      AV_CODEC_ID_NONE },
951 952
};

953 954 955 956 957 958 959 960 961 962 963
static int mxf_get_sorted_table_segments(MXFContext *mxf, int *nb_sorted_segments, MXFIndexTableSegment ***sorted_segments)
{
    int i, j, nb_segments = 0;
    MXFIndexTableSegment **unsorted_segments;
    int last_body_sid = -1, last_index_sid = -1, last_index_start = -1;

    /* count number of segments, allocate arrays and copy unsorted segments */
    for (i = 0; i < mxf->metadata_sets_count; i++)
        if (mxf->metadata_sets[i]->type == IndexTableSegment)
            nb_segments++;

964 965 966
    if (!nb_segments)
        return AVERROR_INVALIDDATA;

967 968
    if (!(unsorted_segments = av_calloc(nb_segments, sizeof(*unsorted_segments))) ||
        !(*sorted_segments  = av_calloc(nb_segments, sizeof(**sorted_segments)))) {
969
        av_freep(sorted_segments);
970 971 972 973 974 975 976 977 978 979 980 981 982
        av_free(unsorted_segments);
        return AVERROR(ENOMEM);
    }

    for (i = j = 0; i < mxf->metadata_sets_count; i++)
        if (mxf->metadata_sets[i]->type == IndexTableSegment)
            unsorted_segments[j++] = (MXFIndexTableSegment*)mxf->metadata_sets[i];

    *nb_sorted_segments = 0;

    /* sort segments by {BodySID, IndexSID, IndexStartPosition}, remove duplicates while we're at it */
    for (i = 0; i < nb_segments; i++) {
        int best = -1, best_body_sid = -1, best_index_sid = -1, best_index_start = -1;
983
        uint64_t best_index_duration = 0;
984 985 986 987 988 989

        for (j = 0; j < nb_segments; j++) {
            MXFIndexTableSegment *s = unsorted_segments[j];

            /* Require larger BosySID, IndexSID or IndexStartPosition then the previous entry. This removes duplicates.
             * We want the smallest values for the keys than what we currently have, unless this is the first such entry this time around.
990
             * If we come across an entry with the same IndexStartPosition but larger IndexDuration, then we'll prefer it over the one we currently have.
991 992
             */
            if ((i == 0     || s->body_sid > last_body_sid || s->index_sid > last_index_sid || s->index_start_position > last_index_start) &&
993 994
                (best == -1 || s->body_sid < best_body_sid || s->index_sid < best_index_sid || s->index_start_position < best_index_start ||
                (s->index_start_position == best_index_start && s->index_duration > best_index_duration))) {
995 996 997 998
                best             = j;
                best_body_sid    = s->body_sid;
                best_index_sid   = s->index_sid;
                best_index_start = s->index_start_position;
999
                best_index_duration = s->index_duration;
1000 1001 1002 1003 1004 1005 1006 1007 1008 1009 1010 1011 1012 1013 1014 1015 1016 1017
            }
        }

        /* no suitable entry found -> we're done */
        if (best == -1)
            break;

        (*sorted_segments)[(*nb_sorted_segments)++] = unsorted_segments[best];
        last_body_sid    = best_body_sid;
        last_index_sid   = best_index_sid;
        last_index_start = best_index_start;
    }

    av_free(unsorted_segments);

    return 0;
}

1018 1019 1020 1021 1022 1023 1024 1025 1026 1027 1028 1029 1030 1031 1032 1033 1034 1035 1036 1037 1038 1039
/**
 * Computes the absolute file offset of the given essence container offset
 */
static int mxf_absolute_bodysid_offset(MXFContext *mxf, int body_sid, int64_t offset, int64_t *offset_out)
{
    int x;
    int64_t offset_in = offset;     /* for logging */

    for (x = 0; x < mxf->partitions_count; x++) {
        MXFPartition *p = &mxf->partitions[x];

        if (p->body_sid != body_sid)
            continue;

        if (offset < p->essence_length || !p->essence_length) {
            *offset_out = p->essence_offset + offset;
            return 0;
        }

        offset -= p->essence_length;
    }

1040 1041
    av_log(mxf->fc, AV_LOG_ERROR,
           "failed to find absolute offset of %"PRIX64" in BodySID %i - partial file?\n",
1042 1043 1044 1045 1046
           offset_in, body_sid);

    return AVERROR_INVALIDDATA;
}

1047 1048 1049 1050 1051 1052 1053 1054 1055 1056 1057 1058 1059 1060 1061 1062 1063 1064 1065 1066 1067 1068 1069
/**
 * Returns the end position of the essence container with given BodySID, or zero if unknown
 */
static int64_t mxf_essence_container_end(MXFContext *mxf, int body_sid)
{
    int x;
    int64_t ret = 0;

    for (x = 0; x < mxf->partitions_count; x++) {
        MXFPartition *p = &mxf->partitions[x];

        if (p->body_sid != body_sid)
            continue;

        if (!p->essence_length)
            return 0;

        ret = p->essence_offset + p->essence_length;
    }

    return ret;
}

1070 1071 1072 1073
/* EditUnit -> absolute offset */
static int mxf_edit_unit_absolute_offset(MXFContext *mxf, MXFIndexTable *index_table, int64_t edit_unit, int64_t *edit_unit_out, int64_t *offset_out, int nag)
{
    int i;
1074
    int64_t offset_temp = 0;
1075 1076 1077 1078 1079 1080 1081 1082 1083 1084 1085 1086 1087 1088 1089

    for (i = 0; i < index_table->nb_segments; i++) {
        MXFIndexTableSegment *s = index_table->segments[i];

        edit_unit = FFMAX(edit_unit, s->index_start_position);  /* clamp if trying to seek before start */

        if (edit_unit < s->index_start_position + s->index_duration) {
            int64_t index = edit_unit - s->index_start_position;

            if (s->edit_unit_byte_count)
                offset_temp += s->edit_unit_byte_count * index;
            else if (s->nb_index_entries) {
                if (s->nb_index_entries == 2 * s->index_duration + 1)
                    index *= 2;     /* Avid index */

1090
                if (index < 0 || index >= s->nb_index_entries) {
1091 1092 1093 1094 1095 1096 1097 1098 1099 1100 1101 1102 1103 1104 1105 1106 1107 1108 1109 1110 1111 1112 1113 1114 1115 1116 1117 1118
                    av_log(mxf->fc, AV_LOG_ERROR, "IndexSID %i segment at %"PRId64" IndexEntryArray too small\n",
                           index_table->index_sid, s->index_start_position);
                    return AVERROR_INVALIDDATA;
                }

                offset_temp = s->stream_offset_entries[index];
            } else {
                av_log(mxf->fc, AV_LOG_ERROR, "IndexSID %i segment at %"PRId64" missing EditUnitByteCount and IndexEntryArray\n",
                       index_table->index_sid, s->index_start_position);
                return AVERROR_INVALIDDATA;
            }

            if (edit_unit_out)
                *edit_unit_out = edit_unit;

            return mxf_absolute_bodysid_offset(mxf, index_table->body_sid, offset_temp, offset_out);
        } else {
            /* EditUnitByteCount == 0 for VBR indexes, which is fine since they use explicit StreamOffsets */
            offset_temp += s->edit_unit_byte_count * s->index_duration;
        }
    }

    if (nag)
        av_log(mxf->fc, AV_LOG_ERROR, "failed to map EditUnit %"PRId64" in IndexSID %i to an offset\n", edit_unit, index_table->index_sid);

    return AVERROR_INVALIDDATA;
}

1119
static int mxf_compute_ptses_fake_index(MXFContext *mxf, MXFIndexTable *index_table)
1120
{
1121
    int i, j, x;
1122 1123 1124
    int8_t max_temporal_offset = -128;

    /* first compute how many entries we have */
1125 1126
    for (i = 0; i < index_table->nb_segments; i++) {
        MXFIndexTableSegment *s = index_table->segments[i];
1127

1128 1129
        if (!s->nb_index_entries) {
            index_table->nb_ptses = 0;
1130
            return 0;                               /* no TemporalOffsets */
1131
        }
1132

1133
        index_table->nb_ptses += s->index_duration;
1134 1135 1136
    }

    /* paranoid check */
1137
    if (index_table->nb_ptses <= 0)
1138 1139
        return 0;

1140 1141 1142
    if (!(index_table->ptses      = av_calloc(index_table->nb_ptses, sizeof(int64_t))) ||
        !(index_table->fake_index = av_calloc(index_table->nb_ptses, sizeof(AVIndexEntry)))) {
        av_freep(&index_table->ptses);
1143
        return AVERROR(ENOMEM);
1144
    }
1145 1146 1147

    /* we may have a few bad TemporalOffsets
     * make sure the corresponding PTSes don't have the bogus value 0 */
1148 1149
    for (x = 0; x < index_table->nb_ptses; x++)
        index_table->ptses[x] = AV_NOPTS_VALUE;
1150 1151 1152 1153 1154 1155 1156 1157 1158 1159 1160 1161 1162 1163 1164 1165 1166 1167 1168 1169 1170 1171 1172 1173 1174 1175 1176 1177

    /**
     * We have this:
     *
     * x  TemporalOffset
     * 0:  0
     * 1:  1
     * 2:  1
     * 3: -2
     * 4:  1
     * 5:  1
     * 6: -2
     *
     * We want to transform it into this:
     *
     * x  DTS PTS
     * 0: -1   0
     * 1:  0   3
     * 2:  1   1
     * 3:  2   2
     * 4:  3   6
     * 5:  4   4
     * 6:  5   5
     *
     * We do this by bucket sorting x by x+TemporalOffset[x] into mxf->ptses,
     * then settings mxf->first_dts = -max(TemporalOffset[x]).
     * The latter makes DTS <= PTS.
     */
1178 1179
    for (i = x = 0; i < index_table->nb_segments; i++) {
        MXFIndexTableSegment *s = index_table->segments[i];
1180
        int index_delta = 1;
1181
        int n = s->nb_index_entries;
1182

1183
        if (s->nb_index_entries == 2 * s->index_duration + 1) {
1184
            index_delta = 2;    /* Avid index */
1185 1186 1187 1188 1189
            /* ignore the last entry - it's the size of the essence container */
            n--;
        }

        for (j = 0; j < n; j += index_delta, x++) {
1190 1191 1192
            int offset = s->temporal_offset_entries[j] / index_delta;
            int index  = x + offset;

1193
            if (x >= index_table->nb_ptses) {
1194 1195
                av_log(mxf->fc, AV_LOG_ERROR,
                       "x >= nb_ptses - IndexEntryCount %i < IndexDuration %"PRId64"?\n",
1196 1197 1198 1199
                       s->nb_index_entries, s->index_duration);
                break;
            }

1200 1201 1202 1203
            index_table->fake_index[x].timestamp = x;
            index_table->fake_index[x].flags = !(s->flag_entries[j] & 0x30) ? AVINDEX_KEYFRAME : 0;

            if (index < 0 || index >= index_table->nb_ptses) {
1204 1205 1206 1207 1208 1209
                av_log(mxf->fc, AV_LOG_ERROR,
                       "index entry %i + TemporalOffset %i = %i, which is out of bounds\n",
                       x, offset, index);
                continue;
            }

1210
            index_table->ptses[index] = x;
1211 1212 1213 1214
            max_temporal_offset = FFMAX(max_temporal_offset, offset);
        }
    }

1215
    index_table->first_dts = -max_temporal_offset;
1216 1217 1218 1219

    return 0;
}

1220 1221 1222 1223 1224 1225
/**
 * Sorts and collects index table segments into index tables.
 * Also computes PTSes if possible.
 */
static int mxf_compute_index_tables(MXFContext *mxf)
{
1226
    int i, j, k, ret, nb_sorted_segments;
1227 1228 1229 1230 1231 1232 1233 1234 1235 1236 1237 1238 1239 1240 1241 1242 1243 1244 1245 1246 1247 1248 1249 1250 1251 1252 1253 1254 1255 1256 1257 1258 1259 1260 1261 1262 1263 1264 1265 1266 1267 1268 1269 1270 1271 1272 1273 1274 1275 1276 1277
    MXFIndexTableSegment **sorted_segments = NULL;

    if ((ret = mxf_get_sorted_table_segments(mxf, &nb_sorted_segments, &sorted_segments)) ||
        nb_sorted_segments <= 0) {
        av_log(mxf->fc, AV_LOG_WARNING, "broken or empty index\n");
        return 0;
    }

    /* sanity check and count unique BodySIDs/IndexSIDs */
    for (i = 0; i < nb_sorted_segments; i++) {
        if (i == 0 || sorted_segments[i-1]->index_sid != sorted_segments[i]->index_sid)
            mxf->nb_index_tables++;
        else if (sorted_segments[i-1]->body_sid != sorted_segments[i]->body_sid) {
            av_log(mxf->fc, AV_LOG_ERROR, "found inconsistent BodySID\n");
            ret = AVERROR_INVALIDDATA;
            goto finish_decoding_index;
        }
    }

    if (!(mxf->index_tables = av_calloc(mxf->nb_index_tables, sizeof(MXFIndexTable)))) {
        av_log(mxf->fc, AV_LOG_ERROR, "failed to allocate index tables\n");
        ret = AVERROR(ENOMEM);
        goto finish_decoding_index;
    }

    /* distribute sorted segments to index tables */
    for (i = j = 0; i < nb_sorted_segments; i++) {
        if (i != 0 && sorted_segments[i-1]->index_sid != sorted_segments[i]->index_sid) {
            /* next IndexSID */
            j++;
        }

        mxf->index_tables[j].nb_segments++;
    }

    for (i = j = 0; j < mxf->nb_index_tables; i += mxf->index_tables[j++].nb_segments) {
        MXFIndexTable *t = &mxf->index_tables[j];

        if (!(t->segments = av_calloc(t->nb_segments, sizeof(MXFIndexTableSegment*)))) {
            av_log(mxf->fc, AV_LOG_ERROR, "failed to allocate IndexTableSegment pointer array\n");
            ret = AVERROR(ENOMEM);
            goto finish_decoding_index;
        }

        if (sorted_segments[i]->index_start_position)
            av_log(mxf->fc, AV_LOG_WARNING, "IndexSID %i starts at EditUnit %"PRId64" - seeking may not work as expected\n",
                   sorted_segments[i]->index_sid, sorted_segments[i]->index_start_position);

        memcpy(t->segments, &sorted_segments[i], t->nb_segments * sizeof(MXFIndexTableSegment*));
        t->index_sid = sorted_segments[i]->index_sid;
        t->body_sid = sorted_segments[i]->body_sid;
1278 1279 1280

        if ((ret = mxf_compute_ptses_fake_index(mxf, t)) < 0)
            goto finish_decoding_index;
1281 1282 1283 1284 1285 1286 1287 1288 1289 1290 1291 1292 1293 1294 1295 1296 1297 1298 1299 1300 1301

        /* fix zero IndexDurations */
        for (k = 0; k < t->nb_segments; k++) {
            if (t->segments[k]->index_duration)
                continue;

            if (t->nb_segments > 1)
                av_log(mxf->fc, AV_LOG_WARNING, "IndexSID %i segment %i has zero IndexDuration and there's more than one segment\n",
                       t->index_sid, k);

            if (mxf->fc->nb_streams <= 0) {
                av_log(mxf->fc, AV_LOG_WARNING, "no streams?\n");
                break;
            }

            /* assume the first stream's duration is reasonable
             * leave index_duration = 0 on further segments in case we have any (unlikely)
             */
            t->segments[k]->index_duration = mxf->fc->streams[0]->duration;
            break;
        }
1302 1303 1304 1305 1306 1307 1308 1309
    }

    ret = 0;
finish_decoding_index:
    av_free(sorted_segments);
    return ret;
}

1310 1311 1312
static int mxf_is_intra_only(MXFDescriptor *descriptor)
{
    return mxf_get_codec_ul(mxf_intra_only_essence_container_uls,
1313
                            &descriptor->essence_container_ul)->id != AV_CODEC_ID_NONE ||
1314
           mxf_get_codec_ul(mxf_intra_only_picture_essence_coding_uls,
1315
                            &descriptor->essence_codec_ul)->id     != AV_CODEC_ID_NONE;
1316 1317
}

1318 1319 1320 1321 1322 1323 1324 1325
static int mxf_add_timecode_metadata(AVDictionary **pm, const char *key, AVTimecode *tc)
{
    char buf[AV_TIMECODE_STR_SIZE];
    av_dict_set(pm, key, av_timecode_make_string(tc, buf, 0), 0);

    return 0;
}

1326 1327 1328 1329
static int mxf_parse_structural_metadata(MXFContext *mxf)
{
    MXFPackage *material_package = NULL;
    MXFPackage *temp_package = NULL;
1330
    int i, j, k, ret;
1331

1332
    av_dlog(mxf->fc, "metadata sets count %d\n", mxf->metadata_sets_count);
1333 1334 1335 1336 1337 1338 1339
    /* TODO: handle multiple material packages (OP3x) */
    for (i = 0; i < mxf->packages_count; i++) {
        material_package = mxf_resolve_strong_ref(mxf, &mxf->packages_refs[i], MaterialPackage);
        if (material_package) break;
    }
    if (!material_package) {
        av_log(mxf->fc, AV_LOG_ERROR, "no material package found\n");
1340
        return AVERROR_INVALIDDATA;
1341 1342 1343 1344 1345 1346 1347 1348 1349
    }

    for (i = 0; i < material_package->tracks_count; i++) {
        MXFPackage *source_package = NULL;
        MXFTrack *material_track = NULL;
        MXFTrack *source_track = NULL;
        MXFTrack *temp_track = NULL;
        MXFDescriptor *descriptor = NULL;
        MXFStructuralComponent *component = NULL;
1350
        MXFTimecodeComponent *mxf_tc = NULL;
1351 1352 1353
        UID *essence_container_ul = NULL;
        const MXFCodecUL *codec_ul = NULL;
        const MXFCodecUL *container_ul = NULL;
1354
        const MXFCodecUL *pix_fmt_ul = NULL;
1355
        AVStream *st;
1356 1357
        AVTimecode tc;
        int flags;
1358 1359 1360 1361 1362 1363

        if (!(material_track = mxf_resolve_strong_ref(mxf, &material_package->tracks_refs[i], Track))) {
            av_log(mxf->fc, AV_LOG_ERROR, "could not resolve material track strong ref\n");
            continue;
        }

1364 1365 1366
        if ((component = mxf_resolve_strong_ref(mxf, &material_track->sequence_ref, TimecodeComponent))) {
            mxf_tc = (MXFTimecodeComponent*)component;
            flags = mxf_tc->drop_frame == 1 ? AV_TIMECODE_FLAG_DROPFRAME : 0;
1367
            if (av_timecode_init(&tc, mxf_tc->rate, flags, mxf_tc->start_frame, mxf->fc) == 0) {
1368 1369 1370 1371
                mxf_add_timecode_metadata(&mxf->fc->metadata, "timecode", &tc);
            }
        }

1372 1373
        if (!(material_track->sequence = mxf_resolve_strong_ref(mxf, &material_track->sequence_ref, Sequence))) {
            av_log(mxf->fc, AV_LOG_ERROR, "could not resolve material track sequence strong ref\n");
1374
            continue;
1375 1376
        }

1377 1378 1379 1380 1381 1382 1383
        for (j = 0; j < material_track->sequence->structural_components_count; j++) {
            component = mxf_resolve_strong_ref(mxf, &material_track->sequence->structural_components_refs[j], TimecodeComponent);
            if (!component)
                continue;

            mxf_tc = (MXFTimecodeComponent*)component;
            flags = mxf_tc->drop_frame == 1 ? AV_TIMECODE_FLAG_DROPFRAME : 0;
1384
            if (av_timecode_init(&tc, mxf_tc->rate, flags, mxf_tc->start_frame, mxf->fc) == 0) {
1385 1386 1387 1388 1389
                mxf_add_timecode_metadata(&mxf->fc->metadata, "timecode", &tc);
                break;
            }
        }

1390 1391 1392 1393 1394 1395 1396 1397 1398 1399 1400 1401 1402 1403 1404 1405
        /* TODO: handle multiple source clips */
        for (j = 0; j < material_track->sequence->structural_components_count; j++) {
            component = mxf_resolve_strong_ref(mxf, &material_track->sequence->structural_components_refs[j], SourceClip);
            if (!component)
                continue;

            for (k = 0; k < mxf->packages_count; k++) {
                temp_package = mxf_resolve_strong_ref(mxf, &mxf->packages_refs[k], SourcePackage);
                if (!temp_package)
                    continue;
                if (!memcmp(temp_package->package_uid, component->source_package_uid, 16)) {
                    source_package = temp_package;
                    break;
                }
            }
            if (!source_package) {
1406
                av_dlog(mxf->fc, "material track %d: no corresponding source package found\n", material_track->track_id);
1407 1408 1409 1410 1411
                break;
            }
            for (k = 0; k < source_package->tracks_count; k++) {
                if (!(temp_track = mxf_resolve_strong_ref(mxf, &source_package->tracks_refs[k], Track))) {
                    av_log(mxf->fc, AV_LOG_ERROR, "could not resolve source track strong ref\n");
1412
                    ret = AVERROR_INVALIDDATA;
1413
                    goto fail_and_free;
1414 1415 1416 1417 1418 1419 1420 1421 1422 1423 1424
                }
                if (temp_track->track_id == component->source_track_id) {
                    source_track = temp_track;
                    break;
                }
            }
            if (!source_track) {
                av_log(mxf->fc, AV_LOG_ERROR, "material track %d: no corresponding source track found\n", material_track->track_id);
                break;
            }
        }
1425
        if (!source_track || !component)
1426 1427
            continue;

1428 1429
        if (!(source_track->sequence = mxf_resolve_strong_ref(mxf, &source_track->sequence_ref, Sequence))) {
            av_log(mxf->fc, AV_LOG_ERROR, "could not resolve source track sequence strong ref\n");
1430
            ret = AVERROR_INVALIDDATA;
1431
            goto fail_and_free;
1432 1433 1434 1435 1436 1437 1438 1439 1440
        }

        /* 0001GL00.MXF.A1.mxf_opatom.mxf has the same SourcePackageID as 0001GL.MXF.V1.mxf_opatom.mxf
         * This would result in both files appearing to have two streams. Work around this by sanity checking DataDefinition */
        if (memcmp(material_track->sequence->data_definition_ul, source_track->sequence->data_definition_ul, 16)) {
            av_log(mxf->fc, AV_LOG_ERROR, "material track %d: DataDefinition mismatch\n", material_track->track_id);
            continue;
        }

1441
        st = avformat_new_stream(mxf->fc, NULL);
1442 1443
        if (!st) {
            av_log(mxf->fc, AV_LOG_ERROR, "could not allocate stream\n");
1444 1445
            ret = AVERROR(ENOMEM);
            goto fail_and_free;
1446
        }
1447
        st->id = source_track->track_id;
1448
        st->priv_data = source_track;
1449
        source_track->original_duration = st->duration = component->duration;
1450 1451 1452
        if (st->duration == -1)
            st->duration = AV_NOPTS_VALUE;
        st->start_time = component->start_position;
1453 1454 1455 1456 1457 1458
        if (material_track->edit_rate.num <= 0 || material_track->edit_rate.den <= 0) {
            av_log(mxf->fc, AV_LOG_WARNING,
                   "invalid edit rate (%d/%d) found on stream #%d, defaulting to 25/1\n",
                   material_track->edit_rate.num, material_track->edit_rate.den, st->index);
            material_track->edit_rate = (AVRational){25, 1};
        }
1459
        avpriv_set_pts_info(st, 64, material_track->edit_rate.den, material_track->edit_rate.num);
1460

1461 1462 1463
        /* ensure SourceTrack EditRate == MaterialTrack EditRate since only the former is accessible via st->priv_data */
        source_track->edit_rate = material_track->edit_rate;

1464
        PRINT_KEY(mxf->fc, "data definition   ul", source_track->sequence->data_definition_ul);
1465 1466
        codec_ul = mxf_get_codec_ul(ff_mxf_data_definition_uls, &source_track->sequence->data_definition_ul);
        st->codec->codec_type = codec_ul->id;
1467 1468 1469 1470 1471 1472 1473 1474 1475 1476 1477 1478 1479 1480 1481 1482 1483 1484 1485 1486 1487 1488 1489 1490 1491 1492 1493 1494 1495 1496 1497 1498 1499 1500 1501 1502 1503 1504

        source_package->descriptor = mxf_resolve_strong_ref(mxf, &source_package->descriptor_ref, AnyType);
        if (source_package->descriptor) {
            if (source_package->descriptor->type == MultipleDescriptor) {
                for (j = 0; j < source_package->descriptor->sub_descriptors_count; j++) {
                    MXFDescriptor *sub_descriptor = mxf_resolve_strong_ref(mxf, &source_package->descriptor->sub_descriptors_refs[j], Descriptor);

                    if (!sub_descriptor) {
                        av_log(mxf->fc, AV_LOG_ERROR, "could not resolve sub descriptor strong ref\n");
                        continue;
                    }
                    if (sub_descriptor->linked_track_id == source_track->track_id) {
                        descriptor = sub_descriptor;
                        break;
                    }
                }
            } else if (source_package->descriptor->type == Descriptor)
                descriptor = source_package->descriptor;
        }
        if (!descriptor) {
            av_log(mxf->fc, AV_LOG_INFO, "source track %d: stream %d, no descriptor found\n", source_track->track_id, st->index);
            continue;
        }
        PRINT_KEY(mxf->fc, "essence codec     ul", descriptor->essence_codec_ul);
        PRINT_KEY(mxf->fc, "essence container ul", descriptor->essence_container_ul);
        essence_container_ul = &descriptor->essence_container_ul;
        /* HACK: replacing the original key with mxf_encrypted_essence_container
         * is not allowed according to s429-6, try to find correct information anyway */
        if (IS_KLV_KEY(essence_container_ul, mxf_encrypted_essence_container)) {
            av_log(mxf->fc, AV_LOG_INFO, "broken encrypted mxf file\n");
            for (k = 0; k < mxf->metadata_sets_count; k++) {
                MXFMetadataSet *metadata = mxf->metadata_sets[k];
                if (metadata->type == CryptoContext) {
                    essence_container_ul = &((MXFCryptoContext *)metadata)->source_container_ul;
                    break;
                }
            }
        }
1505

1506
        /* TODO: drop PictureEssenceCoding and SoundEssenceCompression, only check EssenceContainer */
1507
        codec_ul = mxf_get_codec_ul(ff_mxf_codec_uls, &descriptor->essence_codec_ul);
1508
        st->codec->codec_id = (enum AVCodecID)codec_ul->id;
1509
        if (st->codec->codec_type == AVMEDIA_TYPE_VIDEO) {
1510
            source_track->intra_only = mxf_is_intra_only(descriptor);
1511
            container_ul = mxf_get_codec_ul(mxf_picture_essence_container_uls, essence_container_ul);
1512
            if (st->codec->codec_id == AV_CODEC_ID_NONE)
1513 1514
                st->codec->codec_id = container_ul->id;
            st->codec->width = descriptor->width;
1515 1516 1517 1518 1519 1520 1521
            st->codec->height = descriptor->height; /* Field height, not frame height */
            switch (descriptor->frame_layout) {
                case SegmentedFrame:
                    /* This one is a weird layout I don't fully understand. */
                    av_log(mxf->fc, AV_LOG_INFO, "SegmentedFrame layout isn't currently supported\n");
                    break;
                case FullFrame:
1522
                    st->codec->field_order = AV_FIELD_PROGRESSIVE;
1523 1524 1525 1526 1527 1528 1529 1530
                    break;
                case OneField:
                    /* Every other line is stored and needs to be duplicated. */
                    av_log(mxf->fc, AV_LOG_INFO, "OneField frame layout isn't currently supported\n");
                    break; /* The correct thing to do here is fall through, but by breaking we might be
                              able to decode some streams at half the vertical resolution, rather than not al all.
                              It's also for compatibility with the old behavior. */
                case MixedFields:
1531 1532
                    break;
                case SeparateFields:
1533
                    st->codec->height *= 2; /* Turn field height into frame height. */
1534
                    break;
1535 1536 1537
                default:
                    av_log(mxf->fc, AV_LOG_INFO, "Unknown frame layout type: %d\n", descriptor->frame_layout);
            }
1538
            if (st->codec->codec_id == AV_CODEC_ID_RAWVIDEO) {
1539
                st->codec->pix_fmt = descriptor->pix_fmt;
1540
                if (st->codec->pix_fmt == AV_PIX_FMT_NONE) {
1541 1542
                    pix_fmt_ul = mxf_get_codec_ul(ff_mxf_pixel_format_uls,
                                                  &descriptor->essence_codec_ul);
1543
                    st->codec->pix_fmt = (enum AVPixelFormat)pix_fmt_ul->id;
1544
                    if (st->codec->pix_fmt == AV_PIX_FMT_NONE) {
1545 1546 1547 1548 1549
                        /* support files created before RP224v10 by defaulting to UYVY422
                           if subsampling is 4:2:2 and component depth is 8-bit */
                        if (descriptor->horiz_subsampling == 2 &&
                            descriptor->vert_subsampling == 1 &&
                            descriptor->component_depth == 8) {
1550
                            st->codec->pix_fmt = AV_PIX_FMT_UYVY422;
1551 1552 1553 1554
                        }
                    }
                }
            }
1555
            st->need_parsing = AVSTREAM_PARSE_HEADERS;
1556
        } else if (st->codec->codec_type == AVMEDIA_TYPE_AUDIO) {
1557
            container_ul = mxf_get_codec_ul(mxf_sound_essence_container_uls, essence_container_ul);
1558
            /* Only overwrite existing codec ID if it is unset or A-law, which is the default according to SMPTE RP 224. */
1559 1560
            if (st->codec->codec_id == AV_CODEC_ID_NONE || (st->codec->codec_id == AV_CODEC_ID_PCM_ALAW && (enum AVCodecID)container_ul->id != AV_CODEC_ID_NONE))
                st->codec->codec_id = (enum AVCodecID)container_ul->id;
1561
            st->codec->channels = descriptor->channels;
1562
            st->codec->bits_per_coded_sample = descriptor->bits_per_sample;
1563

1564 1565
            if (descriptor->sample_rate.den > 0) {
                avpriv_set_pts_info(st, 64, descriptor->sample_rate.den, descriptor->sample_rate.num);
1566
                st->codec->sample_rate = descriptor->sample_rate.num / descriptor->sample_rate.den;
1567 1568 1569 1570 1571
            } else {
                av_log(mxf->fc, AV_LOG_WARNING, "invalid sample rate (%d/%d) found for stream #%d, time base forced to 1/48000\n",
                       descriptor->sample_rate.num, descriptor->sample_rate.den, st->index);
                avpriv_set_pts_info(st, 64, 1, 48000);
            }
1572

1573 1574 1575 1576
            /* if duration is set, rescale it from EditRate to SampleRate */
            if (st->duration != AV_NOPTS_VALUE)
                st->duration = av_rescale_q(st->duration, av_inv_q(material_track->edit_rate), st->time_base);

1577 1578
            /* TODO: implement AV_CODEC_ID_RAWAUDIO */
            if (st->codec->codec_id == AV_CODEC_ID_PCM_S16LE) {
1579
                if (descriptor->bits_per_sample > 16 && descriptor->bits_per_sample <= 24)
1580
                    st->codec->codec_id = AV_CODEC_ID_PCM_S24LE;
1581
                else if (descriptor->bits_per_sample == 32)
1582 1583
                    st->codec->codec_id = AV_CODEC_ID_PCM_S32LE;
            } else if (st->codec->codec_id == AV_CODEC_ID_PCM_S16BE) {
1584
                if (descriptor->bits_per_sample > 16 && descriptor->bits_per_sample <= 24)
1585
                    st->codec->codec_id = AV_CODEC_ID_PCM_S24BE;
1586
                else if (descriptor->bits_per_sample == 32)
1587 1588
                    st->codec->codec_id = AV_CODEC_ID_PCM_S32BE;
            } else if (st->codec->codec_id == AV_CODEC_ID_MP2) {
1589 1590 1591
                st->need_parsing = AVSTREAM_PARSE_FULL;
            }
        }
1592 1593 1594 1595 1596 1597 1598
        if (descriptor->extradata) {
            st->codec->extradata = av_mallocz(descriptor->extradata_size + FF_INPUT_BUFFER_PADDING_SIZE);
            if (st->codec->extradata)
                memcpy(st->codec->extradata, descriptor->extradata, descriptor->extradata_size);
        } else if(st->codec->codec_id == CODEC_ID_H264) {
            ff_generate_avci_extradata(st);
        }
1599
        if (st->codec->codec_type != AVMEDIA_TYPE_DATA && (*essence_container_ul)[15] > 0x01) {
1600 1601
            /* TODO: decode timestamps */
            st->need_parsing = AVSTREAM_PARSE_TIMESTAMPS;
1602 1603
        }
    }
1604

1605
    ret = 0;
1606 1607
fail_and_free:
    return ret;
1608 1609 1610 1611
}

static const MXFMetadataReadTableEntry mxf_metadata_read_table[] = {
    { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x05,0x01,0x00 }, mxf_read_primer_pack },
Tomas Härdin's avatar
Tomas Härdin committed
1612 1613 1614 1615 1616 1617 1618 1619 1620 1621
    { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x02,0x01,0x00 }, mxf_read_partition_pack },
    { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x02,0x02,0x00 }, mxf_read_partition_pack },
    { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x02,0x03,0x00 }, mxf_read_partition_pack },
    { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x02,0x04,0x00 }, mxf_read_partition_pack },
    { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x03,0x01,0x00 }, mxf_read_partition_pack },
    { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x03,0x02,0x00 }, mxf_read_partition_pack },
    { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x03,0x03,0x00 }, mxf_read_partition_pack },
    { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x03,0x04,0x00 }, mxf_read_partition_pack },
    { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x04,0x02,0x00 }, mxf_read_partition_pack },
    { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x04,0x04,0x00 }, mxf_read_partition_pack },
1622 1623 1624 1625 1626 1627 1628 1629 1630 1631 1632 1633 1634 1635
    { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x18,0x00 }, mxf_read_content_storage, 0, AnyType },
    { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x37,0x00 }, mxf_read_source_package, sizeof(MXFPackage), SourcePackage },
    { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x36,0x00 }, mxf_read_material_package, sizeof(MXFPackage), MaterialPackage },
    { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x0F,0x00 }, mxf_read_sequence, sizeof(MXFSequence), Sequence },
    { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x11,0x00 }, mxf_read_source_clip, sizeof(MXFStructuralComponent), SourceClip },
    { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x44,0x00 }, mxf_read_generic_descriptor, sizeof(MXFDescriptor), MultipleDescriptor },
    { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x42,0x00 }, mxf_read_generic_descriptor, sizeof(MXFDescriptor), Descriptor }, /* Generic Sound */
    { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x28,0x00 }, mxf_read_generic_descriptor, sizeof(MXFDescriptor), Descriptor }, /* CDCI */
    { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x29,0x00 }, mxf_read_generic_descriptor, sizeof(MXFDescriptor), Descriptor }, /* RGBA */
    { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x51,0x00 }, mxf_read_generic_descriptor, sizeof(MXFDescriptor), Descriptor }, /* MPEG 2 Video */
    { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x48,0x00 }, mxf_read_generic_descriptor, sizeof(MXFDescriptor), Descriptor }, /* Wave */
    { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x47,0x00 }, mxf_read_generic_descriptor, sizeof(MXFDescriptor), Descriptor }, /* AES3 */
    { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x3A,0x00 }, mxf_read_track, sizeof(MXFTrack), Track }, /* Static Track */
    { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x3B,0x00 }, mxf_read_track, sizeof(MXFTrack), Track }, /* Generic Track */
1636
    { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x14,0x00 }, mxf_read_timecode_component, sizeof(MXFTimecodeComponent), TimecodeComponent },
1637
    { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x04,0x01,0x02,0x02,0x00,0x00 }, mxf_read_cryptographic_context, sizeof(MXFCryptoContext), CryptoContext },
1638
    { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x10,0x01,0x00 }, mxf_read_index_table_segment, sizeof(MXFIndexTableSegment), IndexTableSegment },
1639 1640 1641
    { { 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00 }, NULL, 0, AnyType },
};

1642
static int mxf_read_local_tags(MXFContext *mxf, KLVPacket *klv, MXFMetadataReadFunc *read_child, int ctx_size, enum MXFMetadataSetType type)
1643
{
1644
    AVIOContext *pb = mxf->fc->pb;
1645
    MXFMetadataSet *ctx = ctx_size ? av_mallocz(ctx_size) : mxf;
1646
    uint64_t klv_end = avio_tell(pb) + klv->length;
1647 1648

    if (!ctx)
1649
        return AVERROR(ENOMEM);
1650
    while (avio_tell(pb) + 4 < klv_end && !url_feof(pb)) {
1651
        int ret;
1652 1653
        int tag = avio_rb16(pb);
        int size = avio_rb16(pb); /* KLV specified by 0x53 */
1654
        uint64_t next = avio_tell(pb) + size;
1655 1656
        UID uid = {0};

1657
        av_dlog(mxf->fc, "local tag %#04x size %d\n", tag, size);
1658
        if (!size) { /* ignore empty tag, needed for some files with empty UMID tag */
1659
            av_log(mxf->fc, AV_LOG_ERROR, "local tag %#04x with 0 size\n", tag);
1660 1661 1662 1663 1664 1665 1666 1667
            continue;
        }
        if (tag > 0x7FFF) { /* dynamic tag */
            int i;
            for (i = 0; i < mxf->local_tags_count; i++) {
                int local_tag = AV_RB16(mxf->local_tags+i*18);
                if (local_tag == tag) {
                    memcpy(uid, mxf->local_tags+i*18+2, 16);
1668
                    av_dlog(mxf->fc, "local tag %#04x\n", local_tag);
1669 1670 1671 1672 1673
                    PRINT_KEY(mxf->fc, "uid", uid);
                }
            }
        }
        if (ctx_size && tag == 0x3C0A)
1674
            avio_read(pb, ctx->uid, 16);
1675
        else if ((ret = read_child(ctx, pb, tag, size, uid, -1)) < 0)
1676
            return ret;
1677

1678 1679
        /* Accept the 64k local set limit being exceeded (Avid). Don't accept
         * it extending past the end of the KLV though (zzuf5.mxf). */
1680
        if (avio_tell(pb) > klv_end) {
1681 1682 1683
            if (ctx_size)
                av_free(ctx);

1684 1685
            av_log(mxf->fc, AV_LOG_ERROR,
                   "local tag %#04x extends past end of local set @ %#"PRIx64"\n",
1686 1687 1688
                   tag, klv->offset);
            return AVERROR_INVALIDDATA;
        } else if (avio_tell(pb) <= next)   /* only seek forward, else this can loop for a long time */
1689
            avio_seek(pb, next, SEEK_SET);
1690 1691 1692 1693 1694
    }
    if (ctx_size) ctx->type = type;
    return ctx_size ? mxf_add_metadata_set(mxf, ctx) : 0;
}

1695 1696 1697 1698 1699 1700 1701 1702 1703 1704 1705 1706 1707 1708 1709 1710 1711 1712 1713 1714 1715 1716 1717 1718 1719 1720 1721 1722 1723 1724 1725 1726 1727 1728 1729 1730 1731 1732 1733 1734 1735 1736 1737 1738 1739 1740 1741 1742 1743 1744 1745 1746 1747 1748 1749 1750 1751 1752 1753 1754 1755 1756 1757 1758 1759 1760 1761 1762 1763 1764 1765
/**
 * Seeks to the previous partition, if possible
 * @return <= 0 if we should stop parsing, > 0 if we should keep going
 */
static int mxf_seek_to_previous_partition(MXFContext *mxf)
{
    AVIOContext *pb = mxf->fc->pb;

    if (!mxf->current_partition ||
        mxf->run_in + mxf->current_partition->previous_partition <= mxf->last_forward_tell)
        return 0;   /* we've parsed all partitions */

    /* seek to previous partition */
    avio_seek(pb, mxf->run_in + mxf->current_partition->previous_partition, SEEK_SET);
    mxf->current_partition = NULL;

    av_dlog(mxf->fc, "seeking to previous partition\n");

    return 1;
}

/**
 * Called when essence is encountered
 * @return <= 0 if we should stop parsing, > 0 if we should keep going
 */
static int mxf_parse_handle_essence(MXFContext *mxf)
{
    AVIOContext *pb = mxf->fc->pb;
    int64_t ret;

    if (mxf->parsing_backward) {
        return mxf_seek_to_previous_partition(mxf);
    } else {
        if (!mxf->footer_partition) {
            av_dlog(mxf->fc, "no footer\n");
            return 0;
        }

        av_dlog(mxf->fc, "seeking to footer\n");

        /* remember where we were so we don't end up seeking further back than this */
        mxf->last_forward_tell = avio_tell(pb);

        if (!pb->seekable) {
            av_log(mxf->fc, AV_LOG_INFO, "file is not seekable - not parsing footer\n");
            return -1;
        }

        /* seek to footer partition and parse backward */
        if ((ret = avio_seek(pb, mxf->run_in + mxf->footer_partition, SEEK_SET)) < 0) {
            av_log(mxf->fc, AV_LOG_ERROR, "failed to seek to footer @ 0x%"PRIx64" (%"PRId64") - partial file?\n",
                   mxf->run_in + mxf->footer_partition, ret);
            return ret;
        }

        mxf->current_partition = NULL;
        mxf->parsing_backward = 1;
    }

    return 1;
}

/**
 * Called when the next partition or EOF is encountered
 * @return <= 0 if we should stop parsing, > 0 if we should keep going
 */
static int mxf_parse_handle_partition_or_eof(MXFContext *mxf)
{
    return mxf->parsing_backward ? mxf_seek_to_previous_partition(mxf) : 1;
}

1766 1767 1768 1769 1770 1771 1772 1773 1774 1775 1776 1777 1778 1779 1780 1781 1782 1783 1784 1785 1786 1787 1788 1789 1790 1791
/**
 * Figures out the proper offset and length of the essence container in each partition
 */
static void mxf_compute_essence_containers(MXFContext *mxf)
{
    int x;

    /* everything is already correct */
    if (mxf->op == OPAtom)
        return;

    for (x = 0; x < mxf->partitions_count; x++) {
        MXFPartition *p = &mxf->partitions[x];

        if (!p->body_sid)
            continue;       /* BodySID == 0 -> no essence */

        if (x >= mxf->partitions_count - 1)
            break;          /* last partition - can't compute length (and we don't need to) */

        /* essence container spans to the next partition */
        p->essence_length = mxf->partitions[x+1].this_partition - p->essence_offset;

        if (p->essence_length < 0) {
            /* next ThisPartition < essence_offset */
            p->essence_length = 0;
1792 1793
            av_log(mxf->fc, AV_LOG_ERROR,
                   "partition %i: bad ThisPartition = %"PRIX64"\n",
1794 1795 1796 1797 1798 1799 1800 1801 1802 1803 1804 1805 1806
                   x+1, mxf->partitions[x+1].this_partition);
        }
    }
}

static int64_t round_to_kag(int64_t position, int kag_size)
{
    /* TODO: account for run-in? the spec isn't clear whether KAG should account for it */
    /* NOTE: kag_size may be any integer between 1 - 2^10 */
    int64_t ret = (position / kag_size) * kag_size;
    return ret == position ? ret : ret + kag_size;
}

1807
static int is_pcm(enum AVCodecID codec_id)
1808 1809
{
    /* we only care about "normal" PCM codecs until we get samples */
1810
    return codec_id >= AV_CODEC_ID_PCM_S16LE && codec_id < AV_CODEC_ID_PCM_S24DAUD;
1811 1812 1813
}

/**
1814 1815 1816
 * Deal with the case where for some audio atoms EditUnitByteCount is
 * very small (2, 4..). In those cases we should read more than one
 * sample per call to mxf_read_packet().
1817 1818 1819 1820 1821 1822 1823 1824 1825 1826 1827
 */
static void mxf_handle_small_eubc(AVFormatContext *s)
{
    MXFContext *mxf = s->priv_data;

    /* assuming non-OPAtom == frame wrapped
     * no sane writer would wrap 2 byte PCM packets with 20 byte headers.. */
    if (mxf->op != OPAtom)
        return;

    /* expect PCM with exactly one index table segment and a small (< 32) EUBC */
1828 1829 1830 1831 1832
    if (s->nb_streams != 1                                     ||
        s->streams[0]->codec->codec_type != AVMEDIA_TYPE_AUDIO ||
        !is_pcm(s->streams[0]->codec->codec_id)                ||
        mxf->nb_index_tables != 1                              ||
        mxf->index_tables[0].nb_segments != 1                  ||
1833 1834 1835 1836
        mxf->index_tables[0].segments[0]->edit_unit_byte_count >= 32)
        return;

    /* arbitrarily default to 48 kHz PAL audio frame size */
1837 1838 1839
    /* TODO: We could compute this from the ratio between the audio
     *       and video edit rates for 48 kHz NTSC we could use the
     *       1802-1802-1802-1802-1801 pattern. */
1840 1841 1842
    mxf->edit_units_per_packet = 1920;
}

1843
static int mxf_read_header(AVFormatContext *s)
1844 1845 1846
{
    MXFContext *mxf = s->priv_data;
    KLVPacket klv;
1847
    int64_t essence_offset = 0;
1848
    int ret;
1849

1850
    mxf->last_forward_tell = INT64_MAX;
1851
    mxf->edit_units_per_packet = 1;
1852

1853 1854
    if (!mxf_read_sync(s->pb, mxf_header_partition_pack_key, 14)) {
        av_log(s, AV_LOG_ERROR, "could not find header partition pack key\n");
1855
        return AVERROR_INVALIDDATA;
1856
    }
1857
    avio_seek(s->pb, -14, SEEK_CUR);
1858
    mxf->fc = s;
1859 1860
    mxf->run_in = avio_tell(s->pb);

1861 1862 1863
    while (!url_feof(s->pb)) {
        const MXFMetadataReadTableEntry *metadata;

1864 1865 1866 1867 1868 1869 1870 1871
        if (klv_read_packet(&klv, s->pb) < 0) {
            /* EOF - seek to previous partition or stop */
            if(mxf_parse_handle_partition_or_eof(mxf) <= 0)
                break;
            else
                continue;
        }

1872
        PRINT_KEY(s, "read header", klv.key);
1873
        av_dlog(s, "size %"PRIu64" offset %#"PRIx64"\n", klv.length, klv.offset);
1874
        if (IS_KLV_KEY(klv.key, mxf_encrypted_triplet_key) ||
1875
            IS_KLV_KEY(klv.key, mxf_essence_element_key) ||
1876
            IS_KLV_KEY(klv.key, mxf_avid_essence_element_key) ||
1877
            IS_KLV_KEY(klv.key, mxf_system_item_key)) {
1878 1879 1880 1881 1882 1883

            if (!mxf->current_partition) {
                av_log(mxf->fc, AV_LOG_ERROR, "found essence prior to first PartitionPack\n");
                return AVERROR_INVALIDDATA;
            }

1884 1885 1886 1887 1888 1889 1890 1891 1892 1893 1894 1895 1896 1897 1898 1899 1900 1901 1902 1903 1904 1905 1906 1907
            if (!mxf->current_partition->essence_offset) {
                /* for OP1a we compute essence_offset
                 * for OPAtom we point essence_offset after the KL (usually op1a_essence_offset + 20 or 25)
                 * TODO: for OP1a we could eliminate this entire if statement, always stopping parsing at op1a_essence_offset
                 *       for OPAtom we still need the actual essence_offset though (the KL's length can vary)
                 */
                int64_t op1a_essence_offset =
                    round_to_kag(mxf->current_partition->this_partition +
                                 mxf->current_partition->pack_length,       mxf->current_partition->kag_size) +
                    round_to_kag(mxf->current_partition->header_byte_count, mxf->current_partition->kag_size) +
                    round_to_kag(mxf->current_partition->index_byte_count,  mxf->current_partition->kag_size);

                if (mxf->op == OPAtom) {
                    /* point essence_offset to the actual data
                    * OPAtom has all the essence in one big KLV
                    */
                    mxf->current_partition->essence_offset = avio_tell(s->pb);
                    mxf->current_partition->essence_length = klv.length;
                } else {
                    /* NOTE: op1a_essence_offset may be less than to klv.offset (C0023S01.mxf)  */
                    mxf->current_partition->essence_offset = op1a_essence_offset;
                }
            }

1908 1909
            if (!essence_offset)
                essence_offset = klv.offset;
1910 1911 1912 1913

            /* seek to footer, previous partition or stop */
            if (mxf_parse_handle_essence(mxf) <= 0)
                break;
1914
            continue;
1915 1916 1917 1918 1919
        } else if (!memcmp(klv.key, mxf_header_partition_pack_key, 13) &&
                   klv.key[13] >= 2 && klv.key[13] <= 4 && mxf->current_partition) {
            /* next partition pack - keep going, seek to previous partition or stop */
            if(mxf_parse_handle_partition_or_eof(mxf) <= 0)
                break;
1920 1921 1922
            else if (mxf->parsing_backward)
                continue;
            /* we're still parsing forward. proceed to parsing this partition pack */
1923
        }
1924 1925 1926

        for (metadata = mxf_metadata_read_table; metadata->read; metadata++) {
            if (IS_KLV_KEY(klv.key, metadata->key)) {
1927 1928 1929
                int res;
                if (klv.key[5] == 0x53) {
                    res = mxf_read_local_tags(mxf, &klv, metadata->read, metadata->ctx_size, metadata->type);
1930 1931
                } else {
                    uint64_t next = avio_tell(s->pb) + klv.length;
1932
                    res = metadata->read(mxf, s->pb, 0, klv.length, klv.key, klv.offset);
1933 1934 1935 1936 1937 1938 1939 1940

                    /* only seek forward, else this can loop for a long time */
                    if (avio_tell(s->pb) > next) {
                        av_log(s, AV_LOG_ERROR, "read past end of KLV @ %#"PRIx64"\n",
                               klv.offset);
                        return AVERROR_INVALIDDATA;
                    }

1941 1942
                    avio_seek(s->pb, next, SEEK_SET);
                }
1943
                if (res < 0) {
1944
                    av_log(s, AV_LOG_ERROR, "error reading header metadata\n");
1945
                    return res;
1946 1947 1948 1949 1950
                }
                break;
            }
        }
        if (!metadata->read)
1951
            avio_skip(s->pb, klv.length);
1952
    }
1953
    /* FIXME avoid seek */
1954
    if (!essence_offset)  {
1955 1956 1957
        av_log(s, AV_LOG_ERROR, "no essence\n");
        return AVERROR_INVALIDDATA;
    }
1958
    avio_seek(s->pb, essence_offset, SEEK_SET);
1959 1960 1961

    mxf_compute_essence_containers(mxf);

1962 1963 1964
    /* we need to do this before computing the index tables
     * to be able to fill in zero IndexDurations with st->duration */
    if ((ret = mxf_parse_structural_metadata(mxf)) < 0)
1965
        goto fail;
1966

1967
    if ((ret = mxf_compute_index_tables(mxf)) < 0)
1968
        goto fail;
1969 1970 1971 1972 1973 1974 1975

    if (mxf->nb_index_tables > 1) {
        /* TODO: look up which IndexSID to use via EssenceContainerData */
        av_log(mxf->fc, AV_LOG_INFO, "got %i index tables - only the first one (IndexSID %i) will be used\n",
               mxf->nb_index_tables, mxf->index_tables[0].index_sid);
    } else if (mxf->nb_index_tables == 0 && mxf->op == OPAtom) {
        av_log(mxf->fc, AV_LOG_ERROR, "cannot demux OPAtom without an index\n");
1976 1977
        ret = AVERROR_INVALIDDATA;
        goto fail;
1978 1979
    }

1980 1981
    mxf_handle_small_eubc(s);

1982
    return 0;
1983 1984 1985 1986
fail:
    mxf_read_close(s);

    return ret;
1987 1988
}

1989
/**
1990 1991
 * Sets mxf->current_edit_unit based on what offset we're currently at.
 * @return next_ofs if OK, <0 on error
1992
 */
1993
static int64_t mxf_set_current_edit_unit(MXFContext *mxf, int64_t current_offset)
1994
{
1995
    int64_t last_ofs = -1, next_ofs = -1;
1996 1997
    MXFIndexTable *t = &mxf->index_tables[0];

1998 1999
    /* this is called from the OP1a demuxing logic, which means there
     * may be no index tables */
2000
    if (mxf->nb_index_tables <= 0)
2001
        return -1;
2002

2003
    /* find mxf->current_edit_unit so that the next edit unit starts ahead of current_offset */
2004
    while (mxf->current_edit_unit >= 0) {
2005
        if (mxf_edit_unit_absolute_offset(mxf, t, mxf->current_edit_unit + 1, NULL, &next_ofs, 0) < 0)
2006
            return -1;
2007

2008
        if (next_ofs <= last_ofs) {
2009 2010 2011 2012
            /* large next_ofs didn't change or current_edit_unit wrapped
             * around this fixes the infinite loop on zzuf3.mxf */
            av_log(mxf->fc, AV_LOG_ERROR,
                   "next_ofs didn't change. not deriving packet timestamps\n");
2013
            return -1;
2014 2015
        }

2016
        if (next_ofs > current_offset)
2017 2018
            break;

2019
        last_ofs = next_ofs;
2020 2021 2022
        mxf->current_edit_unit++;
    }

2023 2024 2025
    /* not checking mxf->current_edit_unit >= t->nb_ptses here since CBR files may lack IndexEntryArrays */
    if (mxf->current_edit_unit < 0)
        return -1;
2026

2027
    return next_ofs;
2028 2029
}

2030 2031 2032 2033 2034 2035 2036 2037 2038 2039 2040 2041 2042 2043 2044 2045 2046 2047 2048 2049 2050 2051 2052 2053 2054 2055 2056 2057 2058
static int mxf_compute_sample_count(MXFContext *mxf, int stream_index, uint64_t *sample_count)
{
    int i, total = 0, size = 0;
    AVStream *st = mxf->fc->streams[stream_index];
    MXFTrack *track = st->priv_data;
    AVRational time_base = av_inv_q(track->edit_rate);
    AVRational sample_rate = av_inv_q(st->time_base);
    const MXFSamplesPerFrame *spf = NULL;

    if ((sample_rate.num / sample_rate.den) == 48000)
        spf = ff_mxf_get_samples_per_frame(mxf->fc, time_base);
    if (!spf) {
        int remainder = (sample_rate.num * time_base.num) % (time_base.den * sample_rate.den);
        *sample_count = av_q2d(av_mul_q((AVRational){mxf->current_edit_unit, 1},
                                        av_mul_q(sample_rate, time_base)));
        if (remainder)
            av_log(mxf->fc, AV_LOG_WARNING,
                   "seeking detected on stream #%d with time base (%d/%d) and sample rate (%d/%d), audio pts won't be accurate.\n",
                   stream_index, time_base.num, time_base.den, sample_rate.num, sample_rate.den);
        return 0;
    }

    while (spf->samples_per_frame[size]) {
        total += spf->samples_per_frame[size];
        size++;
    }

    av_assert2(size);

2059
    *sample_count = (mxf->current_edit_unit / size) * (uint64_t)total;
2060 2061 2062 2063 2064 2065 2066 2067 2068 2069 2070
    for (i = 0; i < mxf->current_edit_unit % size; i++) {
        *sample_count += spf->samples_per_frame[i];
    }

    return 0;
}

static int mxf_set_audio_pts(MXFContext *mxf, AVCodecContext *codec, AVPacket *pkt)
{
    MXFTrack *track = mxf->fc->streams[pkt->stream_index]->priv_data;
    pkt->pts = track->sample_count;
2071 2072
    if (codec->channels <= 0 || av_get_bits_per_sample(codec->codec_id) <= 0)
        return AVERROR(EINVAL);
2073
    track->sample_count += pkt->size / (codec->channels * (int64_t)av_get_bits_per_sample(codec->codec_id) / 8);
2074 2075 2076
    return 0;
}

2077 2078 2079
static int mxf_read_packet_old(AVFormatContext *s, AVPacket *pkt)
{
    KLVPacket klv;
2080
    MXFContext *mxf = s->priv_data;
2081 2082

    while (!url_feof(s->pb)) {
2083
        int ret;
2084 2085 2086 2087 2088
        if (klv_read_packet(&klv, s->pb) < 0)
            return -1;
        PRINT_KEY(s, "read packet", klv.key);
        av_dlog(s, "size %"PRIu64" offset %#"PRIx64"\n", klv.length, klv.offset);
        if (IS_KLV_KEY(klv.key, mxf_encrypted_triplet_key)) {
2089 2090
            ret = mxf_decrypt_triplet(s, pkt, &klv);
            if (ret < 0) {
2091
                av_log(s, AV_LOG_ERROR, "invalid encoded triplet\n");
2092
                return AVERROR_INVALIDDATA;
2093 2094 2095 2096 2097 2098
            }
            return 0;
        }
        if (IS_KLV_KEY(klv.key, mxf_essence_element_key) ||
            IS_KLV_KEY(klv.key, mxf_avid_essence_element_key)) {
            int index = mxf_get_stream_index(s, &klv);
2099
            int64_t next_ofs, next_klv;
2100 2101
            AVStream *st;
            MXFTrack *track;
2102
            AVCodecContext *codec;
2103

2104 2105 2106 2107
            if (index < 0) {
                av_log(s, AV_LOG_ERROR, "error getting stream index %d\n", AV_RB32(klv.key+12));
                goto skip;
            }
2108 2109 2110 2111

            st = s->streams[index];
            track = st->priv_data;

2112 2113
            if (s->streams[index]->discard == AVDISCARD_ALL)
                goto skip;
2114 2115 2116 2117 2118 2119 2120 2121 2122 2123 2124 2125 2126

            next_klv = avio_tell(s->pb) + klv.length;
            next_ofs = mxf_set_current_edit_unit(mxf, klv.offset);

            if (next_ofs >= 0 && next_klv > next_ofs) {
                /* if this check is hit then it's possible OPAtom was treated as OP1a
                 * truncate the packet since it's probably very large (>2 GiB is common) */
                av_log_ask_for_sample(s,
                    "KLV for edit unit %i extends into next edit unit - OPAtom misinterpreted as OP1a?\n",
                    mxf->current_edit_unit);
                klv.length = next_ofs - avio_tell(s->pb);
            }

2127 2128 2129 2130
            /* check for 8 channels AES3 element */
            if (klv.key[12] == 0x06 && klv.key[13] == 0x01 && klv.key[14] == 0x10) {
                if (mxf_get_d10_aes3_packet(s->pb, s->streams[index], pkt, klv.length) < 0) {
                    av_log(s, AV_LOG_ERROR, "error reading D-10 aes3 frame\n");
2131
                    return AVERROR_INVALIDDATA;
2132 2133
                }
            } else {
2134
                ret = av_get_packet(s->pb, pkt, klv.length);
2135 2136 2137 2138 2139
                if (ret < 0)
                    return ret;
            }
            pkt->stream_index = index;
            pkt->pos = klv.offset;
2140

2141 2142
            codec = s->streams[index]->codec;
            if (codec->codec_type == AVMEDIA_TYPE_VIDEO && next_ofs >= 0) {
2143 2144 2145 2146 2147 2148
                /* mxf->current_edit_unit good - see if we have an index table to derive timestamps from */
                MXFIndexTable *t = &mxf->index_tables[0];

                if (mxf->nb_index_tables >= 1 && mxf->current_edit_unit < t->nb_ptses) {
                    pkt->dts = mxf->current_edit_unit + t->first_dts;
                    pkt->pts = t->ptses[mxf->current_edit_unit];
2149 2150 2151 2152
                } else if (track->intra_only) {
                    /* intra-only -> PTS = EditUnit.
                     * let utils.c figure out DTS since it can be < PTS if low_delay = 0 (Sony IMX30) */
                    pkt->pts = mxf->current_edit_unit;
2153
                }
2154 2155 2156 2157
            } else if (codec->codec_type == AVMEDIA_TYPE_AUDIO) {
                int ret = mxf_set_audio_pts(mxf, codec, pkt);
                if (ret < 0)
                    return ret;
2158 2159 2160 2161
            }

            /* seek for truncated packets */
            avio_seek(s->pb, next_klv, SEEK_SET);
2162

2163 2164 2165 2166 2167 2168 2169 2170 2171 2172 2173
            return 0;
        } else
        skip:
            avio_skip(s->pb, klv.length);
    }
    return AVERROR_EOF;
}

static int mxf_read_packet(AVFormatContext *s, AVPacket *pkt)
{
    MXFContext *mxf = s->priv_data;
2174 2175
    int ret, size;
    int64_t ret64, pos, next_pos;
2176
    AVStream *st;
2177
    MXFIndexTable *t;
2178
    int edit_units;
2179

2180
    if (mxf->op != OPAtom)
2181 2182
        return mxf_read_packet_old(s, pkt);

2183
    /* OPAtom - clip wrapped demuxing */
2184
    /* NOTE: mxf_read_header() makes sure nb_index_tables > 0 for OPAtom */
2185 2186
    st = s->streams[0];
    t = &mxf->index_tables[0];
2187

2188
    if (mxf->current_edit_unit >= st->duration)
2189 2190
        return AVERROR_EOF;

2191 2192
    edit_units = FFMIN(mxf->edit_units_per_packet, st->duration - mxf->current_edit_unit);

2193 2194
    if ((ret = mxf_edit_unit_absolute_offset(mxf, t, mxf->current_edit_unit, NULL, &pos, 1)) < 0)
        return ret;
2195

2196 2197
    /* compute size by finding the next edit unit or the end of the essence container
     * not pretty, but it works */
2198
    if ((ret = mxf_edit_unit_absolute_offset(mxf, t, mxf->current_edit_unit + edit_units, NULL, &next_pos, 0)) < 0 &&
2199 2200 2201 2202
        (next_pos = mxf_essence_container_end(mxf, t->body_sid)) <= 0) {
        av_log(s, AV_LOG_ERROR, "unable to compute the size of the last packet\n");
        return AVERROR_INVALIDDATA;
    }
2203

2204 2205 2206 2207
    if ((size = next_pos - pos) <= 0) {
        av_log(s, AV_LOG_ERROR, "bad size: %i\n", size);
        return AVERROR_INVALIDDATA;
    }
2208

2209 2210
    if ((ret64 = avio_seek(s->pb, pos, SEEK_SET)) < 0)
        return ret64;
2211

2212 2213
    if ((size = av_get_packet(s->pb, pkt, size)) < 0)
        return size;
2214

2215 2216
    pkt->stream_index = 0;

2217 2218 2219 2220
    if (st->codec->codec_type == AVMEDIA_TYPE_VIDEO && t->ptses &&
        mxf->current_edit_unit >= 0 && mxf->current_edit_unit < t->nb_ptses) {
        pkt->dts = mxf->current_edit_unit + t->first_dts;
        pkt->pts = t->ptses[mxf->current_edit_unit];
2221 2222 2223 2224
    } else if (st->codec->codec_type == AVMEDIA_TYPE_AUDIO) {
        int ret = mxf_set_audio_pts(mxf, st->codec, pkt);
        if (ret < 0)
            return ret;
2225 2226
    }

2227
    mxf->current_edit_unit += edit_units;
2228 2229 2230 2231

    return 0;
}

2232 2233 2234
static int mxf_read_close(AVFormatContext *s)
{
    MXFContext *mxf = s->priv_data;
2235
    MXFIndexTableSegment *seg;
2236
    int i;
2237 2238

    av_freep(&mxf->packages_refs);
2239 2240 2241 2242

    for (i = 0; i < s->nb_streams; i++)
        s->streams[i]->priv_data = NULL;

2243 2244
    for (i = 0; i < mxf->metadata_sets_count; i++) {
        switch (mxf->metadata_sets[i]->type) {
2245 2246 2247
        case Descriptor:
            av_freep(&((MXFDescriptor *)mxf->metadata_sets[i])->extradata);
            break;
2248 2249 2250 2251 2252 2253 2254 2255 2256 2257
        case MultipleDescriptor:
            av_freep(&((MXFDescriptor *)mxf->metadata_sets[i])->sub_descriptors_refs);
            break;
        case Sequence:
            av_freep(&((MXFSequence *)mxf->metadata_sets[i])->structural_components_refs);
            break;
        case SourcePackage:
        case MaterialPackage:
            av_freep(&((MXFPackage *)mxf->metadata_sets[i])->tracks_refs);
            break;
2258 2259
        case IndexTableSegment:
            seg = (MXFIndexTableSegment *)mxf->metadata_sets[i];
2260
            av_freep(&seg->temporal_offset_entries);
2261 2262 2263
            av_freep(&seg->flag_entries);
            av_freep(&seg->stream_offset_entries);
            break;
2264 2265 2266 2267 2268
        default:
            break;
        }
        av_freep(&mxf->metadata_sets[i]);
    }
Tomas Härdin's avatar
Tomas Härdin committed
2269
    av_freep(&mxf->partitions);
2270 2271 2272
    av_freep(&mxf->metadata_sets);
    av_freep(&mxf->aesc);
    av_freep(&mxf->local_tags);
2273

2274 2275 2276 2277 2278 2279
    if (mxf->index_tables) {
        for (i = 0; i < mxf->nb_index_tables; i++) {
            av_freep(&mxf->index_tables[i].segments);
            av_freep(&mxf->index_tables[i].ptses);
            av_freep(&mxf->index_tables[i].fake_index);
        }
2280
    }
2281 2282
    av_freep(&mxf->index_tables);

2283 2284 2285 2286
    return 0;
}

static int mxf_probe(AVProbeData *p) {
2287 2288
    const uint8_t *bufp = p->buf;
    const uint8_t *end = p->buf + p->buf_size;
2289 2290 2291 2292 2293 2294 2295 2296 2297 2298 2299 2300 2301 2302 2303 2304 2305 2306 2307

    if (p->buf_size < sizeof(mxf_header_partition_pack_key))
        return 0;

    /* Must skip Run-In Sequence and search for MXF header partition pack key SMPTE 377M 5.5 */
    end -= sizeof(mxf_header_partition_pack_key);
    for (; bufp < end; bufp++) {
        if (IS_KLV_KEY(bufp, mxf_header_partition_pack_key))
            return AVPROBE_SCORE_MAX;
    }
    return 0;
}

/* rudimentary byte seek */
/* XXX: use MXF Index */
static int mxf_read_seek(AVFormatContext *s, int stream_index, int64_t sample_time, int flags)
{
    AVStream *st = s->streams[stream_index];
    int64_t seconds;
2308 2309
    MXFContext* mxf = s->priv_data;
    int64_t seekpos;
2310
    int i, ret;
2311
    MXFIndexTable *t;
2312 2313 2314 2315 2316
    MXFTrack *source_track = st->priv_data;

    /* if audio then truncate sample_time to EditRate */
    if (st->codec->codec_type == AVMEDIA_TYPE_AUDIO)
        sample_time = av_rescale_q(sample_time, st->time_base, av_inv_q(source_track->edit_rate));
2317

2318
    if (mxf->nb_index_tables <= 0) {
2319
    if (!s->bit_rate)
2320
        return AVERROR_INVALIDDATA;
2321 2322 2323
    if (sample_time < 0)
        sample_time = 0;
    seconds = av_rescale(sample_time, st->time_base.num, st->time_base.den);
2324 2325 2326

    if ((ret = avio_seek(s->pb, (s->bit_rate * seconds) >> 3, SEEK_SET)) < 0)
        return ret;
2327
    ff_update_cur_dts(s, st, sample_time);
2328
    } else {
2329
        t = &mxf->index_tables[0];
2330

2331 2332 2333
        /* clamp above zero, else ff_index_search_timestamp() returns negative
         * this also means we allow seeking before the start */
        sample_time = FFMAX(sample_time, 0);
2334

2335 2336 2337 2338 2339 2340 2341
        if (t->fake_index) {
            /* behave as if we have a proper index */
            if ((sample_time = ff_index_search_timestamp(t->fake_index, t->nb_ptses, sample_time, flags)) < 0)
                return sample_time;
        } else {
            /* no IndexEntryArray (one or more CBR segments)
             * make sure we don't seek past the end */
2342
            sample_time = FFMIN(sample_time, source_track->original_duration - 1);
2343 2344
        }

2345 2346 2347
        if ((ret = mxf_edit_unit_absolute_offset(mxf, t, sample_time, &sample_time, &seekpos, 1)) << 0)
            return ret;

2348
        ff_update_cur_dts(s, st, sample_time);
2349
        mxf->current_edit_unit = sample_time;
2350 2351
        avio_seek(s->pb, seekpos, SEEK_SET);
    }
2352 2353 2354 2355 2356 2357 2358 2359 2360 2361 2362 2363 2364 2365

    // Update all tracks sample count
    for (i = 0; i < s->nb_streams; i++) {
        AVStream *cur_st = s->streams[i];
        MXFTrack *cur_track = cur_st->priv_data;
        uint64_t current_sample_count = 0;
        if (cur_st->codec->codec_type == AVMEDIA_TYPE_AUDIO) {
            ret = mxf_compute_sample_count(mxf, i, &current_sample_count);
            if (ret < 0)
                return ret;

            cur_track->sample_count = current_sample_count;
        }
    }
2366 2367 2368
    return 0;
}

2369
AVInputFormat ff_mxf_demuxer = {
2370
    .name           = "mxf",
2371
    .long_name      = NULL_IF_CONFIG_SMALL("MXF (Material eXchange Format)"),
2372 2373 2374 2375 2376 2377
    .priv_data_size = sizeof(MXFContext),
    .read_probe     = mxf_probe,
    .read_header    = mxf_read_header,
    .read_packet    = mxf_read_packet,
    .read_close     = mxf_read_close,
    .read_seek      = mxf_read_seek,
2378
};