def main(): compas_mesh = Mesh.from_obj(os.path.join(DATA, MODEL)) delta = get_param({}, key='delta', defaults_type='gcode') # boolean for delta printers print_volume_x = get_param({}, key='print_volume_x', defaults_type='gcode') # in mm print_volume_y = get_param({}, key='print_volume_y', defaults_type='gcode') # in mm if delta: move_mesh_to_point(compas_mesh, Point(0, 0, 0)) else: move_mesh_to_point(compas_mesh, Point(print_volume_x / 2, print_volume_y / 2, 0)) # ----- slicing slicer = PlanarSlicer(compas_mesh, slicer_type="cgal", layer_height=4.5) slicer.slice_model() generate_brim(slicer, layer_width=3.0, number_of_brim_offsets=4) simplify_paths_rdp_igl(slicer, threshold=0.6) seams_smooth(slicer, smooth_distance=10) slicer.printout_info() save_to_json(slicer.to_data(), OUTPUT_DIR, 'slicer_data.json') # ----- print organization print_organizer = PlanarPrintOrganizer(slicer) print_organizer.create_printpoints() # Set fabrication-related parameters set_extruder_toggle(print_organizer, slicer) print_organizer.printout_info() # create and output gcode gcode_parameters = {} # leave all to default gcode_text = print_organizer.output_gcode(gcode_parameters) utils.save_to_text_file(gcode_text, OUTPUT_DIR, 'my_gcode.gcode')
def main(): compas_mesh = Mesh.from_obj(os.path.join(DATA, MODEL)) move_mesh_to_point(compas_mesh, Point(0, 0, 0)) # Slicing slicer = PlanarSlicer(compas_mesh, slicer_type="cgal", layer_height=5.0) slicer.slice_model() # Sorting into vertical layers and reordering sort_into_vertical_layers(slicer, max_paths_per_layer=10) reorder_vertical_layers(slicer, align_with="x_axis") # Post-processing generate_brim(slicer, layer_width=3.0, number_of_brim_offsets=5) simplify_paths_rdp_igl(slicer, threshold=0.7) seams_smooth(slicer, smooth_distance=10) slicer.printout_info() save_to_json(slicer.to_data(), OUTPUT_DIR, 'slicer_data.json') # PlanarPrintOrganization print_organizer = PlanarPrintOrganizer(slicer) print_organizer.create_printpoints() set_extruder_toggle(print_organizer, slicer) add_safety_printpoints(print_organizer, z_hop=10.0) set_linear_velocity_constant(print_organizer, v=25.0) set_blend_radius(print_organizer, d_fillet=10.0) print_organizer.printout_info() printpoints_data = print_organizer.output_printpoints_dict() utils.save_to_json(printpoints_data, OUTPUT_DIR, 'out_printpoints.json')
def main(): start_time = time.time() ### --- Load stl compas_mesh = Mesh.from_obj(os.path.join(DATA, MODEL)) ### --- Move to origin move_mesh_to_point(compas_mesh, Point(0, 0, 0)) ### --- Slicer # options: 'default' : Both for open and closed paths. But slow # 'cgal' : Very fast. Only for closed paths. Requires additional installation (compas_cgal). slicer = PlanarSlicer(compas_mesh, slicer_type="cgal", layer_height=1.5) slicer.slice_model() ### --- Generate brim generate_brim(slicer, layer_width=3.0, number_of_brim_paths=3) ### --- Simplify the paths by removing points with a certain threshold # change the threshold value to remove more or less points simplify_paths_rdp(slicer, threshold=0.7) ### --- Smooth the seams between layers # change the smooth_distance value to achieve smoother, or more abrupt seams seams_smooth(slicer, smooth_distance=10) ### --- Prints out the info of the slicer slicer.printout_info() viewer = ObjectViewer() viewer.view.use_shaders = False slicer.visualize_on_viewer(viewer) utils.save_to_json(slicer.to_data(), OUTPUT_DIR, 'slicer_data.json') ### --- Fabrication - related information print_organizer = PrintOrganizer(slicer) print_organizer.create_printpoints(compas_mesh) print_organizer.set_extruder_toggle() print_organizer.add_safety_printpoints(z_hop=20) print_organizer.set_linear_velocity("constant", v=25) ### --- Save printpoints dictionary to json file printpoints_data = print_organizer.output_printpoints_dict() utils.save_to_json(printpoints_data, OUTPUT_DIR, 'out_printpoints.json') print_organizer.visualize_on_viewer(viewer, visualize_polyline=True, visualize_printpoints=False) viewer.update() viewer.show() end_time = time.time() print("Total elapsed time", round(end_time - start_time, 2), "seconds")
def create_setup(filename): """ Setting up the stage for testing. """ FILE = os.path.abspath(os.path.join(DATA, filename)) compas_mesh = Mesh.from_obj(FILE) slicer = PlanarSlicer(compas_mesh, slicer_type="default", layer_height=20) slicer.slice_model() generate_brim(slicer, layer_width=3.0, number_of_brim_offsets=3) simplify_paths_rdp(slicer, threshold=1.3) # seams_smooth(slicer, smooth_distance=10) slicer.printout_info() print_organizer = PlanarPrintOrganizer(slicer) print_organizer.create_printpoints() return slicer, print_organizer
def main(): start_time = time.time() ### --- Load stl compas_mesh = Mesh.from_obj(os.path.join(DATA, MODEL)) ### --- Move to origin move_mesh_to_point(compas_mesh, Point(0, 0, 0)) ### --- Slicer # try out different slicers by changing the slicer_type # options: 'default', 'meshcut', 'cgal' slicer = PlanarSlicer(compas_mesh, slicer_type="default", layer_height=1.5) slicer.slice_model() ### --- Generate brim generate_brim(slicer, layer_width=3.0, number_of_brim_paths=3) ### --- Simplify the printpaths by removing points with a certain threshold # change the threshold value to remove more or less points simplify_paths_rdp(slicer, threshold=0.9) ### --- Prints out the info of the slicer slicer.printout_info() viewer = ObjectViewer() viewer.view.use_shaders = False slicer.visualize_on_viewer(viewer) utils.save_to_json(slicer.to_data(), OUTPUT_DIR, 'slicer_data.json') ### --- Fabrication - related information print_organizer = PrintOrganizer(slicer) print_organizer.create_printpoints(compas_mesh) print_organizer.set_extruder_toggle() print_organizer.add_safety_printpoints(z_hop=20) print_organizer.set_linear_velocity("constant", v=25) ### --- Save printpoints dictionary to json file printpoints_data = print_organizer.output_printpoints_dict() utils.save_to_json(printpoints_data, OUTPUT_DIR, 'out_printpoints.json') # # print_organizer.visualize_on_viewer(viewer, visualize_polyline=True, visualize_printpoints=False) viewer.update() viewer.show() end_time = time.time() print("Total elapsed time", round(end_time - start_time, 2), "seconds")
def main(): start_time = time.time() # ========================================================================== # Load mesh # ========================================================================== compas_mesh = Mesh.from_obj(os.path.join(DATA, MODEL)) # ========================================================================== # Move to origin # ========================================================================== move_mesh_to_point(compas_mesh, Point(0, 0, 0)) # ========================================================================== # Slicing # options: 'default': Both for open and closed paths. But slow # 'cgal': Very fast. Only for closed paths. # Requires additional installation (compas_cgal). # ========================================================================== slicer = PlanarSlicer(compas_mesh, slicer_type="cgal", layer_height=1.5) slicer.slice_model() # ========================================================================== # Generate brim / raft # ========================================================================== # NOTE: Typically you would want to use either a brim OR a raft, # however, in this example both are used to explain the functionality generate_brim(slicer, layer_width=3.0, number_of_brim_offsets=4) generate_raft(slicer, raft_offset=20, distance_between_paths=5, direction="xy_diagonal", raft_layers=1) # ========================================================================== # Simplify the paths by removing points with a certain threshold # change the threshold value to remove more or less points # ========================================================================== simplify_paths_rdp_igl(slicer, threshold=0.6) # ========================================================================== # Smooth the seams between layers # change the smooth_distance value to achieve smoother, or more abrupt seams # ========================================================================== seams_smooth(slicer, smooth_distance=10) # ========================================================================== # Prints out the info of the slicer # ========================================================================== slicer.printout_info() # ========================================================================== # Save slicer data to JSON # ========================================================================== save_to_json(slicer.to_data(), OUTPUT_DIR, 'slicer_data.json') # ========================================================================== # Initializes the PlanarPrintOrganizer and creates PrintPoints # ========================================================================== print_organizer = PlanarPrintOrganizer(slicer) print_organizer.create_printpoints() # ========================================================================== # Set fabrication-related parameters # ========================================================================== set_extruder_toggle(print_organizer, slicer) add_safety_printpoints(print_organizer, z_hop=10.0) set_linear_velocity_constant(print_organizer, v=25.0) set_blend_radius(print_organizer, d_fillet=10.0) # ========================================================================== # Prints out the info of the PrintOrganizer # ========================================================================== print_organizer.printout_info() # ========================================================================== # Converts the PrintPoints to data and saves to JSON # ========================================================================= printpoints_data = print_organizer.output_printpoints_dict() utils.save_to_json(printpoints_data, OUTPUT_DIR, 'out_printpoints.json') # ========================================================================== # Initializes the compas_viewer and visualizes results # ========================================================================== viewer = app.App(width=1600, height=1000) # slicer.visualize_on_viewer(viewer, visualize_mesh=False, visualize_paths=True) print_organizer.visualize_on_viewer(viewer, visualize_printpoints=True) viewer.show() end_time = time.time() print("Total elapsed time", round(end_time - start_time, 2), "seconds")
# ========================================================================== move_mesh_to_point(compas_mesh, Point(0, 0, 0)) # ========================================================================== # Slicing # options: 'default': Both for open and closed paths. But slow # 'cgal': Very fast. Only for closed paths. # Requires additional installation (compas_cgal). # ========================================================================== slicer = PlanarSlicer(compas_mesh, slicer_type="cgal", layer_height=1.5) slicer.slice_model() # ========================================================================== # Generate brim # ========================================================================== generate_brim(slicer, layer_width=3.0, number_of_brim_paths=3) # ========================================================================== # Simplify the paths by removing points with a certain threshold # change the threshold value to remove more or less points # ========================================================================== simplify_paths_rdp(slicer, threshold=0.3) ### ======================================================================= ### ----------------------------------------------------------------------- ### INSERT YOUR CUSTOM FUNCTION BELOW ### ----------------------------------------------------------------------- ### ======================================================================= create_overhang_texture(slicer, overhang_distance=15) # ==========================================================================
def main(): start_time = time.time() # --- Load initial_mesh mesh = Mesh.from_obj(os.path.join(DATA_PATH, OBJ_INPUT_NAME)) # --- Load targets (boundaries) low_boundary_vs = utils.load_from_json(DATA_PATH, 'boundaryLOW.json') high_boundary_vs = utils.load_from_json(DATA_PATH, 'boundaryHIGH.json') create_mesh_boundary_attributes(mesh, low_boundary_vs, high_boundary_vs) avg_layer_height = 15.0 parameters = { 'avg_layer_height': avg_layer_height, # controls number of curves that will be generated 'min_layer_height': 0.3, 'max_layer_height': 5.0 # 2.0, } preprocessor = InterpolationSlicingPreprocessor(mesh, parameters, DATA_PATH) preprocessor.create_compound_targets() g_eval = preprocessor.create_gradient_evaluation( norm_filename='gradient_norm.json', g_filename='gradient.json', target_1=preprocessor.target_LOW, target_2=preprocessor.target_HIGH) preprocessor.find_critical_points( g_eval, output_filenames=['minima.json', 'maxima.json', 'saddles.json']) # --- slicing slicer = InterpolationSlicer(mesh, preprocessor, parameters) slicer.slice_model() # compute_norm_of_gradient contours generate_brim(slicer, layer_width=3.0, number_of_brim_offsets=5) seams_smooth(slicer, smooth_distance=10) simplify_paths_rdp_igl(slicer, threshold=0.5) slicer.printout_info() utils.save_to_json(slicer.to_data(), OUTPUT_PATH, 'curved_slicer.json') # --- Print organizer print_organizer = InterpolationPrintOrganizer(slicer, parameters, DATA_PATH) print_organizer.create_printpoints() set_linear_velocity_by_range( print_organizer, param_func=lambda ppt: ppt.layer_height, parameter_range=[avg_layer_height * 0.5, avg_layer_height * 2.0], velocity_range=[150, 70], bound_remapping=False) set_extruder_toggle(print_organizer, slicer) add_safety_printpoints(print_organizer, z_hop=10.0) smooth_printpoints_up_vectors(print_organizer, strength=0.5, iterations=10) smooth_printpoints_layer_heights(print_organizer, strength=0.5, iterations=5) # --- Save printpoints dictionary to json file printpoints_data = print_organizer.output_printpoints_dict() utils.save_to_json(printpoints_data, OUTPUT_PATH, 'out_printpoints.json') # ----- Visualize viewer = app.App(width=1600, height=1000) # slicer.visualize_on_viewer(viewer, visualize_mesh=False, visualize_paths=True) print_organizer.visualize_on_viewer(viewer, visualize_printpoints=True) viewer.show() end_time = time.time() print("Total elapsed time", round(end_time - start_time, 2), "seconds")
def main(): start_time = time.time() # ========================================================================== # Load mesh # ========================================================================== compas_mesh = Mesh.from_obj(os.path.join(DATA, MODEL)) # ========================================================================== # Move to origin # ========================================================================== move_mesh_to_point(compas_mesh, Point(0, 0, 0)) # ========================================================================== # Slicing # options: 'default': Both for open and closed paths. But slow # 'cgal': Very fast. Only for closed paths. # Requires additional installation (compas_cgal). # ========================================================================== slicer = PlanarSlicer(compas_mesh, slicer_type="default", layer_height=1.5) slicer.slice_model() # ========================================================================== # Generate brim # ========================================================================== generate_brim(slicer, layer_width=3.0, number_of_brim_paths=3) # ========================================================================== # Simplify the paths by removing points with a certain threshold # change the threshold value to remove more or less points # ========================================================================== simplify_paths_rdp(slicer, threshold=0.3) # ========================================================================== # Smooth the seams between layers # change the smooth_distance value to achieve smoother, or more abrupt seams # ========================================================================== seams_smooth(slicer, smooth_distance=10) # ========================================================================== # Prints out the info of the slicer # ========================================================================== slicer.printout_info() # ========================================================================== # Save slicer data to JSON # ========================================================================== save_to_json(slicer.to_data(), OUTPUT_DIR, 'slicer_data.json') # ========================================================================== # Initializes the PlanarPrintOrganizer and creates PrintPoints # ========================================================================== print_organizer = PlanarPrintOrganizer(slicer) print_organizer.create_printpoints() # ========================================================================== # Set fabrication-related parameters # ========================================================================== set_extruder_toggle(print_organizer, slicer) add_safety_printpoints(print_organizer, z_hop=10.0) set_linear_velocity(print_organizer, "constant", v=25.0) set_blend_radius(print_organizer, d_fillet=10.0) # ========================================================================== # Prints out the info of the PrintOrganizer # ========================================================================== slicer.printout_info() print_organizer.printout_info() # ========================================================================== # Converts the PrintPoints to data and saves to JSON # ========================================================================= printpoints_data = print_organizer.output_printpoints_dict() utils.save_to_json(printpoints_data, OUTPUT_DIR, 'out_printpoints.json') # ========================================================================== # Initializes the compas_viewer and visualizes results # ========================================================================== """
def main(): start_time = time.time() avg_layer_height = 4.0 parameters = { 'avg_layer_height': avg_layer_height, # controls number of curves that will be generated 'min_layer_height': 0.2, 'max_layer_height': 4.0, 'uneven_upper_targets_offset': 0, 'target_HIGH_smooth_union': [True, [25.0]], # on/off, blend radius } ### --- Load initial_mesh mesh = Mesh.from_obj(os.path.join(DATA_PATH, OBJ_INPUT_NAME)) # --- Load targets (boundaries) low_boundary_vs = utils.load_from_json(DATA_PATH, 'boundaryLOW.json') high_boundary_vs = utils.load_from_json(DATA_PATH, 'boundaryHIGH.json') create_mesh_boundary_attributes(mesh, low_boundary_vs, high_boundary_vs) # --- Create pre-processor preprocessor = InterpolationSlicingPreprocessor(mesh, parameters, DATA_PATH) preprocessor.create_compound_targets() preprocessor.targets_laplacian_smoothing(iterations=4, strength=0.05) ######################################### # --- region split if REGION_SPLIT: # --- ADVANCED slicing with region split g_eval = preprocessor.create_gradient_evaluation( target_1=preprocessor.target_LOW, target_2=preprocessor.target_HIGH, save_output=True) preprocessor.find_critical_points( g_eval, output_filenames=['minima.json', 'maxima.json', 'saddles.json']) preprocessor.region_split( save_split_meshes=True) # split mesh regions on saddle points # utils.interrupt() ######################################### # --- slicing if SLICER: slicers = [] filenames = utils.get_all_files_with_name('split_mesh_', '.json', OUTPUT_PATH) split_meshes = [ Mesh.from_json(os.path.join(OUTPUT_PATH, filename)) for filename in filenames ] for i, split_mesh in enumerate(split_meshes): preprocessor_split = InterpolationSlicingPreprocessor( split_mesh, parameters, DATA_PATH) preprocessor_split.create_compound_targets() preprocessor_split.create_gradient_evaluation( norm_filename='gradient_norm_%d.json' % i, g_filename='gradient_%d.json' % i, target_1=preprocessor_split.target_LOW, target_2=preprocessor_split.target_HIGH) slicer = InterpolationSlicer(split_mesh, preprocessor_split, parameters) if i == 3: slicer.n_multiplier = 0.85 slicer.slice_model() if i == 0: generate_brim(slicer, layer_width=3.0, number_of_brim_offsets=5) seams_smooth(slicer, smooth_distance=0.1) simplify_paths_rdp_igl(slicer, threshold=0.25) utils.save_to_json(slicer.to_data(), OUTPUT_PATH, 'curved_slicer_%d.json' % i) slicers.append(slicer) # utils.interrupt() ######################################### # --- print organization if PRINT_ORGANIZER: filenames = utils.get_all_files_with_name('curved_slicer_', '.json', OUTPUT_PATH) slicers = [ InterpolationSlicer.from_data( utils.load_from_json(OUTPUT_PATH, filename)) for filename in filenames ] for i, slicer in enumerate(slicers): print_organizer = InterpolationPrintOrganizer( slicer, parameters, DATA_PATH) print_organizer.create_printpoints() set_extruder_toggle(print_organizer, slicer) set_blend_radius(print_organizer) add_safety_printpoints(print_organizer, z_hop=10.0) smooth_printpoints_up_vectors(print_organizer, strength=0.5, iterations=10) set_wait_time_on_sharp_corners(print_organizer, threshold=0.5 * math.pi, wait_time=0.2) # --- Save printpoints dictionary to json file printpoints_data = print_organizer.output_printpoints_dict() utils.save_to_json(printpoints_data, OUTPUT_PATH, 'out_printpoints_%d.json' % i) end_time = time.time() print("Total elapsed time", round(end_time - start_time, 2), "seconds")