topic Re: How to efficiently process 10000 individual coordinates to build model. in Fusion API and Scripts Forum

How to efficiently process 10000 individual coordinates to build model.

Anonymous — Fri, 28 Feb 2020 03:11:46 GMT

I am trying to build a 3D solid model using sweep / extrusion feature. I have collection of more than 100000 individual coordinates. I am taking two points for extrusion/ sweep at a time to to avoid self intersecting lines problem. However it's taking forever to finish building all bodies (whole model). Is there any way to speed up this process efficiently. I could not understand the use of BRep bodies. I do not need intermediary bodies to be saved. I am trying to get whole extruded model. I have few questions.

1) How to combine all bodies to make single entity?

2) How to use revolve feature that joins the empty spaces between to extrusion bodies at angle?

3)How to find segments untill the intersecting points to use sweep feature so, it doesn't require to take only two points at a time.

I am very new to scripting and making model in Fusion 360.

I would really appreciate if you could help me out in this.

Here is my code: (Any suggestion for improvement would be highly appreciated.)

import adsk.core, adsk.fusion, traceback
import io

def Extrudeline(spline,distance,radius,rootComp):
path = rootComp.features.createPath(spline)
# create construction plane normal to the spline
planes = rootComp.constructionPlanes
planeInput = planes.createInput()
planeInput.setByDistanceOnPath(path, adsk.core.ValueInput.createByReal(0))
plane = planes.add(planeInput)
extrudes = rootComp.features.extrudeFeatures

sketches = rootComp.sketches
sketch = sketches.add(plane)

center = plane.geometry.origin
center = sketch.modelToSketchSpace(center)
sketch.sketchCurves.sketchCircles.addByCenterRadius(center, radius)
profile = sketch.profiles[0]

#create a extrude input
distances = adsk.core.ValueInput.createByReal(distance)
extrudeInput = extrudes.createInput(profile, adsk.fusion.FeatureOperations.NewBodyFeatureOperation)
extent_distance = adsk.fusion.DistanceExtentDefinition.create(distances)
extrudeInput.setOneSideExtent(extent_distance, adsk.fusion.ExtentDirections.PositiveExtentDirection)
extrude2 = extrudes.add(extrudeInput)
return extrude2

def run(context):
ui = None
try:
app = adsk.core.Application.get()
ui = app.userInterface
# Get all components in the active design.
product = app.activeProduct
design = adsk.fusion.Design.cast(product)
title = 'Import Spline csv'
if not design:
ui.messageBox('No active Fusion design', title)
return

#dlg = ui.createFileDialog()
#dlg.title = 'Open CSV File'
#dlg.filter = 'Comma Separated Values (*.csv);;All Files (*.*)'
#if dlg.showOpen() != adsk.core.DialogResults.DialogOK :
# return

#filename = dlg.filename

filename = '/Users/rizal/Desktop/20mm_cubeCoordinates.txt'

with io.open(filename, 'r',buffering=2000000000,encoding='utf-8-sig') as f:

data = []
line = next(f)

for line in f:
pntStrArr = line.split(',')
try:
data.append([float(pntStr) for pntStr in pntStrArr])
except:
continue

length = len(data)
length=1101
temp =0
radius =0.2

while temp<length - 1:

root = design.rootComponent
sketch = root.sketches.add(root.xYConstructionPlane)
sketch.isComputeDeferred = True
distance = ((data[temp][0]-data[temp+1][0])**2+(data[temp][1]-data[temp+1][1])**2+(data[temp][2]-data[temp+1][2])**2)**0.5
line = sketch.sketchCurves.sketchLines.addByTwoPoints(adsk.core.Point3D.create(data[temp][0], data[temp][1], data[temp][2]),adsk.core.Point3D.create(data[temp+1][0], data[temp+1][1], data[temp+1][2]))

temp += 1

extrude1 = Extrudeline(line,distance,radius,root)
sketch.isComputeDeferred = False

