Discussion :

[jamir89]

Bonjour à tous,

J'ai un soucis sur Matlab, j'ai deux droites non parallèles (en bleu dans l'image) et je cherche à approximer une droite qui passerait au centre des deux première. Au sens des moindres carrés.
J'ai choisi de mettre les abscisses et ordonnées des deux droites dans le même vecteur pour utiliser polyfit mais ça ne donne pas ce que je souhaite (droite verte)...



Avez vous une idée pour que ça fonctionne? Même sans utiliser polyfit...

Mon code:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
p=polyfit(B,A,1); %B contient les ordonnées des deux droites, celle de gauche puis celle de droite 
%A contient les abscisses de la même manière (sans être classés par ordre croissant)
 
c(1,:)=A;c(2,:)=(p(2)+B*p(1));
plot(A,B, 'd', c(2,:), c(1,:));

[Gooby]

Bonjour,

J'entrevois une solution plus géométrique avec des équations de droites, coefficients directeurs, orthogonalité. Où bien si tu travailles sur des matrices, c'est à dire que tu peux donner une représentation matricielle de tes droites, tu peux jeter un oeil sur la fonction bwdist qui pourrait t'être fort utile je pense.

[jamir89]

Bonjour,

Je vais essayer géométriquement en utilisant l'angle entre les droites. L'approximation ne sera pas aux moindres carrés mais pour un premier test ça pourrait aller.

Si quelqu'un a une idée pour améliorer le résultat je suis preneur

Merci Gooby

[Jerome Briot]

Serait-il possible de nous fournir un jeu de données dans un fichier mat (mis dans une archive zip ou rar) ?

[jamir89]

Voila mon jeu de données:

donnees.zip

[Alex3434]

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
load('droite.mat')
load('gauche.mat')
 
y_droit = pixels_droit(:,2);
x_droit = pixels_droit(:,1);
 
y_gauche = pixels_gauche(:,2);
x_gauche = pixels_gauche(:,1);
 
milieu_y_droit = min(y_droit(:))+(max(y_droit(:)) - min(y_droit(:)))/2;
milieu_x_droit = min(x_droit(:))+(max(x_droit(:)) - min(x_droit(:)))/2;
 
milieu_y_gauche = min(y_gauche(:))+(max(y_gauche(:)) - min(y_gauche(:)))/2;
milieu_x_gauche = min(x_gauche(:))+(max(x_gauche(:)) - min(x_gauche(:)))/2;
 
figure
plot(x_droit,y_droit,'g',x_gauche,y_gauche,'r')
line([milieu_x_droit milieu_x_gauche],[milieu_x_gauche milieu_y_gauche]);

Est ce que ça correspond au résultat que tu recherches ?