Optic_flow_RF_estimation_Zhang_et_al_2022/code/OpticflowFitting/transflow_LS_optimization.m

function ssr=transflow_LS_optimization(coef,flowfield)
  flowfield_pos = flowfield(:,1);
  flowfield_vec = flowfield(:,2);
  valid_flowfield_idx = norm(flowfield_vec)>0;
   
  trans_axis = normalize(floQ(coef(2:4)));
  trans_scale = coef(1);
  fit_transField = (trans_axis+flowfield_pos.*trans_axis.*flowfield_pos)./2.*trans_scale; % Translational Field

  residue= fit_transField(valid_flowfield_idx).Qangle(flowfield_vec(valid_flowfield_idx));
  ssr = sum(residue.^2);