Discussion :

[hksa93]

Bonjour,

j'ai un nuage de points (courbe tracée dans un repère 2D), il y a une partie de ces points sont regroupés de manière linéaire (presque droite), comment trouver les points qui appartiennent à une zone linéaire, en vue d'une régression linéaire.
en réalité car je veux prendre que ces points, pas la totalité, comment faire une régression linéaire sur une partie linéaire de points (en Matlab).

Merci

[Jean Dumoncel]

Bonjour,

tu pourrais nous montrer une image de ton nuage de points?

[hksa93]

voila dans l'image joint ce que j'ai envie de faire, c'est régression linéaire aux moindres carrées pour les points de la partie linéaire de la courbe. la partie basse et la partie haute de abscisses ne m’intéresse pas.
Comment automatiser ça, donc détecter surtout les points qui sont presque droit (linéaire).

Merci

[hksa93]

ou pour bien comprendre, je veux avoir ca

[Jean Dumoncel]

Pour le moment ce n'est pas un problème lié à matlab, donc je déplace la discussion.

[pseudocode]

Bonjour,

voila dans l'image joint ce que j'ai envie de faire, c'est régression linéaire aux moindres carrées pour les points de la partie linéaire de la courbe. la partie basse et la partie haute de abscisses ne m’intéresse pas.
Comment automatiser ça, donc détecter surtout les points qui sont presque droit (linéaire).

Et bien je dirais qu'on peut commencer pas approximer la courbe par des segments de droite (algo Douglas–Peucker), puis chercher le segment qui t'intéresse (celui qui regroupe le plus de points ?).

[Aleph69]

Bonsoir,

je ne sais pas si ça marche : faire une régression linéaire, virer le point le plus éloigné de la droite, refaire la régression linéaire, et ainsi de suite.

[hksa93]

salut,
est ce que vous avez un code pour ca?

[pseudocode]

salut,
est ce que vous avez un code pour ca?

En cherchant RANSAC et LINEAR REGRESSION, tu devrais trouver ton bonheur sur le web.

[Alexis.M]

Et une transformé de Hough ? Est-ce que ça ne permetterait pas de détecter la présence de parties linéraires ?

[hksa93]

RANSAC peut etre interessante. Je ne sais pas pour la transformée de Hough car c'est pas un segment de droite (voir les points dans la courbe de la partie lineaire)

[nihao]

Oui, RANSAC semble bien pour toi.
Un exemple de code et des explications :
http://crsouza.blogspot.com/2010/06/...nsac-in-c.html

remarque :

Cette méthode ne fonctionne que si la proportion de points appartenant à la
partie linéaire est suffisante (au moins > 50 % ce qui à l'air d'être le cas sur ton exemple)
mais plus la partie linéaire est importante plus vite tu trouves la solution (moins il y a besoin de tester des valeurs).
Sinon, une extension de cette méthode peut être bien :
http://www.loria.fr/~noury/papers/noury07modeles.pdf (l'article traite d'un autre problème mais le principe est là)


cordialement

[hksa93]

Merci bien nihao, je vais voir de plus près

[Nebulix]

Les raisons pour lesquelles ta courbe s'écarte d'une droite sont "physiques" et ne peuvent pas être devinées par un programme informatique. Les données que tu pourrais en tirer seraient inévitablement faussées.
Le pic à gauche est-il lié à un artefact expérimental, à un autre signal,...? Peux-tu justifier que les n premiers points doivent être exclus de l'analyse ?
La déviation à droite est elle liée à un autre phénomène physique, à la mauvaise soustraction d'une ligne de base,...?

[hksa93]

Tout a fait, il ça des systèmes physiques qui se comportent de cette manière. c'est un comportement de processus en 1/(fréquences).^pente
c'est un spectre de puissance, l’étude sera intéressante dans la partie centrale de linéarité, c'est pour ca qu'on prend cette partie, les autres partie sont exclues.