Rev | Line | |
---|
[37] | 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.