What is the method to transform a Shapely polygon into a numpy pixel array that represents a mask array in Python?

edit

import numpy as np from PIL import Image, ImageDraw from shapely.geometry import Polygon def rasterize_polygon(polygon, width, height): img = Image.new('L', (width, height), 0) draw = ImageDraw.Draw(img) xy = np.array(polygon.exterior.coords.xy).T xy[:, 1] = height - xy[:, 1] xy = tuple(map(tuple, xy)) draw.polygon(xy=xy, outline=1, fill=1) return np.array(img).astype(bool) 

width = 100 height = 100 mask = np.zeros((height, width), dtype=bool) 

polygon = Polygon([(25, 25), (75, 25), (75, 75), (25, 75)]) 

polygon_mask = rasterize_polygon(polygon, width, height) mask = np.logical_or(mask, polygon_mask) 

Image.fromarray(mask.astype(np.uint8)*255).show()