except:
if ui:
ui.messageBox('Failed:\n{}'.format(traceback.format_exc()))

Re: How to efficiently process 10000 individual coordinates to build model.

kandennti — Fri, 28 Feb 2020 08:52:26 GMT

Hi @Anonymous .

I don't understand what we want, but we interpret it as a question about speeding up modeling.

One choice is to use a TemporaryBRepManager Object.

http://help.autodesk.com/view/fusion360/ENU/?guid=GUID-2CA08FA9-8C23-4449-87C4-1CCBA87941C1

I have never used it, but I think fscad created by @JesusFreke is also good.

https://forums.autodesk.com/t5/fusion-360-api-and-scripts/fscad-a-framework-for-programmatic-cad-in-fusion-360/m-p/8601740

If it is not limited to Fusion360, OpenSCAD may be an choice.

https://www.openscad.org/

Re: How to efficiently process 10000 individual coordinates to build model.

Anonymous — Sat, 29 Feb 2020 06:44:05 GMT

Hi @kandennti

Thank you for your response.

Actually I am trying to extrude a circular profile along the path(line) created by the 3D point objects. However, the whole path has multiple self-intersecting places, due to which I am not being able to sweep or extrude at once. To avoid the self-intersecting exceptions, I am taking two (3D point object) coordinates to extrude a body by JoinExtrudeFeature object and running this in the while loop until the end of points. At a minimum, the total number of points is 1800. If I run the 100 points, it gives extrusion in 15 mins but if add further points or try to run the whole collection of points, the Fusion 360 freezes for all day and keeps freezing. I have been able to get an output of 40% data.

But could not build the model out of all the points. Please help me if there is any efficient way to do this. I would really appreciate your help.

#Here is a code I running to build the model

#Any recommendations and feedbacks would help me greatly.

import adsk.core, adsk.fusion, traceback
import io

def Extrudeline(spline,distance,radius,rootComp):

path = rootComp.features.createPath(spline)
# create construction plane normal to the spline
planes = rootComp.constructionPlanes
planeInput = planes.createInput()
planeInput.setByDistanceOnPath(path, adsk.core.ValueInput.createByReal(0))
plane = planes.add(planeInput)
extrudes = rootComp.features.extrudeFeatures

#create a extrude input
distances = adsk.core.ValueInput.createByReal(distance)
extrudeInput = extrudes.createInput(profile, adsk.fusion.FeatureOperations.JoinFeatureOperation)
extent_distance = adsk.fusion.DistanceExtentDefinition.create(distances)
extrudeInput.setOneSideExtent(extent_distance, adsk.fusion.ExtentDirections.PositiveExtentDirection)
extrude2 = extrudes.add(extrudeInput)

return extrude2

filename = '/Desktop/20mm_cubeCoordinates.txt'

with io.open(filename, 'r',buffering=2000000000,encoding='utf-8-sig') as f:

data = []
line = next(f)

for line in f:
pntStrArr = line.split(',')
try:
data.append([float(pntStr) for pntStr in pntStrArr])
except:
continue

length = len(data)
#length=676
temp =0
radius =0.035

while temp<length - 1:

root = design.rootComponent
sketch = root.sketches.add(root.xYConstructionPlane)
sketch.isComputeDeferred = False
distance = ((data[temp][0]/10-data[temp+1][0]/10)**2+(data[temp][1]/10-data[temp+1][1]/10)**2+(data[temp][2]/10-data[temp+1][2]/10)**2)**0.5
line = sketch.sketchCurves.sketchLines.addByTwoPoints(adsk.core.Point3D.create(data[temp][0]/10, data[temp][1]/10, data[temp][2]/10),adsk.core.Point3D.create(data[temp+1][0]/10, data[temp+1][1]/10, data[temp+1][2]/10))

temp += 1

extrude1 = Extrudeline(line,distance,radius,root)

except:
if ui:
ui.messageBox('Failed:\n{}'.format(traceback.format_exc()))

Re: How to efficiently process 10000 individual coordinates to build model.

