Recherche:

La ligne 8 est pas bonne.

Quand tu fais la SVD de la matrice M, la 4e colonne des vecteurs propres à droite correspond à l'équation du plan, comme indiqué dans la ligne 4.

Donc, pour obtenir Z,...

v=[v(1:2);v(4)];

pardon.

ok, une autre façon de faire...


M = [Pts ones(size(Pts,1),1)];
[u s v] = svd(M);
v=v(:,4);
# v(1)*x + v(2)*y + v(3)*z + v(4) = 0
v/=v(3);
v=[v(1:2) v(4)];

oui.

n/=n(3) c'est équivalent à n=n/n(3).

Une autre façon de faire consiste à calculer les valeurs singulières de la matrice formée par tes points, centrés en zéro.

En gros :


M=Pts - mean(Pts)
[u s v]=svd(M)
n=v(:,3) % normale à...