9

私は、ユーザーがマップファイルをアップロードし、ポリゴン内のフィーチャのクエリなどの基本的なマッピング操作を実行できるGeoDjangoアプリケーションを開発しています。

ユーザーがたまたま「Polygon」ではなく「MultiLineString」をアップロードすることを認識しました。これにより、閉じたジオメトリを期待するクエリが失敗します。

PythonでMultiLineStringオブジェクトをポリゴンに変換する最良の方法は何ですか?

4

3 に答える 3

12

ふふ、最初はこれを書いた:

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

しかし、その後、ポリゴン変換を自動的に行う凸包と呼ばれる気の利いた小さなメソッドがあることを発見しました。

>>> 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)),)
于 2010-06-04T14:47:44.840 に答える
2

この小さなコードは多くの時間を節約することができ、おそらく後でジオパンダのより短い形式が組み込まれるでしょう。

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])
于 2019-07-21T19:11:29.580 に答える
2

これはCarlosの回答の変更です。これはより単純で、1つの要素だけでなく、ソースファイルのすべての行を返します。

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

def linestring_to_polygon(fili_shps):
    gdf = gpd.read_file(fili_shps) #LINESTRING
    gdf['geometry'] = [Polygon(mapping(x)['coordinates']) for x in gdf.geometry]
    return gdf
于 2020-05-08T17:08:52.487 に答える