Bonjour,
Je reviens toujours sur les points d'inflexion.
Alors dans mon dernier message je me suis trompé et dire que c'est du code or ce n'est qu'un schéma d'algo , désolé
sinon je reviens sur quelque point dans votre explication:
vous avez parler des angles, arctangentes et atan2 !! or ds votre schéma d'algo y'a que la fonction angle !!!
donc veuillez svp m'expliquer pourquoi vs avez parler de ces 2 fonctions arctangente et atan2 sans que vous les utiliser ?
Quand à moi je possède dans mon package des ces trois fonctions.
donc j'ai réussi à coder la première partie de ce schéma d'algo et j'a obtenu un tableau des angles
mais je suis bloqué dans la partie suivanteoù toujours elle me rend un vecteur vide"NULL")
si angle <= seuil_bas /* On est plutôt vertical */
si seuil = seuil_haut
pt inflexion
seuil = seuil_bas
fin si
sinon
si angle >= seuil_haut /* On est plutôt horizontal */
si seuil = seuil_bas
pt d'inflexion
seuil = seuil_haut
fin si
fin si
veuillez me pardonner si j'ai pas bien compris ceci mais comme je vous ais dit: je possède dans mon package de ces trois fonctions angle, arctangente et atan2.
Aussi ci c'est le cas et je dois utiliser que la fonction angle dans mon code, veuillez m'expliquer un peu cette deuxième partie. Elle me rend toujours un vecteur NULL.
Concernant la premiere partie dont j'ai codé, voici un échantillon des résultats
que j'obtient en utilisant
angle = ( angle(i, i-1) + angle(i+1,i) ) / 2
je stocke tous ces angles dans vecteur :
1 2 3 4 5 6 7 8 9 10 11 12 13
| [1] 3.929276 5.497741 3.925795 7.065522 7.067826 5.497905 3.927240 5.497781
[9] 3.926915 5.497628 3.927009 5.497796 3.927011 7.068561 7.068466 5.497557
[17] 2.356247 3.927081 5.497851 3.927079 5.497697 3.926977 7.067894 5.496508
[25] 3.926993 5.497721 2.356219 2.356196 3.926969 5.497745 2.356210 3.926997
[33] 5.497790 3.926994 7.068566 7.068557 5.497769 3.926940 7.068376 7.068341
[41] 7.068516 7.068528 5.497690 2.356202 3.926990 7.068540 5.497721 3.926966
[49] 7.068511 5.497783 3.926993 5.497785 3.926972 7.068541 5.497777 3.926989
[57] 7.068571 5.497819 2.356269 2.356199 3.926993 7.068544 5.497759 2.356239
[65] 3.926948 5.497702 3.926971 7.068529 5.497758 2.356204 3.927002 5.497782
[73] 3.927005 5.497810 2.356317 3.927053 7.068541 7.067581 5.495789 2.356209
[81] 3.927018 7.068570 5.497783 2.356234 2.356197 3.926995 5.497822 3.927055
[89] 5.497815 3.927041 5.497734 3.926986 5.497778 3.926961 5.497728 3.926980
[97] 5.497763 3.926989 5.497779 3.926964 5.497741 3.927004 7.068517 5.497716 |
et lorsque j'utilise ce code
angle = ( atan2(i, i-1) + atan2(i+1,i) ) / 2
je stocke aussi les résultats dans un vecteur dont je vous montre un échantillon
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
| [1] 1.0861275 1.1870393 1.1662546 1.0509983 1.1592407 1.2996982 1.1481703
[8] 1.1350540 1.1662642 0.9189576 1.1661427 0.9231590 1.1559464 1.3830328
[15] 1.1764127 1.0982413 1.1642893 1.4955434 1.1525557 0.9773623 1.1751766
[22] 1.2986177 1.1749553 1.0018801 1.1724437 1.4391354 1.1758183 0.6563518
[29] 1.1746728 0.3474100 1.1709308 0.9962978 1.1758976 1.9448388 1.1760204
[36] 1.0340820 1.1729605 0.9247127 1.1655750 1.2949155 1.1713335 1.0242407
[43] 1.1773778 1.3365095 1.1727735 1.0870901 1.1656800 1.3399351 1.1631236
[50] 1.3255915 1.1614929 0.8367130 1.1764813 1.4727675 1.1779945 1.3927943
[57] 1.1703661 0.9661680 1.1626525 0.7662359 1.1764523 2.1799955 1.1762079
[64] 1.4085499 1.1770613 0.8122627 1.1762142 1.6398084 1.1764301 0.5861467
[71] 1.1765259 0.3137165 1.1775562 0.8645072 1.1758349 1.9187007 1.1704451
[78] 1.0749685 1.1705362 1.0821037 1.1754082 0.7607101 1.1764380 0.2766385
[85] 1.1763311 1.1195985 1.1769466 2.0684964 1.1776051 1.4777608 1.1744618
[92] 1.0621598 1.1702781 1.1121991 1.1692410 1.2392988 1.1629206 0.4474696
[99] 1.1767896 0.8910010 1.1766828 1.6890174 1.1776029 0.8750174 1.1768248 |
ici la fonction atan2
The arc-tangent of two arguments atan2(y, x) returns the angle between the x-axis and the vector from the origin to (x, y), i.e., for positive arguments atan2(y, x) == atan(y/x)
j'ai voulu vous informer de ces détails pour me bien aider.
ça fait des jours que j'essaye sans arriver toujours et le temps me presse
veuillez me poser vos questions si c'est pas clair.
Merci d'avance
Cordialement B.H
Partager