如何使用For循环将多边形列表显示为一张图

问题描述

所以我正在研究“视线”脚本，其中“多边形”充当建筑物，并根据线是否与“多边形”相交来返回布尔值。

它背后的逻辑有效，但是当我尝试从名为 polycoords 的列表中集成更多polygons并使用For Loop以便可以在Graph上实例化它时，它作为两个单独的图出现而不是一个。

import numpy as np
import matplotlib.pyplot as plt
import shapely.geometry
from shapely.geometry import Linestring
from shapely.geometry import Point,polygon
import descartes

origin = [-1.0,0.0] # Set a point to view from
quality = 7 # Number of grid points squared
polycoords = [[[-1,1],[-1,0.5],[0,1]],[[1,-1],[1,-0.5],[-0,-1]]] # Coordinates of the polygon
fullresults = []
newbool = []

def los (origin,quality,polycoords):
    
    for buildingpoints in range(len(polycoords)):
        x = np.linspace(-1,1,quality)
        y = np.linspace(-1,quality)
        X,Y = np.meshgrid(x,y)
        clip_poly = shapely.geometry.polygon(polycoords[buildingpoints]) 
        fig = plt.figure()
        ax = fig.add_subplot(111)
        polygonbuilding = ax.add_patch(descartes.polygonPatch(clip_poly,fc='pink',alpha=0.3))
        positions = np.vstack([Y.ravel(),X.ravel()])

        for i in range(len(positions)):
            for j in range(len(positions[i])):
                plt.scatter(*positions[::-1])
                x1 = positions[0][j]
                y1 = positions[1][j]
                line = Linestring([origin,(x1,y1)])

                if line.intersects(clip_poly) == True:
                    ax.plot(*np.array(line).T,color='red',linewidth=1,solid_capstyle='round')
                else:
                    ax.plot(*np.array(line).T,color='green',solid_capstyle='round')

                boolintersect = line.intersects(clip_poly)
                listresults = origin,[x1,y1],boolintersect
                fullresults.append(listresults)
                boollist = [x[2] for x in fullresults]
                newbool.append(sum(boollist))
                
    return(fullresults,newbool)

def analysis (losresults):
        percenteq = round((100-(newbool[-1]/len(fullresults))*100))
        print (percenteq,'% of the space is PUBLICALLY visible from point',fullresults[0][0])
        print (100-percenteq,'% of the space is PRIVATE/ISNT visible from point',fullresults[0][0])
        
los(origin,polycoords)
analysis (los)
plt.show()

当我只需要一个同时存在两个多边形的结果时，结果将显示两个单独的图形。我相信这与我有关For循环的代码的结构有关，但是我仍然很新，并且不太确定如何解决此问题。

解决方法

循环中已解决的网格生成已被移出并完成一次，如下所示：

import numpy as np
import matplotlib.pyplot as plt
import shapely.geometry
from shapely.geometry import LineString
from shapely.geometry import Point,Polygon
import descartes

origin = [-1.0,0.0] # Set a point to view from
quality = 7 # Number of grid points squared
polycoords = [[[-1,1],[-1,0.5],[0,1]],[[1,-1],[1,-0.5],[-0,-1]]] # Coordinates of the Polygon
fullresults = []
newbool = []

def init (origin,quality):
    global ax,positions
    x = np.linspace(-1,1,quality)
    y = np.linspace(-1,quality)
    X,Y = np.meshgrid(x,y)
    fig = plt.figure()
    ax = fig.add_subplot(111)
    positions = np.vstack([Y.ravel(),X.ravel()])
        
def los (origin,quality,polycoords):
    global ax,positions
    for buildingpoints in range(len(polycoords)):
        clip_poly = shapely.geometry.Polygon(polycoords[buildingpoints]) 
        polygonbuilding = ax.add_patch(descartes.PolygonPatch(clip_poly,fc='pink',alpha=0.3))

        for i in range(len(positions)):
            for j in range(len(positions[i])):
                plt.scatter(*positions[::-1])
                x1 = positions[0][j]
                y1 = positions[1][j]
                line = LineString([origin,(x1,y1)])

                if line.intersects(clip_poly) == True:
                    ax.plot(*np.array(line).T,color='red',linewidth=1,solid_capstyle='round')
                else:
                    ax.plot(*np.array(line).T,color='green',solid_capstyle='round')

                boolintersect = line.intersects(clip_poly)
                listresults = origin,[x1,y1],boolintersect
                fullresults.append(listresults)
                boollist = [x[2] for x in fullresults]
                newbool.append(sum(boollist))
                
    return(fullresults,newbool)

def analysis (losresults):
        percenteq = round((100-(newbool[-1]/len(fullresults))*100))
        print (percenteq,'% of the space is PUBLICALLY visible from point',fullresults[0][0])
        print (100-percenteq,'% of the space is PRIVATE/ISNT visible from point',fullresults[0][0])
        
init (origin,quality)
los(origin,polycoords)
analysis (los)
plt.show()