Le portail francophone de la géomatique

Bonjour,

La topologie dans un SIG peut-être vue sous l'angle Arc, Face, Noeud. (pour l'aspect vectoriel).
Une face est composée d'arcs (avec le noeud final de l'arc final qui le noeud initial du premier arc)
Un arc est composé de deux noeuds (initial et final) et éventuellement n sommets.
Un noeud est le point d'intersection d'au moins trois arcs.
En utilisant des fonctions de type EstaGaucheDe, Intersection  il est facile d'implanter des algorithmes pour déterminer par exemple quel est le nombre de "polygone" traversés par une "polyligne". De plus cette méthode de stockage permet d'accélérer l'affichage et de minimiser la taille des données.

Concernant les calculs sur un raster, ceci correspond aux études de vectorisation automatique en cours actuellement.

Christophe

Bonjour,

Je ne suis pas un spécialiste, je ne sais donc pas si j'emploie la méthode nommé "buffer" ou autre. Je vous livre mon sentiment au regard des solutions que j'ai adoptées ou étudiées pour ma propre application.

Pour l'instant l'application gére les éléments vectoriels de base (pour les rasters un rectangle):
point (vertex)
Ligne
Polyligne
Polygone (polyligne fermée)
Le rectangle (windows)

Ces objets possèdent des méthodes permettant des opérations sur eux même et entre eux. Ce à l'aide de la géomètrie analytique.

L'élément de base en matière de surface et d'affichage est le rectangle:
Intersection de deux rectangles = un rectangle
Union idem

L'élément de base en "linéaire" : le point
Function
Point est dans Rectangle

Mis à part le point chaque objet est contenu dans un rectangle encadrant. Donc tous les algorithmes cherchant à calculer le résultat de l'intersection de deux objets ou savoir si un objet est contenu dans un autre pourront être optimisés.

Exemple : Point est dans Polygone ?
Si P(x,y) est dans Polygone.RectangleEncadrant Alors
Cherche une solution
Si Non
Solution = Nothing
Fin Si

Cherche une Solution :
Calcul l'intersection de Y=y et de polygone :
Obtient une liste ordonnée des segments(orienté) de polyligne ayant une intersection avec y de P.
Observe la parité de "P est à gauche de segment", les segment étant parcouru selon les X croissant.

Un exemple VB de EstaGaucheDe:
C'est une fonction de l'objet Ligne.

Public Function VertexEstAGauche(v As Vertex) As Boolean
Dim m1 As Double
Dim k1 As Double
Dim Y As Double

If Not mVerticale Then
m1 = (mb.Y - ma.Y) / (mb.X - ma.X)
k1 = mb.Y - m1 * mb.X
Y = v.Y - m1 * v.X - k1
    
Else

If ma.Y <= mb.Y Then
Y = v.X - ma.X
Else
Y = -v.X + ma.X
End If

If Y < 0 Then
VertexEstAGauche = True
Else
VertexEstAGauche = False
End If

End If
End Function

Lorsque l'on passe à une structuration topologique type arc face noeud l'intérêt est de structurer ces données dans une BD.
Par exemple:
Un arc contiendra 2 champs contenant respectivement le numéro de la parcelle à gauche, le numéro de la parcelle à droite.

Donc l'algo pour calculer les parcelles traversées par un polyligne utilisera la combinaison de la géomètrie et des relations entre objets dans la BD.

Polyligne PL,

Dans quel polygone est PL.sommet(0) ? 
         On effectue une requete: Quelles Polylignes ont un rectangle encadrant contenant PL.sommet(0)
          On utilise ensuite l'algo géomètrique sur le résultat (réduit) de la requête.

     Tant que PL.sommet(i&) => PL.sommet.count
     Tant que PL.sommet(i&) appartient au polygone résultat
     i&=i&+1
     Boucle

   Calcul quel arc intersecte PL.sommet(i&-1), PL.sommet(i&)
   En fonction de l'arc obtient le polygone voisin (si PL.sommet(i&-1) est à droite le polygone voisin est à gauche) par une requête

Boucle avec polygone voisin

A+

Christophe