public int CompareTo(Ray3d other)
{
int cmp = Origin.CompareTo(other.Origin);
if (cmp == 0)
{
return(Direction.CompareTo(other.Direction));
}
return(cmp);
}
public int CompareTo(Line3d other)
{
int cmp = A.CompareTo(other.A);
if (cmp != 0)
{
return(cmp);
}
return(B.CompareTo(other.B));
}
public int CompareTo(AABoxd other)
{
int cmp = Min.CompareTo(other.Min);
if (cmp != 0)
{
return(cmp);
}
return(Max.CompareTo(other.Max));
}
public int CompareTo(Sphered other)
{
int cmp = Center.CompareTo(other.Center);
if (cmp != 0)
{
return(cmp);
}
return(Radius.CompareTo(other.Radius));
}
public int CompareTo(Bezier3d other)
{
int cmp = A.CompareTo(other.A);
if (cmp != 0)
{
return(cmp);
}
cmp = B.CompareTo(other.B);
if (cmp != 0)
{
return(cmp);
}
return(C.CompareTo(other.C));
}
public static DMesh3 RemeshMesh(g3.DMesh3 mesh, float minEdgeLength, float maxEdgeLength, float contraintAngle, float smoothSpeed, int smoothPasses, List <Line3d> constrainedLines)
{
// construct mesh projection target
DMesh3 meshCopy = new DMesh3(mesh);
DMeshAABBTree3 tree = new DMeshAABBTree3(meshCopy);
tree.Build();
MeshProjectionTarget target = new MeshProjectionTarget()
{
Mesh = meshCopy,
Spatial = tree
};
MeshConstraints cons = new MeshConstraints();
EdgeRefineFlags useFlags = EdgeRefineFlags.NoFlip;
foreach (int eid in mesh.EdgeIndices())
{
double fAngle = MeshUtil.OpeningAngleD(mesh, eid);
Index2i ev = mesh.GetEdgeV(eid);
if (fAngle > contraintAngle)
{
cons.SetOrUpdateEdgeConstraint(eid, new EdgeConstraint(useFlags));
// TODO Ids based off of ?? What?
int nSetID0 = (mesh.GetVertex(ev[0]).y > 1) ? 1 : 2;
int nSetID1 = (mesh.GetVertex(ev[1]).y > 1) ? 1 : 2;
cons.SetOrUpdateVertexConstraint(ev[0], new VertexConstraint(true, nSetID0));
cons.SetOrUpdateVertexConstraint(ev[1], new VertexConstraint(true, nSetID1));
}
Vector3d p1 = mesh.GetVertex(ev.a);
Vector3d p2 = mesh.GetVertex(ev.b);
foreach (var v in constrainedLines)
{
if (p1.CompareTo(v.Origin) == 0)
{
Vector3d p = v.PointAt(1.0);
if (p2.CompareTo(p) == 0)
{
cons.SetOrUpdateEdgeConstraint(eid, EdgeConstraint.FullyConstrained);
break;
}
}
}
foreach (var v in constrainedLines)
{
if (p2.CompareTo(v.Origin) == 0)
{
Vector3d p = v.PointAt(1.0);
if (p1.CompareTo(p) == 0)
{
cons.SetOrUpdateEdgeConstraint(eid, EdgeConstraint.FullyConstrained);
break;
}
}
}
}
Remesher r = new Remesher(mesh);
r.SetExternalConstraints(cons);
r.SetProjectionTarget(target);
r.Precompute();
r.EnableFlips = r.EnableSplits = r.EnableCollapses = true;
r.MinEdgeLength = minEdgeLength; //0.1f;
r.MaxEdgeLength = maxEdgeLength; // 0.2f;
r.EnableSmoothing = true;
r.SmoothSpeedT = smoothSpeed; // .5;
for (int k = 0; k < smoothPasses; ++k) // smoothPasses = 20
{
r.BasicRemeshPass();
}
return(mesh);
}
public static DMesh3 RemeshMesh(g3.DMesh3 mesh, float minEdgeLength, float maxEdgeLength, float contraintAngle, float smoothSpeed, int smoothPasses, g3.DMesh3 projectMeshInput = null, float projectAmount = 1.0f, float projectedDistance = float.MaxValue, List <Line3d> constrainedLines = null)
{
g3.DMesh3 projectMesh = projectMeshInput;
if (projectMesh == null)
{
projectMesh = mesh;
}
DMesh3 projectMeshCopy = new DMesh3(projectMesh);
DMeshAABBTree3 treeProject = new DMeshAABBTree3(projectMeshCopy);
treeProject.Build();
MeshProjectionTarget targetProject = new MeshProjectionTarget()
{
Mesh = projectMeshCopy,
Spatial = treeProject
};
MeshConstraints cons = new MeshConstraints();
EdgeRefineFlags useFlags = EdgeRefineFlags.NoFlip;
foreach (int eid in mesh.EdgeIndices())
{
double fAngle = MeshUtil.OpeningAngleD(mesh, eid);
if (fAngle > contraintAngle)
{
cons.SetOrUpdateEdgeConstraint(eid, new EdgeConstraint(useFlags));
Index2i ev = mesh.GetEdgeV(eid);
//int nSetID0 = (mesh.GetVertex(ev[0]).y > 1) ? 1 : 2;
//int nSetID1 = (mesh.GetVertex(ev[1]).y > 1) ? 1 : 2;
cons.SetOrUpdateVertexConstraint(ev[0], new VertexConstraint(true));
cons.SetOrUpdateVertexConstraint(ev[1], new VertexConstraint(true));
}
if (constrainedLines != null)
{
Index2i ev = mesh.GetEdgeV(eid);
Vector3d p1 = mesh.GetVertex(ev.a);
Vector3d p2 = mesh.GetVertex(ev.b);
foreach (var v in constrainedLines)
{
if (p1.CompareTo(v.Origin) == 0)
{
Vector3d p = v.PointAt(1.0);
if (p2.CompareTo(p) == 0)
{
cons.SetOrUpdateEdgeConstraint(eid, EdgeConstraint.FullyConstrained);
break;
}
}
}
foreach (var v in constrainedLines)
{
if (p2.CompareTo(v.Origin) == 0)
{
Vector3d p = v.PointAt(1.0);
if (p1.CompareTo(p) == 0)
{
cons.SetOrUpdateEdgeConstraint(eid, EdgeConstraint.FullyConstrained);
break;
}
}
}
}
}
// TODO Constrain Vertices too far away
foreach (int vid in mesh.VertexIndices())
{
var v = mesh.GetVertex(vid);
//v.Distance()
//targetProject.Project()
}
Remesher rProjected = new Remesher(mesh);
rProjected.SetExternalConstraints(cons);
rProjected.SetProjectionTarget(targetProject);
rProjected.Precompute();
rProjected.EnableFlips = rProjected.EnableSplits = rProjected.EnableCollapses = true;
rProjected.MinEdgeLength = minEdgeLength; //0.1f;
rProjected.MaxEdgeLength = maxEdgeLength; // 0.2f;
rProjected.EnableSmoothing = true;
rProjected.SmoothSpeedT = smoothSpeed; // .5;
if (projectMeshInput != null)
{
float bestSmoothPassProjectAmount = projectAmount / smoothPasses;
float testbestSmoothPassProjectAmount = float.MaxValue;
for (float smoothPassProjectAmount = -.1f; smoothPassProjectAmount < 1.1f; smoothPassProjectAmount += 0.005f)
{
double test = 0;
for (int i = 0; i < smoothPasses; i++)
{
test = 1.0 * smoothPassProjectAmount + test * (1 - smoothPassProjectAmount);
}
if (Math.Abs(test - projectAmount) < Math.Abs(testbestSmoothPassProjectAmount - projectAmount))
{
bestSmoothPassProjectAmount = (float)smoothPassProjectAmount;
testbestSmoothPassProjectAmount = (float)test;
}
}
for (int k = 0; k < smoothPasses; ++k) // smoothPasses = 20
{
rProjected.FullProjectionPass(bestSmoothPassProjectAmount, projectedDistance);
}
}
else
{
for (int k = 0; k < smoothPasses; ++k) // smoothPasses = 20
{
rProjected.BasicRemeshPass();
}
}
return(mesh);
}
