Многоугольная треугольная сетка / сетка с фигурным в питоне - PullRequest
Новая
       48

Многоугольная треугольная сетка / сетка с фигурным в питоне

0 голосов
/ 30 ноября 2018

Я хочу сделать сетку из треугольников, которые являются многоугольниками, используя стройные.

Here is picture:

У меня есть список координатных точек (2 координаты для каждой точки) и список соединений. 

import numpy as np
import shapely.geometry as geometry


xlen = 20
ylen = 20
x0=0
y0=0
xPoints = np.arange(x0,xlen+1,1)
yPoints = np.arange(y0,ylen+1,1)

GridPoints = np.array([[[x,y] for x in xPoints] for y in yPoints])

triangles = [[i+j*(ylen+1),
      (i+1)+j*(ylen+1),
      i+(j+1)*(ylen+1)] for i in range(ylen) for j in range(xlen)]

Необходимы полигоны, потому что я будупозже необходимо оптимизировать эту сетку по x и y, чтобы заполнить другой многоугольник как можно большим количеством треугольников.

1 Ответ

0 голосов
/ 07 декабря 2018

Полагаю, вы хотите полностью заполнить область треугольниками.Ваши "треугольники" - это только половина из них.Если вы хотите только свою половину, просто закомментируйте вторую часть цикла for. 

import numpy as np
from shapely.geometry import Polygon
from geopandas import GeoSeries

xlen = 20
ylen = 20
x0 = 0
y0 = 0
xPoints = np.arange(x0, xlen + 1, 1)
yPoints = np.arange(y0, ylen + 1, 1)

GridPoints = list((x, y) for x in xPoints for y in yPoints)

triangles = []  # list of triangles to be populated
for i in range(ylen):
    for j in range(xlen):
        # triangles with perpendicular angle on the bottom left
        triangles.append([i + j * (ylen + 1), (i + 1) + j * (ylen + 1), i + (j + 1) * (ylen + 1)])
        # triangles with perpendicular angle on the top right
        triangles.append([(i + 1) + j * (ylen + 1), i + (j + 1) * (ylen + 1), (i + 1) + (j + 1) * (ylen + 1)])

polygons = []  # list of polygons to be populated
for triangle in triangles:
    polygon = Polygon([GridPoints[triangle[0]], GridPoints[triangle[1]], GridPoints[triangle[2]]])
    polygons.append(polygon)
gs = GeoSeries(polygons)  # save polygons to geopandas GeoSeries

Я сохраняю ваши стройные полигоны в GeoPandas GeoSeries.Результат выглядит так: GeoSeries

...