libmpcodecs/vf_decimate.c

   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 };