openscad/dotSCAD
1
0
Fork 0
mirror of https://github.com/JustinSDK/dotSCAD.git synced 2025-08-17 20:11:50 +02:00
Code Issues Releases Wiki Activity

refactor

Browse Source
This commit is contained in:
Justin Lin
2020-03-22 17:01:38 +08:00
parent c2fad6fb2f
commit 68291ea20c
2 changed files with 39 additions and 28 deletions
Download Patch File Download Diff File Expand all files Collapse all files

48
src/experimental/_impl/_marching_squares_impl.scad
View File

@@ -1,3 +1,17 @@
function interpolated_pt(p0, p1, sigma) =
let(
x0 = p0[0],
y0 = p0[1],
z0 = p0[2],
v = p1 - p0,
t = (sigma - z0) / v[2]
)
[x0 + v[0] * t, y0 + v[1] * t, sigma];
/*
ISOLINES Impl Begin ============================
*/
function _isolines_pn_label(pts, sigma) =
[
for(row = pts)
@@ -13,16 +27,6 @@ function _isolines_corner_value(cell_pts) =
)
c0 + c1 + c2 + c3;
function interpolated_pt(p0, p1, sigma) =
let(
x0 = p0[0],
y0 = p0[1],
z0 = p0[2],
v = p1 - p0,
t = (sigma - z0) / v[2]
)
[x0 + v[0] * t, y0 + v[1] * t, sigma];
function _case1_isolines(cell_pts, sigma) = [
[interpolated_pt(cell_pts[0], cell_pts[1], sigma), interpolated_pt(cell_pts[0], cell_pts[3], sigma)]
];
@@ -102,3 +106,27 @@ function _isolines_of(cell_pts, sigma) =
cv == 6 ? _case12_isolines(cell_pts, sigma) :
cv == 7 ? _case13_isolines(cell_pts, sigma) :
_case14_isolines(cell_pts, sigma);
function _marching_squares_isolines(points, sigma) =
let(labeled_pts = _isolines_pn_label(points, sigma))
[
for(y = [0:len(labeled_pts) - 2])
[
for(x = [0:len(labeled_pts[0]) - 2])
let(
p0 = labeled_pts[y][x],
p1 = labeled_pts[y + 1][x],
p2 = labeled_pts[y + 1][x + 1],
p3 = labeled_pts[y][x + 1],
cell_pts = [p0, p1, p2, p3],
isolines_lt = _isolines_of(cell_pts, sigma)
)
if(isolines_lt != [])
each isolines_lt
]
];
/*
ISOLINES Impl End ============================
*/

19
src/experimental/marching_squares.scad
View File

@@ -1,20 +1,3 @@
use <experimental/_impl/_marching_squares_impl.scad>;
function marching_squares(points, sigma) =
let(labeled_pts = _isolines_pn_label(points, sigma))
[
for(y = [0:len(labeled_pts) - 2])
[
for(x = [0:len(labeled_pts[0]) - 2])
let(
p0 = labeled_pts[y][x],
p1 = labeled_pts[y + 1][x],
p2 = labeled_pts[y + 1][x + 1],
p3 = labeled_pts[y][x + 1],
cell_pts = [p0, p1, p2, p3],
isolines_lt = _isolines_of(cell_pts, sigma)
)
if(isolines_lt != [])
each isolines_lt
]
];
function marching_squares(points, sigma) = _marching_squares_isolines(points, sigma);