Comment convertir un GEOS MultiLineString à Polygon en utilisant Python?

Question

Je développe une application GeoDjango où les utilisateurs peuvent télécharger des fichiers de carte et de faire des opérations de cartographie de base comme l'interrogation des caractéristiques à l'intérieur des polygones.

Je reconnu que les utilisateurs arrivent à télécharger « MultiLineString » s au lieu de « Polygon » est parfois. Cela provoque les requêtes géométries fermées attente à l'échec.

Quelle est la meilleure façon de convertir un objet MultiLineString à un Polygon 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 pour vous automatiquement.

>>> 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)),)

Answer 2

Ce petit code peut économiser beaucoup de temps, peut-être plus tard une forme plus courte geopandas sera incorporé.

import geopandas as gpd
from shapely.geometry import Polygon, mapping

def linestring_to_polygon(fili_shps):
    gdf = gpd.read_file(fili_shps) #LINESTRING
    geom = [x for x in gdf.geometry]
    all_coords = mapping(geom[0])['coordinates']
    lats = [x[1] for x in all_coords]
    lons = [x[0] for x in all_coords]
    polyg = Polygon(zip(lons, lats))
    return gpd.GeoDataFrame(index=[0], crs=gdf.crs, geometry=[polyg])