kandennti — Sat, 29 Feb 2020 16:01:10 GMT

Use this to post your code on the forum. Keep indentation.

Python indentation is important.

I tried using TemporaryBRepManager Object.

#Fusion360API Python script import adsk.core, adsk.fusion, traceback import time def run(context): ui = None try: # data file path path = '/Desktop/20mm_cubeCoordinates.txt' app = adsk.core.Application.get() ui = app.userInterface # new document app.documents.add(adsk.core.DocumentTypes.FusionDesignDocumentType) des = adsk.fusion.Design.cast(app.activeProduct) # direct design des.designType = adsk.fusion.DesignTypes.DirectDesignType # time startTime = time.time() # read file data = readData(path) # point3D pnt3D = adsk.core.Point3D pnts = [pnt3D.create(x, y, z) for (x, y, z) in data] # create cylinder radius =0.035 cylinders = initCylinders(pnts, radius) # add bodies root = des.rootComponent bodies = root.bRepBodies [bodies.add(cyl) for cyl in cylinders] # fin des.designType = adsk.fusion.DesignTypes.ParametricDesignType msg = '-- Done --\nBodyCount : {}\nTime : {} s'.format( len(cylinders),time.time() - startTime) ui.messageBox(msg) except: if ui: ui.messageBox('Failed:\n{}'.format(traceback.format_exc())) def initCylinders( pnts :list, radius :float ) -> list: tmpMgr = adsk.fusion.TemporaryBRepManager.get() vec = adsk.core.Vector3D.create(1.0, 1.0, 1.0) cyls = [] for idx in range(len(pnts) - 1): cyl = tmpMgr.createEllipticalCylinderOrCone( pnts[idx], radius, radius, pnts[idx + 1], radius, vec) if cyl: cyls.append(cyl) return cyls def readData( path :str ) -> list: with open(path, 'r', buffering = -1, encoding = 'utf-8-sig') as f: lst = f.readlines() data = [s.replace('\n', '').split(',') for s in lst] app :adsk.core.Application = adsk.core.Application.get() des :adsk.fusion.Design = app.activeProduct um = des.unitsManager conv = um.convert(1, um.defaultLengthUnits, um.internalUnits) return [(float(x) * conv, float(y) * conv, float(z) * conv) for (x,y,z) in data]

The PC at home was low performance, but could be processed in about 90 seconds.

Re: How to efficiently process 10000 individual coordinates to build model.

Anonymous — Mon, 02 Mar 2020 02:48:13 GMT

Hi @kandennti

Thank you for helping me out with this.

I was able to successfully to get the expected result with your help. However, When I try to take 10,000 point, the Fusion crashes and closes by itself due to segmentation error.

I have got into a new problem with this. I would really appreciate it if you could help to.

1) I am combining the bodies so that I can use the revolve feature. But after combining two bodies I am not being able to select the face that I want to revolve. In the

