]> git.sesse.net Git - ffmpeg/blob - libavfilter/vf_crop.c
Merge remote-tracking branch 'qatar/master'
[ffmpeg] / libavfilter / vf_crop.c
1 /*
2  * Copyright (c) 2007 Bobby Bingham
3  *
4  * This file is part of FFmpeg.
5  *
6  * FFmpeg is free software; you can redistribute it and/or
7  * modify it under the terms of the GNU Lesser General Public
8  * License as published by the Free Software Foundation; either
9  * version 2.1 of the License, or (at your option) any later version.
10  *
11  * FFmpeg is distributed in the hope that it will be useful,
12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14  * Lesser General Public License for more details.
15  *
16  * You should have received a copy of the GNU Lesser General Public
17  * License along with FFmpeg; if not, write to the Free Software
18  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
19  */
20
21 /**
22  * @file
23  * video crop filter
24  */
25
26 /* #define DEBUG */
27
28 #include "avfilter.h"
29 #include "libavutil/eval.h"
30 #include "libavutil/avstring.h"
31 #include "libavutil/libm.h"
32 #include "libavutil/imgutils.h"
33 #include "libavutil/mathematics.h"
34
35 static const char *var_names[] = {
36     "E",
37     "PHI",
38     "PI",
39     "in_w", "iw",   ///< width  of the input video
40     "in_h", "ih",   ///< height of the input video
41     "out_w", "ow",  ///< width  of the cropped video
42     "out_h", "oh",  ///< height of the cropped video
43     "a",
44     "sar",
45     "dar",
46     "hsub",
47     "vsub",
48     "x",
49     "y",
50     "n",            ///< number of frame
51     "pos",          ///< position in the file
52     "t",            ///< timestamp expressed in seconds
53     NULL
54 };
55
56 enum var_name {
57     VAR_E,
58     VAR_PHI,
59     VAR_PI,
60     VAR_IN_W,  VAR_IW,
61     VAR_IN_H,  VAR_IH,
62     VAR_OUT_W, VAR_OW,
63     VAR_OUT_H, VAR_OH,
64     VAR_A,
65     VAR_SAR,
66     VAR_DAR,
67     VAR_HSUB,
68     VAR_VSUB,
69     VAR_X,
70     VAR_Y,
71     VAR_N,
72     VAR_POS,
73     VAR_T,
74     VAR_VARS_NB
75 };
76
77 typedef struct {
78     int  x;             ///< x offset of the non-cropped area with respect to the input area
79     int  y;             ///< y offset of the non-cropped area with respect to the input area
80     int  w;             ///< width of the cropped area
81     int  h;             ///< height of the cropped area
82
83     int max_step[4];    ///< max pixel step for each plane, expressed as a number of bytes
84     int hsub, vsub;     ///< chroma subsampling
85     char x_expr[256], y_expr[256], ow_expr[256], oh_expr[256];
86     AVExpr *x_pexpr, *y_pexpr;  /* parsed expressions for x and y */
87     double var_values[VAR_VARS_NB];
88 } CropContext;
89
90 static int query_formats(AVFilterContext *ctx)
91 {
92     static const enum PixelFormat pix_fmts[] = {
93         PIX_FMT_RGB48BE,      PIX_FMT_RGB48LE,
94         PIX_FMT_BGR48BE,      PIX_FMT_BGR48LE,
95         PIX_FMT_ARGB,         PIX_FMT_RGBA,
96         PIX_FMT_ABGR,         PIX_FMT_BGRA,
97         PIX_FMT_RGB24,        PIX_FMT_BGR24,
98         PIX_FMT_RGB565BE,     PIX_FMT_RGB565LE,
99         PIX_FMT_RGB555BE,     PIX_FMT_RGB555LE,
100         PIX_FMT_BGR565BE,     PIX_FMT_BGR565LE,
101         PIX_FMT_BGR555BE,     PIX_FMT_BGR555LE,
102         PIX_FMT_GRAY16BE,     PIX_FMT_GRAY16LE,
103         PIX_FMT_YUV420P16LE,  PIX_FMT_YUV420P16BE,
104         PIX_FMT_YUV422P16LE,  PIX_FMT_YUV422P16BE,
105         PIX_FMT_YUV444P16LE,  PIX_FMT_YUV444P16BE,
106         PIX_FMT_YUV444P,      PIX_FMT_YUV422P,
107         PIX_FMT_YUV420P,      PIX_FMT_YUV411P,
108         PIX_FMT_YUV410P,      PIX_FMT_YUV440P,
109         PIX_FMT_YUVJ444P,     PIX_FMT_YUVJ422P,
110         PIX_FMT_YUVJ420P,     PIX_FMT_YUVJ440P,
111         PIX_FMT_YUVA420P,
112         PIX_FMT_RGB8,         PIX_FMT_BGR8,
113         PIX_FMT_RGB4_BYTE,    PIX_FMT_BGR4_BYTE,
114         PIX_FMT_PAL8,         PIX_FMT_GRAY8,
115         PIX_FMT_NONE
116     };
117
118     avfilter_set_common_pixel_formats(ctx, avfilter_make_format_list(pix_fmts));
119
120     return 0;
121 }
122
123 static av_cold int init(AVFilterContext *ctx, const char *args, void *opaque)
124 {
125     CropContext *crop = ctx->priv;
126
127     av_strlcpy(crop->ow_expr, "iw", sizeof(crop->ow_expr));
128     av_strlcpy(crop->oh_expr, "ih", sizeof(crop->oh_expr));
129     av_strlcpy(crop->x_expr, "(in_w-out_w)/2", sizeof(crop->x_expr));
130     av_strlcpy(crop->y_expr, "(in_h-out_h)/2", sizeof(crop->y_expr));
131
132     if (args)
133         sscanf(args, "%255[^:]:%255[^:]:%255[^:]:%255[^:]", crop->ow_expr, crop->oh_expr, crop->x_expr, crop->y_expr);
134
135     return 0;
136 }
137
138 static av_cold void uninit(AVFilterContext *ctx)
139 {
140     CropContext *crop = ctx->priv;
141
142     av_expr_free(crop->x_pexpr); crop->x_pexpr = NULL;
143     av_expr_free(crop->y_pexpr); crop->y_pexpr = NULL;
144 }
145
146 static inline int normalize_double(int *n, double d)
147 {
148     int ret = 0;
149
150     if (isnan(d)) {
151         ret = AVERROR(EINVAL);
152     } else if (d > INT_MAX || d < INT_MIN) {
153         *n = d > INT_MAX ? INT_MAX : INT_MIN;
154         ret = AVERROR(EINVAL);
155     } else
156         *n = round(d);
157
158     return ret;
159 }
160
161 static int config_input(AVFilterLink *link)
162 {
163     AVFilterContext *ctx = link->dst;
164     CropContext *crop = ctx->priv;
165     const AVPixFmtDescriptor *pix_desc = &av_pix_fmt_descriptors[link->format];
166     int ret;
167     const char *expr;
168     double res;
169
170     crop->var_values[VAR_E]     = M_E;
171     crop->var_values[VAR_PHI]   = M_PHI;
172     crop->var_values[VAR_PI]    = M_PI;
173     crop->var_values[VAR_IN_W]  = crop->var_values[VAR_IW] = ctx->inputs[0]->w;
174     crop->var_values[VAR_IN_H]  = crop->var_values[VAR_IH] = ctx->inputs[0]->h;
175     crop->var_values[VAR_A]     = (float) link->w / link->h;
176     crop->var_values[VAR_SAR]   = link->sample_aspect_ratio.num ? av_q2d(link->sample_aspect_ratio) : 1;
177     crop->var_values[VAR_DAR]   = crop->var_values[VAR_A] * crop->var_values[VAR_SAR];
178     crop->var_values[VAR_HSUB]  = 1<<pix_desc->log2_chroma_w;
179     crop->var_values[VAR_VSUB]  = 1<<pix_desc->log2_chroma_h;
180     crop->var_values[VAR_X]     = NAN;
181     crop->var_values[VAR_Y]     = NAN;
182     crop->var_values[VAR_OUT_W] = crop->var_values[VAR_OW] = NAN;
183     crop->var_values[VAR_OUT_H] = crop->var_values[VAR_OH] = NAN;
184     crop->var_values[VAR_N]     = 0;
185     crop->var_values[VAR_T]     = NAN;
186     crop->var_values[VAR_POS]   = NAN;
187
188     av_image_fill_max_pixsteps(crop->max_step, NULL, pix_desc);
189     crop->hsub = av_pix_fmt_descriptors[link->format].log2_chroma_w;
190     crop->vsub = av_pix_fmt_descriptors[link->format].log2_chroma_h;
191
192     if ((ret = av_expr_parse_and_eval(&res, (expr = crop->ow_expr),
193                                       var_names, crop->var_values,
194                                       NULL, NULL, NULL, NULL, NULL, 0, ctx)) < 0) goto fail_expr;
195     crop->var_values[VAR_OUT_W] = crop->var_values[VAR_OW] = res;
196     if ((ret = av_expr_parse_and_eval(&res, (expr = crop->oh_expr),
197                                       var_names, crop->var_values,
198                                       NULL, NULL, NULL, NULL, NULL, 0, ctx)) < 0) goto fail_expr;
199     crop->var_values[VAR_OUT_H] = crop->var_values[VAR_OH] = res;
200     /* evaluate again ow as it may depend on oh */
201     if ((ret = av_expr_parse_and_eval(&res, (expr = crop->ow_expr),
202                                       var_names, crop->var_values,
203                                       NULL, NULL, NULL, NULL, NULL, 0, ctx)) < 0) goto fail_expr;
204     crop->var_values[VAR_OUT_W] = crop->var_values[VAR_OW] = res;
205     if (normalize_double(&crop->w, crop->var_values[VAR_OUT_W]) < 0 ||
206         normalize_double(&crop->h, crop->var_values[VAR_OUT_H]) < 0) {
207         av_log(ctx, AV_LOG_ERROR,
208                "Too big value or invalid expression for out_w/ow or out_h/oh. "
209                "Maybe the expression for out_w:'%s' or for out_h:'%s' is self-referencing.\n",
210                crop->ow_expr, crop->oh_expr);
211         return AVERROR(EINVAL);
212     }
213     crop->w &= ~((1 << crop->hsub) - 1);
214     crop->h &= ~((1 << crop->vsub) - 1);
215
216     if ((ret = av_expr_parse(&crop->x_pexpr, crop->x_expr, var_names,
217                              NULL, NULL, NULL, NULL, 0, ctx)) < 0 ||
218         (ret = av_expr_parse(&crop->y_pexpr, crop->y_expr, var_names,
219                              NULL, NULL, NULL, NULL, 0, ctx)) < 0)
220         return AVERROR(EINVAL);
221
222     av_log(ctx, AV_LOG_INFO, "w:%d h:%d -> w:%d h:%d\n",
223            link->w, link->h, crop->w, crop->h);
224
225     if (crop->w <= 0 || crop->h <= 0 ||
226         crop->w > link->w || crop->h > link->h) {
227         av_log(ctx, AV_LOG_ERROR,
228                "Invalid too big or non positive size for width '%d' or height '%d'\n",
229                crop->w, crop->h);
230         return AVERROR(EINVAL);
231     }
232
233     /* set default, required in the case the first computed value for x/y is NAN */
234     crop->x = (link->w - crop->w) / 2;
235     crop->y = (link->h - crop->h) / 2;
236     crop->x &= ~((1 << crop->hsub) - 1);
237     crop->y &= ~((1 << crop->vsub) - 1);
238     return 0;
239
240 fail_expr:
241     av_log(NULL, AV_LOG_ERROR, "Error when evaluating the expression '%s'\n", expr);
242     return ret;
243 }
244
245 static int config_output(AVFilterLink *link)
246 {
247     CropContext *crop = link->src->priv;
248
249     link->w = crop->w;
250     link->h = crop->h;
251
252     return 0;
253 }
254
255 static void start_frame(AVFilterLink *link, AVFilterBufferRef *picref)
256 {
257     AVFilterContext *ctx = link->dst;
258     CropContext *crop = ctx->priv;
259     AVFilterBufferRef *ref2;
260     int i;
261
262     ref2 = avfilter_ref_buffer(picref, ~0);
263     ref2->video->w = crop->w;
264     ref2->video->h = crop->h;
265
266     crop->var_values[VAR_T] = picref->pts == AV_NOPTS_VALUE ?
267         NAN : picref->pts * av_q2d(link->time_base);
268     crop->var_values[VAR_POS] = picref->pos == -1 ? NAN : picref->pos;
269     crop->var_values[VAR_X] = av_expr_eval(crop->x_pexpr, crop->var_values, NULL);
270     crop->var_values[VAR_Y] = av_expr_eval(crop->y_pexpr, crop->var_values, NULL);
271     crop->var_values[VAR_X] = av_expr_eval(crop->x_pexpr, crop->var_values, NULL);
272
273     normalize_double(&crop->x, crop->var_values[VAR_X]);
274     normalize_double(&crop->y, crop->var_values[VAR_Y]);
275
276     if (crop->x < 0) crop->x = 0;
277     if (crop->y < 0) crop->y = 0;
278     if ((unsigned)crop->x + (unsigned)crop->w > link->w) crop->x = link->w - crop->w;
279     if ((unsigned)crop->y + (unsigned)crop->h > link->h) crop->y = link->h - crop->h;
280     crop->x &= ~((1 << crop->hsub) - 1);
281     crop->y &= ~((1 << crop->vsub) - 1);
282
283     av_dlog(ctx, "n:%d t:%f x:%d y:%d x+w:%d y+h:%d\n",
284             (int)crop->var_values[VAR_N], crop->var_values[VAR_T], crop->x,
285             crop->y, crop->x+crop->w, crop->y+crop->h);
286
287     ref2->data[0] += crop->y * ref2->linesize[0];
288     ref2->data[0] += crop->x * crop->max_step[0];
289
290     if (!(av_pix_fmt_descriptors[link->format].flags & PIX_FMT_PAL)) {
291         for (i = 1; i < 3; i ++) {
292             if (ref2->data[i]) {
293                 ref2->data[i] += (crop->y >> crop->vsub) * ref2->linesize[i];
294                 ref2->data[i] += (crop->x * crop->max_step[i]) >> crop->hsub;
295             }
296         }
297     }
298
299     /* alpha plane */
300     if (ref2->data[3]) {
301         ref2->data[3] += crop->y * ref2->linesize[3];
302         ref2->data[3] += crop->x * crop->max_step[3];
303     }
304
305     avfilter_start_frame(link->dst->outputs[0], ref2);
306 }
307
308 static void draw_slice(AVFilterLink *link, int y, int h, int slice_dir)
309 {
310     AVFilterContext *ctx = link->dst;
311     CropContext *crop = ctx->priv;
312
313     if (y >= crop->y + crop->h || y + h <= crop->y)
314         return;
315
316     if (y < crop->y) {
317         h -= crop->y - y;
318         y  = crop->y;
319     }
320     if (y + h > crop->y + crop->h)
321         h = crop->y + crop->h - y;
322
323     avfilter_draw_slice(ctx->outputs[0], y - crop->y, h, slice_dir);
324 }
325
326 static void end_frame(AVFilterLink *link)
327 {
328     CropContext *crop = link->dst->priv;
329
330     crop->var_values[VAR_N] += 1.0;
331     avfilter_unref_buffer(link->cur_buf);
332     avfilter_end_frame(link->dst->outputs[0]);
333 }
334
335 AVFilter avfilter_vf_crop = {
336     .name      = "crop",
337     .description = NULL_IF_CONFIG_SMALL("Crop the input video to width:height:x:y."),
338
339     .priv_size = sizeof(CropContext),
340
341     .query_formats = query_formats,
342     .init          = init,
343     .uninit        = uninit,
344
345     .inputs    = (AVFilterPad[]) {{ .name             = "default",
346                                     .type             = AVMEDIA_TYPE_VIDEO,
347                                     .start_frame      = start_frame,
348                                     .draw_slice       = draw_slice,
349                                     .end_frame        = end_frame,
350                                     .get_video_buffer = avfilter_null_get_video_buffer,
351                                     .config_props     = config_input, },
352                                   { .name = NULL}},
353     .outputs   = (AVFilterPad[]) {{ .name             = "default",
354                                     .type             = AVMEDIA_TYPE_VIDEO,
355                                     .config_props     = config_output, },
356                                   { .name = NULL}},
357 };