private static bool loadModelSTL_binary(SimpleMeshCollection simpleModel, string filename, Matrix4X4 matrix)
{
SimpleMesh vol = new SimpleMesh();
using (FileStream stlStream = File.Open(filename, FileMode.Open, FileAccess.Read, FileShare.ReadWrite))
{
// load it as a binary stl
// skip the first 80 bytes
// read in the number of triangles
stlStream.Position = 0;
BinaryReader br = new BinaryReader(stlStream);
byte[] fileContents = br.ReadBytes((int)stlStream.Length);
int currentPosition = 80;
if (fileContents.Length < currentPosition)
{
return(false);
}
uint numTriangles = System.BitConverter.ToUInt32(fileContents, currentPosition);
long bytesForNormals = numTriangles * 3 * 4;
long bytesForVertices = numTriangles * 3 * 4;
long bytesForAttributs = numTriangles * 2;
currentPosition += 4;
long numBytesRequiredForVertexData = currentPosition + bytesForNormals + bytesForVertices + bytesForAttributs;
if (fileContents.Length < numBytesRequiredForVertexData || numTriangles < 0)
{
stlStream.Close();
return(false);
}
IntPoint[] vector = new IntPoint[3];
for (int i = 0; i < numTriangles; i++)
{
// skip the normal
currentPosition += 3 * 4;
for (int j = 0; j < 3; j++)
{
var vertex = new MatterHackers.VectorMath.Vector3(
System.BitConverter.ToSingle(fileContents, currentPosition + 0 * 4),
System.BitConverter.ToSingle(fileContents, currentPosition + 1 * 4),
System.BitConverter.ToSingle(fileContents, currentPosition + 2 * 4));
var new0 = VectorMath.Vector3.Transform(vertex, matrix);
vector[j] = new IntPoint(new0.X * 1000, new0.Y * 1000, new0.Z * 1000);
currentPosition += 3 * 4;
}
currentPosition += 2; // skip the attribute
vol.addFaceTriangle(vector[2], vector[1], vector[0]);
}
}
if (vol.faceTriangles.Count > 0)
{
simpleModel.SimpleMeshes.Add(vol);
return(true);
}
return(false);
}
public OptimizedMesh(SimpleMesh simpleMesh, OptimizedMeshCollection containingCollection)
{
this.containingCollection = containingCollection;
vertices.Capacity = simpleMesh.faceTriangles.Count * 3;
facesTriangle.Capacity = simpleMesh.faceTriangles.Count;
Dictionary <int, List <int> > indexMap = new Dictionary <int, List <int> >();
Stopwatch t = new Stopwatch();
t.Start();
for (int faceIndex = 0; faceIndex < simpleMesh.faceTriangles.Count; faceIndex++)
{
if (MatterSlice.Canceled)
{
return;
}
OptimizedFace optimizedFace = new OptimizedFace();
if ((faceIndex % 1000 == 0) && t.Elapsed.TotalSeconds > 2)
{
LogOutput.logProgress("optimized", faceIndex + 1, simpleMesh.faceTriangles.Count);
}
for (int vertexIndex = 0; vertexIndex < 3; vertexIndex++)
{
IntPoint p = simpleMesh.faceTriangles[faceIndex].vertices[vertexIndex];
int hash = (int)(((p.X + MELD_DIST / 2) / MELD_DIST) ^ (((p.Y + MELD_DIST / 2) / MELD_DIST) << 10) ^ (((p.Z + MELD_DIST / 2) / MELD_DIST) << 20));
int idx = 0;
bool add = true;
if (indexMap.ContainsKey(hash))
{
for (int n = 0; n < indexMap[hash].Count; n++)
{
if ((vertices[indexMap[hash][n]].position - p).Length() < MELD_DIST)
{
idx = indexMap[hash][n];
add = false;
break;
}
}
}
if (add)
{
if (!indexMap.ContainsKey(hash))
{
indexMap.Add(hash, new List <int>());
}
indexMap[hash].Add(vertices.Count);
idx = vertices.Count;
vertices.Add(new OptimizedPoint3(p));
}
optimizedFace.vertexIndex[vertexIndex] = idx;
}
if (optimizedFace.vertexIndex[0] != optimizedFace.vertexIndex[1] && optimizedFace.vertexIndex[0] != optimizedFace.vertexIndex[2] && optimizedFace.vertexIndex[1] != optimizedFace.vertexIndex[2])
{
//Check if there is a face with the same points
bool duplicate = false;
for (int _idx0 = 0; _idx0 < vertices[optimizedFace.vertexIndex[0]].usedByFacesList.Count; _idx0++)
{
for (int _idx1 = 0; _idx1 < vertices[optimizedFace.vertexIndex[1]].usedByFacesList.Count; _idx1++)
{
for (int _idx2 = 0; _idx2 < vertices[optimizedFace.vertexIndex[2]].usedByFacesList.Count; _idx2++)
{
if (vertices[optimizedFace.vertexIndex[0]].usedByFacesList[_idx0] == vertices[optimizedFace.vertexIndex[1]].usedByFacesList[_idx1] && vertices[optimizedFace.vertexIndex[0]].usedByFacesList[_idx0] == vertices[optimizedFace.vertexIndex[2]].usedByFacesList[_idx2])
{
duplicate = true;
}
}
}
}
if (!duplicate)
{
vertices[optimizedFace.vertexIndex[0]].usedByFacesList.Add(facesTriangle.Count);
vertices[optimizedFace.vertexIndex[1]].usedByFacesList.Add(facesTriangle.Count);
vertices[optimizedFace.vertexIndex[2]].usedByFacesList.Add(facesTriangle.Count);
facesTriangle.Add(optimizedFace);
}
}
}
//fprintf(stdout, "\rAll faces are optimized in %5.1fs.\n",timeElapsed(t));
int openFacesCount = 0;
for (int faceIndex = 0; faceIndex < facesTriangle.Count; faceIndex++)
{
OptimizedFace optimizedFace = facesTriangle[faceIndex];
optimizedFace.touchingFaces[0] = getOtherFaceIndexContainingVertices(optimizedFace.vertexIndex[0], optimizedFace.vertexIndex[1], faceIndex);
optimizedFace.touchingFaces[1] = getOtherFaceIndexContainingVertices(optimizedFace.vertexIndex[1], optimizedFace.vertexIndex[2], faceIndex);
optimizedFace.touchingFaces[2] = getOtherFaceIndexContainingVertices(optimizedFace.vertexIndex[2], optimizedFace.vertexIndex[0], faceIndex);
if (optimizedFace.touchingFaces[0] == -1)
{
openFacesCount++;
}
if (optimizedFace.touchingFaces[1] == -1)
{
openFacesCount++;
}
if (optimizedFace.touchingFaces[2] == -1)
{
openFacesCount++;
}
}
//fprintf(stdout, " Number of open faces: %i\n", openFacesCount);
}
public static bool loadModelSTL_ascii(SimpleMeshCollection simpleModel, string filename, Matrix4X4 matrix)
{
SimpleMesh vol = new SimpleMesh();
using (StreamReader f = new StreamReader(filename))
{
// check for "SOLID"
var vertex = default(MatterHackers.VectorMath.Vector3);
int n = 0;
IntPoint v0 = new IntPoint(0, 0, 0);
IntPoint v1 = new IntPoint(0, 0, 0);
IntPoint v2 = new IntPoint(0, 0, 0);
string line = f.ReadLine();
Regex onlySingleSpaces = new Regex("\\s+", RegexOptions.Compiled);
int lineCount = 0;
while (line != null)
{
if (lineCount++ > 100 && vol.faceTriangles.Count == 0)
{
return(false);
}
line = onlySingleSpaces.Replace(line, " ");
var parts = line.Trim().Split(' ');
if (parts[0].Trim() == "vertex")
{
vertex.X = Convert.ToDouble(parts[1]);
vertex.Y = Convert.ToDouble(parts[2]);
vertex.Z = Convert.ToDouble(parts[3]);
// change the scale from mm to micrometers
n++;
switch (n)
{
case 1:
var new0 = VectorMath.Vector3Ex.Transform(vertex, matrix) * 1000;
v0 = new IntPoint(new0.X, new0.Y, new0.Z);
break;
case 2:
var new1 = VectorMath.Vector3Ex.Transform(vertex, matrix) * 1000;
v1 = new IntPoint(new1.X, new1.Y, new1.Z);
break;
case 3:
var new2 = VectorMath.Vector3Ex.Transform(vertex, matrix) * 1000;
v2 = new IntPoint(new2.X, new2.Y, new2.Z);
vol.addFaceTriangle(v0, v1, v2);
n = 0;
break;
}
}
line = f.ReadLine();
}
}
if (vol.faceTriangles.Count > 3)
{
simpleModel.SimpleMeshes.Add(vol);
return(true);
}
return(false);
}
private static bool loadModelSTLStream_binary(SimpleMeshCollection simpleModel, Stream stlStream, Matrix4X4 matrix)
{
SimpleMesh vol = new SimpleMesh();
// load it as a binary stl
// skip the first 80 bytes
// read in the number of triangles
stlStream.Position = 0;
BinaryReader br = new BinaryReader(stlStream);
byte[] fileContents = br.ReadBytes((int)stlStream.Length);
int currentPosition = 80;
if (fileContents.Length < currentPosition)
{
return(false);
}
uint numTriangles = System.BitConverter.ToUInt32(fileContents, currentPosition);
long bytesForNormals = numTriangles * 3 * 4;
long bytesForVertices = numTriangles * 3 * 4;
long bytesForAttributs = numTriangles * 2;
currentPosition += 4;
long numBytesRequiredForVertexData = currentPosition + bytesForNormals + bytesForVertices + bytesForAttributs;
if (fileContents.Length < numBytesRequiredForVertexData || numTriangles < 0)
{
stlStream.Close();
return(false);
}
IntPoint[] vector = new IntPoint[3];
for (int i = 0; i < numTriangles; i++)
{
// skip the normal
currentPosition += 3 * 4;
for (int j = 0; j < 3; j++)
{
var vertex = new MatterHackers.VectorMath.Vector3(
System.BitConverter.ToSingle(fileContents, currentPosition + 0 * 4),
System.BitConverter.ToSingle(fileContents, currentPosition + 1 * 4),
System.BitConverter.ToSingle(fileContents, currentPosition + 2 * 4));
var new0 = VectorMath.Vector3Ex.Transform(vertex, matrix);
vector[j] = new IntPoint(new0.X * 1000, new0.Y * 1000, new0.Z * 1000);
currentPosition += 3 * 4;
}
currentPosition += 2; // skip the attribute
vol.addFaceTriangle(vector[2], vector[1], vector[0]);
}
// Detect and skip non-visible mesh
var bounds = vol.maxXYZ_um() - vol.minXYZ_um();
if (vol.faceTriangles.Count > 0)
{
if (bounds.X == 0)
{
vol.faceTriangles = new List <SimpleFace>();
}
simpleModel.SimpleMeshes.Add(vol);
return(true);
}
return(false);
}
public OptimizedMesh(SimpleMesh simpleMesh, OptimizedMeshCollection containingCollection)
{
this.containingCollection = containingCollection;
vertices.Capacity = simpleMesh.faceTriangles.Count * 3;
facesTriangle.Capacity = simpleMesh.faceTriangles.Count;
Dictionary<int, List<int>> indexMap = new Dictionary<int, List<int>>();
Stopwatch t = new Stopwatch();
t.Start();
for (int faceIndex = 0; faceIndex < simpleMesh.faceTriangles.Count; faceIndex++)
{
if (MatterSlice.Canceled)
{
return;
}
OptimizedFace optimizedFace = new OptimizedFace();
if ((faceIndex % 1000 == 0) && t.Elapsed.TotalSeconds > 2)
{
LogOutput.logProgress("optimized", faceIndex + 1, simpleMesh.faceTriangles.Count);
}
for (int vertexIndex = 0; vertexIndex < 3; vertexIndex++)
{
Point3 p = simpleMesh.faceTriangles[faceIndex].vertices[vertexIndex];
int hash = (int)(((p.x + MELD_DIST / 2) / MELD_DIST) ^ (((p.y + MELD_DIST / 2) / MELD_DIST) << 10) ^ (((p.z + MELD_DIST / 2) / MELD_DIST) << 20));
int idx = 0;
bool add = true;
if (indexMap.ContainsKey(hash))
{
for (int n = 0; n < indexMap[hash].Count; n++)
{
if ((vertices[indexMap[hash][n]].position - p).AbsLengthLEQ(MELD_DIST))
{
idx = indexMap[hash][n];
add = false;
break;
}
}
}
if (add)
{
if (!indexMap.ContainsKey(hash))
{
indexMap.Add(hash, new List<int>());
}
indexMap[hash].Add(vertices.Count);
idx = vertices.Count;
vertices.Add(new OptimizedPoint3(p));
}
optimizedFace.vertexIndex[vertexIndex] = idx;
}
if (optimizedFace.vertexIndex[0] != optimizedFace.vertexIndex[1] && optimizedFace.vertexIndex[0] != optimizedFace.vertexIndex[2] && optimizedFace.vertexIndex[1] != optimizedFace.vertexIndex[2])
{
//Check if there is a face with the same points
bool duplicate = false;
for (int _idx0 = 0; _idx0 < vertices[optimizedFace.vertexIndex[0]].usedByFacesList.Count; _idx0++)
{
for (int _idx1 = 0; _idx1 < vertices[optimizedFace.vertexIndex[1]].usedByFacesList.Count; _idx1++)
{
for (int _idx2 = 0; _idx2 < vertices[optimizedFace.vertexIndex[2]].usedByFacesList.Count; _idx2++)
{
if (vertices[optimizedFace.vertexIndex[0]].usedByFacesList[_idx0] == vertices[optimizedFace.vertexIndex[1]].usedByFacesList[_idx1] && vertices[optimizedFace.vertexIndex[0]].usedByFacesList[_idx0] == vertices[optimizedFace.vertexIndex[2]].usedByFacesList[_idx2])
duplicate = true;
}
}
}
if (!duplicate)
{
vertices[optimizedFace.vertexIndex[0]].usedByFacesList.Add(facesTriangle.Count);
vertices[optimizedFace.vertexIndex[1]].usedByFacesList.Add(facesTriangle.Count);
vertices[optimizedFace.vertexIndex[2]].usedByFacesList.Add(facesTriangle.Count);
facesTriangle.Add(optimizedFace);
}
}
}
//fprintf(stdout, "\rAll faces are optimized in %5.1fs.\n",timeElapsed(t));
int openFacesCount = 0;
for (int faceIndex = 0; faceIndex < facesTriangle.Count; faceIndex++)
{
OptimizedFace optimizedFace = facesTriangle[faceIndex];
optimizedFace.touchingFaces[0] = getOtherFaceIndexContainingVertices(optimizedFace.vertexIndex[0], optimizedFace.vertexIndex[1], faceIndex);
optimizedFace.touchingFaces[1] = getOtherFaceIndexContainingVertices(optimizedFace.vertexIndex[1], optimizedFace.vertexIndex[2], faceIndex);
optimizedFace.touchingFaces[2] = getOtherFaceIndexContainingVertices(optimizedFace.vertexIndex[2], optimizedFace.vertexIndex[0], faceIndex);
if (optimizedFace.touchingFaces[0] == -1)
{
openFacesCount++;
}
if (optimizedFace.touchingFaces[1] == -1)
{
openFacesCount++;
}
if (optimizedFace.touchingFaces[2] == -1)
{
openFacesCount++;
}
}
//fprintf(stdout, " Number of open faces: %i\n", openFacesCount);
}
private static bool loadModelSTL_binary(SimpleMeshCollection simpleModel, string filename, FMatrix3x3 matrix)
{
SimpleMesh vol = new SimpleMesh();
using (FileStream stlStream = File.Open(filename, FileMode.Open, FileAccess.Read, FileShare.ReadWrite))
{
// load it as a binary stl
// skip the first 80 bytes
// read in the number of triangles
stlStream.Position = 0;
BinaryReader br = new BinaryReader(stlStream);
byte[] fileContents = br.ReadBytes((int)stlStream.Length);
int currentPosition = 80;
uint numTriangles = System.BitConverter.ToUInt32(fileContents, currentPosition);
long bytesForNormals = numTriangles * 3 * 4;
long bytesForVertices = numTriangles * 3 * 4;
long bytesForAttributs = numTriangles * 2;
currentPosition += 4;
long numBytesRequiredForVertexData = currentPosition + bytesForNormals + bytesForVertices + bytesForAttributs;
if (fileContents.Length < numBytesRequiredForVertexData || numTriangles < 0)
{
stlStream.Close();
return false;
}
Point3[] vector = new Point3[3];
for (int i = 0; i < numTriangles; i++)
{
// skip the normal
currentPosition += 3 * 4;
for (int j = 0; j < 3; j++)
{
vector[j] = new Point3(
System.BitConverter.ToSingle(fileContents, currentPosition + 0 * 4) * 1000,
System.BitConverter.ToSingle(fileContents, currentPosition + 1 * 4) * 1000,
System.BitConverter.ToSingle(fileContents, currentPosition + 2 * 4) * 1000);
currentPosition += 3 * 4;
}
currentPosition += 2; // skip the attribute
vol.addFaceTriangle(vector[2], vector[1], vector[0]);
}
}
if (vol.faceTriangles.Count > 0)
{
simpleModel.SimpleMeshes.Add(vol);
return true;
}
return false;
}
public static bool loadModelSTL_ascii(SimpleMeshCollection simpleModel, string filename, FMatrix3x3 matrix)
{
SimpleMesh vol = new SimpleMesh();
using (StreamReader f = new StreamReader(filename))
{
// check for "SOLID"
Vector3 vertex = new Vector3();
int n = 0;
Point3 v0 = new Point3(0, 0, 0);
Point3 v1 = new Point3(0, 0, 0);
Point3 v2 = new Point3(0, 0, 0);
string line = f.ReadLine();
Regex onlySingleSpaces = new Regex("\\s+", RegexOptions.Compiled);
int lineCount = 0;
while (line != null)
{
if(lineCount++ > 100 && vol.faceTriangles.Count == 0)
{
return false;
}
line = onlySingleSpaces.Replace(line, " ");
var parts = line.Trim().Split(' ');
if (parts[0].Trim() == "vertex")
{
vertex.x = Convert.ToDouble(parts[1]);
vertex.y = Convert.ToDouble(parts[2]);
vertex.z = Convert.ToDouble(parts[3]);
// change the scale from mm to micrometers
vertex *= 1000.0;
n++;
switch (n)
{
case 1:
v0 = matrix.apply(vertex);
break;
case 2:
v1 = matrix.apply(vertex);
break;
case 3:
v2 = matrix.apply(vertex);
vol.addFaceTriangle(v0, v1, v2);
n = 0;
break;
}
}
line = f.ReadLine();
}
}
if (vol.faceTriangles.Count > 3)
{
simpleModel.SimpleMeshes.Add(vol);
return true;
}
return false;
}
public static bool loadModelSTL_ascii(SimpleMeshCollection simpleModel, string filename, FMatrix3x3 matrix)
{
SimpleMesh vol = new SimpleMesh();
using (StreamReader f = new StreamReader(filename))
{
// check for "SOLID"
Vector3 vertex = new Vector3();
int n = 0;
Point3 v0 = new Point3(0, 0, 0);
Point3 v1 = new Point3(0, 0, 0);
Point3 v2 = new Point3(0, 0, 0);
string line = f.ReadLine();
Regex onlySingleSpaces = new Regex("\\s+", RegexOptions.Compiled);
int lineCount = 0;
while (line != null)
{
if (lineCount++ > 100 && vol.faceTriangles.Count == 0)
{
return(false);
}
line = onlySingleSpaces.Replace(line, " ");
var parts = line.Trim().Split(' ');
if (parts[0].Trim() == "vertex")
{
vertex.x = Convert.ToDouble(parts[1]);
vertex.y = Convert.ToDouble(parts[2]);
vertex.z = Convert.ToDouble(parts[3]);
// change the scale from mm to micrometers
vertex *= 1000.0;
n++;
switch (n)
{
case 1:
v0 = matrix.apply(vertex);
break;
case 2:
v1 = matrix.apply(vertex);
break;
case 3:
v2 = matrix.apply(vertex);
vol.addFaceTriangle(v0, v1, v2);
n = 0;
break;
}
}
line = f.ReadLine();
}
}
if (vol.faceTriangles.Count > 3)
{
simpleModel.SimpleMeshes.Add(vol);
return(true);
}
return(false);
}
