# from clipper import Area, Clipper, Point, ClipType, PolyType, PolyFillType from clipper import * import math import re from random import randint # from subprocess import call import os #=============================================================================== #=============================================================================== def LoadFile1(lines): # File type 1: first line is total polygons count and subsequent lines # contain the polygon vertex count followed by its coords try: polygons = [] poly = [] for l in lines: vals = re.split(' |, |,', l.strip()) if len(vals) < 2: if (len(poly) > 2): polygons.append(poly) poly = [] else: poly.append(Point(int(vals[0]), int(vals[1]))) if (len(poly) > 2): polygons.append(poly) return polygons except: return None #=============================================================================== def LoadFile2(lines): # File type 2: vertex coords on consecutive lines for each polygon # where each polygon is separated by an empty line try: polygons = [] poly = [] for l in lines: l = l.strip() if (l == ''): if (len(poly) > 2): polygons.append(poly) poly = [] else: vals = re.split(' |, |,', l) poly.append(Point(int(vals[0]), int(vals[1]))) if (len(poly) > 2): polygons.append(poly) return polygons except: return None #=============================================================================== def LoadFile(filename): try: f = open(filename, 'r') try: lines = f.readlines() finally: f.close() # pick file type from format of first line ... if len(lines) == 0: return [] elif not ',' in lines[0]: return LoadFile1(lines) else: return LoadFile2(lines) except: return None #=============================================================================== def SaveToFile(filename, polys, scale = 1.0): invScale = 1.0 / scale try: f = open(filename, 'w') try: if invScale == 1: for poly in polys: for pt in poly: f.write("{0}, {1}\n".format(pt.x, pt.y)) f.write("\n") else: for poly in polys: for pt in poly: f.write("{0:.4f}, {1:.4f}\n".format(pt.x * invScale, pt.y * invScale)) f.write("\n") finally: f.close() except: return #=============================================================================== def RandomPoly(maxWidth, maxHeight, vertCnt): result = [] for _ in range(vertCnt): result.append(Point(randint(0, maxWidth), randint(0, maxHeight))) return result #=============================================================================== # SVGBuilder #=============================================================================== class SVGBuilder(object): def HtmlColor(self, val): return "#{0:06x}".format(val & 0xFFFFFF) def AlphaClr(self, val): return "{0:.2f}".format(float(val >> 24)/255) class StyleInfo(object): def __init__(self): self.fillType = PolyFillType.EvenOdd self.brushClr = 0 self.penClr = 0 self.penWidth = 0.8 self.showCoords = False class StyleInfoPlus(StyleInfo): def __init__(self): SVGBuilder.StyleInfo.__init__(self) self.polygons = [] self.textlines = [] def __init__(self): self.GlobalStyle = SVGBuilder.StyleInfo() self.PolyInfoList = [] self.PathHeader = " \n\n" self.Header = """ \n \n\n""" def AddPolygon(self, poly, brushColor, penColor): if poly is None or len(poly) == 0: return pi = self.StyleInfoPlus() pi.penWidth = self.GlobalStyle.penWidth pi.fillType = self.GlobalStyle.fillType pi.showCoords = self.GlobalStyle.showCoords pi.brushClr = brushColor pi.penClr = penColor pi.polygons.append(poly) self.PolyInfoList.append(pi) def AddPolygons(self, polys, brushColor, penColor): if polys is None or len(polys) == 0: return pi = self.StyleInfoPlus() pi.penWidth = self.GlobalStyle.penWidth pi.fillType = self.GlobalStyle.fillType pi.showCoords = self.GlobalStyle.showCoords pi.brushClr = brushColor pi.penClr = penColor pi.polygons = polys self.PolyInfoList.append(pi) def SaveToFile(self, filename, invScale = 1.0, margin = 10): if len(self.PolyInfoList) == 0: return False if invScale == 0: invScale = 1.0 if margin < 0: margin = 0 pi = self.PolyInfoList[0] # get bounding rect ... left = right = pi.polygons[0][0].x top = bottom = pi.polygons[0][0].y for pi in self.PolyInfoList: for p in pi.polygons: for ip in p: if ip.x < left: left = ip.x if ip.x > right: right = ip.x if ip.y < top: top = ip.y if ip.y > bottom: bottom = ip.y left *= invScale top *= invScale right *= invScale bottom *= invScale offsetX = -left + margin offsetY = -top + margin f = open(filename, 'w') m2 = margin * 2 f.write(self.Header.format(right - left + m2, bottom - top + m2)) for pi in self.PolyInfoList: f.write(self.PathHeader) for p in pi.polygons: cnt = len(p) if cnt < 3: continue f.write(" M {0:.2f} {1:.2f}".format(p[0].x * invScale + offsetX, p[0].y * invScale + offsetY)) for i in range(1,cnt): f.write(" L {0:.2f} {1:.2f}".format(p[i].x * invScale + offsetX, p[i].y * invScale + offsetY)) f.write(" z") fillRule = "evenodd" if pi.fillType != PolyFillType.EvenOdd: fillRule = "nonzero" f.write(self.PathFooter.format(self.HtmlColor(pi.brushClr), self.AlphaClr(pi.brushClr), fillRule, self.HtmlColor(pi.penClr), self.AlphaClr(pi.penClr), pi.penWidth)) if (pi.showCoords): f.write("\n\n") for p in pi.polygons: cnt = len(p) if cnt < 3: continue for pt in p: x = pt.x * invScale + offsetX y = pt.y * invScale + offsetY f.write("{2},{3}\n".format(x, y, pt.x, pt.y)) f.write("\n") f.write("\n") f.write("\n") f.close() return True #=============================================================================== # Main entry ... #=============================================================================== scaleExp = 0 scale = math.pow(10, scaleExp) invScale = 1.0 / scale subj, clip = [], [] #load saved subject and clip polygons ... #subj = LoadFile('./subj.txt') #clip = LoadFile('./clip.txt') # Generate random subject and clip polygons ... subj.append(RandomPoly(640 * scale, 480 * scale, 100)) clip.append(RandomPoly(640 * scale, 480 * scale, 100)) #SaveToFile('./subj2.txt', subj, scale) #SaveToFile('./clip2.txt', clip, scale) # Load the polygons into Clipper and execute the boolean clip op ... c = Clipper() solution = [] pft = PolyFillType.EvenOdd c.AddPolygons(subj, PolyType.Subject) c.AddPolygons(clip, PolyType.Clip) result = c.Execute(ClipType.Intersection, solution, pft, pft) SaveToFile('./solution2.txt', solution, scale) # Create an SVG file to display what's happened ... svgBuilder = SVGBuilder() #svgBuilder.GlobalStyle.showCoords = True svgBuilder.GlobalStyle.fillType = pft svgBuilder.AddPolygons(subj, 0x402020FF, 0x802020FF) #svgBuilder.GlobalStyle.showCoords = False svgBuilder.AddPolygons(clip, 0x40FFFF20, 0x80FF2020) svgBuilder.GlobalStyle.penWidth = 0.6 svgBuilder.AddPolygons(solution, 0x60138013, 0xFF003300) holes = [] for poly in solution: if Area(poly) < 0: holes.append(poly) svgBuilder.AddPolygons(holes, 0x0, 0xFFFF0000) svgBuilder.SaveToFile('./test.svg', invScale, 100) if result: os.startfile('test.svg') # call(('open', 'test.svg')) # print("finished") # else: print("failed")