git.sesse.net Git - ffmpeg/blob - libavfilter/vf_stereo3d.c

   1 /*
   2  * Copyright (c) 2010 Gordon Schmidt <gordon.schmidt <at> s2000.tu-chemnitz.de>
   3  * Copyright (c) 2013 Paul B Mahol
   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 General Public
   9  * License as published by the Free Software Foundation; either
  10  * version 2 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
  15  * GNU General Public License for more details.
  16  *
  17  * You should have received a copy of the GNU General Public License along
  18  * 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 "libavutil/avassert.h"
  23 #include "libavutil/imgutils.h"
  24 #include "libavutil/opt.h"
  25 #include "libavutil/parseutils.h"
  26 #include "libavutil/pixdesc.h"
  27 #include "avfilter.h"
  28 #include "drawutils.h"
  29 #include "formats.h"
  30 #include "internal.h"
  31 #include "video.h"
  32
  33 enum StereoCode {
  34     ANAGLYPH_RC_GRAY,   // anaglyph red/cyan gray
  35     ANAGLYPH_RC_HALF,   // anaglyph red/cyan half colored
  36     ANAGLYPH_RC_COLOR,  // anaglyph red/cyan colored
  37     ANAGLYPH_RC_DUBOIS, // anaglyph red/cyan dubois
  38     ANAGLYPH_GM_GRAY,   // anaglyph green/magenta gray
  39     ANAGLYPH_GM_HALF,   // anaglyph green/magenta half colored
  40     ANAGLYPH_GM_COLOR,  // anaglyph green/magenta colored
  41     ANAGLYPH_GM_DUBOIS, // anaglyph green/magenta dubois
  42     ANAGLYPH_YB_GRAY,   // anaglyph yellow/blue gray
  43     ANAGLYPH_YB_HALF,   // anaglyph yellow/blue half colored
  44     ANAGLYPH_YB_COLOR,  // anaglyph yellow/blue colored
  45     ANAGLYPH_YB_DUBOIS, // anaglyph yellow/blue dubois
  46     ANAGLYPH_RB_GRAY,   // anaglyph red/blue gray
  47     ANAGLYPH_RG_GRAY,   // anaglyph red/green gray
  48     MONO_L,             // mono output for debugging (left eye only)
  49     MONO_R,             // mono output for debugging (right eye only)
  50     INTERLEAVE_ROWS_LR, // row-interleave (left eye has top row)
  51     INTERLEAVE_ROWS_RL, // row-interleave (right eye has top row)
  52     SIDE_BY_SIDE_LR,    // side by side parallel (left eye left, right eye right)
  53     SIDE_BY_SIDE_RL,    // side by side crosseye (right eye left, left eye right)
  54     SIDE_BY_SIDE_2_LR,  // side by side parallel with half width resolution
  55     SIDE_BY_SIDE_2_RL,  // side by side crosseye with half width resolution
  56     ABOVE_BELOW_LR,     // above-below (left eye above, right eye below)
  57     ABOVE_BELOW_RL,     // above-below (right eye above, left eye below)
  58     ABOVE_BELOW_2_LR,   // above-below with half height resolution
  59     ABOVE_BELOW_2_RL,   // above-below with half height resolution
  60     ALTERNATING_LR,     // alternating frames (left eye first, right eye second)
  61     ALTERNATING_RL,     // alternating frames (right eye first, left eye second)
  62     STEREO_CODE_COUNT   // TODO: needs autodetection
  63 };
  64
  65 typedef struct StereoComponent {
  66     int format;                 ///< StereoCode
  67     int width, height;
  68     int off_left, off_right;
  69     int off_lstep, off_rstep;
  70     int row_left, row_right;
  71 } StereoComponent;
  72
  73 static const int ana_coeff[][3][6] = {
  74   [ANAGLYPH_RB_GRAY]   =
  75     {{19595, 38470,  7471,     0,     0,     0},
  76      {    0,     0,     0,     0,     0,     0},
  77      {    0,     0,     0, 19595, 38470,  7471}},
  78   [ANAGLYPH_RG_GRAY]   =
  79     {{19595, 38470,  7471,     0,     0,     0},
  80      {    0,     0,     0, 19595, 38470,  7471},
  81      {    0,     0,     0,     0,     0,     0}},
  82   [ANAGLYPH_RC_GRAY]   =
  83     {{19595, 38470,  7471,     0,     0,     0},
  84      {    0,     0,     0, 19595, 38470,  7471},
  85      {    0,     0,     0, 19595, 38470,  7471}},
  86   [ANAGLYPH_RC_HALF]   =
  87     {{19595, 38470,  7471,     0,     0,     0},
  88      {    0,     0,     0,     0, 65536,     0},
  89      {    0,     0,     0,     0,     0, 65536}},
  90   [ANAGLYPH_RC_COLOR]  =
  91     {{65536,     0,     0,     0,     0,     0},
  92      {    0,     0,     0,     0, 65536,     0},
  93      {    0,     0,     0,     0,     0, 65536}},
  94   [ANAGLYPH_RC_DUBOIS] =
  95     {{29891, 32800, 11559, -2849, -5763,  -102},
  96      {-2627, -2479, -1033, 24804, 48080, -1209},
  97      { -997, -1350,  -358, -4729, -7403, 80373}},
  98   [ANAGLYPH_GM_GRAY]   =
  99     {{    0,     0,     0, 19595, 38470,  7471},
 100      {19595, 38470,  7471,     0,     0,     0},
 101      {    0,     0,     0, 19595, 38470,  7471}},
 102   [ANAGLYPH_GM_HALF]   =
 103     {{    0,     0,     0, 65536,     0,     0},
 104      {19595, 38470,  7471,     0,     0,     0},
 105      {    0,     0,     0,     0,     0, 65536}},
 106   [ANAGLYPH_GM_COLOR]  =
 107     {{    0,     0,     0, 65536,     0,     0},
 108      {    0, 65536,     0,     0,     0,     0},
 109      {    0,     0,     0,     0,     0, 65536}},
 110   [ANAGLYPH_GM_DUBOIS]  =
 111     {{-4063,-10354, -2556, 34669, 46203,  1573},
 112      {18612, 43778,  9372, -1049,  -983, -4260},
 113      { -983, -1769,  1376,   590,  4915, 61407}},
 114   [ANAGLYPH_YB_GRAY]   =
 115     {{    0,     0,     0, 19595, 38470,  7471},
 116      {    0,     0,     0, 19595, 38470,  7471},
 117      {19595, 38470,  7471,     0,     0,     0}},
 118   [ANAGLYPH_YB_HALF]   =
 119     {{    0,     0,     0, 65536,     0,     0},
 120      {    0,     0,     0,     0, 65536,     0},
 121      {19595, 38470,  7471,     0,     0,     0}},
 122   [ANAGLYPH_YB_COLOR]  =
 123     {{    0,     0,     0, 65536,     0,     0},
 124      {    0,     0,     0,     0, 65536,     0},
 125      {    0,     0, 65536,     0,     0,     0}},
 126   [ANAGLYPH_YB_DUBOIS] =
 127     {{65535,-12650,18451,   -987, -7590, -1049},
 128      {-1604, 56032, 4196,    370,  3826, -1049},
 129      {-2345,-10676, 1358,   5801, 11416, 56217}},
 130 };
 131
 132 typedef struct Stereo3DContext {
 133     const AVClass *class;
 134     StereoComponent in, out;
 135     int width, height;
 136     int row_step;
 137     const int *ana_matrix[3];
 138     int nb_planes;
 139     int linesize[4];
 140     int pheight[4];
 141     int hsub, vsub;
 142     int pixstep[4];
 143     AVFrame *prev;
 144     double ts_unit;
 145 } Stereo3DContext;
 146
 147 #define OFFSET(x) offsetof(Stereo3DContext, x)
 148 #define FLAGS AV_OPT_FLAG_FILTERING_PARAM|AV_OPT_FLAG_VIDEO_PARAM
 149
 150 static const AVOption stereo3d_options[] = {
 151     { "in",    "set input format",  OFFSET(in.format),  AV_OPT_TYPE_INT, {.i64=SIDE_BY_SIDE_LR}, SIDE_BY_SIDE_LR, STEREO_CODE_COUNT-1, FLAGS, "in"},
 152     { "ab2l",  "above below half height left first",  0, AV_OPT_TYPE_CONST, {.i64=ABOVE_BELOW_2_LR},  0, 0, FLAGS, "in" },
 153     { "ab2r",  "above below half height right first", 0, AV_OPT_TYPE_CONST, {.i64=ABOVE_BELOW_2_RL},  0, 0, FLAGS, "in" },
 154     { "abl",   "above below left first",              0, AV_OPT_TYPE_CONST, {.i64=ABOVE_BELOW_LR},    0, 0, FLAGS, "in" },
 155     { "abr",   "above below right first",             0, AV_OPT_TYPE_CONST, {.i64=ABOVE_BELOW_RL},    0, 0, FLAGS, "in" },
 156     { "al",    "alternating frames left first",       0, AV_OPT_TYPE_CONST, {.i64=ALTERNATING_LR},    0, 0, FLAGS, "in" },
 157     { "ar",    "alternating frames right first",      0, AV_OPT_TYPE_CONST, {.i64=ALTERNATING_RL},    0, 0, FLAGS, "in" },
 158     { "sbs2l", "side by side half width left first",  0, AV_OPT_TYPE_CONST, {.i64=SIDE_BY_SIDE_2_LR}, 0, 0, FLAGS, "in" },
 159     { "sbs2r", "side by side half width right first", 0, AV_OPT_TYPE_CONST, {.i64=SIDE_BY_SIDE_2_RL}, 0, 0, FLAGS, "in" },
 160     { "sbsl",  "side by side left first",             0, AV_OPT_TYPE_CONST, {.i64=SIDE_BY_SIDE_LR},   0, 0, FLAGS, "in" },
 161     { "sbsr",  "side by side right first",            0, AV_OPT_TYPE_CONST, {.i64=SIDE_BY_SIDE_RL},   0, 0, FLAGS, "in" },
 162     { "out",   "set output format", OFFSET(out.format), AV_OPT_TYPE_INT, {.i64=ANAGLYPH_RC_DUBOIS}, 0, STEREO_CODE_COUNT-1, FLAGS, "out"},
 163     { "ab2l",  "above below half height left first",  0, AV_OPT_TYPE_CONST, {.i64=ABOVE_BELOW_2_LR},   0, 0, FLAGS, "out" },
 164     { "ab2r",  "above below half height right first", 0, AV_OPT_TYPE_CONST, {.i64=ABOVE_BELOW_2_RL},   0, 0, FLAGS, "out" },
 165     { "abl",   "above below left first",              0, AV_OPT_TYPE_CONST, {.i64=ABOVE_BELOW_LR},     0, 0, FLAGS, "out" },
 166     { "abr",   "above below right first",             0, AV_OPT_TYPE_CONST, {.i64=ABOVE_BELOW_RL},     0, 0, FLAGS, "out" },
 167     { "agmc",  "anaglyph green magenta color",        0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_GM_COLOR},  0, 0, FLAGS, "out" },
 168     { "agmd",  "anaglyph green magenta dubois",       0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_GM_DUBOIS}, 0, 0, FLAGS, "out" },
 169     { "agmg",  "anaglyph green magenta gray",         0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_GM_GRAY},   0, 0, FLAGS, "out" },
 170     { "agmh",  "anaglyph green magenta half color",   0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_GM_HALF},   0, 0, FLAGS, "out" },
 171     { "al",    "alternating frames left first",       0, AV_OPT_TYPE_CONST, {.i64=ALTERNATING_LR},     0, 0, FLAGS, "out" },
 172     { "ar",    "alternating frames right first",      0, AV_OPT_TYPE_CONST, {.i64=ALTERNATING_RL},     0, 0, FLAGS, "out" },
 173     { "arbg",  "anaglyph red blue gray",              0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_RB_GRAY},   0, 0, FLAGS, "out" },
 174     { "arcc",  "anaglyph red cyan color",             0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_RC_COLOR},  0, 0, FLAGS, "out" },
 175     { "arcd",  "anaglyph red cyan dubois",            0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_RC_DUBOIS}, 0, 0, FLAGS, "out" },
 176     { "arcg",  "anaglyph red cyan gray",              0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_RC_GRAY},   0, 0, FLAGS, "out" },
 177     { "arch",  "anaglyph red cyan half color",        0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_RC_HALF},   0, 0, FLAGS, "out" },
 178     { "argg",  "anaglyph red green gray",             0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_RG_GRAY},   0, 0, FLAGS, "out" },
 179     { "aybc",  "anaglyph yellow blue color",          0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_YB_COLOR},  0, 0, FLAGS, "out" },
 180     { "aybd",  "anaglyph yellow blue dubois",         0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_YB_DUBOIS}, 0, 0, FLAGS, "out" },
 181     { "aybg",  "anaglyph yellow blue gray",           0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_YB_GRAY},   0, 0, FLAGS, "out" },
 182     { "aybh",  "anaglyph yellow blue half color",     0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_YB_HALF},   0, 0, FLAGS, "out" },
 183     { "irl",   "interleave rows left first",          0, AV_OPT_TYPE_CONST, {.i64=INTERLEAVE_ROWS_LR}, 0, 0, FLAGS, "out" },
 184     { "irr",   "interleave rows right first",         0, AV_OPT_TYPE_CONST, {.i64=INTERLEAVE_ROWS_RL}, 0, 0, FLAGS, "out" },
 185     { "ml",    "mono left",                           0, AV_OPT_TYPE_CONST, {.i64=MONO_L},             0, 0, FLAGS, "out" },
 186     { "mr",    "mono right",                          0, AV_OPT_TYPE_CONST, {.i64=MONO_R},             0, 0, FLAGS, "out" },
 187     { "sbs2l", "side by side half width left first",  0, AV_OPT_TYPE_CONST, {.i64=SIDE_BY_SIDE_2_LR},  0, 0, FLAGS, "out" },
 188     { "sbs2r", "side by side half width right first", 0, AV_OPT_TYPE_CONST, {.i64=SIDE_BY_SIDE_2_RL},  0, 0, FLAGS, "out" },
 189     { "sbsl",  "side by side left first",             0, AV_OPT_TYPE_CONST, {.i64=SIDE_BY_SIDE_LR},    0, 0, FLAGS, "out" },
 190     { "sbsr",  "side by side right first",            0, AV_OPT_TYPE_CONST, {.i64=SIDE_BY_SIDE_RL},    0, 0, FLAGS, "out" },
 191     { NULL }
 192 };
 193
 194 AVFILTER_DEFINE_CLASS(stereo3d);
 195
 196 static const enum AVPixelFormat anaglyph_pix_fmts[] = {
 197     AV_PIX_FMT_RGB24, AV_PIX_FMT_BGR24,
 198     AV_PIX_FMT_NONE
 199 };
 200
 201 static const enum AVPixelFormat other_pix_fmts[] = {
 202     AV_PIX_FMT_RGB24, AV_PIX_FMT_BGR24,
 203     AV_PIX_FMT_RGB48BE, AV_PIX_FMT_BGR48BE,
 204     AV_PIX_FMT_RGB48LE, AV_PIX_FMT_BGR48LE,
 205     AV_PIX_FMT_RGBA64BE, AV_PIX_FMT_BGRA64BE,
 206     AV_PIX_FMT_RGBA64LE, AV_PIX_FMT_BGRA64LE,
 207     AV_PIX_FMT_RGBA,  AV_PIX_FMT_BGRA,
 208     AV_PIX_FMT_ARGB,  AV_PIX_FMT_ABGR,
 209     AV_PIX_FMT_RGB0,  AV_PIX_FMT_BGR0,
 210     AV_PIX_FMT_0RGB,  AV_PIX_FMT_0BGR,
 211     AV_PIX_FMT_GBRP,
 212     AV_PIX_FMT_GBRP9BE,  AV_PIX_FMT_GBRP9LE,
 213     AV_PIX_FMT_GBRP10BE, AV_PIX_FMT_GBRP10LE,
 214     AV_PIX_FMT_GBRP12BE, AV_PIX_FMT_GBRP12LE,
 215     AV_PIX_FMT_GBRP14BE, AV_PIX_FMT_GBRP14LE,
 216     AV_PIX_FMT_GBRP16BE, AV_PIX_FMT_GBRP16LE,
 217     AV_PIX_FMT_YUV410P,
 218     AV_PIX_FMT_YUV411P,
 219     AV_PIX_FMT_YUV420P, AV_PIX_FMT_YUVA420P,
 220     AV_PIX_FMT_YUV422P, AV_PIX_FMT_YUVA422P,
 221     AV_PIX_FMT_YUV440P,
 222     AV_PIX_FMT_YUV444P, AV_PIX_FMT_YUVA444P,
 223     AV_PIX_FMT_YUVJ411P,
 224     AV_PIX_FMT_YUVJ420P,
 225     AV_PIX_FMT_YUVJ422P,
 226     AV_PIX_FMT_YUVJ440P,
 227     AV_PIX_FMT_YUVJ444P,
 228     AV_PIX_FMT_YUV420P9LE,  AV_PIX_FMT_YUVA420P9LE,
 229     AV_PIX_FMT_YUV420P9BE,  AV_PIX_FMT_YUVA420P9BE,
 230     AV_PIX_FMT_YUV422P9LE,  AV_PIX_FMT_YUVA422P9LE,
 231     AV_PIX_FMT_YUV422P9BE,  AV_PIX_FMT_YUVA422P9BE,
 232     AV_PIX_FMT_YUV444P9LE,  AV_PIX_FMT_YUVA444P9LE,
 233     AV_PIX_FMT_YUV444P9BE,  AV_PIX_FMT_YUVA444P9BE,
 234     AV_PIX_FMT_YUV420P10LE, AV_PIX_FMT_YUVA420P10LE,
 235     AV_PIX_FMT_YUV420P10BE, AV_PIX_FMT_YUVA420P10BE,
 236     AV_PIX_FMT_YUV422P10LE, AV_PIX_FMT_YUVA422P10LE,
 237     AV_PIX_FMT_YUV422P10BE, AV_PIX_FMT_YUVA422P10BE,
 238     AV_PIX_FMT_YUV444P10LE, AV_PIX_FMT_YUVA444P10LE,
 239     AV_PIX_FMT_YUV444P10BE, AV_PIX_FMT_YUVA444P10BE,
 240     AV_PIX_FMT_YUV420P12BE,  AV_PIX_FMT_YUV420P12LE,
 241     AV_PIX_FMT_YUV422P12BE,  AV_PIX_FMT_YUV422P12LE,
 242     AV_PIX_FMT_YUV444P12BE,  AV_PIX_FMT_YUV444P12LE,
 243     AV_PIX_FMT_YUV420P14BE,  AV_PIX_FMT_YUV420P14LE,
 244     AV_PIX_FMT_YUV422P14BE,  AV_PIX_FMT_YUV422P14LE,
 245     AV_PIX_FMT_YUV444P14BE,  AV_PIX_FMT_YUV444P14LE,
 246     AV_PIX_FMT_YUV420P16LE, AV_PIX_FMT_YUVA420P16LE,
 247     AV_PIX_FMT_YUV420P16BE, AV_PIX_FMT_YUVA420P16BE,
 248     AV_PIX_FMT_YUV422P16LE, AV_PIX_FMT_YUVA422P16LE,
 249     AV_PIX_FMT_YUV422P16BE, AV_PIX_FMT_YUVA422P16BE,
 250     AV_PIX_FMT_YUV444P16LE, AV_PIX_FMT_YUVA444P16LE,
 251     AV_PIX_FMT_YUV444P16BE, AV_PIX_FMT_YUVA444P16BE,
 252     AV_PIX_FMT_NONE
 253 };
 254
 255 static int query_formats(AVFilterContext *ctx)
 256 {
 257     Stereo3DContext *s = ctx->priv;
 258     const enum AVPixelFormat *pix_fmts;
 259     AVFilterFormats *fmts_list;
 260
 261     switch (s->out.format) {
 262     case ANAGLYPH_GM_COLOR:
 263     case ANAGLYPH_GM_DUBOIS:
 264     case ANAGLYPH_GM_GRAY:
 265     case ANAGLYPH_GM_HALF:
 266     case ANAGLYPH_RB_GRAY:
 267     case ANAGLYPH_RC_COLOR:
 268     case ANAGLYPH_RC_DUBOIS:
 269     case ANAGLYPH_RC_GRAY:
 270     case ANAGLYPH_RC_HALF:
 271     case ANAGLYPH_RG_GRAY:
 272     case ANAGLYPH_YB_COLOR:
 273     case ANAGLYPH_YB_DUBOIS:
 274     case ANAGLYPH_YB_GRAY:
 275     case ANAGLYPH_YB_HALF:
 276         pix_fmts = anaglyph_pix_fmts;
 277         break;
 278     default:
 279         pix_fmts = other_pix_fmts;
 280     }
 281
 282     fmts_list = ff_make_format_list(pix_fmts);
 283     if (!fmts_list)
 284         return AVERROR(ENOMEM);
 285     return ff_set_common_formats(ctx, fmts_list);
 286 }
 287
 288 static int config_output(AVFilterLink *outlink)
 289 {
 290     AVFilterContext *ctx = outlink->src;
 291     AVFilterLink *inlink = ctx->inputs[0];
 292     Stereo3DContext *s = ctx->priv;
 293     AVRational aspect = inlink->sample_aspect_ratio;
 294     AVRational fps = inlink->frame_rate;
 295     AVRational tb = inlink->time_base;
 296     const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(outlink->format);
 297     int ret;
 298
 299     switch (s->in.format) {
 300     case SIDE_BY_SIDE_2_LR:
 301     case SIDE_BY_SIDE_LR:
 302     case SIDE_BY_SIDE_2_RL:
 303     case SIDE_BY_SIDE_RL:
 304         if (inlink->w & 1) {
 305             av_log(ctx, AV_LOG_ERROR, "width must be even\n");
 306             return AVERROR_INVALIDDATA;
 307         }
 308         break;
 309     case ABOVE_BELOW_2_LR:
 310     case ABOVE_BELOW_LR:
 311     case ABOVE_BELOW_2_RL:
 312     case ABOVE_BELOW_RL:
 313         if (s->out.format == INTERLEAVE_ROWS_LR ||
 314             s->out.format == INTERLEAVE_ROWS_RL) {
 315             if (inlink->h & 3) {
 316                 av_log(ctx, AV_LOG_ERROR, "height must be multiple of 4\n");
 317                 return AVERROR_INVALIDDATA;
 318             }
 319         }
 320         if (inlink->h & 1) {
 321             av_log(ctx, AV_LOG_ERROR, "height must be even\n");
 322             return AVERROR_INVALIDDATA;
 323         }
 324         break;
 325     }
 326
 327     s->in.width     =
 328     s->width        = inlink->w;
 329     s->in.height    =
 330     s->height       = inlink->h;
 331     s->row_step     = 1;
 332     s->in.off_lstep =
 333     s->in.off_rstep =
 334     s->in.off_left  =
 335     s->in.off_right =
 336     s->in.row_left  =
 337     s->in.row_right = 0;
 338
 339     switch (s->in.format) {
 340     case SIDE_BY_SIDE_2_LR:
 341         aspect.num     *= 2;
 342     case SIDE_BY_SIDE_LR:
 343         s->width        = inlink->w / 2;
 344         s->in.off_right = s->width;
 345         break;
 346     case SIDE_BY_SIDE_2_RL:
 347         aspect.num     *= 2;
 348     case SIDE_BY_SIDE_RL:
 349         s->width        = inlink->w / 2;
 350         s->in.off_left  = s->width;
 351         break;
 352     case ABOVE_BELOW_2_LR:
 353         aspect.den     *= 2;
 354     case ABOVE_BELOW_LR:
 355         s->in.row_right =
 356         s->height       = inlink->h / 2;
 357         break;
 358     case ABOVE_BELOW_2_RL:
 359         aspect.den     *= 2;
 360     case ABOVE_BELOW_RL:
 361         s->in.row_left  =
 362         s->height       = inlink->h / 2;
 363         break;
 364     case ALTERNATING_RL:
 365     case ALTERNATING_LR:
 366         outlink->flags |= FF_LINK_FLAG_REQUEST_LOOP;
 367         fps.den        *= 2;
 368         tb.num         *= 2;
 369         break;
 370     default:
 371         av_log(ctx, AV_LOG_ERROR, "input format %d is not supported\n", s->in.format);
 372         return AVERROR(EINVAL);
 373     }
 374
 375     s->out.width     = s->width;
 376     s->out.height    = s->height;
 377     s->out.off_lstep =
 378     s->out.off_rstep =
 379     s->out.off_left  =
 380     s->out.off_right =
 381     s->out.row_left  =
 382     s->out.row_right = 0;
 383
 384     switch (s->out.format) {
 385     case ANAGLYPH_RB_GRAY:
 386     case ANAGLYPH_RG_GRAY:
 387     case ANAGLYPH_RC_GRAY:
 388     case ANAGLYPH_RC_HALF:
 389     case ANAGLYPH_RC_COLOR:
 390     case ANAGLYPH_RC_DUBOIS:
 391     case ANAGLYPH_GM_GRAY:
 392     case ANAGLYPH_GM_HALF:
 393     case ANAGLYPH_GM_COLOR:
 394     case ANAGLYPH_GM_DUBOIS:
 395     case ANAGLYPH_YB_GRAY:
 396     case ANAGLYPH_YB_HALF:
 397     case ANAGLYPH_YB_COLOR:
 398     case ANAGLYPH_YB_DUBOIS: {
 399         uint8_t rgba_map[4];
 400
 401         ff_fill_rgba_map(rgba_map, outlink->format);
 402         s->ana_matrix[rgba_map[0]] = &ana_coeff[s->out.format][0][0];
 403         s->ana_matrix[rgba_map[1]] = &ana_coeff[s->out.format][1][0];
 404         s->ana_matrix[rgba_map[2]] = &ana_coeff[s->out.format][2][0];
 405         break;
 406     }
 407     case SIDE_BY_SIDE_2_LR:
 408         aspect.den      *= 2;
 409     case SIDE_BY_SIDE_LR:
 410         s->out.width     = s->width * 2;
 411         s->out.off_right = s->width;
 412         break;
 413     case SIDE_BY_SIDE_2_RL:
 414         aspect.den      *= 2;
 415     case SIDE_BY_SIDE_RL:
 416         s->out.width     = s->width * 2;
 417         s->out.off_left  = s->width;
 418         break;
 419     case ABOVE_BELOW_2_LR:
 420         aspect.num      *= 2;
 421     case ABOVE_BELOW_LR:
 422         s->out.height    = s->height * 2;
 423         s->out.row_right = s->height;
 424         break;
 425     case ABOVE_BELOW_2_RL:
 426         aspect.num      *= 2;
 427     case ABOVE_BELOW_RL:
 428         s->out.height    = s->height * 2;
 429         s->out.row_left  = s->height;
 430         break;
 431     case INTERLEAVE_ROWS_LR:
 432         s->row_step      = 2;
 433         s->height        = s->height / 2;
 434         s->out.off_rstep =
 435         s->in.off_rstep  = 1;
 436         break;
 437     case INTERLEAVE_ROWS_RL:
 438         s->row_step      = 2;
 439         s->height        = s->height / 2;
 440         s->out.off_lstep =
 441         s->in.off_lstep  = 1;
 442         break;
 443     case MONO_R:
 444         s->in.off_left   = s->in.off_right;
 445         s->in.row_left   = s->in.row_right;
 446     case MONO_L:
 447         break;
 448     case ALTERNATING_RL:
 449     case ALTERNATING_LR:
 450         fps.num         *= 2;
 451         tb.den          *= 2;
 452         break;
 453     default:
 454         av_log(ctx, AV_LOG_ERROR, "output format %d is not supported\n", s->out.format);
 455         return AVERROR(EINVAL);
 456     }
 457
 458     outlink->w = s->out.width;
 459     outlink->h = s->out.height;
 460     outlink->frame_rate = fps;
 461     outlink->time_base = tb;
 462     outlink->sample_aspect_ratio = aspect;
 463
 464     if ((ret = av_image_fill_linesizes(s->linesize, outlink->format, s->width)) < 0)
 465         return ret;
 466     s->nb_planes = av_pix_fmt_count_planes(outlink->format);
 467     av_image_fill_max_pixsteps(s->pixstep, NULL, desc);
 468     s->ts_unit = av_q2d(av_inv_q(av_mul_q(outlink->frame_rate, outlink->time_base)));
 469     s->pheight[1] = s->pheight[2] = FF_CEIL_RSHIFT(s->height, desc->log2_chroma_h);
 470     s->pheight[0] = s->pheight[3] = s->height;
 471     s->hsub = desc->log2_chroma_w;
 472     s->vsub = desc->log2_chroma_h;
 473
 474     return 0;
 475 }
 476
 477 static inline uint8_t ana_convert(const int *coeff, const uint8_t *left, const uint8_t *right)
 478 {
 479     int sum;
 480
 481     sum  = coeff[0] * left[0] + coeff[3] * right[0]; //red in
 482     sum += coeff[1] * left[1] + coeff[4] * right[1]; //green in
 483     sum += coeff[2] * left[2] + coeff[5] * right[2]; //blue in
 484
 485     return av_clip_uint8(sum >> 16);
 486 }
 487
 488 static int filter_frame(AVFilterLink *inlink, AVFrame *inpicref)
 489 {
 490     AVFilterContext *ctx  = inlink->dst;
 491     Stereo3DContext *s = ctx->priv;
 492     AVFilterLink *outlink = ctx->outputs[0];
 493     AVFrame *out, *oleft, *oright, *ileft, *iright;
 494     int out_off_left[4], out_off_right[4];
 495     int in_off_left[4], in_off_right[4];
 496     int i;
 497
 498     switch (s->in.format) {
 499     case ALTERNATING_LR:
 500     case ALTERNATING_RL:
 501         if (!s->prev) {
 502             s->prev = inpicref;
 503             return 0;
 504         }
 505         ileft  = s->prev;
 506         iright = inpicref;
 507         if (s->in.format == ALTERNATING_RL)
 508             FFSWAP(AVFrame *, ileft, iright);
 509         break;
 510     default:
 511         ileft = iright = inpicref;
 512     };
 513
 514     out = oleft = oright = ff_get_video_buffer(outlink, outlink->w, outlink->h);
 515     if (!out) {
 516         av_frame_free(&s->prev);
 517         av_frame_free(&inpicref);
 518         return AVERROR(ENOMEM);
 519     }
 520     av_frame_copy_props(out, inpicref);
 521
 522     if (s->out.format == ALTERNATING_LR ||
 523         s->out.format == ALTERNATING_RL) {
 524         oright = ff_get_video_buffer(outlink, outlink->w, outlink->h);
 525         if (!oright) {
 526             av_frame_free(&oleft);
 527             av_frame_free(&s->prev);
 528             av_frame_free(&inpicref);
 529             return AVERROR(ENOMEM);
 530         }
 531         av_frame_copy_props(oright, inpicref);
 532     }
 533
 534     for (i = 0; i < 4; i++) {
 535         int hsub = i == 1 || i == 2 ? s->hsub : 0;
 536         int vsub = i == 1 || i == 2 ? s->vsub : 0;
 537         in_off_left[i]   = (FF_CEIL_RSHIFT(s->in.row_left,   vsub) + s->in.off_lstep)  * ileft->linesize[i]  + FF_CEIL_RSHIFT(s->in.off_left   * s->pixstep[i], hsub);
 538         in_off_right[i]  = (FF_CEIL_RSHIFT(s->in.row_right,  vsub) + s->in.off_rstep)  * iright->linesize[i] + FF_CEIL_RSHIFT(s->in.off_right  * s->pixstep[i], hsub);
 539         out_off_left[i]  = (FF_CEIL_RSHIFT(s->out.row_left,  vsub) + s->out.off_lstep) * oleft->linesize[i]  + FF_CEIL_RSHIFT(s->out.off_left  * s->pixstep[i], hsub);
 540         out_off_right[i] = (FF_CEIL_RSHIFT(s->out.row_right, vsub) + s->out.off_rstep) * oright->linesize[i] + FF_CEIL_RSHIFT(s->out.off_right * s->pixstep[i], hsub);
 541     }
 542
 543     switch (s->out.format) {
 544     case ALTERNATING_LR:
 545     case ALTERNATING_RL:
 546     case SIDE_BY_SIDE_LR:
 547     case SIDE_BY_SIDE_RL:
 548     case SIDE_BY_SIDE_2_LR:
 549     case SIDE_BY_SIDE_2_RL:
 550     case ABOVE_BELOW_LR:
 551     case ABOVE_BELOW_RL:
 552     case ABOVE_BELOW_2_LR:
 553     case ABOVE_BELOW_2_RL:
 554     case INTERLEAVE_ROWS_LR:
 555     case INTERLEAVE_ROWS_RL:
 556         for (i = 0; i < s->nb_planes; i++) {
 557             av_image_copy_plane(oleft->data[i] + out_off_left[i],
 558                                 oleft->linesize[i] * s->row_step,
 559                                 ileft->data[i] + in_off_left[i],
 560                                 ileft->linesize[i] * s->row_step,
 561                                 s->linesize[i], s->pheight[i]);
 562             av_image_copy_plane(oright->data[i] + out_off_right[i],
 563                                 oright->linesize[i] * s->row_step,
 564                                 iright->data[i] + in_off_right[i],
 565                                 iright->linesize[i] * s->row_step,
 566                                 s->linesize[i], s->pheight[i]);
 567         }
 568         break;
 569     case MONO_L:
 570         iright = ileft;
 571     case MONO_R:
 572         for (i = 0; i < s->nb_planes; i++) {
 573             av_image_copy_plane(out->data[i], out->linesize[i],
 574                                 iright->data[i] + in_off_left[i],
 575                                 iright->linesize[i],
 576                                 s->linesize[i], s->pheight[i]);
 577         }
 578         break;
 579     case ANAGLYPH_RB_GRAY:
 580     case ANAGLYPH_RG_GRAY:
 581     case ANAGLYPH_RC_GRAY:
 582     case ANAGLYPH_RC_HALF:
 583     case ANAGLYPH_RC_COLOR:
 584     case ANAGLYPH_RC_DUBOIS:
 585     case ANAGLYPH_GM_GRAY:
 586     case ANAGLYPH_GM_HALF:
 587     case ANAGLYPH_GM_COLOR:
 588     case ANAGLYPH_GM_DUBOIS:
 589     case ANAGLYPH_YB_GRAY:
 590     case ANAGLYPH_YB_HALF:
 591     case ANAGLYPH_YB_COLOR:
 592     case ANAGLYPH_YB_DUBOIS: {
 593         int x, y, il, ir, o;
 594         const uint8_t *lsrc = ileft->data[0];
 595         const uint8_t *rsrc = iright->data[0];
 596         uint8_t *dst = out->data[0];
 597         int out_width = s->out.width;
 598         const int **ana_matrix = s->ana_matrix;
 599
 600         for (y = 0; y < s->out.height; y++) {
 601             o   = out->linesize[0] * y;
 602             il  = in_off_left[0]  + y * ileft->linesize[0];
 603             ir  = in_off_right[0] + y * iright->linesize[0];
 604             for (x = 0; x < out_width; x++, il += 3, ir += 3, o+= 3) {
 605                 dst[o    ] = ana_convert(ana_matrix[0], lsrc + il, rsrc + ir);
 606                 dst[o + 1] = ana_convert(ana_matrix[1], lsrc + il, rsrc + ir);
 607                 dst[o + 2] = ana_convert(ana_matrix[2], lsrc + il, rsrc + ir);
 608             }
 609         }
 610         break;
 611     }
 612     default:
 613         av_assert0(0);
 614     }
 615
 616     av_frame_free(&inpicref);
 617     av_frame_free(&s->prev);
 618     if (oright != oleft) {
 619         if (s->out.format == ALTERNATING_LR)
 620             FFSWAP(AVFrame *, oleft, oright);
 621         oright->pts = outlink->frame_count * s->ts_unit;
 622         ff_filter_frame(outlink, oright);
 623         out = oleft;
 624         oleft->pts = outlink->frame_count * s->ts_unit;
 625     } else if (s->in.format == ALTERNATING_LR ||
 626                s->in.format == ALTERNATING_RL) {
 627         out->pts = outlink->frame_count * s->ts_unit;
 628     }
 629     return ff_filter_frame(outlink, out);
 630 }
 631
 632 static av_cold void uninit(AVFilterContext *ctx)
 633 {
 634     Stereo3DContext *s = ctx->priv;
 635
 636     av_frame_free(&s->prev);
 637 }
 638
 639 static const AVFilterPad stereo3d_inputs[] = {
 640     {
 641         .name         = "default",
 642         .type         = AVMEDIA_TYPE_VIDEO,
 643         .filter_frame = filter_frame,
 644     },
 645     { NULL }
 646 };
 647
 648 static const AVFilterPad stereo3d_outputs[] = {
 649     {
 650         .name         = "default",
 651         .type         = AVMEDIA_TYPE_VIDEO,
 652         .config_props = config_output,
 653     },
 654     { NULL }
 655 };
 656
 657 AVFilter ff_vf_stereo3d = {
 658     .name          = "stereo3d",
 659     .description   = NULL_IF_CONFIG_SMALL("Convert video stereoscopic 3D view."),
 660     .priv_size     = sizeof(Stereo3DContext),
 661     .uninit        = uninit,
 662     .query_formats = query_formats,
 663     .inputs        = stereo3d_inputs,
 664     .outputs       = stereo3d_outputs,
 665     .priv_class    = &stereo3d_class,
 666 };