FFmpeg
 All Data Structures Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Groups Pages
psymodel.c
Go to the documentation of this file.
1 /*
2  * audio encoder psychoacoustic model
3  * Copyright (C) 2008 Konstantin Shishkov
4  *
5  * This file is part of FFmpeg.
6  *
7  * FFmpeg 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  * FFmpeg 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 FFmpeg; if not, write to the Free Software
19  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
20  */
21 
22 #include <string.h>
23 
24 #include "avcodec.h"
25 #include "psymodel.h"
26 #include "iirfilter.h"
27 #include "libavutil/mem.h"
28 
29 extern const FFPsyModel ff_aac_psy_model;
30 
31 av_cold int ff_psy_init(FFPsyContext *ctx, AVCodecContext *avctx, int num_lens,
32  const uint8_t **bands, const int* num_bands,
33  int num_groups, const uint8_t *group_map)
34 {
35  int i, j, k = 0;
36 
37  ctx->avctx = avctx;
38  ctx->ch = av_mallocz_array(sizeof(ctx->ch[0]), avctx->channels * 2);
39  ctx->group = av_mallocz_array(sizeof(ctx->group[0]), num_groups);
40  ctx->bands = av_malloc_array (sizeof(ctx->bands[0]), num_lens);
41  ctx->num_bands = av_malloc_array (sizeof(ctx->num_bands[0]), num_lens);
42  ctx->cutoff = avctx->cutoff;
43 
44  if (!ctx->ch || !ctx->group || !ctx->bands || !ctx->num_bands) {
45  ff_psy_end(ctx);
46  return AVERROR(ENOMEM);
47  }
48 
49  memcpy(ctx->bands, bands, sizeof(ctx->bands[0]) * num_lens);
50  memcpy(ctx->num_bands, num_bands, sizeof(ctx->num_bands[0]) * num_lens);
51 
52  /* assign channels to groups (with virtual channels for coupling) */
53  for (i = 0; i < num_groups; i++) {
54  /* NOTE: Add 1 to handle the AAC chan_config without modification.
55  * This has the side effect of allowing an array of 0s to map
56  * to one channel per group.
57  */
58  ctx->group[i].num_ch = group_map[i] + 1;
59  for (j = 0; j < ctx->group[i].num_ch * 2; j++)
60  ctx->group[i].ch[j] = &ctx->ch[k++];
61  }
62 
63  switch (ctx->avctx->codec_id) {
64  case AV_CODEC_ID_AAC:
65  ctx->model = &ff_aac_psy_model;
66  break;
67  }
68  if (ctx->model->init)
69  return ctx->model->init(ctx);
70  return 0;
71 }
72 
74 {
75  int i = 0, ch = 0;
76 
77  while (ch <= channel)
78  ch += ctx->group[i++].num_ch;
79 
80  return &ctx->group[i-1];
81 }
82 
84 {
85  if (ctx->model && ctx->model->end)
86  ctx->model->end(ctx);
87  av_freep(&ctx->bands);
88  av_freep(&ctx->num_bands);
89  av_freep(&ctx->group);
90  av_freep(&ctx->ch);
91 }
92 
93 typedef struct FFPsyPreprocessContext{
95  float stereo_att;
100 
101 #define FILT_ORDER 4
102 
104 {
106  int i;
107  float cutoff_coeff = 0;
108  ctx = av_mallocz(sizeof(FFPsyPreprocessContext));
109  if (!ctx)
110  return NULL;
111  ctx->avctx = avctx;
112 
113  /* AAC has its own LP method */
114  if (avctx->codec_id != AV_CODEC_ID_AAC) {
115  if (avctx->cutoff > 0)
116  cutoff_coeff = 2.0 * avctx->cutoff / avctx->sample_rate;
117 
118  if (cutoff_coeff && cutoff_coeff < 0.98)
121  cutoff_coeff, 0.0, 0.0);
122  if (ctx->fcoeffs) {
123  ctx->fstate = av_mallocz(sizeof(ctx->fstate[0]) * avctx->channels);
124  for (i = 0; i < avctx->channels; i++)
126  }
127  }
128 
129  ff_iir_filter_init(&ctx->fiir);
130 
131  return ctx;
132 }
133 
134 void ff_psy_preprocess(struct FFPsyPreprocessContext *ctx, float **audio, int channels)
135 {
136  int ch;
137  int frame_size = ctx->avctx->frame_size;
138  FFIIRFilterContext *iir = &ctx->fiir;
139 
140  if (ctx->fstate) {
141  for (ch = 0; ch < channels; ch++)
142  iir->filter_flt(ctx->fcoeffs, ctx->fstate[ch], frame_size,
143  &audio[ch][frame_size], 1, &audio[ch][frame_size], 1);
144  }
145 }
146 
148 {
149  int i;
151  if (ctx->fstate)
152  for (i = 0; i < ctx->avctx->channels; i++)
154  av_freep(&ctx->fstate);
155  av_free(ctx);
156 }
#define NULL
Definition: coverity.c:32
void(* filter_flt)(const struct FFIIRFilterCoeffs *coeffs, struct FFIIRFilterState *state, int size, const float *src, int sstep, float *dst, int dstep)
Perform IIR filtering on floating-point input samples.
Definition: iirfilter.h:62
uint8_t ** bands
scalefactor band sizes for possible frame sizes
Definition: psymodel.h:98
AVFormatContext * ctx
Definition: movenc-test.c:48
FFPsyChannelGroup * group
channel group information
Definition: psymodel.h:94
memory handling functions
av_cold void ff_psy_preprocess_end(struct FFPsyPreprocessContext *ctx)
Cleanup audio preprocessing module.
Definition: psymodel.c:147
psychoacoustic information for an arbitrary group of channels
Definition: psymodel.h:68
av_cold int ff_psy_init(FFPsyContext *ctx, AVCodecContext *avctx, int num_lens, const uint8_t **bands, const int *num_bands, int num_groups, const uint8_t *group_map)
Initialize psychoacoustic model.
Definition: psymodel.c:31
int * num_bands
number of scalefactor bands for possible frame sizes
Definition: psymodel.h:99
av_cold struct FFIIRFilterState * ff_iir_filter_init_state(int order)
Create new filter state.
Definition: iirfilter.c:203
av_cold struct FFIIRFilterCoeffs * ff_iir_filter_init_coeffs(void *avc, enum IIRFilterType filt_type, enum IIRFilterMode filt_mode, int order, float cutoff_ratio, float stopband, float ripple)
Initialize filter coefficients.
Definition: iirfilter.c:161
uint8_t
#define av_cold
Definition: attributes.h:82
struct FFIIRFilterCoeffs * fcoeffs
Definition: psymodel.c:96
context used by psychoacoustic model
Definition: psymodel.h:89
const FFPsyModel ff_aac_psy_model
Definition: aacpsy.c:1047
AVCodecContext * avctx
Definition: psymodel.c:94
FFPsyChannel * ch[PSY_MAX_CHANS]
pointers to the individual channels in the group
Definition: psymodel.h:69
#define AVERROR(e)
Definition: error.h:43
av_cold void ff_iir_filter_free_statep(struct FFIIRFilterState **state)
Free and zero filter state.
Definition: iirfilter.c:302
codec-specific psychoacoustic model implementation
Definition: psymodel.h:114
IIR filter state.
Definition: iirfilter.c:45
uint8_t num_ch
number of channels in this group
Definition: psymodel.h:70
int frame_size
Number of samples per channel in an audio frame.
Definition: avcodec.h:2307
int frame_size
Definition: mxfenc.c:1821
Libavcodec external API header.
enum AVCodecID codec_id
Definition: avcodec.h:1549
int sample_rate
samples per second
Definition: avcodec.h:2287
FFPsyChannelGroup * ff_psy_find_group(FFPsyContext *ctx, int channel)
Determine what group a channel belongs to.
Definition: psymodel.c:73
main external API structure.
Definition: avcodec.h:1532
void(* end)(FFPsyContext *apc)
Definition: psymodel.h:141
int cutoff
lowpass frequency cutoff for analysis
Definition: psymodel.h:96
const struct FFPsyModel * model
encoder-specific model functions
Definition: psymodel.h:91
IIR filter global parameters.
Definition: iirfilter.c:35
void ff_iir_filter_init(FFIIRFilterContext *f)
Initialize FFIIRFilterContext.
Definition: iirfilter.c:317
struct FFIIRFilterState ** fstate
Definition: psymodel.c:97
av_cold struct FFPsyPreprocessContext * ff_psy_preprocess_init(AVCodecContext *avctx)
psychoacoustic model audio preprocessing initialization
Definition: psymodel.c:103
void ff_psy_preprocess(struct FFPsyPreprocessContext *ctx, float **audio, int channels)
Preprocess several channel in audio frame in order to compress it better.
Definition: psymodel.c:134
#define FILT_ORDER
Definition: psymodel.c:101
struct FFIIRFilterContext fiir
Definition: psymodel.c:98
int cutoff
Audio cutoff bandwidth (0 means "automatic")
Definition: avcodec.h:2331
#define av_free(p)
IIR filter interface.
int channels
number of audio channels
Definition: avcodec.h:2288
FFPsyChannel * ch
single channel information
Definition: psymodel.h:93
int(* init)(FFPsyContext *apc)
Definition: psymodel.h:116
static void * av_mallocz_array(size_t nmemb, size_t size)
Definition: mem.h:229
#define av_freep(p)
av_cold void ff_iir_filter_free_coeffsp(struct FFIIRFilterCoeffs **coeffsp)
Free filter coefficients.
Definition: iirfilter.c:307
#define av_malloc_array(a, b)
AVCodecContext * avctx
encoder context
Definition: psymodel.h:90
av_cold void ff_psy_end(FFPsyContext *ctx)
Cleanup model context at the end.
Definition: psymodel.c:83
void * av_mallocz(size_t size)
Allocate a block of size bytes with alignment suitable for all memory accesses (including vectors if ...
Definition: mem.c:252