Fiona访问shapefile中多个多边形的坐标

极地XYZ

我正在处理一个包含4个多边形的shapefile,我需要将它们的坐标提取为元组或numpy数组,以便以后进行操作。虽然fiona可以识别4个多边形,但使用shapely来获得面积仅能为我提供第一个多边形的值。

我想为所有多边形选择xy坐标,并将它们分别放入数组和元组中进行处理。任何帮助,将不胜感激。

import fiona
import pprint
from shapely.geometry import shape, Polygon

c = fiona.open("myShapeFile.shp", 'r')
pprint.pprint(len(c))
poly = Polygon(next(iter(c))["geometry"]["coordinates"][:])
print('Area  :: ',poly.area)

执行给我以下内容:

4

区域:: 483.6379565269977

米兰·切尔马克

尝试这个:

import fiona
import numpy as np

c = fiona.open('myShapeFile.shp')
coords = [np.array(poly['geometry']['coordinates'])
          for poly in c.values()]

coords 是一个4项列表,其中包含根据shapefile中多边形的坐标构建的numpy数组:

[array([[ 3.67739738e+05,  6.35697708e+06, -1.50000000e+03],
       [ 3.67709140e+05,  6.35875947e+06, -1.50000000e+03],
       [ 3.55633222e+05,  6.35932405e+06, -1.50000000e+03],
       [ 3.55340980e+05,  6.35724482e+06, -1.50000000e+03],
       [ 3.54012992e+05,  6.35713924e+06, -1.50000000e+03],
       [ 3.53169984e+05,  6.33306237e+06, -1.50000000e+03],
       [ 3.72936914e+05,  6.33219322e+06, -1.50000000e+03],
       [ 3.73638262e+05,  6.34799928e+06, -1.50000000e+03],
       [ 3.70536492e+05,  6.34865912e+06, -1.50000000e+03],
       [ 3.71030683e+05,  6.35672004e+06, -1.50000000e+03],
       [ 3.67739738e+05,  6.35697708e+06, -1.50000000e+03]]),
 array([[ 3.60619254e+05,  6.35486104e+06, -1.50000000e+03],
       [ 3.61702258e+05,  6.35460318e+06, -1.50000000e+03],
       [ 3.62578977e+05,  6.35362332e+06, -1.50000000e+03],
       [ 3.61805402e+05,  6.35192145e+06, -1.50000000e+03],
       [ 3.60412965e+05,  6.35114787e+06, -1.50000000e+03],
       [ 3.59072102e+05,  6.35171516e+06, -1.50000000e+03],
       [ 3.59020527e+05,  6.35279817e+06, -1.50000000e+03],
       [ 3.60309824e+05,  6.35295288e+06, -1.50000000e+03],
       [ 3.60000391e+05,  6.35408746e+06, -1.50000000e+03],
       [ 3.60619254e+05,  6.35486104e+06, -1.50000000e+03]]),
 array([[ 3.59072102e+05,  6.34439198e+06, -1.50000000e+03],
       [ 3.60928684e+05,  6.34537184e+06, -1.50000000e+03],
       [ 3.62424262e+05,  6.34480455e+06, -1.50000000e+03],
       [ 3.63868270e+05,  6.34361841e+06, -1.50000000e+03],
       [ 3.63146266e+05,  6.34155554e+06, -1.50000000e+03],
       [ 3.61341258e+05,  6.34036939e+06, -1.50000000e+03],
       [ 3.60258250e+05,  6.34109139e+06, -1.50000000e+03],
       [ 3.61650688e+05,  6.34227754e+06, -1.50000000e+03],
       [ 3.61289688e+05,  6.34356684e+06, -1.50000000e+03],
       [ 3.59897250e+05,  6.34294798e+06, -1.50000000e+03],
       [ 3.57989094e+05,  6.34346370e+06, -1.50000000e+03],
       [ 3.59072102e+05,  6.34439198e+06, -1.50000000e+03]]),
 array([[ 3.69489590e+05,  6.34294798e+06, -1.50000000e+03],
       [ 3.71397742e+05,  6.34103982e+06, -1.50000000e+03],
       [ 3.69438016e+05,  6.33928638e+06, -1.50000000e+03],
       [ 3.67581438e+05,  6.34052411e+06, -1.50000000e+03],
       [ 3.67787723e+05,  6.34222597e+06, -1.50000000e+03],
       [ 3.69489590e+05,  6.34294798e+06, -1.50000000e+03]])]

本文收集自互联网,转载请注明来源。

如有侵权,请联系 [email protected] 删除。

编辑于
0

我来说两句

0 条评论
登录 后参与评论

相关文章