110 #define LICENSE_PREFIX "libpostproc license: "
114 #define GET_MODE_BUFFER_SIZE 500
115 #define OPTIONS_ARRAY_SIZE 10
117 #define TEMP_STRIDE 8
120 #if ARCH_X86 && HAVE_INLINE_ASM
144 {
"dr",
"dering", 1, 5, 6,
DERING},
145 {
"al",
"autolevels", 0, 1, 2,
LEVEL_FIX},
154 {
"be",
"bitexact", 1, 0, 0,
BITEXACT},
161 "default",
"hb:a,vb:a,dr:a",
162 "de",
"hb:a,vb:a,dr:a",
163 "fast",
"h1:a,v1:a,dr:a",
164 "fa",
"h1:a,v1:a,dr:a",
165 "ac",
"ha:a:128:7,va:a,dr:a",
179 const int dcOffset= ((
c->nonBQP*
c->ppMode.baseDcDiff)>>8) + 1;
180 const int dcThreshold= dcOffset*2 + 1;
183 numEq += ((unsigned)(
src[0] -
src[1] + dcOffset)) < dcThreshold;
184 numEq += ((unsigned)(
src[1] -
src[2] + dcOffset)) < dcThreshold;
185 numEq += ((unsigned)(
src[2] -
src[3] + dcOffset)) < dcThreshold;
186 numEq += ((unsigned)(
src[3] -
src[4] + dcOffset)) < dcThreshold;
187 numEq += ((unsigned)(
src[4] -
src[5] + dcOffset)) < dcThreshold;
188 numEq += ((unsigned)(
src[5] -
src[6] + dcOffset)) < dcThreshold;
189 numEq += ((unsigned)(
src[6] -
src[7] + dcOffset)) < dcThreshold;
192 return numEq >
c->ppMode.flatnessThreshold;
202 const int dcOffset= ((
c->nonBQP*
c->ppMode.baseDcDiff)>>8) + 1;
203 const int dcThreshold= dcOffset*2 + 1;
207 numEq += ((unsigned)(
src[0] -
src[0+
stride] + dcOffset)) < dcThreshold;
208 numEq += ((unsigned)(
src[1] -
src[1+
stride] + dcOffset)) < dcThreshold;
209 numEq += ((unsigned)(
src[2] -
src[2+
stride] + dcOffset)) < dcThreshold;
210 numEq += ((unsigned)(
src[3] -
src[3+
stride] + dcOffset)) < dcThreshold;
211 numEq += ((unsigned)(
src[4] -
src[4+
stride] + dcOffset)) < dcThreshold;
212 numEq += ((unsigned)(
src[5] -
src[5+
stride] + dcOffset)) < dcThreshold;
213 numEq += ((unsigned)(
src[6] -
src[6+
stride] + dcOffset)) < dcThreshold;
214 numEq += ((unsigned)(
src[7] -
src[7+
stride] + dcOffset)) < dcThreshold;
217 return numEq >
c->ppMode.flatnessThreshold;
224 if((
unsigned)(
src[0] -
src[5] + 2*
QP) > 4*
QP)
return 0;
226 if((
unsigned)(
src[2] -
src[7] + 2*
QP) > 4*
QP)
return 0;
228 if((
unsigned)(
src[4] -
src[1] + 2*
QP) > 4*
QP)
return 0;
230 if((
unsigned)(
src[6] -
src[3] + 2*
QP) > 4*
QP)
return 0;
271 const int middleEnergy= 5*(dst[4] - dst[3]) + 2*(dst[2] - dst[5]);
273 if(
FFABS(middleEnergy) < 8*
c->QP){
274 const int q=(dst[3] - dst[4])/2;
275 const int leftEnergy= 5*(dst[2] - dst[1]) + 2*(dst[0] - dst[3]);
276 const int rightEnergy= 5*(dst[6] - dst[5]) + 2*(dst[4] - dst[7]);
282 d*=
FFSIGN(-middleEnergy);
310 const int first=
FFABS(dst[-1] - dst[0]) <
c->QP ? dst[-1] : dst[0];
311 const int last=
FFABS(dst[8] - dst[7]) <
c->QP ? dst[8] : dst[7];
314 sums[0] = 4*first + dst[0] + dst[1] + dst[2] + 4;
315 sums[1] = sums[0] - first + dst[3];
316 sums[2] = sums[1] - first + dst[4];
317 sums[3] = sums[2] - first + dst[5];
318 sums[4] = sums[3] - first + dst[6];
319 sums[5] = sums[4] - dst[0] + dst[7];
320 sums[6] = sums[5] - dst[1] + last;
321 sums[7] = sums[6] - dst[2] + last;
322 sums[8] = sums[7] - dst[3] + last;
323 sums[9] = sums[8] - dst[4] + last;
325 dst[0]= (sums[0] + sums[2] + 2*dst[0])>>4;
326 dst[1]= (sums[1] + sums[3] + 2*dst[1])>>4;
327 dst[2]= (sums[2] + sums[4] + 2*dst[2])>>4;
328 dst[3]= (sums[3] + sums[5] + 2*dst[3])>>4;
329 dst[4]= (sums[4] + sums[6] + 2*dst[4])>>4;
330 dst[5]= (sums[5] + sums[7] + 2*dst[5])>>4;
331 dst[6]= (sums[6] + sums[8] + 2*dst[6])>>4;
332 dst[7]= (sums[7] + sums[9] + 2*dst[7])>>4;
349 static uint64_t lut[256];
355 int v=
i < 128 ? 2*
i : 2*(
i-256);
364 uint64_t
a= (v/16) & 0xFF;
365 uint64_t
b= (v*3/16) & 0xFF;
366 uint64_t
c= (v*5/16) & 0xFF;
367 uint64_t d= (7*v/16) & 0xFF;
368 uint64_t
A= (0x100 -
a)&0xFF;
369 uint64_t
B= (0x100 -
b)&0xFF;
370 uint64_t
C= (0x100 -
c)&0xFF;
371 uint64_t
D= (0x100 -
c)&0xFF;
373 lut[
i] = (
a<<56) | (
b<<48) | (
c<<40) | (d<<32) |
374 (
D<<24) | (
C<<16) | (
B<<8) | (
A);
408 const int dcOffset= ((
c->nonBQP*
c->ppMode.baseDcDiff)>>8) + 1;
409 const int dcThreshold= dcOffset*2 + 1;
415 numEq += ((unsigned)(
src[-1*step] -
src[0*step] + dcOffset)) < dcThreshold;
416 numEq += ((unsigned)(
src[ 0*step] -
src[1*step] + dcOffset)) < dcThreshold;
417 numEq += ((unsigned)(
src[ 1*step] -
src[2*step] + dcOffset)) < dcThreshold;
418 numEq += ((unsigned)(
src[ 2*step] -
src[3*step] + dcOffset)) < dcThreshold;
419 numEq += ((unsigned)(
src[ 3*step] -
src[4*step] + dcOffset)) < dcThreshold;
420 numEq += ((unsigned)(
src[ 4*step] -
src[5*step] + dcOffset)) < dcThreshold;
421 numEq += ((unsigned)(
src[ 5*step] -
src[6*step] + dcOffset)) < dcThreshold;
422 numEq += ((unsigned)(
src[ 6*step] -
src[7*step] + dcOffset)) < dcThreshold;
423 numEq += ((unsigned)(
src[ 7*step] -
src[8*step] + dcOffset)) < dcThreshold;
424 if(numEq >
c->ppMode.flatnessThreshold){
435 if(
src[x*step] >
src[(x+1)*step]){
448 sums[0] = 4*first +
src[0*step] +
src[1*step] +
src[2*step] + 4;
449 sums[1] = sums[0] - first +
src[3*step];
450 sums[2] = sums[1] - first +
src[4*step];
451 sums[3] = sums[2] - first +
src[5*step];
452 sums[4] = sums[3] - first +
src[6*step];
453 sums[5] = sums[4] -
src[0*step] +
src[7*step];
454 sums[6] = sums[5] -
src[1*step] + last;
455 sums[7] = sums[6] -
src[2*step] + last;
456 sums[8] = sums[7] -
src[3*step] + last;
457 sums[9] = sums[8] -
src[4*step] + last;
469 src[0*step]= (sums[0] + sums[2] + 2*
src[0*step])>>4;
470 src[1*step]= (sums[1] + sums[3] + 2*
src[1*step])>>4;
471 src[2*step]= (sums[2] + sums[4] + 2*
src[2*step])>>4;
472 src[3*step]= (sums[3] + sums[5] + 2*
src[3*step])>>4;
473 src[4*step]= (sums[4] + sums[6] + 2*
src[4*step])>>4;
474 src[5*step]= (sums[5] + sums[7] + 2*
src[5*step])>>4;
475 src[6*step]= (sums[6] + sums[8] + 2*
src[6*step])>>4;
476 src[7*step]= (sums[7] + sums[9] + 2*
src[7*step])>>4;
479 const int middleEnergy= 5*(
src[4*step] -
src[3*step]) + 2*(
src[2*step] -
src[5*step]);
481 if(
FFABS(middleEnergy) < 8*
QP){
482 const int q=(
src[3*step] -
src[4*step])/2;
483 const int leftEnergy= 5*(
src[2*step] -
src[1*step]) + 2*(
src[0*step] -
src[3*step]);
484 const int rightEnergy= 5*(
src[6*step] -
src[5*step]) + 2*(
src[4*step] -
src[7*step]);
490 d*=
FFSIGN(-middleEnergy);
501 d= (d < 0) ? 32 : -32;
519 #define TEMPLATE_PP_C 1
523 # define TEMPLATE_PP_ALTIVEC 1
528 #if ARCH_X86 && HAVE_INLINE_ASM
529 # if CONFIG_RUNTIME_CPUDETECT
530 # define TEMPLATE_PP_MMX 1
532 # define TEMPLATE_PP_MMXEXT 1
534 # define TEMPLATE_PP_3DNOW 1
536 # define TEMPLATE_PP_SSE2 1
539 # if HAVE_SSE2_INLINE
540 # define TEMPLATE_PP_SSE2 1
542 # elif HAVE_MMXEXT_INLINE
543 # define TEMPLATE_PP_MMXEXT 1
545 # elif HAVE_AMD3DNOW_INLINE
546 # define TEMPLATE_PP_3DNOW 1
548 # elif HAVE_MMX_INLINE
549 # define TEMPLATE_PP_MMX 1
556 const int8_t QPs[],
int QPStride,
int isColor,
PPContext *
c2);
561 pp_fn pp = postProcess_C;
567 #if CONFIG_RUNTIME_CPUDETECT
568 #if ARCH_X86 && HAVE_INLINE_ASM
579 pp = postProcess_SSE2;
580 #elif HAVE_MMXEXT_INLINE
581 pp = postProcess_MMX2;
582 #elif HAVE_AMD3DNOW_INLINE
583 pp = postProcess_3DNow;
584 #elif HAVE_MMX_INLINE
585 pp = postProcess_MMX;
587 pp = postProcess_altivec;
592 pp(
src, srcStride, dst, dstStride,
width,
height, QPs, QPStride, isColor,
c);
598 "Available postprocessing filters:\n"
600 "short long name short long option Description\n"
601 "* * a autoq CPU power dependent enabler\n"
602 " c chrom chrominance filtering enabled\n"
603 " y nochrom chrominance filtering disabled\n"
604 " n noluma luma filtering disabled\n"
605 "hb hdeblock (2 threshold) horizontal deblocking filter\n"
606 " 1. difference factor: default=32, higher -> more deblocking\n"
607 " 2. flatness threshold: default=39, lower -> more deblocking\n"
608 " the h & v deblocking filters share these\n"
609 " so you can't set different thresholds for h / v\n"
610 "vb vdeblock (2 threshold) vertical deblocking filter\n"
611 "ha hadeblock (2 threshold) horizontal deblocking filter\n"
612 "va vadeblock (2 threshold) vertical deblocking filter\n"
613 "h1 x1hdeblock experimental h deblock filter 1\n"
614 "v1 x1vdeblock experimental v deblock filter 1\n"
615 "dr dering deringing filter\n"
616 "al autolevels automatic brightness / contrast\n"
617 " f fullyrange stretch luminance to (0..255)\n"
618 "lb linblenddeint linear blend deinterlacer\n"
619 "li linipoldeint linear interpolating deinterlace\n"
620 "ci cubicipoldeint cubic interpolating deinterlacer\n"
621 "md mediandeint median deinterlacer\n"
622 "fd ffmpegdeint ffmpeg deinterlacer\n"
623 "l5 lowpass5 FIR lowpass deinterlacer\n"
624 "de default hb:a,vb:a,dr:a\n"
625 "fa fast h1:a,v1:a,dr:a\n"
626 "ac ha:a:128:7,va:a,dr:a\n"
627 "tn tmpnoise (3 threshold) temporal noise reducer\n"
628 " 1. <= 2. <= 3. larger -> stronger filtering\n"
629 "fq forceQuant <quantizer> force quantizer\n"
631 "<filterName>[:<option>[:<option>...]][[,|/][-]<filterName>[:<option>...]]...\n"
632 "long form example:\n"
633 "vdeblock:autoq/hdeblock:autoq/linblenddeint default,-vdeblock\n"
634 "short form example:\n"
635 "vb:a/hb:a/lb de,-vb\n"
645 static const char filterDelimiters[] =
",/";
646 static const char optionDelimiters[] =
":|";
655 if (!strcmp(
name,
"help")) {
657 for (p =
pp_help; strchr(p,
'\n'); p = strchr(p,
'\n') + 1) {
686 const char *filterName;
694 int numOfUnknownOptions=0;
698 filterToken=
av_strtok(p, filterDelimiters, &tokstate);
699 if(!filterToken)
break;
700 p+= strlen(filterToken) + 1;
701 filterName=
av_strtok(filterToken, optionDelimiters, &tokstate);
708 if(*filterName ==
'-'){
718 if(!strcmp(
"autoq",
option) || !strcmp(
"a",
option)) q= quality;
719 else if(!strcmp(
"nochrom",
option) || !strcmp(
"y",
option)) chrom=0;
720 else if(!strcmp(
"chrom",
option) || !strcmp(
"c",
option)) chrom=1;
721 else if(!strcmp(
"noluma",
option) || !strcmp(
"n",
option)) luma=0;
724 numOfUnknownOptions++;
740 spaceLeft= p -
temp + plen;
745 memmove(p + newlen, p, plen+1);
752 if( !strcmp(
filters[
i].longName, filterName)
753 || !strcmp(
filters[
i].shortName, filterName)){
760 if(q >=
filters[
i].minLumQuality && luma)
762 if(chrom==1 || (chrom==-1 &&
filters[
i].chromDefault))
771 if( !strcmp(
options[o],
"fullyrange")
775 numOfUnknownOptions--;
790 numOfUnknownOptions--;
791 if(numOfNoises >= 3)
break;
799 for(o=0;
options[o] && o<2; o++){
804 numOfUnknownOptions--;
813 for(o=0;
options[o] && o<1; o++){
818 numOfUnknownOptions--;
824 if(!filterNameOk) ppMode->
error++;
825 ppMode->
error += numOfUnknownOptions;
847 int mbWidth = (
width+15)>>4;
848 int mbHeight= (
height+15)>>4;
852 c->qpStride= qpStride;
868 reallocAlign((
void **)&
c->nonBQPTable, qpStride*mbHeight*
sizeof(int8_t));
869 reallocAlign((
void **)&
c->stdQPTable, qpStride*mbHeight*
sizeof(int8_t));
870 reallocAlign((
void **)&
c->forcedQPTable, mbWidth*
sizeof(int8_t));
882 int qpStride= (
width+15)/16 + 2;
889 c->hChromaSubSample= cpuCaps&0x3;
890 c->vChromaSubSample= (cpuCaps>>4)&0x3;
892 c->hChromaSubSample= 1;
893 c->vChromaSubSample= 1;
936 uint8_t * dst[3],
const int dstStride[3],
938 const int8_t *QP_store,
int QPStride,
939 pp_mode *vm,
void *vc,
int pict_type)
941 int mbWidth = (
width+15)>>4;
942 int mbHeight= (
height+15)>>4;
946 int absQPStride =
FFABS(QPStride);
949 if(
c->stride < minStride ||
c->qpStride < absQPStride)
951 FFMAX(minStride,
c->stride),
952 FFMAX(
c->qpStride, absQPStride));
956 QP_store=
c->forcedQPTable;
957 absQPStride = QPStride = 0;
959 for(
i=0;
i<mbWidth;
i++)
c->forcedQPTable[
i]=
mode->forcedQuant;
961 for(
i=0;
i<mbWidth;
i++)
c->forcedQPTable[
i]= 1;
966 const int count=
FFMAX(mbHeight * absQPStride, mbWidth);
967 for(
i=0;
i<(count>>2);
i++){
970 for(
i<<=2;
i<count;
i++){
971 c->stdQPTable[
i] = QP_store[
i]>>1;
973 QP_store=
c->stdQPTable;
974 QPStride= absQPStride;
979 for(y=0; y<mbHeight; y++){
980 for(x=0; x<mbWidth; x++){
988 if((pict_type&7)!=3){
991 const int count=
FFMAX(mbHeight * QPStride, mbWidth);
992 for(
i=0;
i<(count>>2);
i++){
995 for(
i<<=2;
i<count;
i++){
996 c->nonBQPTable[
i] = QP_store[
i] & 0x3F;
1000 for(
i=0;
i<mbHeight;
i++) {
1001 for(j=0; j<absQPStride; j++) {
1002 c->nonBQPTable[
i*absQPStride+j] = QP_store[
i*QPStride+j] & 0x3F;
1014 if (!(
src[1] &&
src[2] && dst[1] && dst[2]))
1020 if(
mode->chromMode){
1026 else if(srcStride[1] == dstStride[1] && srcStride[2] == dstStride[2]){
1032 memcpy(&(dst[1][y*dstStride[1]]), &(
src[1][y*srcStride[1]]),
width);
1033 memcpy(&(dst[2][y*dstStride[2]]), &(
src[2][y*srcStride[2]]),
width);
static double val(void *priv, double ch)
simple assert() macros that are a bit more flexible than ISO C assert().
#define av_assert0(cond)
assert() equivalent, that is always enabled.
Convenience header that includes libavutil's core.
#define FFABS(a)
Absolute value, Note, INT_MIN / INT64_MIN result in undefined behavior as they are not representable ...
#define FFMPEG_CONFIGURATION
int av_get_cpu_flags(void)
Return the flags which specify extensions supported by the CPU.
#define AV_CPU_FLAG_SSE2
PIV SSE2 functions.
#define AV_CPU_FLAG_MMXEXT
SSE integer functions or AMD MMX ext.
#define AV_CPU_FLAG_ALTIVEC
standard
#define AV_CPU_FLAG_MMX
standard MMX
#define AV_CPU_FLAG_3DNOW
AMD 3DNOW.
mode
Use these values in ebur128_init (or'ed).
const OptionDef options[]
#define AV_LOG_DEBUG
Stuff which is only useful for libav* developers.
#define AV_LOG_INFO
Standard information.
#define AV_LOG_ERROR
Something went wrong and cannot losslessly be recovered.
void * av_mallocz(size_t size)
Allocate a memory block with alignment suitable for all memory accesses (including vectors if availab...
char * av_strtok(char *s, const char *delim, char **saveptr)
Split the string into several tokens which can be accessed by successive calls to av_strtok().
size_t av_strlcpy(char *dst, const char *src, size_t size)
Copy the string src to dst, but no more than size - 1 bytes, and null-terminate dst.
const char pp_help[]
a simple help text
#define PP_CPU_CAPS_3DNOW
const char * postproc_configuration(void)
Return the libpostproc build-time configuration.
const char * postproc_license(void)
Return the libpostproc license.
av_cold pp_context * pp_get_context(int width, int height, int cpuCaps)
#define PP_CPU_CAPS_ALTIVEC
unsigned postproc_version(void)
Return the LIBPOSTPROC_VERSION_INT constant.
av_cold void pp_free_context(void *vc)
#define PP_PICT_TYPE_QP2
MPEG2 style QScale.
pp_mode * pp_get_mode_by_name_and_quality(const char *name, int quality)
Return a pp_mode or NULL if an error occurred.
void pp_postprocess(const uint8_t *src[3], const int srcStride[3], uint8_t *dst[3], const int dstStride[3], int width, int height, const int8_t *QP_store, int QPStride, pp_mode *vm, void *vc, int pict_type)
void pp_free_mode(pp_mode *mode)
#define LIBPOSTPROC_VERSION_INT
#define LIBPOSTPROC_VERSION_MICRO
static const uint16_t mask[17]
static const char * context_to_name(void *ptr)
void(* pp_fn)(const uint8_t src[], int srcStride, uint8_t dst[], int dstStride, int width, int height, const int8_t QPs[], int QPStride, int isColor, PPContext *c2)
static const char *const replaceTable[]
static const struct PPFilter filters[]
static av_always_inline void do_a_deblock_C(uint8_t *src, int step, int stride, const PPContext *c, int mode)
accurate deblock filter
static int isHorizMinMaxOk_C(const uint8_t src[], int stride, int QP)
static void postProcess(const uint8_t src[], int srcStride, uint8_t dst[], int dstStride, int width, int height, const int8_t QPs[], int QPStride, int isColor, pp_mode *vm, pp_context *vc)
static int horizClassify_C(const uint8_t src[], int stride, const PPContext *c)
DECLARE_ASM_CONST(8, int, deringThreshold)
static int vertClassify_C(const uint8_t src[], int stride, const PPContext *c)
const char postproc_ffversion[]
#define OPTIONS_ARRAY_SIZE
static const AVClass av_codec_context_class
static void horizX1Filter(uint8_t *src, int stride, int QP)
Experimental Filter 1 (Horizontal) will not damage linear gradients Flat blocks should look like they...
static void reallocAlign(void **p, int size)
static int isHorizDC_C(const uint8_t src[], int stride, const PPContext *c)
Check if the given 8x8 Block is mostly "flat".
static void reallocBuffers(PPContext *c, int width, int height, int stride, int qpStride)
static void doHorizDefFilter_C(uint8_t dst[], int stride, const PPContext *c)
static void doHorizLowPass_C(uint8_t dst[], int stride, const PPContext *c)
Do a horizontal low pass filter on the 10x8 block (dst points to middle 8x8 Block) using the 9-Tap Fi...
static int isVertDC_C(const uint8_t src[], int stride, const PPContext *c)
Check if the middle 8x8 Block in the given 8x16 block is flat.
static int isVertMinMaxOk_C(const uint8_t src[], int stride, int QP)
#define GET_MODE_BUFFER_SIZE
#define FFMPEG_DEINT_FILTER
#define LINEAR_IPOL_DEINT_FILTER
#define CUBIC_IPOL_DEINT_FILTER
#define TEMP_NOISE_FILTER
static void linecpy(void *dest, const void *src, int lines, int stride)
#define MEDIAN_DEINT_FILTER
#define LEVEL_FIX
Brightness & Contrast.
#define LINEAR_BLEND_DEINT_FILTER
#define LOWPASS5_DEINT_FILTER
mmx/mmx2/3dnow postprocess code.
typedef void(RENAME(mix_any_func_type))
#define FF_ARRAY_ELEMS(a)
Describe the class of an AVClass context structure.
Rational number (pair of numerator and denominator).
int mask
Bitmask to turn this filter on.
int maxTmpNoise[3]
for Temporal Noise Reducing filter (Maximal sum of abs differences)
int chromMode
activates filters for chrominance
int maxAllowedY
for brightness correction
AVRational maxClippedThreshold
amount of "black" you are willing to lose to get a brightness-corrected picture
int error
non zero on error
int forcedQuant
quantizer if FORCE_QUANT is used
int lumMode
activates filters for luminance
int minAllowedY
for brightness correction
Contains misc utility macros and inline functions.