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():
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():
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")
for v_key, data in mesh.vertices(data=True):
v_coord = mesh.vertex_coordinates(v_key, axes='xyz')
data['dist_from_plane'] = distance_point_plane(v_coord, plane)
# direction towards point - Vector value (per vertex)
mesh.update_default_vertex_attributes({'direction_to_pt': 0.0})
pt = Point(4.0, 1.0, 0.0)
for v_key, data in mesh.vertices(data=True):
v_coord = mesh.vertex_coordinates(v_key, axes='xyz')
data['direction_to_pt'] = np.array(
normalize_vector(Vector.from_start_end(v_coord, pt)))
# --------------- Slice mesh
slicer = PlanarSlicer(mesh, slicer_type="default", layer_height=5.0)
slicer.slice_model()
simplify_paths_rdp_igl(slicer, threshold=1.0)
slicer_utils.save_to_json(slicer.to_data(), OUTPUT_PATH,
'slicer_data.json')
# --------------- Create printpoints
print_organizer = PlanarPrintOrganizer(slicer)
print_organizer.create_printpoints()
# --------------- Transfer mesh attributes to printpoints
transfer_mesh_attributes_to_printpoints(mesh,
print_organizer.printpoints_dict)
# --------------- Save printpoints to json (only json-serializable attributes are saved)
printpoints_data = print_organizer.output_printpoints_dict()
utils.save_to_json(printpoints_data, OUTPUT_PATH, 'out_printpoints.json')
distance_point_point(pt, proj_pt)
for pt, proj_pt in zip(pts, projected_pts)
] # list with distance between initial+projected pts (one per vertex)
utils.save_to_json(
u, OUTPUT_PATH,
'distance_field.json') # save distance field to json for visualization
# --- assign the scalar field to the mesh's attributes, under the name 'scalar_field'
mesh.update_default_vertex_attributes({'scalar_field': 0.0})
for i, (v_key, data) in enumerate(mesh.vertices(data=True)):
data['scalar_field'] = u[i]
# --- Slice model by generating contours of scalar field
slicer = ScalarFieldSlicer(mesh, u, no_of_isocurves=50)
slicer.slice_model()
# simplify_paths_rdp(slicer, threshold=0.3)
slicer_utils.save_to_json(slicer.to_data(), OUTPUT_PATH,
'isocontours.json') # save results to json
# --- Print organization calculations (i.e. generation of printpoints with fabrication-related information)
simplify_paths_rdp_igl(slicer, threshold=0.3)
print_organizer = ScalarFieldPrintOrganizer(slicer,
parameters={},
DATA_PATH=DATA_PATH)
print_organizer.create_printpoints()
print_organizer.printout_info()
printpoints_data = print_organizer.output_printpoints_dict()
utils.save_to_json(printpoints_data, OUTPUT_PATH,
'out_printpoints.json') # save results to json
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()
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")