git.sesse.net Git - ffmpeg/blob - libavformat/mxfdec.c

   1 /*
   2  * MXF demuxer.
   3  * Copyright (c) 2006 SmartJog S.A., Baptiste Coudurier <baptiste dot coudurier at smartjog dot com>
   4  *
   5  * This file is part of Libav.
   6  *
   7  * Libav is free software; you can redistribute it and/or
   8  * modify it under the terms of the GNU Lesser General Public
   9  * License as published by the Free Software Foundation; either
  10  * version 2.1 of the License, or (at your option) any later version.
  11  *
  12  * Libav is distributed in the hope that it will be useful,
  13  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  14  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  15  * Lesser General Public License for more details.
  16  *
  17  * You should have received a copy of the GNU Lesser General Public
  18  * License along with Libav; if not, write to the Free Software
  19  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
  20  */
  21
  22 /*
  23  * References
  24  * SMPTE 336M KLV Data Encoding Protocol Using Key-Length-Value
  25  * SMPTE 377M MXF File Format Specifications
  26  * SMPTE 378M Operational Pattern 1a
  27  * SMPTE 379M MXF Generic Container
  28  * SMPTE 381M Mapping MPEG Streams into the MXF Generic Container
  29  * SMPTE 382M Mapping AES3 and Broadcast Wave Audio into the MXF Generic Container
  30  * SMPTE 383M Mapping DV-DIF Data to the MXF Generic Container
  31  *
  32  * Principle
  33  * Search for Track numbers which will identify essence element KLV packets.
  34  * Search for SourcePackage which define tracks which contains Track numbers.
  35  * Material Package contains tracks with reference to SourcePackage tracks.
  36  * Search for Descriptors (Picture, Sound) which contains codec info and parameters.
  37  * Assign Descriptors to correct Tracks.
  38  *
  39  * Metadata reading functions read Local Tags, get InstanceUID(0x3C0A) then add MetaDataSet to MXFContext.
  40  * Metadata parsing resolves Strong References to objects.
  41  *
  42  * Simple demuxer, only OP1A supported and some files might not work at all.
  43  * Only tracks with associated descriptors will be decoded. "Highly Desirable" SMPTE 377M D.1
  44  */
  45
  46 #include <stdint.h>
  47
  48 #include "libavutil/aes.h"
  49 #include "libavutil/mathematics.h"
  50 #include "libavcodec/bytestream.h"
  51 #include "avformat.h"
  52 #include "internal.h"
  53 #include "mxf.h"
  54
  55 typedef enum {
  56     Header,
  57     BodyPartition,
  58     Footer
  59 } MXFPartitionType;
  60
  61 typedef enum {
  62     OP1a = 1,
  63     OP1b,
  64     OP1c,
  65     OP2a,
  66     OP2b,
  67     OP2c,
  68     OP3a,
  69     OP3b,
  70     OP3c,
  71     OPAtom,
  72     OPSonyOpt,  /* FATE sample, violates the spec in places */
  73 } MXFOP;
  74
  75 typedef struct {
  76     int closed;
  77     int complete;
  78     MXFPartitionType type;
  79     uint64_t previous_partition;
  80     int index_sid;
  81     int body_sid;
  82     int64_t this_partition;
  83     int64_t essence_offset;         ///< absolute offset of essence
  84     int64_t essence_length;
  85     int32_t kag_size;
  86     int64_t header_byte_count;
  87     int64_t index_byte_count;
  88     int pack_length;
  89 } MXFPartition;
  90
  91 typedef struct {
  92     UID uid;
  93     enum MXFMetadataSetType type;
  94     UID source_container_ul;
  95 } MXFCryptoContext;
  96
  97 typedef struct {
  98     UID uid;
  99     enum MXFMetadataSetType type;
 100     UID source_package_uid;
 101     UID data_definition_ul;
 102     int64_t duration;
 103     int64_t start_position;
 104     int source_track_id;
 105 } MXFStructuralComponent;
 106
 107 typedef struct {
 108     UID uid;
 109     enum MXFMetadataSetType type;
 110     UID data_definition_ul;
 111     UID *structural_components_refs;
 112     int structural_components_count;
 113     int64_t duration;
 114 } MXFSequence;
 115
 116 typedef struct {
 117     UID uid;
 118     enum MXFMetadataSetType type;
 119     MXFSequence *sequence; /* mandatory, and only one */
 120     UID sequence_ref;
 121     int track_id;
 122     uint8_t track_number[4];
 123     AVRational edit_rate;
 124     int intra_only;
 125 } MXFTrack;
 126
 127 typedef struct {
 128     UID uid;
 129     enum MXFMetadataSetType type;
 130     UID essence_container_ul;
 131     UID essence_codec_ul;
 132     AVRational sample_rate;
 133     AVRational aspect_ratio;
 134     int width;
 135     int height; /* Field height, not frame height */
 136     int frame_layout; /* See MXFFrameLayout enum */
 137     int channels;
 138     int bits_per_sample;
 139     unsigned int component_depth;
 140     unsigned int horiz_subsampling;
 141     unsigned int vert_subsampling;
 142     UID *sub_descriptors_refs;
 143     int sub_descriptors_count;
 144     int linked_track_id;
 145     uint8_t *extradata;
 146     int extradata_size;
 147     enum AVPixelFormat pix_fmt;
 148 } MXFDescriptor;
 149
 150 typedef struct {
 151     UID uid;
 152     enum MXFMetadataSetType type;
 153     int edit_unit_byte_count;
 154     int index_sid;
 155     int body_sid;
 156     AVRational index_edit_rate;
 157     uint64_t index_start_position;
 158     uint64_t index_duration;
 159     int8_t *temporal_offset_entries;
 160     int *flag_entries;
 161     uint64_t *stream_offset_entries;
 162     int nb_index_entries;
 163 } MXFIndexTableSegment;
 164
 165 typedef struct {
 166     UID uid;
 167     enum MXFMetadataSetType type;
 168     UID package_uid;
 169     UID *tracks_refs;
 170     int tracks_count;
 171     MXFDescriptor *descriptor; /* only one */
 172     UID descriptor_ref;
 173 } MXFPackage;
 174
 175 typedef struct {
 176     UID uid;
 177     enum MXFMetadataSetType type;
 178 } MXFMetadataSet;
 179
 180 /* decoded index table */
 181 typedef struct {
 182     int index_sid;
 183     int body_sid;
 184     int nb_ptses;               /* number of PTSes or total duration of index */
 185     int64_t first_dts;          /* DTS = EditUnit + first_dts */
 186     int64_t *ptses;             /* maps EditUnit -> PTS */
 187     int nb_segments;
 188     MXFIndexTableSegment **segments;    /* sorted by IndexStartPosition */
 189     AVIndexEntry *fake_index;   /* used for calling ff_index_search_timestamp() */
 190 } MXFIndexTable;
 191
 192 typedef struct {
 193     MXFPartition *partitions;
 194     unsigned partitions_count;
 195     MXFOP op;
 196     UID *packages_refs;
 197     int packages_count;
 198     MXFMetadataSet **metadata_sets;
 199     int metadata_sets_count;
 200     AVFormatContext *fc;
 201     struct AVAES *aesc;
 202     uint8_t *local_tags;
 203     int local_tags_count;
 204     uint64_t footer_partition;
 205     KLVPacket current_klv_data;
 206     int current_klv_index;
 207     int run_in;
 208     MXFPartition *current_partition;
 209     int parsing_backward;
 210     int64_t last_forward_tell;
 211     int last_forward_partition;
 212     int current_edit_unit;
 213     int nb_index_tables;
 214     MXFIndexTable *index_tables;
 215     int edit_units_per_packet;      ///< how many edit units to read at a time (PCM, OPAtom)
 216 } MXFContext;
 217
 218 enum MXFWrappingScheme {
 219     Frame,
 220     Clip,
 221 };
 222
 223 /* NOTE: klv_offset is not set (-1) for local keys */
 224 typedef int MXFMetadataReadFunc(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset);
 225
 226 typedef struct {
 227     const UID key;
 228     MXFMetadataReadFunc *read;
 229     int ctx_size;
 230     enum MXFMetadataSetType type;
 231 } MXFMetadataReadTableEntry;
 232
 233 /* partial keys to match */
 234 static const uint8_t mxf_header_partition_pack_key[]       = { 0x06,0x0e,0x2b,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x02 };
 235 static const uint8_t mxf_essence_element_key[]             = { 0x06,0x0e,0x2b,0x34,0x01,0x02,0x01,0x01,0x0d,0x01,0x03,0x01 };
 236 static const uint8_t mxf_avid_essence_element_key[]        = { 0x06,0x0e,0x2b,0x34,0x01,0x02,0x01,0x01,0x0e,0x04,0x03,0x01 };
 237 static const uint8_t mxf_system_item_key[]                 = { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0D,0x01,0x03,0x01,0x04 };
 238 static const uint8_t mxf_klv_key[]                         = { 0x06,0x0e,0x2b,0x34 };
 239 /* complete keys to match */
 240 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 };
 241 static const uint8_t mxf_encrypted_triplet_key[]           = { 0x06,0x0e,0x2b,0x34,0x02,0x04,0x01,0x07,0x0d,0x01,0x03,0x01,0x02,0x7e,0x01,0x00 };
 242 static const uint8_t mxf_encrypted_essence_container[]     = { 0x06,0x0e,0x2b,0x34,0x04,0x01,0x01,0x07,0x0d,0x01,0x03,0x01,0x02,0x0b,0x01,0x00 };
 243 static const uint8_t mxf_sony_mpeg4_extradata[]            = { 0x06,0x0e,0x2b,0x34,0x04,0x01,0x01,0x01,0x0e,0x06,0x06,0x02,0x02,0x01,0x00,0x00 };
 244
 245 #define IS_KLV_KEY(x, y) (!memcmp(x, y, sizeof(y)))
 246
 247 static int64_t klv_decode_ber_length(AVIOContext *pb)
 248 {
 249     uint64_t size = avio_r8(pb);
 250     if (size & 0x80) { /* long form */
 251         int bytes_num = size & 0x7f;
 252         /* SMPTE 379M 5.3.4 guarantee that bytes_num must not exceed 8 bytes */
 253         if (bytes_num > 8)
 254             return AVERROR_INVALIDDATA;
 255         size = 0;
 256         while (bytes_num--)
 257             size = size << 8 | avio_r8(pb);
 258     }
 259     return size;
 260 }
 261
 262 static int mxf_read_sync(AVIOContext *pb, const uint8_t *key, unsigned size)
 263 {
 264     int i, b;
 265     for (i = 0; i < size && !pb->eof_reached; i++) {
 266         b = avio_r8(pb);
 267         if (b == key[0])
 268             i = 0;
 269         else if (b != key[i])
 270             i = -1;
 271     }
 272     return i == size;
 273 }
 274
 275 static int klv_read_packet(KLVPacket *klv, AVIOContext *pb)
 276 {
 277     if (!mxf_read_sync(pb, mxf_klv_key, 4))
 278         return AVERROR_INVALIDDATA;
 279     klv->offset = avio_tell(pb) - 4;
 280     memcpy(klv->key, mxf_klv_key, 4);
 281     avio_read(pb, klv->key + 4, 12);
 282     klv->length = klv_decode_ber_length(pb);
 283     return klv->length == -1 ? -1 : 0;
 284 }
 285
 286 static int mxf_get_stream_index(AVFormatContext *s, KLVPacket *klv)
 287 {
 288     int i;
 289
 290     for (i = 0; i < s->nb_streams; i++) {
 291         MXFTrack *track = s->streams[i]->priv_data;
 292         /* SMPTE 379M 7.3 */
 293         if (!memcmp(klv->key + sizeof(mxf_essence_element_key), track->track_number, sizeof(track->track_number)))
 294             return i;
 295     }
 296     /* return 0 if only one stream, for OP Atom files with 0 as track number */
 297     return s->nb_streams == 1 ? 0 : -1;
 298 }
 299
 300 /* XXX: use AVBitStreamFilter */
 301 static int mxf_get_d10_aes3_packet(AVIOContext *pb, AVStream *st, AVPacket *pkt, int64_t length)
 302 {
 303     const uint8_t *buf_ptr, *end_ptr;
 304     uint8_t *data_ptr;
 305     int i;
 306
 307     if (length > 61444) /* worst case PAL 1920 samples 8 channels */
 308         return AVERROR_INVALIDDATA;
 309     length = av_get_packet(pb, pkt, length);
 310     if (length < 0)
 311         return length;
 312     data_ptr = pkt->data;
 313     end_ptr = pkt->data + length;
 314     buf_ptr = pkt->data + 4; /* skip SMPTE 331M header */
 315     for (; end_ptr - buf_ptr >= st->codec->channels * 4; ) {
 316         for (i = 0; i < st->codec->channels; i++) {
 317             uint32_t sample = bytestream_get_le32(&buf_ptr);
 318             if (st->codec->bits_per_coded_sample == 24)
 319                 bytestream_put_le24(&data_ptr, (sample >> 4) & 0xffffff);
 320             else
 321                 bytestream_put_le16(&data_ptr, (sample >> 12) & 0xffff);
 322         }
 323         buf_ptr += 32 - st->codec->channels*4; // always 8 channels stored SMPTE 331M
 324     }
 325     av_shrink_packet(pkt, data_ptr - pkt->data);
 326     return 0;
 327 }
 328
 329 static int mxf_decrypt_triplet(AVFormatContext *s, AVPacket *pkt, KLVPacket *klv)
 330 {
 331     static const uint8_t checkv[16] = {0x43, 0x48, 0x55, 0x4b, 0x43, 0x48, 0x55, 0x4b, 0x43, 0x48, 0x55, 0x4b, 0x43, 0x48, 0x55, 0x4b};
 332     MXFContext *mxf = s->priv_data;
 333     AVIOContext *pb = s->pb;
 334     int64_t end = avio_tell(pb) + klv->length;
 335     int64_t size;
 336     uint64_t orig_size;
 337     uint64_t plaintext_size;
 338     uint8_t ivec[16];
 339     uint8_t tmpbuf[16];
 340     int index;
 341
 342     if (!mxf->aesc && s->key && s->keylen == 16) {
 343         mxf->aesc = av_aes_alloc();
 344         if (!mxf->aesc)
 345             return AVERROR(ENOMEM);
 346         av_aes_init(mxf->aesc, s->key, 128, 1);
 347     }
 348     // crypto context
 349     avio_skip(pb, klv_decode_ber_length(pb));
 350     // plaintext offset
 351     klv_decode_ber_length(pb);
 352     plaintext_size = avio_rb64(pb);
 353     // source klv key
 354     klv_decode_ber_length(pb);
 355     avio_read(pb, klv->key, 16);
 356     if (!IS_KLV_KEY(klv, mxf_essence_element_key))
 357         return AVERROR_INVALIDDATA;
 358     index = mxf_get_stream_index(s, klv);
 359     if (index < 0)
 360         return AVERROR_INVALIDDATA;
 361     // source size
 362     klv_decode_ber_length(pb);
 363     orig_size = avio_rb64(pb);
 364     if (orig_size < plaintext_size)
 365         return AVERROR_INVALIDDATA;
 366     // enc. code
 367     size = klv_decode_ber_length(pb);
 368     if (size < 32 || size - 32 < orig_size)
 369         return AVERROR_INVALIDDATA;
 370     avio_read(pb, ivec, 16);
 371     avio_read(pb, tmpbuf, 16);
 372     if (mxf->aesc)
 373         av_aes_crypt(mxf->aesc, tmpbuf, tmpbuf, 1, ivec, 1);
 374     if (memcmp(tmpbuf, checkv, 16))
 375         av_log(s, AV_LOG_ERROR, "probably incorrect decryption key\n");
 376     size -= 32;
 377     size = av_get_packet(pb, pkt, size);
 378     if (size < 0)
 379         return size;
 380     else if (size < plaintext_size)
 381         return AVERROR_INVALIDDATA;
 382     size -= plaintext_size;
 383     if (mxf->aesc)
 384         av_aes_crypt(mxf->aesc, &pkt->data[plaintext_size],
 385                      &pkt->data[plaintext_size], size >> 4, ivec, 1);
 386     av_shrink_packet(pkt, orig_size);
 387     pkt->stream_index = index;
 388     avio_skip(pb, end - avio_tell(pb));
 389     return 0;
 390 }
 391
 392 static int mxf_read_primer_pack(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
 393 {
 394     MXFContext *mxf = arg;
 395     int item_num = avio_rb32(pb);
 396     int item_len = avio_rb32(pb);
 397
 398     if (item_len != 18) {
 399         avpriv_request_sample(pb, "Primer pack item length %d", item_len);
 400         return AVERROR_PATCHWELCOME;
 401     }
 402     if (item_num > UINT_MAX / item_len)
 403         return AVERROR_INVALIDDATA;
 404     mxf->local_tags_count = item_num;
 405     mxf->local_tags = av_malloc(item_num*item_len);
 406     if (!mxf->local_tags)
 407         return AVERROR(ENOMEM);
 408     avio_read(pb, mxf->local_tags, item_num*item_len);
 409     return 0;
 410 }
 411
 412 static int mxf_read_partition_pack(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
 413 {
 414     MXFContext *mxf = arg;
 415     MXFPartition *partition;
 416     UID op;
 417     uint64_t footer_partition;
 418     uint32_t nb_essence_containers;
 419     int err;
 420
 421     if ((err = av_reallocp_array(&mxf->partitions, mxf->partitions_count + 1,
 422                                  sizeof(*mxf->partitions))) < 0) {
 423         mxf->partitions_count = 0;
 424         return err;
 425     }
 426
 427     if (mxf->parsing_backward) {
 428         /* insert the new partition pack in the middle
 429          * this makes the entries in mxf->partitions sorted by offset */
 430         memmove(&mxf->partitions[mxf->last_forward_partition+1],
 431                 &mxf->partitions[mxf->last_forward_partition],
 432                 (mxf->partitions_count - mxf->last_forward_partition)*sizeof(*mxf->partitions));
 433         partition = mxf->current_partition = &mxf->partitions[mxf->last_forward_partition];
 434     } else {
 435         mxf->last_forward_partition++;
 436         partition = mxf->current_partition = &mxf->partitions[mxf->partitions_count];
 437     }
 438
 439     memset(partition, 0, sizeof(*partition));
 440     mxf->partitions_count++;
 441     partition->pack_length = avio_tell(pb) - klv_offset + size;
 442
 443     switch(uid[13]) {
 444     case 2:
 445         partition->type = Header;
 446         break;
 447     case 3:
 448         partition->type = BodyPartition;
 449         break;
 450     case 4:
 451         partition->type = Footer;
 452         break;
 453     default:
 454         av_log(mxf->fc, AV_LOG_ERROR, "unknown partition type %i\n", uid[13]);
 455         return AVERROR_INVALIDDATA;
 456     }
 457
 458     /* consider both footers to be closed (there is only Footer and CompleteFooter) */
 459     partition->closed = partition->type == Footer || !(uid[14] & 1);
 460     partition->complete = uid[14] > 2;
 461     avio_skip(pb, 4);
 462     partition->kag_size = avio_rb32(pb);
 463     partition->this_partition = avio_rb64(pb);
 464     partition->previous_partition = avio_rb64(pb);
 465     footer_partition = avio_rb64(pb);
 466     partition->header_byte_count = avio_rb64(pb);
 467     partition->index_byte_count = avio_rb64(pb);
 468     partition->index_sid = avio_rb32(pb);
 469     avio_skip(pb, 8);
 470     partition->body_sid = avio_rb32(pb);
 471     avio_read(pb, op, sizeof(UID));
 472     nb_essence_containers = avio_rb32(pb);
 473
 474     /* some files don'thave FooterPartition set in every partition */
 475     if (footer_partition) {
 476         if (mxf->footer_partition && mxf->footer_partition != footer_partition) {
 477             av_log(mxf->fc, AV_LOG_ERROR,
 478                    "inconsistent FooterPartition value: %"PRIu64" != %"PRIu64"\n",
 479                    mxf->footer_partition, footer_partition);
 480         } else {
 481             mxf->footer_partition = footer_partition;
 482         }
 483     }
 484
 485     av_dlog(mxf->fc,
 486             "PartitionPack: ThisPartition = 0x%"PRIX64
 487             ", PreviousPartition = 0x%"PRIX64", "
 488             "FooterPartition = 0x%"PRIX64", IndexSID = %i, BodySID = %i\n",
 489             partition->this_partition,
 490             partition->previous_partition, footer_partition,
 491             partition->index_sid, partition->body_sid);
 492
 493     /* sanity check PreviousPartition if set */
 494     if (partition->previous_partition &&
 495         mxf->run_in + partition->previous_partition >= klv_offset) {
 496         av_log(mxf->fc, AV_LOG_ERROR,
 497                "PreviousPartition points to this partition or forward\n");
 498         return AVERROR_INVALIDDATA;
 499     }
 500
 501     if      (op[12] == 1 && op[13] == 1) mxf->op = OP1a;
 502     else if (op[12] == 1 && op[13] == 2) mxf->op = OP1b;
 503     else if (op[12] == 1 && op[13] == 3) mxf->op = OP1c;
 504     else if (op[12] == 2 && op[13] == 1) mxf->op = OP2a;
 505     else if (op[12] == 2 && op[13] == 2) mxf->op = OP2b;
 506     else if (op[12] == 2 && op[13] == 3) mxf->op = OP2c;
 507     else if (op[12] == 3 && op[13] == 1) mxf->op = OP3a;
 508     else if (op[12] == 3 && op[13] == 2) mxf->op = OP3b;
 509     else if (op[12] == 3 && op[13] == 3) mxf->op = OP3c;
 510     else if (op[12] == 64&& op[13] == 1) mxf->op = OPSonyOpt;
 511     else if (op[12] == 0x10) {
 512         /* SMPTE 390m: "There shall be exactly one essence container"
 513          * The following block deals with files that violate this, namely:
 514          * 2011_DCPTEST_24FPS.V.mxf - two ECs, OP1a
 515          * abcdefghiv016f56415e.mxf - zero ECs, OPAtom, output by Avid AirSpeed */
 516         if (nb_essence_containers != 1) {
 517             MXFOP op = nb_essence_containers ? OP1a : OPAtom;
 518
 519             /* only nag once */
 520             if (!mxf->op)
 521                 av_log(mxf->fc, AV_LOG_WARNING,
 522                        "\"OPAtom\" with %u ECs - assuming %s\n",
 523                        nb_essence_containers,
 524                        op == OP1a ? "OP1a" : "OPAtom");
 525
 526             mxf->op = op;
 527         } else
 528             mxf->op = OPAtom;
 529     } else {
 530         av_log(mxf->fc, AV_LOG_ERROR, "unknown operational pattern: %02xh %02xh - guessing OP1a\n", op[12], op[13]);
 531         mxf->op = OP1a;
 532     }
 533
 534     if (partition->kag_size <= 0 || partition->kag_size > (1 << 20)) {
 535         av_log(mxf->fc, AV_LOG_WARNING, "invalid KAGSize %i - guessing ", partition->kag_size);
 536
 537         if (mxf->op == OPSonyOpt)
 538             partition->kag_size = 512;
 539         else
 540             partition->kag_size = 1;
 541
 542         av_log(mxf->fc, AV_LOG_WARNING, "%i\n", partition->kag_size);
 543     }
 544
 545     return 0;
 546 }
 547
 548 static int mxf_add_metadata_set(MXFContext *mxf, void *metadata_set)
 549 {
 550     int err;
 551
 552     if ((err = av_reallocp_array(&mxf->metadata_sets, mxf->metadata_sets_count + 1,
 553                                  sizeof(*mxf->metadata_sets))) < 0) {
 554         mxf->metadata_sets_count = 0;
 555         return err;
 556     }
 557     mxf->metadata_sets[mxf->metadata_sets_count] = metadata_set;
 558     mxf->metadata_sets_count++;
 559     return 0;
 560 }
 561
 562 static int mxf_read_cryptographic_context(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
 563 {
 564     MXFCryptoContext *cryptocontext = arg;
 565     if (size != 16)
 566         return AVERROR_INVALIDDATA;
 567     if (IS_KLV_KEY(uid, mxf_crypto_source_container_ul))
 568         avio_read(pb, cryptocontext->source_container_ul, 16);
 569     return 0;
 570 }
 571
 572 static int mxf_read_content_storage(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
 573 {
 574     MXFContext *mxf = arg;
 575     switch (tag) {
 576     case 0x1901:
 577         mxf->packages_count = avio_rb32(pb);
 578         if (mxf->packages_count >= UINT_MAX / sizeof(UID))
 579             return AVERROR_INVALIDDATA;
 580         mxf->packages_refs = av_malloc(mxf->packages_count * sizeof(UID));
 581         if (!mxf->packages_refs)
 582             return AVERROR(ENOMEM);
 583         avio_skip(pb, 4); /* useless size of objects, always 16 according to specs */
 584         avio_read(pb, (uint8_t *)mxf->packages_refs, mxf->packages_count * sizeof(UID));
 585         break;
 586     }
 587     return 0;
 588 }
 589
 590 static int mxf_read_source_clip(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
 591 {
 592     MXFStructuralComponent *source_clip = arg;
 593     switch(tag) {
 594     case 0x0202:
 595         source_clip->duration = avio_rb64(pb);
 596         break;
 597     case 0x1201:
 598         source_clip->start_position = avio_rb64(pb);
 599         break;
 600     case 0x1101:
 601         /* UMID, only get last 16 bytes */
 602         avio_skip(pb, 16);
 603         avio_read(pb, source_clip->source_package_uid, 16);
 604         break;
 605     case 0x1102:
 606         source_clip->source_track_id = avio_rb32(pb);
 607         break;
 608     }
 609     return 0;
 610 }
 611
 612 static int mxf_read_material_package(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
 613 {
 614     MXFPackage *package = arg;
 615     switch(tag) {
 616     case 0x4403:
 617         package->tracks_count = avio_rb32(pb);
 618         if (package->tracks_count >= UINT_MAX / sizeof(UID))
 619             return AVERROR_INVALIDDATA;
 620         package->tracks_refs = av_malloc(package->tracks_count * sizeof(UID));
 621         if (!package->tracks_refs)
 622             return AVERROR(ENOMEM);
 623         avio_skip(pb, 4); /* useless size of objects, always 16 according to specs */
 624         avio_read(pb, (uint8_t *)package->tracks_refs, package->tracks_count * sizeof(UID));
 625         break;
 626     }
 627     return 0;
 628 }
 629
 630 static int mxf_read_track(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
 631 {
 632     MXFTrack *track = arg;
 633     switch(tag) {
 634     case 0x4801:
 635         track->track_id = avio_rb32(pb);
 636         break;
 637     case 0x4804:
 638         avio_read(pb, track->track_number, 4);
 639         break;
 640     case 0x4B01:
 641         track->edit_rate.num = avio_rb32(pb);
 642         track->edit_rate.den = avio_rb32(pb);
 643         break;
 644     case 0x4803:
 645         avio_read(pb, track->sequence_ref, 16);
 646         break;
 647     }
 648     return 0;
 649 }
 650
 651 static int mxf_read_sequence(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
 652 {
 653     MXFSequence *sequence = arg;
 654     switch(tag) {
 655     case 0x0202:
 656         sequence->duration = avio_rb64(pb);
 657         break;
 658     case 0x0201:
 659         avio_read(pb, sequence->data_definition_ul, 16);
 660         break;
 661     case 0x1001:
 662         sequence->structural_components_count = avio_rb32(pb);
 663         if (sequence->structural_components_count >= UINT_MAX / sizeof(UID))
 664             return AVERROR_INVALIDDATA;
 665         sequence->structural_components_refs = av_malloc(sequence->structural_components_count * sizeof(UID));
 666         if (!sequence->structural_components_refs)
 667             return AVERROR(ENOMEM);
 668         avio_skip(pb, 4); /* useless size of objects, always 16 according to specs */
 669         avio_read(pb, (uint8_t *)sequence->structural_components_refs, sequence->structural_components_count * sizeof(UID));
 670         break;
 671     }
 672     return 0;
 673 }
 674
 675 static int mxf_read_source_package(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
 676 {
 677     MXFPackage *package = arg;
 678     switch(tag) {
 679     case 0x4403:
 680         package->tracks_count = avio_rb32(pb);
 681         if (package->tracks_count >= UINT_MAX / sizeof(UID))
 682             return AVERROR_INVALIDDATA;
 683         package->tracks_refs = av_malloc(package->tracks_count * sizeof(UID));
 684         if (!package->tracks_refs)
 685             return AVERROR(ENOMEM);
 686         avio_skip(pb, 4); /* useless size of objects, always 16 according to specs */
 687         avio_read(pb, (uint8_t *)package->tracks_refs, package->tracks_count * sizeof(UID));
 688         break;
 689     case 0x4401:
 690         /* UMID, only get last 16 bytes */
 691         avio_skip(pb, 16);
 692         avio_read(pb, package->package_uid, 16);
 693         break;
 694     case 0x4701:
 695         avio_read(pb, package->descriptor_ref, 16);
 696         break;
 697     }
 698     return 0;
 699 }
 700
 701 static int mxf_read_index_entry_array(AVIOContext *pb, MXFIndexTableSegment *segment)
 702 {
 703     int i, length;
 704
 705     segment->nb_index_entries = avio_rb32(pb);
 706     if (!segment->nb_index_entries)
 707         return 0;
 708     else if (segment->nb_index_entries < 0 ||
 709              segment->nb_index_entries >
 710              (INT_MAX / sizeof(*segment->stream_offset_entries)))
 711         return AVERROR(ENOMEM);
 712
 713     length = avio_rb32(pb);
 714
 715     segment->temporal_offset_entries = av_mallocz(segment->nb_index_entries *
 716                                  sizeof(*segment->temporal_offset_entries));
 717     segment->flag_entries            = av_mallocz(segment->nb_index_entries *
 718                                  sizeof(*segment->flag_entries));
 719     segment->stream_offset_entries   = av_mallocz(segment->nb_index_entries *
 720                                  sizeof(*segment->stream_offset_entries));
 721
 722     if (!segment->flag_entries || !segment->stream_offset_entries ||
 723         !segment->temporal_offset_entries) {
 724         av_freep(&segment->flag_entries);
 725         av_freep(&segment->stream_offset_entries);
 726         av_freep(&segment->temporal_offset_entries);
 727         return AVERROR(ENOMEM);
 728     }
 729
 730     for (i = 0; i < segment->nb_index_entries; i++) {
 731         segment->temporal_offset_entries[i] = avio_r8(pb);
 732         avio_r8(pb);                                        /* KeyFrameOffset */
 733         segment->flag_entries[i] = avio_r8(pb);
 734         segment->stream_offset_entries[i] = avio_rb64(pb);
 735         avio_skip(pb, length - 11);
 736     }
 737     return 0;
 738 }
 739
 740 static int mxf_read_index_table_segment(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
 741 {
 742     MXFIndexTableSegment *segment = arg;
 743     switch(tag) {
 744     case 0x3F05:
 745         segment->edit_unit_byte_count = avio_rb32(pb);
 746         av_dlog(NULL, "EditUnitByteCount %d\n", segment->edit_unit_byte_count);
 747         break;
 748     case 0x3F06:
 749         segment->index_sid = avio_rb32(pb);
 750         av_dlog(NULL, "IndexSID %d\n", segment->index_sid);
 751         break;
 752     case 0x3F07:
 753         segment->body_sid = avio_rb32(pb);
 754         av_dlog(NULL, "BodySID %d\n", segment->body_sid);
 755         break;
 756     case 0x3F0A:
 757         av_dlog(NULL, "IndexEntryArray found\n");
 758         return mxf_read_index_entry_array(pb, segment);
 759     case 0x3F0B:
 760         segment->index_edit_rate.num = avio_rb32(pb);
 761         segment->index_edit_rate.den = avio_rb32(pb);
 762         av_dlog(NULL, "IndexEditRate %d/%d\n", segment->index_edit_rate.num,
 763                 segment->index_edit_rate.den);
 764         break;
 765     case 0x3F0C:
 766         segment->index_start_position = avio_rb64(pb);
 767         av_dlog(NULL, "IndexStartPosition %"PRId64"\n", segment->index_start_position);
 768         break;
 769     case 0x3F0D:
 770         segment->index_duration = avio_rb64(pb);
 771         av_dlog(NULL, "IndexDuration %"PRId64"\n", segment->index_duration);
 772         break;
 773     }
 774     return 0;
 775 }
 776
 777 static void mxf_read_pixel_layout(AVIOContext *pb, MXFDescriptor *descriptor)
 778 {
 779     int code, value, ofs = 0;
 780     char layout[16] = {0};
 781
 782     do {
 783         code = avio_r8(pb);
 784         value = avio_r8(pb);
 785         av_dlog(NULL, "pixel layout: code %#x\n", code);
 786
 787         if (ofs < 16) {
 788             layout[ofs++] = code;
 789             layout[ofs++] = value;
 790         }
 791     } while (code != 0); /* SMPTE 377M E.2.46 */
 792
 793     ff_mxf_decode_pixel_layout(layout, &descriptor->pix_fmt);
 794 }
 795
 796 static int mxf_read_generic_descriptor(void *arg, AVIOContext *pb, int tag, int size, UID uid, int64_t klv_offset)
 797 {
 798     MXFDescriptor *descriptor = arg;
 799     descriptor->pix_fmt = AV_PIX_FMT_NONE;
 800     switch(tag) {
 801     case 0x3F01:
 802         descriptor->sub_descriptors_count = avio_rb32(pb);
 803         if (descriptor->sub_descriptors_count >= UINT_MAX / sizeof(UID))
 804             return AVERROR_INVALIDDATA;
 805         descriptor->sub_descriptors_refs = av_malloc(descriptor->sub_descriptors_count * sizeof(UID));
 806         if (!descriptor->sub_descriptors_refs)
 807             return AVERROR(ENOMEM);
 808         avio_skip(pb, 4); /* useless size of objects, always 16 according to specs */
 809         avio_read(pb, (uint8_t *)descriptor->sub_descriptors_refs, descriptor->sub_descriptors_count * sizeof(UID));
 810         break;
 811     case 0x3004:
 812         avio_read(pb, descriptor->essence_container_ul, 16);
 813         break;
 814     case 0x3006:
 815         descriptor->linked_track_id = avio_rb32(pb);
 816         break;
 817     case 0x3201: /* PictureEssenceCoding */
 818         avio_read(pb, descriptor->essence_codec_ul, 16);
 819         break;
 820     case 0x3203:
 821         descriptor->width = avio_rb32(pb);
 822         break;
 823     case 0x3202:
 824         descriptor->height = avio_rb32(pb);
 825         break;
 826     case 0x320C:
 827         descriptor->frame_layout = avio_r8(pb);
 828         break;
 829     case 0x320E:
 830         descriptor->aspect_ratio.num = avio_rb32(pb);
 831         descriptor->aspect_ratio.den = avio_rb32(pb);
 832         break;
 833     case 0x3301:
 834         descriptor->component_depth = avio_rb32(pb);
 835         break;
 836     case 0x3302:
 837         descriptor->horiz_subsampling = avio_rb32(pb);
 838         break;
 839     case 0x3308:
 840         descriptor->vert_subsampling = avio_rb32(pb);
 841         break;
 842     case 0x3D03:
 843         descriptor->sample_rate.num = avio_rb32(pb);
 844         descriptor->sample_rate.den = avio_rb32(pb);
 845         break;
 846     case 0x3D06: /* SoundEssenceCompression */
 847         avio_read(pb, descriptor->essence_codec_ul, 16);
 848         break;
 849     case 0x3D07:
 850         descriptor->channels = avio_rb32(pb);
 851         break;
 852     case 0x3D01:
 853         descriptor->bits_per_sample = avio_rb32(pb);
 854         break;
 855     case 0x3401:
 856         mxf_read_pixel_layout(pb, descriptor);
 857         break;
 858     default:
 859         /* Private uid used by SONY C0023S01.mxf */
 860         if (IS_KLV_KEY(uid, mxf_sony_mpeg4_extradata)) {
 861             av_free(descriptor->extradata);
 862             descriptor->extradata_size = 0;
 863             descriptor->extradata = av_malloc(size + FF_INPUT_BUFFER_PADDING_SIZE);
 864             if (!descriptor->extradata)
 865                 return AVERROR(ENOMEM);
 866             descriptor->extradata_size = size;
 867             avio_read(pb, descriptor->extradata, size);
 868         }
 869         break;
 870     }
 871     return 0;
 872 }
 873
 874 /*
 875  * Match an uid independently of the version byte and up to len common bytes
 876  * Returns: boolean
 877  */
 878 static int mxf_match_uid(const UID key, const UID uid, int len)
 879 {
 880     int i;
 881     for (i = 0; i < len; i++) {
 882         if (i != 7 && key[i] != uid[i])
 883             return 0;
 884     }
 885     return 1;
 886 }
 887
 888 static const MXFCodecUL *mxf_get_codec_ul(const MXFCodecUL *uls, UID *uid)
 889 {
 890     while (uls->uid[0]) {
 891         if(mxf_match_uid(uls->uid, *uid, uls->matching_len))
 892             break;
 893         uls++;
 894     }
 895     return uls;
 896 }
 897
 898 static void *mxf_resolve_strong_ref(MXFContext *mxf, UID *strong_ref, enum MXFMetadataSetType type)
 899 {
 900     int i;
 901
 902     if (!strong_ref)
 903         return NULL;
 904     for (i = 0; i < mxf->metadata_sets_count; i++) {
 905         if (!memcmp(*strong_ref, mxf->metadata_sets[i]->uid, 16) &&
 906             (type == AnyType || mxf->metadata_sets[i]->type == type)) {
 907             return mxf->metadata_sets[i];
 908         }
 909     }
 910     return NULL;
 911 }
 912
 913 static const MXFCodecUL mxf_picture_essence_container_uls[] = {
 914     // video essence container uls
 915     { { 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 */
 916     { { 0x06,0x0E,0x2B,0x34,0x04,0x01,0x01,0x01,0x0D,0x01,0x03,0x01,0x02,0x02,0x41,0x01 }, 14,    AV_CODEC_ID_DVVIDEO }, /* DV 625 25mbps */
 917     { { 0x06,0x0E,0x2B,0x34,0x04,0x01,0x01,0x01,0x0D,0x01,0x03,0x01,0x02,0x05,0x00,0x00 }, 14,   AV_CODEC_ID_RAWVIDEO }, /* Uncompressed Picture */
 918     { { 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00 },  0,      AV_CODEC_ID_NONE },
 919 };
 920
 921 /* EC ULs for intra-only formats */
 922 static const MXFCodecUL mxf_intra_only_essence_container_uls[] = {
 923     { { 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 */
 924     { { 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00 },  0,       AV_CODEC_ID_NONE },
 925 };
 926
 927 /* intra-only PictureEssenceCoding ULs, where no corresponding EC UL exists */
 928 static const MXFCodecUL mxf_intra_only_picture_essence_coding_uls[] = {
 929     { { 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 */
 930     { { 0x06,0x0E,0x2B,0x34,0x04,0x01,0x01,0x07,0x04,0x01,0x02,0x02,0x03,0x01,0x01,0x00 }, 14,   AV_CODEC_ID_JPEG2000 }, /* JPEG2000 Codestream */
 931     { { 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00 },  0,       AV_CODEC_ID_NONE },
 932 };
 933
 934 static const MXFCodecUL mxf_sound_essence_container_uls[] = {
 935     // sound essence container uls
 936     { { 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 */
 937     { { 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 */
 938     { { 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 */
 939     { { 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 */
 940     { { 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00 },  0,      AV_CODEC_ID_NONE },
 941 };
 942
 943 static int mxf_get_sorted_table_segments(MXFContext *mxf, int *nb_sorted_segments, MXFIndexTableSegment ***sorted_segments)
 944 {
 945     int i, j, nb_segments = 0;
 946     MXFIndexTableSegment **unsorted_segments;
 947     int last_body_sid = -1, last_index_sid = -1, last_index_start = -1;
 948
 949     /* count number of segments, allocate arrays and copy unsorted segments */
 950     for (i = 0; i < mxf->metadata_sets_count; i++)
 951         if (mxf->metadata_sets[i]->type == IndexTableSegment)
 952             nb_segments++;
 953
 954     if (!nb_segments)
 955         return AVERROR_INVALIDDATA;
 956
 957     *sorted_segments  = av_mallocz(nb_segments * sizeof(**sorted_segments));
 958     unsorted_segments = av_mallocz(nb_segments * sizeof(*unsorted_segments));
 959     if (!*sorted_segments || !unsorted_segments) {
 960         av_freep(sorted_segments);
 961         av_free(unsorted_segments);
 962         return AVERROR(ENOMEM);
 963     }
 964
 965     for (i = j = 0; i < mxf->metadata_sets_count; i++)
 966         if (mxf->metadata_sets[i]->type == IndexTableSegment)
 967             unsorted_segments[j++] = (MXFIndexTableSegment*)mxf->metadata_sets[i];
 968
 969     *nb_sorted_segments = 0;
 970
 971     /* sort segments by {BodySID, IndexSID, IndexStartPosition}, remove duplicates while we're at it */
 972     for (i = 0; i < nb_segments; i++) {
 973         int best = -1, best_body_sid = -1, best_index_sid = -1, best_index_start = -1;
 974         uint64_t best_index_duration = 0;
 975
 976         for (j = 0; j < nb_segments; j++) {
 977             MXFIndexTableSegment *s = unsorted_segments[j];
 978
 979             /* Require larger BosySID, IndexSID or IndexStartPosition then the previous entry. This removes duplicates.
 980              * We want the smallest values for the keys than what we currently have, unless this is the first such entry this time around.
 981              * If we come across an entry with the same IndexStartPosition but larger IndexDuration, then we'll prefer it over the one we currently have.
 982              */
 983             if ((i == 0     || s->body_sid > last_body_sid || s->index_sid > last_index_sid || s->index_start_position > last_index_start) &&
 984                 (best == -1 || s->body_sid < best_body_sid || s->index_sid < best_index_sid || s->index_start_position < best_index_start ||
 985                 (s->index_start_position == best_index_start && s->index_duration > best_index_duration))) {
 986                 best             = j;
 987                 best_body_sid    = s->body_sid;
 988                 best_index_sid   = s->index_sid;
 989                 best_index_start = s->index_start_position;
 990                 best_index_duration = s->index_duration;
 991             }
 992         }
 993
 994         /* no suitable entry found -> we're done */
 995         if (best == -1)
 996             break;
 997
 998         (*sorted_segments)[(*nb_sorted_segments)++] = unsorted_segments[best];
 999         last_body_sid    = best_body_sid;
1000         last_index_sid   = best_index_sid;
1001         last_index_start = best_index_start;
1002     }
1003
1004     av_free(unsorted_segments);
1005
1006     return 0;
1007 }
1008
1009 /**
1010  * Computes the absolute file offset of the given essence container offset
1011  */
1012 static int mxf_absolute_bodysid_offset(MXFContext *mxf, int body_sid, int64_t offset, int64_t *offset_out)
1013 {
1014     int x;
1015     int64_t offset_in = offset;     /* for logging */
1016
1017     for (x = 0; x < mxf->partitions_count; x++) {
1018         MXFPartition *p = &mxf->partitions[x];
1019
1020         if (p->body_sid != body_sid)
1021             continue;
1022
1023         if (offset < p->essence_length || !p->essence_length) {
1024             *offset_out = p->essence_offset + offset;
1025             return 0;
1026         }
1027
1028         offset -= p->essence_length;
1029     }
1030
1031     av_log(mxf->fc, AV_LOG_ERROR,
1032            "failed to find absolute offset of %"PRIX64" in BodySID %i - partial file?\n",
1033            offset_in, body_sid);
1034
1035     return AVERROR_INVALIDDATA;
1036 }
1037
1038 /**
1039  * Returns the end position of the essence container with given BodySID, or zero if unknown
1040  */
1041 static int64_t mxf_essence_container_end(MXFContext *mxf, int body_sid)
1042 {
1043     int x;
1044     int64_t ret = 0;
1045
1046     for (x = 0; x < mxf->partitions_count; x++) {
1047         MXFPartition *p = &mxf->partitions[x];
1048
1049         if (p->body_sid != body_sid)
1050             continue;
1051
1052         if (!p->essence_length)
1053             return 0;
1054
1055         ret = p->essence_offset + p->essence_length;
1056     }
1057
1058     return ret;
1059 }
1060
1061 /* EditUnit -> absolute offset */
1062 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)
1063 {
1064     int i;
1065     int64_t offset_temp = 0;
1066
1067     for (i = 0; i < index_table->nb_segments; i++) {
1068         MXFIndexTableSegment *s = index_table->segments[i];
1069
1070         edit_unit = FFMAX(edit_unit, s->index_start_position);  /* clamp if trying to seek before start */
1071
1072         if (edit_unit < s->index_start_position + s->index_duration) {
1073             int64_t index = edit_unit - s->index_start_position;
1074
1075             if (s->edit_unit_byte_count)
1076                 offset_temp += s->edit_unit_byte_count * index;
1077             else if (s->nb_index_entries) {
1078                 if (s->nb_index_entries == 2 * s->index_duration + 1)
1079                     index *= 2;     /* Avid index */
1080
1081                 if (index < 0 || index >= s->nb_index_entries) {
1082                     av_log(mxf->fc, AV_LOG_ERROR, "IndexSID %i segment at %"PRId64" IndexEntryArray too small\n",
1083                            index_table->index_sid, s->index_start_position);
1084                     return AVERROR_INVALIDDATA;
1085                 }
1086
1087                 offset_temp = s->stream_offset_entries[index];
1088             } else {
1089                 av_log(mxf->fc, AV_LOG_ERROR, "IndexSID %i segment at %"PRId64" missing EditUnitByteCount and IndexEntryArray\n",
1090                        index_table->index_sid, s->index_start_position);
1091                 return AVERROR_INVALIDDATA;
1092             }
1093
1094             if (edit_unit_out)
1095                 *edit_unit_out = edit_unit;
1096
1097             return mxf_absolute_bodysid_offset(mxf, index_table->body_sid, offset_temp, offset_out);
1098         } else {
1099             /* EditUnitByteCount == 0 for VBR indexes, which is fine since they use explicit StreamOffsets */
1100             offset_temp += s->edit_unit_byte_count * s->index_duration;
1101         }
1102     }
1103
1104     if (nag)
1105         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);
1106
1107     return AVERROR_INVALIDDATA;
1108 }
1109
1110 static int mxf_compute_ptses_fake_index(MXFContext *mxf, MXFIndexTable *index_table)
1111 {
1112     int i, j, x;
1113     int8_t max_temporal_offset = -128;
1114
1115     /* first compute how many entries we have */
1116     for (i = 0; i < index_table->nb_segments; i++) {
1117         MXFIndexTableSegment *s = index_table->segments[i];
1118
1119         if (!s->nb_index_entries) {
1120             index_table->nb_ptses = 0;
1121             return 0;                               /* no TemporalOffsets */
1122         }
1123
1124         index_table->nb_ptses += s->index_duration;
1125     }
1126
1127     /* paranoid check */
1128     if (index_table->nb_ptses <= 0)
1129         return 0;
1130
1131     if (index_table->nb_ptses > INT_MAX / sizeof(AVIndexEntry))
1132         return AVERROR(ENOMEM);
1133
1134     index_table->ptses      = av_mallocz(index_table->nb_ptses *
1135                                          sizeof(int64_t));
1136     index_table->fake_index = av_mallocz(index_table->nb_ptses *
1137                                          sizeof(AVIndexEntry));
1138     if (!index_table->ptses || !index_table->fake_index) {
1139         av_freep(&index_table->ptses);
1140         return AVERROR(ENOMEM);
1141     }
1142
1143     /* we may have a few bad TemporalOffsets
1144      * make sure the corresponding PTSes don't have the bogus value 0 */
1145     for (x = 0; x < index_table->nb_ptses; x++)
1146         index_table->ptses[x] = AV_NOPTS_VALUE;
1147
1148     /**
1149      * We have this:
1150      *
1151      * x  TemporalOffset
1152      * 0:  0
1153      * 1:  1
1154      * 2:  1
1155      * 3: -2
1156      * 4:  1
1157      * 5:  1
1158      * 6: -2
1159      *
1160      * We want to transform it into this:
1161      *
1162      * x  DTS PTS
1163      * 0: -1   0
1164      * 1:  0   3
1165      * 2:  1   1
1166      * 3:  2   2
1167      * 4:  3   6
1168      * 5:  4   4
1169      * 6:  5   5
1170      *
1171      * We do this by bucket sorting x by x+TemporalOffset[x] into mxf->ptses,
1172      * then settings mxf->first_dts = -max(TemporalOffset[x]).
1173      * The latter makes DTS <= PTS.
1174      */
1175     for (i = x = 0; i < index_table->nb_segments; i++) {
1176         MXFIndexTableSegment *s = index_table->segments[i];
1177         int index_delta = 1;
1178         int n = s->nb_index_entries;
1179
1180         if (s->nb_index_entries == 2 * s->index_duration + 1) {
1181             index_delta = 2;    /* Avid index */
1182             /* ignore the last entry - it's the size of the essence container */
1183             n--;
1184         }
1185
1186         for (j = 0; j < n; j += index_delta, x++) {
1187             int offset = s->temporal_offset_entries[j] / index_delta;
1188             int index  = x + offset;
1189
1190             if (x >= index_table->nb_ptses) {
1191                 av_log(mxf->fc, AV_LOG_ERROR,
1192                        "x >= nb_ptses - IndexEntryCount %i < IndexDuration %"PRId64"?\n",
1193                        s->nb_index_entries, s->index_duration);
1194                 break;
1195             }
1196
1197             index_table->fake_index[x].timestamp = x;
1198             index_table->fake_index[x].flags = !(s->flag_entries[j] & 0x30) ? AVINDEX_KEYFRAME : 0;
1199
1200             if (index < 0 || index >= index_table->nb_ptses) {
1201                 av_log(mxf->fc, AV_LOG_ERROR,
1202                        "index entry %i + TemporalOffset %i = %i, which is out of bounds\n",
1203                        x, offset, index);
1204                 continue;
1205             }
1206
1207             index_table->ptses[index] = x;
1208             max_temporal_offset = FFMAX(max_temporal_offset, offset);
1209         }
1210     }
1211
1212     index_table->first_dts = -max_temporal_offset;
1213
1214     return 0;
1215 }
1216
1217 /**
1218  * Sorts and collects index table segments into index tables.
1219  * Also computes PTSes if possible.
1220  */
1221 static int mxf_compute_index_tables(MXFContext *mxf)
1222 {
1223     int i, j, k, ret, nb_sorted_segments;
1224     MXFIndexTableSegment **sorted_segments = NULL;
1225
1226     if ((ret = mxf_get_sorted_table_segments(mxf, &nb_sorted_segments, &sorted_segments)) ||
1227         nb_sorted_segments <= 0) {
1228         av_log(mxf->fc, AV_LOG_WARNING, "broken or empty index\n");
1229         return 0;
1230     }
1231
1232     /* sanity check and count unique BodySIDs/IndexSIDs */
1233     for (i = 0; i < nb_sorted_segments; i++) {
1234         if (i == 0 || sorted_segments[i-1]->index_sid != sorted_segments[i]->index_sid)
1235             mxf->nb_index_tables++;
1236         else if (sorted_segments[i-1]->body_sid != sorted_segments[i]->body_sid) {
1237             av_log(mxf->fc, AV_LOG_ERROR, "found inconsistent BodySID\n");
1238             ret = AVERROR_INVALIDDATA;
1239             goto finish_decoding_index;
1240         }
1241     }
1242
1243     mxf->index_tables = av_mallocz_array(mxf->nb_index_tables,
1244                                          sizeof(*mxf->index_tables));
1245     if (!mxf->index_tables) {
1246         av_log(mxf->fc, AV_LOG_ERROR, "failed to allocate index tables\n");
1247         ret = AVERROR(ENOMEM);
1248         goto finish_decoding_index;
1249     }
1250
1251     /* distribute sorted segments to index tables */
1252     for (i = j = 0; i < nb_sorted_segments; i++) {
1253         if (i != 0 && sorted_segments[i-1]->index_sid != sorted_segments[i]->index_sid) {
1254             /* next IndexSID */
1255             j++;
1256         }
1257
1258         mxf->index_tables[j].nb_segments++;
1259     }
1260
1261     for (i = j = 0; j < mxf->nb_index_tables; i += mxf->index_tables[j++].nb_segments) {
1262         MXFIndexTable *t = &mxf->index_tables[j];
1263
1264         t->segments = av_mallocz_array(t->nb_segments,
1265                                        sizeof(*t->segments));
1266
1267         if (!t->segments) {
1268             av_log(mxf->fc, AV_LOG_ERROR, "failed to allocate IndexTableSegment"
1269                    " pointer array\n");
1270             ret = AVERROR(ENOMEM);
1271             goto finish_decoding_index;
1272         }
1273
1274         if (sorted_segments[i]->index_start_position)
1275             av_log(mxf->fc, AV_LOG_WARNING, "IndexSID %i starts at EditUnit %"PRId64" - seeking may not work as expected\n",
1276                    sorted_segments[i]->index_sid, sorted_segments[i]->index_start_position);
1277
1278         memcpy(t->segments, &sorted_segments[i], t->nb_segments * sizeof(MXFIndexTableSegment*));
1279         t->index_sid = sorted_segments[i]->index_sid;
1280         t->body_sid = sorted_segments[i]->body_sid;
1281
1282         if ((ret = mxf_compute_ptses_fake_index(mxf, t)) < 0)
1283             goto finish_decoding_index;
1284
1285         /* fix zero IndexDurations */
1286         for (k = 0; k < t->nb_segments; k++) {
1287             if (t->segments[k]->index_duration)
1288                 continue;
1289
1290             if (t->nb_segments > 1)
1291                 av_log(mxf->fc, AV_LOG_WARNING, "IndexSID %i segment %i has zero IndexDuration and there's more than one segment\n",
1292                        t->index_sid, k);
1293
1294             if (mxf->fc->nb_streams <= 0) {
1295                 av_log(mxf->fc, AV_LOG_WARNING, "no streams?\n");
1296                 break;
1297             }
1298
1299             /* assume the first stream's duration is reasonable
1300              * leave index_duration = 0 on further segments in case we have any (unlikely)
1301              */
1302             t->segments[k]->index_duration = mxf->fc->streams[0]->duration;
1303             break;
1304         }
1305     }
1306
1307     ret = 0;
1308 finish_decoding_index:
1309     av_free(sorted_segments);
1310     return ret;
1311 }
1312
1313 static int mxf_is_intra_only(MXFDescriptor *d)
1314 {
1315     return mxf_get_codec_ul(mxf_intra_only_essence_container_uls,
1316                             &d->essence_container_ul)->id != AV_CODEC_ID_NONE ||
1317            mxf_get_codec_ul(mxf_intra_only_picture_essence_coding_uls,
1318                             &d->essence_codec_ul)->id     != AV_CODEC_ID_NONE;
1319 }
1320
1321 static int mxf_parse_structural_metadata(MXFContext *mxf)
1322 {
1323     MXFPackage *material_package = NULL;
1324     MXFPackage *temp_package = NULL;
1325     int i, j, k, ret;
1326
1327     av_dlog(mxf->fc, "metadata sets count %d\n", mxf->metadata_sets_count);
1328     /* TODO: handle multiple material packages (OP3x) */
1329     for (i = 0; i < mxf->packages_count; i++) {
1330         material_package = mxf_resolve_strong_ref(mxf, &mxf->packages_refs[i], MaterialPackage);
1331         if (material_package) break;
1332     }
1333     if (!material_package) {
1334         av_log(mxf->fc, AV_LOG_ERROR, "no material package found\n");
1335         return AVERROR_INVALIDDATA;
1336     }
1337
1338     for (i = 0; i < material_package->tracks_count; i++) {
1339         MXFPackage *source_package = NULL;
1340         MXFTrack *material_track = NULL;
1341         MXFTrack *source_track = NULL;
1342         MXFTrack *temp_track = NULL;
1343         MXFDescriptor *descriptor = NULL;
1344         MXFStructuralComponent *component = NULL;
1345         UID *essence_container_ul = NULL;
1346         const MXFCodecUL *codec_ul = NULL;
1347         const MXFCodecUL *container_ul = NULL;
1348         const MXFCodecUL *pix_fmt_ul = NULL;
1349         AVStream *st;
1350
1351         if (!(material_track = mxf_resolve_strong_ref(mxf, &material_package->tracks_refs[i], Track))) {
1352             av_log(mxf->fc, AV_LOG_ERROR, "could not resolve material track strong ref\n");
1353             continue;
1354         }
1355
1356         if (!(material_track->sequence = mxf_resolve_strong_ref(mxf, &material_track->sequence_ref, Sequence))) {
1357             av_log(mxf->fc, AV_LOG_ERROR, "could not resolve material track sequence strong ref\n");
1358             continue;
1359         }
1360
1361         /* TODO: handle multiple source clips */
1362         for (j = 0; j < material_track->sequence->structural_components_count; j++) {
1363             /* TODO: handle timecode component */
1364             component = mxf_resolve_strong_ref(mxf, &material_track->sequence->structural_components_refs[j], SourceClip);
1365             if (!component)
1366                 continue;
1367
1368             for (k = 0; k < mxf->packages_count; k++) {
1369                 temp_package = mxf_resolve_strong_ref(mxf, &mxf->packages_refs[k], SourcePackage);
1370                 if (!temp_package)
1371                     continue;
1372                 if (!memcmp(temp_package->package_uid, component->source_package_uid, 16)) {
1373                     source_package = temp_package;
1374                     break;
1375                 }
1376             }
1377             if (!source_package) {
1378                 av_dlog(mxf->fc, "material track %d: no corresponding source package found\n", material_track->track_id);
1379                 break;
1380             }
1381             for (k = 0; k < source_package->tracks_count; k++) {
1382                 if (!(temp_track = mxf_resolve_strong_ref(mxf, &source_package->tracks_refs[k], Track))) {
1383                     av_log(mxf->fc, AV_LOG_ERROR, "could not resolve source track strong ref\n");
1384                     ret = AVERROR_INVALIDDATA;
1385                     goto fail_and_free;
1386                 }
1387                 if (temp_track->track_id == component->source_track_id) {
1388                     source_track = temp_track;
1389                     break;
1390                 }
1391             }
1392             if (!source_track) {
1393                 av_log(mxf->fc, AV_LOG_ERROR, "material track %d: no corresponding source track found\n", material_track->track_id);
1394                 break;
1395             }
1396         }
1397         if (!source_track || !component)
1398             continue;
1399
1400         if (!(source_track->sequence = mxf_resolve_strong_ref(mxf, &source_track->sequence_ref, Sequence))) {
1401             av_log(mxf->fc, AV_LOG_ERROR, "could not resolve source track sequence strong ref\n");
1402             ret = AVERROR_INVALIDDATA;
1403             goto fail_and_free;
1404         }
1405
1406         /* 0001GL00.MXF.A1.mxf_opatom.mxf has the same SourcePackageID as 0001GL.MXF.V1.mxf_opatom.mxf
1407          * This would result in both files appearing to have two streams. Work around this by sanity checking DataDefinition */
1408         if (memcmp(material_track->sequence->data_definition_ul, source_track->sequence->data_definition_ul, 16)) {
1409             av_log(mxf->fc, AV_LOG_ERROR, "material track %d: DataDefinition mismatch\n", material_track->track_id);
1410             continue;
1411         }
1412
1413         st = avformat_new_stream(mxf->fc, NULL);
1414         if (!st) {
1415             av_log(mxf->fc, AV_LOG_ERROR, "could not allocate stream\n");
1416             ret = AVERROR(ENOMEM);
1417             goto fail_and_free;
1418         }
1419         st->id = source_track->track_id;
1420         st->priv_data = source_track;
1421         st->duration = component->duration;
1422         if (st->duration == -1)
1423             st->duration = AV_NOPTS_VALUE;
1424         st->start_time = component->start_position;
1425         avpriv_set_pts_info(st, 64, material_track->edit_rate.den, material_track->edit_rate.num);
1426
1427         PRINT_KEY(mxf->fc, "data definition   ul", source_track->sequence->data_definition_ul);
1428         codec_ul = mxf_get_codec_ul(ff_mxf_data_definition_uls, &source_track->sequence->data_definition_ul);
1429         st->codec->codec_type = codec_ul->id;
1430
1431         source_package->descriptor = mxf_resolve_strong_ref(mxf, &source_package->descriptor_ref, AnyType);
1432         if (source_package->descriptor) {
1433             if (source_package->descriptor->type == MultipleDescriptor) {
1434                 for (j = 0; j < source_package->descriptor->sub_descriptors_count; j++) {
1435                     MXFDescriptor *sub_descriptor = mxf_resolve_strong_ref(mxf, &source_package->descriptor->sub_descriptors_refs[j], Descriptor);
1436
1437                     if (!sub_descriptor) {
1438                         av_log(mxf->fc, AV_LOG_ERROR, "could not resolve sub descriptor strong ref\n");
1439                         continue;
1440                     }
1441                     if (sub_descriptor->linked_track_id == source_track->track_id) {
1442                         descriptor = sub_descriptor;
1443                         break;
1444                     }
1445                 }
1446             } else if (source_package->descriptor->type == Descriptor)
1447                 descriptor = source_package->descriptor;
1448         }
1449         if (!descriptor) {
1450             av_log(mxf->fc, AV_LOG_INFO, "source track %d: stream %d, no descriptor found\n", source_track->track_id, st->index);
1451             continue;
1452         }
1453         PRINT_KEY(mxf->fc, "essence codec     ul", descriptor->essence_codec_ul);
1454         PRINT_KEY(mxf->fc, "essence container ul", descriptor->essence_container_ul);
1455         essence_container_ul = &descriptor->essence_container_ul;
1456         /* HACK: replacing the original key with mxf_encrypted_essence_container
1457          * is not allowed according to s429-6, try to find correct information anyway */
1458         if (IS_KLV_KEY(essence_container_ul, mxf_encrypted_essence_container)) {
1459             av_log(mxf->fc, AV_LOG_INFO, "broken encrypted mxf file\n");
1460             for (k = 0; k < mxf->metadata_sets_count; k++) {
1461                 MXFMetadataSet *metadata = mxf->metadata_sets[k];
1462                 if (metadata->type == CryptoContext) {
1463                     essence_container_ul = &((MXFCryptoContext *)metadata)->source_container_ul;
1464                     break;
1465                 }
1466             }
1467         }
1468
1469         /* TODO: drop PictureEssenceCoding and SoundEssenceCompression, only check EssenceContainer */
1470         codec_ul = mxf_get_codec_ul(ff_mxf_codec_uls, &descriptor->essence_codec_ul);
1471         st->codec->codec_id = codec_ul->id;
1472
1473         if (st->codec->codec_type == AVMEDIA_TYPE_VIDEO) {
1474             source_track->intra_only = mxf_is_intra_only(descriptor);
1475             container_ul = mxf_get_codec_ul(mxf_picture_essence_container_uls, essence_container_ul);
1476             if (st->codec->codec_id == AV_CODEC_ID_NONE)
1477                 st->codec->codec_id = container_ul->id;
1478             st->codec->width = descriptor->width;
1479             /* Field height, not frame height */
1480             st->codec->height = descriptor->height;
1481             switch (descriptor->frame_layout) {
1482                 case SegmentedFrame:
1483                     /* This one is a weird layout I don't fully understand. */
1484                     av_log(mxf->fc, AV_LOG_INFO,
1485                            "SegmentedFrame layout isn't currently supported\n");
1486                     break;
1487                 case FullFrame:
1488                     st->codec->field_order = AV_FIELD_PROGRESSIVE;
1489                     break;
1490                 case OneField:
1491                     /* Every other line is stored and needs to be duplicated. */
1492                     av_log(mxf->fc, AV_LOG_INFO,
1493                            "OneField frame layout isn't currently supported\n");
1494                     break;
1495                     /* The correct thing to do here is fall through, but by
1496                      * breaking we might be able to decode some streams at half
1497                      * the vertical resolution, rather than not al all.
1498                      * It's also for compatibility with the old behavior. */
1499                 case SeparateFields:
1500                 case MixedFields:
1501                     /* Turn field height into frame height. */
1502                     st->codec->height *= 2;
1503                 default:
1504                     av_log(mxf->fc, AV_LOG_INFO,
1505                            "Unknown frame layout type: %d\n",
1506                            descriptor->frame_layout);
1507             }
1508             if (st->codec->codec_id == AV_CODEC_ID_RAWVIDEO) {
1509                 st->codec->pix_fmt = descriptor->pix_fmt;
1510                 if (st->codec->pix_fmt == AV_PIX_FMT_NONE) {
1511                     pix_fmt_ul = mxf_get_codec_ul(ff_mxf_pixel_format_uls,
1512                                                   &descriptor->essence_codec_ul);
1513                     st->codec->pix_fmt = pix_fmt_ul->id;
1514                     if (st->codec->pix_fmt == AV_PIX_FMT_NONE) {
1515                         /* support files created before RP224v10 by defaulting to UYVY422
1516                            if subsampling is 4:2:2 and component depth is 8-bit */
1517                         if (descriptor->horiz_subsampling == 2 &&
1518                             descriptor->vert_subsampling == 1 &&
1519                             descriptor->component_depth == 8) {
1520                             st->codec->pix_fmt = AV_PIX_FMT_UYVY422;
1521                         }
1522                     }
1523                 }
1524             }
1525             st->need_parsing = AVSTREAM_PARSE_HEADERS;
1526         } else if (st->codec->codec_type == AVMEDIA_TYPE_AUDIO) {
1527             container_ul = mxf_get_codec_ul(mxf_sound_essence_container_uls, essence_container_ul);
1528             if (st->codec->codec_id == AV_CODEC_ID_NONE)
1529                 st->codec->codec_id = container_ul->id;
1530             st->codec->channels = descriptor->channels;
1531             st->codec->bits_per_coded_sample = descriptor->bits_per_sample;
1532
1533             if (descriptor->sample_rate.den > 0) {
1534                 st->codec->sample_rate = descriptor->sample_rate.num / descriptor->sample_rate.den;
1535                 avpriv_set_pts_info(st, 64, descriptor->sample_rate.den, descriptor->sample_rate.num);
1536             } else {
1537                 av_log(mxf->fc, AV_LOG_WARNING, "invalid sample rate (%d/%d) "
1538                        "found for stream #%d, time base forced to 1/48000\n",
1539                        descriptor->sample_rate.num, descriptor->sample_rate.den,
1540                        st->index);
1541                 avpriv_set_pts_info(st, 64, 1, 48000);
1542             }
1543
1544             /* TODO: implement AV_CODEC_ID_RAWAUDIO */
1545             if (st->codec->codec_id == AV_CODEC_ID_PCM_S16LE) {
1546                 if (descriptor->bits_per_sample > 16 && descriptor->bits_per_sample <= 24)
1547                     st->codec->codec_id = AV_CODEC_ID_PCM_S24LE;
1548                 else if (descriptor->bits_per_sample == 32)
1549                     st->codec->codec_id = AV_CODEC_ID_PCM_S32LE;
1550             } else if (st->codec->codec_id == AV_CODEC_ID_PCM_S16BE) {
1551                 if (descriptor->bits_per_sample > 16 && descriptor->bits_per_sample <= 24)
1552                     st->codec->codec_id = AV_CODEC_ID_PCM_S24BE;
1553                 else if (descriptor->bits_per_sample == 32)
1554                     st->codec->codec_id = AV_CODEC_ID_PCM_S32BE;
1555             } else if (st->codec->codec_id == AV_CODEC_ID_MP2) {
1556                 st->need_parsing = AVSTREAM_PARSE_FULL;
1557             }
1558         }
1559         if (descriptor->extradata) {
1560             st->codec->extradata = av_mallocz(descriptor->extradata_size + FF_INPUT_BUFFER_PADDING_SIZE);
1561             if (st->codec->extradata) {
1562                 memcpy(st->codec->extradata, descriptor->extradata, descriptor->extradata_size);
1563                 st->codec->extradata_size = descriptor->extradata_size;
1564             }
1565         } else if (st->codec->codec_id == AV_CODEC_ID_H264) {
1566             ret = ff_generate_avci_extradata(st);
1567             if (ret < 0)
1568                 return ret;
1569         }
1570         if (st->codec->codec_type != AVMEDIA_TYPE_DATA && (*essence_container_ul)[15] > 0x01) {
1571             /* TODO: decode timestamps */
1572             st->need_parsing = AVSTREAM_PARSE_TIMESTAMPS;
1573         }
1574     }
1575
1576     ret = 0;
1577 fail_and_free:
1578     return ret;
1579 }
1580
1581 static const MXFMetadataReadTableEntry mxf_metadata_read_table[] = {
1582     { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x05,0x01,0x00 }, mxf_read_primer_pack },
1583     { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x02,0x01,0x00 }, mxf_read_partition_pack },
1584     { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x02,0x02,0x00 }, mxf_read_partition_pack },
1585     { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x02,0x03,0x00 }, mxf_read_partition_pack },
1586     { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x02,0x04,0x00 }, mxf_read_partition_pack },
1587     { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x03,0x01,0x00 }, mxf_read_partition_pack },
1588     { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x03,0x02,0x00 }, mxf_read_partition_pack },
1589     { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x03,0x03,0x00 }, mxf_read_partition_pack },
1590     { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x03,0x04,0x00 }, mxf_read_partition_pack },
1591     { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x04,0x02,0x00 }, mxf_read_partition_pack },
1592     { { 0x06,0x0E,0x2B,0x34,0x02,0x05,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x04,0x04,0x00 }, mxf_read_partition_pack },
1593     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x18,0x00 }, mxf_read_content_storage, 0, AnyType },
1594     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x37,0x00 }, mxf_read_source_package, sizeof(MXFPackage), SourcePackage },
1595     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x36,0x00 }, mxf_read_material_package, sizeof(MXFPackage), MaterialPackage },
1596     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x0F,0x00 }, mxf_read_sequence, sizeof(MXFSequence), Sequence },
1597     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x11,0x00 }, mxf_read_source_clip, sizeof(MXFStructuralComponent), SourceClip },
1598     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x44,0x00 }, mxf_read_generic_descriptor, sizeof(MXFDescriptor), MultipleDescriptor },
1599     { { 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 */
1600     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x28,0x00 }, mxf_read_generic_descriptor, sizeof(MXFDescriptor), Descriptor }, /* CDCI */
1601     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x29,0x00 }, mxf_read_generic_descriptor, sizeof(MXFDescriptor), Descriptor }, /* RGBA */
1602     { { 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 */
1603     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x48,0x00 }, mxf_read_generic_descriptor, sizeof(MXFDescriptor), Descriptor }, /* Wave */
1604     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x47,0x00 }, mxf_read_generic_descriptor, sizeof(MXFDescriptor), Descriptor }, /* AES3 */
1605     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x3A,0x00 }, mxf_read_track, sizeof(MXFTrack), Track }, /* Static Track */
1606     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x01,0x01,0x01,0x01,0x3B,0x00 }, mxf_read_track, sizeof(MXFTrack), Track }, /* Generic Track */
1607     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x04,0x01,0x02,0x02,0x00,0x00 }, mxf_read_cryptographic_context, sizeof(MXFCryptoContext), CryptoContext },
1608     { { 0x06,0x0E,0x2B,0x34,0x02,0x53,0x01,0x01,0x0d,0x01,0x02,0x01,0x01,0x10,0x01,0x00 }, mxf_read_index_table_segment, sizeof(MXFIndexTableSegment), IndexTableSegment },
1609     { { 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00 }, NULL, 0, AnyType },
1610 };
1611
1612 static int mxf_read_local_tags(MXFContext *mxf, KLVPacket *klv, MXFMetadataReadFunc *read_child, int ctx_size, enum MXFMetadataSetType type)
1613 {
1614     AVIOContext *pb = mxf->fc->pb;
1615     MXFMetadataSet *ctx = ctx_size ? av_mallocz(ctx_size) : mxf;
1616     uint64_t klv_end = avio_tell(pb) + klv->length;
1617
1618     if (!ctx)
1619         return AVERROR(ENOMEM);
1620     while (avio_tell(pb) + 4 < klv_end && !pb->eof_reached) {
1621         int ret;
1622         int tag = avio_rb16(pb);
1623         int size = avio_rb16(pb); /* KLV specified by 0x53 */
1624         uint64_t next = avio_tell(pb) + size;
1625         UID uid = {0};
1626
1627         av_dlog(mxf->fc, "local tag %#04x size %d\n", tag, size);
1628         if (!size) { /* ignore empty tag, needed for some files with empty UMID tag */
1629             av_log(mxf->fc, AV_LOG_ERROR, "local tag %#04x with 0 size\n", tag);
1630             continue;
1631         }
1632         if (tag > 0x7FFF) { /* dynamic tag */
1633             int i;
1634             for (i = 0; i < mxf->local_tags_count; i++) {
1635                 int local_tag = AV_RB16(mxf->local_tags+i*18);
1636                 if (local_tag == tag) {
1637                     memcpy(uid, mxf->local_tags+i*18+2, 16);
1638                     av_dlog(mxf->fc, "local tag %#04x\n", local_tag);
1639                     PRINT_KEY(mxf->fc, "uid", uid);
1640                 }
1641             }
1642         }
1643         if (ctx_size && tag == 0x3C0A)
1644             avio_read(pb, ctx->uid, 16);
1645         else if ((ret = read_child(ctx, pb, tag, size, uid, -1)) < 0)
1646             return ret;
1647
1648         /* Accept the 64k local set limit being exceeded (Avid). Don't accept
1649          * it extending past the end of the KLV though (zzuf5.mxf). */
1650         if (avio_tell(pb) > klv_end) {
1651             if (ctx_size)
1652                 av_free(ctx);
1653
1654             av_log(mxf->fc, AV_LOG_ERROR,
1655                    "local tag %#04x extends past end of local set @ %#"PRIx64"\n",
1656                    tag, klv->offset);
1657             return AVERROR_INVALIDDATA;
1658         } else if (avio_tell(pb) <= next)   /* only seek forward, else this can loop for a long time */
1659             avio_seek(pb, next, SEEK_SET);
1660     }
1661     if (ctx_size) ctx->type = type;
1662     return ctx_size ? mxf_add_metadata_set(mxf, ctx) : 0;
1663 }
1664
1665 /**
1666  * Seeks to the previous partition, if possible
1667  * @return <= 0 if we should stop parsing, > 0 if we should keep going
1668  */
1669 static int mxf_seek_to_previous_partition(MXFContext *mxf)
1670 {
1671     AVIOContext *pb = mxf->fc->pb;
1672
1673     if (!mxf->current_partition ||
1674         mxf->run_in + mxf->current_partition->previous_partition <= mxf->last_forward_tell)
1675         return 0;   /* we've parsed all partitions */
1676
1677     /* seek to previous partition */
1678     avio_seek(pb, mxf->run_in + mxf->current_partition->previous_partition, SEEK_SET);
1679     mxf->current_partition = NULL;
1680
1681     av_dlog(mxf->fc, "seeking to previous partition\n");
1682
1683     return 1;
1684 }
1685
1686 /**
1687  * Called when essence is encountered
1688  * @return <= 0 if we should stop parsing, > 0 if we should keep going
1689  */
1690 static int mxf_parse_handle_essence(MXFContext *mxf)
1691 {
1692     AVIOContext *pb = mxf->fc->pb;
1693     int64_t ret;
1694
1695     if (mxf->parsing_backward) {
1696         return mxf_seek_to_previous_partition(mxf);
1697     } else {
1698         if (!mxf->footer_partition) {
1699             av_dlog(mxf->fc, "no footer\n");
1700             return 0;
1701         }
1702
1703         av_dlog(mxf->fc, "seeking to footer\n");
1704
1705         /* remember where we were so we don't end up seeking further back than this */
1706         mxf->last_forward_tell = avio_tell(pb);
1707
1708         if (!pb->seekable) {
1709             av_log(mxf->fc, AV_LOG_INFO, "file is not seekable - not parsing footer\n");
1710             return -1;
1711         }
1712
1713         /* seek to footer partition and parse backward */
1714         if ((ret = avio_seek(pb, mxf->run_in + mxf->footer_partition, SEEK_SET)) < 0) {
1715             av_log(mxf->fc, AV_LOG_ERROR, "failed to seek to footer @ 0x%"PRIx64" (%"PRId64") - partial file?\n",
1716                    mxf->run_in + mxf->footer_partition, ret);
1717             return ret;
1718         }
1719
1720         mxf->current_partition = NULL;
1721         mxf->parsing_backward = 1;
1722     }
1723
1724     return 1;
1725 }
1726
1727 /**
1728  * Called when the next partition or EOF is encountered
1729  * @return <= 0 if we should stop parsing, > 0 if we should keep going
1730  */
1731 static int mxf_parse_handle_partition_or_eof(MXFContext *mxf)
1732 {
1733     return mxf->parsing_backward ? mxf_seek_to_previous_partition(mxf) : 1;
1734 }
1735
1736 /**
1737  * Figure out the proper offset and length of the essence container
1738  * in each partition
1739  */
1740 static void mxf_compute_essence_containers(MXFContext *mxf)
1741 {
1742     int x;
1743
1744     /* everything is already correct */
1745     if (mxf->op == OPAtom)
1746         return;
1747
1748     for (x = 0; x < mxf->partitions_count; x++) {
1749         MXFPartition *p = &mxf->partitions[x];
1750
1751         if (!p->body_sid)
1752             continue;       /* BodySID == 0 -> no essence */
1753
1754         if (x >= mxf->partitions_count - 1)
1755             break;          /* last partition - can't compute length (and we don't need to) */
1756
1757         /* essence container spans to the next partition */
1758         p->essence_length = mxf->partitions[x+1].this_partition - p->essence_offset;
1759
1760         if (p->essence_length < 0) {
1761             /* next ThisPartition < essence_offset */
1762             p->essence_length = 0;
1763             av_log(mxf->fc, AV_LOG_ERROR,
1764                    "partition %i: bad ThisPartition = %"PRIX64"\n",
1765                    x+1, mxf->partitions[x+1].this_partition);
1766         }
1767     }
1768 }
1769
1770 static int64_t round_to_kag(int64_t position, int kag_size)
1771 {
1772     /* TODO: account for run-in? the spec isn't clear whether KAG should account for it */
1773     /* NOTE: kag_size may be any integer between 1 - 2^10 */
1774     int64_t ret = (position / kag_size) * kag_size;
1775     return ret == position ? ret : ret + kag_size;
1776 }
1777
1778 static inline void compute_partition_essence_offset(AVFormatContext *s,
1779                                                     MXFContext *mxf,
1780                                                     KLVPacket *klv)
1781 {
1782     MXFPartition *cur_part = mxf->current_partition;
1783     /* for OP1a we compute essence_offset
1784      * for OPAtom we point essence_offset after the KL
1785      *     (usually op1a_essence_offset + 20 or 25)
1786      * TODO: for OP1a we could eliminate this entire if statement, always
1787      *       stopping parsing at op1a_essence_offset
1788      *       for OPAtom we still need the actual essence_offset though
1789      *       (the KL's length can vary)
1790      */
1791     int64_t op1a_essence_offset =
1792         round_to_kag(cur_part->this_partition + cur_part->pack_length,
1793                      cur_part->kag_size) +
1794         round_to_kag(cur_part->header_byte_count, cur_part->kag_size) +
1795         round_to_kag(cur_part->index_byte_count, cur_part->kag_size);
1796
1797     if (mxf->op == OPAtom) {
1798         /* point essence_offset to the actual data
1799          * OPAtom has all the essence in one big KLV
1800          */
1801         cur_part->essence_offset = avio_tell(s->pb);
1802         cur_part->essence_length = klv->length;
1803     } else {
1804         /* NOTE: op1a_essence_offset may be less than to klv.offset
1805          * (C0023S01.mxf)  */
1806         cur_part->essence_offset = op1a_essence_offset;
1807     }
1808 }
1809
1810 static int is_pcm(enum AVCodecID codec_id)
1811 {
1812     /* we only care about "normal" PCM codecs until we get samples */
1813     return codec_id >= AV_CODEC_ID_PCM_S16LE && codec_id < AV_CODEC_ID_PCM_S24DAUD;
1814 }
1815
1816 /**
1817  * Deal with the case where for some audio atoms EditUnitByteCount is
1818  * very small (2, 4..). In those cases we should read more than one
1819  * sample per call to mxf_read_packet().
1820  */
1821 static void mxf_handle_small_eubc(AVFormatContext *s)
1822 {
1823     MXFContext *mxf = s->priv_data;
1824
1825     /* assuming non-OPAtom == frame wrapped
1826      * no sane writer would wrap 2 byte PCM packets with 20 byte headers.. */
1827     if (mxf->op != OPAtom)
1828         return;
1829
1830     /* expect PCM with exactly one index table segment and a small (< 32) EUBC */
1831     if (s->nb_streams != 1                                     ||
1832         s->streams[0]->codec->codec_type != AVMEDIA_TYPE_AUDIO ||
1833         !is_pcm(s->streams[0]->codec->codec_id)                ||
1834         mxf->nb_index_tables != 1                              ||
1835         mxf->index_tables[0].nb_segments != 1                  ||
1836         mxf->index_tables[0].segments[0]->edit_unit_byte_count >= 32)
1837         return;
1838
1839     /* arbitrarily default to 48 kHz PAL audio frame size */
1840     /* TODO: We could compute this from the ratio between the audio
1841      *       and video edit rates for 48 kHz NTSC we could use the
1842      *       1802-1802-1802-1802-1801 pattern. */
1843     mxf->edit_units_per_packet = 1920;
1844 }
1845
1846 static int mxf_read_header(AVFormatContext *s)
1847 {
1848     MXFContext *mxf = s->priv_data;
1849     KLVPacket klv;
1850     int64_t essence_offset = 0;
1851     int ret;
1852
1853     mxf->last_forward_tell = INT64_MAX;
1854     mxf->edit_units_per_packet = 1;
1855
1856     if (!mxf_read_sync(s->pb, mxf_header_partition_pack_key, 14)) {
1857         av_log(s, AV_LOG_ERROR, "could not find header partition pack key\n");
1858         return AVERROR_INVALIDDATA;
1859     }
1860     avio_seek(s->pb, -14, SEEK_CUR);
1861     mxf->fc = s;
1862     mxf->run_in = avio_tell(s->pb);
1863
1864     while (!s->pb->eof_reached) {
1865         const MXFMetadataReadTableEntry *metadata;
1866
1867         if (klv_read_packet(&klv, s->pb) < 0) {
1868             /* EOF - seek to previous partition or stop */
1869             if(mxf_parse_handle_partition_or_eof(mxf) <= 0)
1870                 break;
1871             else
1872                 continue;
1873         }
1874
1875         PRINT_KEY(s, "read header", klv.key);
1876         av_dlog(s, "size %"PRIu64" offset %#"PRIx64"\n", klv.length, klv.offset);
1877         if (IS_KLV_KEY(klv.key, mxf_encrypted_triplet_key) ||
1878             IS_KLV_KEY(klv.key, mxf_essence_element_key) ||
1879             IS_KLV_KEY(klv.key, mxf_avid_essence_element_key) ||
1880             IS_KLV_KEY(klv.key, mxf_system_item_key)) {
1881
1882             if (!mxf->current_partition) {
1883                 av_log(mxf->fc, AV_LOG_ERROR,
1884                        "found essence prior to first PartitionPack\n");
1885                 return AVERROR_INVALIDDATA;
1886             }
1887
1888             if (!mxf->current_partition->essence_offset) {
1889                 compute_partition_essence_offset(s, mxf, &klv);
1890             }
1891
1892             if (!essence_offset)
1893                 essence_offset = klv.offset;
1894
1895             /* seek to footer, previous partition or stop */
1896             if (mxf_parse_handle_essence(mxf) <= 0)
1897                 break;
1898             continue;
1899         } else if (!memcmp(klv.key, mxf_header_partition_pack_key, 13) &&
1900                    klv.key[13] >= 2 && klv.key[13] <= 4 && mxf->current_partition) {
1901             /* next partition pack - keep going, seek to previous partition or stop */
1902             if(mxf_parse_handle_partition_or_eof(mxf) <= 0)
1903                 break;
1904             else if (mxf->parsing_backward)
1905                 continue;
1906             /* we're still parsing forward. proceed to parsing this partition pack */
1907         }
1908
1909         for (metadata = mxf_metadata_read_table; metadata->read; metadata++) {
1910             if (IS_KLV_KEY(klv.key, metadata->key)) {
1911                 int res;
1912                 if (klv.key[5] == 0x53) {
1913                     res = mxf_read_local_tags(mxf, &klv, metadata->read, metadata->ctx_size, metadata->type);
1914                 } else {
1915                     uint64_t next = avio_tell(s->pb) + klv.length;
1916                     res = metadata->read(mxf, s->pb, 0, klv.length, klv.key, klv.offset);
1917
1918                     /* only seek forward, else this can loop for a long time */
1919                     if (avio_tell(s->pb) > next) {
1920                         av_log(s, AV_LOG_ERROR, "read past end of KLV @ %#"PRIx64"\n",
1921                                klv.offset);
1922                         return AVERROR_INVALIDDATA;
1923                     }
1924
1925                     avio_seek(s->pb, next, SEEK_SET);
1926                 }
1927                 if (res < 0) {
1928                     av_log(s, AV_LOG_ERROR, "error reading header metadata\n");
1929                     return res;
1930                 }
1931                 break;
1932             }
1933         }
1934         if (!metadata->read)
1935             avio_skip(s->pb, klv.length);
1936     }
1937     /* FIXME avoid seek */
1938     if (!essence_offset)  {
1939         av_log(s, AV_LOG_ERROR, "no essence\n");
1940         return AVERROR_INVALIDDATA;
1941     }
1942     avio_seek(s->pb, essence_offset, SEEK_SET);
1943
1944     mxf_compute_essence_containers(mxf);
1945
1946     /* we need to do this before computing the index tables
1947      * to be able to fill in zero IndexDurations with st->duration */
1948     if ((ret = mxf_parse_structural_metadata(mxf)) < 0)
1949         return ret;
1950
1951     if ((ret = mxf_compute_index_tables(mxf)) < 0)
1952         return ret;
1953
1954     if (mxf->nb_index_tables > 1) {
1955         /* TODO: look up which IndexSID to use via EssenceContainerData */
1956         av_log(mxf->fc, AV_LOG_INFO, "got %i index tables - only the first one (IndexSID %i) will be used\n",
1957                mxf->nb_index_tables, mxf->index_tables[0].index_sid);
1958     } else if (mxf->nb_index_tables == 0 && mxf->op == OPAtom) {
1959         av_log(mxf->fc, AV_LOG_ERROR, "cannot demux OPAtom without an index\n");
1960         return AVERROR_INVALIDDATA;
1961     }
1962
1963     mxf_handle_small_eubc(s);
1964
1965     return 0;
1966 }
1967
1968 /**
1969  * Sets mxf->current_edit_unit based on what offset we're currently at.
1970  * @return next_ofs if OK, <0 on error
1971  */
1972 static int64_t mxf_set_current_edit_unit(MXFContext *mxf, int64_t current_offset)
1973 {
1974     int64_t last_ofs = -1, next_ofs = -1;
1975     MXFIndexTable *t = &mxf->index_tables[0];
1976
1977     /* this is called from the OP1a demuxing logic, which means there
1978      * may be no index tables */
1979     if (mxf->nb_index_tables <= 0)
1980         return -1;
1981
1982     /* find mxf->current_edit_unit so that the next edit unit starts ahead
1983      * of current_offset */
1984     while (mxf->current_edit_unit >= 0) {
1985         if (mxf_edit_unit_absolute_offset(mxf, t, mxf->current_edit_unit + 1,
1986                                           NULL, &next_ofs, 0) < 0)
1987             return -1;
1988
1989         if (next_ofs <= last_ofs) {
1990             /* large next_ofs didn't change or current_edit_unit wrapped
1991              * around this fixes the infinite loop on zzuf3.mxf */
1992             av_log(mxf->fc, AV_LOG_ERROR,
1993                    "next_ofs didn't change. not deriving packet timestamps\n");
1994             return -1;
1995         }
1996
1997         if (next_ofs > current_offset)
1998             break;
1999
2000         last_ofs = next_ofs;
2001         mxf->current_edit_unit++;
2002     }
2003
2004     /* not checking mxf->current_edit_unit >= t->nb_ptses here since CBR files
2005      * may lack IndexEntryArrays */
2006     if (mxf->current_edit_unit < 0)
2007         return -1;
2008
2009     return next_ofs;
2010 }
2011
2012 static int mxf_read_packet_old(AVFormatContext *s, AVPacket *pkt)
2013 {
2014     KLVPacket klv;
2015     MXFContext *mxf = s->priv_data;
2016
2017     while (!s->pb->eof_reached) {
2018         if (klv_read_packet(&klv, s->pb) < 0)
2019             return -1;
2020         PRINT_KEY(s, "read packet", klv.key);
2021         av_dlog(s, "size %"PRIu64" offset %#"PRIx64"\n", klv.length, klv.offset);
2022         if (IS_KLV_KEY(klv.key, mxf_encrypted_triplet_key)) {
2023             int res = mxf_decrypt_triplet(s, pkt, &klv);
2024             if (res < 0) {
2025                 av_log(s, AV_LOG_ERROR, "invalid encoded triplet\n");
2026                 return -1;
2027             }
2028             return 0;
2029         }
2030         if (IS_KLV_KEY(klv.key, mxf_essence_element_key) ||
2031             IS_KLV_KEY(klv.key, mxf_avid_essence_element_key)) {
2032             int index = mxf_get_stream_index(s, &klv);
2033             int64_t next_ofs, next_klv;
2034             AVStream *st;
2035             MXFTrack *track;
2036
2037             if (index < 0) {
2038                 av_log(s, AV_LOG_ERROR, "error getting stream index %d\n", AV_RB32(klv.key+12));
2039                 goto skip;
2040             }
2041
2042             st = s->streams[index];
2043             track = st->priv_data;
2044
2045             if (s->streams[index]->discard == AVDISCARD_ALL)
2046                 goto skip;
2047
2048             next_klv = avio_tell(s->pb) + klv.length;
2049             next_ofs = mxf_set_current_edit_unit(mxf, klv.offset);
2050
2051             if (next_ofs >= 0 && next_klv > next_ofs) {
2052                 /* if this check is hit then it's possible OPAtom was treated
2053                  * as OP1a truncate the packet since it's probably very large
2054                  * (>2 GiB is common) */
2055                 avpriv_request_sample(s,
2056                                       "OPAtom misinterpreted as OP1a?"
2057                                       "KLV for edit unit %i extending into "
2058                                       "next edit unit",
2059                                       mxf->current_edit_unit);
2060                 klv.length = next_ofs - avio_tell(s->pb);
2061             }
2062
2063             /* check for 8 channels AES3 element */
2064             if (klv.key[12] == 0x06 && klv.key[13] == 0x01 && klv.key[14] == 0x10) {
2065                 if (mxf_get_d10_aes3_packet(s->pb, s->streams[index], pkt, klv.length) < 0) {
2066                     av_log(s, AV_LOG_ERROR, "error reading D-10 aes3 frame\n");
2067                     return -1;
2068                 }
2069             } else {
2070                 int ret = av_get_packet(s->pb, pkt, klv.length);
2071                 if (ret < 0)
2072                     return ret;
2073             }
2074             pkt->stream_index = index;
2075             pkt->pos = klv.offset;
2076
2077             if (s->streams[index]->codec->codec_type == AVMEDIA_TYPE_VIDEO && next_ofs >= 0) {
2078                 /* mxf->current_edit_unit good - see if we have an
2079                  * index table to derive timestamps from */
2080                 MXFIndexTable *t = &mxf->index_tables[0];
2081
2082                 if (mxf->nb_index_tables >= 1 &&
2083                     mxf->current_edit_unit < t->nb_ptses) {
2084                     pkt->dts = mxf->current_edit_unit + t->first_dts;
2085                     pkt->pts = t->ptses[mxf->current_edit_unit];
2086                 } else if (track->intra_only) {
2087                     /* intra-only -> PTS = EditUnit.
2088                      * let utils.c figure out DTS since it can be
2089                      * < PTS if low_delay = 0 (Sony IMX30) */
2090                     pkt->pts = mxf->current_edit_unit;
2091                 }
2092             }
2093
2094             /* seek for truncated packets */
2095             avio_seek(s->pb, next_klv, SEEK_SET);
2096
2097             return 0;
2098         } else
2099         skip:
2100             avio_skip(s->pb, klv.length);
2101     }
2102     return AVERROR_EOF;
2103 }
2104
2105 static int mxf_read_packet(AVFormatContext *s, AVPacket *pkt)
2106 {
2107     MXFContext *mxf = s->priv_data;
2108     int ret, size;
2109     int64_t ret64, pos, next_pos;
2110     AVStream *st;
2111     MXFIndexTable *t;
2112     int edit_units;
2113
2114     if (mxf->op != OPAtom)
2115         return mxf_read_packet_old(s, pkt);
2116
2117     /* OPAtom - clip wrapped demuxing */
2118     /* NOTE: mxf_read_header() makes sure nb_index_tables > 0 for OPAtom */
2119     st = s->streams[0];
2120     t = &mxf->index_tables[0];
2121
2122     if (mxf->current_edit_unit >= st->duration)
2123         return AVERROR_EOF;
2124
2125     edit_units = FFMIN(mxf->edit_units_per_packet, st->duration - mxf->current_edit_unit);
2126
2127     if ((ret = mxf_edit_unit_absolute_offset(mxf, t, mxf->current_edit_unit, NULL, &pos, 1)) < 0)
2128         return ret;
2129
2130     /* compute size by finding the next edit unit or the end of the essence container
2131      * not pretty, but it works */
2132     if ((ret = mxf_edit_unit_absolute_offset(mxf, t, mxf->current_edit_unit + edit_units, NULL, &next_pos, 0)) < 0 &&
2133         (next_pos = mxf_essence_container_end(mxf, t->body_sid)) <= 0) {
2134         av_log(s, AV_LOG_ERROR, "unable to compute the size of the last packet\n");
2135         return AVERROR_INVALIDDATA;
2136     }
2137
2138     if ((size = next_pos - pos) <= 0) {
2139         av_log(s, AV_LOG_ERROR, "bad size: %i\n", size);
2140         return AVERROR_INVALIDDATA;
2141     }
2142
2143     if ((ret64 = avio_seek(s->pb, pos, SEEK_SET)) < 0)
2144         return ret64;
2145
2146         if ((ret = av_get_packet(s->pb, pkt, size)) != size)
2147             return ret < 0 ? ret : AVERROR_EOF;
2148
2149     if (st->codec->codec_type == AVMEDIA_TYPE_VIDEO && t->ptses &&
2150         mxf->current_edit_unit >= 0 && mxf->current_edit_unit < t->nb_ptses) {
2151         pkt->dts = mxf->current_edit_unit + t->first_dts;
2152         pkt->pts = t->ptses[mxf->current_edit_unit];
2153     }
2154
2155     pkt->stream_index = 0;
2156     mxf->current_edit_unit += edit_units;
2157
2158     return 0;
2159 }
2160
2161
2162 static int mxf_read_close(AVFormatContext *s)
2163 {
2164     MXFContext *mxf = s->priv_data;
2165     MXFIndexTableSegment *seg;
2166     int i;
2167
2168     av_freep(&mxf->packages_refs);
2169
2170     for (i = 0; i < s->nb_streams; i++)
2171         s->streams[i]->priv_data = NULL;
2172
2173     for (i = 0; i < mxf->metadata_sets_count; i++) {
2174         switch (mxf->metadata_sets[i]->type) {
2175         case MultipleDescriptor:
2176             av_freep(&((MXFDescriptor *)mxf->metadata_sets[i])->sub_descriptors_refs);
2177             break;
2178         case Sequence:
2179             av_freep(&((MXFSequence *)mxf->metadata_sets[i])->structural_components_refs);
2180             break;
2181         case SourcePackage:
2182         case MaterialPackage:
2183             av_freep(&((MXFPackage *)mxf->metadata_sets[i])->tracks_refs);
2184             break;
2185         case IndexTableSegment:
2186             seg = (MXFIndexTableSegment *)mxf->metadata_sets[i];
2187             av_freep(&seg->temporal_offset_entries);
2188             av_freep(&seg->flag_entries);
2189             av_freep(&seg->stream_offset_entries);
2190             break;
2191         default:
2192             break;
2193         }
2194         av_freep(&mxf->metadata_sets[i]);
2195     }
2196     av_freep(&mxf->partitions);
2197     av_freep(&mxf->metadata_sets);
2198     av_freep(&mxf->aesc);
2199     av_freep(&mxf->local_tags);
2200
2201     for (i = 0; i < mxf->nb_index_tables; i++) {
2202         av_freep(&mxf->index_tables[i].segments);
2203         av_freep(&mxf->index_tables[i].ptses);
2204         av_freep(&mxf->index_tables[i].fake_index);
2205     }
2206     av_freep(&mxf->index_tables);
2207
2208     return 0;
2209 }
2210
2211 static int mxf_probe(AVProbeData *p) {
2212     uint8_t *bufp = p->buf;
2213     uint8_t *end = p->buf + p->buf_size;
2214
2215     if (p->buf_size < sizeof(mxf_header_partition_pack_key))
2216         return 0;
2217
2218     /* Must skip Run-In Sequence and search for MXF header partition pack key SMPTE 377M 5.5 */
2219     end -= sizeof(mxf_header_partition_pack_key);
2220     for (; bufp < end; bufp++) {
2221         if (IS_KLV_KEY(bufp, mxf_header_partition_pack_key))
2222             return AVPROBE_SCORE_MAX;
2223     }
2224     return 0;
2225 }
2226
2227 /* rudimentary byte seek */
2228 /* XXX: use MXF Index */
2229 static int mxf_read_seek(AVFormatContext *s, int stream_index, int64_t sample_time, int flags)
2230 {
2231     AVStream *st = s->streams[stream_index];
2232     int64_t seconds;
2233     MXFContext* mxf = s->priv_data;
2234     int64_t seekpos;
2235     int ret;
2236     MXFIndexTable *t;
2237
2238     if (mxf->nb_index_tables <= 0) {
2239     if (!s->bit_rate)
2240         return AVERROR_INVALIDDATA;
2241     if (sample_time < 0)
2242         sample_time = 0;
2243     seconds = av_rescale(sample_time, st->time_base.num, st->time_base.den);
2244
2245     seekpos = avio_seek(s->pb, (s->bit_rate * seconds) >> 3, SEEK_SET);
2246     if (seekpos < 0)
2247         return seekpos;
2248
2249     ff_update_cur_dts(s, st, sample_time);
2250     mxf->current_edit_unit = sample_time;
2251     } else {
2252         t = &mxf->index_tables[0];
2253
2254         /* clamp above zero, else ff_index_search_timestamp() returns negative
2255          * this also means we allow seeking before the start */
2256         sample_time = FFMAX(sample_time, 0);
2257
2258         if (t->fake_index) {
2259             /* behave as if we have a proper index */
2260             if ((sample_time = ff_index_search_timestamp(t->fake_index, t->nb_ptses, sample_time, flags)) < 0)
2261                 return sample_time;
2262         } else {
2263             /* no IndexEntryArray (one or more CBR segments)
2264              * make sure we don't seek past the end */
2265             sample_time = FFMIN(sample_time, st->duration - 1);
2266         }
2267
2268         if ((ret = mxf_edit_unit_absolute_offset(mxf, t, sample_time, &sample_time, &seekpos, 1)) << 0)
2269             return ret;
2270
2271         ff_update_cur_dts(s, st, sample_time);
2272         mxf->current_edit_unit = sample_time;
2273         avio_seek(s->pb, seekpos, SEEK_SET);
2274     }
2275     return 0;
2276 }
2277
2278 AVInputFormat ff_mxf_demuxer = {
2279     .name           = "mxf",
2280     .long_name      = NULL_IF_CONFIG_SMALL("MXF (Material eXchange Format)"),
2281     .priv_data_size = sizeof(MXFContext),
2282     .read_probe     = mxf_probe,
2283     .read_header    = mxf_read_header,
2284     .read_packet    = mxf_read_packet,
2285     .read_close     = mxf_read_close,
2286     .read_seek      = mxf_read_seek,
2287 };