#Fusion360API Python script import adsk.core, adsk.fusion, traceback import time def run(context): ui = None try: # data file path path = '/Users/rizal/Desktop/20mm_cubeCoordinates.txt' app = adsk.core.Application.get() ui = app.userInterface # new document # app.documents.add(adsk.core.DocumentTypes.FusionDesignDocumentType) des = adsk.fusion.Design.cast(app.activeProduct) # direct design des.designType = adsk.fusion.DesignTypes.DirectDesignType # time startTime = time.time() # read file data = readData(path) # point3D pnt3D = adsk.core.Point3D pnts = [pnt3D.create(x, y, z) for (x, y, z) in data] # create cylinder radius =0.035 cylinders = initCylinders(pnts, radius) # add bodies root = des.rootComponent bodies = root.bRepBodies # for i in range(bodies.count-1): [bodies.add(cyl) for cyl in cylinders] ################# #Sketching lines to make axis for revolve sketches = root.sketches xyPlane = root.xYConstructionPlane sketch = sketches.add(xyPlane) lines = sketch.sketchCurves.sketchLines axisline = lines.addByTwoPoints( adsk.core.Point3D.create(pnts[1].x,pnts[1].y,pnts[1].z), adsk.core.Point3D.create(pnts[1].x,pnts[1].y,pnts[1].z+0.5)) #Collection object to store bodies for combine operation coll = adsk.core.ObjectCollection.create() coll.add(bodies.item(1)) combine = root.features.combineFeatures combineInput = combine.createInput(bodies.item(0),coll) body = combine.add(combineInput) # face1 = body.faces[2] # face2 = body.faces[3] # revolves = root.features.revolveFeatures # revInput = revolves.createInput(face1,axisline,adsk.fusion.FeatureOperations.JoinFeatureOperation) # revInput.setOneSideToExtent(face2) # ext = revolves.add(revInput) # # for i in range(bodies.count-1): ############ # fin des.designType = adsk.fusion.DesignTypes.ParametricDesignType msg = '-- Done --\nBodyCount : {}\nTime : {}'.format( len(cylinders),time.time() - startTime) ui.messageBox(msg) except: if ui: ui.messageBox('Failed:\n{}'.format(traceback.format_exc())) def initCylinders( pnts :list, radius :float ) -> list: tmpMgr = adsk.fusion.TemporaryBRepManager.get() cyls = [] for idx in range(2): # creating vector to orient epllipse along the path vec = adsk.core.Vector3D.create(pnts[idx].x,pnts[idx].y,pnts[idx].z) #creating elliptical cylinder with major radius 0.02 and minor radius 0.01 and direction of vector cyl = tmpMgr.createEllipticalCylinderOrCone( pnts[idx],0.02, 0.01, pnts[idx + 1], 0.02, vec) if cyl: cyls.append(cyl) return cyls def readData( path :str ) -> list: with open(path, 'r', buffering = -1, encoding = 'utf-8-sig') as f: lst = f.readlines() data = [s.replace('\n', '').split(',') for s in lst] app :adsk.core.Application = adsk.core.Application.get() des :adsk.fusion.Design = app.activeProduct um = des.unitsManager conv = um.convert(1, um.defaultLengthUnits, um.internalUnits) return [(float(x) * conv, float(y) * conv, float(z) * conv) for (x,y,z) in data]

Two combined bodiesSelecting the face (I want to select with Script )Using revolve (needs to be done with scripting)user interface, i can select the face and revolve it to starting face of the next body. However, with the script, I am not being able to select the end face of 1st body to revolve.

Could you please look into the updated code and help me to figure out how to solve this issue.

Re: How to efficiently process 10000 individual coordinates to build model.

kandennti — Mon, 02 Mar 2020 06:01:10 GMT

Using sketches is slow.

Therefore, I decided to create a sphere at the end of the cylinder.

Creating a sphere is everything except the beginning and end.

The combination of all these cylinders and spheres was also handled by the API.

