Statistics
| Branch: | Revision:

ffmpeg / libavfilter / libmpcodecs / vf_decimate.c @ e4852fb3

History | View | Annotate | Download (5.68 KB)

1
/*
2
 * This file is part of MPlayer.
3
 *
4
 * MPlayer is free software; you can redistribute it and/or modify
5
 * it under the terms of the GNU General Public License as published by
6
 * the Free Software Foundation; either version 2 of the License, or
7
 * (at your option) any later version.
8
 *
9
 * MPlayer is distributed in the hope that it will be useful,
10
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
11
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
12
 * GNU General Public License for more details.
13
 *
14
 * You should have received a copy of the GNU General Public License along
15
 * with MPlayer; if not, write to the Free Software Foundation, Inc.,
16
 * 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
17
 */
18

    
19
#include <stdio.h>
20
#include <stdlib.h>
21
#include <string.h>
22

    
23
#include "config.h"
24
#include "mp_msg.h"
25
#include "cpudetect.h"
26

    
27
#include "img_format.h"
28
#include "mp_image.h"
29
#include "vf.h"
30

    
31
#include "libvo/fastmemcpy.h"
32

    
33

    
34
struct vf_priv_s {
35
    int hi, lo;
36
    float frac;
37
    int max, last, cnt;
38
};
39

    
40
#if HAVE_MMX && HAVE_EBX_AVAILABLE
41
static int diff_MMX(unsigned char *old, unsigned char *new, int os, int ns)
42
{
43
    volatile short out[4];
44
    __asm__ (
45
        "movl $8, %%ecx \n\t"
46
        "pxor %%mm4, %%mm4 \n\t"
47
        "pxor %%mm7, %%mm7 \n\t"
48

    
49
        ASMALIGN(4)
50
        "1: \n\t"
51

    
52
        "movq (%%"REG_S"), %%mm0 \n\t"
53
        "movq (%%"REG_S"), %%mm2 \n\t"
54
        "add %%"REG_a", %%"REG_S" \n\t"
55
        "movq (%%"REG_D"), %%mm1 \n\t"
56
        "add %%"REG_b", %%"REG_D" \n\t"
57
        "psubusb %%mm1, %%mm2 \n\t"
58
        "psubusb %%mm0, %%mm1 \n\t"
59
        "movq %%mm2, %%mm0 \n\t"
60
        "movq %%mm1, %%mm3 \n\t"
61
        "punpcklbw %%mm7, %%mm0 \n\t"
62
        "punpcklbw %%mm7, %%mm1 \n\t"
63
        "punpckhbw %%mm7, %%mm2 \n\t"
64
        "punpckhbw %%mm7, %%mm3 \n\t"
65
        "paddw %%mm0, %%mm4 \n\t"
66
        "paddw %%mm1, %%mm4 \n\t"
67
        "paddw %%mm2, %%mm4 \n\t"
68
        "paddw %%mm3, %%mm4 \n\t"
69

    
70
        "decl %%ecx \n\t"
71
        "jnz 1b \n\t"
72
        "movq %%mm4, (%%"REG_d") \n\t"
73
        "emms \n\t"
74
        :
75
        : "S" (old), "D" (new), "a" ((long)os), "b" ((long)ns), "d" (out)
76
        : "%ecx", "memory"
77
        );
78
    return out[0]+out[1]+out[2]+out[3];
79
}
80
#endif
81

    
82
static int diff_C(unsigned char *old, unsigned char *new, int os, int ns)
83
{
84
    int x, y, d=0;
85
    for (y = 8; y; y--) {
86
        for (x = 8; x; x--) {
87
            d += abs(new[x] - old[x]);
88
        }
89
        new += ns;
90
        old += os;
91
    }
92
    return d;
93
}
94

    
95
static int (*diff)(unsigned char *, unsigned char *, int, int);
96

    
97
static int diff_to_drop_plane(int hi, int lo, float frac, unsigned char *old, unsigned char *new, int w, int h, int os, int ns)
98
{
99
    int x, y;
100
    int d, c=0;
101
    int t = (w/16)*(h/16)*frac;
102
    for (y = 0; y < h-7; y += 4) {
103
        for (x = 8; x < w-7; x += 4) {
104
            d = diff(old+x+y*os, new+x+y*ns, os, ns);
105
            if (d > hi) return 0;
106
            if (d > lo) {
107
                c++;
108
                if (c > t) return 0;
109
            }
110
        }
111
    }
112
    return 1;
113
}
114

    
115
static int diff_to_drop(int hi, int lo, float frac, mp_image_t *old, mp_image_t *new)
116
{
117
    if (new->flags & MP_IMGFLAG_PLANAR) {
118
        return diff_to_drop_plane(hi,lo,frac, old->planes[0], new->planes[0],
119
            new->w, new->h, old->stride[0], new->stride[0])
120
            && diff_to_drop_plane(hi,lo,frac, old->planes[1], new->planes[1],
121
            new->chroma_width, new->chroma_height,
122
            old->stride[1], new->stride[1])
123
            && diff_to_drop_plane(hi,lo,frac, old->planes[2], new->planes[2],
124
            new->chroma_width, new->chroma_height,
125
            old->stride[2], new->stride[2]);
126
    }
127
    return diff_to_drop_plane(hi,lo,frac, old->planes[0], new->planes[0],
128
        new->w*(new->bpp/8), new->h, old->stride[0], new->stride[0]);
129
}
130

    
131
static int put_image(struct vf_instance *vf, mp_image_t *mpi, double pts)
132
{
133
    mp_image_t *dmpi;
134

    
135
    dmpi = vf_get_image(vf->next, mpi->imgfmt,
136
        MP_IMGTYPE_STATIC, MP_IMGFLAG_ACCEPT_STRIDE |
137
        MP_IMGFLAG_PRESERVE | MP_IMGFLAG_READABLE,
138
        mpi->width, mpi->height);
139
    dmpi->qscale = mpi->qscale;
140
    dmpi->qstride = mpi->qstride;
141
    dmpi->qscale_type = mpi->qscale_type;
142

    
143
    if (diff_to_drop(vf->priv->hi, vf->priv->lo, vf->priv->frac, dmpi, mpi)) {
144
        if (vf->priv->max == 0)
145
            return 0;
146
        else if ((vf->priv->max > 0) && (vf->priv->cnt++ < vf->priv->max))
147
            return 0;
148
        else if ((vf->priv->max < 0) && (vf->priv->last+1 >= -vf->priv->max))
149
            return vf->priv->last=0;
150
    }
151
    vf->priv->last++;
152
    vf->priv->cnt=0;
153

    
154
    memcpy_pic(dmpi->planes[0], mpi->planes[0], mpi->w, mpi->h,
155
        dmpi->stride[0], mpi->stride[0]);
156
    if (mpi->flags & MP_IMGFLAG_PLANAR) {
157
        memcpy_pic(dmpi->planes[1], mpi->planes[1],
158
            mpi->chroma_width, mpi->chroma_height,
159
            dmpi->stride[1], mpi->stride[1]);
160
        memcpy_pic(dmpi->planes[2], mpi->planes[2],
161
            mpi->chroma_width, mpi->chroma_height,
162
            dmpi->stride[2], mpi->stride[2]);
163
    }
164
    return vf_next_put_image(vf, dmpi, pts);
165
}
166

    
167
static void uninit(struct vf_instance *vf)
168
{
169
    free(vf->priv);
170
}
171

    
172
static int vf_open(vf_instance_t *vf, char *args)
173
{
174
    struct vf_priv_s *p;
175
    vf->put_image = put_image;
176
    vf->uninit = uninit;
177
    vf->default_reqs = VFCAP_ACCEPT_STRIDE;
178
    vf->priv = p = calloc(1, sizeof(struct vf_priv_s));
179
    p->max = 0;
180
    p->hi = 64*12;
181
    p->lo = 64*5;
182
    p->frac = 0.33;
183
    if (args) sscanf(args, "%d:%d:%d:%f", &p->max, &p->hi, &p->lo, &p->frac);
184
    diff = diff_C;
185
#if HAVE_MMX && HAVE_EBX_AVAILABLE
186
    if(gCpuCaps.hasMMX) diff = diff_MMX;
187
#endif
188
    return 1;
189
}
190

    
191
const vf_info_t vf_info_decimate = {
192
    "near-duplicate frame remover",
193
    "decimate",
194
    "Rich Felker",
195
    "",
196
    vf_open,
197
    NULL
198
};