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 ?

[jamir89]

Salut Alex3434

Non ce n'est pas ça. Plutot la droite rouge dans l'image ci dessous. (Je l'ai fait sur paint je n'y arrive toujours pas avec matlab)

[Gooby]

Il y a plusieurs solution. Tu peux utiliser la fonction bwdist pour chercher les points les plus distants de tes deux droites, situés entre ces deux droites. (c'est la méthode qui me semble la plus facile à mettre en place ainsi que la plus fiable)

Tu peux également calculer l'angle que forme ces deux droites et calculer sa bissectrice.

Tu peux aussi prendre le segments que forment les deux extrémités inférieures de tes droites, et calculer la médiatrice. Cela va marcher dans ton cas mais attention aux autres cas.

Il y a probablement une multitude d'autres méthodes, mais ces trois pistes devraient déjà te donner de quoi chercher si tu n'as pas d'autres idées.

[le fab]

salut

ce que tu faisais étais pas trop mal

Code :

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

1
2
3
4
5
6
7
8
9
10
figure
plot(pixels_droit(:,1),pixels_droit(:,2))
hold on
plot(pixels_gauche(:,1),pixels_gauche(:,2))
tout_pixel = [pixels_droit;pixels_gauche];
p = polyfit(tout_pixel(:,2),tout_pixel(:,1),1)
solution = polyval(p,tout_pixel(:,2))
solution = [solution tout_pixel(:,2)];
solution =sortrows(solution,1)
plot(solution(:,1),solution(:,2))

Fabien