Discussion :

[isangoma]

Bonjour,

Desole pour le titre, je n'arrivais pas a en trouver un qui decrive bien mon probleme.

J'ai des courbes qui en principe ressemblent aux deux courbes ci-dessous :



La droite est la moyenne des points de la courbe.
Mon but est de trouver les ordonees des points p1 et p2 pour lesquels "tous" les points entre p1 et p2 se trouvent au dessus de la moyenne.
Dans le cas ideal pas de probleme. Mais parfois, mes courbes sont bruitees comme les deux courbes ci-dessous :



Dans le premier cas, je veux "oublier" la partie centrale de la courbe qui repasse sous la moyenne.
Dans le second cas, ne pas considerer les pics a droite et a gauche.

Cela fait un moment que je cherche, mais je ne trouve pas d'algorithme simple pour regler ce probleme qui pourtant ne me parait pas bien complexe.

Merci par avance de votre aide.

[edit]
J'ai clarifie un peu le probleme et mis a jours les figures.

[Flodelarab]

Bonjour,

Code :

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

1
2
3
Balaie de gauche à droite. 
Si le sens de variation est toujours le même, c'est que tu n'as pas un extrémum. 
Sinon, tu notes la valeur de l'ordonnée.

En quoi cette méthode ne résout pas ton problème?

[isangoma]

Merci de la reponse.
J'avais mal pose ma problematique. Je ne cherche pas les maximum, mais les points p1 et p2 (que je viens d'ajouter).
Desole, mais avec les drogues anti-douleur, mon cerveau a un peu du mal ces jours-ci.

[lautrec1]

La recherche ne doit pas porter sur des extrema, mais sur les deux plus grands maximums locaux au-dessus de la moyenne.

Code :

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

1
2
3
4
5
6
7
8
9
10
maximum1=maximum2=0
en balayant de gauche à droite
si valeur_preprécédente <= valeur_précédente et valeur_courante <= valeur_précédente et valeur_précédente > moyenne 
       //la courbe est un maximum local au-dessus de la moyenne
       si valeur_courante > min(maximum1,maximum2) //valeur_courante doit remplacer le minimum des deux
            si maximum1 < maximum2
                  maximum1=valeur_précédente
            sinon
                  maximum2=valeur_précédente

EDIT: modification de commentaire

[isangoma]

Merci,

Je dois vraiment etre a l'ouest dans mes explications.

Mon but est de trouver les points p1 et p2.
Tous les points compris entre p1 et p2 doivent etre au dessus de la moyenne.
Il peut arriver que des faux couples de p1 et p2 existent comme sur la derniere figure. C'est du bruit et ils doivent etre ignores.
De meme, du bruit peut faire passer un ou deux points entre p1 et p2 en dessous de la moyenne et il faut ignorer ca.

J'espere etre plus clair.

[lautrec1]

Je pars de hypothèses suivantes:
- on ne doit pas tenir compte des intersections de moyenne "autour" d'un minimum de valeur supérieure à bandebruit
- p1 est une intersection "montante", p2 est une intersection non montante
- on recherche les p1 et p2, intersections consécutives montante puis descendante qui maximisent p2-p1

Code :

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

1
2
3
 
fonction EST_MINIMUM (val_cour, val_pre, val_pre_pre)
  retourne (val_pre_pre >= val_pre ET val_cour >= val_pre

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
19
20
Première partie :identifie les p1 et p2 potentiels (élimine implicitement les intersections entourant un minimum de valeur > bandebruit

montant=false
i=0
flag=0
val_cour=val_pre=val_pre_pre=0

Tant que (lire(val_cour) != fin de fichier)
     si (val_cour == moyenne)
         si flag==0 // l'intersection ne succède pas à un minimum à oublier
              pile_intersections[i++]=val_cour
              pile_sens[i++]=(val_cour >=val_pre)
         sinon
             flag=0 //l'intersection succède à minimum à oublier, on ne tient donc pas compte de l'intersection en cours
      si (EST_MINIMUM(val_cour, val_pre, val_pre_pre ET val_cour > bandebruit)
             décrémenter i //on oublie l'intersection précédente
             flag=1
      val_pre_pre=val_pre
      val_pre=val_cour

Code :

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

1
2
3
4
5
6
7
8
9
     
Deuxième partie : identifier p1 montant et p2 descendant qui maximiset p2-p1
p1=p2=max=0

pour j=0; j< i-1; j++
   si ( pile_intersections[j+1] - pile_intersections[j] > max ET pile_sens[j+1]==false ET pile_sens[j] ==true)
      max = pile_intersections[j+1] - pile_intersections[j]
      p1=pile_intersections[j]
      p2=pile_intersections[j+1]

Je pense que cela correspond mieux. Au plaisir!
Bien à toi.

EDIT : pile_sens au lieu de pile_intersections dans la seconde partie de la ligne 6 (dernier block code)

[isangoma]

Merci pour ton aide.

Je teste ca et je te dis si cela me donne de bon resultat.
Mais cela me semble prometteur.

[lautrec1]

Merci pour ton aide.

Je teste ca et je te dis si cela me donne de bon resultat.
Mais cela me semble prometteur.

A la ligne 15 du 2è bloc, il faut ajouter la condition val_cour < moyenne parce que sinon, on oublie les intersections montantes suivant un minimum local > moyenne, qui suit un maximum local.

Or ne doit oublier un minimum local que si bandebruit<minimum local<moyenne.

[isangoma]

Les premiers resulat semble bon.
Mais je pourrais confirmer une fois que j'aurais trouver pourquoi p2 semble toujours un ou deux crans trop a gauche et ainsi faire tous les tests possibles.

[lautrec1]

Les premiers resulat semble bon.
Mais je pourrais confirmer une fois que j'aurais trouver pourquoi p2 semble toujours un ou deux crans trop a gauche et ainsi faire tous les tests possibles.

Il y a une énorme erreur ici:

Code :

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

1
2
3
            pile_intersections[i++]=val_cour
            pile_sens[i++]=(val_cour >=val_pre

-->

Code :

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

1
2
3
4
 
            pile_intersections[i]=val_cour
            pile_sens[i]=(val_cour >=val_pre)
            incrémenter i

[isangoma]

Bon alors les tests sont plus ou moins fini.

Dans 95% des images (les courbes viennent d'image en niveau de gris) cela fonctionne.
Mais sur quelques images, le bruit est trop important et passe pour du signal legitime.

Je vais donc continuer a reflechir a mon probleme.

Merci pour l'aide en tous cas.