SIngle Python script in scripts/poly2raster. Test case is
Exterior ring in blue, red dots indicate crossings on exterior ring. Green lines are inner lines, inner crossings are not represented.
Steps:
Split exterior ring. This gives a list of splits for the exterior ring (blue splits) but also crossings on the exterior ring (exterior_crossings)
Split along horizontal gridlines. For each y level inside bounding box (y = 1,2 and 3 in example), split segments between two crossings that are inside the polygon. This relies on helper split_along_gridlines
Split along vertical gridlines (same as 2)
Concatenate lists of splits (of type LineString) and build split polygons using shapely.ops.polygonize
SIngle Python script in
scripts/poly2raster. Test case isExterior ring in blue, red dots indicate crossings on exterior ring. Green lines are inner lines, inner crossings are not represented.
Steps:
exterior_crossings)ylevel inside bounding box (y = 1,2 and 3 in example), split segments between two crossings that are inside the polygon. This relies on helpersplit_along_gridlinesLineString) and build split polygons usingshapely.ops.polygonizeResult: