Discussion :

[Invité]

Salut, j'essaye de calculer le signe de l'angle entre 2 vecteurs en 3D mais je rencontre quelques difficultés.

Voici une image qui montre ce que je veux faire :



Je voudrais le signe des 2 angles. (l'un est négatif, l'autre est positif)

J'ai essayer cet algorithme mais ça ne fonctionne pas, mes 2 angles sont de même signe :

Code C++ :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
 
float Vec3f::getAngleBetween (const Vec3f &other) {
    if(isNulVector() || other.isNulVector())
        return 0;
    float dotProduct = dot(other);
    Vec3f n3 = normalize().cross(other.normalize());
    Vec3f n1 = normalize();
    Vec3f n2 = other.normalize();
    Matrix4f matrix (n1.x, n1.y, n1.z, 0, n2.x, n2.y, n2.z, 0, n3.x, n3.y, n3.z, 0, 0, 0, 0, 1);
    float fDet = matrix.getDet();
    if (fDet >= 0)
        return Math::acosinus(dotProduct);
    else
        return -Math::acosinus(dotProduct);
}

Code C++ :

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

1
2
3
4
5
 
Vec3f v1(-50, 33.3333, 0);
Vec3f v2(50, 33.3333, 0);
Vec3f v3 (0, 33.3333, 0);
std::cout<<v.getAngleBetween(v1)<<" "<<v.getAngleBetween(v2)<<std::endl;

Voilà donc si quelqu'un pourrait m'aider ça serait sympa.

[Invité]

J'ai trouvé une solution mais il me faut définir les normales manuellement comme ceci :

Code c++ :

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

1
2
3
4
5
6
7
8
9
float Vec3f::getAngleBetween (const Vec3f &other, const Vec3f &n) {
    if(isNulVector() || other.isNulVector())
        return 0;
    float cosinus = dot(other);
    Vec3f crs = cross(other);
    if (n.dot(crs) < 0)
        return -Math::acosinus(cosinus);
    return Math::acosinus(cosinus);
}

Code c++ :

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

1
2
3
4
Vec3f v1(-50, 33.3333, 0);
Vec3f v2(50, 33.3333, 0);
Vec3f v3 (0, 33.3333, 0);
std::cout<<v.getAngleBetween(v1, v.cross(v1))<<" "<<v.getAngleBetween(v2, v2.cross(v)<<std::endl;

[FR119492]

Salut!

calculer le signe de l'angle entre 2 vecteurs en 3D

L'angle entre 2 vecteur n'a pas de signe: si on le regarde d'un côté, il est positif et vu de l'autre il est négatif.
Jean-Marc Blanc

[bertry]

L'angle entre 2 vecteur n'a pas de signe: si on le regarde d'un côté, il est positif et vu de l'autre il est négatif.
Jean-Marc Blanc

En effet, l'angle entre deux vecteurs n'a un signe qu'en 2D, car dans ce cas, le point de vue est imposé. Dès que tu passe à un nombre de dimensions supérieur à 2, il n'est pas possible d'y donner "spontannément" un signe.

[plegat]

Salut,

J'ai trouvé une solution mais il me faut définir les normales manuellement comme ceci :

Plutôt que de les définir manuellement, définis les automatiquement.
Tu prends ton triangle ABC, tu définis sa normale (genre AB*AC en produit vectoriel, et tu normes). C'est elle qui va "définir" le sens positif de tes angles, pour tous les vecteurs dans le plan de ton triangle (en gros, on se ramène en 2D).

Tu calcules ensuite ton produit vectoriel vecteur 0 * vecteur 1, tu fais le produit scalaire avec la normale, si c'est positif => angle positif, si c'est négatif => angle négatif.

[Invité]

Ouais erreur de ma part désolé, j'ai pris la normale du triangle (qui est v1.cross(v2)) et j'ai calculé le signe en fonction de cette normale et maintenant ça marche.

Code c++ :

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

1
2
3
4
Vec3f v1(-50, 33.3333, 0);
Vec3f v2(50, 33.3333, 0);
Vec3f v3 (0, 33.3333, 0);
std::cout<<v.getAngleBetween(v1, v1.cross(v2))<<" "<<v.getAngleBetween(v2, v1.cross(v2)<<std::endl;

Code c++ :

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

1
2
3
4
5
6
7
8
9
10
11
float Vec3f::getAngleBetween (const Vec3f &other, const Vec3f &n) {
    if(isNulVector() || other.isNulVector())
        return 0;
    float cosinus = dot(other);
    if (cosinus == 1)
        return 0;
    Vec3f crs = cross(other);
    if (n.dot(crs) < 0)
        return Math::acosinus(cosinus);
    return -Math::acosinus(cosinus);
}

Sujet résolu.