From 57307528863e2fc8c1c1c7d0489603ce0686e5f0 Mon Sep 17 00:00:00 2001 From: jaseg Date: Tue, 28 Dec 2021 21:40:39 +0100 Subject: WIP --- gerbonara/gerber/aperture_macros/parse.py | 153 +++++++++++++++++++++++++++ gerbonara/gerber/gerber_primitives.py | 14 +++ gerbonara/gerber/graphic_objects.py | 165 ++++++++++++++++++++++++++++++ 3 files changed, 332 insertions(+) create mode 100644 gerbonara/gerber/aperture_macros/parse.py create mode 100644 gerbonara/gerber/gerber_primitives.py create mode 100644 gerbonara/gerber/graphic_objects.py diff --git a/gerbonara/gerber/aperture_macros/parse.py b/gerbonara/gerber/aperture_macros/parse.py new file mode 100644 index 0000000..47f45d1 --- /dev/null +++ b/gerbonara/gerber/aperture_macros/parse.py @@ -0,0 +1,153 @@ +#!/usr/bin/env python +# -*- coding: utf-8 -*- + +# Copyright 2021 Jan Götte + +import operator +import re +import ast +import copy +import math + +import primitive as ap +from expression import * + +from .. import apertures + +def rad_to_deg(x): + return (x / math.pi) * 180 + +def _map_expression(node): + if isinstance(node, ast.Num): + return ConstantExpression(node.n) + + elif isinstance(node, ast.BinOp): + op_map = {ast.Add: operator.add, ast.Sub: operator.sub, ast.Mult: operator.mul, ast.Div: operator.truediv} + return OperatorExpression(op_map[type(node.op)], _map_expression(node.left), _map_expression(node.right)) + + elif isinstance(node, ast.UnaryOp): + if type(node.op) == ast.UAdd: + return _map_expression(node.operand) + else: + return OperatorExpression(operator.sub, ConstantExpression(0), _map_expression(node.operand)) + + elif isinstance(node, ast.Name): + return VariableExpression(int(node.id[3:])) # node.id has format var[0-9]+ + + else: + raise SyntaxError('Invalid aperture macro expression') + +def _parse_expression(expr): + expr = expr.lower().replace('x', '*') + expr = re.sub(r'\$([0-9]+)', r'var\1', expr) + try: + parsed = ast.parse(expr, mode='eval').body + except SyntaxError as e: + raise SyntaxError('Invalid aperture macro expression') from e + return _map_expression(parsed) + +class ApertureMacro: + def __init__(self, name=None, primitives=None, variables=None): + self._name = name + self.comments = [] + self.variables = variables or {} + self.primitives = primitives or [] + + @classmethod + def parse_macro(cls, name, macro, unit): + macro = cls(name) + + blocks = re.sub(r'\s', '', macro).split('*') + for block in blocks: + if not (block := block.strip()): # empty block + continue + + if block[0:1] == '0 ': # comment + macro.comments.append(Comment(block[2:])) + + if block[0] == '$': # variable definition + name, expr = block.partition('=') + number = int(name[1:]) + if number in macro.variables: + raise SyntaxError(f'Re-definition of aperture macro variable {number} inside macro') + macro.variables[number] = _parse_expression(expr) + + else: # primitive + primitive, *args = block.split(',') + args = [_parse_expression(arg) for arg in args] + primitive = ap.PRIMITIVE_CLASSES[int(primitive)](unit=unit, args=args + macro.primitives.append(primitive) + + @property + def name(self): + if self.name is not None: + return self.name + else: + return f'gn_{hash(self)}' + + def __str__(self): + return f'' + + def __eq__(self, other): + return hasattr(other, to_gerber) and self.to_gerber() == other.to_gerber() + + def __hash__(self): + return hash(self.to_gerber()) + + def to_gerber(self, unit=None): + comments = [ c.to_gerber() for c in self.comments ] + variable_defs = [ f'${var.to_gerber(unit)}={expr}' for var, expr in self.variables.items() ] + primitive_defs = [ prim.to_gerber(unit) for prim in self.primitives ] + return '*\n'.join(comments + variable_defs + primitive_defs) + + def to_graphic_primitives(self, offset, rotation:'radians', parameters : [float], unit=None): + variables = dict(self.variables) + for number, value in enumerate(parameters): + if i in variables: + raise SyntaxError(f'Re-definition of aperture macro variable {i} through parameter {value}') + variables[i] = value + + return [ primitive.to_graphic_primitives(offset, rotation, variables, unit) for primitive in self.primitives ] + + def rotated(self, angle): + copy = copy.deepcopy(self) + for primitive in copy.primitives: + primitive.rotation += rad_to_deg(angle) + return copy + + +class GenericMacros: + deg_per_rad = 180 / math.pi + cons, var = VariableExpression + _generic_hole = lambda n: [ + ap.Circle(exposure=0, diameter=var(n), x=0, y=0), + ap.Rectangle(exposure=0, w=var(n), h=var(n+1), x=0, y=0, rotation=var(n+2) * deg_per_rad)] + + circle = ApertureMacro([ + ap.Circle(exposure=1, diameter=var(1), x=0, y=0, rotation=var(4) * deg_per_rad), + *_generic_hole(2)]) + + rect = ApertureMacro([ + ap.Rectangle(exposure=1, w=var(1), h=var(2), x=0, y=0, rotation=var(5) * deg_per_rad), + *_generic_hole(3) ]) + + # w must be larger than h + obround = ApertureMacro([ + ap.Rectangle(exposure=1, w=var(1), h=var(2), x=0, y=0, rotation=var(5) * deg_per_rad), + ap.Circle(exposure=1, diameter=var(2), x=+var(1)/2, y=0, rotation=var(5) * deg_per_rad), + ap.Circle(exposure=1, diameter=var(2), x=-var(1)/2, y=0, rotation=var(5) * deg_per_rad), + *_generic_hole(3) ]) + + polygon = ApertureMacro([ + ap.Polygon(exposure=1, n_vertices=var(2), x=0, y=0, diameter=var(1), rotation=var(3) * deg_per_rad), + pa.Circle(exposure=0, diameter=var(4), x=0, y=0)]) + + +if __name__ == '__main__': + import sys + #for line in sys.stdin: + #expr = _parse_expression(line.strip()) + #print(expr, '->', expr.optimized()) + + for primitive in parse_macro(sys.stdin.read(), 'mm'): + print(primitive) diff --git a/gerbonara/gerber/gerber_primitives.py b/gerbonara/gerber/gerber_primitives.py new file mode 100644 index 0000000..1a40948 --- /dev/null +++ b/gerbonara/gerber/gerber_primitives.py @@ -0,0 +1,14 @@ + +from apertures import * + +class GerberPrimitive: + pass + + def to_graphic_primitives(self): + pass + +class Line(GerberPrimitive): + pass + +class Arc(GerberPrimitive): + pass diff --git a/gerbonara/gerber/graphic_objects.py b/gerbonara/gerber/graphic_objects.py new file mode 100644 index 0000000..55d1b9c --- /dev/null +++ b/gerbonara/gerber/graphic_objects.py @@ -0,0 +1,165 @@ + +import graphic_primitives as gp + +class GerberObject: + _ : KW_ONLY + polarity_dark : bool = True + + def to_primitives(self): + raise NotImplementedError() + +@dataclass +class Flash(GerberObject): + x : float + y : float + aperture : object + + def with_offset(self, dx, dy): + return replace(self, x=self.x+dx, y=self.y+dy) + + def rotate(self, rotation, cx=None, cy=None): + self.x, self.y = gp.rotate_point(self.x, self.y, rotation, cx, cy) + + def to_primitives(self): + yield from self.aperture.flash(self.x, self.y) + + def to_statements(self, gs): + yield from gs.set_polarity(self.polarity_dark) + yield from gs.set_aperture(self.aperture) + yield FlashStmt(self.x, self.y) + +class Region(GerberObject): + def __init__(self, outline=[], arc_centers=None, *, polarity_dark): + super().__init__(self, polarity_dark=polarity_dark) + self.poly = gp.ArcPoly() + + def with_offset(self, dx, dy): + return Region([ (x+dx, y+dy) for x, y in outline ], radii, polarity_dark=self.polarity_dark) + + def rotate(self, angle, cx=0, cy=0): + self.poly.outline = [ gp.rotate_point(x, y, angle, cx, cy) for x, y in self.poly.outline ] + self.poly.arc_centers = [ gp.rotate_point(x, y, angle, cx, cy) for x, y in self.poly.arc_centers ] + + def append(self, obj): + if not self.outline: + self.poly.outline.append(obj.p1) + self.poly.outline.append(obj.p2) + + if isinstance(obj, Arc): + self.poly.arc_centers.append(obj.center) + else: + self.poly.arc_centers.append(None) + + def to_primitives(self): + self.poly.polarity_dark = polarity_dark + yield self.poly + + def to_statements(self, gs): + yield RegionStartStmt() + + yield from gs.set_current_point(self.poly.outline[0]) + + for point, arc_center in zip(self.poly.outline, self.poly.arc_centers): + if arc_center is None: + yield from gs.set_interpolation_mode(LinearModeStmt) + yield InterpolateStmt(*point) + + else: + cx, cy = arc_center + x2, y2 = point + yield from gs.set_interpolation_mode(CircularCCWModeStmt) + yield InterpolateStmt(x2, y2, cx-x2, cy-y2) + + yield RegionEndStmt() + + +class Line(GerberObject): + # Line with *round* end caps. + x1 : float + y1 : float + x2 : float + y2 : float + aperture : object + + def with_offset(self, dx, dy): + return replace(self, x1=self.x1+dx, y1=self.y1+dy, x2=self.x2+dx, y2=self.y2+dy) + + def rotate(self, rotation, cx=None, cy=None): + if cx is None: + cx = (self.x1 + self.x2) / 2 + cy = (self.y1 + self.y2) / 2 + self.x1, self.y1 = gp.rotate_point(self.x1, self.y1, rotation, cx, cy) + self.x2, self.y2 = gp.rotate_point(self.x2, self.y2, rotation, cx, cy) + + @property + def p1(self): + return self.x1, self.y1 + + @property + def p2(self): + return self.x2, self.y2 + + def to_primitives(self): + yield gp.Line(*self.p1, *self.p2, self.aperture.equivalent_width, polarity_dark=self.polarity_dark) + + def to_statements(self, gs): + yield from gs.set_aperture(self.aperture) + yield from gs.set_interpolation_mode(LinearModeStmt) + yield from gs.set_current_point(self.p1) + yield InterpolateStmt(*self.p2) + + +class Arc(GerberObject): + x : float + y : float + r : float + angle1 : float # radians! + angle2 : float # radians! + aperture : object + + @classmethod + def from_coords(kls, start, end, center_delta, aperture, flipped=False, polarity_dark=True): + x0, y0 = start + x1, y1 = end + dx, dy = center_delta + cx, cy = x0+dx, y0+dy + angle1 = math.atan2(y0-cy, x0-cx) + angle2 = math.atan2(y1-cy, x1-cx) + aperture = self.aperture + if flipped: + angle1, angle2 = angle2, angle1 + r = math.sqrt(dx**2 + dy**2) + # r should be approximately (depending on coordinate resolution) equal for center->start and center->end + return kls(cx, cy, r, angle1, angle2, polarity_dark=polarity_dark) + + def with_offset(self, dx, dy): + return replace(self, x=self.x+dx, y=self.y+dy) + + @property + def p1(self): + return self.x + self.r*sin(self.angle1), self.y + self.r*cos(self.angle1) + + @property + def p2(self): + return self.x + self.r*sin(self.angle2), self.y + self.r*cos(self.angle2) + + @property + def center(self): + return (self.x, self.y) + + def rotate(self, rotation, cx=None, cy=None): + self.x, self.y = gp.rotate_point(self.x, self.y, rotation, cx, cy) + self.angle1 = (self.angle1+rotation) % (2*math.pi) + self.angle2 = (self.angle2+rotation) % (2*math.pi) + + def to_primitives(self): + yield gp.Arc(self.x, self.y, self.r, self.angle1, self.angle2, self.aperture.equivalent_width, polarity_dark=self.polarity_dark) + + def to_statements(self, gs): + yield from gs.set_aperture(self.aperture) + yield from gs.set_interpolation_mode(CircularCCWModeStmt) + yield from gs.set_current_point(self.p1) + x2, y2 = self.p2 + yield InterpolateStmt(x2, y2, self.x-x2, self.y-y2) + + -- cgit