Discussion :

[diam's]

voila le problème :

J'ai un plan (de maison, d'immeuble, de champ, ...) sout forme de fichier JPEG

Sur ce plan, je place des points, lesquels définissent des segments de droites permettant de délimiter une surface CLOSE (pièce, appartement, parcelle, etc...), cad le premier point et le dernier point de la poylignes définie par les segments de droite sont les mêmes.

Je définis aussi les échelles (horizontale et verticale) par rapport des élements du plan dont je connais les dimensions.

Là où je sèche, c'est pour calculer la surface : la forme peut être quelconque (ceux qui ont fait de l'arpentage ou qui voient régulièrement certains plans d'architectes me comprendront, il n'y a qu'a voir l'image ci jointe).


J'ai déjà regardé les forums, notamment celui "Algo". Mais très souvent on parle de définir un tableau qu'on va "colorier", de découper en triangles, ou d'utiliser des outils mathématiques que j'ai un peu de mal à comprendre (Pour le découpage en triangle à partir d'un sommet fixe, je vois franchement pas comment faire pour l'image ci dessus). La solution la plus 'envisageable' que j'ai trouvé se trouve là :Code de j.p.mignot (en C++)

Pour mon grand malheur, en plus d'être largué en math, je ne comprends pas le C++ (je comprends en gros comment fonctionne le programme, mais je n'arrive pas à le traduire en Pascal). En fait que représente le * dans

Code :

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

float *DX

S'agit il d'une forme de pointeur ? et que représente de signe de pourcentage dans la partie suivante :

Code :

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

X[(i+1)%Npts]

je suppose qu'il ne s'agit pas d'une division, mais alors, qu'est-ce ?

Donc, quelqu'un pourrait-il m'aider à traduire cet algo en pascal ?
Ou, quelqu'un connait-il une méthode pour calculer la surface d'un élement (comme sur l'image ci-dessus) sans sortir la grosse artillerie (intégrale et autres joyeusetés )?

[sjrd]

En fait que représente le * dans

Code :

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

float *DX

S'agit il d'une forme de pointeur ?

Ce n'est pas une forme de pointeur, c'est un pointeur. DX est un pointeur sur float. En Pascal ça donne :

Code :

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

var DX : ^Single;

et que représente de signe de pourcentage dans la partie suivante :

Code :

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

X[(i+1)%Npts]

je suppose qu'il ne s'agit pas d'une division, mais alors, qu'est-ce ?

C'est stricutement équivalent au mod du Pascal.

Pour le calcul d'aire, je ne connais aucun de ces algorithmes, mais j'associerait le nom "triangularisation à un sommet fixe" à ce type de résolution :

• Tu choisis un des sommets de ta forme. Notons-le O. Les autres sommets seront, dans l'ordre, A, B, C, etc. jusque N. Tu peux tourner dans un sens ou dans l'autre, peu importe.
• Pour chaque paire d'autres sommets consécutifs : A-B, B-C, ..., M-N. On nomme cette paire X-Y.
• Tu calcules l'aire du triangle OXY par la formule de Héron :
  p = demi-périmètre = ( |OX| + |XY| + |YO| ) / 2
  aire = \/¨( p (p-|OX|) (p-|XY|) (p-|YO|) )
• Tu détermines si l'angle XOY est positif ou négatif (tu peux choisir le sens, mais par logique on prendra le sens trigonométrique).
• Si l'angle est positif, alors tu ajoutes l'aire du triangle OXY à l'aire totale de ta forme ; si l'angle est négatif, tu l'en retranches.
• Si à la fin, tu te retrouves avec une aire négative, tu prends sa valeur absolue - c'est juste que tu auras tourné dans le "mauvais" sens.

Si un connaisseur d'algo peut confirmer... Là je viens d'inventer ça sur le vif

[diam's]

Merci Sjrd pour les infos sur les conventions en C++

Par contre, pour ce qui est de la triangulation à sommet fixe, je ne vois pas comment faire. C'est effectivement simple avec une forme géométrique simple, mais avec la forme ci-dessous, c'est plus compliqué. Je m'explique :

Prenons ce dessin, tiré d'un cas réel (et oui, y'a parfois des architectes 'tordus')


Il est aisé, avec ta méthode, de calculer la surface du triangle BOC. Mais il est plus délicat de calculer celle du triangle IOJ. En effet, IOJ n'est pas 'plein' et en plus, une partie de sa surface se recoupe avec celle des triangles QOP, SOT, PON et TOU. A moins que je n'ai mal compris

Alors, j'ai pensé à déplacer le point O à l'intérieur de la forme. Mais quelque soit l'endroit où je le place, je retombe forcément sur le cas 'IOJ ' ci dessus.

N'y aurait-il pas quelqu'un qui aurait développé une fonction/un compo permettant de connaitre la surface d'une forme géométrique fermée dont un connait TOUS les sommets ?

[Pedro]

(et oui, y'a parfois des architectes 'tordus')

En tant que géomètre, je te réponds +100000

Par contre, si c'est bien la surface que tu veux, il existe une méthode pour la calculer à partir des coordonnées

[diam's]

Par contre, si c'est bien la surface que tu veux, il existe une méthode pour la calculer à partir des coordonnées

ça tombe bien, c'est ce que je cherche

par contre, cette méthode, tu l'aurais pas sous le coude par hasard ?

C'est pas un truc du style remplissage de tableau j'espère ? Par ce que, les tableaux, c'est va pas etre facile car :
- si cotation en cm, une case doit représenté 1cm sur 1 cm
- pour un batiment dont la surface projetée peut être inscrite dans un carré de 100 m x 100 m, cela fait 100 000 000 de cases. En admettant que chaque case contienne un booleen (la case est DANS la forme ou AU DEHORS), cela fait 100 000 000 d'octets (un bool est, je crois, codé sur un octet), soit environ 95 mO en mémoire...

Bon d'accord, c'est vrai qu'un bâtiment de 100/100 m, c'est pas fréquent, jusqu'à présent j'ai jamais dépassé le 120 x 50 m, mais quand même...

[LadyWasky]

Va vite voir ici si ça répond à ta question, parce que personnellement, c'est ce que j'utilise
http://www.efg2.com/Lab/Graphics/PolygonArea.htm

[Pedro]

ça tombe bien, c'est ce que je cherche

par contre, cette méthode, tu l'aurais pas sous le coude par hasard ?

Je pense que l'exemple que donne Waskol est le bon

C'est pas un truc du style remplissage de tableau j'espère ? Par ce que, les tableaux, c'est va pas etre facile car :
- si cotation en cm, une case doit représenté 1cm sur 1 cm
- pour un batiment dont la surface projetée peut être inscrite dans un carré de 100 m x 100 m, cela fait 100 000 000 de cases. En admettant que chaque case contienne un booleen (la case est DANS la forme ou AU DEHORS), cela fait 100 000 000 d'octets (un bool est, je crois, codé sur un octet), soit environ 95 mO en mémoire...

Bon d'accord, c'est vrai qu'un bâtiment de 100/100 m, c'est pas fréquent, jusqu'à présent j'ai jamais dépassé le 120 x 50 m, mais quand même...

Non rien de tout cela, ça utilise les coordonnées des extrémités Et puis c'est quoi cette manie qu'ont les architectes de travailler en cm... L'unité SI est le mètre!