{ "spline16", "spline16 interpolation", 0, AV_OPT_TYPE_CONST, {.i64=SPLINE16}, 0, 0, FLAGS, "interp" },
{ "gauss", "gaussian interpolation", 0, AV_OPT_TYPE_CONST, {.i64=GAUSSIAN}, 0, 0, FLAGS, "interp" },
{ "gaussian", "gaussian interpolation", 0, AV_OPT_TYPE_CONST, {.i64=GAUSSIAN}, 0, 0, FLAGS, "interp" },
+ { "mitchell", "mitchell interpolation", 0, AV_OPT_TYPE_CONST, {.i64=MITCHELL}, 0, 0, FLAGS, "interp" },
{ "w", "output width", OFFSET(width), AV_OPT_TYPE_INT, {.i64=0}, 0, INT16_MAX, FLAGS, "w"},
{ "h", "output height", OFFSET(height), AV_OPT_TYPE_INT, {.i64=0}, 0, INT16_MAX, FLAGS, "h"},
{ "in_stereo", "input stereo format", OFFSET(in_stereo), AV_OPT_TYPE_INT, {.i64=STEREO_2D}, 0, NB_STEREO_FMTS-1, FLAGS, "stereo" },
{ "pitch", "pitch rotation", OFFSET(pitch), AV_OPT_TYPE_FLOAT, {.dbl=0.f}, -180.f, 180.f,TFLAGS, "pitch"},
{ "roll", "roll rotation", OFFSET(roll), AV_OPT_TYPE_FLOAT, {.dbl=0.f}, -180.f, 180.f,TFLAGS, "roll"},
{ "rorder", "rotation order", OFFSET(rorder), AV_OPT_TYPE_STRING, {.str="ypr"}, 0, 0,TFLAGS, "rorder"},
- { "h_fov", "output horizontal field of view",OFFSET(h_fov), AV_OPT_TYPE_FLOAT, {.dbl=90.f}, 0.00001f, 360.f,TFLAGS, "h_fov"},
- { "v_fov", "output vertical field of view", OFFSET(v_fov), AV_OPT_TYPE_FLOAT, {.dbl=45.f}, 0.00001f, 360.f,TFLAGS, "v_fov"},
+ { "h_fov", "output horizontal field of view",OFFSET(h_fov), AV_OPT_TYPE_FLOAT, {.dbl=0.f}, 0.f, 360.f,TFLAGS, "h_fov"},
+ { "v_fov", "output vertical field of view", OFFSET(v_fov), AV_OPT_TYPE_FLOAT, {.dbl=0.f}, 0.f, 360.f,TFLAGS, "v_fov"},
{ "d_fov", "output diagonal field of view", OFFSET(d_fov), AV_OPT_TYPE_FLOAT, {.dbl=0.f}, 0.f, 360.f,TFLAGS, "d_fov"},
{ "h_flip", "flip out video horizontally", OFFSET(h_flip), AV_OPT_TYPE_BOOL, {.i64=0}, 0, 1,TFLAGS, "h_flip"},
{ "v_flip", "flip out video vertically", OFFSET(v_flip), AV_OPT_TYPE_BOOL, {.i64=0}, 0, 1,TFLAGS, "v_flip"},
{ "iv_flip", "flip in video vertically", OFFSET(iv_flip), AV_OPT_TYPE_BOOL, {.i64=0}, 0, 1,TFLAGS, "iv_flip"},
{ "in_trans", "transpose video input", OFFSET(in_transpose), AV_OPT_TYPE_BOOL, {.i64=0}, 0, 1, FLAGS, "in_transpose"},
{ "out_trans", "transpose video output", OFFSET(out_transpose), AV_OPT_TYPE_BOOL, {.i64=0}, 0, 1, FLAGS, "out_transpose"},
- { "ih_fov", "input horizontal field of view",OFFSET(ih_fov), AV_OPT_TYPE_FLOAT, {.dbl=90.f}, 0.00001f, 360.f,TFLAGS, "ih_fov"},
- { "iv_fov", "input vertical field of view", OFFSET(iv_fov), AV_OPT_TYPE_FLOAT, {.dbl=45.f}, 0.00001f, 360.f,TFLAGS, "iv_fov"},
+ { "ih_fov", "input horizontal field of view",OFFSET(ih_fov), AV_OPT_TYPE_FLOAT, {.dbl=0.f}, 0.f, 360.f,TFLAGS, "ih_fov"},
+ { "iv_fov", "input vertical field of view", OFFSET(iv_fov), AV_OPT_TYPE_FLOAT, {.dbl=0.f}, 0.f, 360.f,TFLAGS, "iv_fov"},
{ "id_fov", "input diagonal field of view", OFFSET(id_fov), AV_OPT_TYPE_FLOAT, {.dbl=0.f}, 0.f, 360.f,TFLAGS, "id_fov"},
{"alpha_mask", "build mask in alpha plane", OFFSET(alpha), AV_OPT_TYPE_BOOL, {.i64=0}, 0, 1, FLAGS, "alpha"},
{ NULL }
case LANCZOS:
case SPLINE16:
case GAUSSIAN:
+ case MITCHELL:
s->remap_line = depth <= 8 ? remap4_8bit_line_c : remap4_16bit_line_c;
break;
}
}
}
+/**
+ * Calculate 1-dimensional cubic_bc_spline coefficients.
+ *
+ * @param t relative coordinate
+ * @param coeffs coefficients
+ */
+static void calculate_cubic_bc_coeffs(float t, float *coeffs,
+ float b, float c)
+{
+ float sum = 0.f;
+ float p0 = (6.f - 2.f * b) / 6.f,
+ p2 = (-18.f + 12.f * b + 6.f * c) / 6.f,
+ p3 = (12.f - 9.f * b - 6.f * c) / 6.f,
+ q0 = (8.f * b + 24.f * c) / 6.f,
+ q1 = (-12.f * b - 48.f * c) / 6.f,
+ q2 = (6.f * b + 30.f * c) / 6.f,
+ q3 = (-b - 6.f * c) / 6.f;
+
+ for (int i = 0; i < 4; i++) {
+ const float x = fabsf(t - i + 1.f);
+ if (x < 1.f) {
+ coeffs[i] = (p0 + x * x * (p2 + x * p3)) *
+ (p0 + x * x * (p2 + x * p3 / 2.f) / 4.f);
+ } else if (x < 2.f) {
+ coeffs[i] = (q0 + x * (q1 + x * (q2 + x * q3))) *
+ (q0 + x * (q1 + x * (q2 + x / 2.f * q3) / 2.f) / 2.f);
+ } else {
+ coeffs[i] = 0.f;
+ }
+ sum += coeffs[i];
+ }
+
+ for (int i = 0; i < 4; i++) {
+ coeffs[i] /= sum;
+ }
+}
+
+/**
+ * Calculate kernel for mitchell interpolation.
+ *
+ * @param du horizontal relative coordinate
+ * @param dv vertical relative coordinate
+ * @param rmap calculated 4x4 window
+ * @param u u remap data
+ * @param v v remap data
+ * @param ker ker remap data
+ */
+static void mitchell_kernel(float du, float dv, const XYRemap *rmap,
+ int16_t *u, int16_t *v, int16_t *ker)
+{
+ float du_coeffs[4];
+ float dv_coeffs[4];
+
+ calculate_cubic_bc_coeffs(du, du_coeffs, 1.f / 3.f, 1.f / 3.f);
+ calculate_cubic_bc_coeffs(dv, dv_coeffs, 1.f / 3.f, 1.f / 3.f);
+
+ for (int i = 0; i < 4; i++) {
+ for (int j = 0; j < 4; j++) {
+ u[i * 4 + j] = rmap->u[i][j];
+ v[i * 4 + j] = rmap->v[i][j];
+ ker[i * 4 + j] = lrintf(du_coeffs[j] * dv_coeffs[i] * 16385.f);
+ }
+ }
+}
+
/**
* Modulo operation with only positive remainders.
*
static inline int reflecty(int y, int h)
{
if (y < 0) {
- return -y;
+ y = -y;
} else if (y >= h) {
- return 2 * h - 1 - y;
+ y = 2 * h - 1 - y;
}
- return y;
+ return av_clip(y, 0, h - 1);
}
/**
return 1;
}
+/**
+ * Prepare data for processing equirectangular output format.
+ *
+ * @param ctx filter context
+ *
+ * @return error code
+ */
+static int prepare_equirect_out(AVFilterContext *ctx)
+{
+ V360Context *s = ctx->priv;
+
+ s->flat_range[0] = s->h_fov * M_PI / 360.f;
+ s->flat_range[1] = s->v_fov * M_PI / 360.f;
+
+ return 0;
+}
+
/**
* Calculate 3D coordinates on sphere for corresponding frame position in equirectangular format.
*
int i, int j, int width, int height,
float *vec)
{
- const float phi = ((2.f * i + 0.5f) / width - 1.f) * M_PI;
- const float theta = ((2.f * j + 0.5f) / height - 1.f) * M_PI_2;
+ const float phi = ((2.f * i + 0.5f) / width - 1.f) * s->flat_range[0];
+ const float theta = ((2.f * j + 0.5f) / height - 1.f) * s->flat_range[1];
const float sin_phi = sinf(phi);
const float cos_phi = cosf(phi);
return visible;
}
+/**
+ * Prepare data for processing equirectangular input format.
+ *
+ * @param ctx filter context
+ *
+ * @return error code
+ */
+static int prepare_equirect_in(AVFilterContext *ctx)
+{
+ V360Context *s = ctx->priv;
+
+ s->iflat_range[0] = s->ih_fov * M_PI / 360.f;
+ s->iflat_range[1] = s->iv_fov * M_PI / 360.f;
+
+ return 0;
+}
+
/**
* Calculate frame position in equirectangular format for corresponding 3D coordinates on sphere.
*
const float phi = atan2f(vec[0], vec[2]);
const float theta = asinf(vec[1]);
- const float uf = (phi / M_PI + 1.f) * width / 2.f;
- const float vf = (theta / M_PI_2 + 1.f) * height / 2.f;
+ const float uf = (phi / s->iflat_range[0] + 1.f) * width / 2.f;
+ const float vf = (theta / s->iflat_range[1] + 1.f) * height / 2.f;
const int ui = floorf(uf);
const int vi = floorf(vf);
+ int visible;
*du = uf - ui;
*dv = vf - vi;
+ visible = vi >= 0 && vi < height && ui >= 0 && ui < width;
+
for (int i = 0; i < 4; i++) {
for (int j = 0; j < 4; j++) {
us[i][j] = ereflectx(ui + j - 1, vi + i - 1, width, height);
}
}
- return 1;
+ return visible;
}
/**
const float cos_theta = cosf(theta);
vec[0] = cos_theta * sin_phi;
- vec[1] = sin_theta;
- vec[2] = cos_theta * cos_phi;
+ vec[1] = cos_theta * cos_phi;
+ vec[2] = sin_theta;
} else {
vec[0] = 0.f;
vec[1] = 1.f;
return 0;
}
- normalize_vector(vec);
return 1;
}
return 1;
}
-static void multiply_matrix(float c[3][3], const float a[3][3], const float b[3][3])
+static void multiply_quaternion(float c[4], const float a[4], const float b[4])
{
- for (int i = 0; i < 3; i++) {
- for (int j = 0; j < 3; j++) {
- float sum = 0.f;
-
- for (int k = 0; k < 3; k++)
- sum += a[i][k] * b[k][j];
+ c[0] = a[0] * b[0] - a[1] * b[1] - a[2] * b[2] - a[3] * b[3];
+ c[1] = a[1] * b[0] + a[0] * b[1] + a[2] * b[3] - a[3] * b[2];
+ c[2] = a[2] * b[0] + a[0] * b[2] + a[3] * b[1] - a[1] * b[3];
+ c[3] = a[3] * b[0] + a[0] * b[3] + a[1] * b[2] - a[2] * b[1];
+}
- c[i][j] = sum;
- }
- }
+static void conjugate_quaternion(float d[4], const float q[4])
+{
+ d[0] = q[0];
+ d[1] = -q[1];
+ d[2] = -q[2];
+ d[3] = -q[3];
}
/**
- * Calculate rotation matrix for yaw/pitch/roll angles.
+ * Calculate rotation quaternion for yaw/pitch/roll angles.
*/
-static inline void calculate_rotation_matrix(float yaw, float pitch, float roll,
- float rot_mat[3][3],
- const int rotation_order[3])
+static inline void calculate_rotation(float yaw, float pitch, float roll,
+ float rot_quaternion[2][4],
+ const int rotation_order[3])
{
const float yaw_rad = yaw * M_PI / 180.f;
const float pitch_rad = pitch * M_PI / 180.f;
const float roll_rad = roll * M_PI / 180.f;
- const float sin_yaw = sinf(yaw_rad);
- const float cos_yaw = cosf(yaw_rad);
- const float sin_pitch = sinf(pitch_rad);
- const float cos_pitch = cosf(pitch_rad);
- const float sin_roll = sinf(roll_rad);
- const float cos_roll = cosf(roll_rad);
+ const float sin_yaw = sinf(yaw_rad * 0.5f);
+ const float cos_yaw = cosf(yaw_rad * 0.5f);
+ const float sin_pitch = sinf(pitch_rad * 0.5f);
+ const float cos_pitch = cosf(pitch_rad * 0.5f);
+ const float sin_roll = sinf(roll_rad * 0.5f);
+ const float cos_roll = cosf(roll_rad * 0.5f);
- float m[3][3][3];
- float temp[3][3];
+ float m[3][4];
+ float tmp[2][4];
- m[0][0][0] = cos_yaw; m[0][0][1] = 0; m[0][0][2] = sin_yaw;
- m[0][1][0] = 0; m[0][1][1] = 1; m[0][1][2] = 0;
- m[0][2][0] = -sin_yaw; m[0][2][1] = 0; m[0][2][2] = cos_yaw;
+ m[0][0] = cos_yaw; m[0][1] = 0.f; m[0][2] = sin_yaw; m[0][3] = 0.f;
+ m[1][0] = cos_pitch; m[1][1] = sin_pitch; m[1][2] = 0.f; m[1][3] = 0.f;
+ m[2][0] = cos_roll; m[2][1] = 0.f; m[2][2] = 0.f; m[2][3] = sin_roll;
- m[1][0][0] = 1; m[1][0][1] = 0; m[1][0][2] = 0;
- m[1][1][0] = 0; m[1][1][1] = cos_pitch; m[1][1][2] = -sin_pitch;
- m[1][2][0] = 0; m[1][2][1] = sin_pitch; m[1][2][2] = cos_pitch;
+ multiply_quaternion(tmp[0], rot_quaternion[0], m[rotation_order[0]]);
+ multiply_quaternion(tmp[1], tmp[0], m[rotation_order[1]]);
+ multiply_quaternion(rot_quaternion[0], tmp[1], m[rotation_order[2]]);
- m[2][0][0] = cos_roll; m[2][0][1] = -sin_roll; m[2][0][2] = 0;
- m[2][1][0] = sin_roll; m[2][1][1] = cos_roll; m[2][1][2] = 0;
- m[2][2][0] = 0; m[2][2][1] = 0; m[2][2][2] = 1;
-
- multiply_matrix(temp, m[rotation_order[0]], m[rotation_order[1]]);
- multiply_matrix(rot_mat, temp, m[rotation_order[2]]);
+ conjugate_quaternion(rot_quaternion[1], rot_quaternion[0]);
}
/**
- * Rotate vector with given rotation matrix.
+ * Rotate vector with given rotation quaternion.
*
- * @param rot_mat rotation matrix
+ * @param rot_quaternion rotation quaternion
* @param vec vector
*/
-static inline void rotate(const float rot_mat[3][3],
+static inline void rotate(const float rot_quaternion[2][4],
float *vec)
{
- const float x_tmp = vec[0] * rot_mat[0][0] + vec[1] * rot_mat[0][1] + vec[2] * rot_mat[0][2];
- const float y_tmp = vec[0] * rot_mat[1][0] + vec[1] * rot_mat[1][1] + vec[2] * rot_mat[1][2];
- const float z_tmp = vec[0] * rot_mat[2][0] + vec[1] * rot_mat[2][1] + vec[2] * rot_mat[2][2];
+ float qv[4], temp[4], rqv[4];
+
+ qv[0] = 0.f;
+ qv[1] = vec[0];
+ qv[2] = vec[1];
+ qv[3] = vec[2];
- vec[0] = x_tmp;
- vec[1] = y_tmp;
- vec[2] = z_tmp;
+ multiply_quaternion(temp, rot_quaternion[0], qv);
+ multiply_quaternion(rqv, temp, rot_quaternion[1]);
+
+ vec[0] = rqv[1];
+ vec[1] = rqv[2];
+ vec[2] = rqv[3];
}
static inline void set_mirror_modifier(int h_flip, int v_flip, int d_flip,
static void fov_from_dfov(int format, float d_fov, float w, float h, float *h_fov, float *v_fov)
{
switch (format) {
+ case EQUIRECTANGULAR:
+ *h_fov = d_fov;
+ *v_fov = d_fov * 0.5f;
+ break;
case ORTHOGRAPHIC:
{
const float d = 0.5f * hypotf(w, h);
break;
case DUAL_FISHEYE:
{
- const float d = 0.5f * hypotf(w * 0.5f, h);
+ const float d = hypotf(w * 0.5f, h);
- *h_fov = d / w * 2.f * d_fov;
- *v_fov = d / h * d_fov;
+ *h_fov = 0.5f * w / d * d_fov;
+ *v_fov = h / d * d_fov;
}
break;
case FISHEYE:
{
- const float d = 0.5f * hypotf(w, h);
+ const float d = hypotf(w, h);
- *h_fov = d / w * d_fov;
- *v_fov = d / h * d_fov;
+ *h_fov = w / d * d_fov;
+ *v_fov = h / d * d_fov;
}
break;
case FLAT:
else
out_mask = s->out_transform(s, i, j, width, height, vec);
av_assert1(!isnan(vec[0]) && !isnan(vec[1]) && !isnan(vec[2]));
- rotate(s->rot_mat, vec);
+ rotate(s->rot_quaternion, vec);
av_assert1(!isnan(vec[0]) && !isnan(vec[1]) && !isnan(vec[2]));
normalize_vector(vec);
mirror(s->output_mirror_modifier, vec);
const AVPixFmtDescriptor *desc = av_pix_fmt_desc_get(inlink->format);
const int depth = desc->comp[0].depth;
const int sizeof_mask = s->mask_size = (depth + 7) >> 3;
+ float default_h_fov = 360.f;
+ float default_v_fov = 180.f;
+ float default_ih_fov = 360.f;
+ float default_iv_fov = 180.f;
int sizeof_uv;
int sizeof_ker;
int err;
sizeof_uv = sizeof(int16_t) * s->elements;
sizeof_ker = sizeof(int16_t) * s->elements;
break;
+ case MITCHELL:
+ s->calculate_kernel = mitchell_kernel;
+ s->remap_slice = depth <= 8 ? remap4_8bit_slice : remap4_16bit_slice;
+ s->elements = 4 * 4;
+ sizeof_uv = sizeof(int16_t) * s->elements;
+ sizeof_ker = sizeof(int16_t) * s->elements;
+ break;
default:
av_assert0(0);
}
s->in_width = s->inplanewidth[0];
s->in_height = s->inplaneheight[0];
+ switch (s->in) {
+ case CYLINDRICAL:
+ case FLAT:
+ default_ih_fov = 90.f;
+ default_iv_fov = 45.f;
+ break;
+ case EQUISOLID:
+ case ORTHOGRAPHIC:
+ case STEREOGRAPHIC:
+ case DUAL_FISHEYE:
+ case FISHEYE:
+ default_ih_fov = 180.f;
+ default_iv_fov = 180.f;
+ default:
+ break;
+ }
+
+ if (s->ih_fov == 0.f)
+ s->ih_fov = default_ih_fov;
+
+ if (s->iv_fov == 0.f)
+ s->iv_fov = default_iv_fov;
+
if (s->id_fov > 0.f)
fov_from_dfov(s->in, s->id_fov, w, h, &s->ih_fov, &s->iv_fov);
switch (s->in) {
case EQUIRECTANGULAR:
s->in_transform = xyz_to_equirect;
- err = 0;
+ err = prepare_equirect_in(ctx);
wf = w;
hf = h;
break;
switch (s->out) {
case EQUIRECTANGULAR:
s->out_transform = equirect_to_xyz;
- prepare_out = NULL;
+ prepare_out = prepare_equirect_out;
w = lrintf(wf);
h = lrintf(hf);
break;
s->width = w;
s->height = h;
+ switch (s->out) {
+ case CYLINDRICAL:
+ case FLAT:
+ default_h_fov = 90.f;
+ default_v_fov = 45.f;
+ break;
+ case EQUISOLID:
+ case ORTHOGRAPHIC:
+ case STEREOGRAPHIC:
+ case DUAL_FISHEYE:
+ case FISHEYE:
+ default_h_fov = 180.f;
+ default_v_fov = 180.f;
+ break;
+ default:
+ break;
+ }
+
+ if (s->h_fov == 0.f)
+ s->h_fov = default_h_fov;
+
+ if (s->v_fov == 0.f)
+ s->v_fov = default_v_fov;
+
if (s->d_fov > 0.f)
fov_from_dfov(s->out, s->d_fov, w, h, &s->h_fov, &s->v_fov);
return err;
}
- calculate_rotation_matrix(s->yaw, s->pitch, s->roll, s->rot_mat, s->rotation_order);
+ calculate_rotation(s->yaw, s->pitch, s->roll,
+ s->rot_quaternion, s->rotation_order);
+
set_mirror_modifier(s->h_flip, s->v_flip, s->d_flip, s->output_mirror_modifier);
ctx->internal->execute(ctx, v360_slice, NULL, NULL, s->nb_threads);
static int process_command(AVFilterContext *ctx, const char *cmd, const char *args,
char *res, int res_len, int flags)
{
+ V360Context *s = ctx->priv;
int ret;
+ s->yaw = s->pitch = s->roll = 0.f;
+
ret = ff_filter_process_command(ctx, cmd, args, res, res_len, flags);
if (ret < 0)
return ret;
return config_output(ctx->outputs[0]);
}
+static av_cold int init(AVFilterContext *ctx)
+{
+ V360Context *s = ctx->priv;
+
+ s->rot_quaternion[0][0] = 1.f;
+ s->rot_quaternion[0][1] = s->rot_quaternion[0][2] = s->rot_quaternion[0][3] = 0.f;
+
+ return 0;
+}
+
static av_cold void uninit(AVFilterContext *ctx)
{
V360Context *s = ctx->priv;
{ NULL }
};
-AVFilter ff_vf_v360 = {
+const AVFilter ff_vf_v360 = {
.name = "v360",
.description = NULL_IF_CONFIG_SMALL("Convert 360 projection of video."),
.priv_size = sizeof(V360Context),
+ .init = init,
.uninit = uninit,
.query_formats = query_formats,
.inputs = inputs,
projects / ffmpeg / blobdiff