def iges_144_000():
import pyiges.IGESGeomLib as IGES
from pyiges.IGESGeomLib import IGESPoint
from pyiges.IGESCore import IGEStorage
system = IGEStorage()
examples.benchmarks.standard_iges_setup(system, "144-000-benchmark.igs")
para_spline_surface = IGES.IGESTestSplineSurf()
system.Commit(para_spline_surface)
circle = IGES.IGESGeomCircle(IGESPoint(6, 7.25, 0), IGESPoint(6.25, 7.25))
system.Commit(circle)
para_spline_curve = IGES.IGESTestSplineCurve()
system.Commit(para_spline_curve)
curve_on_parametric_surface = IGES.IGESCurveOnParametricSurface(
para_spline_surface, circle, para_spline_curve, 2)
system.Commit(curve_on_parametric_surface)
trimmed_parametric_surface = IGES.IGESTrimmedParaSurface(
para_spline_surface, curve_on_parametric_surface)
system.Commit(trimmed_parametric_surface)
system.save("144-000-benchmark.igs")
if not os.environ.get('READTHEDOCS', None):
print(system)
os.startfile("144-000-benchmark.igs")
def testrun_random_surface(filename="IGESFile_random_surface.igs"):
"""draw a randomiced surface consisting of b spline surface segments"""
system = IGEStorage()
system.StartSection.Prolog = " "
system.GlobalSection.IntegerBits = int(32)
system.GlobalSection.SPMagnitude = int(38)
system.GlobalSection.SPSignificance = int(6)
system.GlobalSection.DPMagnitude = int(38)
system.GlobalSection.DPSignificance = int(15)
system.GlobalSection.MaxNumberLineWeightGrads = int(8)
system.GlobalSection.WidthMaxLineWeightUnits = float(0.016)
system.GlobalSection.MaxCoordValue = float(71)
scale = 2.5
# data = numpy.zeros(shape=(3,5), dtype=int)
data = numpy.random.rand(100, 100)
# print("random data generated:", data.shape)
for line_nr in range(0, data.shape[0]-1):
# print("line number:", line_nr)
for row_nr in range(0, data.shape[1]-1):
x1 = line_nr * scale
x2 = (line_nr + 1) * scale
y1 = row_nr * scale
y2 = (row_nr + 1) * scale
p1 = IGESPoint(x1, y1, float(data[line_nr][row_nr]))
p2 = IGESPoint(x2, y1, float(data[line_nr+1][row_nr]))
p3 = IGESPoint(x1, y2, float(data[line_nr][row_nr+1]))
p4 = IGESPoint(x2, y2, float(data[line_nr+1][row_nr+1]))
square = IGES.IGESRationalBSplineSurface(p1, p2, p3, p4)
system.Commit(square)
# print("save file:", line_nr)
system.save(filename)
def testrun(filename="IGESFile.igs"):
system = IGEStorage()
system.StartSection.Prolog = " "
system.GlobalSection.IntegerBits = int(32)
system.GlobalSection.SPMagnitude = int(38)
system.GlobalSection.SPSignificance = int(6)
system.GlobalSection.DPMagnitude = int(38)
system.GlobalSection.DPSignificance = int(15)
system.GlobalSection.MaxNumberLineWeightGrads = int(8)
system.GlobalSection.WidthMaxLineWeightUnits = float(0.016)
system.GlobalSection.MaxCoordValue = float(71)
P = [[0, 0], [numpy.divide(1, 3),
numpy.divide(numpy.pi, 6)], [numpy.divide(2, 3), 1], [1, 1],
[1 + numpy.divide(1, 3), 1],
[1 + numpy.divide(2, 3),
numpy.divide(numpy.pi, 6)], [2, 0],
[2 + numpy.divide(1, 3), 0 - numpy.divide(numpy.pi, 6)],
[2 + numpy.divide(2, 3), 0 - 1], [3, 0 - 1],
[3 + numpy.divide(1, 3), 0 - 1],
[3 + numpy.divide(2, 3), 0 - numpy.divide(numpy.pi, 6)], [4, 0]]
for i in range(0, 1):
P.extend
P = numpy.transpose(P)
bezi = bezier_curve(P, nTimes=50)
polyln = IGES.IGESGeomPolyline()
for n in range(0, 5):
for i in range(0, len(bezi[0])):
polyln.AddPoint(IGESPoint(bezi[0][i], bezi[1][i], 10))
bezi[0] = bezi[0] + 4
system.Commit(polyln)
line = IGES.IGESGeomLine(IGESPoint(-2, -5, 0), IGESPoint(22, -5, 0))
system.Commit(line)
#system.Commit(IGES.IGESRevolve(polyln, line))
#system.Commit(IGES.IGESExtrude(polyln.DirectoryDataPointer.data, IGESPoint(0,0,10)))
print(system)
system.save(filename)
if not environ.get('READTHEDOCS', None) == 'True':
startfile(filename)
def testrun_rotate_polyline(filename="IGESFile_rotate_polyline.igs"):
"""draw a poly line and rotate it around a line to for a cylindrical object"""
system = IGEStorage()
system.StartSection.Prolog = " "
system.GlobalSection.IntegerBits = int(32)
system.GlobalSection.SPMagnitude = int(38)
system.GlobalSection.SPSignificance = int(6)
system.GlobalSection.DPMagnitude = int(38)
system.GlobalSection.DPSignificance = int(15)
system.GlobalSection.MaxNumberLineWeightGrads = int(8)
system.GlobalSection.WidthMaxLineWeightUnits = float(0.016)
system.GlobalSection.MaxCoordValue = float(71)
# add poly line from numpy array
P = [[0, 0],
[numpy.divide(1, 3), numpy.divide(numpy.pi, 6)],
[numpy.divide(2, 3), 1],
[1, 1],
[1 + numpy.divide(1, 3), 1],
[1 + numpy.divide(2, 3), numpy.divide(numpy.pi, 6)],
[2, 0],
[2 + numpy.divide(1, 3), 0 - numpy.divide(numpy.pi, 6)],
[2 + numpy.divide(2, 3), 0 - 1],
[3, 0 - 1],
[3 + numpy.divide(1, 3), 0 - 1],
[3 + numpy.divide(2, 3), 0 - numpy.divide(numpy.pi, 6)],
[4, 0]]
for i in range(0, 1):
P.extend
P = numpy.transpose(P)
bezi = bezier_curve(P, nTimes = 50)
polyln = IGES.IGESGeomPolyline()
for n in range(0, 5):
for i in range(0, len(bezi[0])):
polyln.AddPoint(IGESPoint(bezi[0][i], bezi[1][i], 10))
bezi[0] = bezi[0] + 4
system.Commit(polyln)
# add straight line
line = IGES.IGESGeomLine(IGESPoint(-2, -5, 0), IGESPoint(22, -5, 0))
system.Commit(line)
# revolve polyline around line
system.Commit(IGES.IGESRevolve(polyln, line))
system.save(filename)
def testrun_extrude_polyline(filename="IGESFile_extrude_polyline.igs"):
"""draw a poly line and extrude it to a target points"""
system = IGEStorage()
system.StartSection.Prolog = " "
system.GlobalSection.IntegerBits = int(32)
system.GlobalSection.SPMagnitude = int(38)
system.GlobalSection.SPSignificance = int(6)
system.GlobalSection.DPMagnitude = int(38)
system.GlobalSection.DPSignificance = int(15)
system.GlobalSection.MaxNumberLineWeightGrads = int(8)
system.GlobalSection.WidthMaxLineWeightUnits = float(0.016)
system.GlobalSection.MaxCoordValue = float(71)
# add poly line from numpy array
P = [[0, 0],
[numpy.divide(1, 3), numpy.divide(numpy.pi, 6)],
[numpy.divide(2, 3), 1],
[1, 1],
[1 + numpy.divide(1, 3), 1],
[1 + numpy.divide(2, 3), numpy.divide(numpy.pi, 6)],
[2, 0],
[2 + numpy.divide(1, 3), 0 - numpy.divide(numpy.pi, 6)],
[2 + numpy.divide(2, 3), 0 - 1],
[3, 0 - 1],
[3 + numpy.divide(1, 3), 0 - 1],
[3 + numpy.divide(2, 3), 0 - numpy.divide(numpy.pi, 6)],
[4, 0]]
for i in range(0, 1):
P.extend
P = numpy.transpose(P)
bezi = bezier_curve(P, nTimes = 50)
polyln = IGES.IGESGeomPolyline()
for n in range(0, 5):
for i in range(0, len(bezi[0])):
polyln.AddPoint(IGESPoint(bezi[0][i], bezi[1][i], 10))
bezi[0] = bezi[0] + 4
system.Commit(polyln)
# extrude polyline to endpoint
system.Commit(IGES.IGESExtrude(polyln, IGESPoint(0, 0, 15)))
system.save(filename)
def testrun_arc_circle(filename="IGESFile_arc_circle.igs"):
"""draw an arc and a circle, both in a plane paralell to the X-Y-plane"""
system = IGEStorage()
system.StartSection.Prolog = " "
system.GlobalSection.IntegerBits = int(32)
system.GlobalSection.SPMagnitude = int(38)
system.GlobalSection.SPSignificance = int(6)
system.GlobalSection.DPMagnitude = int(38)
system.GlobalSection.DPSignificance = int(15)
system.GlobalSection.MaxNumberLineWeightGrads = int(8)
system.GlobalSection.WidthMaxLineWeightUnits = float(0.016)
system.GlobalSection.MaxCoordValue = float(71)
arc = IGES.IGESGeomArc(-3, IGESPoint(2, 5, 0), IGESPoint(-3, 0, 0), IGESPoint(2, 10, 0))
system.Commit(arc)
circ = IGES.IGESGeomCircle(IGESPoint(5, 5, 5), 5)
system.Commit(circ)
system.save(filename)
def testrun_extrude_arc(filename="IGESFile_extrude_arc.igs"):
"""draw arc and extrude it to a target points"""
system = IGEStorage()
system.StartSection.Prolog = " "
system.GlobalSection.IntegerBits = int(32)
system.GlobalSection.SPMagnitude = int(38)
system.GlobalSection.SPSignificance = int(6)
system.GlobalSection.DPMagnitude = int(38)
system.GlobalSection.DPSignificance = int(15)
system.GlobalSection.MaxNumberLineWeightGrads = int(8)
system.GlobalSection.WidthMaxLineWeightUnits = float(0.016)
system.GlobalSection.MaxCoordValue = float(71)
arc = IGES.IGESGeomArc(0, IGESPoint(20, 11, 0), IGESPoint(20.25, 11, 0), IGESPoint(19.75, 11, 0))
system.Commit(arc)
ext = IGES.IGESExtrude(arc, IGESPoint(20.4857, 11.2357, -0.9428))
system.Commit(ext)
system.save(filename)
def iges_144_000():
import pyiges.IGESGeomLib as IGES
from pyiges.IGESCore import IGEStorage
from pyiges.IGESGeomLib import IGESPoint
filename = "144-000-example-1.igs"
system = IGEStorage()
examples.benchmarks.standard_iges_setup(system, filename)
# Circular flat surface profile
circular_flat_surface_profile = IGES.IGESGeomCircle(
IGESPoint(0, 0, 0), IGESPoint(10, 10, 0))
system.Commit(circular_flat_surface_profile)
# Make the surface from out outer profile, this is to be trimmed
circular_flat_surface = IGES.IGESGeomPlane(circular_flat_surface_profile)
system.Commit(circular_flat_surface)
# We're going to put a bunch of holes into the surface
hole = [] # Hole profiles
hole_copse = [] # Curve on Parametric Surface entries (required!)
for i in range(0, 6):
hole.append(
IGES.IGESGeomCircle(IGESPoint(-5 + i * 3, 0, 0),
IGESPoint(-5 + i * 3, 1, 0)))
system.Commit(hole[-1])
# Make those tubes, given we are not using them as a reference,
# Lets directly add them to the system.
system.Commit(IGES.IGESExtrude(hole[-1], IGESPoint(0, 0, i + 1)))
system.Commit(IGES.IGESExtrude(hole[-1], IGESPoint(0, 0, -i - 1)))
hole_copse.append(
IGES.IGESCurveOnParametricSurface(circular_flat_surface, hole[-1],
hole[-1], 1))
system.Commit(hole_copse[-1])
# Setup the surface we are to put holes into, we're using the whole surface,
# as indicated by the 0 in the outer profile parameter
trimmed_surface = IGES.IGESTrimmedParaSurface(circular_flat_surface, 0)
trimmed_surface.N1 = 0
# Lets put those holes in!
for _hole in hole_copse:
trimmed_surface.add_bounding_profile(_hole)
system.Commit(trimmed_surface)
system.save(filename)
if not os.environ.get('READTHEDOCS', None):
print(system)
os.startfile(filename)
def testrun_torus(filename="IGESFile_torus.igs"):
"""Dont know if this works, neither FreeCAD not IGES Viewer will display anyting for a torus element"""
system = IGEStorage()
system.StartSection.Prolog = " "
system.GlobalSection.IntegerBits = int(32)
system.GlobalSection.SPMagnitude = int(38)
system.GlobalSection.SPSignificance = int(6)
system.GlobalSection.DPMagnitude = int(38)
system.GlobalSection.DPSignificance = int(15)
system.GlobalSection.MaxNumberLineWeightGrads = int(8)
system.GlobalSection.WidthMaxLineWeightUnits = float(0.016)
system.GlobalSection.MaxCoordValue = float(71)
torus = IGES.IGESGeomTorus(50, 10, IGESPoint(0,0,0), IGESVector(0,0,1))
system.Commit(torus)
system.save(filename)
def iges_144_000():
import pyiges.IGESGeomLib as IGES
from pyiges.IGESCore import IGEStorage
from pyiges.IGESGeomLib import IGESPoint
filename = "144-000-example-2.igs"
system = IGEStorage()
examples.benchmarks.standard_iges_setup(system, filename)
# Setup a cone like surface, the cone is boring but whatever!
cone_surface_profile = IGES.IGESGeomPolyline(IGESPoint(-10, 0, 0),
IGESPoint(0, 0, 10))
system.Commit(cone_surface_profile)
# The centre of the revolve
center_line = IGES.IGESGeomLine(IGESPoint(0, 0, 0), IGESPoint(0, 0, 10))
system.Commit(center_line)
# Make the surface from the profile, this is to be trimmed
cone_surface = IGES.IGESRevolve(cone_surface_profile, center_line)
system.Commit(cone_surface)
#Setup the surface we are to put holes into, we're using the whole surface,
# as indicated by the 0 in the outer profile parameter
trimmed_surface = IGES.IGESTrimmedParaSurface(cone_surface, 0)
trimmed_surface.N1 = 0
hole = []
hole_copse = []
# We're going to put a bunch of holes into the surface, because we don't
# care about the index of the hole we can get away with specifying -8 to 8
for holenum in range(-8, 8):
hole.append(IGES.IGESGeomPolyline())
hole_data = [[], [], []]
# Top line following the cone profile
for i in range(21 * holenum, 20 + 21 * holenum):
hole_data[0].append(3)
hole_data[1].append(i)
hole_data[2].append(7)
# Line down the cone profile
hole_data[0].append(7)
hole_data[1].append(hole_data[1][-1])
hole_data[2].append(3)
# Bottom line following the cone profile
for i in range(20 + 21 * holenum, 21 * holenum, -1):
hole_data[0].append(7)
hole_data[1].append(i)
hole_data[2].append(3)
# Line back to start of the cone profile
hole_data[0].append(hole_data[0][0])
hole_data[1].append(hole_data[1][0])
hole_data[2].append(hole_data[2][0])
# Convert to rectangular coordinate system
hole_rect_cs = cyl_to_rect(hole_data)
# Push the points into the common format
for i in range(0, len(hole_rect_cs[0])):
hole[-1].AddPoint(
IGESPoint(hole_rect_cs[0][i], hole_rect_cs[1][i],
hole_rect_cs[2][i]))
# Commit this hole
system.Commit(hole[-1])
# Lets make a really very crazy looking thing!
# And revolve the cutting profile by a bit
# Setup Revolve Line data
random_line_data = cyl_to_rect([[-10, 10],
[(20 * holenum) + 90,
(20 * holenum) + 90], [10, 10]])
# Make revolve line
revolve_line = IGES.IGESGeomLine(
IGESPoint(random_line_data[0][0], random_line_data[1][0],
random_line_data[2][0]),
IGESPoint(random_line_data[0][1], random_line_data[1][1],
random_line_data[2][1]))
# Commit revolve line
system.Commit(revolve_line)
# Make the revolve
system.Commit(IGES.IGESRevolve(hole[-1], revolve_line, -1, 0))
# Create the Curve on Parametric Surface
hole_copse.append(
IGES.IGESCurveOnParametricSurface(cone_surface, hole[-1], hole[-1],
2))
system.Commit(hole_copse[-1])
# Lets put that holes in!
trimmed_surface.add_bounding_profile(hole_copse[-1])
# commit that surface!
system.Commit(trimmed_surface)
system.save(filename)
if not os.environ.get('READTHEDOCS', None):
print(system)
os.startfile(filename)
def iges_414_000():
import pyiges.IGESGeomLib as IGES
from pyiges.IGESCore import IGEStorage
from pyiges.IGESGeomLib import IGESPoint
system = IGEStorage()
examples.benchmarks.standard_iges_setup(system, "414-000-benchmark.igs")
line_1 = IGES.IGESGeomLine(IGESPoint(-0.13, 0, 0), IGESPoint(0, 0.25, 0))
line_2 = IGES.IGESGeomLine(IGESPoint(0, 0.25, 0), IGESPoint(0.13, 0, 0))
line_3 = IGES.IGESGeomLine(IGESPoint(0, 0.25, 0), IGESPoint(0, 0.44, 0))
line_4 = IGES.IGESGeomLine(IGESPoint(0, 0.38, 0), IGESPoint(0.09, 0.28, 0))
line_5 = IGES.IGESGeomLine(IGESPoint(0, 0.38, 0),
IGESPoint(-0.09, 0.28, 0))
circle_1 = IGES.IGESGeomArc(0, IGESPoint(0, 0.5), IGESPoint(0.06, 0.5),
IGESPoint(0.06, 0.5))
system.Commit(line_1)
system.Commit(line_2)
system.Commit(line_3)
system.Commit(line_4)
system.Commit(line_5)
system.Commit(circle_1)
person_group = IGES.IGESGroup("PERSON", line_1, line_2, line_3, line_4,
line_5, circle_1)
system.Commit(person_group)
person_array = IGES.IGESCircularArray(person_group, 3,
IGESPoint(5.5, 1.5, 0), 0.5, 0.52,
2.09)
system.Commit(person_array)
system.save("414-000-benchmark.igs")
if not os.environ.get('READTHEDOCS', None):
print(system)
os.startfile("414-000-benchmark.igs")
numpy.divide(numpy.pi, 6)], [2, 0],
[2 + numpy.divide(1, 3), 0 - numpy.divide(numpy.pi, 6)],
[2 + numpy.divide(2, 3), 0 - 1], [3, 0 - 1],
[3 + numpy.divide(1, 3), 0 - 1],
[3 + numpy.divide(2, 3), 0 - numpy.divide(numpy.pi, 6)], [4, 0]]
for i in range(0, 1):
P.extend
P = numpy.transpose(P)
bezi = bezier_curve(P, nTimes=50)
polyln = IGES.IGESGeomPolyline()
system.Commit(polyln)
line = IGES.IGESGeomLine(IGESPoint(-2, -5, 0), IGESPoint(22, -5, 0))
print(line)
system.Commit(line)
line = IGES.IGESGeomLine(IGESPoint(-2, -1, 0), IGESPoint(22, -3, 3))
print(line)
system.Commit(line)
system.Commit(IGES.IGESRevolve(polyln, line))
# system.Commit(IGES.IGESExtrude(polyln.DirectoryDataPointer.data, IGESPoint(0,0,10)))
# system.Commit(IGES.IGESRevolve(line0, line))
#system.Commit(IGES.IGESRevolve(polyln, line))
#system.Commit(IGES.IGESExtrude(polyln.DirectoryDataPointer.data, IGESPoint(0,0,10)))
def test_using_IGESPoint_radius(self):
circle = IGESGeomCircle(IGESPoint(0, 0, 0), IGESPoint(5, 5, 0))
self.assertEqual(
circle.ParameterData, [0, 0, 0, 5, 5, 5, 5],
msg='IGESGeomCircle convert from [IGESPoint, IGESPoint] list fail')
def test_using_IGESPoint(self):
point = IGESGeomPoint(IGESPoint(1, 2, 3)) # correct usage
self.assertEqual(point.ParameterData, [1, 2, 3, 0],
msg='IGESGeomPoint convert from IGESPoint fail')