Comment convertir un GEOS MultiLineString en polygone en utilisant Python?

Question

Je suis en train de développer une application GeoDjango dans laquelle les utilisateurs peuvent télécharger des fichiers de cartes et effectuer des opérations de mappage de base comme l'interrogation de fonctions à l'intérieur de polygones.

J'ai reconnu que les utilisateurs peuvent télécharger des "MultiLineString" au lieu de "Polygon" parfois. Cela provoque l'échec des requêtes qui attendent des géométries fermées.

Quelle est la meilleure façon de convertir un objet MultiLineString en un Polygone en Python?

Merci.

- OMat

Answer 1

Hehe, tout d'abord je l'ai écrit ceci:

def close_geometry(self, geometry): 
    if geometry.empty or geometry[0].empty: 
     return geometry # empty 

    if(geometry[-1][-1] == geometry[0][0]): 
     return geometry # already closed 

    result = None 
    for linestring in geom: 
     if result is None: 
      resultstring = linestring.clone() 
     else: 
      resultstring.extend(linestring.coords) 

    geom = Polygon(resultstring) 

    return geom 

mais je découvre qu'il ya un peu nifty méthode appelée convex_hull qui fait la conversion de polygones automatiquement pour vous .

>>> s1 = LineString((0, 0), (1, 1), (1, 2), (0, 1)) 
>>> s1.convex_hull 
<Polygon object at ...> 
>>> s1.convex_hull.coords 
(((0.0, 0.0), (0.0, 1.0), (1.0, 2.0), (1.0, 1.0), (0.0, 0.0)),) 

>>> m1=MultiLineString(s1) 
>>> m1.convex_hull 
<Polygon object at...> 
>>> m1.convex_hull.coords 
(((0.0, 0.0), (0.0, 1.0), (1.0, 2.0), (1.0, 1.0), (0.0, 0.0)),)