getinitpospara.m

cootes提出的ASM算法的matlab实现。在cootes用于课堂使用的源代码的基础上改写。改正了原来代码中的一些错误

function [s,Theta,cx,cy]=GetInitPosPara(XInputLeft,XInputRight,XBaseLeft,XBaseRight)
% function [s,Theta,cx,cy]=GetInitPosPara(XInputLeft,XInputRight,XBaseLeft,XBaseRight)
% XInputLeft,XInputRight是训练得到的平均形状模型中的左右嘴角位置
% XBaseLeft,XBaseRight是用Boosting方法得到的图上嘴角的预测位置
% 返回使平均模型接近图上嘴位置的变换参数


%                                 [ s*cos(Theta)  s*sin(Theta)  0]
% [x* y* 1] = [x y 1]*T = [x y 1]*[-s*sin(Theta)  s*cos(Theta)  0]
%                                 [     cx            cy        1]


Input=[XInputLeft;XInputRight];
Base=[XBaseLeft;XBaseRight];

mytform = cp2tform(Input,Base,'linear conformal');
T=mytform.tdata.T;
cx=T(3,1);
cy=T(3,2);

s=norm(XBaseLeft-XBaseRight)/norm(XInputLeft-XInputRight);

Theta=asin(T(1,2)/s);
if abs(Theta)<1e-6
    Theta=0;
end