libavcodec/hq_hqa.c

   1 /*
   2  * Canopus HQ/HQA decoder
   3  *
   4  * This file is part of Libav.
   5  *
   6  * Libav is free software; you can redistribute it and/or
   7  * modify it under the terms of the GNU Lesser General Public
   8  * License as published by the Free Software Foundation; either
   9  * version 2.1 of the License, or (at your option) any later version.
  10  *
  11  * Libav is distributed in the hope that it will be useful,
  12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  14  * Lesser General Public License for more details.
  15  *
  16  * You should have received a copy of the GNU Lesser General Public
  17  * License along with Libav; if not, write to the Free Software
  18  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
  19  */
  20
  21 #include <stdint.h>
  22
  23 #include "libavutil/attributes.h"
  24 #include "libavutil/intreadwrite.h"
  25
  26 #include "avcodec.h"
  27 #include "bitstream.h"
  28 #include "canopus.h"
  29 #include "internal.h"
  30
  31 #include "hq_hqa.h"
  32 #include "hq_hqadsp.h"
  33
  34 /* HQ/HQA slices are a set of macroblocks belonging to a frame, and
  35  * they usually form a pseudorandom pattern (probably because it is
  36  * nicer to display on partial decode).
  37  *
  38  * For HQA it just happens that each slice is on every 8th macroblock,
  39  * but they can be on any frame width like
  40  *   X.......X.
  41  *   ......X...
  42  *   ....X.....
  43  *   ..X.......
  44  * etc.
  45  *
  46  * The original decoder has special handling for edge macroblocks,
  47  * while lavc simply aligns coded_width and coded_height.
  48  */
  49
  50 static inline void put_blocks(HQContext *c, AVFrame *pic,
  51                               int plane, int x, int y, int ilace,
  52                               int16_t *block0, int16_t *block1)
  53 {
  54     uint8_t *p = pic->data[plane] + x;
  55
  56     c->hqhqadsp.idct_put(p + y * pic->linesize[plane],
  57                          pic->linesize[plane] << ilace, block0);
  58     c->hqhqadsp.idct_put(p + (y + (ilace ? 1 : 8)) * pic->linesize[plane],
  59                          pic->linesize[plane] << ilace, block1);
  60 }
  61
  62 static int hq_decode_block(HQContext *c, BitstreamContext *bc, int16_t block[64],
  63                            int qsel, int is_chroma, int is_hqa)
  64 {
  65     const int32_t *q;
  66     int val, pos = 1;
  67
  68     memset(block, 0, 64 * sizeof(*block));
  69
  70     if (!is_hqa) {
  71         block[0] = bitstream_read_signed(bc, 9) << 6;
  72         q = ff_hq_quants[qsel][is_chroma][bitstream_read(bc, 2)];
  73     } else {
  74         q = ff_hq_quants[qsel][is_chroma][bitstream_read(bc, 2)];
  75         block[0] = bitstream_read_signed(bc, 9) << 6;
  76     }
  77
  78     for (;;) {
  79         val = bitstream_read_vlc(bc, c->hq_ac_vlc.table, 9, 2);
  80         if (val < 0)
  81             return AVERROR_INVALIDDATA;
  82
  83         pos += ff_hq_ac_skips[val];
  84         if (pos >= 64)
  85             break;
  86         block[ff_zigzag_direct[pos]] = (ff_hq_ac_syms[val] * q[pos]) >> 12;
  87         pos++;
  88     }
  89
  90     return 0;
  91 }
  92
  93 static int hq_decode_mb(HQContext *c, AVFrame *pic,
  94                         BitstreamContext *bc, int x, int y)
  95 {
  96     int qgroup, flag;
  97     int i, ret;
  98
  99     qgroup = bitstream_read(bc, 4);
 100     flag   = bitstream_read_bit(bc);
 101
 102     for (i = 0; i < 8; i++) {
 103         ret = hq_decode_block(c, bc, c->block[i], qgroup, i >= 4, 0);
 104         if (ret < 0)
 105             return ret;
 106     }
 107
 108     put_blocks(c, pic, 0, x,      y, flag, c->block[0], c->block[2]);
 109     put_blocks(c, pic, 0, x + 8,  y, flag, c->block[1], c->block[3]);
 110     put_blocks(c, pic, 2, x >> 1, y, flag, c->block[4], c->block[5]);
 111     put_blocks(c, pic, 1, x >> 1, y, flag, c->block[6], c->block[7]);
 112
 113     return 0;
 114 }
 115
 116 static int hq_decode_frame(HQContext *ctx, AVFrame *pic,
 117                            int prof_num, size_t data_size)
 118 {
 119     const HQProfile *profile;
 120     BitstreamContext bc;
 121     const uint8_t *perm, *src = ctx->gbc.buffer;
 122     uint32_t slice_off[21];
 123     int slice, start_off, next_off, i, ret;
 124
 125     if (prof_num >= NUM_HQ_PROFILES) {
 126         profile = &ff_hq_profile[0];
 127         avpriv_request_sample(ctx->avctx, "HQ Profile %d", prof_num);
 128     } else {
 129         profile = &ff_hq_profile[prof_num];
 130         av_log(ctx->avctx, AV_LOG_VERBOSE, "HQ Profile %d\n", prof_num);
 131     }
 132
 133     ctx->avctx->coded_width         = FFALIGN(profile->width,  16);
 134     ctx->avctx->coded_height        = FFALIGN(profile->height, 16);
 135     ctx->avctx->width               = profile->width;
 136     ctx->avctx->height              = profile->height;
 137     ctx->avctx->bits_per_raw_sample = 8;
 138     ctx->avctx->pix_fmt             = AV_PIX_FMT_YUV422P;
 139
 140     ret = ff_get_buffer(ctx->avctx, pic, 0);
 141     if (ret < 0) {
 142         av_log(ctx->avctx, AV_LOG_ERROR, "Could not allocate buffer.\n");
 143         return ret;
 144     }
 145
 146     /* Offsets are stored from CUV position, so adjust them accordingly. */
 147     for (i = 0; i < profile->num_slices + 1; i++)
 148         slice_off[i] = bytestream2_get_be24(&ctx->gbc) - 4;
 149
 150     next_off = 0;
 151     for (slice = 0; slice < profile->num_slices; slice++) {
 152         start_off = next_off;
 153         next_off  = profile->tab_h * (slice + 1) / profile->num_slices;
 154         perm = profile->perm_tab + start_off * profile->tab_w * 2;
 155
 156         if (slice_off[slice] < (profile->num_slices + 1) * 3 ||
 157             slice_off[slice] >= slice_off[slice + 1] ||
 158             slice_off[slice + 1] > data_size) {
 159             av_log(ctx->avctx, AV_LOG_ERROR,
 160                    "Invalid slice size %zu.\n", data_size);
 161             break;
 162         }
 163         bitstream_init8(&bc, src + slice_off[slice],
 164                         slice_off[slice + 1] - slice_off[slice]);
 165
 166         for (i = 0; i < (next_off - start_off) * profile->tab_w; i++) {
 167             ret = hq_decode_mb(ctx, pic, &bc, perm[0] * 16, perm[1] * 16);
 168             if (ret < 0) {
 169                 av_log(ctx->avctx, AV_LOG_ERROR,
 170                        "Error decoding macroblock %d at slice %d.\n", i, slice);
 171                 return ret;
 172             }
 173             perm += 2;
 174         }
 175     }
 176
 177     return 0;
 178 }
 179
 180 static int hqa_decode_mb(HQContext *c, AVFrame *pic, int qgroup,
 181                          BitstreamContext *bc, int x, int y)
 182 {
 183     int flag = 0;
 184     int i, ret, cbp;
 185
 186     cbp = bitstream_read_vlc(bc, c->hqa_cbp_vlc.table, 5, 1);
 187
 188     for (i = 0; i < 12; i++)
 189         memset(c->block[i], 0, sizeof(*c->block));
 190     for (i = 0; i < 12; i++)
 191         c->block[i][0] = -128 * (1 << 6);
 192
 193     if (cbp) {
 194         flag = bitstream_read_bit(bc);
 195
 196         cbp |= cbp << 4;
 197         if (cbp & 0x3)
 198             cbp |= 0x500;
 199         if (cbp & 0xC)
 200             cbp |= 0xA00;
 201         for (i = 0; i < 12; i++) {
 202             if (!(cbp & (1 << i)))
 203                 continue;
 204             ret = hq_decode_block(c, bc, c->block[i], qgroup, i >= 8, 1);
 205             if (ret < 0)
 206                 return ret;
 207         }
 208     }
 209
 210     put_blocks(c, pic, 3, x,      y, flag, c->block[ 0], c->block[ 2]);
 211     put_blocks(c, pic, 3, x + 8,  y, flag, c->block[ 1], c->block[ 3]);
 212     put_blocks(c, pic, 0, x,      y, flag, c->block[ 4], c->block[ 6]);
 213     put_blocks(c, pic, 0, x + 8,  y, flag, c->block[ 5], c->block[ 7]);
 214     put_blocks(c, pic, 2, x >> 1, y, flag, c->block[ 8], c->block[ 9]);
 215     put_blocks(c, pic, 1, x >> 1, y, flag, c->block[10], c->block[11]);
 216
 217     return 0;
 218 }
 219
 220 static int hqa_decode_slice(HQContext *ctx, AVFrame *pic, BitstreamContext *bc,
 221                             int quant, int slice_no, int w, int h)
 222 {
 223     int i, j, off;
 224     int ret;
 225
 226     for (i = 0; i < h; i += 16) {
 227         off = (slice_no * 16 + i * 3) & 0x70;
 228         for (j = off; j < w; j += 128) {
 229             ret = hqa_decode_mb(ctx, pic, quant, bc, j, i);
 230             if (ret < 0) {
 231                 av_log(ctx->avctx, AV_LOG_ERROR,
 232                        "Error decoding macroblock at %dx%d.\n", i, j);
 233                 return ret;
 234             }
 235         }
 236     }
 237
 238     return 0;
 239 }
 240
 241 static int hqa_decode_frame(HQContext *ctx, AVFrame *pic, size_t data_size)
 242 {
 243     BitstreamContext bc;
 244     const int num_slices = 8;
 245     uint32_t slice_off[9];
 246     int i, slice, ret;
 247     int width, height, quant;
 248     const uint8_t *src = ctx->gbc.buffer;
 249
 250     width  = bytestream2_get_be16(&ctx->gbc);
 251     height = bytestream2_get_be16(&ctx->gbc);
 252
 253     ctx->avctx->coded_width         = FFALIGN(width,  16);
 254     ctx->avctx->coded_height        = FFALIGN(height, 16);
 255     ctx->avctx->width               = width;
 256     ctx->avctx->height              = height;
 257     ctx->avctx->bits_per_raw_sample = 8;
 258     ctx->avctx->pix_fmt             = AV_PIX_FMT_YUVA422P;
 259
 260     av_log(ctx->avctx, AV_LOG_VERBOSE, "HQA Profile\n");
 261
 262     quant = bytestream2_get_byte(&ctx->gbc);
 263     bytestream2_skip(&ctx->gbc, 3);
 264     if (quant >= NUM_HQ_QUANTS) {
 265         av_log(ctx->avctx, AV_LOG_ERROR,
 266                "Invalid quantization matrix %d.\n", quant);
 267         return AVERROR_INVALIDDATA;
 268     }
 269
 270     ret = ff_get_buffer(ctx->avctx, pic, 0);
 271     if (ret < 0) {
 272         av_log(ctx->avctx, AV_LOG_ERROR, "Could not allocate buffer.\n");
 273         return ret;
 274     }
 275
 276     /* Offsets are stored from HQA1 position, so adjust them accordingly. */
 277     for (i = 0; i < num_slices + 1; i++)
 278         slice_off[i] = bytestream2_get_be32(&ctx->gbc) - 4;
 279
 280     for (slice = 0; slice < num_slices; slice++) {
 281         if (slice_off[slice] < (num_slices + 1) * 3 ||
 282             slice_off[slice] >= slice_off[slice + 1] ||
 283             slice_off[slice + 1] > data_size) {
 284             av_log(ctx->avctx, AV_LOG_ERROR,
 285                    "Invalid slice size %zu.\n", data_size);
 286             break;
 287         }
 288         bitstream_init8(&bc, src + slice_off[slice],
 289                         slice_off[slice + 1] - slice_off[slice]);
 290
 291         ret = hqa_decode_slice(ctx, pic, &bc, quant, slice, width, height);
 292         if (ret < 0)
 293             return ret;
 294     }
 295
 296     return 0;
 297 }
 298
 299 static int hq_hqa_decode_frame(AVCodecContext *avctx, void *data,
 300                                int *got_frame, AVPacket *avpkt)
 301 {
 302     HQContext *ctx = avctx->priv_data;
 303     AVFrame *pic = data;
 304     uint32_t info_tag;
 305     unsigned int data_size;
 306     int tag, ret;
 307
 308     bytestream2_init(&ctx->gbc, avpkt->data, avpkt->size);
 309     if (bytestream2_get_bytes_left(&ctx->gbc) < 4 + 4) {
 310         av_log(avctx, AV_LOG_ERROR, "Frame is too small (%d).\n", avpkt->size);
 311         return AVERROR_INVALIDDATA;
 312     }
 313
 314     info_tag = bytestream2_peek_le32(&ctx->gbc);
 315     if (info_tag == MKTAG('I', 'N', 'F', 'O')) {
 316         int info_size;
 317         bytestream2_skip(&ctx->gbc, 4);
 318         info_size = bytestream2_get_le32(&ctx->gbc);
 319         if (bytestream2_get_bytes_left(&ctx->gbc) < info_size) {
 320             av_log(avctx, AV_LOG_ERROR, "Invalid INFO size (%d).\n", info_size);
 321             return AVERROR_INVALIDDATA;
 322         }
 323         ff_canopus_parse_info_tag(avctx, ctx->gbc.buffer, info_size);
 324
 325         bytestream2_skip(&ctx->gbc, info_size);
 326     }
 327
 328     data_size = bytestream2_get_bytes_left(&ctx->gbc);
 329     if (data_size < 4) {
 330         av_log(avctx, AV_LOG_ERROR, "Frame is too small (%d).\n", data_size);
 331         return AVERROR_INVALIDDATA;
 332     }
 333
 334     /* HQ defines dimensions and number of slices, and thus slice traversal
 335      * order. HQA has no size constraint and a fixed number of slices, so it
 336      * needs a separate scheme for it. */
 337     tag = bytestream2_get_le32(&ctx->gbc);
 338     if ((tag & 0x00FFFFFF) == (MKTAG('U', 'V', 'C', ' ') & 0x00FFFFFF)) {
 339         ret = hq_decode_frame(ctx, pic, tag >> 24, data_size);
 340     } else if (tag == MKTAG('H', 'Q', 'A', '1')) {
 341         ret = hqa_decode_frame(ctx, pic, data_size);
 342     } else {
 343         av_log(avctx, AV_LOG_ERROR, "Not a HQ/HQA frame.\n");
 344         return AVERROR_INVALIDDATA;
 345     }
 346     if (ret < 0) {
 347         av_log(avctx, AV_LOG_ERROR, "Error decoding frame.\n");
 348         return ret;
 349     }
 350
 351     pic->key_frame = 1;
 352     pic->pict_type = AV_PICTURE_TYPE_I;
 353
 354     *got_frame = 1;
 355
 356     return avpkt->size;
 357 }
 358
 359 static av_cold int hq_hqa_decode_init(AVCodecContext *avctx)
 360 {
 361     HQContext *ctx = avctx->priv_data;
 362     ctx->avctx = avctx;
 363
 364     ff_hqdsp_init(&ctx->hqhqadsp);
 365
 366     return ff_hq_init_vlcs(ctx);
 367 }
 368
 369 static av_cold int hq_hqa_decode_close(AVCodecContext *avctx)
 370 {
 371     HQContext *ctx = avctx->priv_data;
 372
 373     ff_free_vlc(&ctx->hq_ac_vlc);
 374     ff_free_vlc(&ctx->hqa_cbp_vlc);
 375
 376     return 0;
 377 }
 378
 379 AVCodec ff_hq_hqa_decoder = {
 380     .name           = "hq_hqa",
 381     .long_name      = NULL_IF_CONFIG_SMALL("Canopus HQ/HQA"),
 382     .type           = AVMEDIA_TYPE_VIDEO,
 383     .id             = AV_CODEC_ID_HQ_HQA,
 384     .priv_data_size = sizeof(HQContext),
 385     .init           = hq_hqa_decode_init,
 386     .decode         = hq_hqa_decode_frame,
 387     .close          = hq_hqa_decode_close,
 388     .capabilities   = AV_CODEC_CAP_DR1,
 389     .caps_internal  = FF_CODEC_CAP_INIT_THREADSAFE |
 390                       FF_CODEC_CAP_INIT_CLEANUP,
 391 };