#Fusion360API Python script import adsk.core, adsk.fusion, traceback import time def run(context): ui = None try: # data file path path = '/Desktop/20mm_cubeCoordinates.txt' app = adsk.core.Application.get() ui = app.userInterface # new document app.documents.add(adsk.core.DocumentTypes.FusionDesignDocumentType) des = adsk.fusion.Design.cast(app.activeProduct) # direct design des.designType = adsk.fusion.DesignTypes.DirectDesignType # time startTime = time.time() # read file data = readData(path) # point3D pnt3D = adsk.core.Point3D pnts = [pnt3D.create(x, y, z) for (x, y, z) in data] # create cylinder radius = 0.035 cylinders = initCylinders(pnts, radius) # create Sphere spheres = initSpheres(pnts, radius) # union body = unionBodies(unionBodies(cylinders[0], cylinders[1:]),spheres) # add bodies root = des.rootComponent bodies = root.bRepBodies bodies.add(body) # fin des.designType = adsk.fusion.DesignTypes.ParametricDesignType msg = '-- Done --\nPointCount : {}\nTime : {} s'.format( len(pnts),time.time() - startTime) ui.messageBox(msg) except: if ui: ui.messageBox('Failed:\n{}'.format(traceback.format_exc())) def unionBodies( body :adsk.fusion.BRepBody, bodies :list, ) -> adsk.fusion.BRepBody: tmpMgr = adsk.fusion.TemporaryBRepManager.get() unionBool = adsk.fusion.BooleanTypes.UnionBooleanType target :adsk.fusion.BRepBody = body for tool in bodies: tmpMgr.booleanOperation(target, tool, unionBool) return target def initSpheres( pnts :list, radius :float ) -> list: tmpMgr = adsk.fusion.TemporaryBRepManager.get() sphs = [] for pnt in pnts[1:-1]: sph = tmpMgr.createSphere(pnt, radius) if sph: sphs.append(sph) return sphs def initCylinders( pnts :list, radius :float ) -> list: tmpMgr = adsk.fusion.TemporaryBRepManager.get() vec = adsk.core.Vector3D.create(1.0, 1.0, 1.0) cyls = [] for idx in range(len(pnts) - 1): cyl = tmpMgr.createEllipticalCylinderOrCone( pnts[idx], radius, radius, pnts[idx + 1], radius, vec) if cyl: cyls.append(cyl) return cyls def readData( path :str ) -> list: with open(path, 'r', buffering = -1, encoding = 'utf-8-sig') as f: lst = f.readlines() data = [s.replace('\n', '').split(',') for s in lst] app :adsk.core.Application = adsk.core.Application.get() des :adsk.fusion.Design = app.activeProduct um = des.unitsManager conv = um.convert(1, um.defaultLengthUnits, um.internalUnits) return [(float(x) * conv, float(y) * conv, float(z) * conv) for (x,y,z) in data]

However, When I try to take 10,000 point, the Fusion crashes and closes by itself due to segmentation error.

No more ideas. Why not split and process the data?

Re: How to efficiently process 10000 individual coordinates to build model.

Anonymous — Mon, 02 Mar 2020 06:41:48 GMT

Hi @kandennti ,

Thats a great solution.I tried it and gives a good result if the cylinder is circluar.

I am wodering how would you approach if the cylinder is elliptical shape.

For the data, i tried doing it in fractions and adding to component.

I have another file with 5075 points. My fusion app works untill 5000 points (at single or doing multiple time and combining bodies) however if i add another faraction of 20 points after 5000 my fusion crashes and closes. I dont know why its not being able to 75 more points.

Re: How to efficiently process 10000 individual coordinates to build model.

Anonymous — Mon, 02 Mar 2020 06:45:55 GMT

Hi @kandennti ,

It's my humble request if you could give an idea if the shape is elliptical.

Would you mind trying another file? One file is a whole object.

I would really appreciate it if you could give me feedback on it.

new.txt has 5075 points and

shape.txt has hell loads of points.

Re: How to efficiently process 10000 individual coordinates to build model.

kandennti — Mon, 02 Mar 2020 10:45:49 GMT

Attach f3d file and script file which processed "new.txt".

However, it cannot be determined whether the processing is correct.

"shape.txt" can create BRepBody,

bodies.add(body)

Is a difficult state.

From tomorrow, my work may be so busy that I can't afford the time to tackle this problem.

Re: How to efficiently process 10000 individual coordinates to build model.

kandennti — Mon, 09 Mar 2020 06:24:02 GMT

The last one was incorrect.

Attach the slightly improved data.

I tried "shape.txt", but the data is big and Fusion360 goes down. It can be processed halfway, but ...

This is the limit for me.

Re: How to efficiently process 10000 individual coordinates to build model.

Anonymous — Sun, 15 Mar 2020 21:39:19 GMT

Hi @kandennti ,

I would like to express my sincere thanks for your help and time you have spent in finding a solution to the problem I encountered.

Your solution helped me greatly to move forward in my project.

Thanks,

Khem Raj Rijal