Line | |
---|
1 | function [u2,T] = pointnormiso(u); |
---|
2 | % pointnormiso Isotropic point normalization |
---|
3 | % |
---|
4 | % [u2,T] = pointnormiso(u); |
---|
5 | % u ... 3xN input data homogenous |
---|
6 | % |
---|
7 | % u2 ... 3xN normalized data homogenous |
---|
8 | % T ... 3x3 transformation matrix the does the tranformation |
---|
9 | % |
---|
10 | % $Id: pointnormiso.m,v 1.1 2005/05/23 16:16:00 svoboda Exp $ |
---|
11 | |
---|
12 | n=size(u,2); |
---|
13 | |
---|
14 | xmean = mean(u(1,:)); |
---|
15 | ymean = mean(u(2,:)); |
---|
16 | |
---|
17 | u2 = u; |
---|
18 | u2(1:2,:) = u(1:2,:) - repmat([xmean;ymean],1,n); |
---|
19 | |
---|
20 | scale = sqrt(2)/mean(sqrt(sum(u2(1:2,:).^2))); |
---|
21 | |
---|
22 | u2(1:2,:) = scale*u2(1:2,:); |
---|
23 | |
---|
24 | T = diag([scale,scale,1]); |
---|
25 | T(1,3) = -scale*xmean; |
---|
26 | T(2,3) = -scale*ymean; |
---|
27 | |
---|
28 | return; |
---|
Note: See
TracBrowser
for help on using the repository browser.