void
schro_encoder_global_estimation (SchroMotionEst *me)
{
SchroParams *params = me->params;
SchroMotionField *mf, *mf_orig;
int i;
SCHRO_ERROR("Global prediction is broken. Please try again later");
for(i=0;i<params->num_refs;i++) {
mf_orig = me->downsampled_mf[i][0];
mf = schro_motion_field_new (mf_orig->x_num_blocks, mf_orig->y_num_blocks);
memcpy (mf->motion_vectors, mf_orig->motion_vectors,
sizeof(SchroMotionVector)*mf->x_num_blocks*mf->y_num_blocks);
schro_motion_field_global_estimation (mf,
&me->encoder_frame->params.global_motion[i],
params->mv_precision);
if (i == 0) {
schro_motion_global_metric (mf,
me->encoder_frame->filtered_frame,
me->encoder_frame->ref_frame[0]->filtered_frame);
} else {
schro_motion_global_metric (mf, me->encoder_frame->filtered_frame,
me->encoder_frame->ref_frame[1]->filtered_frame);
}
schro_list_append (me->encoder_frame->motion_field_list, mf);
}
}
int
test_full_field (int width, int height, double *a, double *b, int r, int hole)
{
SchroMotionField *mf;
SchroMotionVector *mv;
SchroGlobalMotion gm;
double mult;
int i,j;
mf = schro_motion_field_new (100, 80);
printf("test_full_field: r=%d hole=%d\n"
"[%6.4f %6.4f %6.4f %6.4f] [%g %g]\n",
r, hole,
a[0], a[1], a[2], a[3], b[0], b[1]);
for(j=0;j<mf->y_num_blocks;j++){
for(i=0;i<mf->x_num_blocks;i++){
mv = mf->motion_vectors + j*mf->x_num_blocks + i;
mv->dx[0] = rint((a[0]-1)*8*i + a[1]*8*j + b[0] + r * rand_f64());
mv->dy[0] = rint(a[2]*8*i + (a[3]-1)*8*j + b[1] + r * rand_f64());
if (hole && abs(mf->y_num_blocks/2 - j) < 10 &&
abs(mf->x_num_blocks/2 - i) < 10) {
mv->dx[0] = 0;
mv->dy[0] = 0;
}
}
}
schro_motion_field_global_estimation (mf, &gm, 0);
mult = (1<<gm.a_exp);
printf("[%6.4f %6.4f %6.4f %6.4f] [%d %d]\n", gm.a00/mult, gm.a01/mult,
gm.a10/mult, gm.a11/mult, gm.b0, gm.b1);
if (fabs(gm.a00/mult - a[0]) > 0.01) return 0;
if (fabs(gm.a01/mult - a[1]) > 0.01) return 0;
if (fabs(gm.a10/mult - a[2]) > 0.01) return 0;
if (fabs(gm.a11/mult - a[3]) > 0.01) return 0;
if (fabs(gm.b0 - b[0]) > 1) return 0;
if (fabs(gm.b1 - b[1]) > 1) return 0;
return 1;
}