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-2015 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/intreadwrite.h"
  25 #include "libavutil/opt.h"
  26 #include "libavutil/parseutils.h"
  27 #include "libavutil/pixdesc.h"
  28 #include "avfilter.h"
  29 #include "drawutils.h"
  30 #include "formats.h"
  31 #include "internal.h"
  32 #include "video.h"
  33 #include "stereo3d.h"
  34
  35 enum StereoCode {
  36     ANAGLYPH_RC_GRAY,   // anaglyph red/cyan gray
  37     ANAGLYPH_RC_HALF,   // anaglyph red/cyan half colored
  38     ANAGLYPH_RC_COLOR,  // anaglyph red/cyan colored
  39     ANAGLYPH_RC_DUBOIS, // anaglyph red/cyan dubois
  40     ANAGLYPH_GM_GRAY,   // anaglyph green/magenta gray
  41     ANAGLYPH_GM_HALF,   // anaglyph green/magenta half colored
  42     ANAGLYPH_GM_COLOR,  // anaglyph green/magenta colored
  43     ANAGLYPH_GM_DUBOIS, // anaglyph green/magenta dubois
  44     ANAGLYPH_YB_GRAY,   // anaglyph yellow/blue gray
  45     ANAGLYPH_YB_HALF,   // anaglyph yellow/blue half colored
  46     ANAGLYPH_YB_COLOR,  // anaglyph yellow/blue colored
  47     ANAGLYPH_YB_DUBOIS, // anaglyph yellow/blue dubois
  48     ANAGLYPH_RB_GRAY,   // anaglyph red/blue gray
  49     ANAGLYPH_RG_GRAY,   // anaglyph red/green gray
  50     MONO_L,             // mono output for debugging (left eye only)
  51     MONO_R,             // mono output for debugging (right eye only)
  52     INTERLEAVE_ROWS_LR, // row-interleave (left eye has top row)
  53     INTERLEAVE_ROWS_RL, // row-interleave (right eye has top row)
  54     SIDE_BY_SIDE_LR,    // side by side parallel (left eye left, right eye right)
  55     SIDE_BY_SIDE_RL,    // side by side crosseye (right eye left, left eye right)
  56     SIDE_BY_SIDE_2_LR,  // side by side parallel with half width resolution
  57     SIDE_BY_SIDE_2_RL,  // side by side crosseye with half width resolution
  58     ABOVE_BELOW_LR,     // above-below (left eye above, right eye below)
  59     ABOVE_BELOW_RL,     // above-below (right eye above, left eye below)
  60     ABOVE_BELOW_2_LR,   // above-below with half height resolution
  61     ABOVE_BELOW_2_RL,   // above-below with half height resolution
  62     ALTERNATING_LR,     // alternating frames (left eye first, right eye second)
  63     ALTERNATING_RL,     // alternating frames (right eye first, left eye second)
  64     CHECKERBOARD_LR,    // checkerboard pattern (left eye first, right eye second)
  65     CHECKERBOARD_RL,    // checkerboard pattern (right eye first, left eye second)
  66     INTERLEAVE_COLS_LR, // column-interleave (left eye first, right eye second)
  67     INTERLEAVE_COLS_RL, // column-interleave (right eye first, left eye second)
  68     HDMI,               // HDMI frame pack (left eye first, right eye second)
  69     STEREO_CODE_COUNT   // TODO: needs autodetection
  70 };
  71
  72 typedef struct StereoComponent {
  73     int format;                 ///< StereoCode
  74     int width, height;
  75     int off_left, off_right;
  76     int off_lstep, off_rstep;
  77     int row_left, row_right;
  78     int row_step;
  79 } StereoComponent;
  80
  81 static const int ana_coeff[][3][6] = {
  82   [ANAGLYPH_RB_GRAY]   =
  83     {{19595, 38470,  7471,     0,     0,     0},
  84      {    0,     0,     0,     0,     0,     0},
  85      {    0,     0,     0, 19595, 38470,  7471}},
  86   [ANAGLYPH_RG_GRAY]   =
  87     {{19595, 38470,  7471,     0,     0,     0},
  88      {    0,     0,     0, 19595, 38470,  7471},
  89      {    0,     0,     0,     0,     0,     0}},
  90   [ANAGLYPH_RC_GRAY]   =
  91     {{19595, 38470,  7471,     0,     0,     0},
  92      {    0,     0,     0, 19595, 38470,  7471},
  93      {    0,     0,     0, 19595, 38470,  7471}},
  94   [ANAGLYPH_RC_HALF]   =
  95     {{19595, 38470,  7471,     0,     0,     0},
  96      {    0,     0,     0,     0, 65536,     0},
  97      {    0,     0,     0,     0,     0, 65536}},
  98   [ANAGLYPH_RC_COLOR]  =
  99     {{65536,     0,     0,     0,     0,     0},
 100      {    0,     0,     0,     0, 65536,     0},
 101      {    0,     0,     0,     0,     0, 65536}},
 102   [ANAGLYPH_RC_DUBOIS] =
 103     {{29891, 32800, 11559, -2849, -5763,  -102},
 104      {-2627, -2479, -1033, 24804, 48080, -1209},
 105      { -997, -1350,  -358, -4729, -7403, 80373}},
 106   [ANAGLYPH_GM_GRAY]   =
 107     {{    0,     0,     0, 19595, 38470,  7471},
 108      {19595, 38470,  7471,     0,     0,     0},
 109      {    0,     0,     0, 19595, 38470,  7471}},
 110   [ANAGLYPH_GM_HALF]   =
 111     {{    0,     0,     0, 65536,     0,     0},
 112      {19595, 38470,  7471,     0,     0,     0},
 113      {    0,     0,     0,     0,     0, 65536}},
 114   [ANAGLYPH_GM_COLOR]  =
 115     {{    0,     0,     0, 65536,     0,     0},
 116      {    0, 65536,     0,     0,     0,     0},
 117      {    0,     0,     0,     0,     0, 65536}},
 118   [ANAGLYPH_GM_DUBOIS]  =
 119     {{-4063,-10354, -2556, 34669, 46203,  1573},
 120      {18612, 43778,  9372, -1049,  -983, -4260},
 121      { -983, -1769,  1376,   590,  4915, 61407}},
 122   [ANAGLYPH_YB_GRAY]   =
 123     {{    0,     0,     0, 19595, 38470,  7471},
 124      {    0,     0,     0, 19595, 38470,  7471},
 125      {19595, 38470,  7471,     0,     0,     0}},
 126   [ANAGLYPH_YB_HALF]   =
 127     {{    0,     0,     0, 65536,     0,     0},
 128      {    0,     0,     0,     0, 65536,     0},
 129      {19595, 38470,  7471,     0,     0,     0}},
 130   [ANAGLYPH_YB_COLOR]  =
 131     {{    0,     0,     0, 65536,     0,     0},
 132      {    0,     0,     0,     0, 65536,     0},
 133      {    0,     0, 65536,     0,     0,     0}},
 134   [ANAGLYPH_YB_DUBOIS] =
 135     {{65535,-12650,18451,   -987, -7590, -1049},
 136      {-1604, 56032, 4196,    370,  3826, -1049},
 137      {-2345,-10676, 1358,   5801, 11416, 56217}},
 138 };
 139
 140 typedef struct Stereo3DContext {
 141     const AVClass *class;
 142     StereoComponent in, out;
 143     int width, height;
 144     const int *ana_matrix[3];
 145     int nb_planes;
 146     int linesize[4];
 147     int pheight[4];
 148     int hsub, vsub;
 149     int pixstep[4];
 150     AVFrame *prev;
 151     int blanks;
 152     int in_off_left[4], in_off_right[4];
 153     Stereo3DDSPContext dsp;
 154 } Stereo3DContext;
 155
 156 #define OFFSET(x) offsetof(Stereo3DContext, x)
 157 #define FLAGS AV_OPT_FLAG_FILTERING_PARAM|AV_OPT_FLAG_VIDEO_PARAM
 158
 159 static const AVOption stereo3d_options[] = {
 160     { "in",    "set input format",  OFFSET(in.format),   AV_OPT_TYPE_INT,   {.i64=SIDE_BY_SIDE_LR}, INTERLEAVE_ROWS_LR, STEREO_CODE_COUNT-1, FLAGS, "in"},
 161     { "ab2l",  "above below half height left first",  0, AV_OPT_TYPE_CONST, {.i64=ABOVE_BELOW_2_LR},   0, 0, FLAGS, "in" },
 162     { "ab2r",  "above below half height right first", 0, AV_OPT_TYPE_CONST, {.i64=ABOVE_BELOW_2_RL},   0, 0, FLAGS, "in" },
 163     { "abl",   "above below left first",              0, AV_OPT_TYPE_CONST, {.i64=ABOVE_BELOW_LR},     0, 0, FLAGS, "in" },
 164     { "abr",   "above below right first",             0, AV_OPT_TYPE_CONST, {.i64=ABOVE_BELOW_RL},     0, 0, FLAGS, "in" },
 165     { "al",    "alternating frames left first",       0, AV_OPT_TYPE_CONST, {.i64=ALTERNATING_LR},     0, 0, FLAGS, "in" },
 166     { "ar",    "alternating frames right first",      0, AV_OPT_TYPE_CONST, {.i64=ALTERNATING_RL},     0, 0, FLAGS, "in" },
 167     { "sbs2l", "side by side half width left first",  0, AV_OPT_TYPE_CONST, {.i64=SIDE_BY_SIDE_2_LR},  0, 0, FLAGS, "in" },
 168     { "sbs2r", "side by side half width right first", 0, AV_OPT_TYPE_CONST, {.i64=SIDE_BY_SIDE_2_RL},  0, 0, FLAGS, "in" },
 169     { "sbsl",  "side by side left first",             0, AV_OPT_TYPE_CONST, {.i64=SIDE_BY_SIDE_LR},    0, 0, FLAGS, "in" },
 170     { "sbsr",  "side by side right first",            0, AV_OPT_TYPE_CONST, {.i64=SIDE_BY_SIDE_RL},    0, 0, FLAGS, "in" },
 171     { "irl",   "interleave rows left first",          0, AV_OPT_TYPE_CONST, {.i64=INTERLEAVE_ROWS_LR}, 0, 0, FLAGS, "in" },
 172     { "irr",   "interleave rows right first",         0, AV_OPT_TYPE_CONST, {.i64=INTERLEAVE_ROWS_RL}, 0, 0, FLAGS, "in" },
 173     { "out",   "set output format", OFFSET(out.format),  AV_OPT_TYPE_INT,   {.i64=ANAGLYPH_RC_DUBOIS}, 0, STEREO_CODE_COUNT-1, FLAGS, "out"},
 174     { "ab2l",  "above below half height left first",  0, AV_OPT_TYPE_CONST, {.i64=ABOVE_BELOW_2_LR},   0, 0, FLAGS, "out" },
 175     { "ab2r",  "above below half height right first", 0, AV_OPT_TYPE_CONST, {.i64=ABOVE_BELOW_2_RL},   0, 0, FLAGS, "out" },
 176     { "abl",   "above below left first",              0, AV_OPT_TYPE_CONST, {.i64=ABOVE_BELOW_LR},     0, 0, FLAGS, "out" },
 177     { "abr",   "above below right first",             0, AV_OPT_TYPE_CONST, {.i64=ABOVE_BELOW_RL},     0, 0, FLAGS, "out" },
 178     { "agmc",  "anaglyph green magenta color",        0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_GM_COLOR},  0, 0, FLAGS, "out" },
 179     { "agmd",  "anaglyph green magenta dubois",       0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_GM_DUBOIS}, 0, 0, FLAGS, "out" },
 180     { "agmg",  "anaglyph green magenta gray",         0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_GM_GRAY},   0, 0, FLAGS, "out" },
 181     { "agmh",  "anaglyph green magenta half color",   0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_GM_HALF},   0, 0, FLAGS, "out" },
 182     { "al",    "alternating frames left first",       0, AV_OPT_TYPE_CONST, {.i64=ALTERNATING_LR},     0, 0, FLAGS, "out" },
 183     { "ar",    "alternating frames right first",      0, AV_OPT_TYPE_CONST, {.i64=ALTERNATING_RL},     0, 0, FLAGS, "out" },
 184     { "arbg",  "anaglyph red blue gray",              0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_RB_GRAY},   0, 0, FLAGS, "out" },
 185     { "arcc",  "anaglyph red cyan color",             0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_RC_COLOR},  0, 0, FLAGS, "out" },
 186     { "arcd",  "anaglyph red cyan dubois",            0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_RC_DUBOIS}, 0, 0, FLAGS, "out" },
 187     { "arcg",  "anaglyph red cyan gray",              0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_RC_GRAY},   0, 0, FLAGS, "out" },
 188     { "arch",  "anaglyph red cyan half color",        0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_RC_HALF},   0, 0, FLAGS, "out" },
 189     { "argg",  "anaglyph red green gray",             0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_RG_GRAY},   0, 0, FLAGS, "out" },
 190     { "aybc",  "anaglyph yellow blue color",          0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_YB_COLOR},  0, 0, FLAGS, "out" },
 191     { "aybd",  "anaglyph yellow blue dubois",         0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_YB_DUBOIS}, 0, 0, FLAGS, "out" },
 192     { "aybg",  "anaglyph yellow blue gray",           0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_YB_GRAY},   0, 0, FLAGS, "out" },
 193     { "aybh",  "anaglyph yellow blue half color",     0, AV_OPT_TYPE_CONST, {.i64=ANAGLYPH_YB_HALF},   0, 0, FLAGS, "out" },
 194     { "irl",   "interleave rows left first",          0, AV_OPT_TYPE_CONST, {.i64=INTERLEAVE_ROWS_LR}, 0, 0, FLAGS, "out" },
 195     { "irr",   "interleave rows right first",         0, AV_OPT_TYPE_CONST, {.i64=INTERLEAVE_ROWS_RL}, 0, 0, FLAGS, "out" },
 196     { "ml",    "mono left",                           0, AV_OPT_TYPE_CONST, {.i64=MONO_L},             0, 0, FLAGS, "out" },
 197     { "mr",    "mono right",                          0, AV_OPT_TYPE_CONST, {.i64=MONO_R},             0, 0, FLAGS, "out" },
 198     { "sbs2l", "side by side half width left first",  0, AV_OPT_TYPE_CONST, {.i64=SIDE_BY_SIDE_2_LR},  0, 0, FLAGS, "out" },
 199     { "sbs2r", "side by side half width right first", 0, AV_OPT_TYPE_CONST, {.i64=SIDE_BY_SIDE_2_RL},  0, 0, FLAGS, "out" },
 200     { "sbsl",  "side by side left first",             0, AV_OPT_TYPE_CONST, {.i64=SIDE_BY_SIDE_LR},    0, 0, FLAGS, "out" },
 201     { "sbsr",  "side by side right first",            0, AV_OPT_TYPE_CONST, {.i64=SIDE_BY_SIDE_RL},    0, 0, FLAGS, "out" },
 202     { "chl",   "checkerboard left first",             0, AV_OPT_TYPE_CONST, {.i64=CHECKERBOARD_LR},    0, 0, FLAGS, "out" },
 203     { "chr",   "checkerboard right first",            0, AV_OPT_TYPE_CONST, {.i64=CHECKERBOARD_RL},    0, 0, FLAGS, "out" },
 204     { "icl",   "interleave columns left first",       0, AV_OPT_TYPE_CONST, {.i64=INTERLEAVE_COLS_LR}, 0, 0, FLAGS, "out" },
 205     { "icr",   "interleave columns right first",      0, AV_OPT_TYPE_CONST, {.i64=INTERLEAVE_COLS_RL}, 0, 0, FLAGS, "out" },
 206     { "hdmi",  "HDMI frame pack",                     0, AV_OPT_TYPE_CONST, {.i64=HDMI},               0, 0, FLAGS, "out" },
 207     { NULL }
 208 };
 209
 210 AVFILTER_DEFINE_CLASS(stereo3d);
 211
 212 static const enum AVPixelFormat anaglyph_pix_fmts[] = {
 213     AV_PIX_FMT_RGB24, AV_PIX_FMT_BGR24,
 214     AV_PIX_FMT_NONE
 215 };
 216
 217 static const enum AVPixelFormat other_pix_fmts[] = {
 218     AV_PIX_FMT_RGB24, AV_PIX_FMT_BGR24,
 219     AV_PIX_FMT_RGB48BE, AV_PIX_FMT_BGR48BE,
 220     AV_PIX_FMT_RGB48LE, AV_PIX_FMT_BGR48LE,
 221     AV_PIX_FMT_RGBA64BE, AV_PIX_FMT_BGRA64BE,
 222     AV_PIX_FMT_RGBA64LE, AV_PIX_FMT_BGRA64LE,
 223     AV_PIX_FMT_RGBA,  AV_PIX_FMT_BGRA,
 224     AV_PIX_FMT_ARGB,  AV_PIX_FMT_ABGR,
 225     AV_PIX_FMT_RGB0,  AV_PIX_FMT_BGR0,
 226     AV_PIX_FMT_0RGB,  AV_PIX_FMT_0BGR,
 227     AV_PIX_FMT_GBRP,
 228     AV_PIX_FMT_GBRP9BE,  AV_PIX_FMT_GBRP9LE,
 229     AV_PIX_FMT_GBRP10BE, AV_PIX_FMT_GBRP10LE,
 230     AV_PIX_FMT_GBRP12BE, AV_PIX_FMT_GBRP12LE,
 231     AV_PIX_FMT_GBRP14BE, AV_PIX_FMT_GBRP14LE,
 232     AV_PIX_FMT_GBRP16BE, AV_PIX_FMT_GBRP16LE,
 233     AV_PIX_FMT_YUV410P,
 234     AV_PIX_FMT_YUV411P,
 235     AV_PIX_FMT_YUV420P, AV_PIX_FMT_YUVA420P,
 236     AV_PIX_FMT_YUV422P, AV_PIX_FMT_YUVA422P,
 237     AV_PIX_FMT_YUV440P,
 238     AV_PIX_FMT_YUV444P, AV_PIX_FMT_YUVA444P,
 239     AV_PIX_FMT_YUVJ411P,
 240     AV_PIX_FMT_YUVJ420P,
 241     AV_PIX_FMT_YUVJ422P,
 242     AV_PIX_FMT_YUVJ440P,
 243     AV_PIX_FMT_YUVJ444P,
 244     AV_PIX_FMT_YUV420P9LE,  AV_PIX_FMT_YUVA420P9LE,
 245     AV_PIX_FMT_YUV420P9BE,  AV_PIX_FMT_YUVA420P9BE,
 246     AV_PIX_FMT_YUV422P9LE,  AV_PIX_FMT_YUVA422P9LE,
 247     AV_PIX_FMT_YUV422P9BE,  AV_PIX_FMT_YUVA422P9BE,
 248     AV_PIX_FMT_YUV444P9LE,  AV_PIX_FMT_YUVA444P9LE,
 249     AV_PIX_FMT_YUV444P9BE,  AV_PIX_FMT_YUVA444P9BE,
 250     AV_PIX_FMT_YUV420P10LE, AV_PIX_FMT_YUVA420P10LE,
 251     AV_PIX_FMT_YUV420P10BE, AV_PIX_FMT_YUVA420P10BE,
 252     AV_PIX_FMT_YUV422P10LE, AV_PIX_FMT_YUVA422P10LE,
 253     AV_PIX_FMT_YUV422P10BE, AV_PIX_FMT_YUVA422P10BE,
 254     AV_PIX_FMT_YUV444P10LE, AV_PIX_FMT_YUVA444P10LE,
 255     AV_PIX_FMT_YUV444P10BE, AV_PIX_FMT_YUVA444P10BE,
 256     AV_PIX_FMT_YUV420P12BE,  AV_PIX_FMT_YUV420P12LE,
 257     AV_PIX_FMT_YUV422P12BE,  AV_PIX_FMT_YUV422P12LE,
 258     AV_PIX_FMT_YUV444P12BE,  AV_PIX_FMT_YUV444P12LE,
 259     AV_PIX_FMT_YUV420P14BE,  AV_PIX_FMT_YUV420P14LE,
 260     AV_PIX_FMT_YUV422P14BE,  AV_PIX_FMT_YUV422P14LE,
 261     AV_PIX_FMT_YUV444P14BE,  AV_PIX_FMT_YUV444P14LE,
 262     AV_PIX_FMT_YUV420P16LE, AV_PIX_FMT_YUVA420P16LE,
 263     AV_PIX_FMT_YUV420P16BE, AV_PIX_FMT_YUVA420P16BE,
 264     AV_PIX_FMT_YUV422P16LE, AV_PIX_FMT_YUVA422P16LE,
 265     AV_PIX_FMT_YUV422P16BE, AV_PIX_FMT_YUVA422P16BE,
 266     AV_PIX_FMT_YUV444P16LE, AV_PIX_FMT_YUVA444P16LE,
 267     AV_PIX_FMT_YUV444P16BE, AV_PIX_FMT_YUVA444P16BE,
 268     AV_PIX_FMT_NONE
 269 };
 270
 271 static int query_formats(AVFilterContext *ctx)
 272 {
 273     Stereo3DContext *s = ctx->priv;
 274     const enum AVPixelFormat *pix_fmts;
 275     AVFilterFormats *fmts_list;
 276
 277     switch (s->out.format) {
 278     case ANAGLYPH_GM_COLOR:
 279     case ANAGLYPH_GM_DUBOIS:
 280     case ANAGLYPH_GM_GRAY:
 281     case ANAGLYPH_GM_HALF:
 282     case ANAGLYPH_RB_GRAY:
 283     case ANAGLYPH_RC_COLOR:
 284     case ANAGLYPH_RC_DUBOIS:
 285     case ANAGLYPH_RC_GRAY:
 286     case ANAGLYPH_RC_HALF:
 287     case ANAGLYPH_RG_GRAY:
 288     case ANAGLYPH_YB_COLOR:
 289     case ANAGLYPH_YB_DUBOIS:
 290     case ANAGLYPH_YB_GRAY:
 291     case ANAGLYPH_YB_HALF:
 292         pix_fmts = anaglyph_pix_fmts;
 293         break;
 294     default:
 295         pix_fmts = other_pix_fmts;
 296     }
 297
 298     fmts_list = ff_make_format_list(pix_fmts);
 299     if (!fmts_list)
 300         return AVERROR(ENOMEM);
 301     return ff_set_common_formats(ctx, fmts_list);
 302 }
 303
 304 static inline uint8_t ana_convert(const int *coeff, const uint8_t *left, const uint8_t *right)
 305 {
 306     int sum;
 307
 308     sum  = coeff[0] * left[0] + coeff[3] * right[0]; //red in
 309     sum += coeff[1] * left[1] + coeff[4] * right[1]; //green in
 310     sum += coeff[2] * left[2] + coeff[5] * right[2]; //blue in
 311
 312     return av_clip_uint8(sum >> 16);
 313 }
 314
 315 static void anaglyph(uint8_t *dst, uint8_t *lsrc, uint8_t *rsrc,
 316                      ptrdiff_t dst_linesize, ptrdiff_t l_linesize, ptrdiff_t r_linesize,
 317                      int width, int height,
 318                      const int *ana_matrix_r, const int *ana_matrix_g, const int *ana_matrix_b)
 319 {
 320     int x, y, o;
 321
 322     for (y = 0; y < height; y++) {
 323         for (o = 0, x = 0; x < width; x++, o+= 3) {
 324             dst[o    ] = ana_convert(ana_matrix_r, lsrc + o, rsrc + o);
 325             dst[o + 1] = ana_convert(ana_matrix_g, lsrc + o, rsrc + o);
 326             dst[o + 2] = ana_convert(ana_matrix_b, lsrc + o, rsrc + o);
 327         }
 328
 329         dst  += dst_linesize;
 330         lsrc += l_linesize;
 331         rsrc += r_linesize;
 332     }
 333 }
 334
 335 static int config_output(AVFilterLink *outlink)
 336 {
 337     AVFilterContext *ctx = outlink->src;
 338     AVFilterLink *inlink = ctx->inputs[0];
 339     Stereo3DContext *s = ctx->priv;
 340     AVRational aspect = inlink->sample_aspect_ratio;
 341     AVRational fps = inlink->frame_rate;
 342     AVRational tb = inlink->time_base;
 343     const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(outlink->format);
 344     int ret;
 345
 346     switch (s->in.format) {
 347     case SIDE_BY_SIDE_2_LR:
 348     case SIDE_BY_SIDE_LR:
 349     case SIDE_BY_SIDE_2_RL:
 350     case SIDE_BY_SIDE_RL:
 351         if (inlink->w & 1) {
 352             av_log(ctx, AV_LOG_ERROR, "width must be even\n");
 353             return AVERROR_INVALIDDATA;
 354         }
 355         break;
 356     case INTERLEAVE_ROWS_LR:
 357     case INTERLEAVE_ROWS_RL:
 358     case ABOVE_BELOW_2_LR:
 359     case ABOVE_BELOW_LR:
 360     case ABOVE_BELOW_2_RL:
 361     case ABOVE_BELOW_RL:
 362         if (inlink->h & 1) {
 363             av_log(ctx, AV_LOG_ERROR, "height must be even\n");
 364             return AVERROR_INVALIDDATA;
 365         }
 366         break;
 367     }
 368
 369     s->in.width     =
 370     s->width        = inlink->w;
 371     s->in.height    =
 372     s->height       = inlink->h;
 373     s->in.off_lstep =
 374     s->in.off_rstep =
 375     s->in.off_left  =
 376     s->in.off_right =
 377     s->in.row_left  =
 378     s->in.row_right = 0;
 379     s->in.row_step  = 1;
 380
 381     switch (s->in.format) {
 382     case SIDE_BY_SIDE_2_LR:
 383         aspect.num     *= 2;
 384     case SIDE_BY_SIDE_LR:
 385         s->width        = inlink->w / 2;
 386         s->in.off_right = s->width;
 387         break;
 388     case SIDE_BY_SIDE_2_RL:
 389         aspect.num     *= 2;
 390     case SIDE_BY_SIDE_RL:
 391         s->width        = inlink->w / 2;
 392         s->in.off_left  = s->width;
 393         break;
 394     case ABOVE_BELOW_2_LR:
 395         aspect.den     *= 2;
 396     case ABOVE_BELOW_LR:
 397         s->in.row_right =
 398         s->height       = inlink->h / 2;
 399         break;
 400     case ABOVE_BELOW_2_RL:
 401         aspect.den     *= 2;
 402     case ABOVE_BELOW_RL:
 403         s->in.row_left  =
 404         s->height       = inlink->h / 2;
 405         break;
 406     case ALTERNATING_RL:
 407     case ALTERNATING_LR:
 408         fps.den        *= 2;
 409         tb.num         *= 2;
 410         break;
 411     case INTERLEAVE_ROWS_LR:
 412     case INTERLEAVE_ROWS_RL:
 413         s->in.row_step  = 2;
 414         if (s->in.format == INTERLEAVE_ROWS_RL)
 415             s->in.off_lstep = 1;
 416         else
 417             s->in.off_rstep = 1;
 418         if (s->out.format != CHECKERBOARD_LR &&
 419             s->out.format != CHECKERBOARD_RL)
 420             s->height   = inlink->h / 2;
 421         break;
 422     default:
 423         av_log(ctx, AV_LOG_ERROR, "input format %d is not supported\n", s->in.format);
 424         return AVERROR(EINVAL);
 425     }
 426
 427     s->out.width     = s->width;
 428     s->out.height    = s->height;
 429     s->out.off_lstep =
 430     s->out.off_rstep =
 431     s->out.off_left  =
 432     s->out.off_right =
 433     s->out.row_left  =
 434     s->out.row_right = 0;
 435     s->out.row_step  = 1;
 436
 437     switch (s->out.format) {
 438     case ANAGLYPH_RB_GRAY:
 439     case ANAGLYPH_RG_GRAY:
 440     case ANAGLYPH_RC_GRAY:
 441     case ANAGLYPH_RC_HALF:
 442     case ANAGLYPH_RC_COLOR:
 443     case ANAGLYPH_RC_DUBOIS:
 444     case ANAGLYPH_GM_GRAY:
 445     case ANAGLYPH_GM_HALF:
 446     case ANAGLYPH_GM_COLOR:
 447     case ANAGLYPH_GM_DUBOIS:
 448     case ANAGLYPH_YB_GRAY:
 449     case ANAGLYPH_YB_HALF:
 450     case ANAGLYPH_YB_COLOR:
 451     case ANAGLYPH_YB_DUBOIS: {
 452         uint8_t rgba_map[4];
 453
 454         ff_fill_rgba_map(rgba_map, outlink->format);
 455         s->ana_matrix[rgba_map[0]] = &ana_coeff[s->out.format][0][0];
 456         s->ana_matrix[rgba_map[1]] = &ana_coeff[s->out.format][1][0];
 457         s->ana_matrix[rgba_map[2]] = &ana_coeff[s->out.format][2][0];
 458         break;
 459     }
 460     case SIDE_BY_SIDE_2_LR:
 461         aspect.den      *= 2;
 462     case SIDE_BY_SIDE_LR:
 463         s->out.width     = s->width * 2;
 464         s->out.off_right = s->width;
 465         break;
 466     case SIDE_BY_SIDE_2_RL:
 467         aspect.den      *= 2;
 468     case SIDE_BY_SIDE_RL:
 469         s->out.width     = s->width * 2;
 470         s->out.off_left  = s->width;
 471         break;
 472     case ABOVE_BELOW_2_LR:
 473         aspect.num      *= 2;
 474     case ABOVE_BELOW_LR:
 475         s->out.height    = s->height * 2;
 476         s->out.row_right = s->height;
 477         break;
 478     case HDMI:
 479         if (s->height != 720 && s->height != 1080) {
 480             av_log(ctx, AV_LOG_ERROR, "Only 720 and 1080 height supported\n");
 481             return AVERROR(EINVAL);
 482         }
 483
 484         s->blanks = s->height / 24;
 485         s->out.height    = s->height * 2 + s->blanks;
 486         s->out.row_right = s->height + s->blanks;
 487         break;
 488     case ABOVE_BELOW_2_RL:
 489         aspect.num      *= 2;
 490     case ABOVE_BELOW_RL:
 491         s->out.height    = s->height * 2;
 492         s->out.row_left  = s->height;
 493         break;
 494     case INTERLEAVE_ROWS_LR:
 495         s->in.row_step   = 1 + (s->in.format == INTERLEAVE_ROWS_RL);
 496         s->out.row_step  = 2;
 497         s->out.height    = s->height * 2;
 498         s->out.off_rstep = 1;
 499         break;
 500     case INTERLEAVE_ROWS_RL:
 501         s->in.row_step   = 1 + (s->in.format == INTERLEAVE_ROWS_LR);
 502         s->out.row_step  = 2;
 503         s->out.height    = s->height * 2;
 504         s->out.off_lstep = 1;
 505         break;
 506     case MONO_R:
 507         s->in.off_left   = s->in.off_right;
 508         s->in.row_left   = s->in.row_right;
 509         if (s->in.format == INTERLEAVE_ROWS_LR)
 510             FFSWAP(int, s->in.off_lstep, s->in.off_rstep);
 511         break;
 512     case MONO_L:
 513         if (s->in.format == INTERLEAVE_ROWS_RL)
 514             FFSWAP(int, s->in.off_lstep, s->in.off_rstep);
 515         break;
 516     case ALTERNATING_RL:
 517     case ALTERNATING_LR:
 518         fps.num         *= 2;
 519         tb.den          *= 2;
 520         break;
 521     case CHECKERBOARD_LR:
 522     case CHECKERBOARD_RL:
 523         s->out.width     = s->width * 2;
 524         break;
 525     case INTERLEAVE_COLS_LR:
 526     case INTERLEAVE_COLS_RL:
 527         s->out.width     = s->width * 2;
 528         break;
 529     default:
 530         av_log(ctx, AV_LOG_ERROR, "output format %d is not supported\n", s->out.format);
 531         return AVERROR(EINVAL);
 532     }
 533
 534     outlink->w = s->out.width;
 535     outlink->h = s->out.height;
 536     outlink->frame_rate = fps;
 537     outlink->time_base = tb;
 538     outlink->sample_aspect_ratio = aspect;
 539
 540     if ((ret = av_image_fill_linesizes(s->linesize, outlink->format, s->width)) < 0)
 541         return ret;
 542     s->nb_planes = av_pix_fmt_count_planes(outlink->format);
 543     av_image_fill_max_pixsteps(s->pixstep, NULL, desc);
 544     s->pheight[1] = s->pheight[2] = FF_CEIL_RSHIFT(s->height, desc->log2_chroma_h);
 545     s->pheight[0] = s->pheight[3] = s->height;
 546     s->hsub = desc->log2_chroma_w;
 547     s->vsub = desc->log2_chroma_h;
 548
 549     s->dsp.anaglyph = anaglyph;
 550     if (ARCH_X86)
 551         ff_stereo3d_init_x86(&s->dsp);
 552
 553     return 0;
 554 }
 555
 556 typedef struct ThreadData {
 557     AVFrame *ileft, *iright;
 558     AVFrame *out;
 559 } ThreadData;
 560
 561 static int filter_slice(AVFilterContext *ctx, void *arg, int jobnr, int nb_jobs)
 562 {
 563     Stereo3DContext *s = ctx->priv;
 564     ThreadData *td = arg;
 565     AVFrame *ileft = td->ileft;
 566     AVFrame *iright = td->iright;
 567     AVFrame *out = td->out;
 568     int height = s->out.height;
 569     int start = (height *  jobnr   ) / nb_jobs;
 570     int end   = (height * (jobnr+1)) / nb_jobs;
 571     const int **ana_matrix = s->ana_matrix;
 572
 573     s->dsp.anaglyph(out->data[0] + out->linesize[0] * start,
 574              ileft ->data[0] + s->in_off_left [0]  + ileft->linesize[0] * start * s->in.row_step,
 575              iright->data[0] + s->in_off_right[0] + iright->linesize[0] * start * s->in.row_step,
 576              out->linesize[0],
 577              ileft->linesize[0] * s->in.row_step,
 578              iright->linesize[0] * s->in.row_step,
 579              s->out.width, end - start,
 580              ana_matrix[0], ana_matrix[1], ana_matrix[2]);
 581
 582     return 0;
 583 }
 584
 585 static int filter_frame(AVFilterLink *inlink, AVFrame *inpicref)
 586 {
 587     AVFilterContext *ctx  = inlink->dst;
 588     Stereo3DContext *s = ctx->priv;
 589     AVFilterLink *outlink = ctx->outputs[0];
 590     AVFrame *out, *oleft, *oright, *ileft, *iright;
 591     int out_off_left[4], out_off_right[4];
 592     int i, ret;
 593
 594     if (s->in.format == s->out.format)
 595         return ff_filter_frame(outlink, inpicref);
 596
 597     switch (s->out.format) {
 598     case ALTERNATING_LR:
 599     case ALTERNATING_RL:
 600         if (!s->prev) {
 601             s->prev = inpicref;
 602             return 0;
 603         }
 604         break;
 605     };
 606
 607     switch (s->in.format) {
 608     case ALTERNATING_LR:
 609     case ALTERNATING_RL:
 610         if (!s->prev) {
 611             s->prev = inpicref;
 612             return 0;
 613         }
 614         ileft  = s->prev;
 615         iright = inpicref;
 616         if (s->in.format == ALTERNATING_RL)
 617             FFSWAP(AVFrame *, ileft, iright);
 618         break;
 619     default:
 620         ileft = iright = inpicref;
 621     };
 622
 623     if ((s->out.format == ALTERNATING_LR ||
 624          s->out.format == ALTERNATING_RL) &&
 625         (s->in.format == SIDE_BY_SIDE_LR ||
 626          s->in.format == SIDE_BY_SIDE_RL ||
 627          s->in.format == SIDE_BY_SIDE_2_LR ||
 628          s->in.format == SIDE_BY_SIDE_2_RL ||
 629          s->in.format == ABOVE_BELOW_LR ||
 630          s->in.format == ABOVE_BELOW_RL ||
 631          s->in.format == ABOVE_BELOW_2_LR ||
 632          s->in.format == ABOVE_BELOW_2_RL ||
 633          s->in.format == INTERLEAVE_ROWS_LR ||
 634          s->in.format == INTERLEAVE_ROWS_RL)) {
 635         oright = av_frame_clone(s->prev);
 636         oleft  = av_frame_clone(s->prev);
 637         if (!oright || !oleft) {
 638             av_frame_free(&oright);
 639             av_frame_free(&oleft);
 640             av_frame_free(&s->prev);
 641             av_frame_free(&inpicref);
 642             return AVERROR(ENOMEM);
 643         }
 644     } else if ((s->out.format == MONO_L ||
 645                 s->out.format == MONO_R) &&
 646         (s->in.format == SIDE_BY_SIDE_LR ||
 647          s->in.format == SIDE_BY_SIDE_RL ||
 648          s->in.format == SIDE_BY_SIDE_2_LR ||
 649          s->in.format == SIDE_BY_SIDE_2_RL ||
 650          s->in.format == ABOVE_BELOW_LR ||
 651          s->in.format == ABOVE_BELOW_RL ||
 652          s->in.format == ABOVE_BELOW_2_LR ||
 653          s->in.format == ABOVE_BELOW_2_RL ||
 654          s->in.format == INTERLEAVE_ROWS_LR ||
 655          s->in.format == INTERLEAVE_ROWS_RL)) {
 656         out = oleft = oright = av_frame_clone(inpicref);
 657         if (!out) {
 658             av_frame_free(&s->prev);
 659             av_frame_free(&inpicref);
 660             return AVERROR(ENOMEM);
 661         }
 662     } else if ((s->out.format == MONO_L && s->in.format == ALTERNATING_LR) ||
 663                (s->out.format == MONO_R && s->in.format == ALTERNATING_RL)) {
 664         s->prev->pts /= 2;
 665         ret = ff_filter_frame(outlink, s->prev);
 666         av_frame_free(&inpicref);
 667         s->prev = NULL;
 668         return ret;
 669     } else if ((s->out.format == MONO_L && s->in.format == ALTERNATING_RL) ||
 670                (s->out.format == MONO_R && s->in.format == ALTERNATING_LR)) {
 671         av_frame_free(&s->prev);
 672         inpicref->pts /= 2;
 673         return ff_filter_frame(outlink, inpicref);
 674     } else if ((s->out.format == ALTERNATING_LR && s->in.format == ALTERNATING_RL) ||
 675                (s->out.format == ALTERNATING_RL && s->in.format == ALTERNATING_LR)) {
 676         FFSWAP(int64_t, s->prev->pts, inpicref->pts);
 677         ff_filter_frame(outlink, inpicref);
 678         ret = ff_filter_frame(outlink, s->prev);
 679         s->prev = NULL;
 680         return ret;
 681     } else {
 682         out = oleft = oright = ff_get_video_buffer(outlink, outlink->w, outlink->h);
 683         if (!out) {
 684             av_frame_free(&s->prev);
 685             av_frame_free(&inpicref);
 686             return AVERROR(ENOMEM);
 687         }
 688         av_frame_copy_props(out, inpicref);
 689
 690         if (s->out.format == ALTERNATING_LR ||
 691             s->out.format == ALTERNATING_RL) {
 692             oright = ff_get_video_buffer(outlink, outlink->w, outlink->h);
 693             if (!oright) {
 694                 av_frame_free(&oleft);
 695                 av_frame_free(&s->prev);
 696                 av_frame_free(&inpicref);
 697                 return AVERROR(ENOMEM);
 698             }
 699             av_frame_copy_props(oright, s->prev);
 700         }
 701     }
 702
 703     for (i = 0; i < 4; i++) {
 704         int hsub = i == 1 || i == 2 ? s->hsub : 0;
 705         int vsub = i == 1 || i == 2 ? s->vsub : 0;
 706         s->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);
 707         s->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);
 708         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);
 709         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);
 710     }
 711
 712     switch (s->out.format) {
 713     case ALTERNATING_LR:
 714     case ALTERNATING_RL:
 715         switch (s->in.format) {
 716         case INTERLEAVE_ROWS_LR:
 717         case INTERLEAVE_ROWS_RL:
 718             for (i = 0; i < s->nb_planes; i++) {
 719                 oleft->linesize[i]  *= 2;
 720                 oright->linesize[i] *= 2;
 721             }
 722         case ABOVE_BELOW_LR:
 723         case ABOVE_BELOW_RL:
 724         case ABOVE_BELOW_2_LR:
 725         case ABOVE_BELOW_2_RL:
 726         case SIDE_BY_SIDE_LR:
 727         case SIDE_BY_SIDE_RL:
 728         case SIDE_BY_SIDE_2_LR:
 729         case SIDE_BY_SIDE_2_RL:
 730             oleft->width   = outlink->w;
 731             oright->width  = outlink->w;
 732             oleft->height  = outlink->h;
 733             oright->height = outlink->h;
 734
 735             for (i = 0; i < s->nb_planes; i++) {
 736                 oleft->data[i]  += s->in_off_left[i];
 737                 oright->data[i] += s->in_off_right[i];
 738             }
 739             break;
 740         default:
 741             goto copy;
 742             break;
 743         }
 744         break;
 745     case HDMI:
 746         for (i = 0; i < s->nb_planes; i++) {
 747             int j, h = s->height >> ((i == 1 || i == 2) ? s->vsub : 0);
 748             int b = (s->blanks) >> ((i == 1 || i == 2) ? s->vsub : 0);
 749
 750             for (j = h; j < h + b; j++)
 751                 memset(oleft->data[i] + j * s->linesize[i], 0, s->linesize[i]);
 752         }
 753     case SIDE_BY_SIDE_LR:
 754     case SIDE_BY_SIDE_RL:
 755     case SIDE_BY_SIDE_2_LR:
 756     case SIDE_BY_SIDE_2_RL:
 757     case ABOVE_BELOW_LR:
 758     case ABOVE_BELOW_RL:
 759     case ABOVE_BELOW_2_LR:
 760     case ABOVE_BELOW_2_RL:
 761     case INTERLEAVE_ROWS_LR:
 762     case INTERLEAVE_ROWS_RL:
 763 copy:
 764         for (i = 0; i < s->nb_planes; i++) {
 765             av_image_copy_plane(oleft->data[i] + out_off_left[i],
 766                                 oleft->linesize[i] * s->out.row_step,
 767                                 ileft->data[i] + s->in_off_left[i],
 768                                 ileft->linesize[i] * s->in.row_step,
 769                                 s->linesize[i], s->pheight[i]);
 770             av_image_copy_plane(oright->data[i] + out_off_right[i],
 771                                 oright->linesize[i] * s->out.row_step,
 772                                 iright->data[i] + s->in_off_right[i],
 773                                 iright->linesize[i] * s->in.row_step,
 774                                 s->linesize[i], s->pheight[i]);
 775         }
 776         break;
 777     case MONO_L:
 778         iright = ileft;
 779     case MONO_R:
 780         switch (s->in.format) {
 781         case INTERLEAVE_ROWS_LR:
 782         case INTERLEAVE_ROWS_RL:
 783             for (i = 0; i < s->nb_planes; i++) {
 784                 out->linesize[i] *= 2;
 785             }
 786         case ABOVE_BELOW_LR:
 787         case ABOVE_BELOW_RL:
 788         case ABOVE_BELOW_2_LR:
 789         case ABOVE_BELOW_2_RL:
 790         case SIDE_BY_SIDE_LR:
 791         case SIDE_BY_SIDE_RL:
 792         case SIDE_BY_SIDE_2_LR:
 793         case SIDE_BY_SIDE_2_RL:
 794             out->width  = outlink->w;
 795             out->height = outlink->h;
 796
 797             for (i = 0; i < s->nb_planes; i++) {
 798                 out->data[i] += s->in_off_left[i];
 799             }
 800             break;
 801         default:
 802             for (i = 0; i < s->nb_planes; i++) {
 803                 av_image_copy_plane(out->data[i], out->linesize[i],
 804                                     iright->data[i] + s->in_off_left[i],
 805                                     iright->linesize[i] * s->in.row_step,
 806                                     s->linesize[i], s->pheight[i]);
 807             }
 808             break;
 809         }
 810         break;
 811     case ANAGLYPH_RB_GRAY:
 812     case ANAGLYPH_RG_GRAY:
 813     case ANAGLYPH_RC_GRAY:
 814     case ANAGLYPH_RC_HALF:
 815     case ANAGLYPH_RC_COLOR:
 816     case ANAGLYPH_RC_DUBOIS:
 817     case ANAGLYPH_GM_GRAY:
 818     case ANAGLYPH_GM_HALF:
 819     case ANAGLYPH_GM_COLOR:
 820     case ANAGLYPH_GM_DUBOIS:
 821     case ANAGLYPH_YB_GRAY:
 822     case ANAGLYPH_YB_HALF:
 823     case ANAGLYPH_YB_COLOR:
 824     case ANAGLYPH_YB_DUBOIS: {
 825         ThreadData td;
 826
 827         td.ileft = ileft; td.iright = iright; td.out = out;
 828         ctx->internal->execute(ctx, filter_slice, &td, NULL,
 829                                FFMIN(s->out.height, ctx->graph->nb_threads));
 830         break;
 831     }
 832     case CHECKERBOARD_RL:
 833     case CHECKERBOARD_LR:
 834         for (i = 0; i < s->nb_planes; i++) {
 835             int x, y;
 836
 837             for (y = 0; y < s->pheight[i]; y++) {
 838                 uint8_t *dst = out->data[i] + out->linesize[i] * y;
 839                 uint8_t *left = ileft->data[i] + ileft->linesize[i] * y + s->in_off_left[i];
 840                 uint8_t *right = iright->data[i] + iright->linesize[i] * y + s->in_off_right[i];
 841                 int p, b;
 842
 843                 if (s->out.format == CHECKERBOARD_RL)
 844                     FFSWAP(uint8_t*, left, right);
 845                 switch (s->pixstep[i]) {
 846                 case 1:
 847                     for (x = 0, b = 0, p = 0; x < s->linesize[i] * 2; x+=2, p++, b++) {
 848                         dst[x  ] = (b&1) == (y&1) ? left[p] : right[p];
 849                         dst[x+1] = (b&1) != (y&1) ? left[p] : right[p];
 850                     }
 851                     break;
 852                 case 2:
 853                     for (x = 0, b = 0, p = 0; x < s->linesize[i] * 2; x+=4, p+=2, b++) {
 854                         AV_WN16(&dst[x  ], (b&1) == (y&1) ? AV_RN16(&left[p]) : AV_RN16(&right[p]));
 855                         AV_WN16(&dst[x+2], (b&1) != (y&1) ? AV_RN16(&left[p]) : AV_RN16(&right[p]));
 856                     }
 857                     break;
 858                 case 3:
 859                     for (x = 0, b = 0, p = 0; x < s->linesize[i] * 2; x+=6, p+=3, b++) {
 860                         AV_WB24(&dst[x  ], (b&1) == (y&1) ? AV_RB24(&left[p]) : AV_RB24(&right[p]));
 861                         AV_WB24(&dst[x+3], (b&1) != (y&1) ? AV_RB24(&left[p]) : AV_RB24(&right[p]));
 862                     }
 863                     break;
 864                 case 4:
 865                     for (x = 0, b = 0, p = 0; x < s->linesize[i] * 2; x+=8, p+=4, b++) {
 866                         AV_WN32(&dst[x  ], (b&1) == (y&1) ? AV_RN32(&left[p]) : AV_RN32(&right[p]));
 867                         AV_WN32(&dst[x+4], (b&1) != (y&1) ? AV_RN32(&left[p]) : AV_RN32(&right[p]));
 868                     }
 869                     break;
 870                 case 6:
 871                     for (x = 0, b = 0, p = 0; x < s->linesize[i] * 2; x+=12, p+=6, b++) {
 872                         AV_WB48(&dst[x  ], (b&1) == (y&1) ? AV_RB48(&left[p]) : AV_RB48(&right[p]));
 873                         AV_WB48(&dst[x+6], (b&1) != (y&1) ? AV_RB48(&left[p]) : AV_RB48(&right[p]));
 874                     }
 875                     break;
 876                 case 8:
 877                     for (x = 0, b = 0, p = 0; x < s->linesize[i] * 2; x+=16, p+=8, b++) {
 878                         AV_WN64(&dst[x  ], (b&1) == (y&1) ? AV_RN64(&left[p]) : AV_RN64(&right[p]));
 879                         AV_WN64(&dst[x+8], (b&1) != (y&1) ? AV_RN64(&left[p]) : AV_RN64(&right[p]));
 880                     }
 881                     break;
 882                 }
 883             }
 884         }
 885         break;
 886     case INTERLEAVE_COLS_LR:
 887     case INTERLEAVE_COLS_RL:
 888         for (i = 0; i < s->nb_planes; i++) {
 889             int x, y;
 890
 891             for (y = 0; y < s->pheight[i]; y++) {
 892                 uint8_t *dst = out->data[i] + out->linesize[i] * y;
 893                 uint8_t *left = ileft->data[i] + ileft->linesize[i] * y * s->in.row_step + s->in_off_left[i];
 894                 uint8_t *right = iright->data[i] + iright->linesize[i] * y * s->in.row_step + s->in_off_right[i];
 895                 int p, b;
 896
 897                 if (s->out.format == INTERLEAVE_COLS_LR)
 898                     FFSWAP(uint8_t*, left, right);
 899
 900                 switch (s->pixstep[i]) {
 901                 case 1:
 902                     for (x = 0, b = 0, p = 0; x < s->linesize[i] * 2; x+=2, p++, b++) {
 903                         dst[x  ] =   b&1  ? left[p] : right[p];
 904                         dst[x+1] = !(b&1) ? left[p] : right[p];
 905                     }
 906                     break;
 907                 case 2:
 908                     for (x = 0, b = 0, p = 0; x < s->linesize[i] * 2; x+=4, p+=2, b++) {
 909                         AV_WN16(&dst[x  ],   b&1  ? AV_RN16(&left[p]) : AV_RN16(&right[p]));
 910                         AV_WN16(&dst[x+2], !(b&1) ? AV_RN16(&left[p]) : AV_RN16(&right[p]));
 911                     }
 912                     break;
 913                 case 3:
 914                     for (x = 0, b = 0, p = 0; x < s->linesize[i] * 2; x+=6, p+=3, b++) {
 915                         AV_WB24(&dst[x  ],   b&1  ? AV_RB24(&left[p]) : AV_RB24(&right[p]));
 916                         AV_WB24(&dst[x+3], !(b&1) ? AV_RB24(&left[p]) : AV_RB24(&right[p]));
 917                     }
 918                     break;
 919                 case 4:
 920                     for (x = 0, b = 0, p = 0; x < s->linesize[i] * 2; x+=8, p+=4, b++) {
 921                         AV_WN32(&dst[x  ],   b&1  ? AV_RN32(&left[p]) : AV_RN32(&right[p]));
 922                         AV_WN32(&dst[x+4], !(b&1) ? AV_RN32(&left[p]) : AV_RN32(&right[p]));
 923                     }
 924                     break;
 925                 case 6:
 926                     for (x = 0, b = 0, p = 0; x < s->linesize[i] * 2; x+=12, p+=6, b++) {
 927                         AV_WB48(&dst[x  ],   b&1  ? AV_RB48(&left[p]) : AV_RB48(&right[p]));
 928                         AV_WB48(&dst[x+6], !(b&1) ? AV_RB48(&left[p]) : AV_RB48(&right[p]));
 929                     }
 930                     break;
 931                 case 8:
 932                     for (x = 0, b = 0, p = 0; x < s->linesize[i] * 2; x+=16, p+=8, b++) {
 933                         AV_WN64(&dst[x  ],   b&1 ?  AV_RN64(&left[p]) : AV_RN64(&right[p]));
 934                         AV_WN64(&dst[x+8], !(b&1) ? AV_RN64(&left[p]) : AV_RN64(&right[p]));
 935                     }
 936                     break;
 937                 }
 938             }
 939         }
 940         break;
 941     default:
 942         av_assert0(0);
 943     }
 944
 945     if (oright != oleft) {
 946         if (s->out.format == ALTERNATING_LR)
 947             FFSWAP(AVFrame *, oleft, oright);
 948         oright->pts = s->prev->pts * 2;
 949         ff_filter_frame(outlink, oright);
 950         out = oleft;
 951         oleft->pts = s->prev->pts + inpicref->pts;
 952         av_frame_free(&s->prev);
 953         s->prev = inpicref;
 954     } else if (s->in.format == ALTERNATING_LR ||
 955                s->in.format == ALTERNATING_RL) {
 956         out->pts = s->prev->pts / 2;
 957         av_frame_free(&s->prev);
 958         av_frame_free(&inpicref);
 959     } else {
 960         av_frame_free(&s->prev);
 961         av_frame_free(&inpicref);
 962     }
 963     return ff_filter_frame(outlink, out);
 964 }
 965
 966 static av_cold void uninit(AVFilterContext *ctx)
 967 {
 968     Stereo3DContext *s = ctx->priv;
 969
 970     av_frame_free(&s->prev);
 971 }
 972
 973 static const AVFilterPad stereo3d_inputs[] = {
 974     {
 975         .name         = "default",
 976         .type         = AVMEDIA_TYPE_VIDEO,
 977         .filter_frame = filter_frame,
 978     },
 979     { NULL }
 980 };
 981
 982 static const AVFilterPad stereo3d_outputs[] = {
 983     {
 984         .name         = "default",
 985         .type         = AVMEDIA_TYPE_VIDEO,
 986         .config_props = config_output,
 987     },
 988     { NULL }
 989 };
 990
 991 AVFilter ff_vf_stereo3d = {
 992     .name          = "stereo3d",
 993     .description   = NULL_IF_CONFIG_SMALL("Convert video stereoscopic 3D view."),
 994     .priv_size     = sizeof(Stereo3DContext),
 995     .uninit        = uninit,
 996     .query_formats = query_formats,
 997     .inputs        = stereo3d_inputs,
 998     .outputs       = stereo3d_outputs,
 999     .priv_class    = &stereo3d_class,
1000     .flags         = AVFILTER_FLAG_SLICE_THREADS,
1001 };