summaryrefslogtreecommitdiff
path: root/libav/mpegvideo.c
diff options
context:
space:
mode:
Diffstat (limited to 'libav/mpegvideo.c')
-rw-r--r--libav/mpegvideo.c1098
1 files changed, 0 insertions, 1098 deletions
diff --git a/libav/mpegvideo.c b/libav/mpegvideo.c
deleted file mode 100644
index 4987b38af4..0000000000
--- a/libav/mpegvideo.c
+++ /dev/null
@@ -1,1098 +0,0 @@
-/*
- * The simplest mpeg encoder
- * Copyright (c) 2000 Gerard Lantau.
- *
- * This program is free software; you can redistribute it and/or modify
- * it under the terms of the GNU General Public License as published by
- * the Free Software Foundation; either version 2 of the License, or
- * (at your option) any later version.
- *
- * This program is distributed in the hope that it will be useful,
- * but WITHOUT ANY WARRANTY; without even the implied warranty of
- * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
- * GNU General Public License for more details.
- *
- * You should have received a copy of the GNU General Public License
- * along with this program; if not, write to the Free Software
- * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
- */
-#include <stdlib.h>
-#include <stdio.h>
-#include <netinet/in.h>
-#include <math.h>
-#include "avcodec.h"
-#include "mpegvideo.h"
-
-//#define DEBUG
-
-/* depends on JPEG librarie */
-extern void jpeg_fdct_ifast (DCTELEM * data);
-
-/* depends on mpeg */
-extern void j_rev_dct (DCTELEM *data);
-
-/* for jpeg fast DCT */
-#define CONST_BITS 14
-
-static const unsigned short aanscales[64] = {
- /* precomputed values scaled up by 14 bits */
- 16384, 22725, 21407, 19266, 16384, 12873, 8867, 4520,
- 22725, 31521, 29692, 26722, 22725, 17855, 12299, 6270,
- 21407, 29692, 27969, 25172, 21407, 16819, 11585, 5906,
- 19266, 26722, 25172, 22654, 19266, 15137, 10426, 5315,
- 16384, 22725, 21407, 19266, 16384, 12873, 8867, 4520,
- 12873, 17855, 16819, 15137, 12873, 10114, 6967, 3552,
- 8867, 12299, 11585, 10426, 8867, 6967, 4799, 2446,
- 4520, 6270, 5906, 5315, 4520, 3552, 2446, 1247
-};
-
-static UINT8 cropTbl[256 + 2 * MAX_NEG_CROP];
-static UINT32 squareTbl[512];
-
-static void encode_picture(MpegEncContext *s, int picture_number);
-static void rate_control_init(MpegEncContext *s);
-static int rate_estimate_qscale(MpegEncContext *s);
-static void mpeg1_skip_picture(MpegEncContext *s, int pict_num);
-
-#include "mpegencodevlc.h"
-
-static void put_header(MpegEncContext *s, int header)
-{
- align_put_bits(&s->pb);
- put_bits(&s->pb, 32, header);
-}
-
-static void convert_matrix(int *qmat, const UINT8 *quant_matrix, int qscale)
-{
- int i;
-
- for(i=0;i<64;i++) {
- qmat[i] = (int)((1 << 22) * 16384.0 / (aanscales[i] * qscale * quant_matrix[i]));
- }
-}
-
-
-int MPV_encode_init(AVEncodeContext *avctx)
-{
- MpegEncContext *s = avctx->priv_data;
- int pict_size, c_size, i;
- UINT8 *pict;
-
- s->bit_rate = avctx->bit_rate;
- s->frame_rate = avctx->rate;
- s->width = avctx->width;
- s->height = avctx->height;
- s->gop_size = avctx->gop_size;
- if (s->gop_size <= 1) {
- s->intra_only = 1;
- s->gop_size = 12;
- } else {
- s->intra_only = 0;
- }
-
- switch(avctx->codec->id) {
- case CODEC_ID_MPEG1VIDEO:
- s->out_format = FMT_MPEG1;
- break;
- case CODEC_ID_MJPEG:
- s->out_format = FMT_MJPEG;
- s->intra_only = 1; /* force intra only for jpeg */
- if (mjpeg_init(s) < 0)
- return -1;
- break;
- case CODEC_ID_H263:
- s->out_format = FMT_H263;
- break;
- case CODEC_ID_RV10:
- s->out_format = FMT_H263;
- s->h263_rv10 = 1;
- break;
- default:
- return -1;
- }
-
- switch(s->frame_rate) {
- case 24:
- s->frame_rate_index = 2;
- break;
- case 25:
- s->frame_rate_index = 3;
- break;
- case 30:
- s->frame_rate_index = 5;
- break;
- case 50:
- s->frame_rate_index = 6;
- break;
- case 60:
- s->frame_rate_index = 8;
- break;
- default:
- /* we accept lower frame rates than 24 for low bit rate mpeg */
- if (s->frame_rate >= 1 && s->frame_rate < 24) {
- s->frame_rate_index = 2;
- } else {
- return -1;
- }
- break;
- }
-
- /* init */
- s->mb_width = s->width / 16;
- s->mb_height = s->height / 16;
-
- c_size = s->width * s->height;
- pict_size = (c_size * 3) / 2;
- pict = malloc(pict_size);
- if (pict == NULL)
- return -1;
- s->last_picture[0] = pict;
- s->last_picture[1] = pict + c_size;
- s->last_picture[2] = pict + c_size + (c_size / 4);
-
- pict = malloc(pict_size);
- if (pict == NULL)
- return -1;
- s->last_picture[0] = pict;
- s->last_picture[1] = pict + c_size;
- s->last_picture[2] = pict + c_size + (c_size / 4);
-
- pict = malloc(pict_size);
- if (pict == NULL) {
- free(s->last_picture[0]);
- return -1;
- }
- s->current_picture[0] = pict;
- s->current_picture[1] = pict + c_size;
- s->current_picture[2] = pict + c_size + (c_size / 4);
-
- for(i=0;i<256;i++) cropTbl[i + MAX_NEG_CROP] = i;
- for(i=0;i<MAX_NEG_CROP;i++) {
- cropTbl[i] = 0;
- cropTbl[i + MAX_NEG_CROP + 256] = 255;
- }
-
- for(i=0;i<512;i++) {
- squareTbl[i] = (i - 256) * (i - 256);
- }
-
- /* rate control init */
- rate_control_init(s);
-
- s->picture_number = 0;
- s->fake_picture_number = 0;
-
- return 0;
-}
-
-int MPV_encode_end(AVEncodeContext *avctx)
-{
- MpegEncContext *s = avctx->priv_data;
-#if 0
- /* end of sequence */
- if (s->out_format == FMT_MPEG1) {
- put_header(s, SEQ_END_CODE);
- }
-
- if (!s->flush_frames)
- flush_put_bits(&s->pb);
-#endif
- free(s->last_picture[0]);
- free(s->current_picture[0]);
- if (s->out_format == FMT_MJPEG)
- mjpeg_close(s);
- return 0;
-}
-
-int MPV_encode_picture(AVEncodeContext *avctx,
- unsigned char *buf, int buf_size, void *data)
-{
- MpegEncContext *s = avctx->priv_data;
- int i;
-
- memcpy(s->new_picture, data, 3 * sizeof(UINT8 *));
-
- init_put_bits(&s->pb, buf, buf_size, NULL, NULL);
-
- /* group of picture */
- if (s->out_format == FMT_MPEG1) {
- unsigned int vbv_buffer_size;
- unsigned int time_code, fps, n;
-
- if ((s->picture_number % s->gop_size) == 0) {
- /* mpeg1 header repeated every gop */
- put_header(s, SEQ_START_CODE);
-
- put_bits(&s->pb, 12, s->width);
- put_bits(&s->pb, 12, s->height);
- put_bits(&s->pb, 4, 1); /* 1/1 aspect ratio */
- put_bits(&s->pb, 4, s->frame_rate_index);
- put_bits(&s->pb, 18, 0x3ffff);
- put_bits(&s->pb, 1, 1); /* marker */
- /* vbv buffer size: slightly greater than an I frame. We add
- some margin just in case */
- vbv_buffer_size = (3 * s->I_frame_bits) / (2 * 8);
- put_bits(&s->pb, 10, (vbv_buffer_size + 16383) / 16384);
- put_bits(&s->pb, 1, 1); /* constrained parameter flag */
- put_bits(&s->pb, 1, 0); /* no custom intra matrix */
- put_bits(&s->pb, 1, 0); /* no custom non intra matrix */
-
- put_header(s, GOP_START_CODE);
- put_bits(&s->pb, 1, 0); /* do drop frame */
- /* time code : we must convert from the real frame rate to a
- fake mpeg frame rate in case of low frame rate */
- fps = frame_rate_tab[s->frame_rate_index];
- time_code = s->fake_picture_number;
- s->gop_picture_number = time_code;
- put_bits(&s->pb, 5, (time_code / (fps * 3600)) % 24);
- put_bits(&s->pb, 6, (time_code / (fps * 60)) % 60);
- put_bits(&s->pb, 1, 1);
- put_bits(&s->pb, 6, (time_code / fps) % 60);
- put_bits(&s->pb, 6, (time_code % fps));
- put_bits(&s->pb, 1, 1); /* closed gop */
- put_bits(&s->pb, 1, 0); /* broken link */
- }
-
- if (s->frame_rate < 24 && s->picture_number > 0) {
- /* insert empty P pictures to slow down to the desired
- frame rate. Each fake pictures takes about 20 bytes */
- fps = frame_rate_tab[s->frame_rate_index];
- n = ((s->picture_number * fps) / s->frame_rate) - 1;
- while (s->fake_picture_number < n) {
- mpeg1_skip_picture(s, s->fake_picture_number -
- s->gop_picture_number);
- s->fake_picture_number++;
- }
-
- }
- s->fake_picture_number++;
- }
-
-
- if (!s->intra_only) {
- /* first picture of GOP is intra */
- if ((s->picture_number % s->gop_size) == 0)
- s->pict_type = I_TYPE;
- else
- s->pict_type = P_TYPE;
- } else {
- s->pict_type = I_TYPE;
- }
- avctx->key_frame = (s->pict_type == I_TYPE);
-
- encode_picture(s, s->picture_number);
-
- /* swap current and last picture */
- for(i=0;i<3;i++) {
- UINT8 *tmp;
-
- tmp = s->last_picture[i];
- s->last_picture[i] = s->current_picture[i];
- s->current_picture[i] = tmp;
- }
- s->picture_number++;
-
- if (s->out_format == FMT_MJPEG)
- mjpeg_picture_trailer(s);
-
- flush_put_bits(&s->pb);
- s->total_bits += (s->pb.buf_ptr - s->pb.buf) * 8;
- return s->pb.buf_ptr - s->pb.buf;
-}
-
-/* insert a fake P picture */
-static void mpeg1_skip_picture(MpegEncContext *s, int pict_num)
-{
- unsigned int mb_incr;
-
- /* mpeg1 picture header */
- put_header(s, PICTURE_START_CODE);
- /* temporal reference */
- put_bits(&s->pb, 10, pict_num & 0x3ff);
-
- put_bits(&s->pb, 3, P_TYPE);
- put_bits(&s->pb, 16, 0xffff); /* non constant bit rate */
-
- put_bits(&s->pb, 1, 1); /* integer coordinates */
- put_bits(&s->pb, 3, 1); /* forward_f_code */
-
- put_bits(&s->pb, 1, 0); /* extra bit picture */
-
- /* only one slice */
- put_header(s, SLICE_MIN_START_CODE);
- put_bits(&s->pb, 5, 1); /* quantizer scale */
- put_bits(&s->pb, 1, 0); /* slice extra information */
-
- mb_incr = 1;
- put_bits(&s->pb, mbAddrIncrTable[mb_incr][1],
- mbAddrIncrTable[mb_incr][0]);
-
- /* empty macroblock */
- put_bits(&s->pb, 3, 1); /* motion only */
-
- /* zero motion x & y */
- put_bits(&s->pb, 1, 1);
- put_bits(&s->pb, 1, 1);
-
- /* output a number of empty slice */
- mb_incr = s->mb_width * s->mb_height - 1;
- while (mb_incr > 33) {
- put_bits(&s->pb, 11, 0x008);
- mb_incr -= 33;
- }
- put_bits(&s->pb, mbAddrIncrTable[mb_incr][1],
- mbAddrIncrTable[mb_incr][0]);
-
- /* empty macroblock */
- put_bits(&s->pb, 3, 1); /* motion only */
-
- /* zero motion x & y */
- put_bits(&s->pb, 1, 1);
- put_bits(&s->pb, 1, 1);
-}
-
-static int pix_sum(UINT8 *pix, int line_size)
-{
- int s, i, j;
-
- s = 0;
- for(i=0;i<16;i++) {
- for(j=0;j<16;j+=8) {
- s += pix[0];
- s += pix[1];
- s += pix[2];
- s += pix[3];
- s += pix[4];
- s += pix[5];
- s += pix[6];
- s += pix[7];
- pix += 8;
- }
- pix += line_size - 16;
- }
- return s;
-}
-
-static int pix_norm1(UINT8 *pix, int line_size)
-{
- int s, i, j;
- UINT32 *sq = squareTbl + 256;
-
- s = 0;
- for(i=0;i<16;i++) {
- for(j=0;j<16;j+=8) {
- s += sq[pix[0]];
- s += sq[pix[1]];
- s += sq[pix[2]];
- s += sq[pix[3]];
- s += sq[pix[4]];
- s += sq[pix[5]];
- s += sq[pix[6]];
- s += sq[pix[7]];
- pix += 8;
- }
- pix += line_size - 16;
- }
- return s;
-}
-
-static int pix_norm(UINT8 *pix1, UINT8 *pix2, int line_size)
-{
- int s, i, j;
- UINT32 *sq = squareTbl + 256;
-
- s = 0;
- for(i=0;i<16;i++) {
- for(j=0;j<16;j+=8) {
- s += sq[pix1[0] - pix2[0]];
- s += sq[pix1[1] - pix2[1]];
- s += sq[pix1[2] - pix2[2]];
- s += sq[pix1[3] - pix2[3]];
- s += sq[pix1[4] - pix2[4]];
- s += sq[pix1[5] - pix2[5]];
- s += sq[pix1[6] - pix2[6]];
- s += sq[pix1[7] - pix2[7]];
- pix1 += 8;
- pix2 += 8;
- }
- pix1 += line_size - 16;
- pix2 += line_size - 16;
- }
- return s;
-}
-
-
-static int estimate_motion(MpegEncContext *s,
- int mb_x, int mb_y,
- int *mx_ptr, int *my_ptr)
-{
- UINT8 *pix, *ppix;
- int sum, varc, vard;
-
- pix = s->new_picture[0] + (mb_y * 16 * s->width) + mb_x * 16;
- ppix = s->last_picture[0] + (mb_y * 16 * s->width) + mb_x * 16;
-
- sum = pix_sum(pix, s->width);
- varc = pix_norm1(pix, s->width);
- vard = pix_norm(pix, ppix, s->width);
-
- vard = vard >> 8;
- sum = sum >> 8;
- varc = (varc >> 8) - sum * sum;
-
- *mx_ptr = 0;
- *my_ptr = 0;
- if (vard <= 64) {
- return 0;
- } else if (vard < varc) {
- return 0;
- } else {
- return 1;
- }
-}
-
-static void get_pixels(DCTELEM *block, const UINT8 *pixels, int line_size);
-static void put_pixels(const DCTELEM *block, UINT8 *pixels, int line_size);
-static void sub_pixels(DCTELEM *block, const UINT8 *pixels, int line_size);
-static void add_pixels(DCTELEM *block, const UINT8 *pixels, int line_size);
-static int dct_quantize(MpegEncContext *s, DCTELEM *block, int qscale);
-static void encode_block(MpegEncContext *s,
- DCTELEM *block,
- int component);
-static void dct_unquantize(MpegEncContext *s, DCTELEM *block, int qscale);
-static void mpeg1_encode_mb(MpegEncContext *s, int mb_x, int mb_y,
- DCTELEM block[6][64],
- int motion_x, int motion_y);
-
-static void encode_picture(MpegEncContext *s, int picture_number)
-{
- int mb_x, mb_y;
- UINT8 *ptr;
- DCTELEM block[6][64];
- int i, motion_x, motion_y;
-
- s->picture_number = picture_number;
- s->qscale = rate_estimate_qscale(s);
-
- /* precompute matrix */
- if (s->out_format == FMT_MJPEG) {
- /* for mjpeg, we do include qscale in the matrix */
- s->init_intra_matrix[0] = default_intra_matrix[0];
- for(i=1;i<64;i++)
- s->init_intra_matrix[i] = (default_intra_matrix[i] * s->qscale) >> 3;
- convert_matrix(s->intra_matrix, s->init_intra_matrix, 8);
- } else {
- convert_matrix(s->intra_matrix, default_intra_matrix, s->qscale);
- convert_matrix(s->non_intra_matrix, default_non_intra_matrix, s->qscale);
- }
-
- switch(s->out_format) {
- case FMT_MJPEG:
- mjpeg_picture_header(s);
- break;
- case FMT_H263:
- if (s->h263_rv10)
- rv10_encode_picture_header(s, picture_number);
- else
- h263_picture_header(s, picture_number);
- break;
- case FMT_MPEG1:
- /* mpeg1 picture header */
- put_header(s, PICTURE_START_CODE);
- /* temporal reference */
- put_bits(&s->pb, 10, (s->fake_picture_number -
- s->gop_picture_number) & 0x3ff);
-
- put_bits(&s->pb, 3, s->pict_type);
- put_bits(&s->pb, 16, 0xffff); /* non constant bit rate */
-
- if (s->pict_type == P_TYPE) {
- put_bits(&s->pb, 1, 1); /* integer coordinates */
- put_bits(&s->pb, 3, 1); /* forward_f_code */
- }
-
- put_bits(&s->pb, 1, 0); /* extra bit picture */
-
- /* only one slice */
- put_header(s, SLICE_MIN_START_CODE);
- put_bits(&s->pb, 5, s->qscale); /* quantizer scale */
- put_bits(&s->pb, 1, 0); /* slice extra information */
- break;
- }
-
- /* init last dc values */
- /* XXX: quant matrix value is implied here */
- s->last_dc[0] = 128;
- s->last_dc[1] = 128;
- s->last_dc[2] = 128;
- s->mb_incr = 1;
-
- for(mb_y=0; mb_y < s->mb_height; mb_y++) {
- for(mb_x=0; mb_x < s->mb_width; mb_x++) {
- /* compute motion vector and macro block type (intra or non intra) */
- motion_x = 0;
- motion_y = 0;
- if (s->pict_type == P_TYPE) {
- s->mb_intra = estimate_motion(s, mb_x, mb_y,
- &motion_x,
- &motion_y);
- } else {
- s->mb_intra = 1;
- }
-
- /* reset intra predictors if non intra mb */
- if (!s->mb_intra) {
- s->last_dc[0] = 128;
- s->last_dc[1] = 128;
- s->last_dc[2] = 128;
- }
-
- /* get the pixels */
- ptr = s->new_picture[0] + (mb_y * 16 * s->width) + mb_x * 16;
- get_pixels(block[0], ptr, s->width);
- get_pixels(block[1], ptr + 8, s->width);
- get_pixels(block[2], ptr + 8 * s->width, s->width);
- get_pixels(block[3], ptr + 8 * s->width + 8, s->width);
- ptr = s->new_picture[1] + (mb_y * 8 * (s->width >> 1)) + mb_x * 8;
- get_pixels(block[4],ptr, s->width >> 1);
-
- ptr = s->new_picture[2] + (mb_y * 8 * (s->width >> 1)) + mb_x * 8;
- get_pixels(block[5],ptr, s->width >> 1);
-
- /* subtract previous frame if non intra */
- if (!s->mb_intra) {
- ptr = s->last_picture[0] +
- ((mb_y * 16 + motion_y) * s->width) + (mb_x * 16 + motion_x);
-
- sub_pixels(block[0], ptr, s->width);
- sub_pixels(block[1], ptr + 8, s->width);
- sub_pixels(block[2], ptr + s->width * 8, s->width);
- sub_pixels(block[3], ptr + 8 + s->width * 8, s->width);
- ptr = s->last_picture[1] +
- ((mb_y * 8 + (motion_y >> 1)) * (s->width >> 1)) +
- (mb_x * 8 + (motion_x >> 1));
- sub_pixels(block[4], ptr, s->width >> 1);
- ptr = s->last_picture[2] +
- ((mb_y * 8 + (motion_y >> 1)) * (s->width >> 1)) +
- (mb_x * 8 + (motion_x >> 1));
- sub_pixels(block[5], ptr, s->width >> 1);
- }
-
- /* DCT & quantize */
- for(i=0;i<6;i++) {
- int last_index;
- last_index = dct_quantize(s, block[i], s->qscale);
- s->block_last_index[i] = last_index;
- }
-
- /* huffman encode */
- switch(s->out_format) {
- case FMT_MPEG1:
- mpeg1_encode_mb(s, mb_x, mb_y, block, motion_x, motion_y);
- break;
- case FMT_H263:
- h263_encode_mb(s, block, motion_x, motion_y);
- break;
- case FMT_MJPEG:
- mjpeg_encode_mb(s, block);
- break;
- }
-
- /* decompress blocks so that we keep the state of the decoder */
- if (!s->intra_only) {
- for(i=0;i<6;i++) {
- if (s->block_last_index[i] >= 0) {
- dct_unquantize(s, block[i], s->qscale);
- }
- }
-
- if (!s->mb_intra) {
- ptr = s->last_picture[0] +
- ((mb_y * 16 + motion_y) * s->width) + (mb_x * 16 + motion_x);
-
- add_pixels(block[0], ptr, s->width);
- add_pixels(block[1], ptr + 8, s->width);
- add_pixels(block[2], ptr + s->width * 8, s->width);
- add_pixels(block[3], ptr + 8 + s->width * 8, s->width);
- ptr = s->last_picture[1] +
- ((mb_y * 8 + (motion_y >> 1)) * (s->width >> 1)) +
- (mb_x * 8 + (motion_x >> 1));
- add_pixels(block[4], ptr, s->width >> 1);
- ptr = s->last_picture[2] +
- ((mb_y * 8 + (motion_y >> 1)) * (s->width >> 1)) +
- (mb_x * 8 + (motion_x >> 1));
- add_pixels(block[5], ptr, s->width >> 1);
- }
-
- /* write the pixels */
- ptr = s->current_picture[0] + (mb_y * 16 * s->width) + mb_x * 16;
- put_pixels(block[0], ptr, s->width);
- put_pixels(block[1], ptr + 8, s->width);
- put_pixels(block[2], ptr + 8 * s->width, s->width);
- put_pixels(block[3], ptr + 8 * s->width + 8, s->width);
- ptr = s->current_picture[1] + (mb_y * 8 * (s->width >> 1)) + mb_x * 8;
- put_pixels(block[4],ptr, s->width >> 1);
-
- ptr = s->current_picture[2] + (mb_y * 8 * (s->width >> 1)) + mb_x * 8;
- put_pixels(block[5],ptr, s->width >> 1);
- }
- }
- }
-}
-
-static void mpeg1_encode_mb(MpegEncContext *s, int mb_x, int mb_y,
- DCTELEM block[6][64],
- int motion_x, int motion_y)
-{
- int mb_incr, i, cbp;
-
- /* compute cbp */
- cbp = 0;
- for(i=0;i<6;i++) {
- if (s->block_last_index[i] >= 0)
- cbp |= 1 << (5 - i);
- }
-
- /* skip macroblock, except if first or last macroblock of a slice */
- if ((cbp | motion_x | motion_y) == 0 &&
- (!((mb_x | mb_y) == 0 ||
- (mb_x == s->mb_width - 1 && mb_y == s->mb_height - 1)))) {
- s->mb_incr++;
- } else {
- /* output mb incr */
- mb_incr = s->mb_incr;
-
- while (mb_incr > 33) {
- put_bits(&s->pb, 11, 0x008);
- mb_incr -= 33;
- }
- put_bits(&s->pb, mbAddrIncrTable[mb_incr][1],
- mbAddrIncrTable[mb_incr][0]);
-
- if (s->pict_type == I_TYPE) {
- put_bits(&s->pb, 1, 1); /* macroblock_type : macroblock_quant = 0 */
- } else {
- if (s->mb_intra) {
- put_bits(&s->pb, 5, 0x03);
- } else {
- if (motion_x == 0 && motion_y == 0) {
- if (cbp != 0) {
- put_bits(&s->pb, 2, 1); /* macroblock_pattern only */
- put_bits(&s->pb, mbPatTable[cbp][1], mbPatTable[cbp][0]);
- } else {
- put_bits(&s->pb, 3, 1); /* motion only & zero motion vectors */
- /* zero motion x & y */
- put_bits(&s->pb, 1, 1);
- put_bits(&s->pb, 1, 1);
- }
- } else {
- /* XXX: not used yet */
- put_bits(&s->pb, mbPatTable[cbp][1], mbPatTable[cbp][0]);
- }
- }
-
- }
-
- for(i=0;i<6;i++) {
- if (cbp & (1 << (5 - i))) {
- encode_block(s, block[i], i);
- }
- }
- s->mb_incr = 1;
- }
-}
-
-static void get_pixels(DCTELEM *block, const UINT8 *pixels, int line_size)
-{
- DCTELEM *p;
- const UINT8 *pix;
- int i;
-
- /* read the pixels */
- p = block;
- pix = pixels;
- for(i=0;i<8;i++) {
- p[0] = pix[0];
- p[1] = pix[1];
- p[2] = pix[2];
- p[3] = pix[3];
- p[4] = pix[4];
- p[5] = pix[5];
- p[6] = pix[6];
- p[7] = pix[7];
- pix += line_size;
- p += 8;
- }
-}
-
-static void put_pixels(const DCTELEM *block, UINT8 *pixels, int line_size)
-{
- const DCTELEM *p;
- UINT8 *pix;
- int i;
- UINT8 *cm = cropTbl + MAX_NEG_CROP;
-
- /* read the pixels */
- p = block;
- pix = pixels;
- for(i=0;i<8;i++) {
- pix[0] = cm[p[0]];
- pix[1] = cm[p[1]];
- pix[2] = cm[p[2]];
- pix[3] = cm[p[3]];
- pix[4] = cm[p[4]];
- pix[5] = cm[p[5]];
- pix[6] = cm[p[6]];
- pix[7] = cm[p[7]];
- pix += line_size;
- p += 8;
- }
-}
-
-static void sub_pixels(DCTELEM *block, const UINT8 *pixels, int line_size)
-{
- DCTELEM *p;
- const UINT8 *pix;
- int i;
-
- /* read the pixels */
- p = block;
- pix = pixels;
- for(i=0;i<8;i++) {
- p[0] -= pix[0];
- p[1] -= pix[1];
- p[2] -= pix[2];
- p[3] -= pix[3];
- p[4] -= pix[4];
- p[5] -= pix[5];
- p[6] -= pix[6];
- p[7] -= pix[7];
- pix += line_size;
- p += 8;
- }
-}
-
-static void add_pixels(DCTELEM *block, const UINT8 *pixels, int line_size)
-{
- DCTELEM *p;
- const UINT8 *pix;
- int i;
-
- /* read the pixels */
- p = block;
- pix = pixels;
- for(i=0;i<8;i++) {
- p[0] += pix[0];
- p[1] += pix[1];
- p[2] += pix[2];
- p[3] += pix[3];
- p[4] += pix[4];
- p[5] += pix[5];
- p[6] += pix[6];
- p[7] += pix[7];
- pix += line_size;
- p += 8;
- }
-}
-
-#define USE_FAST_MUL
-
-static int dct_quantize(MpegEncContext *s,
- DCTELEM *block,
- int qscale)
-{
- int i, j, level, last_non_zero;
-#ifdef USE_FAST_MUL
- const int *qmat;
-#else
- const UINT8 *qmat;
-#endif
-
- jpeg_fdct_ifast (block);
-
- if (s->mb_intra) {
- block[0] = (block[0] + 4 * 8) >> 6;
- i = 1;
- last_non_zero = 0;
- if (s->out_format == FMT_H263) {
-#ifdef USE_FAST_MUL
- qmat = s->non_intra_matrix;
-#else
- qmat = default_non_intra_matrix;
-#endif
- } else {
-#ifdef USE_FAST_MUL
- qmat = s->intra_matrix;
-#else
- qmat = default_intra_matrix;
-#endif
- }
- } else {
- i = 0;
- last_non_zero = -1;
-#ifdef USE_FAST_MUL
- qmat = s->non_intra_matrix;
-#else
- qmat = default_non_intra_matrix;
-#endif
- }
-
- for(;i<64;i++) {
- j = zigzag_direct[i];
- level = block[j];
-#ifdef USE_FAST_MUL
- level = (level * qmat[j]) / (1 << 22);
-#else
- /* post dct normalization */
- level = (level << 11) / aanscales[j];
- /* quantification */
- level = (8 * level) / (qscale * qmat[j]);
-#endif
- block[j] = level;
- if (level)
- last_non_zero = i;
- }
- return last_non_zero;
-}
-
-static void dct_unquantize(MpegEncContext *s,
- DCTELEM *block, int qscale)
-{
- int i, level, coeff;
- const UINT8 *quant_matrix;
-
- if (s->mb_intra) {
- block[0] = block[0] << 3;
- if (s->out_format == FMT_H263) {
- i = 1;
- goto unquant_even;
- }
- quant_matrix = default_intra_matrix;
- for(i=1;i<64;i++) {
- block[i] = (block[i] * qscale * quant_matrix[i]) >> 3;
- }
- } else {
- i = 0;
- unquant_even:
- quant_matrix = default_non_intra_matrix;
- for(;i<64;i++) {
- level = block[i];
- if (level) {
- if (level < 0) {
- coeff = (((level << 1) - 1) * qscale *
- ((int) (quant_matrix[i]))) >> 4;
- coeff += (coeff & 1);
- } else {
- coeff = (((level << 1) + 1) * qscale *
- ((int) (quant_matrix[i]))) >> 4;
- coeff -= (coeff & 1);
- }
- block[i] = coeff;
- }
- }
- }
-
- j_rev_dct(block);
-}
-
-
-static inline void encode_dc(MpegEncContext *s, int diff, int component)
-{
- int adiff, index;
-
- // printf("dc=%d c=%d\n", diff, component);
- adiff = abs(diff);
- index = vlc_dc_table[adiff];
- if (component == 0) {
- put_bits(&s->pb, vlc_dc_lum_bits[index], vlc_dc_lum_code[index]);
- } else {
- put_bits(&s->pb, vlc_dc_chroma_bits[index], vlc_dc_chroma_code[index]);
- }
- if (diff > 0) {
- put_bits(&s->pb, index, (diff & ((1 << index) - 1)));
- } else if (diff < 0) {
- put_bits(&s->pb, index, ((diff - 1) & ((1 << index) - 1)));
- }
-}
-
-static void encode_block(MpegEncContext *s,
- DCTELEM *block,
- int n)
-{
- int alevel, level, last_non_zero, dc, diff, i, j, run, last_index;
- int code, nbits, component;
-
- last_index = s->block_last_index[n];
-
- /* DC coef */
- if (s->mb_intra) {
- component = (n <= 3 ? 0 : n - 4 + 1);
- dc = block[0]; /* overflow is impossible */
- diff = dc - s->last_dc[component];
- encode_dc(s, diff, component);
- s->last_dc[component] = dc;
- i = 1;
- } else {
- /* encode the first coefficient : needs to be done here because
- it is handled slightly differently */
- level = block[0];
- if (abs(level) == 1) {
- code = ((UINT32)level >> 31); /* the sign bit */
- put_bits(&s->pb, 2, code | 0x02);
- i = 1;
- } else {
- i = 0;
- last_non_zero = -1;
- goto next_coef;
- }
- }
-
- /* now quantify & encode AC coefs */
- last_non_zero = i - 1;
- for(;i<=last_index;i++) {
- j = zigzag_direct[i];
- level = block[j];
- next_coef:
-#if 0
- if (level != 0)
- printf("level[%d]=%d\n", i, level);
-#endif
- /* encode using VLC */
- if (level != 0) {
- run = i - last_non_zero - 1;
- alevel = abs(level);
- // printf("run=%d level=%d\n", run, level);
- if ( (run < HUFF_MAXRUN) && (alevel < huff_maxlevel[run])) {
- /* encode using the Huffman tables */
- code = (huff_table[run])[alevel];
- nbits = (huff_bits[run])[alevel];
- code |= ((UINT32)level >> 31); /* the sign bit */
-
- put_bits(&s->pb, nbits, code);
- } else {
- /* escape: only clip in this case */
- if (level > 255)
- level = 255;
- else if (level < -255)
- level = -255;
- put_bits(&s->pb, 6, 0x1);
- put_bits(&s->pb, 6, run);
- if (alevel < 128) {
- put_bits(&s->pb, 8, level & 0xff);
- } else {
- if (level < 0) {
- put_bits(&s->pb, 16, 0x8001 + level + 255);
- } else {
- put_bits(&s->pb, 16, level & 0xffff);
- }
- }
- }
- last_non_zero = i;
- }
- }
- /* end of block */
- put_bits(&s->pb, 2, 0x2);
-}
-
-
-/* rate control */
-
-/* an I frame is I_FRAME_SIZE_RATIO bigger than a P frame */
-#define I_FRAME_SIZE_RATIO 1.5
-#define QSCALE_K 20
-
-static void rate_control_init(MpegEncContext *s)
-{
- s->wanted_bits = 0;
-
- if (s->intra_only) {
- s->I_frame_bits = s->bit_rate / s->frame_rate;
- s->P_frame_bits = s->I_frame_bits;
- } else {
- s->P_frame_bits = (int) ((float)(s->gop_size * s->bit_rate) /
- (float)(s->frame_rate * (I_FRAME_SIZE_RATIO + s->gop_size - 1)));
- s->I_frame_bits = (int)(s->P_frame_bits * I_FRAME_SIZE_RATIO);
- }
-
-#if defined(DEBUG)
- printf("I_frame_size=%d P_frame_size=%d\n",
- s->I_frame_bits, s->P_frame_bits);
-#endif
-}
-
-
-/*
- * This heuristic is rather poor, but at least we do not have to
- * change the qscale at every macroblock.
- */
-static int rate_estimate_qscale(MpegEncContext *s)
-{
- long long total_bits = s->total_bits;
- float q;
- int qscale, diff;
-
- if (s->pict_type == I_TYPE) {
- s->wanted_bits += s->I_frame_bits;
- } else {
- s->wanted_bits += s->P_frame_bits;
- }
- diff = s->wanted_bits - total_bits;
- q = 31.0 - (float)diff / (QSCALE_K * s->mb_height * s->mb_width);
- /* adjust for I frame */
- if (s->pict_type == I_TYPE && !s->intra_only) {
- q /= I_FRAME_SIZE_RATIO;
- }
-
- if (q < 1)
- q = 1;
- else if (q > 31)
- q = 31;
- qscale = (int)(q + 0.5);
-#if defined(DEBUG)
- printf("%d: total=%Ld br=%0.1f diff=%d qest=%0.1f\n",
- s->picture_number,
- total_bits, (float)s->frame_rate * total_bits / s->picture_number,
- diff, q);
-#endif
- return qscale;
-}
-
-AVEncoder mpeg1video_encoder = {
- "mpeg1video",
- CODEC_TYPE_VIDEO,
- CODEC_ID_MPEG1VIDEO,
- sizeof(MpegEncContext),
- MPV_encode_init,
- MPV_encode_picture,
- MPV_encode_end,
-};
-
-AVEncoder h263_encoder = {
- "h263",
- CODEC_TYPE_VIDEO,
- CODEC_ID_H263,
- sizeof(MpegEncContext),
- MPV_encode_init,
- MPV_encode_picture,
- MPV_encode_end,
-};
-
-AVEncoder rv10_encoder = {
- "rv10",
- CODEC_TYPE_VIDEO,
- CODEC_ID_RV10,
- sizeof(MpegEncContext),
- MPV_encode_init,
- MPV_encode_picture,
- MPV_encode_end,
-};
-
-AVEncoder mjpeg_encoder = {
- "mjpeg",
- CODEC_TYPE_VIDEO,
- CODEC_ID_MJPEG,
- sizeof(MpegEncContext),
- MPV_encode_init,
- MPV_encode_picture,
- MPV_encode_end,
-};