def __init__(self, cgal_elem, swoop_elem):
from CGAL.CGAL_Kernel import \
Segment_2,\
Polygon_2,\
Point_2,\
Bbox_2,\
Iso_rectangle_2,\
do_intersect
self.cgal_elem = cgal_elem
self.swoop_elem = swoop_elem
bbox = cgal_elem.bbox()
self.rect = Rectangle.from_cgal_bbox(bbox, check=False)
#do_intersect only works with iso_rectangle, not bbox
self.iso_rect = Iso_rectangle_2(bbox.xmin(), bbox.ymin(), bbox.xmax(), bbox.ymax())
def __init__(self, cgal_elem, swoop_elem):
from CGAL.CGAL_Kernel import \
Segment_2,\
Polygon_2,\
Point_2,\
Bbox_2,\
Iso_rectangle_2,\
do_intersect
self.cgal_elem = cgal_elem
self.swoop_elem = swoop_elem
bbox = cgal_elem.bbox()
self.rect = Rectangle.from_cgal_bbox(bbox, check=False)
#do_intersect only works with iso_rectangle, not bbox
self.iso_rect = Iso_rectangle_2(bbox.xmin(), bbox.ymin(), bbox.xmax(),
bbox.ymax())
def get_bounding_box(self, layer=None, type=None):
"""
Get the minimum bounding box enclosing this list of primitive elements
More accurate than self._get_cgal_elem().bbox(), because it accounts for segment width
:param layer: Swoop layer to filter on
"""
def max_min(vertex_list, width=None):
max = np.maximum.reduce(vertex_list)
min = np.minimum.reduce(vertex_list)
if width is not None:
radius = np.array([width, width]) / 2.0
max += radius
min -= radius
return max,min
if isinstance(self, Swoop.Polygon):
vertices = [v.get_point() for v in self.get_vertices()]
return Rectangle(*max_min(vertices, self.get_width()))
elif isinstance(self, Swoop.Wire):
if self.get_curve() is not None:
theta = math.radians(self.get_curve()) # angle swept by arc
theta = math.fmod(theta, 2*math.pi)
p1 = self.get_point(0) # 2 points on the circle
p2 = self.get_point(1)
return arc_bounding_box(p1, p2, theta).pad(self.get_width() / 2.0)
else:
vertices = [self.get_point(0), self.get_point(1)]
return Rectangle(*max_min(vertices, self.get_width()), check=False)
elif isinstance(self, Swoop.Rectangle):
#get_cgal_elem already handles rotation
bbox = self._get_cgal_elem().bbox()
return Rectangle.from_cgal_bbox(bbox, check=False)
elif isinstance(self, Swoop.Via) or isinstance(self, Swoop.Pad):
#These assume default settings
#Unfortunately, restring can change the sizes of things after import
center = self.get_point()
if self.get_diameter() is None:
# Gotta figure it out from drill
if self.get_drill() < 1.0:
diameter = self.get_drill() + 0.5
else:
diameter = self.get_drill()*1.5
else:
diameter = self.get_diameter()
radius = diameter/2.0 * np.ones(2)
angle = 0.0
if isinstance(self, Swoop.Pad) and self.get_rot() is not None:
angle_m = angle_match(self.get_rot())
angle = angle_m['angle']
if angle_m['mirrored']:
angle = 180.0 - angle
rotate_matrix = Rectangle.rotation_matrix(math.radians(angle))
if self.get_shape() in ['long', 'offset']:
# This is basically a square pad with 180º arc endcaps
if self.get_shape()=='offset':
center[0] += diameter/2.0 # Center is offset to the left, correct it
rect = Rectangle(center - radius, center + radius)
verts = list(rect.vertices())
verts = map(lambda v: rotate_matrix.dot(v - center) + center, verts)
left_cap = arc_bounding_box(verts[0], verts[3], pi)
right_cap = arc_bounding_box(verts[2], verts[1], pi)
rect = Rectangle.union(rect, left_cap)
rect = Rectangle.union(rect, right_cap)
return rect
elif self.get_shape()=='octagon':
vertex = np.array([diameter/2.0,
diameter/2.0 * math.tan(2*pi / 16 )])
vertex = rotate_matrix.dot(vertex)
vertices = []
rot = Rectangle.rotation_matrix(2*pi/8)
for i in xrange(8):
vertices.append(vertex.copy() + center)
vertex = rot.dot(vertex)
return Rectangle.from_vertices(vertices, check=True)
else:
rect = Rectangle(center - radius, center + radius)
if self.get_shape()=='square':
rect.rotate_resize(angle,origin=center)
return rect
elif isinstance(self, Swoop.Smd):
center = self.get_point()
radius = np.array([self.get_dx(), self.get_dy()]) / 2.0
if self.get_rot() is not None:
angle = angle_match(self.get_rot())
radius = abs(Rectangle.rotation_matrix(math.radians(angle['angle'])).dot(radius))
return Rectangle(center - radius, center + radius)
elif isinstance(self, Swoop.Circle):
center = self.get_point()
radius = np.ones(2) * (self.get_radius() + self.get_width()/2.0)
return Rectangle(center - radius, center + radius)
elif isinstance(self, Swoop.Hole):
center = self.get_point()
radius = np.ones(2) * self.get_drill()/2.0
return Rectangle(center - radius, center + radius)
elif isinstance(self, Swoop.Element):
# Element objects do not have enough information for the bounding box
# That gets set in the constructor
if hasattr(self,"_extension_geo_elem"):
return self._extension_geo_elem.rect
else:
tform = self.get_transform()
bbox = self.find_package().get_bounding_box(layer=layer, type=type)
if bbox is None:
return None
return tform.apply(bbox)
elif isinstance(self, Swoop.Package):
rect = None
for c in self.get_children():
if ((layer is None) or (hasattr(c,"get_layer") and c.get_layer()==layer))\
and ((type is None) or isinstance(c, type)):
r = c.get_bounding_box()
if r is not None:
if rect is None:
rect = r
rect = Rectangle.union(rect, r)
return rect
else:
return None
def get_bounding_box(self, layer=None, type=None):
"""
Get the minimum bounding box enclosing this list of primitive elements
More accurate than self._get_cgal_elem().bbox(), because it accounts for segment width
:param layer: Swoop layer to filter on
"""
def max_min(vertex_list, width=None):
max = np.maximum.reduce(vertex_list)
min = np.minimum.reduce(vertex_list)
if width is not None:
radius = np.array([width, width]) / 2.0
max += radius
min -= radius
return max, min
if isinstance(self, Swoop.Polygon):
vertices = [v.get_point() for v in self.get_vertices()]
return Rectangle(*max_min(vertices, self.get_width()))
elif isinstance(self, Swoop.Wire):
if self.get_curve() is not None:
theta = math.radians(self.get_curve()) # angle swept by arc
theta = math.fmod(theta, 2 * math.pi)
p1 = self.get_point(0) # 2 points on the circle
p2 = self.get_point(1)
return arc_bounding_box(p1, p2,
theta).pad(self.get_width() / 2.0)
else:
vertices = [self.get_point(0), self.get_point(1)]
return Rectangle(*max_min(vertices, self.get_width()),
check=False)
elif isinstance(self, Swoop.Rectangle):
#get_cgal_elem already handles rotation
bbox = self._get_cgal_elem().bbox()
return Rectangle.from_cgal_bbox(bbox, check=False)
elif isinstance(self, Swoop.Via) or isinstance(self, Swoop.Pad):
#These assume default settings
#Unfortunately, restring can change the sizes of things after import
center = self.get_point()
if self.get_diameter() is None:
# Gotta figure it out from drill
if self.get_drill() < 1.0:
diameter = self.get_drill() + 0.5
else:
diameter = self.get_drill() * 1.5
else:
diameter = self.get_diameter()
radius = diameter / 2.0 * np.ones(2)
angle = 0.0
if isinstance(self, Swoop.Pad) and self.get_rot() is not None:
angle_m = angle_match(self.get_rot())
angle = angle_m['angle']
if angle_m['mirrored']:
angle = 180.0 - angle
rotate_matrix = Rectangle.rotation_matrix(math.radians(angle))
if self.get_shape() in ['long', 'offset']:
# This is basically a square pad with 180º arc endcaps
if self.get_shape() == 'offset':
center[
0] += diameter / 2.0 # Center is offset to the left, correct it
rect = Rectangle(center - radius, center + radius)
verts = list(rect.vertices())
verts = [rotate_matrix.dot(v - center) + center for v in verts]
left_cap = arc_bounding_box(verts[0], verts[3], pi)
right_cap = arc_bounding_box(verts[2], verts[1], pi)
rect = Rectangle.union(rect, left_cap)
rect = Rectangle.union(rect, right_cap)
return rect
elif self.get_shape() == 'octagon':
vertex = np.array(
[diameter / 2.0, diameter / 2.0 * math.tan(2 * pi / 16)])
vertex = rotate_matrix.dot(vertex)
vertices = []
rot = Rectangle.rotation_matrix(2 * pi / 8)
for i in range(8):
vertices.append(vertex.copy() + center)
vertex = rot.dot(vertex)
return Rectangle.from_vertices(vertices, check=True)
else:
rect = Rectangle(center - radius, center + radius)
if self.get_shape() == 'square':
rect.rotate_resize(angle, origin=center)
return rect
elif isinstance(self, Swoop.Smd):
center = self.get_point()
radius = np.array([self.get_dx(), self.get_dy()]) / 2.0
if self.get_rot() is not None:
angle = angle_match(self.get_rot())
radius = abs(
Rectangle.rotation_matrix(math.radians(
angle['angle'])).dot(radius))
return Rectangle(center - radius, center + radius)
elif isinstance(self, Swoop.Circle):
center = self.get_point()
radius = np.ones(2) * (self.get_radius() + self.get_width() / 2.0)
return Rectangle(center - radius, center + radius)
elif isinstance(self, Swoop.Hole):
center = self.get_point()
radius = np.ones(2) * self.get_drill() / 2.0
return Rectangle(center - radius, center + radius)
elif isinstance(self, Swoop.Element):
# Element objects do not have enough information for the bounding box
# That gets set in the constructor
if hasattr(self, "_extension_geo_elem"):
return self._extension_geo_elem.rect
else:
tform = self.get_transform()
bbox = self.find_package().get_bounding_box(layer=layer,
type=type)
if bbox is None:
return None
return tform.apply(bbox)
elif isinstance(self, Swoop.Package):
rect = None
for c in self.get_children():
if ((layer is None) or (hasattr(c,"get_layer") and c.get_layer()==layer))\
and ((type is None) or isinstance(c, type)):
r = c.get_bounding_box()
if r is not None:
if rect is None:
rect = r
rect = Rectangle.union(rect, r)
return rect
else:
return None