Skip to content

Commit

Permalink
Slicer less strict to detect inner volume of a mesh
Browse files Browse the repository at this point in the history
  • Loading branch information
janper committed Apr 24, 2021
1 parent 1e17948 commit 91f0695
Show file tree
Hide file tree
Showing 2 changed files with 11 additions and 10 deletions.
21 changes: 11 additions & 10 deletions Components/SliceGeometry.cs
View file
Original file line number Diff line number Diff line change
Expand Up @@ -218,9 +218,9 @@ protected override void SolveInstance(IGH_DataAccess DA) {
}

private static bool IsOnEdgeUnitized(Point3d geometryPoint) {
return Math.Abs(Math.Abs(geometryPoint.X % 1) - 0.5) <= RhinoMath.SqrtEpsilon
|| Math.Abs(Math.Abs(geometryPoint.Y % 1) - 0.5) <= RhinoMath.SqrtEpsilon
|| Math.Abs(Math.Abs(geometryPoint.Z % 1) - 0.5) <= RhinoMath.SqrtEpsilon;
return Math.Abs(Math.Abs(geometryPoint.X % 1) - 0.5) <= Config.EPSILON
|| Math.Abs(Math.Abs(geometryPoint.Y % 1) - 0.5) <= Config.EPSILON
|| Math.Abs(Math.Abs(geometryPoint.Z % 1) - 0.5) <= Config.EPSILON;
}

/// <summary>
Expand Down Expand Up @@ -281,17 +281,18 @@ private static IEnumerable<Point3d> PopulateCurve(double distance, Curve curve)
/// entirely inside mesh volume.</returns>
private static List<Point3i> CentersFromMeshVolume(Mesh mesh) {
var pointsInsideMesh = new List<Point3i>();
var boundingBox = mesh.GetBoundingBox(false);
var boundingBox = mesh.GetBoundingBox(true);
var slotInterval = new Interval(-0.5, 0.5);
var plane = Plane.WorldXY;
for (var z = Math.Floor(boundingBox.Min.Z - 1); z < Math.Ceiling(boundingBox.Max.Z + 1); z++) {
for (var y = Math.Floor(boundingBox.Min.Y - 1); y < Math.Ceiling(boundingBox.Max.Y + 1); y++) {
for (var x = Math.Floor(boundingBox.Min.X - 1); x < Math.Ceiling(boundingBox.Max.X + 1); x++) {
var testSlotCenter = new Point3d(x, y, z);
var plane = Plane.WorldXY;
plane.Origin = testSlotCenter;
var box = new Box(plane, slotInterval, slotInterval, slotInterval);
var boxPoints = box.GetCorners();
if (boxPoints.All(point => mesh.IsPointInside(point, RhinoMath.SqrtEpsilon, false))) {
if (boxPoints.All(point => mesh.IsPointInside(point, Config.EPSILON, false)
|| point.DistanceTo(mesh.ClosestPoint(point)) < Config.EPSILON)) {
pointsInsideMesh.Add(new Point3i(testSlotCenter));
}
}
Expand All @@ -312,15 +313,15 @@ private static List<Point3i> CentersFromMeshVolumeAndSurface(Mesh mesh) {
var pointsInsideMesh = new List<Point3i>();
var boundingBox = mesh.GetBoundingBox(false);
var slotInterval = new Interval(-0.5, 0.5);
var plane = Plane.WorldXY;
for (var z = Math.Floor(boundingBox.Min.Z - 1); z < Math.Ceiling(boundingBox.Max.Z + 1); z++) {
for (var y = Math.Floor(boundingBox.Min.Y - 1); y < Math.Ceiling(boundingBox.Max.Y + 1); y++) {
for (var x = Math.Floor(boundingBox.Min.X - 1); x < Math.Ceiling(boundingBox.Max.X + 1); x++) {
var testSlotCenter = new Point3d(x, y, z);
var plane = Plane.WorldXY;
plane.Origin = testSlotCenter;
var box = new Box(plane, slotInterval, slotInterval, slotInterval);
var boxPoints = box.GetCorners();
if (boxPoints.Any(point => mesh.IsPointInside(point, Rhino.RhinoMath.SqrtEpsilon, false))) {
if (boxPoints.Any(point => mesh.IsPointInside(point, Config.EPSILON, false) || point.DistanceTo(mesh.ClosestPoint(point)) < Config.EPSILON)) {
pointsInsideMesh.Add(new Point3i(testSlotCenter));
}
}
Expand All @@ -340,16 +341,16 @@ private static List<Point3i> CentersFromMeshSurface(Mesh mesh) {
var pointsInsideMesh = new List<Point3i>();
var boundingBox = mesh.GetBoundingBox(false);
var slotInterval = new Interval(-0.5, 0.5);
var plane = Plane.WorldXY;
for (var z = Math.Floor(boundingBox.Min.Z - 1); z < Math.Ceiling(boundingBox.Max.Z + 1); z++) {
for (var y = Math.Floor(boundingBox.Min.Y - 1); y < Math.Ceiling(boundingBox.Max.Y + 1); y++) {
for (var x = Math.Floor(boundingBox.Min.X - 1); x < Math.Ceiling(boundingBox.Max.X + 1); x++) {
var testSlotCenter = new Point3d(x, y, z);
var plane = Plane.WorldXY;
plane.Origin = testSlotCenter;
var box = new Box(plane, slotInterval, slotInterval, slotInterval);
var boxPoints = box.GetCorners();
var boxcornersInside = boxPoints
.Count(point => mesh.IsPointInside(point, RhinoMath.SqrtEpsilon, false));
.Count(point => mesh.IsPointInside(point, Config.EPSILON, false) || point.DistanceTo(mesh.ClosestPoint(point)) < Config.EPSILON);
if (boxcornersInside > 0 && boxcornersInside < 6) {
pointsInsideMesh.Add(new Point3i(testSlotCenter));
}
Expand Down
Binary file modified bin/Monoceros.pdb
View file
Binary file not shown.

0 comments on commit 91f0695

Please sign in to comment.