Kendall’s Shape Coordinates

Kendall coordinates是除了Bookstein coordinates以外的另一种去除Object位置信息(对应于Euclidean Similarity Transformation中的移动)的方法。Kendall coordinates利用Helmert submatrix去除Object的Location信息。

Helmert submatrix H k − 1 × k H_{k-1 \times k} Hk1×k是去掉了第一行的Helmert matrix。Full Helmert matrix H F H^{F} HF k × k k \times k k×k的正交方阵。其第一行中每个元素的值为 1 / k 1/\sqrt{k} 1/k ,其后每行元素的定义见下文。而且,每行vector的模等于1,任意两行彼此相互正交。因此,Full Helmert matrix即相当于是一个 k k k维的直角坐标系,每行vector即是对应方向上的单位向量。数据 X k × m X_{k \times m} Xk×m经过没有第一行向量的Helmert submatix的变换,让变换后的数据 ( H X ) k − 1 × m (HX)_{k-1 \times m} (HX)k1×m不再依赖数据 X X X的原始位置信息。

注: H T H = C H^{T}H = C HTH=C C C C是对Object的Size进行测量的Centroid size方法中的centring matrix, C = I k − 1 k 1 k 1 k T C = I_{k}-\frac{1}{k}\bf{1}_{k}\bf{1}_{k}^{T} C=Ikk11k1kT

Definition 1: The j-th row of the Helmert submatrix H H H is given by: ( h j , … , h j , − j h j , 0 , … , 0 ) , h j = − j ( j + 1 ) − 1 / 2 , (h_{j}, \dots, h_{j}, -jh_{j}, 0, \dots, 0), h_{j} = -{j(j+1)}^{-1/2}, (hj,,hj,jhj,0,,0),hj=j(j+1)1/2, and so the j-th row consists of h j h_{j} hj repeated j j j times, followed by − j h j -jh_{j} jhj and then k − j − 1 k-j-1 kj1 zeros, j = 1 , … , k − 1 j = 1, \dots, k-1 j=1,,k1.

示例

对于 k = 3 k=3 k=3的full Helmert matrix:
H F = [ 1 / 3 1 / 3 1 / 3 − 1 / 2 − 1 / 2 0 − 1 / 6 − 1 / 6 2 / 6 ] H^{F} = \begin{bmatrix} 1/\sqrt{3} & 1/\sqrt{3} & 1/\sqrt{3} \\ -1/\sqrt{2} & -1/\sqrt{2} & 0 \\ -1/\sqrt{6} & -1/\sqrt{6} & 2/\sqrt{6} \end{bmatrix} HF=1/3 1/2 1/6 1/3 1/2 1/6 1/3 02/6
它的Helmert submatrix:
H = [ − 1 / 2 − 1 / 2 0 − 1 / 6 − 1 / 6 2 / 6 ] H = \begin{bmatrix} -1/\sqrt{2} & -1/\sqrt{2} & 0 \\ -1/\sqrt{6} & -1/\sqrt{6} & 2/\sqrt{6} \end{bmatrix} H=[1/2 1/6 1/2 1/6 02/6 ]
对于 k = 4 k=4 k=4的full Helmert matrix:
H F = [ 1 / 2 1 / 2 1 / 2 1 / 2 − 1 / 2 1 / 2 0 0 − 1 / 6 − 1 / 6 2 / 6 0 − 1 / 12 − 1 / 12 − 1 / 12 3 / 12 ] H^{F} = \begin{bmatrix} 1/2 & 1/2 & 1/2 & 1/2 \\ -1/\sqrt{2} & 1/\sqrt{2} & 0 & 0 \\ -1/\sqrt{6} & -1/\sqrt{6} & 2/\sqrt{6} & 0 \\ -1/\sqrt{12} & -1/\sqrt{12} & -1/\sqrt{12} & 3/\sqrt{12} \end{bmatrix} HF=1/21/2 1/6 1/12 1/21/2 1/6 1/12 1/202/6 1/12 1/2003/12
它的Helmert submatrix:
H = [ − 1 / 2 1 / 2 0 0 − 1 / 6 − 1 / 6 2 / 6 0 − 1 / 12 − 1 / 12 − 1 / 12 3 / 12 ] H = \begin{bmatrix} -1/\sqrt{2} & 1/\sqrt{2} & 0 & 0 \\ -1/\sqrt{6} & -1/\sqrt{6} & 2/\sqrt{6} & 0 \\ -1/\sqrt{12} & -1/\sqrt{12} & -1/\sqrt{12} & 3/\sqrt{12} \end{bmatrix} H=1/2 1/6 1/12 1/2 1/6 1/12 02/6 1/12 003/12

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值