Discussion :

[romrom45]

Bonjour à tous,

me voilà depuis un petit moment face à un sérieux problème !
Je voudrais utiliser la méthode de ransac, pour exclure automatiquement des valeurs aberrantes lors de l'acquisition de données.

Cependant, avant d'en arriver là, je bute sur la génération d'un nuage de points suivant une fonction simple (de la forme ax+b) mais possèdant des outliers, c'est à dire des points qui ne suivent pas la fonction (des valeurs aberrantes en gros).

Je sais donc que je dois affecter à mes coordonnées des valeurs de la forme :

m_px[i] = i*(x_max-x_min)/N + x_min //ou N est le nombre total de points

m_py[i] = A*m_px[i] + B + "outliers"

De plus il me faut 3 paramètres:
g1 qui est la distance minimum d'un outlier à la courbe
g2 qui est la distance maximum d'un outlier à la courbe
P qui va etre le pourcentage d'outliers.

Voila, j'espere être assez clair vu que c'est encore assez confus pour moi même..
Dans l'attente d'une aide de votre part !

[therwald]

Petites précisions: que veux tu?

• Générer un/des jeux de test?
• Implémenter l'algo de résolution?

[romrom45]

les deux en fait !
Mais pour l'instant ce que je cherche à faire c'est tester mon programme, en insérant des valeurs aléatoires abérantes. c'est cette partie qui me chiffone, je n'arrive pas à générer mes coordonnées de points.

[Ekleog]

Code :

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

1
2
3
double rand_double(double min, double max) {
    return static_cast<double>(rand()) / RAND_MAX * (max - min) + min;
}

D'où : (avec i le point que tu veux "aberriser", [min|max]_[x|y] les bornes)

Code :

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

1
2
3
 
m_px[i] = rand_double(min_x, max_x);
m_py[i] = rand_double(min_y, max_y);

C'est pour mettre un point aléatoire. Est-ce que tu voulais bien ça ?

[therwald]

Ca lui génèrera sa population de base, mais je pense qu'il manque un terme supplémentaire qui ajouterait le terme (également aléatoire) 'outlier' tous les 100/p points de manière à créer les points aberrants. Quelque chose comme:

Code :

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

1
2
3
4
5
6
7
 
yBase = rand_double(min_y, max_y);
if(rand()<RAND_MAX*p){
m_py[i] = yBase+rand_double(yBase+g1,yBase+g2);
}else{
m_py[i] = yBase;
}

PS: romrom45 j'espère que ceci était une situation professionnelle et non un TP, sinon tu viens de te priver de l'essentiel du bénéfice de l'exercice (la réflexion personnelle, qui permet d'apprendre vraiement...)